CFD中文网

@qingdong_wang Info << phi即可

86f6fd36-007a-42a5-b105-de5d4fcefffa-image.png
cited:http://openfoamwiki.net/index.php/See_the_MRF_development

图为在旋转坐标系内解绝对速度的控制方程，其中 \vec \omega \times \vec U_i 这项称为科氏力，旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程，其中的通量是相对通量（\vec U_R \vec U_I）,同时由于动量是一个矢量，因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献，和由于运输的物理量U_I的扩散贡献。

因此，我在思考是不是诸如湍动能，湍流耗散率这些标量，也可以认为在旋转坐标系内，运输物理量为湍动能，湍流耗散率这些标量，通量为相对速度的通量，但是由于这些量都是标量，因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小，此处的voricity是不是也应该用相对速度？

如果你认为动量方程增量很小，可以尝试把动量方程的松弛因子设为1，看看稳定性:duang: 好像问题不是这么简单

我找着我第一次发的方程改得，就是第二个block的代码，算出来的结果和MPPICInterfoam的基本一样。当然有可能我的粒子直径比较小，对计算结果影响不大？

我没明白您下面写的这个是什么意思。如果不看mrf的话， 你这个和interfoam相比就是减上了一项

fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

这一项就是65293cb6-a527-49eb-bdc4-2355ec82c507-image.png
我不知道怎么推出的这一项？

pEqn里面用到了rAUc

就是全场吧， 因为1-alpha1 是对于水为0， 所以

gsum（（1-alpha1） × mesh.V()）

就是bubble的volume；下面的那个速度也是一样的， 但是那个速度只是一个平均速度。

Double check if the lagrangian particle is two-way coupled with the fluid phase. If it is one-way coupled, the restriction can be abandoned.

positionfiles可以用matlab或者python处理

@东岳 但是那个貌似是在v1906这种版本，我现在用的是OpenFoam 7, 貌似没有MPPICInterFoam

WechatIMG47.png

请问一下, 如果得到这个bubble的projection area和horizontal diameter呢?
我现在尝试两种方法, 第一种对气泡旁的alpha数据在python中进行3D griddata处理, 得到alpha < 0.5的数值, 然后将所有坐标投影到一个平面(比如 z=0), 然后使用CV2进行图像分析,得到半径, 圆心和投影面积
但是因为需要计算量巨大, 我们在griddata中的space number取30-50, 这样造成误差比较大.

第二种办法是使用openfoam的后处理sample获取 z=z_bubble的constant plane然后进行2D griddata, 得到的数据比较好, 但是问题是当遇到图2中的气泡时, 这个方法就不好用了, 请问openfoam有没有一种后处理方法可以将气泡表面捕获然后给出projection的坐标呢?

好像foam_extend里面已经有了Rusche这个代码？你看看？
这个算法我看过但没自己植入过。因此这些细节问题也不好说，没什么建议了

@李东岳 哦哦，我想问问老师，那对于追踪大气泡（尺度和网格单元差不多）的LPT方法有相关的论文吗？或者相关的方法？

patch和wall基本一样，只不过wall可以用壁面函数，对结果没影响

[2] [1] ##00 [38] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) at ??:? [38] #1 Foam::sigFpe::sigHandler(int) at at ??:???:? [1] [2] ##11 Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHandler(int) at ??:? [38] #2 ? at ??:? [2] #2 ? at ??:? [1] #2 ? in "/lib64/libpthread.so.0" [38] #3 ? in "/lib64/libpthread.so.0" [2] #3 ? in "/lib64/libpthread.so.0" [1] #3 ? in "/lib64/libm.so.6" [38] #4 ? in "/lib64/libm.so.6" [2] #4 ? in "/lib64/libm.so.6" [1] #4 ? in "/lib64/libm.so.6" [38] #5 ? in "/lib64/libm.so.6" [2] #5 ? in "/lib64/libm.so.6" [1] #5 ? in "/lib64/libm.so.6" in "/lib64/libm.so.6" in "/lib64/libm.so.6" _pmiu_daemon(SIGCHLD): [NID 00126] [c0-0c1s15n2] [Tue Jul 9 10:43:27 2019] PE RANK 38 exit signal Floating point exception _pmiu_daemon(SIGCHLD): [NID 00082] [c0-0c1s4n2] [Tue Jul 9 10:43:27 2019] PE RANK 1 exit signal Floating point exception [NID 00126] 2019-07-09 10:43:27 Apid 8697814: initiated application termination

I have tried it again, the same problem occurs.
Last time i have used shear gradient boundary, and this time i just use zeroGraident. However, ....

I guess. the problem is the signal floating exception???

Tks

带摩擦的运输方程是 u_t = lambda * u_x = eps * u_xx, 通过变化, 我们得到了, 0_1545168405877_3bd01f52-6748-4ece-bac9-e55eeef96e15-image.png ,其中D是导数矩阵, D2 是二阶导数矩阵
然后, 我们对M进行求特征值(EW), 得到了这种图片, 我想知道这是为什么呢?
0_1545168538892_dba5eac1-7474-4ce6-b089-77d3d8903c03-image.png