44 : NonlinearSWESolver(pSession)
65 static NekDouble g = m_params[
"gravity"]();
78 NekDouble hC, huC, hvC, SL, SR, hstar, ustar, Sstar;
79 hstar = 0.5 * (cL + cR) + 0.25 * (uL - uR);
82 ustar = 0.5 * (uL + uR) + cL - cR;
87 SL = uL - cL *
sqrt(0.5 * ((hstar * hstar + hstar * hL) / (hL * hL)));
97 SR = uR + cR *
sqrt(0.5 * ((hstar * hstar + hstar * hR) / (hR * hR)));
104 if (fabs(hR * (uR - SR) - hL * (uL - SL)) <= 1.0e-10)
110 Sstar = (SL * hR * (uR - SR) - SR * hL * (uL - SL)) /
111 (hR * (uR - SR) - hL * (uL - SL));
117 huf = uL * uL * hL + 0.5 * g * hL * hL;
123 huf = uR * uR * hR + 0.5 * g * hR * hR;
126 else if ((SL < 0) && (Sstar >= 0))
128 hC = hL * ((SL - uL) / (SL - Sstar));
132 hf = hL * uL + SL * (hC - hL);
133 huf = (uL * uL * hL + 0.5 * g * hL * hL) + SL * (huC - hL * uL);
134 hvf = (uL * vL * hL) + SL * (hvC - hL * vL);
136 else if ((SR > 0) && (Sstar <= 0))
138 hC = hR * ((SR - uR) / (SR - Sstar));
142 hf = hR * uR + SR * (hC - hR);
143 huf = (uR * uR * hR + 0.5 * g * hR * hR) + SR * (huC - hR * uR);
144 hvf = (uR * vR * hR) + SR * (hvC - hR * vR);
149 "Error in HLLC solver -- non physical combination of "
#define NEKERROR(type, msg)
Assert Level 0 – Fundamental assert which is used whether in FULLDEBUG, DEBUG or OPT compilation mode...
static std::string solverName
static RiemannSolverSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession)
HLLCSolver(const LibUtilities::SessionReaderSharedPtr &pSession)
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
std::shared_ptr< SessionReader > SessionReaderSharedPtr
RiemannSolverFactory & GetRiemannSolverFactory()
scalarT< T > sqrt(scalarT< T > in)