[Kr_kim]

I've successfully simulated kOmegaSST in OF-1.5
but I got following error message using kOmegaSSTSAS in OF-1.6
Is there anybody that know what this problem is and how to resolve it?

Thanks in advance


#0 Foam::error:rintStack(Foam::Ostream&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#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/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#5 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 "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#6 Foam::incompressible::LESModels::kOmegaSSTSAS::F2( ) const in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so"
#7 Foam::incompressible::LESModels::kOmegaSSTSAS::upd ateSubGridScaleFields(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so"
#8 Foam::incompressible::LESModels::kOmegaSSTSAS::kOm egaSSTSAS(Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so"
#9 Foam::incompressible::LESModel::adddictionaryConst ructorToTable<Foam::incompressible::LESModels::kOm egaSSTSAS>::New(Foam::GeometricField<Foam::Vector< double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so"
#10 Foam::incompressible::LESModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so"
#11 Foam::incompressible::turbulenceModel::addturbulen ceModelConstructorToTable<Foam::incompressible::LE SModel>::NewturbulenceModel(Foam::GeometricField<F oam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so"
#12 Foam::incompressible::turbulenceModel::New(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleTurbulenceModel.so"
#13 main in "/home/kkim64/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/pisoFoam"
#14 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#15 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122

[ngj]

Hi Kim

You divide by 0. This imply that you have put zero on either boundaries or internal volumes for k or omega. These needs to be finite, e.g. 1e-11.

Best regards,

Niels

[Kr_kim]

Thaks for the reply
Although I changed the values which you recommended but the same problem occurred. Is there any problem in code in itself?

http://openfoamwiki.net/index.php/HowTo_debugging

If there is some problem in code itself, how can I debug it(for specific files)

Best regards,

Kyungrock Kim

[Kr_kim]

Originally Posted by ngj

Hi Kim

You divide by 0. This imply that you have put zero on either boundaries or internal volumes for k or omega. These needs to be finite, e.g. 1e-11.

Best regards,

Niels

I have one more question on kOmegaSSTSAS. SAS mode was worked correctly in pisoFoam by quantifying correct initial and boudary data for omega and k at pisoFoam.

Question) But is it viable at channelFoam. It looks possible, I guess, because I have the same error with in pisoFoam(dividing by zero blah blah...) but I couldn't overcome it yet. To make sure it. is it possible in your opinion and if it is, what's the difference to pisoFoam?

[Kr_kim]

Originally Posted by ngj

Hi Kim

You divide by 0. This imply that you have put zero on either boundaries or internal volumes for k or omega. These needs to be finite, e.g. 1e-11.

Best regards,

Niels

Ubar is bulk velocity and it is used at pisoFoam and channelFoam with the same token. So k, epsilon and omega should be calculated based on Ubar values and asssign these values at internal Field, and the boundary condition at the both wall for k and omega should be zeroGradient. Are they right? I belive so. I'm waiting for your reply

Thanks in advance.

[podallaire]

moved my post to Bugs - thanks