[rishibigghe]

Hello,

I have been trying to simulate a jet flow of methane in air. I am using reactingFoam solver with chemistry turned off as I am only interested in mass fraction without reacting. I am running the solver to achieve a steady state i.e. increasing time continuously to achieve a steady state result with a transient solver.

I am using a coarse mesh just in case to determine a rough approximate of steady state time and arrive to a solution to sample set-up. I am running it on openfoam 1.5 in one machine and in openfoam 1.7 in another machine.

When i run in openfoam 1.5, the case continues without any error. But when i run in oepnfoam 1.7, the cases crashes unexpected with floating point exception. Following is the error message that i get.

__________________________________________________ _____________________
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::LimitedScheme<double, Foam::limitedLinearLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::fv::gaussConvectionScheme<double>::fvmDiv(Fo am::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) const in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam"
#7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam"
#8 Foam::compressible::RASModels::kEpsilon::correct() in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#9
in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam"
#10 __libc_start_main in "/lib64/libc.so.6"
#11
in "/home/ddovizio/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam"
Floating point exception (core dumped)
__________________________________________________ ________________

Can anyone explain as to why I am experiencing this error?

I have another query. Before performing the above simulation i performed a jet of air in air just to check the geometry,mesh, velocity profiles. The velocity profile that i got with simpleFoam (for air jet) should be same as reactingFoam velocity profile (for methane jet). But, in actual it is coming completely different instead of all the boundary conditions remaining same.

In my openfoam 1.5 the UEqn.H file was missing from the reactingFoam solver. Is it there in yours.

Kindly reply if you can help in resolving the problem

Regards
Rishi