CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problem running IDDES with pimpleFoam (http://www.cfd-online.com/Forums/openfoam-solving/71243-problem-running-iddes-pimplefoam.html)

charlie December 21, 2009 09:28

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

braennstroem December 22, 2009 12:14

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

charlie December 22, 2009 18:35

Hi Fabian,

Thanks very much - that seems to have solved the problem!

Merry Christmas,

Charlie.

gmc_unsa March 10, 2010 10:00

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

braennstroem March 11, 2010 15:57

Hi,

I do not know this model, but I assume, that you have an initial field for B, which is set to 0...

Fabian

gmc_unsa March 17, 2010 09:57

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

Joanna Huang April 30, 2013 05:05

Quote:

Originally Posted by charlie (Post 240771)
Hi Fabian,

Thanks very much - that seems to have solved the problem!

Merry Christmas,

Charlie.


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

菊爆大队总队长 July 11, 2013 11:45

Quote:

Originally Posted by Joanna Huang (Post 424110)
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


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 12:38.