[jmart7]

Hello everyone,

So I am trying to simulate a multiple spray injector case based on the aachenBomb for sprayFoam but using a different geometry, but when the case is being initialized I get a floating point exception after the kEpsilon RAS model is selected. Looking at the error, I can see that it has do do with calculating the kEpsilon coefficients when using mutkWallFunction in the boundaries.

I have tried dieselFoam and I get the same error. Here is the output, any help would be greatly appreciated! Thanks!

Create mesh for time = 0

Reading chemistry properties


Reading g
Reading thermophysical properties

Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hsPsiMixtureThermo<reactingMixture<gasThermoPhysic s>>
Selecting chemistryReader chemkinReader
ODEChemistryModel: Number of species = 5 and reactions = 1
Selecting ODE solver SIBS
Creating component thermo properties:
multi-component carrier - 5 species
liquids - 1 components
solids - 0 components

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
#0 Foam::error:rintStack(Foam::Ostream&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam::compressible::RASModels::mutkWallFunctionFvP atchScalarField::calcMut() const in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#4 Foam::compressible::RASModels::mutkWallFunctionFvP atchScalarField::updateCoeffs() in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#5 Foam::fvPatchField<double>::evaluate(Foam::UPstrea m::commsTypes) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"
#7 Foam::compressible::RASModels::kEpsilon::kEpsilon( Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#8 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<Foam::compressible::RASModels::kEpsilo n>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#9 Foam::compressible::RASModel::New(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#10 Foam::compressible::turbulenceModel::addturbulence ModelConstructorToTable<Foam::compressible::RASMod el>::NewturbulenceModel(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#11 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#12 main in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"
#13 __libc_start_main in "/lib64/tls/libc.so.6"
#14 Foam::UOPstream::write(char) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"

[Smed]

Any luck resolving this? I'm getting a similar error.

[kornickel]

---POST DELETED---

Sorry, wrong thread

[dkxls]

sprayFoam has basically been unusable ever since it was introduced in 2.0.x.

Even if in some cases it would have produced some results (i.e. didn't crash), these results were most likely wrong - at least for engine spray simulation!

I fixed most of the issues concerning the thermodynamics and contribute them to upstream OpenFOAM (2.2.x), but there are still several issues, especially concerning the breakup models.
For more information, see the list of bugs I reported/fixed in the last couple of weeks (and there are still more to come...):
http://www.openfoam.org/mantisbt/vie...eporter_id=294

Also have a look here:
http://www.cfd-online.com/Forums/ope...parameter.html

To get started with sprayFoam, get a recent git checkout, set the atomization model to "none" (it's broken altogether right now) and use maybe the ReitzDiwakar or PilchErdman breakup model (at least don't use ReitzKHRT or TAB/ETAB).

I would very much appreciate more people testing this stuff...

[kornickel]

---POST DELETED---

Sorry, wrong thread

[dkxls]

Oho!
I just read the header and didn't exactly read your message, so yes you should probably start a new thread.

And then again, my post is not related to RANS turbulence models at all, so ... also questionable if it really belongs here.