[wangle]

I am a new user of OpenFOAM.I had write a new turbulenceModel,but there are a lot of things wrong when I used it to calculate a case.There are some error the compiler told me.But i don't know what's wrong with it.I need your help.Thank you very much!



Selecting RAS turbulence model kEpsilonNew1

1. 0 Foam::errorrintStack(Foam::Ostream&) in "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
2. 1 Foam::sigFpe::sigFpeHandler(int) in "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
3. 2 Uninterpreted: [0xffffe420]
4. 3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
5. 4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libfiniteVolume.so"
6. 5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foamperator/<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 "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libfiniteVolume.so"
7. 6 Foam::compressible::RASModels::kEpsilonNew1::kEpsi lonNew1(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&) in "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcompressibleRASModels.so"
8. 7 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<Foam::compressible::RASModels::kEpsilo nNew1>::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&) in "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcompressibleRASModels.so"
9. 8 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&) in "/home/lele/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcompressibleRASModels.so"
10. 9 main in "/home/lele/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/sonicTurbFoam"
11. 10 __libc_start_main in "/lib/libc.so.6"
12. 11 Foam::regIOobject::readIfModified() in "/home/lele/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/sonicTurbFoam"

floating number exception

[ngj]

Hi

You are having a 0 as boundary condition/initial condition on either k or epsilon. This is not allowed as you will divide by zero. Choose a physical sound and finite value.

/ Niels