~$ 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::T(const ThermoType&, Foam::scalar, Foam::scalar, Foam::scalar, FType, dFdTType, LimitType, bool) [with ThermoType = Foam::species::thermo >, Foam::sensibleInternalEnergy>; FType = double (Foam::hConstThermo >::*)(double, double) const; dFdTType = double (Foam::hConstThermo >::*)(double, double) const; LimitType = double (Foam::hConstThermo >::*)(double) const; Thermo = Foam::hConstThermo >; 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::sensibleInternalEnergy>::T >, Foam::sensibleInternalEnergy>, double (Foam::hConstThermo >::*)(double, double) const, double (Foam::hConstThermo >::*)(double, double) const, double (Foam::hConstThermo >::*)(double) const>(Foam::species::thermo >, Foam::sensibleInternalEnergy> const&, double, double, double, double (Foam::hConstThermo >::*)(double, double) const, double (Foam::hConstThermo >::*)(double, double) const, double (Foam::hConstThermo >::*)(double) const, bool) at ??:? #3 Foam::hePsiThermo >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:? #4 virtual thunk to Foam::hePsiThermo >, 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$