~$ 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$