44 RegisterCreatorFunction(
"PressureInflowFile",
45 PressureInflowFileBC::create,
46 "Pressure inflow (file) boundary condition.");
54 :
CFSBndCond(pSession, pFields, pTraceNormals, pSpaceDim, bcRegion, cnt)
56 int nvariables =
m_fields.num_elements();
61 GetBndCondExpansions()[
m_bcRegion]->GetNpoints();
62 for (
int i = 0; i < nvariables; ++i)
78 int nTracePts =
m_fields[0]->GetTrace()->GetNpoints();
79 int nVariables = physarray.num_elements();
86 NekDouble gammaMinusOneInv = 1.0 / gammaMinusOne;
92 for (i = 0; i < nDimensions; ++i)
94 Vmath::Vdiv(nTracePts, Fwd[i+1], 1, Fwd[0], 1, Vel, 1);
101 m_varConv->GetAbsoluteVelocity(Fwd, absVel);
108 m_varConv->GetSoundSpeed(Fwd, pressure, soundSpeed);
112 Vmath::Vdiv(nTracePts, Vn, 1, soundSpeed, 1, Mach, 1);
116 int e, id1, id2,
npts, pnt;
124 GetExp(e)->GetTotPoints();
130 for (i = 0; i <
npts; i++)
135 if (Mach[pnt] < 0.99)
138 for (j = 0; j < nVariables-1; ++j)
146 for (j = 1; j < nVariables-1; ++j)
152 rhoeb = gammaMinusOneInv * pressure[pnt] + Ek;
155 UpdatePhys())[id1+i] = rhoeb;
160 for (j = 0; j < nVariables; ++j)
164 UpdatePhys())[id1+i] = Fwd[j][pnt];
int m_spacedim
Space dimension.
virtual void v_Apply(Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time)
void Vvtvp(int n, const T *w, const int incw, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
vvtvp (vector times vector plus vector): z = w*x + y
void Vdiv(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Multiply vector z = x/y.
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array of fields.
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
void Vabs(int n, const T *x, const int incx, T *y, const int incy)
vabs: y = |x|
CFSBndCondFactory & GetCFSBndCondFactory()
Declaration of the boundary condition factory singleton.
Encapsulates the user-defined boundary conditions for compressible flow solver.
NekDouble m_gamma
Parameters of the flow.
VariableConverterSharedPtr m_varConv
Auxiliary object to convert variables.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
Trace normals.
int m_bcRegion
Id of the boundary region.
Array< OneD, Array< OneD, NekDouble > > m_fieldStorage