zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Compressible solvers: Improved transonic option in the pressure equation

Browse Source 
The handling of the div(phid,p) term for transonic support in the pressure
equation is now consistent such that conservation is achieved at convergence of
the pressure system irrespective of the scheme chosen for div(phid,p) and the
relaxation of the pressure equation.

The rhoSimpleFoam tutorials have been updated and improved.
This commit is contained in:
Henry Weller
2022-02-22 13:05:40 +00:00
parent 69858a80ec
commit b827c0d740
17 changed files with 77 additions and 87 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
tutorials/compressible/rhoSimpleFoam/aerofoilNACA0012/constant/physicalProperties
View File

@ -22,7 +22,7 @@ thermoType
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
energy sensibleEnthalpy;
}
mixture // air at room temperature (293 K)
@ -31,11 +31,13 @@ mixture // air at room temperature (293 K)
{
molWeight 28.9;
}
thermodynamics
{
Cp 1005;
Hf 0;
}
transport
{
mu 1.82e-05;

4
tutorials/compressible/rhoSimpleFoam/aerofoilNACA0012/system/controlDict
View File

@ -15,13 +15,13 @@ FoamFile
application rhoSimpleFoam;
startFrom latestTime;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 2000;
endTime 5000;
deltaT 1;

11
tutorials/compressible/rhoSimpleFoam/aerofoilNACA0012/system/fvSchemes
View File

@ -34,19 +34,14 @@ divSchemes
div(phi,U) bounded Gauss linearUpwind limited;
turbulence bounded Gauss upwind;
energy bounded Gauss linearUpwind limited;
div(phi,h) $energy;
div(phi,K) $energy;
turbulence bounded Gauss upwind;
div(phi,k) $turbulence;
div(phi,omega) $turbulence;
div(phi,e) $energy;
div(phi,K) $energy;
div(phi,Ekp) $energy;
div(phid,p) Gauss upwind;
div((phi|interpolate(rho)),p) bounded Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

10
tutorials/compressible/rhoSimpleFoam/aerofoilNACA0012/system/fvSolution
View File

@ -18,12 +18,13 @@ solvers
p
{
solver GAMG;
smoother GaussSeidel;
smoother DIC;
tolerance 1e-8;
relTol 0.01;
}
"(U|k|omega|e)"
"(U|k|omega|h)"
{
solver PBiCGStab;
preconditioner DILU;
@ -38,7 +39,7 @@ SIMPLE
{
p 1e-6;
U 1e-5;
"(k|omega|e)" 1e-5;
"(k|omega|h)" 1e-5;
}
nNonOrthogonalCorrectors 0;
@ -51,10 +52,11 @@ relaxationFactors
p 0.3;
rho 0.01;
}
equations
{
U 0.7;
e 0.7;
h 0.7;
"(k|omega)" 0.7;
}
}

2
tutorials/compressible/rhoSimpleFoam/angledDuctExplicitFixedCoeff/constant/physicalProperties
View File

@ -22,7 +22,7 @@ thermoType
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
energy sensibleEnthalpy;
}
mixture

9
tutorials/compressible/rhoSimpleFoam/angledDuctExplicitFixedCoeff/system/fvSchemes
View File

@ -29,11 +29,14 @@ divSchemes
default none;
div(phi,U) bounded Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div(phi,e) bounded Gauss upwind;
div(phi,h) bounded Gauss upwind;
div(phi,K) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,Ekp) bounded Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes

17
tutorials/compressible/rhoSimpleFoam/angledDuctExplicitFixedCoeff/system/fvSolution
View File

