su2code · Cristopher-Morales · Sep 15, 2023 · Sep 18, 2023 · Sep 20, 2023 · Sep 20, 2023
diff --git a/Common/include/CConfig.hpp b/Common/include/CConfig.hpp
@@ -821,6 +821,7 @@ class CConfig {
   su2double Gamma,      /*!< \brief Ratio of specific heats of the gas. */
   Bulk_Modulus,         /*!< \brief Value of the bulk modulus for incompressible flows. */
   Beta_Factor,          /*!< \brief Value of the epsilon^2 multiplier for Beta for the incompressible preconditioner. */
+  Density_Constant,     /*!< \brief Constant Density for ConstantDensity model (only for incompressible flows). */
   Gas_Constant,         /*!< \brief Specific gas constant. */
   Gas_ConstantND,       /*!< \brief Non-dimensional specific gas constant. */
   *Molecular_Weight;    /*!< \brief Molecular weight of an incompressible ideal gas (g/mol). */
@@ -831,8 +832,7 @@ class CConfig {
   Thermal_Expansion_CoeffND,  /*!< \brief Non-dimensional thermal expansion coefficient. */
   Inc_Density_Ref,       /*!< \brief Reference density for custom incompressible non-dim. */
   Inc_Velocity_Ref,      /*!< \brief Reference velocity for custom incompressible non-dim. */
-  Inc_Temperature_Ref,   /*!< \brief Reference temperature for custom incompressible non-dim. */
-  Inc_Density_Init,      /*!< \brief Initial density for incompressible flows. */
+  Inc_Temperature_Ref,   /*!< \brief Reference temperature for custom incompressible non-dim. */ 
   Inc_Temperature_Init,  /*!< \brief Initial temperature for incompressible flows w/ heat transfer. */
   Heat_Flux_Ref,         /*!< \brief Reference heat flux for non-dim. */
   Gas_Constant_Ref,      /*!< \brief Reference specific gas constant. */
@@ -1668,6 +1668,12 @@ class CConfig {
    */
   su2double GetMolecular_Weight(unsigned short val_index = 0) const { return Molecular_Weight[val_index]; }
 
+   /*!
+   * \brief Get the value of the constant density for incompressible flows.
+   * \return Constant density for incompressible flows when CONSTANT_DENSITY fluid model is used.
+   */
+  su2double GetInc_Constant_Density(void) const { return Density_Constant; }
+
   /*!
    * \brief Get the value of specific heat at constant pressure.
    * \return Value of the constant: Cp
@@ -2080,12 +2086,6 @@ class CConfig {
    */
   su2double GetInc_Temperature_Ref(void) const { return Inc_Temperature_Ref; }
 
-  /*!
-   * \brief Get the value of the initial density for incompressible flows.
-   * \return Initial density for incompressible flows.
-   */
-  su2double GetInc_Density_Init(void) const { return Inc_Density_Init; }
-
   /*!
    * \brief Get the value of the initial velocity for incompressible flows.
    * \return Initial velocity for incompressible flows.

diff --git a/Common/src/CConfig.cpp b/Common/src/CConfig.cpp
@@ -1246,6 +1246,9 @@ void CConfig::SetConfig_Options() {
   /* DESCRIPTION: default value for AIR */
   addDoubleListOption("MU_CONSTANT", nMu_Constant, Mu_Constant);
 
+  /*!\brief CONSTANT_DENSITY \n DESCRIPTION: Constant density for incompressible flows (1.2886 kg/m^3 by default) \ingroup Config*/
+  addDoubleOption("DENSITY_CONSTANT", Density_Constant, 1.2886);
+
   /*--- Options related to Sutherland Viscosity Model ---*/
 
   /* DESCRIPTION: Sutherland Viscosity Ref default value for AIR SI */
@@ -1334,8 +1337,6 @@ void CConfig::SetConfig_Options() {
   addDoubleOption("INC_VELOCITY_REF", Inc_Velocity_Ref, 1.0);
   /*!\brief INC_TEMPERATURE_REF \n DESCRIPTION: Reference temperature for incompressible flows with the energy equation (1.0 by default) \ingroup Config*/
   addDoubleOption("INC_TEMPERATURE_REF", Inc_Temperature_Ref, 1.0);
-  /*!\brief INC_DENSITY_INIT \n DESCRIPTION: Initial density for incompressible flows (1.2886 kg/m^3 by default) \ingroup Config*/
-  addDoubleOption("INC_DENSITY_INIT", Inc_Density_Init, 1.2886);
   /*!\brief INC_VELOCITY_INIT \n DESCRIPTION: Initial velocity for incompressible flows (1.0,0,0 m/s by default) \ingroup Config*/
   vel_init[0] = 1.0; vel_init[1] = 0.0; vel_init[2] = 0.0;
   addDoubleArrayOption("INC_VELOCITY_INIT", 3, vel_init);
@@ -3022,6 +3023,11 @@ void CConfig::SetConfig_Parsing(istream& config_buffer){
             newString.append("DYNAMIC_ANALYSIS is deprecated. Use TIME_DOMAIN instead.\n\n");
           else if (!option_name.compare("SPECIES_USE_STRONG_BC"))
             newString.append("SPECIES_USE_STRONG_BC is deprecated. Use MARKER_SPECIES_STRONG_BC= (marker1, ...) instead.\n\n");
+          else if (!option_name.compare("INC_DENSITY_INIT"))
+            newString.append(
+                "INC_DENSITY_INIT is deprecated. The initial density is computed using the ideal gas law based on the "
+                "THERMODYNAMICS_PRESSURE and INC_TEMPERATURE_INIT. For Using the FLUID_MODEL=CONSTANT_DENSITY, the "
+                "constant density must be given using the following option DENSITY_CONSTANT.\n\n");
           else {
             /*--- Find the most likely candidate for the unrecognized option, based on the length
              of start and end character sequences shared by candidates and the option. ---*/

diff --git a/SU2_CFD/src/solvers/CIncEulerSolver.cpp b/SU2_CFD/src/solvers/CIncEulerSolver.cpp
@@ -265,7 +265,6 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i
 
   /*--- Compute dimensional free-stream values. ---*/
 
-  Density_FreeStream     = config->GetInc_Density_Init();     config->SetDensity_FreeStream(Density_FreeStream);
   Temperature_FreeStream = config->GetInc_Temperature_Init(); config->SetTemperature_FreeStream(Temperature_FreeStream);
   Pressure_FreeStream    = 0.0; config->SetPressure_FreeStream(Pressure_FreeStream);
 
@@ -291,24 +290,24 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i
 
     case CONSTANT_DENSITY:
 
-      auxFluidModel = new CConstantDensity(Density_FreeStream, config->GetSpecific_Heat_Cp());
+      auxFluidModel = new CConstantDensity(config->GetInc_Constant_Density(), config->GetSpecific_Heat_Cp());
       auxFluidModel->SetTDState_T(Temperature_FreeStream);
+      Density_FreeStream= auxFluidModel->GetDensity(); config->SetDensity_FreeStream(Density_FreeStream);
       break;
 
     case INC_IDEAL_GAS:
 
       config->SetGas_Constant(UNIVERSAL_GAS_CONSTANT/(config->GetMolecular_Weight()/1000.0));
-      Pressure_Thermodynamic = Density_FreeStream*Temperature_FreeStream*config->GetGas_Constant();
+      Pressure_Thermodynamic = config->GetPressure_Thermodynamic();
       auxFluidModel = new CIncIdealGas(config->GetSpecific_Heat_Cp(), config->GetGas_Constant(), Pressure_Thermodynamic);
       auxFluidModel->SetTDState_T(Temperature_FreeStream);
-      Pressure_Thermodynamic = auxFluidModel->GetPressure();
-      config->SetPressure_Thermodynamic(Pressure_Thermodynamic);
+      Density_FreeStream= auxFluidModel->GetDensity(); config->SetDensity_FreeStream(Density_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();
+      Pressure_Thermodynamic = config->GetPressure_Thermodynamic();
       auxFluidModel = new CIncIdealGasPolynomial<N_POLY_COEFFS>(config->GetGas_Constant(), Pressure_Thermodynamic);
       if (viscous) {
         /*--- Variable Cp model via polynomial. ---*/
@@ -317,8 +316,7 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i
         auxFluidModel->SetCpModel(config);
       }
       auxFluidModel->SetTDState_T(Temperature_FreeStream);
-      Pressure_Thermodynamic = auxFluidModel->GetPressure();
-      config->SetPressure_Thermodynamic(Pressure_Thermodynamic);
+      Density_FreeStream= auxFluidModel->GetDensity(); config->SetDensity_FreeStream(Density_FreeStream);
       break;
 
     case FLUID_MIXTURE:
@@ -327,15 +325,16 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i
       Pressure_Thermodynamic = config->GetPressure_Thermodynamic();
       auxFluidModel = new CFluidScalar(config->GetSpecific_Heat_Cp(), config->GetGas_Constant(), Pressure_Thermodynamic, config);
       auxFluidModel->SetTDState_T(Temperature_FreeStream, config->GetSpecies_Init());
+      Density_FreeStream= auxFluidModel->GetDensity(); config->SetDensity_FreeStream(Density_FreeStream);
       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());
+      Density_FreeStream= auxFluidModel->GetDensity(); config->SetDensity_FreeStream(Density_FreeStream);
       break;
 
     default:

diff --git a/TestCases/coupled_cht/disc_adj_incomp_2d/flow_cylinder.cfg b/TestCases/coupled_cht/disc_adj_incomp_2d/flow_cylinder.cfg
@@ -27,7 +27,6 @@ MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
 %
 INC_DENSITY_MODEL= VARIABLE
 INC_ENERGY_EQUATION = YES
-INC_DENSITY_INIT= 0.000210322
 INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
 INC_NONDIM= DIMENSIONAL
@@ -37,6 +36,7 @@ INC_NONDIM= DIMENSIONAL
 FLUID_MODEL= INC_IDEAL_GAS
 SPECIFIC_HEAT_CP= 1004.703
 MOLECULAR_WEIGHT= 28.96
+THERMODYNAMIC_PRESSURE = 17.39861946
 
 % --------------------------- VISCOSITY MODEL ---------------------------------%
 %

diff --git a/TestCases/coupled_cht/disc_adj_unsteadyCHT_cylinder/fluid.cfg b/TestCases/coupled_cht/disc_adj_unsteadyCHT_cylinder/fluid.cfg
@@ -10,7 +10,7 @@ OBJECTIVE_WEIGHT= 0.0
 INC_NONDIM= DIMENSIONAL
 % ~Water at room temperature
 INC_DENSITY_MODEL= CONSTANT
-INC_DENSITY_INIT= 1.0e3
+DENSITY_CONSTANT = 1.0e3
 %
 % Re = rho*v*d/mu_dyn = 1e3*v*1/1e-3 = 200 => v = 2.0e-4
 INC_VELOCITY_INIT= ( 2.0e-4, 0.0, 0.0 )

diff --git a/TestCases/coupled_cht/incomp_2d/flow_cylinder.cfg b/TestCases/coupled_cht/incomp_2d/flow_cylinder.cfg
@@ -25,7 +25,6 @@ MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
 %
 INC_DENSITY_MODEL= VARIABLE
 INC_ENERGY_EQUATION = YES
-INC_DENSITY_INIT= 0.000210322
 INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
 INC_NONDIM= DIMENSIONAL
@@ -35,6 +34,7 @@ INC_NONDIM= DIMENSIONAL
 FLUID_MODEL= INC_IDEAL_GAS
 SPECIFIC_HEAT_CP= 1004.703
 MOLECULAR_WEIGHT= 28.96
+THERMODYNAMIC_PRESSURE = 17.39861946
 
 % --------------------------- VISCOSITY MODEL ---------------------------------%
 %

diff --git a/TestCases/coupled_cht/incomp_2d_unsteady/flow_cylinder.cfg b/TestCases/coupled_cht/incomp_2d_unsteady/flow_cylinder.cfg
@@ -25,7 +25,6 @@ MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
 %
 INC_DENSITY_MODEL= VARIABLE
 INC_ENERGY_EQUATION = YES
-INC_DENSITY_INIT= 0.0210322
 INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
 INC_NONDIM= DIMENSIONAL
@@ -35,6 +34,7 @@ INC_NONDIM= DIMENSIONAL
 FLUID_MODEL= INC_IDEAL_GAS
 SPECIFIC_HEAT_CP= 1004.703
 MOLECULAR_WEIGHT= 28.96
+THERMODYNAMIC_PRESSURE = 1739.861946
 
 % --------------------------- VISCOSITY MODEL ---------------------------------%
 %

diff --git a/TestCases/deformation/intersection_prevention/def_intersect.cfg b/TestCases/deformation/intersection_prevention/def_intersect.cfg
@@ -29,8 +29,8 @@ AXISYMMETRIC = NO
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
 INC_DENSITY_MODEL= CONSTANT
+DENSITY_CONSTANT= 1.1728
 INC_ENERGY_EQUATION = NO
-INC_DENSITY_INIT= 1.1728
 INC_VELOCITY_INIT= (1.0, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 300.0
 INC_INLET_TYPE= VELOCITY_INLET

diff --git a/TestCases/disc_adj_heat/disc_adj_heat.cfg b/TestCases/disc_adj_heat/disc_adj_heat.cfg
@@ -33,9 +33,9 @@ ITER = 11
 
 INC_DENSITY_MODEL= CONSTANT
 
-INC_ENERGY_EQUATION= NO
+DENSITY_CONSTANT= 998.2
 
-INC_DENSITY_INIT= 998.2
+INC_ENERGY_EQUATION= NO
 
 INC_VELOCITY_INIT= ( 0.25, 0.0, 0.0 )
 

diff --git a/TestCases/disc_adj_heat/primal.cfg b/TestCases/disc_adj_heat/primal.cfg
@@ -20,8 +20,8 @@ OBJECTIVE_FUNCTION= TOTAL_HEATFLUX
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 
 INC_DENSITY_MODEL= CONSTANT
+DENSITY_CONSTANT= 998.2
 INC_ENERGY_EQUATION= NO
-INC_DENSITY_INIT= 998.2
 INC_VELOCITY_INIT= ( 0.25, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 300.0
 INC_NONDIM= INITIAL_VALUES

diff --git a/TestCases/disc_adj_incomp_euler/naca0012/incomp_NACA0012_disc.cfg b/TestCases/disc_adj_incomp_euler/naca0012/incomp_NACA0012_disc.cfg
@@ -17,7 +17,7 @@ READ_BINARY_RESTART= NO
 
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
-INC_DENSITY_INIT= 1.2886
+DENSITY_CONSTANT= 1.2886
 INC_VELOCITY_INIT= ( 33.1588, 1.15793, 0.0 )
 
 % ---------------------- REFERENCE VALUE DEFINITION ---------------------------%

diff --git a/TestCases/disc_adj_incomp_navierstokes/cylinder/heated_cylinder.cfg b/TestCases/disc_adj_incomp_navierstokes/cylinder/heated_cylinder.cfg
@@ -19,7 +19,6 @@ RESTART_SOL= NO
 %
 INC_DENSITY_MODEL= VARIABLE
 INC_ENERGY_EQUATION = YES
-INC_DENSITY_INIT= 0.000210322
 INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
 
@@ -28,6 +27,7 @@ INC_TEMPERATURE_INIT= 288.15
 FLUID_MODEL= INC_IDEAL_GAS
 SPECIFIC_HEAT_CP= 1004.703
 MOLECULAR_WEIGHT= 28.96
+THERMODYNAMIC_PRESSURE = 17.39861946
 
 % --------------------------- VISCOSITY MODEL ---------------------------------%
 %

diff --git a/TestCases/disc_adj_incomp_rans/naca0012/turb_naca0012_sa.cfg b/TestCases/disc_adj_incomp_rans/naca0012/turb_naca0012_sa.cfg
@@ -18,7 +18,7 @@ RESTART_SOL= NO
 
 % -------------------- INCOMPRESSIBLE FREE-STREAM DEFINITION ------------------%
 %
-INC_DENSITY_INIT= 2.13163
+DENSITY_CONSTANT= 2.13163
 % AoA 10 deg
 INC_VELOCITY_INIT= ( 51.36481493540834, 9.0570027322096198, 0.0 )
 

diff --git a/TestCases/disc_adj_incomp_rans/naca0012/turb_naca0012_sst.cfg b/TestCases/disc_adj_incomp_rans/naca0012/turb_naca0012_sst.cfg
@@ -17,7 +17,7 @@ RESTART_SOL= NO
 
 % -------------------- INCOMPRESSIBLE FREE-STREAM DEFINITION ------------------%
 %
-INC_DENSITY_INIT= 2.13163
+DENSITY_CONSTANT= 2.13163
 % AoA 10 deg
 INC_VELOCITY_INIT= ( 51.36481493540834, 9.0570027322096198, 0.0 )
 

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
@@ -19,7 +19,6 @@ RESTART_SOL = YES
 %
 INC_DENSITY_MODEL= FLAMELET
 INC_ENERGY_EQUATION = YES
-INC_DENSITY_INIT= 1.00
 INC_VELOCITY_INIT= (0.5, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 300.0
 INC_NONDIM= DIMENSIONAL

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
@@ -28,7 +28,6 @@ CUSTOM_OUTPUTS= 'avg_NOx : MassFlowAvg{SPECIES[3]}[outlet]'
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
 INC_DENSITY_MODEL= FLAMELET
-INC_DENSITY_INIT= 1.1766
 INC_VELOCITY_INIT= ( 0.2, 0, 0 )
 INC_ENERGY_EQUATION= YES
 INC_TEMPERATURE_INIT= 300.0

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
@@ -29,7 +29,6 @@ CUSTOM_OUTPUTS= 'avg_NOx : MassFlowAvg{SPECIES[3]}[outlet]'
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
 INC_DENSITY_MODEL= FLAMELET
-INC_DENSITY_INIT= 1.1766
 INC_VELOCITY_INIT= ( 0.2, 0, 0 )
 INC_ENERGY_EQUATION= YES
 INC_TEMPERATURE_INIT= 300.0

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
@@ -43,7 +43,6 @@ MARKER_SPECIES_STRONG_BC= (inlet, outlet, wall_burner)
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 
 INC_DENSITY_MODEL= FLAMELET
-INC_DENSITY_INIT= 1.1766
 INC_VELOCITY_INIT= ( 0.2, 0, 0 )
 INC_ENERGY_EQUATION= YES
 INC_TEMPERATURE_INIT= 300.0

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
@@ -45,7 +45,6 @@ MARKER_SPECIES_STRONG_BC= (inlet, wall_burner)
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 
 INC_DENSITY_MODEL= FLAMELET
-INC_DENSITY_INIT= 1.1766
 INC_VELOCITY_INIT= ( 0.2, 0, 0 )
 INC_ENERGY_EQUATION= YES
 INC_TEMPERATURE_INIT= 300.0

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
@@ -20,7 +20,6 @@ RESTART_SOL= YES
 %
 INC_DENSITY_MODEL= FLAMELET
 INC_ENERGY_EQUATION= YES
-INC_DENSITY_INIT= 1.00
 INC_VELOCITY_INIT= (0.5, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 300.0
 INC_NONDIM= DIMENSIONAL

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
@@ -19,7 +19,6 @@ RESTART_SOL= YES
 %
 INC_DENSITY_MODEL= VARIABLE
 INC_ENERGY_EQUATION = YES
-INC_DENSITY_INIT= 1.00
 INC_VELOCITY_INIT= (0.7, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 300.0
 INC_NONDIM= DIMENSIONAL

diff --git a/TestCases/incomp_euler/naca0012/incomp_NACA0012.cfg b/TestCases/incomp_euler/naca0012/incomp_NACA0012.cfg
@@ -17,7 +17,7 @@ RESTART_SOL= NO
 
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
-INC_DENSITY_INIT= 998.2
+DENSITY_CONSTANT= 998.2
 INC_VELOCITY_INIT= ( 1.775, 0.0, 0.0 )
 INC_INLET_TYPE= VELOCITY_INLET
 INC_INLET_DAMPING= 0.1

diff --git a/TestCases/incomp_euler/nozzle/inv_nozzle.cfg b/TestCases/incomp_euler/nozzle/inv_nozzle.cfg
@@ -17,8 +17,8 @@ RESTART_SOL= NO
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
 INC_DENSITY_MODEL= CONSTANT
+DENSITY_CONSTANT= 1.0
 INC_ENERGY_EQUATION = NO
-INC_DENSITY_INIT= 1.0
 INC_VELOCITY_INIT= ( 1.0, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
 INC_INLET_TYPE= PRESSURE_INLET

diff --git a/TestCases/incomp_navierstokes/bend/lam_bend.cfg b/TestCases/incomp_navierstokes/bend/lam_bend.cfg
@@ -19,8 +19,8 @@ RESTART_SOL= NO
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
 INC_DENSITY_MODEL= CONSTANT
+DENSITY_CONSTANT= 1.2886
 INC_ENERGY_EQUATION = NO
-INC_DENSITY_INIT= 1.2886
 INC_VELOCITY_INIT= ( 0.1, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
 INC_NONDIM= INITIAL_VALUES

diff --git a/TestCases/incomp_navierstokes/buoyancy_cavity/lam_buoyancy_cavity.cfg b/TestCases/incomp_navierstokes/buoyancy_cavity/lam_buoyancy_cavity.cfg
@@ -24,7 +24,7 @@ INC_ENERGY_EQUATION = YES
 %INC_DENSITY_INIT= 0.00018903539834  % Ra ~ 1e3
 %INC_DENSITY_INIT= 0.00059778241716  % Ra ~ 1e4
 %INC_DENSITY_INIT= 0.00189035398341  % Ra ~ 1e5
-INC_DENSITY_INIT= 0.00597782417156  % Ra ~ 1e6
+%INC_DENSITY_INIT= 0.00597782417156  % Ra ~ 1e6
 %
 INC_VELOCITY_INIT= ( 1.0, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
@@ -34,6 +34,7 @@ INC_TEMPERATURE_INIT= 288.15
 FLUID_MODEL= INC_IDEAL_GAS
 SPECIFIC_HEAT_CP= 1004.703
 MOLECULAR_WEIGHT= 28.96
+THERMODYNAMIC_PRESSURE = 494.5078878
 
 % --------------------------- VISCOSITY MODEL ---------------------------------%
 %

diff --git a/TestCases/incomp_navierstokes/cylinder/incomp_cylinder.cfg b/TestCases/incomp_navierstokes/cylinder/incomp_cylinder.cfg
@@ -18,7 +18,7 @@ RESTART_SOL= NO
 
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
-INC_DENSITY_INIT= 998.2
+DENSITY_CONSTANT= 998.2
 INC_VELOCITY_INIT= ( 0.000008, 0.0, 0.0 )
 
 % --------------------------- VISCOSITY MODEL ---------------------------------%

diff --git a/TestCases/incomp_navierstokes/cylinder/poly_cylinder.cfg b/TestCases/incomp_navierstokes/cylinder/poly_cylinder.cfg
@@ -20,7 +20,6 @@ RESTART_SOL= NO
 %
 INC_DENSITY_MODEL= VARIABLE
 INC_ENERGY_EQUATION = YES
-INC_DENSITY_INIT= 0.000210322
 INC_VELOCITY_INIT= ( 3.40297, 0.0, 0.0 )
 INC_TEMPERATURE_INIT= 288.15
 INC_NONDIM= DIMENSIONAL
@@ -31,6 +30,7 @@ FLUID_MODEL= INC_IDEAL_GAS_POLY
 SPECIFIC_HEAT_CP= 1004.703
 MOLECULAR_WEIGHT= 28.96
 CP_POLYCOEFFS= ( 1004.703, 0.1, 0.0, 0.0, 0.0)
+THERMODYNAMIC_PRESSURE = 17.39861946
 
 % --------------------------- VISCOSITY MODEL ---------------------------------%
 %