~$ run
run$ cp -r $FOAM_TUTORIALS/compressible/rhoPimpleFoam/laminar/forwardStep/ pimpleForwardStep
run$ cd pimpleForwardStep
pimpleForwardStep$ blockMesh
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10-b2f0c22c5002
Exec   : blockMesh
Date   : Sep 19 2022
Time   : 11:19:38
Host   : "VM"
PID    : 2724
I/O    : uncollated
Case   : /home/vagrant/OpenFOAM/vagrant-10/run/pimpleForwardStep
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Reading "blockMeshDict"

Creating block mesh from
    "system/blockMeshDict"
No non-linear block edges defined
No non-planar block faces defined
Creating topology blocks
Creating topology patches

Creating block mesh topology

Check topology

	Basic statistics
		Number of internal faces : 2
		Number of boundary faces : 14
		Number of defined boundary faces : 14
		Number of undefined boundary faces : 0
	Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1
    Block 0 cell size :
        i : 0.024 .. 0.024
        j : 0.02 .. 0.02
        k : 0.1
    Block 1 cell size :
        i : 0.024 .. 0.024
        j : 0.02 .. 0.02
        k : 0.1
    Block 2 cell size :
        i : 0.024 .. 0.024
        j : 0.02 .. 0.02
        k : 0.1

There are no merge patch pairs edges

Writing polyMesh
----------------
Mesh Information
----------------
  boundingBox: (0 0 -0.05) (3 1 0.05)
  nPoints: 10852
  nCells: 5250
  nFaces: 21175
  nInternalFaces: 10325
----------------
Patches
----------------
  patch 0 (start: 10325 size: 50) name: inlet
  patch 1 (start: 10375 size: 40) name: outlet
  patch 2 (start: 10415 size: 25) name: bottom
  patch 3 (start: 10440 size: 125) name: top
  patch 4 (start: 10565 size: 110) name: obstacle
  patch 5 (start: 10675 size: 10500) name: defaultFaces

End

pimpleForwardStep$ checkMesh
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10-b2f0c22c5002
Exec   : checkMesh
Date   : Sep 19 2022
Time   : 12:39:27
Host   : "cfdVM"
PID    : 3555
I/O    : uncollated
Case   : /home/vagrant/OpenFOAM/vagrant-10/run/pfs
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create polyMesh for time = 0

Time = 0s

Mesh stats
    points:           10852
    internal points:  0
    faces:            21175
    internal faces:   10325
    cells:            5250
    faces per cell:   6
    boundary patches: 6
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     5250
    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.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    inlet               50       102      ok (non-closed singly connected)  
    outlet              40       82       ok (non-closed singly connected)  
    bottom              25       52       ok (non-closed singly connected)  
    top                 125      252      ok (non-closed singly connected)  
    obstacle            110      222      ok (non-closed singly connected)  
    defaultFaces        10500    10852    ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (0 0 -0.05) (3 1 0.05)
    Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
    Mesh has 2 solution (non-empty) directions (1 1 0)
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (5.04971e-18 1.13618e-17 -1.80824e-16) OK.
    Max cell openness = 1.12938e-16 OK.
    Max aspect ratio = 1.2 OK.
    Minimum face area = 0.00048. Maximum face area = 0.0024.  Face area magnitudes OK.
    Min volume = 4.8e-05. Max volume = 4.8e-05.  Total volume = 0.252.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.09978e-13 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

