[saturn_K]

Hi

I'm trying to simulate gas dispersion in OpenFoam with rhoReacringBuoyantFoam.

Finally, I can solve all the obvious mistakes and errors; but still my program doesn't work and even I couldn't understand what's going on this time.

Is there anybody who can help me to solve this?

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |

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

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Reading thermophysical properties

Selecting thermodynamics package
{
type heRhoThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

Selecting chemistryReader foamChemistryReader
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
bounding k, min: 0 max: 0 average: 0
bounding omega, min: 0 max: 0 average: 0
RAS
{
RASModel kOmegaSST;
turbulence on;
printCoeffs on;
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.5555555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}

Creating reaction model

Selecting combustion model none

Reading g

Reading hRef
Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

Creating finite volume options from "constant/fvOptions"

Selecting finite volume options model type radiation
Source: radiation
Selecting radiationModel P1
Selecting absorptionEmissionModel greyMeanCombustion
Selecting scatterModel none
Selecting sootModel none
#0 Foam::error:rintStack(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::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoReactingBuoyantFoam"
#5 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Fo am::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foa m::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::F2() const at ??:?
#6 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Fo am::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foa m::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::F23() const at ??:?
#7 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Fo am::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foa m::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correctNut() at ??:?
#8 ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoReactingBuoyantFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoReactingBuoyantFoam"
Floating point exception (core dumped)