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

#include <AUSM3Solver.h>

Inheritance diagram for Nektar::AUSM3Solver:
[legend]

Static Public Member Functions

static RiemannSolverSharedPtr create (const LibUtilities::SessionReaderSharedPtr &pSession)
 
- Static Public Member Functions inherited from Nektar::AUSM0Solver
static RiemannSolverSharedPtr create (const LibUtilities::SessionReaderSharedPtr &pSession)
 

Static Public Attributes

static std::string solverName
 
- Static Public Attributes inherited from Nektar::AUSM0Solver
static std::string solverName
 

Protected Member Functions

 AUSM3Solver (const LibUtilities::SessionReaderSharedPtr &pSession)
 
void v_PointSolve (double rhoL, double rhouL, double rhovL, double rhowL, double EL, double rhoR, double rhouR, double rhovR, double rhowR, double ER, double &rhof, double &rhouf, double &rhovf, double &rhowf, double &Ef) override
 AUSM3 Riemann solver. More...
 
- Protected Member Functions inherited from Nektar::AUSM0Solver
 AUSM0Solver (const LibUtilities::SessionReaderSharedPtr &pSession)
 
void v_PointSolve (double rhoL, double rhouL, double rhovL, double rhowL, double EL, double rhoR, double rhouR, double rhovR, double rhowR, double ER, double &rhof, double &rhouf, double &rhovf, double &rhowf, double &Ef) override
 AUSM0 Riemann solver. More...
 
double M1Function (int A, double M)
 
double M2Function (int A, double M)
 
double M4Function (int A, double beta, double M)
 
double P5Function (int A, double alpha, double M)
 
- Protected Member Functions inherited from Nektar::CompressibleSolver
 CompressibleSolver (const LibUtilities::SessionReaderSharedPtr &pSession)
 Session ctor. More...
 
 CompressibleSolver ()
 Programmatic ctor. More...
 
void v_Solve (const int nDim, const Array< OneD, const Array< OneD, ND > > &Fwd, const Array< OneD, const Array< OneD, ND > > &Bwd, Array< OneD, Array< OneD, ND > > &flux) override
 
virtual void v_ArraySolve (const Array< OneD, const Array< OneD, ND > > &Fwd, const Array< OneD, const Array< OneD, ND > > &Bwd, Array< OneD, Array< OneD, ND > > &flux)
 
virtual void v_PointSolve (ND rhoL, ND rhouL, ND rhovL, ND rhowL, ND EL, ND rhoR, ND rhouR, ND rhovR, ND rhowR, ND ER, ND &rhof, ND &rhouf, ND &rhovf, ND &rhowf, ND &Ef)
 
ND GetRoeSoundSpeed (ND rhoL, ND pL, ND eL, ND HL, ND srL, ND rhoR, ND pR, ND eR, ND HR, ND srR, ND HRoe, ND URoe2, ND srLR)
 
- 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)
 

Protected Attributes

NekDouble m_Mco
 
- Protected Attributes inherited from Nektar::CompressibleSolver
bool m_pointSolve
 
EquationOfStateSharedPtr m_eos
 
bool m_idealGas
 
- 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...
 

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 Types inherited from Nektar::CompressibleSolver
using ND = NekDouble
 

Detailed Description

Definition at line 43 of file AUSM3Solver.h.

Constructor & Destructor Documentation

◆ AUSM3Solver()

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

Definition at line 44 of file AUSM3Solver.cpp.

45 : AUSM0Solver(pSession)
46{
47 pSession->LoadParameter("Mco", m_Mco, 0.01);
48}
AUSM0Solver(const LibUtilities::SessionReaderSharedPtr &pSession)
Definition: AUSM0Solver.cpp:44

References m_Mco.

Referenced by create().

Member Function Documentation

◆ create()

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

Definition at line 46 of file AUSM3Solver.h.

48 {
49 return RiemannSolverSharedPtr(new AUSM3Solver(pSession));
50 }
AUSM3Solver(const LibUtilities::SessionReaderSharedPtr &pSession)
Definition: AUSM3Solver.cpp:44
std::shared_ptr< RiemannSolver > RiemannSolverSharedPtr
A shared pointer to an EquationSystem object.

References AUSM3Solver().

◆ v_PointSolve()

void Nektar::AUSM3Solver::v_PointSolve ( double  rhoL,
double  rhouL,
double  rhovL,
double  rhowL,
double  EL,
double  rhoR,
double  rhouR,
double  rhovR,
double  rhowR,
double  ER,
double &  rhof,
double &  rhouf,
double &  rhovf,
double &  rhowf,
double &  Ef 
)
overrideprotectedvirtual

AUSM3 Riemann solver.

