Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Static Public Member Functions | Static Public Attributes | Protected Member Functions | Private Member Functions | Friends | List of all members
Nektar::IsentropicVortexBC Class Reference

Wall boundary conditions for compressible flow problems. More...

#include <IsentropicVortexBC.h>

Inheritance diagram for Nektar::IsentropicVortexBC:
Inheritance graph
[legend]
Collaboration diagram for Nektar::IsentropicVortexBC:
Collaboration graph
[legend]

Static Public Member Functions

static CFSBndCondSharedPtr create (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const int pSpaceDim, const int bcRegion, const int cnt)
 Creates an instance of this class. More...
 

Static Public Attributes

static std::string className
 Name of the class. More...
 

Protected Member Functions

virtual void v_Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time)
 
- Protected Member Functions inherited from Nektar::CFSBndCond
 CFSBndCond (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const int pSpaceDim, const int bcRegion, const int cnt)
 Constructor. More...
 

Private Member Functions

 IsentropicVortexBC (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const int pSpaceDim, const int bcRegion, const int cnt)
 
void EvaluateIsentropicVortex (const Array< OneD, NekDouble > &x, const Array< OneD, NekDouble > &y, const Array< OneD, NekDouble > &z, Array< OneD, Array< OneD, NekDouble > > &u, NekDouble time, const int o=0)
 
virtual ~IsentropicVortexBC (void)
 

Friends

class MemoryManager< IsentropicVortexBC >
 

Additional Inherited Members

- Public Member Functions inherited from Nektar::CFSBndCond
virtual ~CFSBndCond ()
 
void Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time=0)
 Apply the boundary condition. More...
 
- Protected Attributes inherited from Nektar::CFSBndCond
LibUtilities::SessionReaderSharedPtr m_session
 Session reader. More...
 
Array< OneD,
MultiRegions::ExpListSharedPtr
m_fields
 Array of fields. More...
 
Array< OneD, Array< OneD,
NekDouble > > 
m_traceNormals
 Trace normals. More...
 
int m_spacedim
 Space dimension. More...
 
VariableConverterSharedPtr m_varConv
 Auxiliary object to convert variables. More...
 
NekDouble m_gamma
 Parameters of the flow. More...
 
NekDouble m_rhoInf
 
NekDouble m_pInf
 
Array< OneD, NekDoublem_velInf
 
int m_bcRegion
 Id of the boundary region. More...
 
int m_offset
 Offset. More...
 

Detailed Description

Wall boundary conditions for compressible flow problems.

Definition at line 48 of file IsentropicVortexBC.h.

Constructor & Destructor Documentation

Nektar::IsentropicVortexBC::IsentropicVortexBC ( const LibUtilities::SessionReaderSharedPtr pSession,
const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const Array< OneD, Array< OneD, NekDouble > > &  pTraceNormals,
const int  pSpaceDim,
const int  bcRegion,
const int  cnt 
)
private

Definition at line 48 of file IsentropicVortexBC.cpp.

54  : CFSBndCond(pSession, pFields, pTraceNormals, pSpaceDim, bcRegion, cnt)
55 {
56 }
CFSBndCond(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const int pSpaceDim, const int bcRegion, const int cnt)
Constructor.
Definition: CFSBndCond.cpp:51
virtual Nektar::IsentropicVortexBC::~IsentropicVortexBC ( void  )
inlineprivatevirtual

Definition at line 93 of file IsentropicVortexBC.h.

93 {};

Member Function Documentation

static CFSBndCondSharedPtr Nektar::IsentropicVortexBC::create ( const LibUtilities::SessionReaderSharedPtr pSession,
const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const Array< OneD, Array< OneD, NekDouble > > &  pTraceNormals,
const int  pSpaceDim,
const int  bcRegion,
const int  cnt 
)
inlinestatic

Creates an instance of this class.

Definition at line 55 of file IsentropicVortexBC.h.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and CellMLToNektar.cellml_metadata::p.

60  {
62  AllocateSharedPtr(pSession, pFields,
63  pTraceNormals, pSpaceDim, bcRegion, cnt);
64  return p;
65  }
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
boost::shared_ptr< CFSBndCond > CFSBndCondSharedPtr
A shared pointer to a boundary condition object.
Definition: CFSBndCond.h:49
void Nektar::IsentropicVortexBC::EvaluateIsentropicVortex ( const Array< OneD, NekDouble > &  x,
const Array< OneD, NekDouble > &  y,
const Array< OneD, NekDouble > &  z,
Array< OneD, Array< OneD, NekDouble > > &  u,
NekDouble  time,
const int  o = 0 
)
private

Definition at line 97 of file IsentropicVortexBC.cpp.

References Nektar::CFSBndCond::m_gamma, Nektar::CFSBndCond::m_spacedim, and Vmath::Zero().

