[Edoardo1993]

Hello,
I am having a foam fatal error with reactingFoam, only when using a slightly finer grid. With coarse mesh it's totally fine. I don't know why min/max(T) diverges.. sometimes negative values or enormous values like 1e6 have occurred.

Courant Number mean: 0.001135298157 max: 0.9182607171
deltaT = 3.299359615e-05
Time = 0.044934

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
Chemistry time solved max/min : 0.5333121972 / 5.463038001e-06
DILUPBiCGStab:  Solving for H2, Initial residual = 7.293812575e-06, Final residual = 2.220911406e-08, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 7.290405096e-06, Final residual = 2.219849253e-08, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 7.225032898e-06, Final residual = 2.199064069e-08, No Iterations 1
min/max(T) = 200, 309.4691418
GAMG:  Solving for p, Initial residual = 0.03586210971, Final residual = 5.248603349e-05, No Iterations 2
GAMG:  Solving for p, Initial residual = 5.573778688e-05, Final residual = 2.5091275e-07, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.136252712e-12, global = 1.018235007e-13, cumulative = 4.010593461e-11
GAMG:  Solving for p, Initial residual = 0.0004663544021, Final residual = 3.278041203e-06, No Iterations 2
GAMG:  Solving for p, Initial residual = 3.288315574e-06, Final residual = 3.64851378e-09, No Iterations 4
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.659801872e-14, global = -1.50631397e-15, cumulative = 4.010442829e-11
PIMPLE: iteration 2
Chemistry time solved max/min : 0.5333121981 / 5.482100715e-06
DILUPBiCGStab:  Solving for H2, Initial residual = 5.971544133e-07, Final residual = 4.894734739e-09, No Iterations 1
DILUPBiCGStab:  Solving for O2, Initial residual = 5.974402205e-07, Final residual = 4.897210375e-09, No Iterations 1
DILUPBiCGStab:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab:  Solving for h, Initial residual = 6.013661987e-07, Final residual = 4.709978953e-09, No Iterations 1
[1] 
[1] 
[1] --> FOAM FATAL ERROR: (openfoam-2212 patch=230110)
[1] Maximum number of iterations exceeded: 100 when starting from T0:200 old T:1101.685303 new T:200 f:-296826.0785 p:101403.8258 tol:0.02
[1] 
[1]     From Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::janafThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::janafThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>]
[1]     in file ./src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.
[1] 
FOAM parallel run aborting
[1] 
[1] #0  Foam::error::printStack(Foam::Ostream&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
[1] #1  Foam::error::simpleExit(int, bool) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
[1] #2  ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
[1] #3  ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
[1] #4  ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
[1] #5  ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/reactingFoam
[1] #6  ? in /lib/x86_64-linux-gnu/libc.so.6
[1] #7  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[1] #8  ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/reactingFoam
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

This is checkMesh -allTopology -allGeometry:

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

Create mesh for time = 0

Enabling all (cell, face, edge, point) topology checks.

Enabling all geometry checks.

Time = 0

Mesh stats
    points:           262701
    faces:            762500
    internal faces:   737500
    cells:            250000
    faces per cell:   6
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       0

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

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Topological cell zip-up check OK.
    Face-face connectivity OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology Bounding box
    inlet               2484     2592     ok (non-closed singly connected)   (0 0 0) (0 1.5 1.5)
    outlet              2500     2601     ok (non-closed singly connected)   (3 0 0) (3 1.5 1.5)
    freestream          15000    15251    ok (non-closed singly connected)   (0 0 0) (3 1.5 1.5)
    walls               5000     5151     ok (non-closed singly connected)   (0 0 0) (3 0 1.5)
    inletFuel           16       25       ok (non-closed singly connected)   (0 0.7401806535 0.7401806535) (0 0.7598193465 0.7598193465)

Checking faceZone topology for multiply connected surfaces...
    No faceZones found.

Checking basic cellZone addressing...
    No cellZones found.

Checking basic pointZone addressing...
    No pointZones found.

Checking geometry...
    Overall domain bounding box (0 0 0) (3 1.5 1.5)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (1.253985221e-16 1.344603441e-16 9.136210307e-17) OK.
    Max cell openness = 3.297977928e-16 OK.
    Max aspect ratio = 20.41398827 OK.
    Minimum face area = 2.119361571e-05. Maximum face area = 0.008652924449.  Face area magnitudes OK.
    Min volume = 9.95876417e-08. Max volume = 0.0007967011336.  Total volume = 6.75.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.353287224e-13 OK.
    Coupled point location match (average 0) OK.
    Face tets OK.
    Min/max edge length = 0.004603652432 0.09397890673 OK.
    All angles in faces OK.
    Face flatness (1 = flat, 0 = butterfly) : min = 1  average = 1
    All face flatness OK.
    Cell determinant (wellposedness) : minimum: 7.183349046e-05 average: 0.3523742609
 ***Cells with small determinant (< 0.001) found, number of cells: 368
  <<Writing 368 under-determined cells to set underdeterminedCells
    Concave cell check OK.
    Face interpolation weight : minimum: 0.4688349092 average: 0.4771777174
    Face interpolation weight check OK.
    Face volume ratio : minimum: 0.8826538438 average: 0.913622635
    Face volume ratio check OK.

Failed 1 mesh checks.

End