I compared myLTSPimpleFoam with simpleFoam. Ahmed body in previous articles was used for the test case. The computational grid is presented here.
The calculation setup for simpleFoam is almost identical to the previous case, but I made slight changes for fvSchemes and fvSolution.
The calculation setup for simpleFoam is almost identical to the previous case, but I made slight changes for fvSchemes and fvSolution.
fvSchemes (differences from the previous article are colored in red):
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(U) cellLimited Gauss linear 1;
}
divSchemes
{
default none;
div(phi,U) bounded Gauss GammaV 1;
div(phi,k) bounded Gauss Gamma 1;
div(phi,epsilon) bounded Gauss Gamma 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited corrected 0.333;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited corrected 0.333;
}
fluxRequired
{
default no;
p ;
}
fvSolution (differences from the previous article are colored in red):
solvers
{
"p.*"
{
solver GAMG;
tolerance 1e-8;
relTol 0;
smoother GaussSeidel;
nPreSweeps 1;
nPostSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 500;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U.*|k.*|epsilon.*)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0;
nSweeps 1;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 2;
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
epsilon 0.7;
}
}
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(U) cellLimited Gauss linear 1;
}
{
default none;
div(phi,U) bounded Gauss GammaV 1;
div(phi,k) bounded Gauss Gamma 1;
div(phi,epsilon) bounded Gauss Gamma 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited corrected 0.333;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited corrected 0.333;
}
fluxRequired
{
default no;
p ;
}
fvSolution (differences from the previous article are colored in red):
solvers
{
"p.*"
{
solver GAMG;
tolerance 1e-8;
relTol 0;
smoother GaussSeidel;
nPreSweeps 1;
nPostSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 500;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U.*|k.*|epsilon.*)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0;
nSweeps 1;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 2;
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
epsilon 0.7;
}
}
Comments
Post a Comment