/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.1-262087cdf8db
Exec : PBCavitatingFoam
Date : Jan 04 2021
Time : 21:20:13
Host : "mgt"
PID : 292963
Case : /home/test/liubo/liubo/PBC
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading g
PIMPLE: max iterations = 20
field "(U|p_rgh|p)" : relTol 0.01, tolerance 1e-06
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Constructing phaseChangeTwoPhaseMixture
Selecting phaseChange model SchnerrSauer
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleInternalEnergy;
}
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleInternalEnergy;
}
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}
leave from phaseChangeTwoPhaseMixture Constructor
Reading thermophysical properties
Calculating field g.h
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}
Creating field kinetic energy K
Courant Number mean: 0.02677338955 max: 0.2119866949
No finite volume options present
Starting time loop
Courant Number mean: 0.02677338955 max: 0.2119866949
Interface Courant Number mean: 0 max: 0
^o^ deltaT = 1e-05
Time = 1e-05
PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 0.9999999964
Min(alpha.water) = 0 Max(alpha.water) = 1
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = 0 Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 0.9999999927
Min(alpha.water) = 0 Max(alpha.water) = 1
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = 0 Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 0.999999989
Min(alpha.water) = 0 Max(alpha.water) = 1
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = 0 Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 0.9999999853
Min(alpha.water) = 0 Max(alpha.water) = 1
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = 0 Max(alpha.air) = 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 1.565205693e-09, No Iterations 5
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 1.423403593e-09, No Iterations 5
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 1.232277888e-09, No Iterations 5
smoothSolver: Solving for T, Initial residual = 0.9999963229, Final residual = 0.09144259942, No Iterations 6
min(T) 292.9992852
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 6085871.068, No Iterations 50
GAMG: Solving for p_rgh, Initial residual = 2.170075185e-15, Final residual = 2.170075185e-15, No Iterations 0
max(U) 4889963.244
min(p_rgh) -849012.5866
GAMG: Solving for p_rgh, Initial residual = 2.661173357e-10, Final residual = 2.661173357e-10, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.661173357e-10, Final residual = 2.661173357e-10, No Iterations 0
max(U) 1.142439092e+14
min(p_rgh) -849012.5866
GAMG: Solving for p_rgh, Initial residual = 2.676702599e-10, Final residual = 2.676702599e-10, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.676702599e-10, Final residual = 2.676702599e-10, No Iterations 0
max(U) 1.425091029e+14
min(p_rgh) -849012.5866
GAMG: Solving for p_rgh, Initial residual = 2.68040637e-10, Final residual = 2.68040637e-10, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.68040637e-10, Final residual = 2.68040637e-10, No Iterations 0
max(U) 1.489703193e+14
min(p_rgh) -849012.5866
GAMG: Solving for p_rgh, Initial residual = 2.681313485e-10, Final residual = 2.681313485e-10, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.681313485e-10, Final residual = 2.681313485e-10, No Iterations 0
max(U) 1.504634378e+14
min(p_rgh) -849012.5866
PIMPLE: iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 100322.0152
Min(alpha.water) = 0 Max(alpha.water) = 2.810378902e+11
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = -2.810378902e+11 Max(alpha.air) = 2.858846245e+11
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 3.059042308e+16
Min(alpha.water) = 0 Max(alpha.water) = 8.173458136e+22
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = -8.173458136e+22 Max(alpha.air) = 6.422063805e+21
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 2.539802609e+27
Min(alpha.water) = 0 Max(alpha.water) = 1.817814046e+34
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = -1.817814046e+34 Max(alpha.air) = 2.285753964e+34
MULES: Solving for alpha.air
MULES: Solving for alpha.vapour
Liquid phase volume fraction = 1.050092687e+39
Min(alpha.water) = 0 Max(alpha.water) = 6.532567778e+45
Min(alpha.vapour) = 0 Max(alpha.vapour) = 0
Min(alpha.air) = -6.532567778e+45 Max(alpha.air) = 5.600607821e+44
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.07551363127, Final residual = 4.375766351e-09, No Iterations 5
smoothSolver: Solving for Uy, Initial residual = 0.08313923201, Final residual = 8.360047038e-09, No Iterations 5
smoothSolver: Solving for Uz, Initial residual = 0.08809270175, Final residual = 8.908948216e-10, No Iterations 6
smoothSolver: Solving for T, Initial residual = 1, Final residual = 6.382659004e+119, No Iterations 1000
--> FOAM FATAL ERROR:
Maximum number of iterations exceeded
From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /opt/ohpc/pub/apps/OpenFOAM/OpenFOAM-2.3.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::TEs(double, double, double) const at ??:?
#3 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 Foam::threePhaseMixtureThermo::correct() at ??:?
#6
at ??:?
#7 __libc_start_main in "/lib64/libc.so.6"
#8
at ??:?
