45 "Stagnation conditions inflow boundary condition.");
51 const int pSpaceDim,
const int bcRegion,
const int cnt)
52 :
CFSBndCond(pSession, pFields, pTraceNormals, pSpaceDim, bcRegion, cnt)
54 const size_t nvariables =
m_fields.size();
55 const int expdim =
m_fields[0]->GetGraph()->GetMeshDimension();
56 const int spacedim =
m_fields[0]->GetGraph()->GetSpaceDimension();
62 for (
int i = 0; i < nvariables; ++i)
71 for (
int j = 0; j < numBCPts; ++j)
75 for (
int i = 0; i < spacedim; ++i)
84 for (
int i = 0; i < spacedim; ++i)
91 for (
int i = 0; i < expdim; ++i)
95 for (
int i = expdim; i < spacedim; ++i)
107 const size_t nTracePts = Fwd[0].size();
108 const size_t nVariables = physarray.size();
111 "Number of trace points does not match in "
112 "StagnationInflowBC::v_Apply()");
118 m_varConv->GetAbsoluteVelocity(Fwd, absVel);
135 for (
int i = 0; i < npts; ++i)
145 NekDouble hStat = hStag - 0.5 * pow(absVel[id2 + i], 2);
154 ->UpdatePhys())[id1 + i] = rho;
161 ->UpdatePhys())[id1 + i] =
168 ->UpdatePhys())[id1 + i] =
169 rho * (hStat /
m_gamma + 0.5 * pow(absVel[id2 + i], 2));
#define ASSERTL0(condition, msg)
Encapsulates the user-defined boundary conditions for compressible flow solver.
int m_spacedim
Space dimension.
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
Trace normals.
NekDouble m_gamma
Parameters of the flow.
int m_bcRegion
Id of the boundary region.
VariableConverterSharedPtr m_varConv
Auxiliary object to convert variables.
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array of fields.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
StagnationInflowBC(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 v_Apply(Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time) override
static std::string className
Name of the class.
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.
Array< OneD, Array< OneD, NekDouble > > m_fieldStorage
std::shared_ptr< SessionReader > SessionReaderSharedPtr
CFSBndCondFactory & GetCFSBndCondFactory()
Declaration of the boundary condition factory singleton.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
scalarT< T > sqrt(scalarT< T > in)