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/)
-   -   1.6.x: Adaptive wall function: BoundaryFoam Tutorial (http://www.cfd-online.com/Forums/openfoam-solving/74317-1-6-x-adaptive-wall-function-boundaryfoam-tutorial.html)

bastil March 30, 2010 04:43

1.6.x: Adaptive wall function: BoundaryFoam Tutorial
 
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:
Code:

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 May 18, 2010 18:20

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 May 19, 2010 18:20

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 May 21, 2010 03:49

Quote:

Originally Posted by louisgag (Post 259518)
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 May 25, 2010 13:11

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
Code:

(1e-10 0 0)
Best regards,


-Louis


All times are GMT -4. The time now is 13:49.