From 187153d7b3ca7b4584e483ba0b08a35df802f3a1 Mon Sep 17 00:00:00 2001
From: Dylan Cormet <cormet.dylan@gmail.com>
Date: Mon, 28 Apr 2025 16:24:06 +0200
Subject: [PATCH] correct calculus of DL[Dim+1] and DR[Dim+1]
---
src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp b/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp
index 373dc1c2..1a30528b 100644
--- a/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp
@@ -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]);
--
GitLab