Referenced by v_Apply().

104 {
105  int nq = x.num_elements();
106 
107  // Flow parameters
108  const NekDouble x0 = 5.0;
109  const NekDouble y0 = 0.0;
110  const NekDouble beta = 5.0;
111  const NekDouble u0 = 1.0;
112  const NekDouble v0 = 0.5;
113  const NekDouble gamma = m_gamma;
114  NekDouble r, xbar, ybar, tmp;
115  NekDouble fac = 1.0/(16.0*gamma*M_PI*M_PI);
116 
117  // In 3D zero rhow field.
118  if (m_spacedim == 3)
119  {
120  Vmath::Zero(nq, &u[3][o], 1);
121  }
122 
123  // Fill storage
124  for (int i = 0; i < nq; ++i)
125  {
126  xbar = x[i] - u0*time - x0;
127  ybar = y[i] - v0*time - y0;
128  r = sqrt(xbar*xbar + ybar*ybar);
129  tmp = beta*exp(1-r*r);
130  u[0][i+o] = pow(1.0 - (gamma-1.0)*tmp*tmp*fac, 1.0/(gamma-1.0));
131  u[1][i+o] = u[0][i+o]*(u0 - tmp*ybar/(2*M_PI));
132  u[2][i+o] = u[0][i+o]*(v0 + tmp*xbar/(2*M_PI));
133  u[m_spacedim+1][i+o] = pow(u[0][i+o], gamma)/(gamma-1.0) +
134  0.5*(u[1][i+o]*u[1][i+o] + u[2][i+o]*u[2][i+o]) / u[0][i+o];
135  }
136 }
int m_spacedim
Space dimension.
Definition: CFSBndCond.h:90
double NekDouble
NekDouble m_gamma
Parameters of the flow.
Definition: CFSBndCond.h:95
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:373
void Nektar::IsentropicVortexBC::v_Apply ( Array< OneD, Array< OneD, NekDouble > > &  Fwd,
Array< OneD, Array< OneD, NekDouble > > &  physarray,
const NekDouble time 
)
protectedvirtual

Implements Nektar::CFSBndCond.

Definition at line 58 of file IsentropicVortexBC.cpp.

References EvaluateIsentropicVortex(), Nektar::CFSBndCond::m_bcRegion, Nektar::CFSBndCond::m_fields, Nektar::CFSBndCond::m_offset, and Vmath::Vcopy().

62 {
63  int nvariables = physarray.num_elements();
64 
65  const Array<OneD, const int> &bndTraceMap = m_fields[0]->GetTraceBndMap();
66  // loop over Boundary Regions
67  int npoints, id1, id2, e_max;
68 
69  e_max = m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExpSize();
70 
71  for(int e = 0; e < e_max; ++e)
72  {
73  npoints = m_fields[0]->
74  GetBndCondExpansions()[m_bcRegion]->GetExp(e)->GetTotPoints();
75  id1 = m_fields[0]->
76  GetBndCondExpansions()[m_bcRegion]->GetPhys_Offset(e);
77  id2 = m_fields[0]->GetTrace()->GetPhys_Offset(bndTraceMap[m_offset+e]);
78 
79  Array<OneD,NekDouble> x(npoints, 0.0);
80  Array<OneD,NekDouble> y(npoints, 0.0);
81  Array<OneD,NekDouble> z(npoints, 0.0);
82 
83  m_fields[0]->GetBndCondExpansions()[m_bcRegion]->
84  GetExp(e)->GetCoords(x, y, z);
85 
86  EvaluateIsentropicVortex(x, y, z, Fwd, time, id2);
87 
88  for (int i = 0; i < nvariables; ++i)
89  {
90  Vmath::Vcopy(npoints, &Fwd[i][id2], 1,
91  &(m_fields[i]->GetBndCondExpansions()[m_bcRegion]->
92  UpdatePhys())[id1], 1);
93  }
94  }
95 }
void EvaluateIsentropicVortex(const Array< OneD, NekDouble > &x, const Array< OneD, NekDouble > &y, const Array< OneD, NekDouble > &z, Array< OneD, Array< OneD, NekDouble > > &u, NekDouble time, const int o=0)
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array of fields.
Definition: CFSBndCond.h:86
int m_offset
Offset.
Definition: CFSBndCond.h:103
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1061
int m_bcRegion
Id of the boundary region.
Definition: CFSBndCond.h:101

Friends And Related Function Documentation

friend class MemoryManager< IsentropicVortexBC >
friend

Definition at line 52 of file IsentropicVortexBC.h.

Member Data Documentation

std::string Nektar::IsentropicVortexBC::className
static
Initial value:
RegisterCreatorFunction("IsentropicVortex",
"Isentropic vortex boundary condition.")

Name of the class.

Definition at line 68 of file IsentropicVortexBC.h.