[bastil]

Hi,

following this thread:
http://www.cfd-online.com/Forums/ope...megasst-2.html

I am trying to run the BoundaryFoam tutorial with Adaptive wall functions:

• constant/polyMesh from launderSharma (low-Re Mesh)
• all the others from WallFunction

This combination gives a running tutorial. When modifing: nutWallFunction to nutSpalaratAllmarasWallFunction in o/nut it crashes:

Selecting RAS turbulence model kEpsilon
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/opt/OpenFOAM/OpenFOAM-1.6.x/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 "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#4  Foam::operator/(Foam::tmp<Foam::Field<double> > const&, Foam::UList<double> const&) in "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#5  Foam::incompressible::RASModels::nutSpalartAllmarasWallFunctionFvPatchScalarField::calcNut() const in "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#6  Foam::incompressible::RASModels::nutWallFunctionFvPatchScalarField::updateCoeffs() in "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#7  Foam::fvPatchField<double>::evaluate(Foam::Pstream::commsTypes) in "/opt/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/boundaryFoam"
#8  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/opt/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/boundaryFoam"
#9  Foam::incompressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#10  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kEpsilon>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#11  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/opt/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#12  main in "/opt/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/boundaryFoam"
#13  __libc_start_main in "/lib64/libc.so.6"
#14  __gxx_personality_v0 in "/opt/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/boundaryFoam"
Floating point exception

Can anybody reproduce this, please? What is going wrong.

The latest 1.6.x also have pure LowRe Nut Wall functions. Changing nutWallFunction to nutLowReWallFunction everything is running fine...

Ideas?

Regards Bastian

[louisgag]

Hi Bastian,

I have the same issue when running my own case with the nutSpallartAll..WallFunction. And setting to nutLowReWallFunction everything runs fine. Thus for now, I will simply use the LowRe function...

However, I don't know what is the issue.

Good luck with that,


-Louis

[louisgag]

Ok I solved it. It works now.

EDIT: false alarm, I had made a typo in the name and wrote it as nutSpallartAllmarasWallFunction, which made it run, but.. improperly...

[bastil]

Originally Posted by louisgag

Ok I solved it. It works now.

EDIT: false alarm, I had made a typo in the name and wrote it as nutSpallartAllmarasWallFunction, which made it run, but.. improperly...

Ok otherwise I would be interested in the solution.
I did not have the time to check but there were changes of the wall functions pushed into git after I stated this thread. Which version did you use?

Regards Bastian

[louisgag]

Hi Bastian.

- For my version number I couldn't find it; however, it dates from around early april.

- For the typo in "SpallartAllmaras" I found that if you give any weird name the same model is used, probably some default wall model. It is not nutWallFunction.

- For the working nutSpalartAllmarasWallFunction, try setting your wall U value to something like

(1e-10 0 0)

Best regards,


-Louis