[vizvaz]

Hello,


I'm trying to simulate supersonic flow in nozzle using rhoCentralFoam for couple of weeks now. I have gone through the multiple posts on this forum and nothing helped in solving my issue. So, reaching to you guys at least few could help me in addressing my issue.



Problem Statement:- Supersonic flow in nozzle to understand shock waves



Inputs:-Nozzle is designed with stagnation pressure of 2.26e+06 pa, stagnation temperature of 1200k where the flow exits to pressure of 38521 pa.



Following are the Boundary conditions used



p

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

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

internalField   uniform 38521;

boundaryField
{
    inlet
    {
    //type        fixedValue;
    //value        uniform 2.267e+6;
    type           totalPressure;
    p0        uniform 2.267e+6;
    value        uniform 2.267e+6;
    }

    outlet
    {
    type         waveTransmissive;
    value         uniform 38521;        
    field        p;
    gamma         1.4;
    lInf         100;
    fieldInf     38521;
    }

    nozSym
    {
        type            symmetryPlane;
    }
    
    wall1
    {
        type            zeroGradient;
    }
    
    wall2
    {
        type            zeroGradient;
    }
    
    nozWall
    {
        type            zeroGradient;
    }

    frontAndBackPlanes
    {
        type            empty;
    }
}

// ************************************************************************* //

T

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

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 1200;

boundaryField
{
    inlet
    {
    type        totalTemperature;
    gamma        1.4;
    T0        uniform 1200;
    }

    outlet
    {
    type            zeroGradient;
    }

    nozSym
    {
        type            symmetryPlane;
    }
    
    wall1
    {
        type            zeroGradient;
    }
    
    wall2
    {
        type            zeroGradient;
    }
    
    nozWall
    {
        type            zeroGradient;
    }

    frontAndBackPlanes
    {
        type            empty;
    }
}

// ************************************************************************* //

u

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

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
    type           pressureInletOutletVelocity;
    phi        phi;
    value        uniform (0 0 0);
    }

    outlet
    {
    type            zeroGradient;
    }

    nozSym
    {
        type            symmetryPlane;
    }
    
    wall1
    {
        type            slip;
    }
    
    wall2
    {
        type            slip;
    }
    
    nozWall
    {
        type            slip;
    }

    frontAndBackPlanes
    {
        type            empty;
    }
}

// ************************************************************************* //

thermophysicalProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    //transport       sutherland;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        nMoles      1;
        molWeight   28.9;
    }
    thermodynamics
    {
        Cp          1005;
        Hf          0;
    }
    transport
    {
        //mu          1.84e-05;
        mu          0;
        Pr          0.713;

        As              1.4792e-06;
        Ts              116;
    }
}

//dpdt no;

// ************************************************************************* //

turbulenceProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType  laminar;

/*
RAS
{
    RASModel        kOmega;

    turbulence      on;

    printCoeffs     on;
}
*/

 // ************************************************************************* //

Issue:- Floating Point Exception. Solution runs for some time and gets terminated after some time step.



I have attached controlDict, fvSchemes and fvSolutions if anyone really wants to have a look at it. It doesn't have any issue with mesh and checkMesh returned OK. I understand the problem has ill-posed BC or setup.

openFoam experts seeking your advise to direct me to right approach.


Regards,
VV

[Fluid_]

Hey, I am also working on the same problems with similar boundary conditions. If you were able to solve the problem, can you share the details here?

[limengjie]

Hello, I am in trouble with 2D RANS supersonic jet fire. Similar boundary conditions and same error (float exception). Did you solve this problem? If you did, would you like help me? Thanks very much.