[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.