47 pSession->LoadParameter(
"Mco",
m_Mco, 0.01);
70 double rhowL,
double EL,
double rhoR,
71 double rhouR,
double rhovR,
double rhowR,
72 double ER,
double &rhof,
double &rhouf,
73 double &rhovf,
double &rhowf,
double &Ef)
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;
102 NekDouble Mtilde = 0.5 * (ML * ML + MR * MR);
111 NekDouble Mp = -(Kp / fa) * ((pR - pL) / (rhoA * cA * cA)) *
112 std::max(1.0 - sigma * Mtilde, 0.0);
116 NekDouble pu = -2.0 * Ku * rhoA * cA * cA * (MR - ML) *
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);
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);
double P5Function(int A, double alpha, double M)
double M4Function(int A, double beta, double M)
static RiemannSolverSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession)
AUSM3Solver(const LibUtilities::SessionReaderSharedPtr &pSession)
static std::string solverName
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.
EquationOfStateSharedPtr m_eos
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
std::shared_ptr< SessionReader > SessionReaderSharedPtr
@ beta
Gauss Radau pinned at x=-1,.
RiemannSolverFactory & GetRiemannSolverFactory()
scalarT< T > sqrt(scalarT< T > in)