[ryanc6]

Hello all,

I am using OpenFOAM v1812 on a high performance computer. I am attempting to model a simple geometric supersonic projectile in a compressible fluid. I chose sonicFOAM as my solver. The projectile is cylindrical with a hemispherical nose shape. I have modeled it within Fusion360 and imported it using the blockMesh, surfaceFeatureExtract, and snappyHexMesh utilities. The mesh generation seems to work fine. As shown by the following output messages of those utilities.

blockmesh.out

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1812 OPENFOAM=1812
Arch   : "LSB;label=32;scalar=64"
Exec   : blockMesh
Date   : Jul 22 2019
Time   : 10:23:36
Host   : nid04231
PID    : 2373
I/O    : uncollated
Case   : /p/work/ryanc6/run/OpenFOAM/SDb
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

Creating block mesh from
    "/p/work/ryanc6/run/OpenFOAM/SDb/system/blockMeshDict"
Creating block edges
No non-planar block faces defined
Creating topology blocks
Creating topology patches

Creating block mesh topology

Check topology

        Basic statistics
                Number of internal faces : 0
                Number of boundary faces : 6
                Number of defined boundary faces : 6
                Number of undefined boundary faces : 0
        Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 0.0254
    Block 0 cell size :
        i : 0.00635 .. 0.00635
        j : 0.00846667 .. 0.00846667
        k : 0.00846667 .. 0.00846667


There are no merge patch pairs

Writing polyMesh with 0 cellZones
----------------
Mesh Information
----------------
  boundingBox: (-0.4572 -0.3048 -0.3048) (0.4572 0.3048 0.3048)
  nPoints: 772705
  nCells: 746496
  nFaces: 2265408
  nInternalFaces: 2213568
----------------
Patches
----------------
  patch 0 (start: 2213568 size: 46656) name: sides
  patch 1 (start: 2260224 size: 5184) name: inlet

End

*Note the snappyHexMesh output code was over 700 lines long so here is just the end segment saying it finished properly.

snappyHexMesh.out

No straight edges simplified and no points removed ...
Snapped mesh : cells:2834400  faces:8580516  points:2912371
Cells per refinement level:
    0   707296
    1   55552
    2   2063872
    3   3584
    4   4096
Writing mesh to time constant
Wrote mesh in = 16.11 s.
Mesh snapped in = 144.13 s.
Checking final mesh ...
Checking faces in error :
    non-orthogonality > 20  degrees                        : 137056
    faces with face pyramid volume < 1e-13                 : 0
    faces with face-decomposition tet quality < 1e+15      : 8443460
    faces with concavity > 50  degrees                     : 0
    faces with skewness > 4   (internal) or 20  (boundary) : 0
    faces with interpolation weights (0..1)  < 0.05        : 0
    faces with volume ratio of neighbour cells < 0.01      : 0
    faces with face twist < 0.02                           : 0
    faces on cells with determinant < 0.001                : 0
Finished meshing with 8580516 illegal faces (concave, zero area or negative cell pyramid volume)
Finished meshing in = 232.46 s.
End

I notice the error comes from the sonicFOAM utility. Here is a copy of the sonicFOAM output file:

sonicFOAM.out

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1812 OPENFOAM=1812
Arch   : "LSB;label=32;scalar=64"
Exec   : sonicFoam
Date   : Jul 22 2019
Time   : 10:28:18
Host   : nid04231
PID    : 2993
I/O    : uncollated
Case   : /p/work/ryanc6/run/OpenFOAM/SDb
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


PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating field kinetic energy K

No MRF models present

No finite volume options present

Starting time loop

Time = 0.2

Courant Number mean: 23317.8 max: 327563
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
Application 13672019 exit codes: 136
Application 13672019 resources: utime ~14s, stime ~20s, Rss ~2810192, inblocks ~0, outblocks ~5394888

Here is the error message generated:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libc.so.6
#3  Foam::symGaussSeidelSmoother::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::symGaussSeidelSmoother::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  ? at ??:?
#10  __libc_start_main in /lib64/libc.so.6
#11  ? at /home/abuild/rpmbuild/BUILD/glibc-2.19/csu/../sysdeps/x86_64/start.S:125
_pmiu_daemon(SIGCHLD): [NID 04231] [c8-1c0s1n3] [Mon Jul 22 10:28:52 2019] PE RANK 0 exit signal Floating point exception

I have no idea what this error message is trying to say and therefore no clue what is wrong with my code or how to fix it. Any help is appreciated, thanks. I also attached a copy of all the files. I appreciate the help!

ryanc6

[giovanni.medici]

Apparently the timestep you are using is way too high, the Courant number shall be at most < than 1 (in principle you could reach values slightly higher than 1 with pimple).
So I would suggest trying to reduce timestep. The sonicDyMFoam solver carries the possibility to adapt the timestep to a specified Courant value (my personal opinion is that this approach works only when some degree of convergence has been reached, otherwise, e.g. in the first part of the computation, it is very likely to obtain oscillations of timestep hence courant that may lead to divergence).
I would strongly suggest you to take a look to (and inspiration from) the tutorials for sonicFoam and sonicDyMFoam.