[colin47]

Hey everyone,


I've been dealing with a persistent error for over two weeks now and could use some help.


I'm working on simulating incompressible flow around a cylinder using pimpleFoam and LES, based on a mesh converted from Salome. I've decomposed the mesh with the scotch method in decomposePar. The cylinder surface is set as a wall, and the rear surfaces (top, bottom, and two sides) are set as cyclicAMI. I've varied the spanwise direction with 1, 2, 10, 25, and 50 layers, resulting in mesh sizes of 1, 2, 10, 25, and 50 million cells, respectively. Strangely, I've successfully run the simulations for 1, 2, and 10 layers without any issues. However, hitting 25 and 50 layers led to simulation crashes with the errors attached. Even after testing the mesh with URANS (k-epsilon) and switching to linear and first-order schemes, the errors persist. It's becoming increasingly perplexing.

Though the mesh is converted from Salome to OpenFOAM, after running checkMesh command it shows mesh is 'ok'.

The minimum surface area is 6e-7 and min. vlolume of cell is 1e-7...

Initially, I thought it might be due to insufficient CPUs or RAM. Despite scaling up from 4 nodes and 128 CPUs for the initial three cases to 32 nodes and 1024 CPUs for the 25-layer case, I'm still puzzled as to why these errors persist. I've attached the files from the 0, constant, system folders and the log file for your review. Any assistance would be immensely appreciated. Thank you in advance!


errors:

Starting time loop

Courant Number mean: 1.25088e-05 max: 0.108801
deltaT = 0.003
Time = 0.003

PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 0.999962, Final residual = 1.91444e-09, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 4.80386e-10, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 1.43133e-09, No Iterations 2
[210] #0 Foam::error:rintStack(Foam::Ostream&)[211] #0
[547] #1 Foam::sigFpe::sigHandler(int) at ??:?
at ??:?

......

[412] #2 ?[414] #2 ?/lib64/libpthread.so.0
at ??:?
[39] #2 ?[875] #3 at Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&)??:?
at in [126] #2 ???:?

.......


[c137-034:38070:0:38070] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
at ??:?
[89] #6 __libc_start_main at ??:?
[588] #6 __libc_start_main at ??:?
[127] #6 __libc_start_main at ??:?

......

= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= RANK 352 PID 243228 RUNNING AT c127-082
= KILLED BY SIGNAL: 8 (Floating point exception)

......

[colin47]

the related system files are attached