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

[Floating Point Exception] Canoncial Cylinder Case meshed with Pointwise

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 7, 2020, 09:13
Default [SOLVED] Floating Point Exception - Canonical Cylinder Case meshed with Pointwise
  #1
New Member
 
Roberto Martinez
Join Date: Mar 2020
Posts: 2
Rep Power: 0
robertomzc is on a distinguished road
Hi all,



I am trying to run a pretty canonical case of a cylinder attached to one wall (which I treat as symPlane in the simulation).

The unstructured mesh is done in Pointwise,with a Trex grid (Prism layer) around the cylinder of y+ ~ 80 (have tried also decreasing to 1 but same problem remained).

checkMesh is giving the okay and in general mesh looks good without any major parameter going to crazy values.

Solver is simpleFoam, with keSST model, basically copied from the simpleFoam/motorcycle tutorial case.

My guess is I am having some problem with the boundary conditions, because when I launch the simulation with simpleFoam I get directly a Floating Point Exception in the smoother (full log below)
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1912                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _f3950763fe-20191219 OPENFOAM=1912
Arch   : "LSB;label=32;scalar=64"
Exec   : simpleFoam
Date   : May 07 2020
Time   : 13:11:06
Host   : ubuntu
PID    : 34652
I/O    : uncollated
Case   : /home/rmartinez/OpenFOAM/rmartinez-v1912/run/cylinder
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: no convergence criteria found. Calculations will run for 500 steps.

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    RASModel        kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
    decayControl    false;
    kInf            0;
    omegaInf        0;
}

No MRF models present

No finite volume options present

Starting time loop

Time = 1

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libc.so.6
#3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#8  Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:?
#9  ? in ~/OpenFOAM/OpenFOAM-v1912/platforms/linux64GccDPInt32Opt/bin/simpleFoam
#10  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11  ? in ~/OpenFOAM/OpenFOAM-v1912/platforms/linux64GccDPInt32Opt/bin/simpleFoam
Floating point exception (core dumped)
I leave the case without the mesh attached as well in case you want to perform your own checks.


What could be causing this problem?



Thanks for your time,
Rober
Attached Files
File Type: zip cylinder_nomesh.zip (7.3 KB, 0 views)

Last edited by robertomzc; May 12, 2020 at 12:37.
robertomzc is offline   Reply With Quote

Old   May 12, 2020, 12:36
Default Solution
  #2
New Member
 
Roberto Martinez
Join Date: Mar 2020
Posts: 2
Rep Power: 0
robertomzc is on a distinguished road
The problem relied in the communication between PointWise and OpenFOAM.
For now the workaround that makes it work is running
Code:
 renumberMesh -overwrite
before running the case
robertomzc is offline   Reply With Quote

Reply

Tags
boundary conditions, floating point exception, simplefoam

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
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 15:53
[blockMesh] error EOF in blockMesh Ahmed Khattab OpenFOAM Meshing & Mesh Conversion 7 May 17, 2012 01:37
Ercoftac centrifugal pump case study bug, Openfoam 1.6-ext (Floating point exception) ArianeJasmin OpenFOAM Programming & Development 1 April 4, 2012 15:11
Problem in unstructured circular cylinder case Tarak OpenFOAM Running, Solving & CFD 3 January 5, 2011 04:13
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 13:24


All times are GMT -4. The time now is 17:56.