Problem running IDDES with pimpleFoam
Dear OpenFOAM users,
I am trying to test the new SpalartAllmarasIDDES implementation in OpenFOAM 1.6. I am using pimpleFoam to test its capabilities to exceed a Courant number of 1. I believe I have set up everything correctly, however when I run my case the following error is output. I don't understand why it wants details of a RAS model, since I specified an LES model. I'm not sure if this is a bug or an error in my usage. I would be very grateful for any help! Best regards, Charlie. ---------------- Error output: Selecting incompressible transport model Newtonian Selecting turbulence model type LESModel Selecting LES turbulence model SpalartAllmarasIDDES request for RASModel RASProperties from objectRegistry region0 failed available objects of type RASModel are 0 ( ) #0 Foam::error::printStack(Foam::Ostream&) in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/pimpleFoam" #3 Foam::incompressible::RASModel const& Foam::objectRegistry::lookupObject<Foam::incompres sible::RASModel>(Foam::word const&) const in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleRASModels.so" #4 Foam::incompressible::RASModels::nutWallFunctionFv PatchScalarField::calcNut() const in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleRASModels.so" #5 Foam::incompressible::RASModels::nutWallFunctionFv PatchScalarField::updateCoeffs() in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleRASModels.so" #6 Foam::fvPatchField<double>::evaluate(Foam::Pstream ::commsTypes) in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/pimpleFoam" #7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/pimpleFoam" #8 Foam::incompressible::LESModels::SpalartAllmaras:: updateSubGridScaleFields() in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so" #9 Foam::incompressible::LESModels::SpalartAllmaras:: SpalartAllmaras(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/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so" #10 Foam::incompressible::LESModels::SpalartAllmarasID DES::SpalartAllmarasIDDES(Foam::GeometricField<Foa m::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so" #11 Foam::incompressible::LESModel::adddictionaryConst ructorToTable<Foam::incompressible::LESModels::Spa lartAllmarasIDDES>::New(Foam::GeometricField<Foam: :Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so" #12 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/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so" #13 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/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleLESModels.so" #14 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/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libincompressibleTurbulenceModel.so" #15 main in "/home/mockett/CFD/Solvers/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/pimpleFoam" #16 __libc_start_main in "/lib/tls/libc.so.6" #17 _start at ../sysdeps/i386/elf/start.S:122 From function objectRegistry::lookupObject<Type>(const word&) const in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.6/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 140. FOAM aborting Aborted |
Hi Charlie,
one chance is, that you used the wrong wall treatment in nuTilda!? It has to be something like nuSgs...WallTreatment. Best Regards! Fabian |
Hi Fabian,
Thanks very much - that seems to have solved the problem! Merry Christmas, Charlie. |
Problem running LRRDiffStress
Dear Fabian,
We are trying to test LRRDiffStress implementation in OpenFOAM 1.6 with LESModel. We are using buoyantPisoFoam for a cavity. When we run our case the following error is output. We don't understand the error message. We don't know if our boundary condition for B, nutilda or nusgs, are correct? We would be very grateful for any help! Ana and Sonia Reading g Reading thermophysical properties Reading field T Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting LES turbulence model LRRDiffStress LRRDiffStressCoeffs { ce 1.048; couplingFactor 0; ck 0.09; c1 1.8; c2 0.6; } Courant Number mean: 0 max: 0 Starting time loop Time = 0.001 Courant Number mean: 0 max: 0 DILUPBiCG: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 1.35638e-08, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 1.35807e-07, No Iterations 1 DICPCG: Solving for p, Initial residual = 1, Final residual = 0.363883, No Iterations 7 time step continuity errors : sum local = 3.56917e-05, global = -7.6233e-22, cumulative = -7.6233e-22 DICPCG: Solving for p, Initial residual = 0.155533, Final residual = 0.155533, No Iterations 0 time step continuity errors : sum local = 3.57045e-05, global = 1.23879e-21, cumulative = 4.76456e-22 DICPCG: Solving for p, Initial residual = 0.155533, Final residual = 8.23312e-07, No Iterations 71 time step continuity errors : sum local = 1.89001e-10, global = 1.71523e-21, cumulative = 2.19168e-21 #0 Foam::error::printStack(Foam::Ostream&) in "/home/foam6/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/foam6/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib64/libc.so.6" #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/foam6/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/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/foam6/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.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::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/foam6/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so" #6 Foam::incompressible::LESModels::LRRDiffStress::co rrect(Foam::tmp<Foam::GeometricField<Foam::Tensor< double>, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/foam6/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so" #7 Foam::incompressible::LESModel::correct() in "/home/foam6/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so" #8 main in "/home/foam6/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/flotacionLESFoam" #9 __libc_start_main in "/lib64/libc.so.6" #10 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/foam6/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/flotacionLESFoam" Floating point exception |
Hi,
I do not know this model, but I assume, that you have an initial field for B, which is set to 0... Fabian |
Dear Fabian
Sorry by our tardy answering, respect to B field it was have actually initialized on 0, we changed for to small value but the problem continued. We found that the problem was boundary condition, too. It changed it to small value and worked fine. Thank for you help Ana and Sonia |
Quote:
Hi Charlie, Can I ask how you sovle the problem. Since it doesn't work for me. I have tried may ways, still have the error, could you pleas kindly to tell me how to set the IDDES simulation. Thank you so much. Best regards, Joanna |
Quote:
Hi I have the same problem with SA DDES simulation. Could you tell me how you solved it? |
All times are GMT -4. The time now is 07:53. |