35 #ifndef NEKTAR_SOLVERUTILS_ADVECTION
36 #define NEKTAR_SOLVERUTILS_ADVECTION
59 typedef std::function<
void(
const Array<OneD, Array<OneD, NekDouble>> &,
60 Array<OneD, Array<OneD, Array<OneD, NekDouble>>> &)>
83 const int nConvectiveFields,
95 const int nConvectiveFields,
101 v_AdvectVolumeFlux(nConvectiveFields, pFields, pAdvVel, pInarray,
107 const int nConvectiveFields,
117 v_AdvectTraceFlux(nConvectiveFields, pFields, pAdvVel, pInarray,
118 pTraceFlux, pTime, pFwd, pBwd);
128 const int nConvectiveFields,
145 template <
typename FuncPo
interT,
typename ObjectPo
interT>
149 std::bind(
func, obj, std::placeholders::_1, std::placeholders::_2);
169 m_fluxVector = fluxVector;
181 v_SetBaseFlow(inarray, fields);
184 template <
typename DataType,
typename TypeNekBlkMatSharedPtr>
186 const int nConvectiveFields,
const int nSpaceDim,
193 template <
typename DataType,
typename TypeNekBlkMatSharedPtr>
194 void CalcJacobTraceInteg(
203 const int &nConvectiveFields,
207 v_AddVolumJacToMat(pFields, nConvectiveFields, ElmtJacArray, gmtxarray);
226 const int nConvectiveFields,
238 const int nConvectiveFields,
246 const int nConvectiveFields,
257 const int nConvectiveFields,
274 const int &nConvectiveFields,
#define SOLVER_UTILS_EXPORT
Provides a generic Factory class.
An abstract base class encapsulating the concept of advection of a vector field.
void SetBaseFlow(const Array< OneD, Array< OneD, NekDouble >> &inarray, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
Set the base flow used for linearised advection objects.
int m_spaceDim
Storage for space dimension. Used for homogeneous extension.
virtual SOLVER_UTILS_EXPORT ~Advection()
AdvectionFluxVecCB m_fluxVector
Callback function to the flux vector (set when advection is in conservative form).
virtual SOLVER_UTILS_EXPORT void v_Advect(const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &advVel, const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray)=0
Advects a vector field.
SOLVER_UTILS_EXPORT void AdvectTraceFlux(const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &pAdvVel, const Array< OneD, Array< OneD, NekDouble >> &pInarray, Array< OneD, Array< OneD, NekDouble >> &pTraceFlux, const NekDouble &pTime, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray)
Interface function to advect the Trace field.
SOLVER_UTILS_EXPORT void AddVolumJacToMat(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int &nConvectiveFields, const TensorOfArray5D< NekDouble > &ElmtJacArray, Array< OneD, Array< OneD, SNekBlkMatSharedPtr >> &gmtxarray)
SOLVER_UTILS_EXPORT void AdvectVolumeFlux(const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &pAdvVel, const Array< OneD, Array< OneD, NekDouble >> &pInarray, TensorOfArray3D< NekDouble > &pVolumeFlux, const NekDouble &pTime)
Interface function to advect the Volume field.
SOLVER_UTILS_EXPORT void AddTraceJacToMat(const int nConvectiveFields, const int nSpaceDim, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacCons, Array< OneD, Array< OneD, TypeNekBlkMatSharedPtr >> &gmtxarray, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacGrad, const Array< OneD, Array< OneD, DataType >> &TracePntJacGradSign)
RiemannSolverSharedPtr m_riemann
Riemann solver for DG-type schemes.
void SetFluxVector(AdvectionFluxVecCB fluxVector)
Set the flux vector callback function.
void SetRiemannSolver(RiemannSolverSharedPtr riemann)
Set a Riemann solver object for this advection object.
void SetFluxVector(FuncPointerT func, ObjectPointerT obj)
Set the flux vector callback function.
std::shared_ptr< SessionReader > SessionReaderSharedPtr
std::function< void(const Array< OneD, Array< OneD, NekDouble >> &, Array< OneD, Array< OneD, Array< OneD, NekDouble >>> &)> AdvectionFluxVecCB
std::shared_ptr< RiemannSolver > RiemannSolverSharedPtr
A shared pointer to an EquationSystem object.
AdvectionFactory & GetAdvectionFactory()
Gets the factory for initialising advection objects.
LibUtilities::NekFactory< std::string, Advection, std::string > AdvectionFactory
Datatype of the NekFactory used to instantiate classes derived from the Advection class.
std::shared_ptr< Advection > AdvectionSharedPtr
A shared pointer to an Advection object.
The above copyright notice and this permission notice shall be included.
static Array< OneD, Array< OneD, NekDouble > > NullNekDoubleArrayOfArray