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

Problem with turbulence model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 8, 2010, 16:47
Default Problem with turbulence model
  #1
New Member
 
Aditya
Join Date: Aug 2010
Posts: 6
Rep Power: 6
akonduri is on a distinguished road
Hi,

I am trying to do LES of flow past a square cylinder. For LES model, I am using locDynOneEqEddy. When I run the code, I get the following error. Can someone tell me want is going wrong?

============================
ERROR:
=============================

Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type LESModel
Selecting LES turbulence model locDynOneEqEddy
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 __kernel_standard at interp.c:0
#4 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#6 Foam::incompressible::LESModels::locDynOneEqEddy:: ck(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#7 Foam::incompressible::LESModels::locDynOneEqEddy:: updateSubGridScaleFields(Foam::GeometricField<Foam ::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#8 Foam::incompressible::LESModels::locDynOneEqEddy:: locDynOneEqEddy(Foam::GeometricField<Foam::Vector< double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#9 Foam::incompressible::LESModel::adddictionaryConst ructorToTable<Foam::incompressible::LESModels::loc DynOneEqEddy>::New(Foam::GeometricField<Foam::Vect or<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/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/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/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/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/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/akonduri/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleTurbulenceModel.so"
#13 main in "/home/akonduri/OpenFOAM/akonduri-1.6/applications/bin/linux64GccDPOpt/pisoFoamAdi"
#14 __libc_start_main in "/lib64/libc.so.6"
#15 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/akonduri/OpenFOAM/akonduri-1.6/applications/bin/linux64GccDPOpt/pisoFoamAdi"
akonduri is offline   Reply With Quote

Old   September 9, 2010, 14:56
Default
  #2
New Member
 
Aditya
Join Date: Aug 2010
Posts: 6
Rep Power: 6
akonduri is on a distinguished road
Can any one please tell what is the problem?
akonduri is offline   Reply With Quote

Old   September 17, 2010, 00:49
Default
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
I would guess a bad initialization or boundary condition setup. You have a division by zero somewhere in a power operator (zero argument of negative exponent power, maybe), given line #4.

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto 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
multiphaseInterFoam for RAS turbulence model chiven OpenFOAM Bugs 7 August 23, 2011 02:52
build your own turbulence model with buoyancy Thomas Baumann OpenFOAM 11 November 23, 2009 09:53
Fan heater model: what turbulence source to use? andy20 Main CFD Forum 0 March 2, 2008 13:46
problem with MFR model for multiphase mixing tanks Srinivas Main CFD Forum 1 November 7, 2005 15:16
multiphase mixing Problem with MRF model in MixSim Srinivas FLUENT 0 October 17, 2005 06:35


All times are GMT -4. The time now is 06:29.