CFD中文网

@李东岳 , 李老师您好，我已经在chtMultiRegionFoam中植入了第一张图所示的拉普拉斯及泊松方程，

实现了对电势在流固区域的耦合求解，但电荷的求解似乎不对，然后电场力的求解结果也和参考文章里面差的太多，但我边界设置理论上是符合文章要求的，第二张图和第三张图是我和参考文章计算结果的对比。参考文章及我的算例和编译的求解器链接放在这里：https://pan.baidu.com/s/1U8NBBon8CFVvDoF9Fsi9UA 提取码：mppz
麻烦李老师能抽空看一下，检查一下问题，不胜感激。一个头俩大

@李东岳 李老师您好，按照您的提示我把div(grad())项改成了laplacian操作符，但是在单独求解图1所示的方程时迭代了几步之后它就挂了，图2 是我的代码，其中epsilonr、ka均是常数，rhoel定义的场是volScalarField,不知是不是自己代码编写的问题，还望李老师指点迷津

@李东岳 好的，感谢李老师

@李东岳 感谢李老师回复 ，问题解决，OF小白一只，多多包涵

@东岳，李老师您好，目前自己在做气动热这块，就想求解第一张图上放的这个方程，自己编译的代码如第二张图所示，epsilonr和rhoel定义的场都是volScalarField,但编译的时候报这种错

@李东岳 在 关于cht的网格划分问题 中说：

writeZones

万分感谢李老师的解答，转化格式（写入命令语句）为
fluentMeshToFoam xxx.msh -writeZones，
本人OpenFoam小白，望老师多多包涵

@李东岳 李老师您好，我的网格是ICEM画的，然后再导入OpenFoam的，网格分了两个区，为什么导入之后PolyMesh里面看不到任何关于网格区域的信息，如果我想在其中一个区域添加源项该如何调用这个区域呢

@我是河滩 在 指定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 ??:?
浮点数例外 (核心已转储)