[mihir1310]

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 .

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;
    }
}

U

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;
    }
}

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

[mihir1310]

this is the output :

diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error:rintStack(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:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::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