求教：Laval 喷管模拟产生激波边界条件设置问题

yingqing

想要设置喷管，入口总压102387.1Pa，出口静压42449.07Pa, 入口温度294.45K。采用kOmegaSST湍流模型。但是运行报告湍流边界条件设置有问题。改了好几次都不行，希望各位老师，能帮我看一下，看看是哪里设置存在问题。以下是我的初始条件：

U T p omega nut Ma k alphat
运行报错结果：

tsing@tsing:~/OpenFOAM/tsing-v2206/run/EosNozzlePerfect$ rhoPimpleFoam 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2206                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _76d719d1e6-20220624 OPENFOAM=2206 version=v2206
Arch   : "LSB;label=32;scalar=64"
Exec   : rhoPimpleFoam
Date   : Feb 24 2025
Time   : 18:29:00
Host   : tsing
PID    : 218467
I/O    : uncollated
Case   : /home/tsing/OpenFOAM/tsing-v2206/run/EosNozzlePerfect
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


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 kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    RASModel        kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
    decayControl    false;
    kInf            0;
    omegaInf        0;
}

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
DILUPBiCGStab:  Solving for Ux, Initial residual = 1, Final residual = 8.37036e-10, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 1, Final residual = 4.3098e-10, No Iterations 1
DILUPBiCGStab:  Solving for e, Initial residual = 1, Final residual = 4.75567e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00613491, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00680296, Final residual = 3.5161e-05, No Iterations 4
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.28493e-08, global = -7.8312e-09, cumulative = -7.8312e-09
GAMG:  Solving for p, Initial residual = 0.147144, Final residual = 0.000685025, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.000646506, Final residual = 5.05846e-07, No Iterations 5
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.88756e-09, global = -4.94043e-10, cumulative = -8.32524e-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::omegaWallFunctionFvPatchScalarField::calculate(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#4  Foam::omegaWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#5  Foam::omegaWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#6  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<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)

李东岳

换一个湍流模型能收敛么
用旧版本的能收敛么？比如OF9

yingqing

换成层流模型可以。但换成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)

李东岳

这两个问题都出现在壁面函数上目前智能判断到这里 OF9你有就试试

yingqing

OF9尝试了，还是避免函数有问题，结果如下：

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

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
Creating field dpdt

Creating field kinetic energy K

No MRF models present

No fvModels present
No fvConstraints present
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::nutkWallFunctionFvPatchScalarField::nut() const at ??:?
#4  Foam::nutWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#5  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/home/tsing/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#6  Foam::RASModels::kEpsilon<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correctNut() at ??:?
#7  ? in "/home/tsing/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  ? in "/home/tsing/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
Floating point exception (core dumped)

李东岳

能公开就把算例网格处理到3万以下，上传我debug

另外你可以看看这个 http://dyfluid.com/openfoam.html#laval

yingqing

网址上的案例看了，使用的是层流模型。我使用层流模型也可以跑通。但使用RAS模型加入壁面函数出问题，麻烦李老师给我看看，debug一下
EosNozzleTestUpdate.zip

李东岳

rhoPimpleFoam_2.tar.xz

CFD中文网

求教：Laval 喷管模拟产生激波边界条件设置问题