|
[Sponsors] |
December 12, 2013, 06:18 |
Problem with dynOneEqEddy
|
#1 |
Senior Member
Join Date: Jan 2013
Posts: 134
Rep Power: 13 |
hi everyone!
I've came across a very strange error in my LES simulation with dynOneEqEddy model. When my velocity at the initial condition is set to 1m/s, it runs smoothly. However, when velocity is set to be 1.1m/s, or 0.9m/s, or even 0.999 and 1.001, the simulation breaks with the following error: 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 dynOneEqEddy Selecting LES delta type cubeRootVol bounding k, min: 0 max: 0.03 average: 0.03 #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so" #5 Foam::incompressible::LESModels::dynOneEqEddy::ck( Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so" #6 Foam::incompressible::LESModels::dynOneEqEddy::upd ateSubGridScaleFields(Foam::GeometricField<Foam::S ymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so" #7 Foam::incompressible::LESModels::dynOneEqEddy::dyn OneEqEddy(Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so" #8 Foam::incompressible::LESModel::adddictionaryConst ructorToTable<Foam::incompressible::LESModels::dyn OneEqEddy>::New(Foam::GeometricField<Foam::Vector< double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so" #9 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&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so" #10 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&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleLESModels.so" #11 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&, Foam::word const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleTurbulenceModel.so" #12 in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/pimpleFoam" #13 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6" #14 in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/pimpleFoam" Floating point exception (core dumped) I tried on several different computers, without any luck. From the message it seems that the error comes from the implementation of the dynOneEqEddy. Can someone explain this? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
UDF compiling problem | Wouter | Fluent UDF and Scheme Programming | 6 | June 6, 2012 04:43 |
Gambit - meshing over airfoil wrapping (?) problem | JFDC | FLUENT | 1 | July 11, 2011 05:59 |
natural convection problem for a CHT problem | Se-Hee | CFX | 2 | June 10, 2007 06:29 |
Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 19:13 |
Is this problem well posed? | Thomas P. Abraham | Main CFD Forum | 5 | September 8, 1999 14:52 |