Code:
thermophysicalProps:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType hPsiThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>>;mixture
{
specie
{
nMoles 1; //Stoffmenge in mol
molWeight 28.9; //Molmasse in g/mol
}
thermodynamics
{
Cp 1004; //Wärmekapazität cp
Hf 0; //Schmelzwärme Hf
}
transport
{
As 1.4584e-06; //Sutherlandkoeffizient
Ts 110.33; //Sutherlandtemperatur
}
}
// ************************************************************************* //
fvSolution:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0.05;
smoother GaussSeidel;
cacheAgglomeration off;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-7;
relTol 0.1;
}
"(k|omega)"
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-07;
relTol 0.1;
}
h
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0.1;
//solver smoothSolver;
//smoother DILUGaussSeidel;
//nSweeps 2;
//tolerance 1e-10;
//relTol 0.01;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 1;
rhoMax rhoMax [1 -3 0 0 0 0 0] 5;
rhoMin rhoMin [1 -3 0 0 0 0 0] 0;
residualControl
{
p 1e-9;
U 1e-6;
T 1e-7;
// possibly check turbulence fields
"(k|omega)" 1e-9;
}
}
relaxationFactors
{
p 0.5;
rho 0.01;
U 0.7;
"(k|omega)" 0.9;
h 0.7;
}
// ************************************************************************* //
fvSchemes:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(U) Gauss linear;
grad(p) Gauss linear;
}
divSchemes
{
default none;
// UEqn
div(phi,U) Gauss upwind;
div((muEff*dev2(T(grad(U))))) Gauss linear;
div(((muEff*(grad(U)+grad(U).T()))&U)) Gauss linear;
div((((0.666667*muEff)*div(U))*U)) Gauss linear;
div(U,p) Gauss upwind;
// hEqn
div(phi,h) Gauss upwind;
// omegaEqn
div(phi,omega) Gauss upwind;
// epsilonEqn
div(phi,epsilon) Gauss upwind;
// kEqn
div(phi,k) Gauss upwind;
// experimental
div(phid,p) Gauss upwind;
div(U) Gauss linear;
}
laplacianSchemes
{
default none;
// UEqn
laplacian(muEff,U) Gauss linear corrected;
// hEqn
laplacian(alphaEff,h) Gauss linear corrected;
// pCorrect
laplacian((rho*rAU),p) Gauss linear corrected;
laplacian((rho*(1|A(U))),p) Gauss linear corrected;
// kEqn
laplacian(DkEff,k) Gauss linear corrected;
// omegaEqn
laplacian(DomegaEff,omega) Gauss linear corrected;
// epsilonEqn
// laplacian(DepsilonEff,epsilon) Gauss linear corrected;
// potential
laplacian(1,p) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
div(U,p) upwind phi;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
}
// ************************************************************************* //