#include <LEESolver.h>
Inheritance diagram for Nektar::LEESolver:
Protected Member Functions | |
| LEESolver (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 |
 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, RSScalarFuncType > | m_scalars |
| Map of scalar function types. More... | |
| std::map< std::string, RSVecFuncType > | m_vectors |
| Map of vector function types. More... | |
| std::map< std::string, RSParamFuncType > | m_params |
| Map of parameter function types. More... | |
| std::map< std::string, RSScalarFuncType > | m_auxScal |
| Map of auxiliary scalar function types. More... | |
| std::map< std::string, RSVecFuncType > | m_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 48 of file LEESolver.h.
Constructor & Destructor Documentation
◆ LEESolver()
|
protected |
Definition at line 46 of file LEESolver.cpp.
47 : AcousticSolver(pSession)
48{
50}
AcousticSolver(const LibUtilities::SessionReaderSharedPtr &pSession)
Definition: RiemannSolvers/AcousticSolver.cpp:46
bool m_requiresRotation
Indicates whether the Riemann solver requires a rotation to be applied to the velocity fields.
Definition: RiemannSolver.h:145
References Nektar::SolverUtils::RiemannSolver::m_requiresRotation.
Member Function Documentation
◆ v_Solve()
|
overrideprotectedvirtual |
Reimplemented from Nektar::AcousticSolver.
Definition at line 55 of file LEESolver.cpp.
59{
60 int nTracePts = Fwd[0].size();
61
62 Array<OneD, Array<OneD, NekDouble>> bfFwd(nDim + 3);
63 Array<OneD, Array<OneD, NekDouble>> bfBwd(nDim + 3);
64 for (int i = 0; i < nDim + 3; i++)
65 {
66 bfFwd[i] = Array<OneD, NekDouble>(nTracePts);
67 bfBwd[i] = Array<OneD, NekDouble>(nTracePts);
68 }
69
70 GetRotBasefield(bfFwd, bfBwd);
71
72 int expDim = nDim;
73 NekDouble vF, wF;
74
75 if (expDim == 1)
76 {
77 for (int i = 0; i < nTracePts; ++i)
78 {
79 v_PointSolve(Fwd[0][i], Fwd[1][i], Fwd[2][i], 0.0, 0.0, Bwd[0][i],
80 Bwd[1][i], Bwd[2][i], 0.0, 0.0, bfFwd[0][i],
81 bfFwd[1][i], bfFwd[2][i], 0.0, 0.0, bfBwd[0][i],
82 bfBwd[1][i], bfBwd[2][i], 0.0, 0.0, flux[0][i],
83 flux[1][i], flux[2][i], vF, wF);
84 }
85 }
86 else if (expDim == 2)
87 {
88 for (int i = 0; i < nTracePts; ++i)
89 {
90 v_PointSolve(Fwd[0][i], Fwd[1][i], Fwd[2][i], Fwd[3][i], 0.0,
91 Bwd[0][i], Bwd[1][i], Bwd[2][i], Bwd[3][i], 0.0,
92 bfFwd[0][i], bfFwd[1][i], bfFwd[2][i], bfFwd[3][i],
93 0.0, bfBwd[0][i], bfBwd[1][i], bfBwd[2][i],
94 bfBwd[3][i], 0.0, flux[0][i], flux[1][i], flux[2][i],
95 flux[3][i], wF);
96 }
97 }
98 else if (expDim == 3)
99 {
100 for (int i = 0; i < nTracePts; ++i)
101 {
102 v_PointSolve(Fwd[0][i], Fwd[1][i], Fwd[2][i], Fwd[3][i], Fwd[4][i],
103 Bwd[0][i], Bwd[1][i], Bwd[2][i], Bwd[3][i], Bwd[4][i],
104 bfFwd[0][i], bfFwd[1][i], bfFwd[2][i], bfFwd[3][i],
105 bfFwd[4][i], bfBwd[0][i], bfBwd[1][i], bfBwd[2][i],
106 bfBwd[3][i], bfBwd[4][i], flux[0][i], flux[1][i],
107 flux[2][i], flux[3][i], flux[4][i]);
108 }
109 }
110}
void GetRotBasefield(Array< OneD, Array< OneD, NekDouble > > &bfFwd, Array< OneD, Array< OneD, NekDouble > > &bfBwd)
Definition: RiemannSolvers/AcousticSolver.cpp:115
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
References Nektar::AcousticSolver::GetRotBasefield(), and Nektar::AcousticSolver::v_PointSolve().
