换成层流模型可以。但换成kEpsilon模型，也是在Epsilon壁面函数错误。输出结果如下所示。
PIMPLE: Operating solver in PISO mode

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            heRhoThermo;
    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 RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 4.30558e-05 max: 0.0192008

Starting time loop

Courant Number mean: 4.30558e-05 max: 0.0192008
Time = 5e-05

PIMPLE: iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 4.88632e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 2.83045e-08, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 1, Final residual = 3.20771e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00612066, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00680187, Final residual = 3.35179e-05, No Iterations 4
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.89239e-08, global = -7.08945e-09, cumulative = -7.08945e-09
GAMG:  Solving for p, Initial residual = 0.147326, Final residual = 0.000684668, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00064622, Final residual = 4.92284e-07, No Iterations 5
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.8365e-09, global = -5.67277e-10, cumulative = -7.65673e-09
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3  Foam::epsilonWallFunctionFvPatchScalarField::calculate(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#4  Foam::epsilonWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#5  Foam::epsilonWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#6  Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:?
#7  ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam
#8  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#9  ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam
Floating point exception (core dumped)
CFD中文网

yingqing
