同问，重心的位置怎么求，，，

@东岳 谢谢老师解答~

@东岳 感谢东岳老师的回答。

我参考了相关文献，low-Mach-number assumption 应该是忽略了压力波动，使得密度、焓等变量与压力波动无关。

那是不是意味着问题还是属于可压缩的呢？论文里面提到：

With these assumptions, the compressible energy equation is greatly simplified so that the resulting low-Machnumber equations are represented in tensorial notations as follows:

然后就是上面的三个控制方程。

那上面有关DNS的问题老师了解吗？

谢谢！

把文件里面的 simulationtype 设置成 laminar，网格足够密，时间步长足够小，是不是就算是DNS呢？

should be quasi-DNS

代码里面turbulence->divDevRhoReff(U)以及turbulence->alphaEff()这些湍流相关的项，会不会对DNS直接求解NS方程的方法造成影响？

如果设置为层流，turbulence->divDevRhoReff(U)直接返回层流项
因此不会有影响

@wwzhao 谢谢，sorry，没有仔细搜索，建议删帖

人找到了

@东岳 在 可压缩状态方程 中说：

可能是液体的压力和压强关系较大

我试过了，计算时这个方程可以考虑液体可压缩性并且结果也可以接受，本来我还以为有什么流体力学的理论给出的这个公式呢，谢谢东岳老师~

@东岳 谢谢东岳老师，确实好使。
只不过我在HPC上不能改标准求解器，还需要自定义一个myDPMFoam才能编译，运行的时候还需要定义div(alphaPhic,k.air)。研究好长时间才明白……

0_1541517377747_e8feb7f7-66e9-4596-8d6b-38fc6f7ad433-image.png
0_1541517320383_77c5f7ac-d93b-4e1a-b7ae-5cc1a483afae-image.png
将进口边界条件设置为大气边界层，运行报错，请大神指教哪里出了问题？

0_1541514111009_33d16365-121d-47bf-b842-32e36006ab8c-image.png
0_1541514148537_3c9f7071-5c62-4b8f-b151-e2d78c41ca7d-image.png 0_1541514164661_e93f6f1f-d619-4aec-87bf-a267cebe8fd3-image.png
0_1541514189389_f2cdec75-6e7b-4fe1-bd94-e6ea2a932467-image.png 0_1541514218249_5d745cd3-2ba3-4023-890d-71ad8358d24c-image.png
请问为什么我在knematicCloudProperties中设置的所有面都是颗粒反弹，但是实际运行过程中颗粒从下边界左边逸出了？

@tidedrinker 0_1541396115574_8a12b935-4737-407d-8a75-79e3bcf4fe98-image.png

老师，颗粒时间步长设置得越小，计算颗粒轨迹就越准确吗？有没有默认的MaxCo值？

Cellcells()不能并行，东岳老师是怎么做的？@东岳

@刘雄国 非常感谢

@j 太感谢了

@lambs 关于湍流的入口边界条件，可以参照《Large Eddy Simulation for Incompressible flow》 by Pierre Sagout 书中的两种方法：

1.Stochastic Recontruction from a statistical one-point description. 就是用随机方法生成湍流场，比较容易上手的有 Smirnov-Shi-Celik 方法，书中有详细介绍，要在流场中形成fully developed turbulent 需要很长的domain 和 time marching.

Deterministic precursor simulation. 在运行真正的算例之前，用循环边界条件加在一个domain 的四周，相当于一个无限长的 doamin， 生成充分发展的湍流场，然后将其中的一个断面map 到 真正算例里的入流边界，亦需要很长的时间达到quasi-steady state.

我觉得你的问题已经解决了，看起来是某个符号遗漏

modest radiative heat transfer
这本书太厚了，去看下周怀春老师的视频吧！

有关李老师后面的内容，硬交互，我不是很理解。“在OpenFOAM直接把一些数据提取出来，硬生成一个矩阵，进入F库计算后返回”，这个过程，是不是等同于让C程序计算完后，将结果数据带入F之中？

是的

我的F代码输入和输出都是个固定内容格式的txt文件。输入螺旋桨的几何特征型值表以及工况，最终拟定格式输出螺旋桨的定常或非定常性能。不是按照OpenFOAM的格式进行的输入输出

我想研究船桨干扰问题，船（C代码）和桨（F代码）两者是互相迭代的过程。这个迭代过程硬交互（Brute-force）能否实现自动，而不是完全孤立一般的手动迭代？

看起来是可以的，你用F代码计算，后生成txt文件，然后通过C语言把txt文件内的数据写成OpenFOAM的格式，然后用OpenFOAM计算船的流场。

这样如何

另外我觉得你这个不是大问题，如果做上述对接操作，纯前端的活，不需要太久。

软件版本：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 ***

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