CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

LES case diverges after 5 flow-through-time

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 5, 2016, 04:19
Default LES case diverges after 5 flow-through-time
  #1
New Member
 
James Zh
Join Date: Mar 2016
Posts: 3
Rep Power: 3
mcshidada is on a distinguished road
Hi Foamers,

This is my first post

I'm doing Ma7 LES on rhoCentralFoam with kEqn model. The geometry is quite simple. The mapped inflow boundary condition is used as in the attached geometry. The back and forth patches are cyclic, while the top one is supersonicfreestream.

details of physical and numerical boundary conditions:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

6
(
    INLET
    {
        type            mappedPatch;
        inGroups        1(mappedPatch);
        nFaces          5355;
        startFace       5836184;
        sampleMode      nearestCell;
        sampleRegion    region0;
        samplePatch     none;
        offsetMode      uniform;
        offset          (0.06 0 0);
    }
    SYM
    {
        type            cyclic;
        inGroups        1(cyclic);
        nFaces          31195;
        startFace       5841539;
        matchTolerance  0.0006;
        transform       unknown;
        neighbourPatch  SYM2;
    }
    SYM2
    {
        type            cyclic;
        inGroups        1(cyclic);
        nFaces          31195;
        startFace       5872734;
        matchTolerance  0.0006;
        transform       unknown;
        neighbourPatch  SYM;
    }
    WALL_D
    {
        type            wall;
        inGroups        1(wall);
        nFaces          23121;
        startFace       5903929;
    }
    OUTLET
    {
        type            patch;
        nFaces          5355;
        startFace       5927050;
    }
    FAROUT
    {
        type            patch;
        nFaces          23121;
        startFace       5932405;
    }
)

// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    INLET
    {
        type            calculated;
        value           uniform 0;
    }
    FAROUT
    {
        type            calculated;
        value           uniform 0;
    }
    OUTLET
    {
        type            calculated;
        value           uniform 0;
    }

    SYM
    {
        type            cyclic;
    }
    SYM2
    {
        type            cyclic;
    }

    WALL_D
    {
        type            zeroGradient;
    }

}


// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 4925.8;

boundaryField
{
    INLET
    {
        type                mapped;
        value               uniform 4925.8;
        interpolationScheme cell;
        setAverage          true;
        average             4925.8;
    }
    FAROUT
    {
        type            zeroGradient;
    }
    OUTLET
    {
        type            zeroGradient;
    }

    SYM
    {
        type            cyclic;
    }
    SYM2
    {
        type            cyclic;
    }


    WALL_D
    {
        type            zeroGradient;
    }
}


// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    INLET
    {
        type            zeroGradient;
    }
    FAROUT
    {
        type            zeroGradient;
    }
    OUTLET
    {
        type            zeroGradient;
    }

    SYM
    {
        type            cyclic;
    }
    SYM2
    {
        type            cyclic;
    }


    WALL_D
    {
        type            zeroGradient;
    }
}


// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 1341.28589;

boundaryField
{
    INLET
    {
        type            fixedValue;
        value           uniform 1341.28589;
    }
    FAROUT
    {
        type            zeroGradient;
    }

    OUTLET
    {
        type            waveTransmissive;
        field           p;
        phi             phi;
        rho             rho;
        psi             thermo:psi;
        gamma           1.4;
        fieldInf        1341.28589;
        lInf            10;
        value           uniform 1341.28589;
    }

    SYM
    {
        type            cyclic;
    }
    SYM2
    {
        type            cyclic;
    }

    WALL_D
    {
        type            zeroGradient;
    }
}

// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 62.9;

boundaryField
{
    INLET
    {
        type                fixedValue;
        value               uniform 62.9;
    }
    FAROUT
    {
        type            zeroGradient;
    }

    OUTLET
    {
        type            zeroGradient;
    }

    SYM
    {
        type            cyclic;
    }
    SYM2
    {
        type            cyclic;
    }


    WALL_D
    {
        type            zeroGradient;
    }
}

// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (1146.1 0 0);

boundaryField
{
    INLET
    {
        type                mapped;
        value               uniform (1146.1 0 0);
        interpolationScheme cell;
        setAverage          false;
        average             (1146.1 0 0);
    }
    FAROUT
    {
        type            supersonicFreestream;
        pInf            1341.28589;
        TInf            62.9;
        UInf            (1146.1 0 0);
        gamma           1.4;
        value           uniform (1146.1 0 0);
    }

    OUTLET
    {
        type            pressureInletOutletVelocity;
        inletValue      uniform (0 0 0);
        value           uniform (0 0 0);
    }

    SYM
    {
        type            cyclic;
    }
    SYM2
    {
        type            cyclic;
    }


    WALL_D
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}

// ************************************************************************* //
Schemes and solution file
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fluxScheme      Kurganov;

ddtSchemes
{
    default         backward;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(tauMC)      Gauss linear;
    div(phi,k)      Gauss Gamma 0.2;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         none;
    reconstruct(rho)    vanLeer;
    reconstruct(U)      vanLeerV;
    reconstruct(T)      vanLeer;
    interpolate(rhoU)   linear;
    interpolate(muEff)  linear;
    interpolate(tauMC)  linear;
    interpolate(rho)    linear;
}

snGradSchemes
{
    default         corrected;
}


// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "(rho|rhoU|rhoE)"
    {
        solver          diagonal;
    }

    U
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-07;
        relTol          0;
    }

    e
    {
        $U;
        tolerance       1e-7;
        relTol          0;
    }
    k
    {
        $U;
        tolerance       1e-7;
        relTol          0;
    }
}


// ************************************************************************* //
There are field contours of k, pressure, velocity magnitude and temperature in the attachment at time step just before it diverged.
I monitored the mass flow rate of all the patches and it seemed OK with the inlet, outlet and the freestream one. The mass flow rate of cyclic patch kept increasing up to around 0.01 of the oulet's, is it common? Maybe because the turbulent boundary layer is still developing, i am not sure.
The case has been running for about two days(around 5 flow-through-time), finally blow up...
Delta t is set to 1e-8 and the maxCo is less than 0.2 when it diverged.
I have stuck here for a week...
Any help will be appreciated.

Regards
James Zh
Attached Images
File Type: png g.png (24.4 KB, 7 views)
Attached Files
File Type: zip k_p_contour.zip (125.1 KB, 1 views)
File Type: zip T_U_contour.zip (117.9 KB, 1 views)

Last edited by mcshidada; October 5, 2016 at 07:01.
mcshidada is offline   Reply With Quote

Old   October 5, 2016, 06:36
Default
  #2
New Member
 
James Zh
Join Date: Mar 2016
Posts: 3
Rep Power: 3
mcshidada is on a distinguished road
Any comment is welcomed
mcshidada is offline   Reply With Quote

Reply

Tags
high-speed, les, rhocentralfoam

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM Running, Solving & CFD 8 March 28, 2016 11:08
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
Same SimpleFOAM Case converges with openFOAM 2.1 but diverges with openFOAM 2.0.1 alsdia OpenFOAM Running, Solving & CFD 3 October 22, 2012 11:25
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


All times are GMT -4. The time now is 05:11.