Skip to content

Commit 7132d99

Browse files
Cristopher-Moralespcarruscag
Cristopher-Morales
and
pcarruscag
authored
BC inlet for SST solver based on local conditions at the Inlet (#1796)
* BC inlet based on fluid properties at the inlets * cleaning up * moving implementations to UniformInlet * addressing previous comments * fixing segmentation faults * let postprocess dealing with muT * fixing inconsistencies * small fix * updating residuals * update restart test cases * updating residuals * small fix * point testcases to develop * generic way to access primitive indices Co-authored-by: Pedro Gomes <[email protected]>
1 parent 7f96fbe commit 7132d99

File tree

11 files changed

+155
-34
lines changed

11 files changed

+155
-34
lines changed

‎Common/src/CConfig.cpp‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1222,7 +1222,7 @@ void CConfig::SetConfig_Options() {
12221222
addDoubleListOption("MU_T_REF", nMu_Temperature_Ref, Mu_Temperature_Ref);
12231223
/* DESCRIPTION: Sutherland constant, default value for AIR SI */
12241224
addDoubleListOption("SUTHERLAND_CONSTANT", nMu_S, Mu_S);
1225-
1225+
12261226
/*--- Options related to Viscosity Model ---*/
12271227
/*!\brief MIXINGVISCOSITY_MODEL \n DESCRIPTION: Mixing model of the viscosity \n OPTIONS: See \link ViscosityModel_Map \endlink \n DEFAULT: DAVIDSON \ingroup Config*/
12281228
addEnumOption("MIXING_VISCOSITY_MODEL", Kind_MixingViscosityModel, MixingViscosityModel_Map, MIXINGVISCOSITYMODEL::DAVIDSON);
@@ -3798,7 +3798,7 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i
37983798

37993799
if (Kind_FluidModel == FLUID_MIXTURE) {
38003800
/*--- Check whether the number of entries of each specified fluid property equals the number of transported scalar
3801-
equations solved + 1. nMolecular_Weight and nSpecific_Heat_Cp are used because it is required for the fluid mixing models.
3801+
equations solved + 1. nMolecular_Weight and nSpecific_Heat_Cp are used because it is required for the fluid mixing models.
38023802
* Cp is required in case of MIXTURE_FLUID_MODEL because the energy equation needs to be active.--- */
38033803
if (nMolecular_Weight != nSpecies_Init + 1 || nSpecific_Heat_Cp != nSpecies_Init + 1) {
38043804
SU2_MPI::Error(
@@ -6054,10 +6054,10 @@ void CConfig::SetOutput(SU2_COMPONENT val_software, unsigned short val_izone) {
60546054
}
60556055
cout << "." << endl;
60566056
break;
6057-
}
6057+
}
60586058
switch (Kind_Trans_Model) {
60596059
case TURB_TRANS_MODEL::NONE: break;
6060-
case TURB_TRANS_MODEL::LM: cout << "Transition model: Langtry and Menter's 4 equation model (2009)" << endl; break;
6060+
case TURB_TRANS_MODEL::LM: cout << "Transition model: Langtry and Menter's 4 equation model (2009)" << endl; break;
60616061
}
60626062
cout << "Hybrid RANS/LES: ";
60636063
switch (Kind_HybridRANSLES) {

‎SU2_CFD/include/fluid/CFluidModel.hpp‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -83,7 +83,7 @@ class CFluidModel {
8383
* \brief Instantiate the right type of conductivity model based on config.
8484
*/
8585
static unique_ptr<CConductivityModel> MakeThermalConductivityModel(const CConfig* config, unsigned short iSpecies);
86-
86+
8787
/*!
8888
* \brief Instantiate the right type of mass diffusivity model based on config.
8989
*/
@@ -167,7 +167,7 @@ class CFluidModel {
167167
mass_diffusivity = MassDiffusivity->GetDiffusivity();
168168
return mass_diffusivity;
169169
}
170-
170+
171171
/*!
172172
* \brief Get fluid pressure partial derivative.
173173
*/
@@ -242,7 +242,7 @@ class CFluidModel {
242242
* \brief Set thermal conductivity model.
243243
*/
244244
virtual void SetThermalConductivityModel(const CConfig* config);
245-
245+
246246
/*!
247247
* \brief Set mass diffusivity model.
248248
*/

‎SU2_CFD/include/variables/CEulerVariable.hpp‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -39,7 +39,7 @@
3939
* \author F. Palacios, T. Economon
4040
*/
4141
class CEulerVariable : public CFlowVariable {
42-
public:
42+
public:
4343
static constexpr size_t MAXNVAR = 12;
4444

4545
template <class IndexType>

‎SU2_CFD/include/variables/CPrimitiveIndices.hpp‎

Lines changed: 92 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,92 @@
1+
/*!
2+
* \file CPrimitiveIndices.hpp
3+
* \brief Abstract representation of flow primitive variable indices that tries to be efficient.
4+
* \version 7.4.0 "Blackbird"
5+
*
6+
* SU2 Project Website: https://su2code.github.io
7+
*
8+
* The SU2 Project is maintained by the SU2 Foundation
9+
* (http://su2foundation.org)
10+
*
11+
* Copyright 2012-2022, SU2 Contributors (cf. AUTHORS.md)
12+
*
13+
* SU2 is free software; you can redistribute it and/or
14+
* modify it under the terms of the GNU Lesser General Public
15+
* License as published by the Free Software Foundation; either
16+
* version 2.1 of the License, or (at your option) any later version.
17+
*
18+
* SU2 is distributed in the hope that it will be useful,
19+
* but WITHOUT ANY WARRANTY; without even the implied warranty of
20+
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
21+
* Lesser General Public License for more details.
22+
*
23+
* You should have received a copy of the GNU Lesser General Public
24+
* License along with SU2. If not, see <http://www.gnu.org/licenses/>.
25+
*/
26+
27+
#pragma once
28+
29+
#include <cstdint>
30+
#include <cstring>
31+
#include "CEulerVariable.hpp"
32+
#include "CIncEulerVariable.hpp"
33+
#include "CNEMOEulerVariable.hpp"
34+
35+
/*!
36+
* \class CPrimitiveIndices
37+
* \brief Abstract wrapper for the CIndices classes of CEulerVariable, CIncEulerVariable, etc.
38+
* This should only be used when the concrete type of solver variable is not known.
39+
* \note We are using some tricks to try to make this more efficient than virtual functions.
40+
*/
41+
template <class IndexType>
42+
struct CPrimitiveIndices {
43+
CPrimitiveIndices(bool incompressible, bool nemo, IndexType nDim, IndexType nSpecies) {
44+
/*--- Instantiate the right type of CIndices, without heap allocations. ---*/
45+
if (incompressible) {
46+
type_ = 0;
47+
Construct<CIncEulerVariable::template CIndices<IndexType>>(nDim, nSpecies);
48+
} else if (nemo) {
49+
type_ = 1;
50+
Construct<CNEMOEulerVariable::template CIndices<IndexType>>(nDim, nSpecies);
51+
} else {
52+
type_ = 2;
53+
Construct<CEulerVariable::template CIndices<IndexType>>(nDim, nSpecies);
54+
}
55+
}
56+
57+
/*--- To use the CIndices object we stored on the buffer, we cast to the
58+
* concrete type using the same logic used during construction. ---*/
59+
#define GET_INDEX_IMPL(NAME) \
60+
inline IndexType NAME() const { \
61+
switch (type_) { \
62+
case 0: \
63+
return reinterpret_cast<const CIncEulerVariable::template CIndices<IndexType>*>(data_)->NAME(); \
64+
case 1: \
65+
return reinterpret_cast<const CNEMOEulerVariable::template CIndices<IndexType>*>(data_)->NAME(); \
66+
default: \
67+
return reinterpret_cast<const CEulerVariable::template CIndices<IndexType>*>(data_)->NAME(); \
68+
} \
69+
}
70+
GET_INDEX_IMPL(Temperature)
71+
GET_INDEX_IMPL(Velocity)
72+
GET_INDEX_IMPL(Pressure)
73+
GET_INDEX_IMPL(Density)
74+
GET_INDEX_IMPL(Enthalpy)
75+
GET_INDEX_IMPL(SoundSpeed)
76+
GET_INDEX_IMPL(LaminarViscosity)
77+
GET_INDEX_IMPL(EddyViscosity)
78+
GET_INDEX_IMPL(ThermalConductivity)
79+
GET_INDEX_IMPL(CpTotal)
80+
GET_INDEX_IMPL(Temperature_ve)
81+
82+
private:
83+
template <class ConcreteIndices>
84+
void Construct(IndexType nDim, IndexType nSpecies) {
85+
/*--- Build the indices object in the static buffer owned by this class. ---*/
86+
static_assert(sizeof(ConcreteIndices) <= 2 * sizeof(IndexType), "");
87+
new(data_) ConcreteIndices(nDim, nSpecies);
88+
}
89+
90+
alignas(sizeof(IndexType)) char data_[2 * sizeof(IndexType)]{}; /*!< \brief Space for the largest CIndices class. */
91+
uint8_t type_;
92+
};

‎SU2_CFD/src/solvers/CTurbSSTSolver.cpp‎

Lines changed: 36 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
#include "../../include/solvers/CTurbSSTSolver.hpp"
2929
#include "../../include/variables/CTurbSSTVariable.hpp"
3030
#include "../../include/variables/CFlowVariable.hpp"
31+
#include "../../include/variables/CPrimitiveIndices.hpp"
3132
#include "../../../Common/include/parallelization/omp_structure.hpp"
3233
#include "../../../Common/include/toolboxes/geometry_toolbox.hpp"
3334

@@ -322,8 +323,8 @@ void CTurbSSTSolver::Source_Residual(CGeometry *geometry, CSolver **solver_conta
322323
numerics->SetCrossDiff(nodes->GetCrossDiff(iPoint));
323324

324325
/*--- Effective Intermittency ---*/
325-
326-
if (TURB_TRANS_MODEL::LM == config->GetKind_Trans_Model()) {
326+
327+
if (TURB_TRANS_MODEL::LM == config->GetKind_Trans_Model()) {
327328
numerics->SetIntermittencyEff(solver_container[TRANS_SOL]->GetNodes()->GetIntermittencyEff(iPoint));
328329
}
329330

@@ -578,6 +579,9 @@ void CTurbSSTSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container, C
578579

579580
const bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
580581

582+
const CPrimitiveIndices<unsigned short> prim_idx(config->GetKind_Regime() == ENUM_REGIME::INCOMPRESSIBLE,
583+
config->GetNEMOProblem(), nDim, config->GetnSpecies());
584+
581585
/*--- Loop over all the vertices on this boundary marker ---*/
582586

583587
SU2_OMP_FOR_STAT(OMP_MIN_SIZE)
@@ -608,13 +612,38 @@ void CTurbSSTSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container, C
608612

609613
conv_numerics->SetPrimitive(V_domain, V_inlet);
610614

611-
/*--- Non-dimensionalize Inlet_TurbVars if Inlet-Files are used. ---*/
612615
su2double Inlet_Vars[MAXNVAR];
613-
Inlet_Vars[0] = Inlet_TurbVars[val_marker][iVertex][0];
614-
Inlet_Vars[1] = Inlet_TurbVars[val_marker][iVertex][1];
615616
if (config->GetInlet_Profile_From_File()) {
616-
Inlet_Vars[0] /= pow(config->GetVelocity_Ref(), 2);
617-
Inlet_Vars[1] *= config->GetViscosity_Ref() / (config->GetDensity_Ref() * pow(config->GetVelocity_Ref(), 2));
617+
/*--- Non-dimensionalize Inlet_TurbVars if Inlet-Files are used. ---*/
618+
Inlet_Vars[0] = Inlet_TurbVars[val_marker][iVertex][0] / pow(config->GetVelocity_Ref(), 2);
619+
Inlet_Vars[1] = Inlet_TurbVars[val_marker][iVertex][1] * config->GetViscosity_Ref() /
620+
(config->GetDensity_Ref() * pow(config->GetVelocity_Ref(), 2));
621+
} else {
622+
/*--- Obtain fluid model for computing the kine and omega to impose at the inlet boundary. ---*/
623+
CFluidModel* FluidModel = solver_container[FLOW_SOL]->GetFluidModel();
624+
625+
/*--- Obtain flow velocity vector at inlet boundary node ---*/
626+
627+
const su2double* Velocity_Inlet = &V_inlet[prim_idx.Velocity()];
628+
su2double Density_Inlet;
629+
if (config->GetKind_Regime() == ENUM_REGIME::COMPRESSIBLE) {
630+
Density_Inlet = V_inlet[prim_idx.Density()];
631+
FluidModel->SetTDState_Prho(V_inlet[prim_idx.Pressure()], Density_Inlet);
632+
} else {
633+
const su2double* Scalar_Inlet = nullptr;
634+
if (config->GetKind_Species_Model() != SPECIES_MODEL::NONE) {
635+
Scalar_Inlet = config->GetInlet_SpeciesVal(config->GetMarker_All_TagBound(val_marker));
636+
}
637+
FluidModel->SetTDState_T(V_inlet[prim_idx.Temperature()], Scalar_Inlet);
638+
Density_Inlet = FluidModel->GetDensity();
639+
}
640+
const su2double Laminar_Viscosity_Inlet = FluidModel->GetLaminarViscosity();
641+
const su2double Intensity = config->GetTurbulenceIntensity_FreeStream();
642+
const su2double viscRatio = config->GetTurb2LamViscRatio_FreeStream();
643+
const su2double VelMag2 = GeometryToolbox::SquaredNorm(nDim, Velocity_Inlet);
644+
645+
Inlet_Vars[0] = 3.0 / 2.0 * (VelMag2 * pow(Intensity, 2));
646+
Inlet_Vars[1] = Density_Inlet * Inlet_Vars[0] / (Laminar_Viscosity_Inlet * viscRatio);
618647
}
619648

620649
/*--- Set the turbulent variable states. Use free-stream SST

‎TestCases/hybrid_regression.py‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -254,7 +254,7 @@ def main():
254254
axi_rans_air_nozzle_restart.cfg_dir = "axisymmetric_rans/air_nozzle"
255255
axi_rans_air_nozzle_restart.cfg_file = "air_nozzle_restart.cfg"
256256
axi_rans_air_nozzle_restart.test_iter = 10
257-
axi_rans_air_nozzle_restart.test_vals = [-12.094405, -6.649212, -8.899548, -2.491859, -1924.900000]
257+
axi_rans_air_nozzle_restart.test_vals = [-12.085749, -7.439410, -8.772091, -4.045091, -1924.800000]
258258
axi_rans_air_nozzle_restart.test_vals_aarch64 = [-12.093539, -6.630357, -8.798732, -2.399130, -1938.200000]
259259
test_list.append(axi_rans_air_nozzle_restart)
260260

@@ -633,7 +633,7 @@ def main():
633633
bars_SST_2D.cfg_dir = "sliding_interface/bars_SST_2D"
634634
bars_SST_2D.cfg_file = "bars.cfg"
635635
bars_SST_2D.test_iter = 13
636-
bars_SST_2D.test_vals = [13.000000, -0.609170, -1.523885]
636+
bars_SST_2D.test_vals = [13.000000, -0.604409, -1.523885]
637637
bars_SST_2D.multizone = True
638638
test_list.append(bars_SST_2D)
639639

‎TestCases/parallel_regression.py‎

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -394,7 +394,7 @@ def main():
394394
axi_rans_air_nozzle_restart.cfg_dir = "axisymmetric_rans/air_nozzle"
395395
axi_rans_air_nozzle_restart.cfg_file = "air_nozzle_restart.cfg"
396396
axi_rans_air_nozzle_restart.test_iter = 10
397-
axi_rans_air_nozzle_restart.test_vals = [-12.098072, -6.657157, -8.872565, -2.491820, -1924.900000]
397+
axi_rans_air_nozzle_restart.test_vals = [-12.089268, -7.493381, -8.716391, -4.021218, -1924.800000]
398398
axi_rans_air_nozzle_restart.tol = 0.0001
399399
test_list.append(axi_rans_air_nozzle_restart)
400400

@@ -919,7 +919,7 @@ def main():
919919
coolprop_nozzle.cfg_dir = "nicf/coolprop"
920920
coolprop_nozzle.cfg_file = "coolprop_nozzle.cfg"
921921
coolprop_nozzle.test_iter = 10
922-
coolprop_nozzle.test_vals = [-4.692515, -3.779307, 3.419174, 0.000000, 0.000000]
922+
coolprop_nozzle.test_vals = [-4.692149, -1.965079, 3.428130, 0.000000, 0.000000]
923923
test_list.append(coolprop_nozzle)
924924

925925
######################################
@@ -1025,7 +1025,7 @@ def main():
10251025
bars_SST_2D.cfg_dir = "sliding_interface/bars_SST_2D"
10261026
bars_SST_2D.cfg_file = "bars.cfg"
10271027
bars_SST_2D.test_iter = 13
1028-
bars_SST_2D.test_vals = [13.000000, -0.609170, -1.523885]
1028+
bars_SST_2D.test_vals = [13.000000, -0.604409, -1.523885]
10291029
bars_SST_2D.command = TestCase.Command(exec = "SU2_CFD")
10301030
bars_SST_2D.multizone = True
10311031
test_list.append(bars_SST_2D)
@@ -1244,7 +1244,7 @@ def main():
12441244
pywrapper_unsteadyCHT.cfg_dir = "py_wrapper/flatPlate_unsteady_CHT"
12451245
pywrapper_unsteadyCHT.cfg_file = "unsteady_CHT_FlatPlate_Conf.cfg"
12461246
pywrapper_unsteadyCHT.test_iter = 5
1247-
pywrapper_unsteadyCHT.test_vals = [-1.614167, 2.245726, -0.001241, 0.175715]
1247+
pywrapper_unsteadyCHT.test_vals = [-1.614167, 2.245725, -0.001241, 0.175713]
12481248
pywrapper_unsteadyCHT.command = TestCase.Command("mpirun -np 2", "python", "launch_unsteady_CHT_FlatPlate.py --parallel -f")
12491249
pywrapper_unsteadyCHT.unsteady = True
12501250
pywrapper_unsteadyCHT.new_output = True
@@ -1330,7 +1330,7 @@ def main():
13301330
species2_primitiveVenturi_mixingmodel.cfg_dir = "species_transport/venturi_primitive_3species"
13311331
species2_primitiveVenturi_mixingmodel.cfg_file = "species2_primitiveVenturi_mixingmodel.cfg"
13321332
species2_primitiveVenturi_mixingmodel.test_iter = 50
1333-
species2_primitiveVenturi_mixingmodel.test_vals = [-5.476757, -4.585053, -4.579453, -5.693611, -0.069507, -5.630568, 5.000000, -1.667485, 5.000000, -4.865481, 5.000000, -1.168603, 0.000425, 0.000405, 0.000020, 0.000000]
1333+
species2_primitiveVenturi_mixingmodel.test_vals = [-5.477173, -4.589710, -4.582462, -5.690334, -0.069005, -5.631182, 5.000000, -1.668297, 5.000000, -4.866834, 5.000000, -1.168385, 0.000425, 0.000404, 0.000021, 0.000000]
13341334
species2_primitiveVenturi_mixingmodel.new_output = True
13351335
test_list.append(species2_primitiveVenturi_mixingmodel)
13361336

@@ -1339,7 +1339,7 @@ def main():
13391339
species2_primitiveVenturi_mixingmodel_viscosity.cfg_dir = "species_transport/venturi_primitive_3species"
13401340
species2_primitiveVenturi_mixingmodel_viscosity.cfg_file = "species2_primitiveVenturi_mixingmodel_viscosity.cfg"
13411341
species2_primitiveVenturi_mixingmodel_viscosity.test_iter = 50
1342-
species2_primitiveVenturi_mixingmodel_viscosity.test_vals = [-4.421969, -4.027438, -4.346703, -5.465691, 0.436661, -4.685235, 5.000000, -1.816533, 5.000000, -5.298965, 5.000000, -1.724682, 2.287735, 0.971641, 0.607256, 0.708838]
1342+
species2_primitiveVenturi_mixingmodel_viscosity.test_vals = [-4.415103, -4.019010, -4.341952, -5.471507, 0.437029, -4.685683, 5.000000, -1.814481, 5.000000, -5.282862, 5.000000, -1.717346, 2.280789, 0.971585, 0.607256, 0.701947]
13431343
species2_primitiveVenturi_mixingmodel_viscosity.new_output = True
13441344
test_list.append(species2_primitiveVenturi_mixingmodel_viscosity)
13451345

@@ -1348,7 +1348,7 @@ def main():
13481348
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg_dir = "species_transport/venturi_primitive_3species"
13491349
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg_file = "species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg"
13501350
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_iter = 50
1351-
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_vals = [-6.123576, -4.998143, -4.888374, -7.379916, 2.465109, -5.639096, 30.000000, -5.719962, 12.000000, -8.161358, 10.000000, -8.602778, 2.084606, 1.000000, 0.600000, 0.484606]
1351+
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_vals = [-6.119300, -4.997334, -4.886952, -7.382047, 2.439615, -5.627780, 30.000000, -5.723010, 12.000000, -8.145296, 9.000000, -8.075817, 2.084563, 1.000000, 0.600000, 0.484563]
13521352
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.su2_exec = "mpirun -n 2 SU2_CFD"
13531353
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.timeout = 1600
13541354
species2_primitiveVenturi_mixingmodel_heatcapacity_H2.new_output = True
@@ -1360,7 +1360,7 @@ def main():
13601360
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg_dir = "species_transport/venturi_primitive_3species"
13611361
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg_file = "species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg"
13621362
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_iter = 50
1363-
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_vals = [-5.729433, -5.302965, -5.193202, -8.383664, 2.160286, -5.244935, 30.000000, -5.719915, 12.000000, -8.161218, 10.000000, -8.602842, 2.084608, 1.000000, 0.600000, 0.484608]
1363+
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_vals = [-5.725157, -5.302154, -5.191782, -8.385791, 2.134792, -5.233618, 30.000000, -5.722966, 12.000000, -8.145150, 9.000000, -8.075867, 2.084565, 1.000000, 0.600000, 0.484565]
13641364
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.su2_exec = "mpirun -n 2 SU2_CFD"
13651365
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.timeout = 1600
13661366
species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.new_output = True

‎TestCases/parallel_regression_AD.py‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -150,7 +150,7 @@ def main():
150150
discadj_axisymmetric_rans_nozzle.cfg_dir = "axisymmetric_rans/air_nozzle"
151151
discadj_axisymmetric_rans_nozzle.cfg_file = "air_nozzle_restart.cfg"
152152
discadj_axisymmetric_rans_nozzle.test_iter = 10
153-
discadj_axisymmetric_rans_nozzle.test_vals = [9.524990, 5.023431, 9.417118, 2.545450, 0.000000, -246129999999999991087104.000000]
153+
discadj_axisymmetric_rans_nozzle.test_vals = [9.523445, 5.023948, 7.509000, 2.841876, 0.000000, -246139999999999992659968.000000]
154154
discadj_axisymmetric_rans_nozzle.no_restart = True
155155
test_list.append(discadj_axisymmetric_rans_nozzle)
156156

‎TestCases/serial_regression.py‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -264,7 +264,7 @@ def main():
264264
turb_wallfunction_flatplate_sst.cfg_dir = "wallfunctions/flatplate/compressible_SST"
265265
turb_wallfunction_flatplate_sst.cfg_file = "turb_SST_flatplate.cfg"
266266
turb_wallfunction_flatplate_sst.test_iter = 10
267-
turb_wallfunction_flatplate_sst.test_vals = [-4.229952, -1.930555, -1.998476, 1.250387, -1.635646, 1.462394, 10.000000, -2.151907, 0.072873, 0.002514] #last 10 columns
267+
turb_wallfunction_flatplate_sst.test_vals = [-4.229999, -1.930714, -1.998514, 1.250334, -1.635256, 1.474343, 10.000000, -2.152260, 0.072873, 0.002514] #last 10 columns
268268
test_list.append(turb_wallfunction_flatplate_sst)
269269

270270
# FLAT PLATE, WALL FUNCTIONS, COMPRESSIBLE SA
@@ -353,7 +353,7 @@ def main():
353353
axi_rans_air_nozzle_restart.cfg_dir = "axisymmetric_rans/air_nozzle"
354354
axi_rans_air_nozzle_restart.cfg_file = "air_nozzle_restart.cfg"
355355
axi_rans_air_nozzle_restart.test_iter = 10
356-
axi_rans_air_nozzle_restart.test_vals = [-12.098379, -6.656873, -8.869081, -2.491818, -1924.900000]
356+
axi_rans_air_nozzle_restart.test_vals = [-12.088204, -7.497812, -8.714688, -4.019578, -1924.800000]
357357
axi_rans_air_nozzle_restart.tol = 0.0001
358358
test_list.append(axi_rans_air_nozzle_restart)
359359

@@ -1023,7 +1023,7 @@ def main():
10231023
bars_SST_2D.cfg_dir = "sliding_interface/bars_SST_2D"
10241024
bars_SST_2D.cfg_file = "bars.cfg"
10251025
bars_SST_2D.test_iter = 13
1026-
bars_SST_2D.test_vals = [13.000000, -0.609170, -1.523885]
1026+
bars_SST_2D.test_vals = [13.000000, -0.604409, -1.523885]
10271027
bars_SST_2D.multizone = True
10281028
test_list.append(bars_SST_2D)
10291029

0 commit comments

Comments
 (0)