[yeharav]

Hi all,

I have installed foam-extend-5.0 on ubuntu 22.04 either with the .deb file, or compiled from sources (git branch 73c2a3d539c20e6b97cd625d5d79fb7f78149f5b)

The simpleFoam tutorials crush - I have tested motorbike, and pitzDaily.
The simpleFoam of the pimpleDynMFoam in the wingMotion tutorial also crushes.

I have attached below the output of the pitzDaily.
To run it:

>> blockMesh
>> simpleFoam

It crushes very fast.

Any ideas?


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
Build : 5.0-73c2a3d539c2
Exec : simpleFoam
Date : Feb 21 2023
Time : 19:58:56
Host : "lt-yehudaa"
PID : 658991
CtrlDict : "/home/yehudaa/Projects/2023/foamExtendTutorials/simpleFoam/pitzDaily/system/controlDict"
Case : /home/yehudaa/Projects/2023/foamExtendTutorials/simpleFoam/pitzDaily
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0


SIMPLE: convergence criteria
field p tolerance 0.0001
field U tolerance 0.0001
field "(k|epsilon)" tolerance 0.0001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
alphaEps 0.76923;
sigmaEps 1.3;
}


Starting time loop

Time = 1

BiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.00464911, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00960309, No Iterations 163
time step continuity errors : sum local = 0.475961, global = -0.00434833, cumulative = -0.00434833
BiCGStab: Solving for epsilon, Initial residual = 0.594985, Final residual = 0.00495, No Iterations 1
BiCGStab: Solving for k, Initial residual = 1, Final residual = 0.0031536, No Iterations 1
ExecutionTime = 0.07 s ClockTime = 0 s

Time = 2

BiCGStab: Solving for Ux, Initial residual = 0.561136, Final residual = 0.00840767, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0.980401, Final residual = 0.0158018, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.905008, Final residual = 0.00863013, No Iterations 57
time step continuity errors : sum local = 27.5255, global = 2.94505, cumulative = 2.9407
BiCGStab: Solving for epsilon, Initial residual = 0.999276, Final residual = 0.00800445, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.996406, Final residual = 0.000380816, No Iterations 1
ExecutionTime = 0.09 s ClockTime = 0 s

Time = 3

BiCGStab: Solving for Ux, Initial residual = 0.516662, Final residual = 0.00621829, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0.998397, Final residual = 0.0176732, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.999941, Final residual = 0.0089741, No Iterations 19
time step continuity errors : sum local = 53029.6, global = -41.905, cumulative = -38.9643
BiCGStab: Solving for epsilon, Initial residual = 1, Final residual = 0.00122703, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.999998, Final residual = 0.00573344, No Iterations 1
ExecutionTime = 0.09 s ClockTime = 0 s

Time = 4

BiCGStab: Solving for Ux, Initial residual = 0.550045, Final residual = 0.0056509, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0.999997, Final residual = 0.0160376, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00862799, No Iterations 58
time step continuity errors : sum local = 5.33155e+10, global = -1.23755e+08, cumulative = -1.23755e+08
BiCGStab: Solving for epsilon, Initial residual = 0.999972, Final residual = 0.00643736, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.999955, Final residual = 0.002415, No Iterations 1
ExecutionTime = 0.11 s ClockTime = 0 s

Time = 5

BiCGStab: Solving for Ux, Initial residual = 0.421288, Final residual = 0.00214245, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.00222639, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00733539, No Iterations 133
time step continuity errors : sum local = 1.12466e+19, global = 1.96772e+15, cumulative = 1.96772e+15
BiCGStab: Solving for epsilon, Initial residual = 0.999994, Final residual = 0.000797934, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.999988, Final residual = 0.0591909, No Iterations 1
ExecutionTime = 0.13 s ClockTime = 0 s

Time = 6

BiCGStab: Solving for Ux, Initial residual = 0.525405, Final residual = 0.00822925, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.059309, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00503789, No Iterations 5
time step continuity errors : sum local = 3.055e+39, global = 3.65535e+27, cumulative = 3.65535e+27
BiCGStab: Solving for epsilon, Initial residual = 1, Final residual = 0.000516864, No Iterations 1
BiCGStab: Solving for k, Initial residual = 1, Final residual = 2.5594e-09, No Iterations 1
ExecutionTime = 0.14 s ClockTime = 0 s

