Skip to content
Snippets Groups Projects
Commit 187153d7 authored by Cormet Dylan's avatar Cormet Dylan
Browse files

correct calculus of DL[Dim+1] and DR[Dim+1]

parent 5893b488
No related branches found
No related tags found
No related merge requests found
......@@ -1230,7 +1230,7 @@ class HybridHLLcRusanovEulerianCompositeSolver_v2
for (size_t d = 2; d < (Dimension + 1); ++d){
DL[d] = Uj[d-1];
}
DL[Dimension + 1] = Ej + (SC - uL) * (SC + (Pressionj / (rhoj * diffVelocityj)));
DL[Dimension + 1] = Ej + rhoj * (SC - uL) * (SC + (Pressionj / (rhoj * diffVelocityj)));
const Rp UCL = rhoj * (diffVelocityj / (MinVpNormjk - SC)) * DL;
const Rp diffStates = MinVpNormjk * l2Norm(Cjf_loc) * (UCL - StateAtFace[j][l]);
......@@ -1250,7 +1250,7 @@ class HybridHLLcRusanovEulerianCompositeSolver_v2
for (size_t d = 2; d < (Dimension + 1); ++d){
DR[d] = UK[d-1];
}
DR[Dimension + 1] = EK + (SC - uR) * (SC + (PressionK / (rhoK * diffVelocityK)));
DR[Dimension + 1] = EK + rhoK * (SC - uR) * (SC + (PressionK / (rhoK * diffVelocityK)));
const Rp UCR = rhoK * (diffVelocityK / (MaxVpNormjk - SC)) * DR;
const Rp diffStates = MaxVpNormjk * l2Norm(Cjf_loc) * (UCR - StateAtFace[K][R]);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment