CFD中文网

求助，目前在使用hy2Foam求解器计算类似乘波体外形，已经计算过钝头算例可以得到收敛解。但目前计算外形相对复杂，二维计算，计算域为1*0.55m，网格数量42000。计算过程中经常报错，求问有没有什么调试建议。下面是报错信息。

Mean and max Courant Numbers = 5.29417e-05 0.0491796
deltaT = 3.99074e-11
Time = 2.017990284e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for N2, Initial residual = 0.000161695, Final residual = 3.0494e-13, No Iterations 1
DILUPBiCG:  Solving for O2, Initial residual = 0.000148597, Final residual = 1.8029e-13, No Iterations 1
DILUPBiCG:  Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 1.07385e-06, Final residual = 1.49309e-06, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 8.83461e-07, Final residual = 1.19425e-06, No Iterations 1000
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 2.94288e-05, Final residual = 3.15326e+222, No Iterations 400
Attempt to use rho2ReactionThermo out of temperature range 11152 times during this iteration.
	-1.64159e+236 < Tlow: 200
Phase no 1.0  ExecutionTime = 3564.14 s  ClockTime = 3564 s  Iteration no 4641 (0.79 s)

Mean and max Courant Numbers = 5.41311e-05 0.0489379
deltaT = 4.07806e-11
Time = 2.018031064e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for N2, Initial residual = 0.000162296, Final residual = 1.67236e-16, No Iterations 2
DILUPBiCG:  Solving for O2, Initial residual = 0.000148383, Final residual = 6.70388e-18, No Iterations 2
DILUPBiCG:  Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
#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::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#8  Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:?
#9  ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
#10  ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
#11  ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"

@biaowai 在 使用mapFields 操作更改网格计算出错 中说：

印象可以修改差值方式或者用mapFieldsPar试试

我用的1706版本，尝试过修改差值方式还是不行。去试试mapFieldsPar。

@学流体的小明 确实是数据插值有问题，手动把右下角网格点数据改掉了，大概两行，一百多个点，现在可以跑起来了

计算域不变，边界条件不变，将12000个网格加密到15000个后，使用 mapFields -consistent ../sourceFolder 命令进行结果映射的时候，映射之后的场是这个样子，感觉分布不合理，是不是 mapFields 操作有问题。

使用hy2Foam求解器做RamC模型验证，在驻点线画图时发现异常，组分质量分数结果使用线性坐标为正常曲线，如下所示。
但当转换为指数坐标时，曲线出现奇怪波动，请问是为什么呢？

但计算过程中会出现CFL数突然增大之后报错的问题Mean and max Courant Numbers = 0.0196702 0.2
deltaT = 1.34917e-09
Time = 1.381229231e-07

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for N2, Initial residual = 0.102714, Final residual = 8.31435e-10, No Iterations 11
DILUPBiCG: Solving for O2, Initial residual = 0.0998228, Final residual = 2.97854e-08, No Iterations 11
DILUPBiCG: Solving for NO, Initial residual = 0.00103043, Final residual = 4.54498e-07, No Iterations 11
DILUPBiCG: Solving for N2+, Initial residual = 0.157236, Final residual = 1.27729e-09, No Iterations 11
DILUPBiCG: Solving for O2+, Initial residual = 0.157241, Final residual = 1.44374e-09, No Iterations 11
DILUPBiCG: Solving for NO+, Initial residual = 0.157279, Final residual = 1.20494e-09, No Iterations 11
DILUPBiCG: Solving for N, Initial residual = 0.00643012, Final residual = 4.19377e-07, No Iterations 11
DILUPBiCG: Solving for O, Initial residual = 0.00193657, Final residual = 2.2919e-07, No Iterations 11
DILUPBiCG: Solving for N+, Initial residual = 0.157174, Final residual = 1.12526e-09, No Iterations 11
DILUPBiCG: Solving for O+, Initial residual = 0.156703, Final residual = 1.23372e-08, No Iterations 11
DILUPBiCG: Solving for e-, Initial residual = 0.0591289, Final residual = 2.85246e-15, No Iterations 5
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 5.41102e-05, Final residual = 5.24963e-09, No Iterations 120
smoothSolver: Solving for Uy, Initial residual = 0.00128686, Final residual = 1.27448e-07, No Iterations 160
smoothSolver: Solving for Uz, Initial residual = 0.000139985, Final residual = 1.28962e-08, No Iterations 70
diagonal: Solving for rhoEv, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for hevel, Initial residual = 0.000257938, Final residual = 3.98118e-07, No Iterations 400
smoothSolver: Solving for e, Initial residual = 0.00686842, Final residual = 6.70799e-06, No Iterations 140
Attempt to use rho2ReactionThermo out of temperature range 139 times during this iteration.
 Thigh: 40000 < 640000
Phase no 1.0 ExecutionTime = 4.97 s ClockTime = 5 s Iteration no 112 (0.04 s)

Mean and max Courant Numbers = 0.0196515 0.552368
deltaT = 4.88502e-10
Time = 1.386114254e-07

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for N2, Initial residual = 0.0831271, Final residual = 4.94891e-09, No Iterations 11
DILUPBiCG: Solving for O2, Initial residual = 0.0830444, Final residual = 6.0302e-09, No Iterations 11
DILUPBiCG: Solving for NO, Initial residual = 0.0225179, Final residual = 4.10796e-07, No Iterations 11
DILUPBiCG: Solving for N2+, Initial residual = 0.070932, Final residual = 1.23497e-08, No Iterations 11
DILUPBiCG: Solving for O2+, Initial residual = 0.0707959, Final residual = 1.16027e-08, No Iterations 11
DILUPBiCG: Solving for NO+, Initial residual = 0.0712528, Final residual = 1.38649e-08, No Iterations 11
DILUPBiCG: Solving for N, Initial residual = 0.117227, Final residual = 8.58776e-08, No Iterations 11
DILUPBiCG: Solving for O, Initial residual = 0.209694, Final residual = 8.63459e-07, No Iterations 11
DILUPBiCG: Solving for N+, Initial residual = 0.240723, Final residual = 1.5458e-05, No Iterations 11
DILUPBiCG: Solving for O+, Initial residual = 0.0709349, Final residual = 3.14266e-08, No Iterations 11
DILUPBiCG: Solving for e-, Initial residual = 0.0330941, Final residual = 7.94061e-14, No Iterations 3
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1.89119e-05, Final residual = 1.65027e-09, No Iterations 70
smoothSolver: Solving for Uy, Initial residual = 0.000264353, Final residual = 2.41377e-08, No Iterations 102
smoothSolver: Solving for Uz, Initial residual = 5.29899e-05, Final residual = 5.01779e-09, No Iterations 70
diagonal: Solving for rhoEv, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for hevel, Initial residual = 0.278772, Final residual = 249.287, No Iterations 400
smoothSolver: Solving for e, Initial residual = 0.00510102, Final residual = 1.56626e+32, No Iterations 400
Attempt to use rho2ReactionThermo out of temperature range 3106 times during this iteration.
-2.31822e+10 < Tlow: 200
 Thigh: 40000 < 3.66402e+10
Phase no 1.0 ExecutionTime = 5.02 s ClockTime = 5 s Iteration no 113 (0.05 s)

[10] #0 Foam::error::printStack(Foam::Ostream&)[9] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
at ??:?
[9] #1 Foam::sigFpe::sigHandler(int)[10] #1 Foam::sigFpe::sigHandler(int) at ??:?
[9] #2 ? at ??:?
[10] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[10] #3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[9] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[9] #5 at ??:?
[10] #4 ?Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[9] #6 in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[10] #5 ?? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[9] #7 __libc_start_main in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[10] #6 in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #8 ?? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[10] #7 __libc_start_main in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[DESKTOP-N36BFE3:23240] *** Process received signal ***
[DESKTOP-N36BFE3:23240] Signal: Floating point exception (8)
[DESKTOP-N36BFE3:23240] Signal code: (-6)
[DESKTOP-N36BFE3:23240] Failing at address: 0x3e800005ac8
[DESKTOP-N36BFE3:23240] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x43090)[0x7f49b1370090]
[DESKTOP-N36BFE3:23240] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f49b137000b]
[DESKTOP-N36BFE3:23240] [ 2] in "/lib/x86_64-linux-gnu/libc./lib/x86_64-linux-gnu/libc.so.6(+0x43090)[0x7f49b1370090]
[DESKTOP-N36BFE3:23240] [ 3] so.6"
[10] #8 /home/sy/OpenFOAM/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x3a)[0x7f49b1f24b7a]
[DESKTOP-N36BFE3:23240] [ 4] hy2Foam(ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8+0x174)[0x55dafcc83f74]
[DESKTOP-N36BFE3:23240] [ 5] hy2Foam(+0x3fed5)[0x55dafcc10ed5]
[DESKTOP-N36BFE3:23240] [ 6] hy2Foam(+0x3464c)[0x55dafcc0564c]
[DESKTOP-N36BFE3:23240] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f49b1351083]
[DESKTOP-N36BFE3:23240] [ 8] hy2Foam(+0x35f1e)[0x55dafcc06f1e]
[DESKTOP-N36BFE3:23240] *** End of error message ***
? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[DESKTOP-N36BFE3:23241] *** Process received signal ***
[DESKTOP-N36BFE3:23241] Signal: Floating point exception (8)
[DESKTOP-N36BFE3:23241] Signal code: (-6)
[DESKTOP-N36BFE3:23241] Failing at address: 0x3e800005ac9
[DESKTOP-N36BFE3:23241] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x43090)[0x7fb8d8d28090]
[DESKTOP-N36BFE3:23241] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7fb8d8d2800b]
[DESKTOP-N36BFE3:23241] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x43090)[0x7fb8d8d28090]
[DESKTOP-N36BFE3:23241] [ 3] /home/sy/OpenFOAM/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x3a)[0x7fb8d98dcb7a]
[DESKTOP-N36BFE3:23241] [ 4] hy2Foam(ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8+0x174)[0x55f634616f74]
[DESKTOP-N36BFE3:23241] [ 5] hy2Foam(+0x3fed5)[0x55f6345a3ed5]
[DESKTOP-N36BFE3:23241] [ 6] hy2Foam(+0x3464c)[0x55f63459864c]
[DESKTOP-N36BFE3:23241] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7fb8d8d09083]
[DESKTOP-N36BFE3:23241] [ 8] hy2Foam(+0x35f1e)[0x55f634599f1e]
[DESKTOP-N36BFE3:23241] *** End of error message ***

Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.

mpirun noticed that process rank 9 with PID 0 on node DESKTOP-N36BFE3 exited on signal 8 (Floating point exception).

已经解决了，需要在controlDic中修改输出精度，至少在12以上

使用pointwise输出二维网格，网格检查结果良好
但是用
makeAxialMesh -case ~/OpenFOAM/sy-v1706/run/00mesh -axis symmetry -wedge BaseAndTop
collapseEdges指令生成楔形网格后检查网格质量，正交性变差，导致计算不断报错。

想请教一下，相同算例直接计算可以得到收敛结果，但并行运算时算到大概一半时间时报错。下面是decomposeParDict文件和报错信息。
报错信息.txt

@ShaneHEEE 设置的是从x方向流入 我再检查一下

@李东岳 好的，谢谢老师，我再检查看看

已经计算至0.001s，按照我设置的初始速度应该已经走过一半距离，但现在数值还是只存在与第一层网格。有点搞不懂到底是哪个部分出现了问题

不是标准算例 是自己建的模型 是可压缩流动

使用hy2Foam求解器计算三维翼型流场过程中出现内场维持初值不变的情况，建立一个简单的长方体算例出现相同情况，仅入口两层网格存在数值，下面是计算至0.001s的情况 采用100x10x10的结构化网格，下面是入口，上面是出口，四周设置为零梯度条件。

仿照https://www.cfd-online.com/Forums/openfoam-solving/60471-call-fortran-program-openfoam.html进行Fortran程序的调用，编译过程中报错： from chemistryModel/rho2ChemistryModel/rho2ChemistryModels.C:37:
/home/sy/OpenFOAM/sy-v1706/src/thermophysicalModels/strath/strathSpecie/lnInclude/CEATransportI.H:152:7: error: expected constructor, destructor, or type conversion before ‘(’ token
phypro(&TG, &TE, &PT, &PE, &CON[0], &CP_FINAL[0][0], &H_FINAL[0][0], &MIU_FINAL, &K_FINAL[0]);
但仿照https://www.cnblogs.com/snake553/p/6962386.html编写C++程序可以实现Fortran程序的调用

增加了一种读取表格计算物性的方法，编译无报错但实际应用时显示下列报错信息

@Elysiumor 在 使用hy2Foam求解器，更换bluntedCone算例网格便计算报错 中说：

@sy2315317 请问您按照官网操作，使用gnuplot可以实时查看残差吗

不能，提示缺了一个库，但装不上

@Elysiumor 直接运算可以得出结果，没试过单线程，并行运算15小时左右

@Elysiumor 我也是照着官网教程安装的~

hy2Foam求解器的bluntedCone算例使用12万网格进行计算，计算时间较长。尝试使用pointwise绘制简化网格进行计算，分别使用400x150、300x100、150x50网格进行计算，正交性较好，均报错。下面是错误信息：
Foam::rarefied<Foam::powerLawEuckenTransport<Foam::species::multiThermo<Foam::decoupledEnergyModesThermo<Foam::perfect2GasFoam::advancedSpecie >, Foam::sensible2InternalEnergy> > >::correct(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[1] #4 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[1] #5 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
[1] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
使用正交性较差的60x100的网格结果可以计算，但激波位置一直发生变化。
想问hy2Foam求解器对网格有什么要求？激波位置为什么会一直变化？