[arun1994]

Hey all, I am trying to solve a simple flat plate film cooling problem. I am attaching the 0, constant and system files.

This is the error I get.

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: Convergence criteria found
        p: tolerance 0.001
        U: tolerance 0.0001
        e: tolerance 0.001
        "(k|epsilon|omega)": tolerance 0.001

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

pressureControl
    pMax 0
    pMin 0

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

No MRF models present

No finite volume options present
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::operator/(Foam::tmp<Foam::Field<double> > const&, Foam::UList<double> const&) at ??:?
#5  Foam::CompressibleTurbulenceModel<Foam::fluidThermo>::nu(int) const at ??:?
#6  Foam::nutkWallFunctionFvPatchScalarField::nut() const at ??:?
#7  Foam::nutWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#8  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
#9  Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correctNut() at ??:?
#10  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
 Floating point exception (core dumped)

In this problem, there is a mainstream flow and a " jet in crossflow". I started to use rhosimplefoam as it is easier to add temperatures to the inlet flows. I have read elsewhere that the above floating point exception error occurs when there is a division by wero for k, but I really cant find the error ( maybe it is in front of my eyes, but I spent enough time finding it before posting here)



Will do wonders if someone could help

[dlahaye]

Can you try to solve isothermal case using simpleFoam first?

[arun1994]

Hello, I have already tried the iso thermal case and it runs well using simpleFoam