46 "Unsteady Advection equation.");
67 m_session->MatchSolverInfo(
"GJPStabilisation",
"False",
77 std::vector<std::string> vel;
106 m_session->LoadSolverInfo(
"AdvectionType", advName,
"WeakDG");
119 m_session->LoadSolverInfo(
"UpwindType", riemName,
"Upwind");
134 m_session->LoadSolverInfo(
"AdvectionType", advName,
152 ASSERTL0(
false,
"Unsupported projection type.");
170 ASSERTL0(
false,
"Implicit unsteady Advection not set up.");
195 for (
int i = 0; i < velfield.size(); ++i)
197 m_fields[0]->ExtractTracePhys(velfield[i], tmp);
218 int nVariables = inarray.size();
233 for (
int i = 0; i < nVariables; ++i)
242 x->Apply(
m_fields, inarray, outarray, time);
258 int nVariables = inarray.size();
270 if (inarray != outarray)
274 for (
int i = 0; i < nVariables; ++i)
285 int ncoeffs =
m_fields[0]->GetNcoeffs();
294 for (
int i = 0; i < nVariables; ++i)
297 std::dynamic_pointer_cast<MultiRegions::ContField>(
300 m_fields[i]->IProductWRTBase(inarray[i], wsp);
303 cfield->GetGJPForcing();
308 if (GJPData->IsSemiImplicit())
322 mtype, cfield->GetLocalToGlobalMap(),
factors);
326 m_fields[i]->BwdTrans(coeffs, outarray[i]);
331 for (
int i = 0; i < nVariables; ++i)
333 m_fields[i]->FwdTrans(inarray[i], coeffs);
334 m_fields[i]->BwdTrans(coeffs, outarray[i]);
340 ASSERTL0(
false,
"Unknown projection scheme");
356 "Dimension of flux array and velocity array do not match");
358 const int nq =
m_fields[0]->GetNpoints();
360 for (
int i = 0; i < flux.size(); ++i)
362 for (
int j = 0; j < flux[0].size(); ++j)
381 "Dimension of flux array and velocity array do not match");
383 int nq = physfield[0].size();
384 int nVariables = physfield.size();
392 nq =
m_fields[0]->Get1DScaledTotPoints(OneDptscale);
400 for (
int i = 0; i < nVariables; ++i)
409 m_fields[0]->PhysInterp1DScaled(OneDptscale, physfield[i],
423 for (
int i = 0; i < flux.size(); ++i)
425 for (
int j = 0; j < flux[0].size(); ++j)
427 Vmath::Vmul(nq, physfieldInterp[i], 1, velocityInterp[j], 1,
428 fluxInterp[i][j], 1);
433 for (
int i = 0; i < nVariables; ++i)
437 m_fields[0]->PhysGalerkinProjection1DScaled(
438 OneDptscale, fluxInterp[i][j], flux[i][j]);
449 s,
"GJP Stab. Impl. ",
450 m_session->GetSolverInfo(
"GJPStabilisation"));
#define ASSERTL0(condition, msg)
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode....
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
tBaseSharedPtr CreateInstance(tKey idKey, tParam... args)
Create an instance of the class referred to by idKey.
void DefineProjection(FuncPointerT func, ObjectPointerT obj)
void DefineOdeRhs(FuncPointerT func, ObjectPointerT obj)
void AccumulateRegion(std::string, int iolevel=0)
Accumulate elapsed time for a region.
Describe a linear system.
A base class for PDEs which include an advection component.
SolverUtils::AdvectionSharedPtr m_advObject
Advection term.
SOLVER_UTILS_EXPORT void v_InitObject(bool DeclareField=true) override
Initialisation object for EquationSystem.
int m_spacedim
Spatial dimension (>= expansion dim).
SOLVER_UTILS_EXPORT int GetTraceNpoints()
int m_expdim
Expansion dimension.
NekDouble m_timestep
Time step size.
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array holding all dependent variables.
bool m_specHP_dealiasing
Flag to determine if dealisising is usde for the Spectral/hp element discretisation.
LibUtilities::SessionReaderSharedPtr m_session
The session reader.
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
Array holding trace normals for DG simulations in the forwards direction.
SOLVER_UTILS_EXPORT int GetNpoints()
enum MultiRegions::ProjectionType m_projectionType
Type of projection; e.g continuous or discontinuous.
SOLVER_UTILS_EXPORT void SetBoundaryConditions(NekDouble time)
Evaluates the boundary conditions at the given time.
SOLVER_UTILS_EXPORT SessionFunctionSharedPtr GetFunction(std::string name, const MultiRegions::ExpListSharedPtr &field=MultiRegions::NullExpListSharedPtr, bool cache=false)
Get a SessionFunction by name.
static SOLVER_UTILS_EXPORT std::vector< ForcingSharedPtr > Load(const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields=0)
Base class for unsteady solvers.
LibUtilities::TimeIntegrationSchemeOperators m_ode
The time integration scheme operators to use.
bool m_explicitAdvection
Indicates if explicit or implicit treatment of advection is used.
SOLVER_UTILS_EXPORT void v_GenerateSummary(SummaryList &s) override
Print a summary of time stepping parameters.
bool m_homoInitialFwd
Flag to determine if simulation should start in homogeneous forward transformed state.
std::vector< SolverUtils::ForcingSharedPtr > m_forcing
Forcing terms.
Array< OneD, NekDouble > m_traceVn
Array< OneD, NekDouble > & GetNormalVel(const Array< OneD, const Array< OneD, NekDouble > > &velfield)
Get the normal velocity based on input velfield.
Array< OneD, NekDouble > & GetNormalVelocity()
Get the normal velocity.
void v_GenerateSummary(SolverUtils::SummaryList &s) override
Print Summary.
Array< OneD, Array< OneD, NekDouble > > m_velocity
Advection velocity.
static SolverUtils::EquationSystemSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
Creates an instance of this class.
void DoOdeRhs(const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
Compute the RHS.
void DoOdeProjection(const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
Compute the projection.
void GetFluxVectorDeAlias(const Array< OneD, Array< OneD, NekDouble > > &physfield, Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &flux)
Evaluate the flux at each solution point using dealiasing.
UnsteadyAdvection(const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
Session reader.
void v_InitObject(bool DeclareFields=true) override
Initialise the object.
bool m_useGJPStabilisation
static std::string className
Name of class.
void GetFluxVector(const Array< OneD, Array< OneD, NekDouble > > &physfield, Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &flux)
Evaluate the flux at each solution point.
SolverUtils::RiemannSolverSharedPtr m_riemannSolver
std::shared_ptr< SessionReader > SessionReaderSharedPtr
@ eMixed_CG_Discontinuous
std::shared_ptr< GJPStabilisation > GJPStabilisationSharedPtr
std::shared_ptr< ContField > ContFieldSharedPtr
AdvectionFactory & GetAdvectionFactory()
Gets the factory for initialising advection objects.
std::vector< std::pair< std::string, std::string > > SummaryList
EquationSystemFactory & GetEquationSystemFactory()
void AddSummaryItem(SummaryList &l, const std::string &name, const std::string &value)
Adds a summary item to the summary info list.
RiemannSolverFactory & GetRiemannSolverFactory()
std::shared_ptr< MeshGraph > MeshGraphSharedPtr
std::map< ConstFactorType, NekDouble > ConstFactorMap
StdRegions::ConstFactorMap factors
static Array< OneD, NekDouble > NullNekDouble1DArray
void Vmul(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Multiply vector z = x*y.
void Neg(int n, T *x, const int incx)
Negate x = -x.
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 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)