求教：Laval 喷管模拟产生激波边界条件设置问题
-
想要设置喷管，入口总压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) -
换成层流模型可以。但换成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尝试了，还是避免函数有问题，结果如下：
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) -
网址上的案例看了，使用的是层流模型。我使用层流模型也可以跑通。但使用RAS模型加入壁面函数出问题，麻烦李老师给我看看，debug一下
EosNozzleTestUpdate.zip
