CFD中文网

仿照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求解器对网格有什么要求？激波位置为什么会一直变化？