From c2d86ea20886fec616bb20fd453be50cccfba2dd Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 1 Dec 2022 10:50:09 +0100 Subject: [PATCH 1/7] adding thermodynamics pressure option --- Common/src/CConfig.cpp | 2 ++ SU2_CFD/include/fluid/CFluidScalar.hpp | 1 - SU2_CFD/src/fluid/CFluidScalar.cpp | 5 ++--- SU2_CFD/src/solvers/CIncEulerSolver.cpp | 3 +-- 4 files changed, 5 insertions(+), 6 deletions(-) diff --git a/Common/src/CConfig.cpp b/Common/src/CConfig.cpp index c5586098980..7514aff4f7e 100644 --- a/Common/src/CConfig.cpp +++ b/Common/src/CConfig.cpp @@ -1161,6 +1161,8 @@ void CConfig::SetConfig_Options() { addDoubleOption("GAS_CONSTANT", Gas_Constant, 287.058); /*!\brief GAMMA_VALUE \n DESCRIPTION: Ratio of specific heats (1.4 (air), only for compressible flows) \ingroup Config*/ addDoubleOption("GAMMA_VALUE", Gamma, 1.4); + /*!\brief THERMODYNAMIC_PRESSURE \n DESCRIPTION: Thermodynamics(operating) Pressure (101325 Pa), only for incompressible flows) \ingroup Config*/ + addDoubleOption("THERMODYNAMIC_PRESSURE", Pressure_Thermodynamic, 101325.0); /*!\brief CP_VALUE \n DESCRIPTION: Specific heat at constant pressure, Cp (1004.703 J/kg*K (air), constant density incompressible fluids only) \ingroup Config*/ addDoubleListOption("SPECIFIC_HEAT_CP", nSpecific_Heat_Cp, Specific_Heat_Cp); /*!\brief THERMAL_EXPANSION_COEFF \n DESCRIPTION: Thermal expansion coefficient (0.00347 K^-1 (air), used for Boussinesq approximation for liquids/non-ideal gases) \ingroup Config*/ diff --git a/SU2_CFD/include/fluid/CFluidScalar.hpp b/SU2_CFD/include/fluid/CFluidScalar.hpp index 619b869c757..354e0a71910 100644 --- a/SU2_CFD/include/fluid/CFluidScalar.hpp +++ b/SU2_CFD/include/fluid/CFluidScalar.hpp @@ -42,7 +42,6 @@ class CFluidScalar final : public CFluidModel { const int n_species_mixture; /*!< \brief Number of species in mixture. */ su2double Gas_Constant; /*!< \brief Specific gas constant. */ const su2double Gamma; /*!< \brief Ratio of specific heats of the gas. */ - const su2double Pressure_Thermodynamic; /*!< \brief Constant pressure thermodynamic. */ const su2double GasConstant_Ref; /*!< \brief Gas constant reference needed for Nondimensional problems. */ const bool wilke; diff --git a/SU2_CFD/src/fluid/CFluidScalar.cpp b/SU2_CFD/src/fluid/CFluidScalar.cpp index eb8f4624aec..9c1ddc6006b 100644 --- a/SU2_CFD/src/fluid/CFluidScalar.cpp +++ b/SU2_CFD/src/fluid/CFluidScalar.cpp @@ -48,14 +48,13 @@ CFluidScalar::CFluidScalar(su2double val_Cp, su2double val_gas_constant, su2doub n_species_mixture(config->GetnSpecies() + 1), Gas_Constant(val_gas_constant), Gamma(config->GetGamma()), - Pressure_Thermodynamic(value_pressure_operating), GasConstant_Ref(config->GetGas_Constant_Ref()), wilke(config->GetKind_MixingViscosityModel() == MIXINGVISCOSITYMODEL::WILKE), davidson(config->GetKind_MixingViscosityModel() == MIXINGVISCOSITYMODEL::DAVIDSON) { if (n_species_mixture > ARRAYSIZE) { SU2_MPI::Error("Too many species, increase ARRAYSIZE", CURRENT_FUNCTION); } - + Pressure = value_pressure_operating; for (int iVar = 0; iVar < n_species_mixture; iVar++) { molarMasses[iVar] = config->GetMolecular_Weight(iVar); specificHeat[iVar] = config->GetSpecific_Heat_CpND(iVar); @@ -216,7 +215,7 @@ void CFluidScalar::SetTDState_T(const su2double val_temperature, const su2double MassToMoleFractions(val_scalars); ComputeGasConstant(); Temperature = val_temperature; - Density = Pressure_Thermodynamic / (Temperature * Gas_Constant); + Density = Pressure / (Temperature * Gas_Constant); Cp = ComputeMeanSpecificHeatCp(val_scalars); Cv = Cp - Gas_Constant; diff --git a/SU2_CFD/src/solvers/CIncEulerSolver.cpp b/SU2_CFD/src/solvers/CIncEulerSolver.cpp index e7f238e7334..a85c66fa11b 100644 --- a/SU2_CFD/src/solvers/CIncEulerSolver.cpp +++ b/SU2_CFD/src/solvers/CIncEulerSolver.cpp @@ -328,10 +328,9 @@ void CIncEulerSolver::SetNondimensionalization(CConfig *config, unsigned short i case FLUID_MIXTURE: 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 CFluidScalar(config->GetSpecific_Heat_Cp(), config->GetGas_Constant(), Pressure_Thermodynamic, config); auxFluidModel->SetTDState_T(Temperature_FreeStream, config->GetSpecies_Init()); - config->SetPressure_Thermodynamic(Pressure_Thermodynamic); break; default: From 39756ec8a56839d69f03b55ea251451c9414ef93 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 1 Dec 2022 12:38:32 +0100 Subject: [PATCH 2/7] updating residuals --- TestCases/parallel_regression.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/TestCases/parallel_regression.py b/TestCases/parallel_regression.py index 75833a37f67..aa36c35788c 100644 --- a/TestCases/parallel_regression.py +++ b/TestCases/parallel_regression.py @@ -1338,7 +1338,7 @@ def main(): species2_primitiveVenturi_mixingmodel.cfg_dir = "species_transport/venturi_primitive_3species" species2_primitiveVenturi_mixingmodel.cfg_file = "species2_primitiveVenturi_mixingmodel.cfg" species2_primitiveVenturi_mixingmodel.test_iter = 50 - 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] + species2_primitiveVenturi_mixingmodel.test_vals = [-5.477236, -4.589782, -4.582531, -5.690416, -0.068943, -5.631254, 5.000000, -1.668303, 5.000000, -4.866978, 5.000000, -1.168368, 0.000425, 0.000404, 0.000021, 0.000000] species2_primitiveVenturi_mixingmodel.new_output = True test_list.append(species2_primitiveVenturi_mixingmodel) @@ -1347,7 +1347,7 @@ def main(): species2_primitiveVenturi_mixingmodel_boundedscalar.cfg_dir = "species_transport/venturi_primitive_3species" species2_primitiveVenturi_mixingmodel_boundedscalar.cfg_file = "species2_primitiveVenturi_mixingmodel_boundedscalar.cfg" species2_primitiveVenturi_mixingmodel_boundedscalar.test_iter = 50 - species2_primitiveVenturi_mixingmodel_boundedscalar.test_vals = [-5.419758, -4.490843, -4.496264, -5.724852, -0.120393, -5.693152, 5.000000, -1.766881, 5.000000, -4.958352, 5.000000, -2.150266, 0.000300, 0.000300, 0.000000, 0.000000] + species2_primitiveVenturi_mixingmodel_boundedscalar.test_vals = [-5.419823, -4.490921, -4.496340, -5.724950, -0.120331, -5.693224, 5.000000, -1.766882, 5.000000, -4.958487, 5.000000, -2.150247, 0.000300, 0.000300, 0.000000, 0.000000] species2_primitiveVenturi_mixingmodel_boundedscalar.new_output = True test_list.append(species2_primitiveVenturi_mixingmodel_boundedscalar) @@ -1365,7 +1365,7 @@ def main(): species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg_dir = "species_transport/venturi_primitive_3species" species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg_file = "species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg" species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_iter = 50 - 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] + species2_primitiveVenturi_mixingmodel_heatcapacity_H2.test_vals = [-6.105905, -4.979417, -4.870704, -7.338344, 2.435662, -5.621107, 30.000000, -5.737172, 12.000000, -8.148164, 9.000000, -8.040415, 2.084589, 1.000000, 0.600000, 0.484589] species2_primitiveVenturi_mixingmodel_heatcapacity_H2.su2_exec = "mpirun -n 2 SU2_CFD" species2_primitiveVenturi_mixingmodel_heatcapacity_H2.timeout = 1600 species2_primitiveVenturi_mixingmodel_heatcapacity_H2.new_output = True @@ -1377,7 +1377,7 @@ def main(): species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg_dir = "species_transport/venturi_primitive_3species" species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg_file = "species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg" species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_iter = 50 - 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] + species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.test_vals = [-5.711760, -5.284237, -5.175532, -8.342082, 2.130839, -5.226944, 30.000000, -5.737157, 12.000000, -8.148020, 9.000000, -8.040460, 2.084591, 1.000000, 0.600000, 0.484591] species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.su2_exec = "mpirun -n 2 SU2_CFD" species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.timeout = 1600 species2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.new_output = True @@ -1427,7 +1427,7 @@ def main(): species3_multizone_restart.cfg_dir = "species_transport/multizone" species3_multizone_restart.cfg_file = "configMaster.cfg" species3_multizone_restart.test_iter = 5 - species3_multizone_restart.test_vals = [-4.773393, -4.254100] + species3_multizone_restart.test_vals = [-2.439389, -2.284570] species3_multizone_restart.new_output = True species3_multizone_restart.multizone = True test_list.append(species3_multizone_restart) From 839ee90073a50d4777061ca79d95d8c1ad2cff6e Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 1 Dec 2022 13:28:51 +0100 Subject: [PATCH 3/7] updating residual species mixingmodel_viscosity --- TestCases/parallel_regression.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/TestCases/parallel_regression.py b/TestCases/parallel_regression.py index aa36c35788c..cabc4e75580 100644 --- a/TestCases/parallel_regression.py +++ b/TestCases/parallel_regression.py @@ -1356,7 +1356,7 @@ def main(): species2_primitiveVenturi_mixingmodel_viscosity.cfg_dir = "species_transport/venturi_primitive_3species" species2_primitiveVenturi_mixingmodel_viscosity.cfg_file = "species2_primitiveVenturi_mixingmodel_viscosity.cfg" species2_primitiveVenturi_mixingmodel_viscosity.test_iter = 50 - 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] + species2_primitiveVenturi_mixingmodel_viscosity.test_vals = [-4.738778, -4.325766, -4.610914, -5.834431, 0.521730, -4.934375, 5.000000, -1.887191, 5.000000, -5.499917, 5.000000, -1.770845, 2.292904, 0.971941, 0.608500, 0.712464] species2_primitiveVenturi_mixingmodel_viscosity.new_output = True test_list.append(species2_primitiveVenturi_mixingmodel_viscosity) From 3221a21d9c9d117e877628be9439a446c5c4316c Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 1 Dec 2022 14:04:25 +0100 Subject: [PATCH 4/7] updating restart test case --- .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 8e2d4414353..5b235e4a6a3 100644 --- a/.github/workflows/regression.yml +++ b/.github/workflows/regression.yml @@ -128,7 +128,7 @@ jobs: uses: docker://ghcr.io/su2code/su2/test-su2:220614-1237 with: # -t -c - args: -b ${{github.ref}} -t develop -c develop -s ${{matrix.testscript}} + args: -b ${{github.ref}} -t develop -c feature_update_restart_multizone -s ${{matrix.testscript}} - name: Cleanup uses: docker://ghcr.io/su2code/su2/test-su2:220614-1237 with: From 9ef9710e4743c90a4561bd9f5e72b427e9ce3384 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 1 Dec 2022 14:50:45 +0100 Subject: [PATCH 5/7] updating residuals restart --- TestCases/parallel_regression.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/TestCases/parallel_regression.py b/TestCases/parallel_regression.py index cabc4e75580..888f91a2912 100644 --- a/TestCases/parallel_regression.py +++ b/TestCases/parallel_regression.py @@ -1427,7 +1427,7 @@ def main(): species3_multizone_restart.cfg_dir = "species_transport/multizone" species3_multizone_restart.cfg_file = "configMaster.cfg" species3_multizone_restart.test_iter = 5 - species3_multizone_restart.test_vals = [-2.439389, -2.284570] + species3_multizone_restart.test_vals = [-7.169637, -6.381576] species3_multizone_restart.new_output = True species3_multizone_restart.multizone = True test_list.append(species3_multizone_restart) From 85c08571f68fa4f064c77e92811479d964f76bb5 Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Sun, 4 Dec 2022 18:05:50 +0100 Subject: [PATCH 6/7] addressing comments --- .github/workflows/regression.yml | 2 +- SU2_CFD/include/fluid/CFluidScalar.hpp | 1 + SU2_CFD/src/fluid/CFluidScalar.cpp | 5 +++-- 3 files changed, 5 insertions(+), 3 deletions(-) diff --git a/.github/workflows/regression.yml b/.github/workflows/regression.yml index 5b235e4a6a3..8e2d4414353 100644 --- a/.github/workflows/regression.yml +++ b/.github/workflows/regression.yml @@ -128,7 +128,7 @@ jobs: uses: docker://ghcr.io/su2code/su2/test-su2:220614-1237 with: # -t -c - args: -b ${{github.ref}} -t develop -c feature_update_restart_multizone -s ${{matrix.testscript}} + args: -b ${{github.ref}} -t develop -c develop -s ${{matrix.testscript}} - name: Cleanup uses: docker://ghcr.io/su2code/su2/test-su2:220614-1237 with: diff --git a/SU2_CFD/include/fluid/CFluidScalar.hpp b/SU2_CFD/include/fluid/CFluidScalar.hpp index 354e0a71910..619b869c757 100644 --- a/SU2_CFD/include/fluid/CFluidScalar.hpp +++ b/SU2_CFD/include/fluid/CFluidScalar.hpp @@ -42,6 +42,7 @@ class CFluidScalar final : public CFluidModel { const int n_species_mixture; /*!< \brief Number of species in mixture. */ su2double Gas_Constant; /*!< \brief Specific gas constant. */ const su2double Gamma; /*!< \brief Ratio of specific heats of the gas. */ + const su2double Pressure_Thermodynamic; /*!< \brief Constant pressure thermodynamic. */ const su2double GasConstant_Ref; /*!< \brief Gas constant reference needed for Nondimensional problems. */ const bool wilke; diff --git a/SU2_CFD/src/fluid/CFluidScalar.cpp b/SU2_CFD/src/fluid/CFluidScalar.cpp index 9c1ddc6006b..eb8f4624aec 100644 --- a/SU2_CFD/src/fluid/CFluidScalar.cpp +++ b/SU2_CFD/src/fluid/CFluidScalar.cpp @@ -48,13 +48,14 @@ CFluidScalar::CFluidScalar(su2double val_Cp, su2double val_gas_constant, su2doub n_species_mixture(config->GetnSpecies() + 1), Gas_Constant(val_gas_constant), Gamma(config->GetGamma()), + Pressure_Thermodynamic(value_pressure_operating), GasConstant_Ref(config->GetGas_Constant_Ref()), wilke(config->GetKind_MixingViscosityModel() == MIXINGVISCOSITYMODEL::WILKE), davidson(config->GetKind_MixingViscosityModel() == MIXINGVISCOSITYMODEL::DAVIDSON) { if (n_species_mixture > ARRAYSIZE) { SU2_MPI::Error("Too many species, increase ARRAYSIZE", CURRENT_FUNCTION); } - Pressure = value_pressure_operating; + for (int iVar = 0; iVar < n_species_mixture; iVar++) { molarMasses[iVar] = config->GetMolecular_Weight(iVar); specificHeat[iVar] = config->GetSpecific_Heat_CpND(iVar); @@ -215,7 +216,7 @@ void CFluidScalar::SetTDState_T(const su2double val_temperature, const su2double MassToMoleFractions(val_scalars); ComputeGasConstant(); Temperature = val_temperature; - Density = Pressure / (Temperature * Gas_Constant); + Density = Pressure_Thermodynamic / (Temperature * Gas_Constant); Cp = ComputeMeanSpecificHeatCp(val_scalars); Cv = Cp - Gas_Constant; From fab4c3566b13c2be9b2b4420d9a83822c6b7e7fd Mon Sep 17 00:00:00 2001 From: Cristopher-Morales Date: Thu, 8 Dec 2022 08:46:51 +0100 Subject: [PATCH 7/7] updating config template and test cases --- .../species2_primitiveVenturi_mixingmodel.cfg | 2 ++ .../species2_primitiveVenturi_mixingmodel_boundedscalar.cfg | 2 ++ .../species2_primitiveVenturi_mixingmodel_heatcapacity_H2.cfg | 2 ++ ...pecies2_primitiveVenturi_mixingmodel_heatcapacity_H2_ND.cfg | 2 ++ .../species2_primitiveVenturi_mixingmodel_viscosity.cfg | 2 ++ config_template.cfg | 3 +++ 6 files changed, 13 insertions(+) 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 5d4595e1a14..53f2764f995 100644 --- a/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel.cfg +++ b/TestCases/species_transport/venturi_primitive_3species/species2_primitiveVenturi_mixingmodel.cfg @@ -32,6 +32,8 @@ INC_NONDIM= DIMENSIONAL % FLUID_MODEL= FLUID_MIXTURE % +THERMODYNAMIC_PRESSURE= 101325.0 +% MOLECULAR_WEIGHT= 28.96, 16.043 % SPECIFIC_HEAT_CP = 1009.39, 1009.39 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 4c668a8ca0a..f51545a1174 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 @@ -32,6 +32,8 @@ INC_NONDIM= DIMENSIONAL % FLUID_MODEL= FLUID_MIXTURE % +THERMODYNAMIC_PRESSURE= 101325.0 +% MOLECULAR_WEIGHT= 28.96, 16.043 % SPECIFIC_HEAT_CP = 1009.39, 1009.39 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 0fce29f4142..24b943f677b 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 @@ -33,6 +33,8 @@ INC_NONDIM= DIMENSIONAL % FLUID_MODEL= FLUID_MIXTURE % +THERMODYNAMIC_PRESSURE= 101325.0 +% MOLECULAR_WEIGHT= 2.01588, 28.960 % SPECIFIC_HEAT_CP = 14310, 1009.39 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 6af9b3ed7c9..0f53f027374 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 @@ -34,6 +34,8 @@ INC_NONDIM= INITIAL_VALUES % FLUID_MODEL= FLUID_MIXTURE % +THERMODYNAMIC_PRESSURE= 101325.0 +% MOLECULAR_WEIGHT= 2.01588, 28.960 % SPECIFIC_HEAT_CP = 14310, 1009.39 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 adc7050dfb6..e92a8c9088c 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 @@ -32,6 +32,8 @@ INC_NONDIM= DIMENSIONAL % FLUID_MODEL= FLUID_MIXTURE % +THERMODYNAMIC_PRESSURE= 101325.0 +% MOLECULAR_WEIGHT= 16.043, 28.960 % SPECIFIC_HEAT_CP = 2224.43, 1009.39 diff --git a/config_template.cfg b/config_template.cfg index 5b76a40e403..9b676406417 100644 --- a/config_template.cfg +++ b/config_template.cfg @@ -337,6 +337,9 @@ CRITICAL_PRESSURE= 3588550.0 % Acentri factor (0.035 (air)) ACENTRIC_FACTOR= 0.035 % +% Thermodynamics(operating) Pressure (101325 Pa default value, only for incompressible flow and FLUID_MIXTURE) +THERMODYNAMIC_PRESSURE= 101325.0 +% % Specific heat at constant pressure, Cp (1004.703 J/kg*K (air)). % Incompressible fluids with energy eqn. (CONSTANT_DENSITY, INC_IDEAL_GAS) and the heat equation. SPECIFIC_HEAT_CP= 1004.703