@我是河滩 在 指定cellzone区域添加动量源项的问题 中说：

遍历了整个网格区域进行搜索

前辈您好，您是如何遍历整个网格区域并获取网格坐标的呢，是利用mesh.C().component()语句吗

李老师您好，请问这个问题解决了吗？我也遇到类似的问题，就是只改湍流模型的名字，在编译的时候也出现了一些error

@Rui 我自己尝试过的办法有检查边界，更改离散格式和时间步长等等，换过求解器为reactingFoam.但都未能解决这个问题，希望大佬们能够指点，谢谢~

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.00441557, Final residual = 1.54512e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00624943, Final residual = 1.18395e-07, No Iterations 1
DILUPBiCG:  Solving for H2, Initial residual = 0.00188285, Final residual = 7.43415e-08, No Iterations 1
DILUPBiCG:  Solving for O2, Initial residual = 0.00188283, Final residual = 1.01514e-07, No Iterations 1
DILUPBiCG:  Solving for H, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for OH, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for H2O, Initial residual = 0.00188283, Final residual = 1.01514e-07, No Iterations 1
DILUPBiCG:  Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for h, Initial residual = 0.00147451, Final residual = 5.51534e-08, No Iterations 1
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 30.769
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = -214.865
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 147.685
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 147.687
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = -215.127
--> FOAM Warning : 
    From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
    attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 28.5007
min/max(T) = 200, 534.025
smoothSolver:  Solving for p, Initial residual = 0.00340835, Final residual = 7.54003e-11, No Iterations 1
smoothSolver:  Solving for p, Initial residual = 2.66055e-10, Final residual = 2.66055e-10, No Iterations 0
smoothSolver:  Solving for p, Initial residual = 2.66055e-10, Final residual = 2.66055e-10, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.58233e-12, global = 8.75409e-14, cumulative = 6.66904e-12
smoothSolver:  Solving for p, Initial residual = 2.46585e-05, Final residual = 5.53512e-13, No Iterations 1
smoothSolver:  Solving for p, Initial residual = 2.76545e-12, Final residual = 2.76545e-12, No Iterations 0
smoothSolver:  Solving for p, Initial residual = 2.76545e-12, Final residual = 2.76545e-12, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.65649e-14, global = 1.10285e-15, cumulative = 6.67015e-12
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  log in "/lib/x86_64-linux-gnu/libm.so.6"
#4  Foam::omegaWallFunctionFvPatchScalarField::calculate(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#5  Foam::omegaWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#6  Foam::omegaWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#7  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::updateCoeffs() at ??:?
#8  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:?
#9  ? at ??:?
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? at ??:?
浮点数例外 (核心已转储)