CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

kOmegaSSTSAS model

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 24, 2009, 21:00
Wink kOmegaSSTSAS model
  #1
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
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
Kr_kim is offline   Reply With Quote

Old   August 25, 2009, 04:01
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,660
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
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
ngj is offline   Reply With Quote

Old   September 1, 2009, 22:59
Default
  #3
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
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 is offline   Reply With Quote

Old   October 19, 2009, 14:20
Default kOmegaSSTSAS on channelFoam
  #4
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
Quote:
Originally Posted by ngj View Post
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 is offline   Reply With Quote

Old   October 19, 2009, 14:39
Default to clarify the problem
  #5
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 8
Kr_kim is on a distinguished road
Quote:
Originally Posted by ngj View Post
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.
Kr_kim is offline   Reply With Quote

Old   February 4, 2010, 14:46
Default
  #6
Senior Member
 
Pierre-Olivier Dallaire
Join Date: Mar 2009
Location: Montreal, Quebec, Canada
Posts: 187
Rep Power: 8
podallaire is on a distinguished road
moved my post to Bugs - thanks

Last edited by podallaire; February 4, 2010 at 14:58. Reason: moved to Bugs
podallaire is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Reynolds Stress model in CFX vs Fluent Tim CFX 1 October 7, 2009 06:19
references about the fan/radiator model Mihai ARGHIR Main CFD Forum 1 January 8, 2001 16:49
references about the fan/radiator model Mihai ARGHIR Main CFD Forum 0 December 21, 2000 04:06
references about the fan/radiator model Mihai ARGHIR Main CFD Forum 1 December 17, 2000 08:01
references about the fan/radiator model Mihai ARGHIR FLUENT 0 December 17, 2000 07:40


All times are GMT -4. The time now is 09:57.