Time = 7

BiCGStab: Solving for Ux, Initial residual = 0.699943, Final residual = 4.29177e-11, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.0335931, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 2.86737e-09, No Iterations 2
time step continuity errors : sum local = 1.02702e+68, global = -2.47902e+51, cumulative = -2.47902e+51
Floating point exception (core dumped)

[orangeDuck]

I also am experiencing this. Were you able to solve it?

[amirhosseintaran]

Hello,

Did you follow the instructions in Allrun file ?
especially for motorbike case, you need to use snappyHexMesh in order to create your geometry.

Make sure to execute Allrun file by ./Allrun

Bests,
Amirhossein.

[PetarT]

I experience same problem, the pitzDaily crashes at 7th iteration, other solvers also fave similar problems. I noticed that all solvers have poor convergence it may be that there is a problem with numerical schemes. Does anyone has answer to this?


Update: running with openfoam2306 it converges successfully. It seems that the solver itself is broken.

[josemm92]

Hi,


I have found out the same issue with all turbulence tutorials. Only the tutorials with laminar flow work correctly.


Any help would be appreciated.

[kol1978]

Originally Posted by josemm92

Hi,


I have found out the same issue with all turbulence tutorials. Only the tutorials with laminar flow work correctly.


Any help would be appreciated.

I just installed it and this is a trial run - everything is the same as yours... HELP! SOS help SOS:сумасшедший::запутался:

kol@kol:~/foam/kol-5.0/run/pitzDaily$ fe50
kol@kol:~/foam/kol-5.0/run/pitzDaily$ blockMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
Build : 5.0-287705b4a589
Exec : blockMesh
Date : Jun 06 2024
Time : 00:26:38
Host : "kol"
PID : 7059
CtrlDict : "/home/kol/foam/kol-5.0/run/pitzDaily/system/controlDict"
Case : /home/kol/foam/kol-5.0/run/pitzDaily
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

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


Creating block mesh from
"/home/kol/foam/kol-5.0/run/pitzDaily/system/blockMeshDict"


Creating block mesh topology

Creating polyMesh from blockMesh

Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (-0.0206 -0.0254 -0.0005) (0.29 0.0254 0.0005)
nPoints: 25012
nCells: 12225
nFaces: 49180
nInternalFaces: 24170
----------------
Patches
----------------
patch 0 (start: 24170 size: 30) name: inlet
patch 1 (start: 24200 size: 57) name: outlet
patch 2 (start: 24257 size: 223) name: upperWall
patch 3 (start: 24480 size: 250) name: lowerWall
patch 4 (start: 24730 size: 24450) name: frontAndBack

End

kol@kol:~/foam/kol-5.0/run/pitzDaily$ checkMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
Build : 5.0-287705b4a589
Exec : checkMesh
Date : Jun 06 2024
Time : 00:27:33
Host : "kol"
PID : 7061
CtrlDict : "/home/kol/foam/kol-5.0/run/pitzDaily/system/controlDict"
Case : /home/kol/foam/kol-5.0/run/pitzDaily
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

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

Create polyMesh for time = 0

Time = 0

Mesh stats
all points: 25012
live points: 25012
all faces: 49180
live faces: 49180
internal faces: 24170
cells: 12225
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 12225
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Area [m^2] Surface topology
inlet 30 62 2.54e-05 ok (non-closed singly connected)
outlet 57 116 3.32e-05 ok (non-closed singly connected)
upperWall 223 448 0.00031106 ok (non-closed singly connected)
lowerWall 250 502 0.00033646 ok (non-closed singly connected)
frontAndBack 24450 25012 0.0290321 ok (non-closed singly connected)

