软件版本：openfoam-v1806
硬件信息：超算（x86_64架构）
系统：centos 6.6
问题描述：openfoam编译完成后在登录结点串行并行计算都无报错，但是通过bsub命令提交任务计算出现下面的报错：

[m5605:28439] Signal: Segmentation fault (11) [m5605:28439] Signal code: Address not mapped (1) [m5605:28439] Failing at address: 0x29c38a0 [m5605:28442] [ 0] /lib64/libc.so.6(+0x326a0)[0x2b582a1406a0] [m5605:28442] [ 1] /home/daiyut/tools/openmpi215/lib/openmpi/mca_btl_openib.so(+0xf380)[0x2b583afbf380] [m5605:28442] [ 2] /home/daiyut/tools/openmpi215/lib/libopen-pal.so.20(mca_btl_base_select+0x156)[0x2b582b7370f6] [m5605:28442] [ 3] /home/daiyut/tools/openmpi215/lib/openmpi/mca_bml_r2.so(mca_bml_r2_component_init+0x11)[0x2b583adaf251] [m5605:28442] [ 4] /home/daiyut/tools/openmpi215/lib/libmpi.so.20(mca_bml_base_init+0x7c)[0x2b582aefc0ec] [m5605:28442] [ 5] /home/daiyut/tools/openmpi215/lib/libmpi.so.20(ompi_mpi_init+0x3ad)[0x2b582aebf59d] [m5605:28442] [ 6] /home/daiyut/tools/openmpi215/lib/libmpi.so.20(MPI_Init_thread+0x46)[0x2b582aede3f6] [m5605:28442] [ 7] /home/daiyut/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so(_ZN4Foam8UPstream4initERiRPPcb+0x67)[0x2b582a4a90c7] [m5605:28442] [ 8] /home/daiyut/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7argListC1ERiRPPcbbb+0x67c)[0x2b5828e2489c] [m5605:28442] [ 9] icoFoam[0x453149] [m5605:28442] [10] /lib64/libc.so.6(__libc_start_main+0xfd)[0x2b582a12cd5d] [m5605:28442] [11] icoFoam[0x41c6c9] [m5605:28442] *** End of error message ***

请问有前辈知道这个是什么原因造成的吗？以及有什么解决办法吗？

好像没有。

给你举个例子，在计算压力p和速度U的时候，并行分解之后会出现一个processor0to1类似这种边界条件。不处理的话，phi在processor0to1是通过U计算出来的。你可以手动处理processor0to1这个边界条件的phi值，比如全部设为0. 然后你的处理器边界就0通量了

@wwzhao 你是指严格小于1吗？不一定的，我们也用过超过2的库朗数，结果也还行。

@东岳 老师，我重新把coneNozzleInjectionModel文件全部复制到一个新的文件夹下，然后重命名为MyInjection然后将模型文件MyInjection.C与MyInjection.H中的模型名称coneNozzleInjection全部替换成MyInjection，编译完成后运行sprayFoam发现关键字myInjection并不能被识别为可用的颗粒喷射模型。下图为我运行求解器的结果：
0_1540376406398_2018-10-24 18-17-46 的屏幕截图.png
其中颗粒喷射子模型修改我是参照下面地址所示的步骤进行的：https://openfoamwiki.net/index.php/Howto_addSpraySubModels
0_1540376624909_2018-10-24 18-22-03 的屏幕截图.png
请问老师，这是不是哪里文件没有修改造成的？

谢谢，我尝试下重新编译IHFOAM看能不能解决

如果是稳态，不需要增加fvm::ddt()，但看起来问题不是来源于fvm::ddt()。更详细的错误信息判断不出了。

@东岳 好的，多谢东岳前辈的意见^_^

想讲提取的自由液面波动的曲线数据输出，利用plot date工具，其中只涉及自由表面波动的问题能得到理想的曲线，如下面图1和2，可是当自由液面出现气泡时，曲线变得波动起伏，不能按理想曲线排列，如下面图3和4，想知道该用哪个filter怎样改变设置条件能得到云图中那样理想的曲线数据，希望懂的大神能指点一二，谢谢！
0_1540196637813_5abc6bca-edf7-48d3-98bb-eea182ac6014-image.png
0_1540196684260_cf5d2c26-19e0-4907-93f2-df12199272d6-image.png
0_1540196828240_58c4864b-8945-4d63-8db6-1765ae3da3ef-image.png
0_1540196846936_b07d6dd2-daf2-462d-84ae-d46ee372bfcf-image.png

为社区做贡献😄

@东岳 哦了，谢啦！

0_1539859784313_69ca8214-5828-4a5a-aea6-836bf99a13c1-image.png 上面的问题已解决，出现新的报错，是因为定义面/点不同吗？

没人有类似的需求？

@cc_cherry 你好，请问一下你知道a1,a2,a3,a4怎么选取的吗？谢谢您，麻烦您了。

大家好，本人最近在模拟一个具有固体区域和液体区域的模型，需要考虑固液之间的换热，所以想请问是否有人用过chtMultiRegionFoam求解器，以及能否提供一下算例学习一下，谢谢

嗯嗯 是使用wavetransmissive边界条件 谢谢

@东岳 老师，是用filedAverage这个函数来获取LES流场的平均p、U、k吗？但是用LES计算时要求时间步很小，计算到RANS稳态结果的时间计算量非常大

这个问题我也遇到过，是不是用了动网格求解器进行了并行计算？

@东岳 谢谢东岳老师，我还想问一下，OpenFOAM自带算例中的topoSet一般都是指定一个方形区域（boxToCell），如果我想指定一个圆柱区域该如何使用呢？关于topoSet，老师有没有什么较为系统和详细的资料推荐呀？

在计算粒子的Eo的时候，OpenFOAM定义

template<class ParcelType> inline Foam::scalar Foam::KinematicParcel<ParcelType>::Eo ( const vector& a, const scalar d, const scalar sigma ) const { vector dir = U_/(mag(U_) + ROOTVSMALL); return mag(a & dir)*(rho_ - rhoc_)*sqr(d)/(sigma + ROOTVSMALL); }

注意上面的(rho_ - rhoc_)，如果粒子密度小于连续相密度，在计算升力的时候：

scalar Eo = p.Eo(g, p.d(), sigma_); scalar dH = p.d()*cbrt(1.0 + 0.163*pow(Eo, 0.757));

对于Eo的负值，pow(Eo, 0.757)为虚数

最近在用interFoam算板上液膜流动，用固定步长计算时，液膜表面观察不到任何变化，用变步长计算到同样的时间时，又可以观察到液膜表面的波动，请问有人遇到过吗？
固定步长计算结果：
0_1539420931486_ae58dccc-52ef-4d14-be19-14a6aeb0090e-图片.png
变步长计算结果：
0_1539421027533_e4c454a3-88ec-4842-a05b-480f07a0f73c-图片.png