@ -29,22 +29,24 @@ solvers
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-06;
tolerance 1e-6;
relTol 0.1;
}
e
h
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-06;
tolerance 1e-6;
relTol 0.1;
}
"(k|epsilon)"
{
$U;
tolerance 1e-07;
relTol 0.1;
}
@ -59,9 +61,7 @@ SIMPLE
p 1e-2;
U 1e-4;
e 1e-3;
// Possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
"(k|epsilon)" 1e-3;
}
}
@ -72,10 +72,11 @@ relaxationFactors
p 0.3;
rho 0.01;
}
equations
{
U 0.7;
e 0.7;
h 0.7;
"(k|epsilon)" 0.7;
}
}

5
tutorials/compressible/rhoSimpleFoam/squareBend/0/U
View File

@ -23,6 +23,7 @@ boundaryField
{
type noSlip;
}
inlet
{
type flowRateInletVelocity;
@ -30,11 +31,11 @@ boundaryField
rhoInlet 0.5; // Guess for rho
profile turbulentBL;
}
outlet
{
type inletOutlet;
type pressureInletOutletVelocity;
value uniform (0 0 0);
inletValue uniform (0 0 0);
}
}

2
tutorials/compressible/rhoSimpleFoam/squareBend/constant/physicalProperties
View File

@ -22,7 +22,7 @@ thermoType
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
energy sensibleEnthalpy;
}
mixture

10
tutorials/compressible/rhoSimpleFoam/squareBend/system/fvSchemes
View File

@ -29,14 +29,16 @@ divSchemes
default none;
div(phi,U) bounded Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div(phi,e) bounded Gauss upwind;
div(phi,h) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,K) bounded Gauss upwind;
div(phid,p) Gauss upwind;
div(phi,Ekp) bounded Gauss upwind;
div((phi|interpolate(rho)),p) Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes

21
tutorials/compressible/rhoSimpleFoam/squareBend/system/fvSolution
View File

@ -19,19 +19,19 @@ solvers
p
{
solver GAMG;
smoother DILU;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;
}
"(U|e|k|epsilon)"
"(U|h|k|epsilon)"
{
solver GAMG;
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;
}
}
@ -46,10 +46,8 @@ SIMPLE
{
p 1e-3;
U 1e-4;
e 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
h 1e-3;
"(k|epsilon)" 1e-3;
}
}
@ -59,11 +57,12 @@ relaxationFactors
{
p 1;
}
equations
{
p 1;
U 0.9;
e 0.8;
h 0.9;
k 0.9;
epsilon 0.9;
}

2
tutorials/compressible/rhoSimpleFoam/squareBendLiq/constant/physicalProperties
View File

@ -19,7 +19,7 @@ thermoType
type heRhoThermo;
mixture pureMixture;
properties liquid;
energy sensibleInternalEnergy;
energy sensibleEnthalpy;
}
mixture

6
tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSchemes
View File

@ -31,12 +31,12 @@ divSchemes
default none;
div(phi,U) bounded Gauss linearUpwindV limited;
div(phi,h) bounded Gauss linearUpwind limited;
div(phi,e) bounded Gauss linearUpwind limited;
div(phi,K) bounded Gauss linearUpwind limited;
div(phi,epsilon) bounded Gauss linearUpwind limited;
div(phi,k) bounded Gauss linearUpwind limited;
div(phi,Ekp) bounded Gauss linearUpwind limited;
div(phi,K) bounded Gauss linearUpwind limited;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

6
tutorials/compressible/rhoSimpleFoam/squareBendLiq/system/fvSolution
View File

@ -25,7 +25,7 @@ solvers
relTol 0.1;
}
"(U|h|e|k|epsilon)"
"(U|h|k|epsilon)"
{
solver smoothSolver;
smoother symGaussSeidel;
@ -46,7 +46,7 @@ SIMPLE
{
p 1e-4;
U 1e-4;
"(k|omega|epsilon|h|e)" 1e-3;
"(k|omega|epsilon|h)" 1e-3;
}
}
@ -56,11 +56,11 @@ relaxationFactors
{
p 0.3;
}
equations
{
U 0.7;
h 0.7;
e 0.7;
k 0.7;
epsilon 0.7;
}