CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

correctors due to the mesh description

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 15, 2022, 10:20
Post correctors due to the mesh description
  #1
Member
 
Said CATALBAS
Join Date: Feb 2020
Location: Türkiye
Posts: 58
Rep Power: 4
saidc. is on a distinguished road
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;
}
Attached Images
File Type: png meshDecsription_FORUM.png (38.7 KB, 9 views)
saidc. is offline   Reply With Quote

Reply

Tags
accuracy, mesh, nonorthofaces

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[swak4Foam] funkyDoCalc with OF2.3 massflow NiFl OpenFOAM Community Contributions 14 November 25, 2020 04:30
hybrid mesh created by pointwise fails mesh check in fluent due to left-handed faces Hamid.de Pointwise & Gridgen 4 November 9, 2016 12:20
Fatal overflow error, singularity, due to local mesh form? Thomasbr CFX 7 July 8, 2016 10:57
error message cuteapathy CFX 14 March 20, 2012 07:45
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55


All times are GMT -4. The time now is 23:06.