pimpleForwardStep$ rhoPimpleFoam
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10-b2f0c22c5002
Exec   : rhoPimpleFoam
Date   : Sep 19 2022
Time   : 11:19:56
Host   : "VM"
PID    : 2725
I/O    : uncollated
Case   : /home/vagrant/OpenFOAM/vagrant-10/run/pimpleForwardStep
nProcs : 1
sigFpe : Enabling floating point exception trapping (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


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calculations will do 2 corrections


PIMPLE: Operating solver in transient mode with 2 outer correctors


Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating thermophysical transport model

Selecting thermophysical transport type laminar
Selecting default laminar thermophysical transport model unityLewisFourier
Creating field dpdt

Creating field kinetic energy K

No MRF models present

No fvModels present
No fvConstraints present
Courant Number mean: 0.249762 max: 0.25

Starting time loop

Courant Number mean: 0.249762 max: 0.25
Time = 0.002s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.835954, Final residual = 0.303956, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.992372, Final residual = 7.00073e-17, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 1, Final residual = 3.99937e-06, No Iterations 3
smoothSolver:  Solving for p, Initial residual = 1, Final residual = 1.29382e-07, No Iterations 4
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.95445e-11, global = 4.59826e-11, cumulative = 4.59826e-11
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.029104, Final residual = 2.88103e-10, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0365878, Final residual = 1.60986e-08, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0964448, Final residual = 1.72866e-08, No Iterations 1
smoothSolver:  Solving for p, Initial residual = 0.0216961, Final residual = 3.69034e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.10126e-10, global = 1.09875e-10, cumulative = 1.55858e-10
ExecutionTime = 0.203425 s  ClockTime = 0 s

Courant Number mean: 0.249764 max: 0.270674
Time = 0.004s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.185341, Final residual = 2.40791e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.191307, Final residual = 1.41438e-07, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.17141, Final residual = 1.09114e-06, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 0.192614, Final residual = 1.40018e-07, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.22113e-10, global = 5.2066e-10, cumulative = 6.76518e-10
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.0113229, Final residual = 3.42484e-10, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.010938, Final residual = 7.62427e-09, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0504379, Final residual = 3.25284e-06, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 0.0234216, Final residual = 1.37557e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.28819e-11, global = 6.4928e-11, cumulative = 7.41446e-10
ExecutionTime = 0.215673 s  ClockTime = 0 s

Courant Number mean: 0.249779 max: 0.295615
Time = 0.006s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.135168, Final residual = 3.93174e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.1365, Final residual = 2.90439e-07, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.116595, Final residual = 6.98715e-06, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 0.121759, Final residual = 8.39205e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.89575e-10, global = 5.7584e-10, cumulative = 1.31729e-09
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.00754074, Final residual = 1.54294e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00709415, Final residual = 1.5106e-08, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0379848, Final residual = 6.52938e-08, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 0.015404, Final residual = 1.14999e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.10532e-10, global = 1.0552e-10, cumulative = 1.42281e-09
ExecutionTime = 0.228095 s  ClockTime = 1 s

Courant Number mean: 0.249817 max: 0.330713
Time = 0.008s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.109998, Final residual = 9.61884e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.108976, Final residual = 6.31217e-07, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0980622, Final residual = 9.65876e-06, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 0.0935742, Final residual = 6.04758e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.48988e-10, global = 5.94325e-10, cumulative = 2.01713e-09
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.00605867, Final residual = 4.89571e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00567702, Final residual = 2.96264e-08, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0424891, Final residual = 9.99417e-06, No Iterations 1
smoothSolver:  Solving for p, Initial residual = 0.0106971, Final residual = 1.1627e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.58234e-10, global = 1.50284e-10, cumulative = 2.16742e-09
ExecutionTime = 0.240562 s  ClockTime = 1 s

Courant Number mean: 0.249888 max: 0.383413
Time = 0.01s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.0949074, Final residual = 1.95134e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0917156, Final residual = 1.24199e-06, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0940598, Final residual = 1.95416e-07, No Iterations 3
smoothSolver:  Solving for p, Initial residual = 0.0800918, Final residual = 5.25456e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.83051e-10, global = 6.5904e-10, cumulative = 2.82646e-09
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.00560314, Final residual = 1.49972e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00517889, Final residual = 6.07081e-08, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0597194, Final residual = 3.11597e-08, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 0.00739937, Final residual = 1.23015e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.22097e-10, global = 2.09978e-10, cumulative = 3.03643e-09
ExecutionTime = 0.252256 s  ClockTime = 1 s

Courant Number mean: 0.250026 max: 0.472331
Time = 0.012s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.085275, Final residual = 3.64692e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0796627, Final residual = 2.26134e-06, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0954682, Final residual = 3.37051e-07, No Iterations 3
smoothSolver:  Solving for p, Initial residual = 0.0725467, Final residual = 5.90646e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.15919e-09, global = 9.37893e-10, cumulative = 3.97433e-09
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.00574302, Final residual = 4.45068e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00505768, Final residual = 1.34611e-07, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.0796528, Final residual = 9.6544e-06, No Iterations 1
smoothSolver:  Solving for p, Initial residual = 0.00686812, Final residual = 1.35242e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.10289e-10, global = 2.82329e-10, cumulative = 4.25666e-09
ExecutionTime = 0.264436 s  ClockTime = 1 s

Courant Number mean: 0.25033 max: 0.650559
Time = 0.014s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.0792685, Final residual = 6.64422e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0707776, Final residual = 4.08768e-06, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.109434, Final residual = 8.11861e-07, No Iterations 3
smoothSolver:  Solving for p, Initial residual = 0.0703919, Final residual = 5.85818e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.48671e-09, global = 1.01439e-09, cumulative = 5.27104e-09
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.0070433, Final residual = 1.61834e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00532753, Final residual = 3.49144e-07, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.117347, Final residual = 2.9863e-06, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 0.00646344, Final residual = 1.12867e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.17861e-10, global = 2.89641e-10, cumulative = 5.56068e-09
ExecutionTime = 0.275743 s  ClockTime = 1 s

Courant Number mean: 0.251275 max: 1.17915
Time = 0.016s

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.0764808, Final residual = 1.29e-06, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.064002, Final residual = 7.94223e-06, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.26368, Final residual = 4.9006e-06, No Iterations 3
smoothSolver:  Solving for p, Initial residual = 0.0715941, Final residual = 6.77228e-08, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.94255e-09, global = 1.06042e-09, cumulative = 6.6211e-09
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.00795199, Final residual = 9.99449e-11, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.021992, Final residual = 1.56364e-09, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.285134, Final residual = 2.3536e-08, No Iterations 2


--> FOAM FATAL ERROR: 
Negative initial temperature T0: -3.09503

    From function static Foam::scalar Foam::species::thermo<Thermo, Type>::T(const ThermoType&, Foam::scalar, Foam::scalar, Foam::scalar, FType, dFdTType, LimitType, bool) [with ThermoType = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>; FType = double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const; dFdTType = double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const; LimitType = double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double) const; Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 299.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  double Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::T<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double) const>(Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> const&, double, double, double, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double, double) const, double (Foam::hConstThermo<Foam::perfectGas<Foam::specie> >::*)(double) const, bool) at ??:?
#3  Foam::hePsiThermo<Foam::psiThermo::composite, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4  virtual thunk to Foam::hePsiThermo<Foam::psiThermo::composite, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5  ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#6  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
Aborted (core dumped)
pimpleForwardStep$