CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

1.6.x: Adaptive wall function: BoundaryFoam Tutorial

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 30, 2010, 05:43
Default 1.6.x: Adaptive wall function: BoundaryFoam Tutorial
  #1
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 530
Rep Power: 20
bastil is on a distinguished road
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
bastil is offline   Reply With Quote

Old   May 18, 2010, 19:20
Default
  #2
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 338
Rep Power: 18
louisgag is on a distinguished road
Send a message via ICQ to 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 is offline   Reply With Quote

Old   May 19, 2010, 19:20
Default
  #3
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 338
Rep Power: 18
louisgag is on a distinguished road
Send a message via ICQ to 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...

Last edited by louisgag; May 20, 2010 at 17:45.
louisgag is offline   Reply With Quote

Old   May 21, 2010, 04:49
Default
  #4
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 530
Rep Power: 20
bastil is on a distinguished road
Quote:
Originally Posted by louisgag View Post
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
bastil is offline   Reply With Quote

Old   May 25, 2010, 14:11
Default
  #5
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 338
Rep Power: 18
louisgag is on a distinguished road
Send a message via ICQ to 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
Code:
(1e-10 0 0)
Best regards,


-Louis
louisgag is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 10:31
Need some wall function approaches! yka8150 Main CFD Forum 0 September 22, 2009 00:08
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 15:00
Wall function tod Phoenics 1 May 19, 2003 06:05
Wall function in adverse pressure gradients stephane baralon Main CFD Forum 11 September 2, 1999 05:05


All times are GMT -4. The time now is 15:27.