CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   kOmegaSSTSAS model (http://www.cfd-online.com/Forums/openfoam-solving/67719-komegasstsas-model.html)

Kr_kim August 24, 2009 21:00

kOmegaSSTSAS model
 
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::printStack(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::operator/<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 August 25, 2009 04:01

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 September 1, 2009 22:59

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 October 19, 2009 14:20

kOmegaSSTSAS on channelFoam
 
Quote:

Originally Posted by ngj (Post 227386)
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 October 19, 2009 14:39

to clarify the problem
 
Quote:

Originally Posted by ngj (Post 227386)
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 February 4, 2010 13:46

moved my post to Bugs - thanks


All times are GMT -4. The time now is 08:20.