From 5ef8edfc0cdb3647f610cd59013a4ccd4957f544 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Thu, 2 Jul 2020 09:42:07 +0100 Subject: [PATCH] tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner: Changed to perfectGas Only perfectGas and real-gas equations of state are consistent with standard Janaf thermo data based on Cp. Using other equations of state is possible but the Janaf Cp data would have to be modified for consistency. --- .../chtMultiRegionFoam/reverseBurner/0/gas/p | 2 +- .../reverseBurner/0/gas/p_rgh | 2 +- .../reverseBurner/0/solid/p | 2 +- .../reverseBurner/constant/gas/fvOptions | 2 +- .../constant/gas/thermo.compressibleGas | 20 +++++-------------- .../constant/gas/thermophysicalProperties | 4 ++-- .../reverseBurner/system/gas/fvSchemes | 2 +- .../reverseBurner/system/gas/fvSolution | 4 ++-- 8 files changed, 14 insertions(+), 24 deletions(-) diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p index 433d372bf8..2e1d4666db 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p @@ -17,7 +17,7 @@ FoamFile dimensions [ 1 -1 -2 0 0 0 0 ]; -internalField uniform 0; +internalField uniform 1e6; boundaryField { diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh index 70c835912b..bbb80ddd00 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh @@ -17,7 +17,7 @@ FoamFile dimensions [ 1 -1 -2 0 0 0 0 ]; -internalField uniform 0; +internalField uniform 1e6; boundaryField { diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p index addefc7a2e..1d33d3f9fa 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p @@ -17,7 +17,7 @@ FoamFile dimensions [ 1 -1 -2 0 0 0 0 ]; -internalField uniform 0; +internalField uniform 1e6; boundaryField { diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions index 45592dc58e..6c73634312 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions @@ -29,7 +29,7 @@ options sources { - e + h { explicit 5e7; // kg/m/s^3 implicit 0; diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas index e07b9a98f5..6f04539719 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas @@ -22,9 +22,7 @@ O2 molWeight 31.9988; } equationOfState - { - pRef 1e6; - } + {} thermodynamics { Tlow 200; @@ -47,9 +45,7 @@ H2O molWeight 18.0153; } equationOfState - { - pRef 1e6; - } + {} thermodynamics { Tlow 200; @@ -72,9 +68,7 @@ CH4 molWeight 16.0428; } equationOfState - { - pRef 1e6; - } + {} thermodynamics { Tlow 200; @@ -97,9 +91,7 @@ CO2 molWeight 44.01; } equationOfState - { - pRef 1e6; - } + {} thermodynamics { Tlow 200; @@ -122,9 +114,7 @@ N2 molWeight 28.0134; } equationOfState - { - pRef 1e6; - } + {} thermodynamics { Tlow 200; diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties index a3b6403128..8779e3be95 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties @@ -21,9 +21,9 @@ thermoType mixture multiComponentMixture; transport sutherland; thermo janaf; - equationOfState incompressiblePerfectGas; + equationOfState perfectGas; specie specie; - energy sensibleInternalEnergy; + energy sensibleEnthalpy; } inertSpecie N2; diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes index 636767cb23..4349ee1971 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes @@ -30,7 +30,7 @@ divSchemes default none; div(phi,U) Gauss linearUpwind grad(U); - div(phi,e) Gauss linearUpwind limitedGrad; + div(phi,h) Gauss linearUpwind limitedGrad; div(phi,Yi_h) Gauss limitedLinear 1; div(phi,k) Gauss linearUpwind limitedGrad; div(phi,epsilon) Gauss linearUpwind limitedGrad; diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution index 7cc406ecc6..c5028bb07a 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution @@ -36,7 +36,7 @@ solvers relTol 0; } - "(U|e|k|epsilon|Yi)" + "(U|h|k|epsilon|Yi)" { solver PBiCGStab; preconditioner DILU; @@ -44,7 +44,7 @@ solvers relTol 0.1; } - "(U|e|k|epsilon|Yi)Final" + "(U|h|k|epsilon|Yi)Final" { $U; tolerance 1e-7;