[charlie]

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:rintStack(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:perator<< <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:bjectRegistry::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]

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]

Hi Fabian,

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

Merry Christmas,

Charlie.

[gmc_unsa]

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:rintStack(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:perator/<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]

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]

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]

Originally Posted by charlie

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

[菊爆大队总队长]

Originally Posted by Joanna Huang

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?