CFD中文网

@东方白杨 我记得前面有话题讨论过，你看这个：https://bitbucket.org/lesituu/postchannelflow/src/master/

@李东岳 @欧阳 我计算ODW，一百万网格，氢氨机理s59r356，cpuLoad（redistributionInterval 20; maxImbalance 0.1;）和DLBFoam（只使用负载均衡，不使用pyJac）效果基本相当。

并行 sample 的还是串行的，我记得并行的时候会漏点

@wangfei9088 赞👍

@wangfei9088 再多问大佬一句，代码是open source的吗？我看到了CPC Library link to program files。我们用的求解器是基于OpenFOAM 12的内部求解器和v2312俄罗斯的一个求解器，都没考虑混合平均多组分扩散。感觉这对算缓燃还是挺重要的。

@wangfei9088 追踪了cpc的文章，看到名字，第一想到了你

@wangfei9088 歪个楼，大佬最近发没发CPC，看到个同名的

@欧阳 测试过，但是没继续使用，我现在用的是org的cpuLoad

@EricLiu 是纯气体

@EricLiu 你好，我也没用过vof那个求解器

也可以按照以前的方法写，PDRfoam你看下

@李东岳 李老师，还想请教下，一般你推荐使用binary还是ascii？binary好像存储更小，读取更快，但vim等直接没法看

@李东岳 是的李老师，这确实像是个bug。前几天我在OpenFOAM Issue Tracking提问了，可能没描述清楚，直接给我删了。刚刚debug发现改写出格式可以避免合成拉格朗日数据错误。

在controlDict中将writeFormat ascii; 改为 binary; 不会报错。。。。

在拉格朗日计算时，使用了org最新的负载平衡功能，但是不能合成数据，一合成就报错

--> FOAM FATAL IO ERROR: 
wrong token type - expected Scalar, found on line 34 the punctuation token '-'

file: /public/home/foamer/case/2P_H2+NH3_H25_M8_θ35°/ERL0.3_d75_sMF/processor1/1e-06/lagrangian/cloud/mass0 at line 34.

    From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)
    in file lnInclude/Scalar.C at line 102.

FOAM exiting

@李玉棠 blastFoam好像就可以

@东方白杨 我也没有很好解决，一年多没用SHM了，忘记了~

@EricLiu 你可能需要先选对应的分支版本，我估计你是直接下载的默认分支的，那个不对

rCF是这样的，库朗数一定要足够小，单纯加密网格没改善。这是我的对比计算结果

一些可能的分析（一家之言）：

注：bF为blastFoam，rCF为rhoCentralFoam，sMF是我自己的求解器

@欧阳

sudo apt install intel-mkl
dpkg -L libmkl-dev
gedit ~/.bashrc
export MKLROOT="/usr/include/mkl"
source ~/.bashrc
**exit**

source MKL，MKL 会污染系统 mpirun，导致 OpenFOAM 计算时并行挂起
oneapi  #这个命令会导致编译器变成 intel 编译器，所以 export MKLROOT 后要退出重进系统

./Allwmake --clean --platform MKL