saidc. |
September 15, 2022 09:20 |
correctors due to the mesh description
1 Attachment(s)
Hi,
I've a really complex geometry and I can't create better mesh than I've already done. Because of this I need to force the foam to get accurate results but it diverges with high nNonOrthogonalCorrectors. I calculated 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999 quantiles of mesh data to describe it as much as possible. I believe 0.9999 is the one that carrying more information. The simulation is laminar, closed natural convection case. All boundries are wall and has adiabatic bc. Solver is buoyantPimpleFoam.
I want to know your ideas how you would set the fvSolution and fvSchemes for this mesh. I attached meshDescription.png, please check it.
Code:
solvers
{
"rho.*"
{
solver PCG;
preconditioner DIC;
tolerance 0;
relTol 0;
}
p_rgh
{
solver GAMG;
tolerance 1e-06;
relTol 0.01;
smoother GaussSeidel;
}
p_rghFinal
{
solver GAMG;
tolerance 1e-07;
relTol 0.0;
smoother GaussSeidel;
}
"(U|h)"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-08;
relTol 0.0;
}
"(U|h)Final"
{
$U;
tolerance 1e-08;
relTol 0.0;
}
}
PIMPLE
{
momentumPredictor no;
nNonOrthogonalCorrectors 2;
nOuterCorrectors 1;
nCorrectors 3;
pRefCell 0;
pRefValue 1e5;
}
Code:
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(U) cellLimited Gauss linear 0.75;
}
divSchemes
{
default none;
div(phi,U) Gauss linearUpwind grad(U);
div(phi,h) Gauss upwind;
div(phi,K) Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited 0.75;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited 0.75;
}
|