From d7bbcc3e5a8c2f0174ad0a29274925034ccab982 Mon Sep 17 00:00:00 2001
From: HOCH PHILIPPE <philippe.hoch@gmail.com>
Date: Tue, 1 Apr 2025 18:57:51 +0200
Subject: [PATCH] Correction for 3D Wall BD CL
---
src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp | 2 ++
src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp | 3 +++
src/scheme/RusanovEulerianCompositeSolver_o2.cpp | 2 ++
src/scheme/RusanovEulerianCompositeSolver_v2.cpp | 2 ++
src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp | 2 ++
5 files changed, 11 insertions(+)
diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
index fd856ffbc..8752bbecf 100644
--- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
@@ -660,6 +660,8 @@ class RoeViscousFormEulerianCompositeSolver_v2
continue;
normal *= 1. / nbnormal;
+ normal *= 1. / l2Norm(normal);
+
const auto& edge_cell_list = edge_to_cell_matrix[edge_id];
const auto& edge_local_number_in_its_cells = edge_local_numbers_in_their_cells.itemArray(edge_id);
diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp
index 934c39b96..b16466349 100644
--- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp
+++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp
@@ -664,6 +664,9 @@ class RoeViscousFormEulerianCompositeSolver_v2_o2
if (nbnormal == 0)
continue;
normal *= 1. / nbnormal;
+
+ normal *= 1. / l2Norm(normal);
+
const auto& edge_cell_list = edge_to_cell_matrix[edge_id];
const auto& edge_local_number_in_its_cells = edge_local_numbers_in_their_cells.itemArray(edge_id);
diff --git a/src/scheme/RusanovEulerianCompositeSolver_o2.cpp b/src/scheme/RusanovEulerianCompositeSolver_o2.cpp
index 9e59e875d..e17e9dc71 100644
--- a/src/scheme/RusanovEulerianCompositeSolver_o2.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver_o2.cpp
@@ -333,6 +333,8 @@ class RusanovEulerianCompositeSolver_o2
continue;
normal *= 1. / nbnormal;
+ normal *= 1. / l2Norm(normal);
+
const auto& edge_cell_list = edge_to_cell_matrix[edge_id];
const auto& edge_local_number_in_its_cells = edge_local_numbers_in_their_cells.itemArray(edge_id);
diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
index 52269034a..a68b66b90 100644
--- a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
@@ -323,6 +323,8 @@ class RusanovEulerianCompositeSolver_v2
continue;
normal *= 1. / nbnormal;
+ normal *= 1. / l2Norm(normal);
+
const auto& edge_cell_list = edge_to_cell_matrix[edge_id];
const auto& edge_local_number_in_its_cells = edge_local_numbers_in_their_cells.itemArray(edge_id);
diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp
index 399c0c230..e1c4f46ab 100644
--- a/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp
@@ -331,6 +331,8 @@ class RusanovEulerianCompositeSolver_v2_o2
continue;
normal *= 1. / nbnormal;
+ normal *= 1. / l2Norm(normal);
+
const auto& edge_cell_list = edge_to_cell_matrix[edge_id];
const auto& edge_local_number_in_its_cells = edge_local_numbers_in_their_cells.itemArray(edge_id);
--
GitLab