Nektar++
Static Public Member Functions | Static Public Attributes | Protected Member Functions | List of all members
Nektar::APEUpwindSolver Class Reference

#include <APEUpwindSolver.h>

Inheritance diagram for Nektar::APEUpwindSolver:
[legend]

Static Public Member Functions

static RiemannSolverSharedPtr create (const LibUtilities::SessionReaderSharedPtr &pSession)
 

Static Public Attributes

static std::string solverName
 

Protected Member Functions

 APEUpwindSolver (const LibUtilities::SessionReaderSharedPtr &pSession)
 
void v_PointSolve (NekDouble pL, NekDouble rhoL, NekDouble uL, NekDouble vL, NekDouble wL, NekDouble pR, NekDouble rhoR, NekDouble uR, NekDouble vR, NekDouble wR, NekDouble c0sqL, NekDouble rho0L, NekDouble u0L, NekDouble v0L, NekDouble w0L, NekDouble c0sqR, NekDouble rho0R, NekDouble u0R, NekDouble v0R, NekDouble w0R, NekDouble &pF, NekDouble &rhoF, NekDouble &uF, NekDouble &vF, NekDouble &wF) override
 Upwind Riemann solver. More...
 
- Protected Member Functions inherited from Nektar::AcousticSolver
 AcousticSolver (const LibUtilities::SessionReaderSharedPtr &pSession)
 
void v_Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux) override
 
virtual void v_PointSolve (NekDouble pL, NekDouble rhoL, NekDouble uL, NekDouble vL, NekDouble wL, NekDouble pR, NekDouble rhoR, NekDouble uR, NekDouble vR, NekDouble wR, NekDouble c0sqL, NekDouble rho0L, NekDouble u0L, NekDouble v0L, NekDouble w0L, NekDouble c0sqR, NekDouble rho0R, NekDouble u0R, NekDouble v0R, NekDouble w0R, NekDouble &pF, NekDouble &rhoF, NekDouble &uF, NekDouble &vF, NekDouble &wF)=0
 
void GetRotBasefield (Array< OneD, Array< OneD, NekDouble > > &bfFwd, Array< OneD, Array< OneD, NekDouble > > &bfBwd)
 
- Protected Member Functions inherited from Nektar::SolverUtils::RiemannSolver
SOLVER_UTILS_EXPORT RiemannSolver ()
 
SOLVER_UTILS_EXPORT RiemannSolver (const LibUtilities::SessionReaderSharedPtr &pSession)
 
virtual SOLVER_UTILS_EXPORT ~RiemannSolver ()
 
virtual void v_Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)=0
 
SOLVER_UTILS_EXPORT void GenerateRotationMatrices (const Array< OneD, const Array< OneD, NekDouble > > &normals)
 Generate rotation matrices for 3D expansions. More...
 
void FromToRotation (Array< OneD, const NekDouble > &from, Array< OneD, const NekDouble > &to, NekDouble *mat)
 A function for creating a rotation matrix that rotates a vector from into another vector to. More...
 
SOLVER_UTILS_EXPORT void rotateToNormal (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &normals, const Array< OneD, const Array< OneD, NekDouble > > &vecLocs, Array< OneD, Array< OneD, NekDouble > > &outarray)
 Rotate a vector field to trace normal. More...
 
SOLVER_UTILS_EXPORT void rotateFromNormal (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &normals, const Array< OneD, const Array< OneD, NekDouble > > &vecLocs, Array< OneD, Array< OneD, NekDouble > > &outarray)
 Rotate a vector field from trace normal. More...
 
SOLVER_UTILS_EXPORT bool CheckScalars (std::string name)
 Determine whether a scalar has been defined in m_scalars. More...
 
SOLVER_UTILS_EXPORT bool CheckVectors (std::string name)
 Determine whether a vector has been defined in m_vectors. More...
 
SOLVER_UTILS_EXPORT bool CheckParams (std::string name)
 Determine whether a parameter has been defined in m_params. More...
 
SOLVER_UTILS_EXPORT bool CheckAuxScal (std::string name)
 Determine whether a scalar has been defined in m_auxScal. More...
 
SOLVER_UTILS_EXPORT bool CheckAuxVec (std::string name)
 Determine whether a vector has been defined in m_auxVec. More...
 
virtual SOLVER_UTILS_EXPORT void v_CalcFluxJacobian (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, const Array< OneD, const Array< OneD, NekDouble > > &normals, DNekBlkMatSharedPtr &FJac, DNekBlkMatSharedPtr &BJac)
 

Additional Inherited Members

- Public Member Functions inherited from Nektar::SolverUtils::RiemannSolver
SOLVER_UTILS_EXPORT void Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)
 Perform the Riemann solve given the forwards and backwards spaces. More...
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetScalar (std::string name, FuncPointerT func, ObjectPointerT obj)
 
void SetScalar (std::string name, RSScalarFuncType fp)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetVector (std::string name, FuncPointerT func, ObjectPointerT obj)
 
void SetVector (std::string name, RSVecFuncType fp)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetParam (std::string name, FuncPointerT func, ObjectPointerT obj)
 
void SetALEFlag (bool &ALE)
 
void SetParam (std::string name, RSParamFuncType fp)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetAuxScal (std::string name, FuncPointerT func, ObjectPointerT obj)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetAuxVec (std::string name, FuncPointerT func, ObjectPointerT obj)
 
void SetAuxVec (std::string name, RSVecFuncType fp)
 
std::map< std::string, RSScalarFuncType > & GetScalars ()
 
std::map< std::string, RSVecFuncType > & GetVectors ()
 
std::map< std::string, RSParamFuncType > & GetParams ()
 
SOLVER_UTILS_EXPORT void CalcFluxJacobian (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, DNekBlkMatSharedPtr &FJac, DNekBlkMatSharedPtr &BJac)
 Calculate the flux jacobian of Fwd and Bwd. More...
 
- Public Attributes inherited from Nektar::SolverUtils::RiemannSolver
int m_spacedim
 
- Protected Attributes inherited from Nektar::SolverUtils::RiemannSolver
bool m_requiresRotation
 Indicates whether the Riemann solver requires a rotation to be applied to the velocity fields. More...
 
std::map< std::string, RSScalarFuncTypem_scalars
 Map of scalar function types. More...
 
std::map< std::string, RSVecFuncTypem_vectors
 Map of vector function types. More...
 
std::map< std::string, RSParamFuncTypem_params
 Map of parameter function types. More...
 
std::map< std::string, RSScalarFuncTypem_auxScal
 Map of auxiliary scalar function types. More...
 
std::map< std::string, RSVecFuncTypem_auxVec
 Map of auxiliary vector function types. More...
 
Array< OneD, Array< OneD, NekDouble > > m_rotMat
 Rotation matrices for each trace quadrature point. More...
 
Array< OneD, Array< OneD, Array< OneD, NekDouble > > > m_rotStorage
 Rotation storage. More...
 
bool m_ALESolver = false
 Flag if using the ALE formulation. More...
 

Detailed Description

Definition at line 47 of file APEUpwindSolver.h.

Constructor & Destructor Documentation

◆ APEUpwindSolver()

Nektar::APEUpwindSolver::APEUpwindSolver ( const LibUtilities::SessionReaderSharedPtr pSession)
protected

Definition at line 46 of file APEUpwindSolver.cpp.

48 : AcousticSolver(pSession)
49{
50}
AcousticSolver(const LibUtilities::SessionReaderSharedPtr &pSession)

Referenced by create().

Member Function Documentation

◆ create()

static RiemannSolverSharedPtr Nektar::APEUpwindSolver::create ( const LibUtilities::SessionReaderSharedPtr pSession)
inlinestatic

Definition at line 50 of file APEUpwindSolver.h.

52 {
53 return RiemannSolverSharedPtr(new APEUpwindSolver(pSession));
54 }
APEUpwindSolver(const LibUtilities::SessionReaderSharedPtr &pSession)
std::shared_ptr< RiemannSolver > RiemannSolverSharedPtr
A shared pointer to an EquationSystem object.

References APEUpwindSolver().

◆ v_PointSolve()

void Nektar::APEUpwindSolver::v_PointSolve ( NekDouble  pL,
NekDouble  rhoL,
NekDouble  uL,
NekDouble  vL,
NekDouble  wL,
NekDouble  pR,
NekDouble  rhoR,
NekDouble  uR,
NekDouble  vR,
NekDouble  wR,
NekDouble  c0sqL,
NekDouble  rho0L,
NekDouble  u0L,
NekDouble  v0L,
NekDouble  w0L,
NekDouble  c0sqR,
NekDouble  rho0R,
NekDouble  u0R,
NekDouble  v0R,
NekDouble  w0R,
NekDouble pF,
NekDouble rhoF,
NekDouble uF,
NekDouble vF,
NekDouble wF 
)
overrideprotectedvirtual

Upwind Riemann solver.

The fluxes are straight out of sympy, so lets just hope the compiler optimizes them for us.

Parameters
pLPerturbation pressure left state
rhoLPerturbation density left state
pRPerturbation pressure right state
rhoRPerturbation density right state
uLx perturbation velocity component left state
uRx perturbation velocity component right state
vLy perturbation velocity component left state
vRy perturbation velocity component right state
wLz perturbation velocity component left state
wRz perturbation velocity component right state
c0sqLBase pressure left state
c0sqRBase pressure right state
rho0LBase density left state
rho0RBase density right state
u0LBase x velocity component left state
u0RBase x velocity component right state
v0LBase y velocity component left state
v0RBase y velocity component right state
w0LBase z velocity component left state
w0RBase z velocity component right state
pFComputed Riemann flux for perturbation pressure
rhoFComputed Riemann flux for perturbation density
uFComputed Riemann flux for x perturbation velocity component
vFComputed Riemann flux for y perturbation velocity component
wFComputed Riemann flux for z perturbation velocity component

Implements Nektar::AcousticSolver.

Definition at line 84 of file APEUpwindSolver.cpp.

92{
93 // Speed of sound
94 NekDouble c0L = sqrt(c0sqL);
95 NekDouble c0R = sqrt(c0sqR);
96 NekDouble c0M = (c0L + c0R) / 2.0;
97
98 NekDouble u0M = (u0L + u0R) / 2.0;
99
100 pF = 0.0;
101 uF = 0.0;
102 vF = 0.0;
103 wF = 0.0;
104
105 // lambda_3
106 if (u0M - c0M > 0)
107 {
108 pF = pF + 0.5 * (c0L - u0L) *
109 (-rho0L * v0L * vL * (c0L * u0L + c0sqL) -
110 rho0L * w0L * wL * (c0L * u0L + c0sqL) +
111 (c0sqL - pow(u0L, 2)) * (rho0L * c0L * uL - pL)) /
112 (c0sqL - pow(u0L, 2));
113
114 uF = uF +
115 0.5 * (c0L - u0L) *
116 (rho0L * c0L * (c0L * u0L + c0sqL) * (v0L * vL + w0L * wL) -
117 rho0L * c0sqL * uL * (c0sqL - pow(u0L, 2)) +
118 c0L * pL * (c0sqL - pow(u0L, 2))) /
119 (rho0L * c0sqL * (c0sqL - pow(u0L, 2)));
120 }
121 else
122 {
123 pF = pF + 0.5 * (c0R - u0R) *
124 (-rho0R * v0R * vR * (c0R * u0R + c0sqR) -
125 rho0R * w0R * wR * (c0R * u0R + c0sqR) +
126 (c0sqR - pow(u0R, 2)) * (rho0R * c0R * uR - pR)) /
127 (c0sqR - pow(u0R, 2));
128
129 uF = uF +
130 0.5 * (c0R - u0R) *
131 (rho0R * c0R * (c0R * u0R + c0sqR) * (v0R * vR + w0R * wR) -
132 rho0R * c0sqR * uR * (c0sqR - pow(u0R, 2)) +
133 c0R * pR * (c0sqR - pow(u0R, 2))) /
134 (rho0R * c0sqR * (c0sqR - pow(u0R, 2)));
135 }
136
137 // lambda_4
138 if (u0M + c0M > 0)
139 {
140 pF = pF + 0.5 * (c0L + u0L) *
141 (-rho0L * v0L * vL * (c0L * u0L - c0sqL) -
142 rho0L * w0L * wL * (c0L * u0L - c0sqL) +
143 (c0sqL - pow(u0L, 2)) * (rho0L * c0L * uL + pL)) /
144 (c0sqL - pow(u0L, 2));
145
146 uF = uF +
147 0.5 * (c0L + u0L) *
148 (-rho0L * c0L * (c0L * u0L - c0sqL) * (v0L * vL + w0L * wL) +
149 rho0L * c0sqL * uL * (c0sqL - pow(u0L, 2)) +
150 c0L * pL * (c0sqL - pow(u0L, 2))) /
151 (rho0L * c0sqL * (c0sqL - pow(u0L, 2)));
152 }
153 else
154 {
155 pF = pF + 0.5 * (c0R + u0R) *
156 (-rho0R * v0R * vR * (c0R * u0R - c0sqR) -
157 rho0R * w0R * wR * (c0R * u0R - c0sqR) +
158 (c0sqR - pow(u0R, 2)) * (rho0R * c0R * uR + pR)) /
159 (c0sqR - pow(u0R, 2));
160
161 uF = uF +
162 0.5 * (c0R + u0R) *
163 (-rho0R * c0R * (c0R * u0R - c0sqR) * (v0R * vR + w0R * wR) +
164 rho0R * c0sqR * uR * (c0sqR - pow(u0R, 2)) +
165 c0R * pR * (c0sqR - pow(u0R, 2))) /
166 (rho0R * c0sqR * (c0sqR - pow(u0R, 2)));
167 }
168}
double NekDouble
scalarT< T > sqrt(scalarT< T > in)
Definition: scalar.hpp:285

References tinysimd::sqrt().

Member Data Documentation

◆ solverName

std::string Nektar::APEUpwindSolver::solverName
static
Initial value:
=
"Upwind solver for the APE equation")
static RiemannSolverSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession)
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
RiemannSolverFactory & GetRiemannSolverFactory()

Definition at line 56 of file APEUpwindSolver.h.