82 NekDouble eL = (EL - 0.5 * (rhouL * uL + rhovL * vL + rhowL * wL)) / rhoL;
83 NekDouble eR = (ER - 0.5 * (rhouR * uR + rhovR * vR + rhowR * wR)) / rhoR;
101 NekDouble uRoe = (srL * uL + srR * uR) / srLR;
102 NekDouble vRoe = (srL * vL + srR * vR) / srLR;
103 NekDouble wRoe = (srL * wL + srR * wR) / srLR;
104 NekDouble URoe2 = uRoe * uRoe + vRoe * vRoe + wRoe * wRoe;
105 NekDouble HRoe = (srL * HL + srR * HR) / srLR;
107 srR, HRoe, URoe2, srLR);
110 NekDouble SL = std::min(uL - cL, uRoe - cRoe);
111 NekDouble SR = std::max(uR + cR, uRoe + cRoe);
117 rhouf = rhouL * uL + pL;
126 rhouf = rhouR * uR + pR;
134 NekDouble SM = (pR - pL + rhouL * (SL - uL) - rhouR * (SR - uR)) /
135 (rhoL * (SL - uL) - rhoR * (SR - uR));
136 NekDouble rhoML = rhoL * (SL - uL) / (SL - SM);
141 rhoML * (EL / rhoL + (SM - uL) * (SM + pL / (rhoL * (SL - uL))));
143 NekDouble rhoMR = rhoR * (SR - uR) / (SR - SM);
148 rhoMR * (ER / rhoR + (SM - uR) * (SM + pR / (rhoR * (SR - uR))));
150 if (SL < 0.0 && SM >= 0.0)
152 rhof = rhouL + SL * (rhoML - rhoL);
153 rhouf = rhouL * uL + pL + SL * (rhouML - rhouL);
154 rhovf = rhouL * vL + SL * (rhovML - rhovL);
155 rhowf = rhouL * wL + SL * (rhowML - rhowL);
156 Ef = uL * (EL + pL) + SL * (EML - EL);
158 else if (SM < 0.0 && SR > 0.0)
160 rhof = rhouR + SR * (rhoMR - rhoR);
161 rhouf = rhouR * uR + pR + SR * (rhouMR - rhouR);
162 rhovf = rhouR * vR + SR * (rhovMR - rhovR);
163 rhowf = rhouR * wR + SR * (rhowMR - rhowR);
164 Ef = uR * (ER + pR) + SR * (EMR - ER);
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)
EquationOfStateSharedPtr m_eos
void v_PointSolve(NekDouble rhoL, NekDouble rhouL, NekDouble rhovL, NekDouble rhowL, NekDouble EL, NekDouble rhoR, NekDouble rhouR, NekDouble rhovR, NekDouble rhowR, NekDouble ER, NekDouble &rhof, NekDouble &rhouf, NekDouble &rhovf, NekDouble &rhowf, NekDouble &Ef) override
HLLC Riemann solver.
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)