45 "No-slip (viscous) wall boundary condition.");
50 "Adiabatic wall boundary condition.");
56 const int pSpaceDim,
const int bcRegion,
const int cnt)
57 :
CFSBndCond(pSession, pFields, pTraceNormals, pSpaceDim, bcRegion, cnt)
69 int nVariables = physarray.size();
79 for (i = 0; i < nVariables; ++i)
84 m_fields[i]->EvaluateBoundaryConditions(time, varName);
96 int e, id1, id2, nBCEdgePts, eMax;
100 for (e = 0; e < eMax; ++e)
114 Vmath::Zero(nBCEdgePts, &Fwd[nVariables - 1][id2], 1);
124 for (i = 0; i < nVariables; ++i)
134 for (i = 0; i < nVariables; ++i)
139 ->UpdatePhys())[id1],
Encapsulates the user-defined boundary conditions for compressible flow solver.
NekDouble m_diffusionAveWeight
Weight for average calculation of diffusion term.
LibUtilities::SessionReaderSharedPtr m_session
Session reader.
int m_spacedim
Space dimension.
int m_bcRegion
Id of the boundary region.
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array of fields.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
static std::string classNameAdiabatic
WallViscousBC(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)
Array< OneD, Array< OneD, NekDouble > > m_bndPhys
void v_Apply(Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time) override
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.
static std::string classNameViscous
Name of the class.
std::shared_ptr< SessionReader > SessionReaderSharedPtr
CFSBndCondFactory & GetCFSBndCondFactory()
Declaration of the boundary condition factory singleton.
void Neg(int n, T *x, const int incx)
Negate x = -x.
void Vadd(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Add vector z = x+y.
void Zero(int n, T *x, const int incx)
Zero vector.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)