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

simpleFoam Floating point error after applying refineMesh

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

LinkBack Thread Tools Display Modes
Old   April 28, 2014, 22:21
Exclamation simpleFoam Floating point error after applying refineMesh
New Member
CoolHand's Avatar
Join Date: Jan 2014
Location: South Florida
Posts: 19
Rep Power: 5
CoolHand is on a distinguished road
I'm posting this to Solving because I think it's a problem with Solver/Mesh interaction. Please let me know if its not the right place for it.

I have a 2D simulation of flow over a half cylinder. It works perfectly until I use refineMesh on it. After using refineMesh, checkMesh reports:

***Number of edges not aligned with or perpendicular to non-empty directions: 1374
  <<Writing 2061 points on non-aligned edges to set nonAlignedEdges
and if I try to run the simulation it gives the error:

Starting time loop

Time = 1

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#2   in "/lib/x86_64-linux-gnu/"
#3   at symmTensorField.C:0
#4  Foam::inv(Foam::UList<Foam::SymmTensor<double> > const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#5  Foam::leastSquaresVectors::makeLeastSquaresVectors() const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#6  Foam::leastSquaresVectors::pVectors() const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#7  Foam::fv::leastSquaresGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#8  Foam::fv::cellLimitedGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#9  Foam::fv::gradScheme<Foam::Vector<double> >::grad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#10  Foam::tmp<Foam::GeometricField<Foam::outerProduct<Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#11  Foam::incompressible::RASModels::kEpsilon::divDevReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
#13  __libc_start_main in "/lib/x86_64-linux-gnu/"
 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)
The refinement goes across the bottom of the domain from the inlet to the outlet. If I reduce the refinement to be only 1 radius upstream and downstream of the cylinder I don't get the errors.

Here is my case file:

The mesh was created in the GMSH format by a MatLab code that I wrote. The refined mesh is in the 1 folder.

Any ideas?
CoolHand is offline   Reply With Quote


floating point exception, refinemesh, simplefoam

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
problem in toposet Ahmed Khattab OpenFOAM Pre-Processing 23 November 2, 2017 13:17
Floating point error and divergence detected aannjj FLUENT 0 July 2, 2013 03:44
sprayEngineFoam - floating point exception error ris OpenFOAM Running, Solving & CFD 3 May 27, 2013 07:01
No convergence; floating point exception all the time Tobi OpenFOAM 2 September 19, 2012 04:01

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