Checking geometry...
This is a 2-D mesh
Overall domain bounding box (-0.0206 -0.0254 -0.0005) (0.29 0.0254 0.0005)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
Mesh (non-empty, non-wedge) dimensions 2
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (1.67338e-19 3.10465e-18 -5.10502e-15) Threshold = 1e-06 OK.
Max cell openness = 3.2572e-16 OK.
Max aspect ratio = 7.74923 OK.
Minumum face area = 1.65978e-07. Maximum face area = 5.14451e-06. Face area magnitudes OK.
Min volume = 1.65978e-10. Max volume = 4.08301e-09. Total volume = 1.4516e-05. Cell volumes OK.
Mesh non-orthogonality Max: 5.94745 average: 1.41505 Threshold = 70
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.260443 OK.

Mesh OK.

End

No VM guests are running outdated hypervisor (qemu) binaries on this host.
kol@kol:~/foam/kol-5.0/run/pitzDaily$ mc

Select an editor. To change later, run 'select-editor'.
1. /bin/nano <---- easiest
2. /usr/bin/mcedit
3. /usr/bin/emacs

Choose 1-3 [1]: 1

kol@kol:~/foam/kol-5.0/run/pitzDaily$ checkMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
Build : 5.0-287705b4a589
Exec : checkMesh
Date : Jun 06 2024
Time : 00:30:09
Host : "kol"
PID : 7248
CtrlDict : "/home/kol/foam/kol-5.0/run/pitzDaily/system/controlDict"
Case : /home/kol/foam/kol-5.0/run/pitzDaily
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

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

Create polyMesh for time = 0

Time = 0

Mesh stats
all points: 25012
live points: 25012
all faces: 49180
live faces: 49180
internal faces: 24170
cells: 12225
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 12225
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Area [m^2] Surface topology
inlet 30 62 2.54e-05 ok (non-closed singly connected)
outlet 57 116 3.32e-05 ok (non-closed singly connected)
upperWall 223 448 0.00031106 ok (non-closed singly connected)
lowerWall 250 502 0.00033646 ok (non-closed singly connected)
frontAndBack 24450 25012 0.0290321 ok (non-closed singly connected)

Checking geometry...
This is a 2-D mesh
Overall domain bounding box (-0.0206 -0.0254 -0.0005) (0.29 0.0254 0.0005)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
Mesh (non-empty, non-wedge) dimensions 2
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (1.67338e-19 3.10465e-18 -5.10502e-15) Threshold = 1e-06 OK.
Max cell openness = 3.2572e-16 OK.
Max aspect ratio = 7.74923 OK.
Minumum face area = 1.65978e-07. Maximum face area = 5.14451e-06. Face area magnitudes OK.
Min volume = 1.65978e-10. Max volume = 4.08301e-09. Total volume = 1.4516e-05. Cell volumes OK.
Mesh non-orthogonality Max: 5.94745 average: 1.41505 Threshold = 70
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.260443 OK.

Mesh OK.

End

kol@kol:~/foam/kol-5.0/run/pitzDaily$ mc

kol@kol:~/foam/kol-5.0/run/pitzDaily$ simpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 5.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
Build : 5.0-287705b4a589
Exec : simpleFoam
Date : Jun 06 2024
Time : 00:36:00
Host : "kol"
PID : 7330
CtrlDict : "/home/kol/foam/kol-5.0/run/pitzDaily/system/controlDict"
Case : /home/kol/foam/kol-5.0/run/pitzDaily
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0


SIMPLE: convergence criteria
field p tolerance 0.0001
field U tolerance 0.0001
field "(k|epsilon)" tolerance 0.0001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
alphaEps 0.76923;
sigmaEps 1.3;
}


Starting time loop

Time = 1

BiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.00464911, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00960309, No Iterations 163
time step continuity errors : sum local = 0.475961, global = -0.00434833, cumulative = -0.00434833
BiCGStab: Solving for epsilon, Initial residual = 0.594985, Final residual = 0.00495, No Iterations 1
BiCGStab: Solving for k, Initial residual = 1, Final residual = 0.0031536, No Iterations 1
ExecutionTime = 0.18 s ClockTime = 1 s

Time = 2

