[arun1994]

Hello,


I have a film cooling problem of a flat plate with RANS (k-ep). The geometry is simple as in figure attached. I am using rhoSimpleFoam to solve the problem but I have a floating point exception always I have tried to play with the boundary conditions quite a bit but I still cannot find where the problem is.

[25] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[25] #1  Foam::sigFpe::sigHandler(int) at ??:?
[25] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[25] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
[25] #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
[25] #5  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
[25] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[25] #7  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
[ARUBUNTU:10065] *** Process received signal ***
[ARUBUNTU:10065] Signal: Floating point exception (8)
[ARUBUNTU:10065] Signal code:  (-6)
[ARUBUNTU:10065] Failing at address: 0x3e800002751
[ARUBUNTU:10065] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3efd0)[0x7fc201030fd0]
[ARUBUNTU:10065] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7fc201030f47]
[ARUBUNTU:10065] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3efd0)[0x7fc201030fd0]
[ARUBUNTU:10065] [ 3] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE9calculateEv+0x1e2)[0x7fc206c0eb02]
[ARUBUNTU:10065] [ 4] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11hePsiThermoINS_9psiThermoENS_11pureMixtureINS_19sutherlandTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE7correctEv+0x2a)[0x7fc206c4458a]
[ARUBUNTU:10065] [ 5] rhoSimpleFoam(+0x2889f)[0x562fd890489f]
[ARUBUNTU:10065] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7fc201013b97]
[ARUBUNTU:10065] [ 7] rhoSimpleFoam(+0x2d1aa)[0x562fd89091aa]
[ARUBUNTU:10065] *** End of error message ***

My 0 files look like this :
----------------------------------------------------------------------------------

alphat

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    Body.outlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_base
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.wall_side2
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.topwall
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.wall_side
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
    Body.inlet_coolant
    {
        type            calculated;
        value           uniform 0;
    }
    Body.inlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_coolant
    {
        type            compressible::alphatWallFunction;
        Prt             0.85;
        value           uniform 0;
    }
}

----------------------------------------------------------------------------------------------------
velocity U

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

internalField   uniform (0 0 0);

boundaryField
{
    Body.outlet
    {
        type            inletOutlet;
        value           uniform (0 0 0);
        inletValue      uniform (0 0 0);
    }
    Body.wall_base
    {
        type            noSlip;
    }
    Body.wall_side2
    {
        type            noSlip;
    }
    Body.topwall
    {
        type            noSlip;
    }
    Body.wall_side
    {
        type            noSlip;
    }
    Body.inlet_coolant
    {
    type surfaceNormalFixedValue;
    refValue uniform -10;
        rhoInlet        1.177;    // Guess for rho
    }
    Body.inlet
    {
    type surfaceNormalFixedValue;
    refValue uniform -10;
        rhoInlet        1;    // Guess for rho
    }
    Body.wall_coolant
    {
        type            noSlip;
    }
}

----------------------------------------------------------------------------------------------





Temperature

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    Body.outlet
    {

        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }
    Body.wall_base
    {
        type            zeroGradient;
    }
    Body.wall_side2
    {
        type            zeroGradient;
    }
    Body.topwall
    {
        type            zeroGradient;
    }
    Body.wall_side
    {
        type            zeroGradient;
    }
    Body.inlet_coolant
    {
        type            fixedValue;
        value           uniform 313;
    }
    Body.inlet
    {
        type            fixedValue;
        value           uniform 293;
    }
    Body.wall_coolant
    {
        type            zeroGradient;
    }
}

----------------------------------------------------------------------------------------------
Pressure

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

internalField   uniform 101325;

boundaryField
{
    Body.outlet
    {
        type            fixedValue;
        value           uniform 101325;

        // type            mixed;
        // refValue        uniform 110000;
        // refGradient     uniform 0;
        // valueFraction   uniform 1;
        // type            transonicOutletPressure;
        // U               U;
        // phi             phi;
        // gamma           1.4;
        // psi             psi;
        // pInf            uniform 110000;
    }
    Body.wall_base
    {
        type            zeroGradient;
    }
    Body.wall_side2
    {
        type            zeroGradient;
    }
    Body.topwall
    {
        type            zeroGradient;
    }
    Body.wall_side
    {
        type            zeroGradient;
    }
    Body.inlet_coolant
    {
        type            zeroGradient;
        // type            mixed;
        refValue        uniform 101325;
        refGradient     uniform 0;
        valueFraction   uniform 0.3;
    }
    Body.inlet
    {
        type            zeroGradient;
        // type            mixed;
        refValue        uniform 101325;
        refGradient     uniform 0;
        valueFraction   uniform 0.3;
    }
    Body.wall_coolant
    {
        type            zeroGradient;
    }
}

----------------------------------------------------------------------------------------------
nut

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

internalField   uniform 0;

boundaryField
{
    Body.outlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_base
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.wall_side2
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.topwall
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.wall_side
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
    Body.inlet_coolant
    {
        type            calculated;
        value           uniform 0;
    }
    Body.inlet
    {
        type            calculated;
        value           uniform 0;
    }
    Body.wall_coolant
    {
        type            nutkWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0;
    }
}

----------------------------------------------------------------------------------------------


k

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

internalField   uniform 1;

boundaryField
{
    Body.outlet
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 1;
    }
    Body.wall_base
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.wall_side2
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.topwall
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.wall_side
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
    Body.inlet_coolant
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;
        value           uniform 1;
    }
    Body.inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;
        value           uniform 1;
    }
    Body.wall_coolant
    {
        type            kqRWallFunction;
        value           uniform 1;
    }
}

----------------------------------------------------------------------------------------------
epsilon

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

internalField   uniform 200;

boundaryField
{
    Body.outlet
    {
        type            inletOutlet;
        inletValue      uniform 200;
        value           uniform 200;
    }
    Body.wall_base
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.wall_side2
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.topwall
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.wall_side
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
    Body.inlet_coolant
    {
        type            turbulentMixingLengthDissipationRateInlet;
        mixingLength    0.005;
        value           uniform 200;
    }
    Body.inlet
    {
        type            turbulentMixingLengthDissipationRateInlet;
        mixingLength    0.005;
        value           uniform 200;
    }
    Body.wall_coolant
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 200;
    }
}

----------------------------------------------------------------------------------------------
I am posting the BCs because I think this is where the problem should be.



If I am missing some details let me know please

Thanks for your time!
Arun