Parameters
rhoLDensity left state.
rhoRDensity right state.
rhouLx-momentum component left state.
rhouRx-momentum component right state.
rhovLy-momentum component left state.
rhovRy-momentum component right state.
rhowLz-momentum component left state.
rhowRz-momentum component right state.
ELEnergy left state.
EREnergy right state.
rhofComputed Riemann flux for density.
rhoufComputed Riemann flux for x-momentum component
rhovfComputed Riemann flux for y-momentum component
rhowfComputed Riemann flux for z-momentum component
EfComputed Riemann flux for energy.

Reimplemented from Nektar::AUSM0Solver.

Definition at line 69 of file AUSM3Solver.cpp.

74{
75 // Left and Right velocities
76 NekDouble uL = rhouL / rhoL;
77 NekDouble vL = rhovL / rhoL;
78 NekDouble wL = rhowL / rhoL;
79 NekDouble uR = rhouR / rhoR;
80 NekDouble vR = rhovR / rhoR;
81 NekDouble wR = rhowR / rhoR;
82
83 // Internal energy (per unit mass)
84 NekDouble eL = (EL - 0.5 * (rhouL * uL + rhovL * vL + rhowL * wL)) / rhoL;
85 NekDouble eR = (ER - 0.5 * (rhouR * uR + rhovR * vR + rhowR * wR)) / rhoR;
86 // Pressure
87 NekDouble pL = m_eos->GetPressure(rhoL, eL);
88 NekDouble pR = m_eos->GetPressure(rhoR, eR);
89 // Speed of sound
90 NekDouble cL = m_eos->GetSoundSpeed(rhoL, eL);
91 NekDouble cR = m_eos->GetSoundSpeed(rhoR, eR);
92
93 // Average speeds of sound
94 NekDouble cA = 0.5 * (cL + cR);
95
96 // Local Mach numbers
97 NekDouble ML = uL / cA;
98 NekDouble MR = uR / cA;
99
100 // Parameters for specify the upwinding
101 // Note: if fa = 1 then AUSM3 = AUSM3
102 NekDouble Mtilde = 0.5 * (ML * ML + MR * MR);
103 NekDouble Mo = std::sqrt(std::min(1.0, std::max(Mtilde, m_Mco * m_Mco)));
104 NekDouble fa = Mo * (2.0 - Mo);
105 NekDouble beta = 0.125;
106 NekDouble alpha = 0.1875;
107 NekDouble sigma = 1.0;
108 NekDouble Kp = 0.25;
109 NekDouble Ku = 0.75;
110 NekDouble rhoA = 0.5 * (rhoL + rhoR);
111 NekDouble Mp = -(Kp / fa) * ((pR - pL) / (rhoA * cA * cA)) *
112 std::max(1.0 - sigma * Mtilde, 0.0);
113
114 NekDouble Mbar = M4Function(0, beta, ML) + M4Function(1, beta, MR) + Mp;
115
116 NekDouble pu = -2.0 * Ku * rhoA * cA * cA * (MR - ML) *
117 P5Function(0, alpha, ML) * P5Function(1, alpha, MR);
118
119 NekDouble pbar =
120 pL * P5Function(0, alpha, ML) + pR * P5Function(1, alpha, MR) + pu;
121
122 if (Mbar >= 0.0)
123 {
124 rhof = cA * Mbar * rhoL;
125 rhouf = cA * Mbar * rhoL * uL + pbar;
126 rhovf = cA * Mbar * rhoL * vL;
127 rhowf = cA * Mbar * rhoL * wL;
128 Ef = cA * Mbar * (EL + pL);
129 }
130 else
131 {
132 rhof = cA * Mbar * rhoR;
133 rhouf = cA * Mbar * rhoR * uR + pbar;
134 rhovf = cA * Mbar * rhoR * vR;
135 rhowf = cA * Mbar * rhoR * wR;
136 Ef = cA * Mbar * (ER + pR);
137 }
138}
double P5Function(int A, double alpha, double M)
double M4Function(int A, double beta, double M)
EquationOfStateSharedPtr m_eos
@ beta
Gauss Radau pinned at x=-1,.
Definition: PointsType.h:59
double NekDouble
scalarT< T > sqrt(scalarT< T > in)
Definition: scalar.hpp:285

References Nektar::LibUtilities::beta, Nektar::AUSM0Solver::M4Function(), Nektar::CompressibleSolver::m_eos, m_Mco, Nektar::AUSM0Solver::P5Function(), and tinysimd::sqrt().

Member Data Documentation

◆ m_Mco

NekDouble Nektar::AUSM3Solver::m_Mco
protected

Definition at line 62 of file AUSM3Solver.h.

Referenced by AUSM3Solver(), and v_PointSolve().

◆ solverName

std::string Nektar::AUSM3Solver::solverName
static
Initial value:
=
"AUSM3", AUSM3Solver::create, "AUSM3 Riemann solver")
static RiemannSolverSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession)
Definition: AUSM3Solver.h:46
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
RiemannSolverFactory & GetRiemannSolverFactory()

Definition at line 52 of file AUSM3Solver.h.