From efcf8a5e703d0469bb73b51c8d2eb64b6c828815 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 15 Jan 2026 22:33:36 +0100 Subject: [PATCH 1/4] initial commit --- Common/include/CConfig.hpp | 8 +++++++ Common/include/option_structure.hpp | 12 ++++++++++ Common/src/CConfig.cpp | 4 ++++ SU2_CFD/src/solvers/CIncEulerSolver.cpp | 24 ++++++++----------- .../lam_prem_ch4_cfd.cfg | 1 + .../lam_prem_ch4_hx_ad.cfg | 1 + .../lam_prem_ch4_hx_dot.cfg | 1 + .../lam_prem_ch4_cht_cfd_fluid.cfg | 1 + .../lam_prem_ch4_cht_ad_fluid.cfg | 1 + .../lam_prem_ch4_cfd_axi.cfg | 1 + .../lam_partial_prem_ch4_cfd.cfg | 1 + .../laminar_premixed_h2_flame_cfd.cfg | 1 + .../species_transport/multizone/zone1.cfg | 3 ++- .../species_transport/multizone/zone2.cfg | 3 ++- .../active_species_transport_temp_limits.cfg | 1 + .../species2_primitiveVenturi_JST.cfg | 1 + ...pecies2_primitiveVenturi_Lax_Friedrich.cfg | 1 + .../species2_primitiveVenturi_mixingmodel.cfg | 1 + ...eVenturi_mixingmodel_TURBULENT_MARKERS.cfg | 1 + ...itiveVenturi_mixingmodel_boundedscalar.cfg | 1 + ...iveVenturi_mixingmodel_heatcapacity_H2.cfg | 1 + ...Venturi_mixingmodel_heatcapacity_H2_ND.cfg | 1 + ...primitiveVenturi_mixingmodel_viscosity.cfg | 1 + .../sandia_jet/validation.cfg | 1 + 24 files changed, 56 insertions(+), 16 deletions(-) diff --git a/Common/include/CConfig.hpp b/Common/include/CConfig.hpp index f198df434761..92d255973aa0 100644 --- a/Common/include/CConfig.hpp +++ b/Common/include/CConfig.hpp @@ -523,6 +523,7 @@ class CConfig { LINEAR_SOLVER_INNER Kind_Linear_Solver_Inner; /*!< \brief Inner solver used in nested Krylov schemes. */ unsigned short Kind_FluidModel, /*!< \brief Kind of the Fluid Model: Ideal, van der Waals, etc. */ Kind_InitOption, /*!< \brief Kind of Init option to choose if initializing with Reynolds number or with thermodynamic conditions */ + Kind_InitOption_Inc, /*!< \brief Kind of Init option for incompressible solver, initializing with initial density or operating pressure */ Kind_GridMovement, /*!< \brief Kind of the static mesh movement. */ *Kind_SurfaceMovement, /*!< \brief Kind of the static mesh movement. */ nKind_SurfaceMovement, /*!< \brief Kind of the dynamic mesh movement. */ @@ -4044,6 +4045,13 @@ class CConfig { * \return free stream option */ unsigned short GetKind_InitOption(void) const { return Kind_InitOption; } + + /*! + * \brief Get kind of initialization option for incompressible flows. + * \return Kind initialization option for incompressible flows. + */ + unsigned short GetKind_InitOptionInc(void) const { return Kind_InitOption_Inc; } + /*! * \brief Get the value of the critical pressure. * \return Critical pressure. diff --git a/Common/include/option_structure.hpp b/Common/include/option_structure.hpp index 690350fae050..1ba3279bc3b8 100644 --- a/Common/include/option_structure.hpp +++ b/Common/include/option_structure.hpp @@ -651,6 +651,18 @@ static const MapType InitOption_Map = { MakePair("TD_CONDITIONS", TD_CONDITIONS) }; +/*! + * \brief Types of initialization option incompressible solver + */ +enum ENUM_INIT_OPTION_INC { + DENSITY_INIT = 0, /*!< \brief Density initalization. */ + OPERATING_PRESSURE = 1 /*!< \brief Operating pressure initalization. */ +}; +static const MapType InitOptionInc_Map = { + MakePair("DENSITY_INIT", DENSITY_INIT) + MakePair("OPERATING_PRESSURE", OPERATING_PRESSURE) +}; + /*! * \brief Types of freestream specification */ diff --git a/Common/src/CConfig.cpp b/Common/src/CConfig.cpp index a9ecb3a94b24..cd3608c8760d 100644 --- a/Common/src/CConfig.cpp +++ b/Common/src/CConfig.cpp @@ -1369,6 +1369,10 @@ void CConfig::SetConfig_Options() { addDoubleOption("MACH_NUMBER", Mach, 0.0); /*!\brief INIT_OPTION \n DESCRIPTION: Init option to choose between Reynolds or thermodynamics quantities for initializing the solution \n OPTIONS: see \link InitOption_Map \endlink \n DEFAULT REYNOLDS \ingroup Config*/ addEnumOption("INIT_OPTION", Kind_InitOption, InitOption_Map, REYNOLDS); + /*!\brief INIT_OPTION_INC \n DESCRIPTION: Init option for incompressible solver to choose between initial density or + * operating pressure initialization of the solution \n OPTIONS: see \link InitOption_Map \endlink \n DEFAULT REYNOLDS + * \ingroup Config*/ + addEnumOption("INIT_OPTION_INC", Kind_InitOption_Inc, InitOptionInc_Map, DENSITY_INIT); /* DESCRIPTION: Free-stream option to choose between density and temperature for initializing the solution */ addEnumOption("FREESTREAM_OPTION", Kind_FreeStreamOption, FreeStreamOption_Map, FREESTREAM_OPTION::TEMPERATURE_FS); /*!\brief FREESTREAM_PRESSURE\n DESCRIPTION: Free-stream pressure (101325.0 N/m^2 by default) \ingroup Config*/ diff --git a/SU2_CFD/src/solvers/CIncEulerSolver.cpp b/SU2_CFD/src/solvers/CIncEulerSolver.cpp index c99f87896c2d..5d2eecddc2ed 100644 --- a/SU2_CFD/src/solvers/CIncEulerSolver.cpp +++ b/SU2_CFD/src/solvers/CIncEulerSolver.cpp @@ -260,6 +260,7 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i bool tkeNeeded = ((turbulent) && ((config->GetKind_Turb_Model() == TURB_MODEL::SST))); bool energy = config->GetEnergy_Equation(); bool boussinesq = (config->GetKind_DensityModel() == INC_DENSITYMODEL::BOUSSINESQ); + bool density_init = (config->GetKind_InitOptionInc() == DENSITY_INIT); /*--- Compute dimensional free-stream values. ---*/ @@ -283,6 +284,14 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i } ModVel_FreeStream = sqrt(ModVel_FreeStream); config->SetModVel_FreeStream(ModVel_FreeStream); + config->SetGas_Constant(UNIVERSAL_GAS_CONSTANT / (config->GetMolecular_Weight() / 1000.0)); + if (density_init) { + Pressure_Thermodynamic = Density_FreeStream * Temperature_FreeStream * config->GetGas_Constant(); + } else { + Pressure_Thermodynamic = config->GetPressure_Thermodynamic(); + } + config->SetPressure_Thermodynamic(Pressure_Thermodynamic); + CFluidModel* auxFluidModel = nullptr; switch (config->GetKind_FluidModel()) { @@ -295,18 +304,12 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i case INC_IDEAL_GAS: - config->SetGas_Constant(UNIVERSAL_GAS_CONSTANT/(config->GetMolecular_Weight()/1000.0)); - Pressure_Thermodynamic = Density_FreeStream*Temperature_FreeStream*config->GetGas_Constant(); auxFluidModel = new CIncIdealGas(config->GetSpecific_Heat_Cp(), config->GetGas_Constant(), Pressure_Thermodynamic, STD_REF_TEMP); - auxFluidModel->SetTDState_T(Temperature_FreeStream); - Pressure_Thermodynamic = auxFluidModel->GetPressure(); - config->SetPressure_Thermodynamic(Pressure_Thermodynamic); + auxFluidModel->SetTDState_T(Temperature_FreeStream); break; case INC_IDEAL_GAS_POLY: - config->SetGas_Constant(UNIVERSAL_GAS_CONSTANT/(config->GetMolecular_Weight()/1000.0)); - Pressure_Thermodynamic = Density_FreeStream*Temperature_FreeStream*config->GetGas_Constant(); auxFluidModel = new CIncIdealGasPolynomial(config->GetGas_Constant(), Pressure_Thermodynamic, STD_REF_TEMP); if (viscous) { /*--- Variable Cp model via polynomial. ---*/ @@ -315,24 +318,17 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i auxFluidModel->SetCpModel(config, Temperature_FreeStream); } auxFluidModel->SetTDState_T(Temperature_FreeStream); - Pressure_Thermodynamic = auxFluidModel->GetPressure(); - config->SetPressure_Thermodynamic(Pressure_Thermodynamic); break; case FLUID_MIXTURE: - config->SetGas_Constant(UNIVERSAL_GAS_CONSTANT / (config->GetMolecular_Weight() / 1000.0)); - Pressure_Thermodynamic = config->GetPressure_Thermodynamic(); auxFluidModel = new CFluidScalar(Pressure_Thermodynamic, config); auxFluidModel->SetTDState_T(Temperature_FreeStream, config->GetSpecies_Init()); break; case FLUID_FLAMELET: - config->SetGas_Constant(UNIVERSAL_GAS_CONSTANT / (config->GetMolecular_Weight() / 1000.0)); - Pressure_Thermodynamic = config->GetPressure_Thermodynamic(); auxFluidModel = new CFluidFlamelet(config, Pressure_Thermodynamic); - config->SetPressure_Thermodynamic(Pressure_Thermodynamic); auxFluidModel->SetTDState_T(Temperature_FreeStream, config->GetSpecies_Init()); break; diff --git a/TestCases/flamelet/01_laminar_premixed_ch4_flame_cfd/lam_prem_ch4_cfd.cfg b/TestCases/flamelet/01_laminar_premixed_ch4_flame_cfd/lam_prem_ch4_cfd.cfg index 12d6c6bf8615..ba36c610390a 100644 --- a/TestCases/flamelet/01_laminar_premixed_ch4_flame_cfd/lam_prem_ch4_cfd.cfg +++ b/TestCases/flamelet/01_laminar_premixed_ch4_flame_cfd/lam_prem_ch4_cfd.cfg @@ -23,6 +23,7 @@ INC_DENSITY_INIT= 1.00 INC_VELOCITY_INIT= (0.5, 0.0, 0.0 ) INC_TEMPERATURE_INIT= 300.0 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID MODEL --------------------------------------- % % diff --git a/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_ad.cfg b/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_ad.cfg index 2c8724cafd46..c54b73896d37 100644 --- a/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_ad.cfg +++ b/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_ad.cfg @@ -33,6 +33,7 @@ INC_VELOCITY_INIT= ( 0.2, 0, 0 ) INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_dot.cfg b/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_dot.cfg index 120171d8e27c..101d0465d4bd 100644 --- a/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_dot.cfg +++ b/TestCases/flamelet/02_laminar_premixed_ch4_flame_hx_ad/lam_prem_ch4_hx_dot.cfg @@ -34,6 +34,7 @@ INC_VELOCITY_INIT= ( 0.2, 0, 0 ) INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/flamelet/03_laminar_premixed_ch4_flame_cht_cfd/lam_prem_ch4_cht_cfd_fluid.cfg b/TestCases/flamelet/03_laminar_premixed_ch4_flame_cht_cfd/lam_prem_ch4_cht_cfd_fluid.cfg index 005bea4123af..af34a6d990f6 100644 --- a/TestCases/flamelet/03_laminar_premixed_ch4_flame_cht_cfd/lam_prem_ch4_cht_cfd_fluid.cfg +++ b/TestCases/flamelet/03_laminar_premixed_ch4_flame_cht_cfd/lam_prem_ch4_cht_cfd_fluid.cfg @@ -50,6 +50,7 @@ INC_VELOCITY_INIT= ( 0.2, 0, 0 ) INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % -------------------- FLUID MODEL --------------------------------------- % diff --git a/TestCases/flamelet/04_laminar_premixed_ch4_flame_cht_ad/lam_prem_ch4_cht_ad_fluid.cfg b/TestCases/flamelet/04_laminar_premixed_ch4_flame_cht_ad/lam_prem_ch4_cht_ad_fluid.cfg index 8b6b374d2046..b51b51f18a5d 100644 --- a/TestCases/flamelet/04_laminar_premixed_ch4_flame_cht_ad/lam_prem_ch4_cht_ad_fluid.cfg +++ b/TestCases/flamelet/04_laminar_premixed_ch4_flame_cht_ad/lam_prem_ch4_cht_ad_fluid.cfg @@ -52,6 +52,7 @@ INC_VELOCITY_INIT= ( 0.2, 0, 0 ) INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % -------------------- FLUID MODEL --------------------------------------- % diff --git a/TestCases/flamelet/05_laminar_premixed_ch4_flame_cfd_axi/lam_prem_ch4_cfd_axi.cfg b/TestCases/flamelet/05_laminar_premixed_ch4_flame_cfd_axi/lam_prem_ch4_cfd_axi.cfg index d8ca68d94fd9..3e403aad7cad 100644 --- a/TestCases/flamelet/05_laminar_premixed_ch4_flame_cfd_axi/lam_prem_ch4_cfd_axi.cfg +++ b/TestCases/flamelet/05_laminar_premixed_ch4_flame_cfd_axi/lam_prem_ch4_cfd_axi.cfg @@ -24,6 +24,7 @@ INC_DENSITY_INIT= 1.00 INC_VELOCITY_INIT= (0.5, 0.0, 0.0 ) INC_TEMPERATURE_INIT= 300.0 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % AXISYMMETRIC= YES diff --git a/TestCases/flamelet/06_laminar_partial_premixed_ch4_flame_cfd/lam_partial_prem_ch4_cfd.cfg b/TestCases/flamelet/06_laminar_partial_premixed_ch4_flame_cfd/lam_partial_prem_ch4_cfd.cfg index 6d5fded99026..ea7b1ab8c696 100644 --- a/TestCases/flamelet/06_laminar_partial_premixed_ch4_flame_cfd/lam_partial_prem_ch4_cfd.cfg +++ b/TestCases/flamelet/06_laminar_partial_premixed_ch4_flame_cfd/lam_partial_prem_ch4_cfd.cfg @@ -23,6 +23,7 @@ INC_DENSITY_INIT= 1.00 INC_VELOCITY_INIT= (0.7, 0.0, 0.0 ) INC_TEMPERATURE_INIT= 300.0 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID MODEL --------------------------------------- % % diff --git a/TestCases/flamelet/07_laminar_premixed_h2_flame_cfd/laminar_premixed_h2_flame_cfd.cfg b/TestCases/flamelet/07_laminar_premixed_h2_flame_cfd/laminar_premixed_h2_flame_cfd.cfg index 93328224dc0b..4d21f38ea7fa 100644 --- a/TestCases/flamelet/07_laminar_premixed_h2_flame_cfd/laminar_premixed_h2_flame_cfd.cfg +++ b/TestCases/flamelet/07_laminar_premixed_h2_flame_cfd/laminar_premixed_h2_flame_cfd.cfg @@ -22,6 +22,7 @@ INC_VELOCITY_INIT= (1.13, 0.0, 0.0 ) INC_TEMPERATURE_INIT= 300.0 THERMODYNAMIC_PRESSURE= 101325 INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % -------------------- FLUID MODEL --------------------------------------- % % diff --git a/TestCases/species_transport/multizone/zone1.cfg b/TestCases/species_transport/multizone/zone1.cfg index c18f0565b1e0..08cf8c8e418d 100644 --- a/TestCases/species_transport/multizone/zone1.cfg +++ b/TestCases/species_transport/multizone/zone1.cfg @@ -31,7 +31,8 @@ INC_DENSITY_INIT= 1.2 INC_VELOCITY_INIT= ( 1.0, 0.0, 0.0 ) INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 -INC_NONDIM= DIMENSIONAL +INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/multizone/zone2.cfg b/TestCases/species_transport/multizone/zone2.cfg index 2ad80c3dd58b..064243e245ed 100644 --- a/TestCases/species_transport/multizone/zone2.cfg +++ b/TestCases/species_transport/multizone/zone2.cfg @@ -31,7 +31,8 @@ INC_DENSITY_INIT= 1.2 INC_VELOCITY_INIT= ( 1.0, 0.0, 0.0 ) INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 -INC_NONDIM= DIMENSIONAL +INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/passive_transport_validation/active_species_transport_temp_limits.cfg b/TestCases/species_transport/passive_transport_validation/active_species_transport_temp_limits.cfg index 12dc4d7bb5d0..e1f7b5f897b3 100644 --- a/TestCases/species_transport/passive_transport_validation/active_species_transport_temp_limits.cfg +++ b/TestCases/species_transport/passive_transport_validation/active_species_transport_temp_limits.cfg @@ -28,6 +28,7 @@ TEMPERATURE_LIMITS= 250, 400 INC_TEMPERATURE_INIT= 350.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_JST.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_JST.cfg index 2bd04a386be3..2c097c77754a 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_JST.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_JST.cfg @@ -27,6 +27,7 @@ INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 350.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_Lax_Friedrich.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_Lax_Friedrich.cfg index dbe7fde8ebfc..df704e4d2f44 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_Lax_Friedrich.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_Lax_Friedrich.cfg @@ -28,6 +28,7 @@ INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 350.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel.cfg index d3fd801cf569..9355d4c53ffc 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel.cfg @@ -27,6 +27,7 @@ INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_TURBULENT_MARKERS.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_TURBULENT_MARKERS.cfg index 8dbb6ff6d47f..1248c4389684 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_TURBULENT_MARKERS.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_TURBULENT_MARKERS.cfg @@ -28,6 +28,7 @@ INC_ENERGY_EQUATION= NO INC_TEMPERATURE_INIT= 300.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_boundedscalar.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_boundedscalar.cfg index c80ec5088350..01284a7476a7 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_boundedscalar.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_boundedscalar.cfg @@ -27,6 +27,7 @@ INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg index 7d55656b64a4..7b488672442e 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg @@ -28,6 +28,7 @@ INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg index 7feb95f422e5..01a8750d3473 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg @@ -29,6 +29,7 @@ INC_ENERGY_EQUATION= YES INC_TEMPERATURE_INIT= 300.0 % INC_NONDIM= INITIAL_VALUES +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_viscosity.cfg b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_viscosity.cfg index a1117a029654..b9413340583b 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_viscosity.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel_viscosity.cfg @@ -28,6 +28,7 @@ INC_ENERGY_EQUATION= NO INC_TEMPERATURE_INIT= 300.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % diff --git a/TestCases/vandv/species_transport/sandia_jet/validation.cfg b/TestCases/vandv/species_transport/sandia_jet/validation.cfg index 42fd54500511..9b5a5c0b0f30 100644 --- a/TestCases/vandv/species_transport/sandia_jet/validation.cfg +++ b/TestCases/vandv/species_transport/sandia_jet/validation.cfg @@ -27,6 +27,7 @@ INC_ENERGY_EQUATION= NO INC_TEMPERATURE_INIT= 294.0 % INC_NONDIM= DIMENSIONAL +INIT_OPTION_INC= OPERATING_PRESSURE % % -------------------- FLUID PROPERTIES ------------------------------------- % % From 0df7d14e974f4cb04cf55be638cf724011f492b4 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Fri, 16 Jan 2026 09:12:07 +0100 Subject: [PATCH 2/4] checking tutorials cases --- .github/workflows/regression.yml | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/.github/workflows/regression.yml b/.github/workflows/regression.yml index c1bebd3928ff..3bd73816fca2 100644 --- a/.github/workflows/regression.yml +++ b/.github/workflows/regression.yml @@ -211,7 +211,7 @@ jobs: uses: docker://ghcr.io/su2code/su2/test-su2:250717-1402 with: # -t -c - args: -b ${{github.ref}} -t develop -c develop -s ${{matrix.testscript}} + args: -b ${{github.ref}} -t feature_initialization -c develop -s ${{matrix.testscript}} - name: Cleanup uses: docker://ghcr.io/su2code/su2/test-su2:250717-1402 with: From a27c73c0efb19ddd281bd2b956c9ff63a90df91c Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Fri, 16 Jan 2026 10:04:07 +0100 Subject: [PATCH 3/4] updating config_template --- config_template.cfg | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/config_template.cfg b/config_template.cfg index e519f5971aaf..061f6a825ab3 100644 --- a/config_template.cfg +++ b/config_template.cfg @@ -276,6 +276,10 @@ INC_TEMPERATURE_INIT= 288.15 % INC_*_REF values are ignored unless REFERENCE_VALUES is chosen. INC_NONDIM= INITIAL_VALUES % +% Init option incompressible solver to choose between initial density (default) or operating pressure +% for initializing the solution (DENSITY_INIT, OPERATING_PRESSURE) +INIT_OPTION_INC= DENSITY_INIT +% % Reference density for incompressible flows (1.0 kg/m^3 by default) INC_DENSITY_REF= 1.0 % From b3e5240041773c43e8d92715b8542d1e8a750fb3 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Tue, 20 Jan 2026 15:53:56 +0100 Subject: [PATCH 4/4] throw an error when INIT_OPTION_INC=OPERATING_PRESSURE is used with CONSTANT_DENSITY fluid model --- Common/src/CConfig.cpp | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/Common/src/CConfig.cpp b/Common/src/CConfig.cpp index a727d396881f..1c5c06185434 100644 --- a/Common/src/CConfig.cpp +++ b/Common/src/CConfig.cpp @@ -3562,6 +3562,12 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i SU2_MPI::Error("CoolProp can not be used with non-dimensionalization.", CURRENT_FUNCTION); } + /*--- Check if CONSTANT_DENSITY model is used with INIT_OPTION_INC=OPERATING_PRESSURE. ---*/ + if (Kind_FluidModel == CONSTANT_DENSITY && Kind_InitOption_Inc == OPERATING_PRESSURE) { + SU2_MPI::Error("CONSTANT_DENSITY fluid model can only be used with INIT_OPTION_INC=DENSITY_INIT.", + CURRENT_FUNCTION); + } + /*--- STL_BINARY output not implemented yet, but already a value in option_structure.hpp---*/ for (unsigned short iVolumeFile = 0; iVolumeFile < nVolumeOutputFiles; iVolumeFile++) { if (VolumeOutputFiles[iVolumeFile] == OUTPUT_TYPE::STL_BINARY){