用compressibleinterfoam好像能量方程有问题，经常出现低速和超音速都出现，求解决？

Liubo

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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 ??:?

bestucan

我猜是网格质量问题
https://cfd-china.com/topic/4357/checkmesh中发现的几个网格问题如何检查修改/4?_=1610591734371

参考这个帖子描述的，四两拨千斤，千斤拨四两

LYT

温度一直往负的走，不知道是这求解器自身有问题，还是我们边界条件给的不行

CFD中文网

用compressibleinterfoam好像能量方程有问题，经常出现低速和超音速都出现，求解决？