@东岳 非常感谢前辈，已经解决。

@xiaohsu2006 请问问题解决了吗？我也遇到了类似的问题。谢谢！

解决了，边界条件设置的问题，

inlet { type uniformFixedValue; uniformValue table ( (0 (1 0 0)) (1000 (2 0 0)) (2000 (3 0 0)) ); value $internalField; }

少设置了 value $internalField;
结果没有生效

有可能是分配的processor过多

我今天在decomposePar时出现了FOAM Warning 也是提示no points in mesh 和 no cells in mesh

我把system下decomposeParDict里面的numberOfSubdomains减小了 同时相应地给x、y、z方向分配的数量也减小了，再运行decomposePar就不会报错了。
:chitang:

谢谢。

谢谢李老师

@星星星星晴 然而也没其他办法了吧:zoule:

@东岳 好的，谢谢东岳老师的回复

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的坐标呢?

好的 谢谢东岳老师！

@东岳 东岳老师高抬了- -我做OpenFOAM才做了不到两年，你才是大佬啊

很古老的OpenFOAM里面存在pd，我已经记不清什么概念了，和p_rgh不是一个东西

好的，我认真看看，非常感谢您的回答，谢谢

谢谢东岳老师。
我的研究目的是看使用毛巾往复擦桌子过程中，由于毛巾的往复运动，引起的周围气流变化，毛巾速度是0.2m/s，如下图，先向左前进1s，再向右运动1s，回来。d49e6d5e-44a9-45c2-add6-7afe01a3bee3-image.png
这种匀速的往复匀速，速度边界条件应该如何给定？

谢谢您。

@zhe 直接设置水体区域，没有另外处理

dyfluid@dyfluid:~/twoWayHyQMOMFoam/tutorials/PTC_2D/testCase2$ oneWayHyQMOMFoam oneWayHyQMOMFoam: symbol lookup error: oneWayHyQMOMFoam: undefined symbol: _ZN4Foam6HyQMOMC1ERKNS_7PtrListINS1_INS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEEEE

原因：当时调整了src下的的构造函数，紧跟着solver那面的构造函数也加了一个参数进去。然后重新编译了src，忘记了编译solver，出错。
方案：重新编译solver后解决。

#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigSegv::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"

原因：在class里面增加了一个私有成员，编译库之后，求解器可以对类进行初始化，随即报错。报错只有3行，没有任何其他信息
方案：编译class之后，尤其这种增加私有成员的时候，记得重新编译solver一下就可以解决

@东岳 ok:huahua:

嗯嗯，好的，谢谢李老师指点:chitang:

边界条件、网格、格式

不可压缩LES模型算法非常成熟了，一般都是设置不正确带来的发散

按照twoPhaseEulerFoam中的EEqns.H写出能量方程如下（算例文件中能量形式是sensibleInternalEnergy）大神们先看看这么写对吗？

fvScalarMatrix E1Eqn ( fvm::ddt(alpha1, rho1, he1) + fvm::div(alphaRhoPhi1, he1) - fvm::Sp(contErr1, he1) + fvc::ddt(alpha1, rho1, K1) + fvc::div(alphaRhoPhi1, K1) - contErr1*K1 + ( he1.name() == thermo1.phasePropertyName("e") ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p) + p*fvc::ddt(alpha1) : -alpha1*dpdt ) - fvm::laplacian ( fvc::interpolate(alpha1) *fvc::interpolate(thermo1.alphaEff(phase1.turbulence().mut())), he1 ) ); E1Eqn.relax(); E1Eqn -= ( Kh*(thermo2.T() - thermo1.T()) + Kh*he1/Cpv1 - fvm::Sp(Kh/Cpv1, he1) + alpha1*rho1*(U1&g) + fvOptions(alpha1, rho1, he1) );

a8753d25-8ea0-4cf2-8fae-10a71aa5cdfd-image.png
那如果不考虑传热的话，只把导热系数设为0和温度设为一样，这样做对吗？
我为了不考虑传热，是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了，我这种做法对吗？