CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

[isoAdvector] Different results with polyhedral and hexahedral mesh

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 21, 2020, 07:31
Default [isoAdvector] Different results with polyhedral and hexahedral mesh
  #1
New Member
 
Liz
Join Date: Apr 2020
Location: Germany
Posts: 10
Rep Power: 6
redTo is on a distinguished road
Hello eveybody,

I'm using compressibleInterIsoFoam (OF v2006) for a nozzle spray simulation and I'm comparing the results between polyMesh and hexaMesh.

For both cases I have defined the same conditions and the results are really different (see attached file).

Any idea why is this happening?

I'm grateful for any help!

Cheers,
Liz

PS: Inlet condition: nozzle bottom 4,6 g/s
outlet condition: outlet box 1 bar

fvSolution
Code:
solvers
{
    "alpha.oil.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 3;
        cAlpha          1;

        reconstructionScheme plicRDF;
        vof2IsoTol 	1e-8;
        surfCellTol 	1e-6; 
        nAlphaBounds    3;
        snapTol         1e-12;
        clip            true;

    }

    "pcorr.*"
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-7;
        relTol          0.01;
    }
    
    ".*(rho|rhoFinal)"
    {
        solver          diagonal;
    }

    p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-8;
        relTol          0.01;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "(U|T|k|epsilon|B).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-08;
        relTol          0;
    }
}

PIMPLE
{
    momentumPredictor no;    
    nOuterCorrectors 1;
    nCorrectors     3;
    nNonOrthogonalCorrectors 0;

}

relaxationFactors
{
     fields
    {
        "p.*"           0.1; 
    }
    equations
    {
        "k.*" 0.7;
	"epsilon.*" 0.7;
	"U.*" 0.8;
    }
}
fvSchemes
Code:
 ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    div(rhoPhi,U)  Gauss upwind;
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
 
    div(phi,thermo:rho.oil) Gauss upwind;
    div(phi,thermo:rho.air) Gauss upwind;

    div(phi,k)      Gauss upwind;
    div(rhoPhi,k)   Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(rhoPhi,epsilon)  Gauss upwind;

    div(rhoPhi,T)  Gauss upwind;
    div(rhoPhi,K)  Gauss upwind;
    div(phi,p)      Gauss upwind;

    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}
Attached Files
File Type: pdf nozzle.pdf (183.2 KB, 28 views)

Last edited by redTo; October 21, 2020 at 08:41.
redTo is offline   Reply With Quote

Reply

Tags
compressibleinterfoam, isoadvector

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
foam-extend-4.1 release hjasak OpenFOAM Announcements from Other Sources 19 July 16, 2021 05:02
Polyhedral Mesh Diego Flores FLUENT 17 August 21, 2020 13:52
Fluent warning Hexahedral or polyhedral cell problems lehoanganh07 FLUENT 1 November 6, 2013 01:49
[ANSYS Meshing] Hexahedral and polyhedral cell problems lehoanganh07 ANSYS Meshing & Geometry 0 November 3, 2013 09:42
Fluent 6.3.16 + Polyehdral Meshes Carlos FLUENT 3 March 15, 2007 03:09


All times are GMT -4. The time now is 15:47.