CFD Online Discussion Forums - rhoCentalFoam error

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - rhoCentalFoam error (https://www.cfd-online.com/Forums/openfoam-solving/72128-rhocentalfoam-error.html)

rhoCentalFoam error

Could anyone please help me with this :

Im trying to solve a free-jet problem by rhoCentralFoam [OF 1.5]. My domain represents a close chamber in which the jet enters at Supersonic conditions.
My code crashes due to a FPE , I am getting -ve pressure & density at the outlet .

this is obviously due to solving the h equation. I do not understand what is causing this error . ill post my BC;s & the thermoPhysicalProperties directory .

Code:

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



internalField   uniform 42000;



boundaryField

{

    wall

    {

        type            slip;

    }

    axis

    {

        type            empty;

    }

    outlet

    {

        type            nonReflective;

        refValue        uniform 6600;

        gamma           1.4;

        phi             phi;

        psi             psi;

        valueFraction   uniform 0.2;

        refGradient     uniform 0.2;

        value           uniform 6600;

    }

    inlet

    {

        type            totalPressure;

        U               U;

        phi             phi;

        rho             none;

        psi             none;

        gamma           1.4;

        p0              uniform 101325;

        value           uniform 101325;

    }

    frontAndBackPlanes

    {

        type            empty;

    }

    front

    {

        type            wedge;

    }

    back

    {

        type            wedge;

    }

}

Code:

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



internalField   uniform (260.39 0 0);



boundaryField

{

    wall

    {

        type            slip;

    }

    axis

    {

        type            empty;

    }

    outlet

    {

        type            zeroGradient;

    }

    inlet

    {

        type            fixedValue;

        value           uniform (315 0 0);

    }

    frontAndBackPlanes

    {

        type            empty;

    }

    front

    {

        type            wedge;

    }

    back

    {

        type            wedge;

    }

}

Code:

thermoType      hThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>>;



mixture         perfectGas 1 28.9 1004.3 200.0 1.458e-6 110.4;



rhoMin          rhoMin     [1 -3 0 0 0]     0.001;

Thank you

this is the output :

diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error::printStack(Foam::Ostream&) in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam::perfectGas> > > > >::calculate() in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam::perfectGas> > > > >::correct() in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5 main in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoCentralFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception