CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

OpenFOAM crashes using high resolution grid at wall

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

LinkBack Thread Tools Display Modes
Old   February 17, 2014, 10:44
Default OpenFOAM crashes using high resolution grid at wall
Senior Member
Join Date: Jun 2012
Location: Germany, Bochum
Posts: 146
Rep Power: 8
Bazinga is on a distinguished road
Dear all,

I am really new to OpenFOAM, and after going through some of the tutorials I found the simpleFoam/turbineSiting tutorial which helped me to create a simulation with an atmospheric boundary layer (ABL).

For my own test simulation I wanna start with an empty domain with an ABL just going through it. This works fine for an uniform mesh of 0.1m in a 1.6x1.8x12m domain.

For my further studies it is necessary to have a high mesh resolution of the boundary layer at the bottom of the domain. So, I tried to put the node of the first cell at the bottom at about y+ = 30 which should be above the buffer layer to use wall functions for the viscous layer at the wall.

In my case I considered the freestream velocity as 6m/s, and the boundary length as 1.2m which gave me - using the cfd-online tool - a y+=30 wall distance of 0.0016m. However, when I make the mesh this dense at the bottom of my domain the simulation crashes in OpenFoam. Here is the errorcode if that helps:

time step continuity errors : sum local = 3.59678e+85, global = 7.97804e+75, cumulative = 7.97804e+75
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/"
#2 Uninterpreted:
#3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/"
#4 void Foam::multiply<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 "/opt/openfoam222/platforms/linuxGccDPOpt/lib/"
#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 "/opt/openfoam222/platforms/linuxGccDPOpt/lib/"
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/"
in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/simpleFoam"
#8 __libc_start_main in "/lib/i386-linux-gnu/"
in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)
The simulation settings are mostly copy pasted from the turbineSiting tutorial although I set the wall function for the bottom of the domain as:
nutkWallfunction for the walls
kqRWallFunction for k
epsilonWallFunction for eps
no-slip condition for U
zeroGradient for p

I also changed the inlet conditions of the ABL using Uref= 4 m/s, Href = 0,25m.

I am thankful for any comments regarding my thoughts about the simulations or hints why my denser mesh doesn't work.
Bazinga is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Enhanced Wall Treatment paduchev FLUENT 21 August 14, 2017 00:54
Grid deformation for high aspect ratio grids praveen SU2 Shape Design 5 February 15, 2013 20:50
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
how to understand high resolution scheme and high order scheme iilw1314 Main CFD Forum 3 September 30, 2010 10:43
Grid refinement with wall functions Mark Render Main CFD Forum 8 May 2, 2000 07:09

All times are GMT -4. The time now is 00:47.