[akonduri]

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]

Can any one please tell what is the problem?

[alberto]

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,