[robertomzc]

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)

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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

[robertomzc]

The problem relied in the communication between PointWise and OpenFOAM.
For now the workaround that makes it work is running

 renumberMesh -overwrite

before running the case

[cosmictrex]

Originally Posted by robertomzc

The problem relied in the communication between PointWise and OpenFOAM.
For now the workaround that makes it work is running

 renumberMesh -overwrite

before running the case

Thank you so very much for this. I have been stuck on this all day. May I ask how you figured this out and what exactly is the source of miscommunication between PW and OpenFOAM?