From 91c3bc19b75dcbb248063ff6f92ae74c2415b4a5 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Tue, 14 Jan 2025 21:28:45 +0100 Subject: [PATCH 1/3] computing tke based on inlet values --- SU2_CFD/src/solvers/CEulerSolver.cpp | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/SU2_CFD/src/solvers/CEulerSolver.cpp b/SU2_CFD/src/solvers/CEulerSolver.cpp index 01626563d40..615039cc9e7 100644 --- a/SU2_CFD/src/solvers/CEulerSolver.cpp +++ b/SU2_CFD/src/solvers/CEulerSolver.cpp @@ -7356,7 +7356,12 @@ void CEulerSolver::BC_Supersonic_Inlet(CGeometry *geometry, CSolver **solver_con const su2double Velocity2 = GeometryToolbox::SquaredNorm(int(MAXNDIM), Velocity); su2double Energy = Pressure / (Density * Gamma_Minus_One) + 0.5 * Velocity2; - if (tkeNeeded) Energy += GetTke_Inf(); + if (tkeNeeded) { + const su2double* Turb_Properties = config->GetInlet_TurbVal(config->GetMarker_All_TagBound(val_marker)); + const su2double Intensity = Turb_Properties[0]; + const su2double Tke = 3.0 / 2.0 * (Velocity2* pow(Intensity, 2)); + Energy += Tke; + } /*--- Primitive variables, using the derived quantities. ---*/ From f837bd5fc0103619de1d810f7b689a62f28141e5 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 23 Jan 2025 10:19:51 +0100 Subject: [PATCH 2/3] small fix --- SU2_CFD/src/solvers/CEulerSolver.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/SU2_CFD/src/solvers/CEulerSolver.cpp b/SU2_CFD/src/solvers/CEulerSolver.cpp index 615039cc9e7..b619ed87622 100644 --- a/SU2_CFD/src/solvers/CEulerSolver.cpp +++ b/SU2_CFD/src/solvers/CEulerSolver.cpp @@ -7357,7 +7357,7 @@ void CEulerSolver::BC_Supersonic_Inlet(CGeometry *geometry, CSolver **solver_con const su2double Velocity2 = GeometryToolbox::SquaredNorm(int(MAXNDIM), Velocity); su2double Energy = Pressure / (Density * Gamma_Minus_One) + 0.5 * Velocity2; if (tkeNeeded) { - const su2double* Turb_Properties = config->GetInlet_TurbVal(config->GetMarker_All_TagBound(val_marker)); + const su2double* Turb_Properties = config->GetInlet_TurbVal(Marker_Tag); const su2double Intensity = Turb_Properties[0]; const su2double Tke = 3.0 / 2.0 * (Velocity2* pow(Intensity, 2)); Energy += Tke; From ad8a784c6fd2b9bd1f5605692a9e034dcfaf4102 Mon Sep 17 00:00:00 2001 From: Cristopher Morales <98025159+Cristopher-Morales@users.noreply.github.com> Date: Wed, 5 Feb 2025 16:00:07 +0100 Subject: [PATCH 3/3] Update SU2_CFD/src/solvers/CEulerSolver.cpp Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com> --- SU2_CFD/src/solvers/CEulerSolver.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/SU2_CFD/src/solvers/CEulerSolver.cpp b/SU2_CFD/src/solvers/CEulerSolver.cpp index b619ed87622..b0e1381fef0 100644 --- a/SU2_CFD/src/solvers/CEulerSolver.cpp +++ b/SU2_CFD/src/solvers/CEulerSolver.cpp @@ -7359,7 +7359,7 @@ void CEulerSolver::BC_Supersonic_Inlet(CGeometry *geometry, CSolver **solver_con if (tkeNeeded) { const su2double* Turb_Properties = config->GetInlet_TurbVal(Marker_Tag); const su2double Intensity = Turb_Properties[0]; - const su2double Tke = 3.0 / 2.0 * (Velocity2* pow(Intensity, 2)); + const su2double Tke = 3.0 / 2.0 * (Velocity2 * pow(Intensity, 2)); Energy += Tke; }