BiCGStab: Solving for Ux, Initial residual = 0.561136, Final residual = 0.00840767, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0.980401, Final residual = 0.0158018, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.905008, Final residual = 0.00863013, No Iterations 57
time step continuity errors : sum local = 27.5255, global = 2.94505, cumulative = 2.9407
BiCGStab: Solving for epsilon, Initial residual = 0.999276, Final residual = 0.00800445, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.996406, Final residual = 0.000380816, No Iterations 1
ExecutionTime = 0.21 s ClockTime = 1 s

Time = 3

BiCGStab: Solving for Ux, Initial residual = 0.516662, Final residual = 0.00621829, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0.998397, Final residual = 0.0176732, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.999941, Final residual = 0.0089741, No Iterations 19
time step continuity errors : sum local = 53029.6, global = -41.905, cumulative = -38.9643
BiCGStab: Solving for epsilon, Initial residual = 1, Final residual = 0.00122703, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.999998, Final residual = 0.00573344, No Iterations 1
ExecutionTime = 0.24 s ClockTime = 1 s

Time = 4

BiCGStab: Solving for Ux, Initial residual = 0.550045, Final residual = 0.0056509, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 0.999997, Final residual = 0.0160376, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00862799, No Iterations 58
time step continuity errors : sum local = 5.33155e+10, global = -1.23755e+08, cumulative = -1.23755e+08
BiCGStab: Solving for epsilon, Initial residual = 0.999972, Final residual = 0.00643736, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.999955, Final residual = 0.002415, No Iterations 1
ExecutionTime = 0.27 s ClockTime = 1 s

Time = 5

BiCGStab: Solving for Ux, Initial residual = 0.421288, Final residual = 0.00214245, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.00222639, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00733539, No Iterations 133
time step continuity errors : sum local = 1.12466e+19, global = 1.96772e+15, cumulative = 1.96772e+15
BiCGStab: Solving for epsilon, Initial residual = 0.999994, Final residual = 0.000797934, No Iterations 1
BiCGStab: Solving for k, Initial residual = 0.999988, Final residual = 0.0591909, No Iterations 1
ExecutionTime = 0.33 s ClockTime = 1 s

Time = 6

BiCGStab: Solving for Ux, Initial residual = 0.525405, Final residual = 0.00822925, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.059309, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00503789, No Iterations 5
time step continuity errors : sum local = 3.055e+39, global = 3.65535e+27, cumulative = 3.65535e+27
BiCGStab: Solving for epsilon, Initial residual = 1, Final residual = 0.000516864, No Iterations 1
BiCGStab: Solving for k, Initial residual = 1, Final residual = 2.5594e-09, No Iterations 1
ExecutionTime = 0.34 s ClockTime = 1 s

Time = 7

BiCGStab: Solving for Ux, Initial residual = 0.699943, Final residual = 4.29177e-11, No Iterations 1
BiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.0335931, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 2.86737e-09, No Iterations 2
time step continuity errors : sum local = 1.02702e+68, global = -2.47902e+51, cumulative = -2.47902e+51
Floating point exception (core dumped)
kol@kol:~/foam/kol-5.0/run/pitzDaily$

[DVSoares]

Hello,

Has anyone been able to find a solution to these crashes, or to identify if the errors were produced by the solver or by something going wrong at compilation?

Thank you

******* Full edit here *******

I gave some time to investigate and it seems the release 5.0 (main branch) is incompatible with GCC 13, the one that comes with Ubuntu 24.04.

To solve, use instead the nextRelease branch. The solvers were updated and should work normally without any different configuration. To do that, use this on your "foam" folder:

git clone --branch nextRelease https://git.code.sf.net/p/foam-extend/foam-extend-5.0 foam-extend-5.1

Check that the address remains '5.0' but the code was updated to '5.1'. Therefore, to source this branch of foam-extend correctly, add this line to your ~/.bashrc and adapt the correct path to your foam folder:

# foam-extend 5.1
 alias fe51='source $HOME/foam/foam-extend-5.1/etc/bashrc'

Proceed normally with compilation and test your tutorials. They should work normally.

[ThierryRochemont]

Switching to foam-extend-5.1 (nextRelease) solve the issue. Thanks DVSoares!