inletOutlet边界条件针对流入和流出的情况需要作区分，比如流出的时候是零法向梯度，流入的时候需要指定某个值。判断流入和流出的信息是通过phi来识别的。OpenFOAM中边界处的phi在大于0的时候，全部流出。因此对于这一类条件，需要指定phi，你的案例中，指定的是phi.air

@东岳 i其实是第一个网格，n是最后一个网格
也就是
先求每一个网格中的alpha * x 再对整个计算域求和
然后求和所有计算域中alpha的值
两者做比值

不好意思 数据在这里
0_1542843045909_48db69f6-3dde-40bc-9b8d-70c1b51caa28-image.png Dd [um] Number
0.00000220 17200.00000
0.00000439 71000.00000
0.00000659 97900.00000
0.00000878 84400.00000
0.00001100 53900.00000
0.00001320 35000.00000
0.00001540 23900.00000
0.00001760 15300.00000
0.00001980 12200.00000
0.00002200 9060.00000
0.00002410 6550.00000
0.00002630 5540.00000
0.00002850 4090.00000
0.00003070 3440.00000
0.00003290 2530.00000
0.00003510 2210.00000
0.00003730 1960.00000
0.00003950 1490.00000
0.00004170 1240.00000
0.00004390 938.00000
0.00004610 860.00000
0.00004830 743.00000
0.00005050 652.00000
0.00005270 550.00000
0.00005490 427.00000
0.00005710 339.00000
0.00005930 282.00000
0.00006150 268.00000
0.00006370 194.00000
0.00006590 208.00000
0.00006810 171.00000
0.00007030 131.00000
0.00007240 106.00000
0.00007460 125.00000
0.00007680 76.00000
0.00007900 66.10000
0.00008120 52.80000
0.00008340 63.80000
0.00008560 23.50000
0.00008780 42.40000
0.00009000 34.10000
0.00009220 33.50000
0.00009440 18.30000
0.00009660 21.60000
0.00009880 17.70000
0.00010100 14.00000
0.00010300 10.40000
0.00010500 10.30000
0.00010800 10.20000
0.00011000 6.71000
0.00011200 0.00000
0.00011400 3.29000
0.00011600 6.51000
0.00011900 0.00000
0.00012100 0.00000
0.00012300 0.00000
0.00012500 0.00000
0.00012700 0.00000
0.00013000 0.00000
0.00013200 0.00000
0.00013400 0.00000
0.00013600 2.98000
0.00013800 2.96000
0.00014100 0.00000
0.00014300 0.00000
0.00014500 0.00000
0.00014700 0.00000
0.00014900 0.00000
0.00015100 0.00000
0.00015400 0.00000
0.00015600 0.00000
0.00015800 0.00000
0.00016000 0.00000
0.00016200 0.00000
0.00016500 0.00000
0.00016700 0.00000
0.00016900 0.00000
0.00017100 0.00000
0.00017300 0.00000
0.00017600 0.00000
0.00017800 0.00000
0.00018000 0.00000
0.00018200 0.00000
0.00018400 0.00000
0.00018700 0.00000
0.00018900 0.00000
0.00019100 0.00000
0.00019300 0.00000
0.00019500 0.00000
0.00019800 0.00000
0.00020000 0.00000
0.00020200 0.00000
0.00020400 0.00000
0.00020600 0.00000
0.00020900 0.00000
0.00021100 0.00000
0.00021300 0.00000
0.00021500 0.00000
0.00021700 0.00000
0.00022000 0.00000
0.00022200 0.00000
0.00022400 0.00000
0.00022600 0.00000
0.00022800 0.00000
0.00023100 0.00000
0.00023300 0.00000
0.00023500 0.00000
0.00023700 0.00000
0.00023900 0.00000
0.00024100 0.00000
0.00024400 0.00000
0.00024600 0.00000
0.00024800 0.00000
0.00025000 0.00000
0.00025200 0.00000
0.00025500 0.00000
0.00025700 0.00000
0.00025900 0.00000
0.00026100 0.00000
0.00026300 0.00000
0.00026600 0.00000
0.00026800 0.00000
0.00027000 0.00000
0.00027200 0.00000
0.00027400 0.00000
0.00027700 0.00000
0.00027900 0.00000
0.00028100 0.00000
0.00028300 1.84000
0.00028500 0.00000
0.00028800 0.00000
0.00029000 0.00000
0.00029200 0.00000
0.00029400 0.00000
0.00029600 0.00000
0.00029900 0.00000
0.00030100 0.00000
0.00030300 0.00000
0.00030500 0.00000
0.00030700 0.00000
0.00031000 0.00000
0.00031200 0.00000
0.00031400 0.00000
0.00031600 0.00000
0.00031800 0.00000
0.00032100 0.00000
0.00032300 0.00000
0.00032500 0.00000
0.00032700 0.00000
0.00032900 0.00000
0.00033200 0.00000
0.00033400 0.00000
0.00033600 0.00000
0.00033800 0.00000
0.00034000 0.00000
0.00034200 0.00000

从代码上解释更好说：

scalarField U(-phi/patch().magSf()); //其中的phi你要给定（指定通量），然后即可计算出速度U的大小（模） vectorField n(patch().nf()); //面的法向 if (phiName.dimensions() == dimVelocity*dimArea) { // volumetric flow-rate operator==(n*U); //速度 } else if (phiName.dimensions() == dimDensity*dimVelocity*dimArea) { const fvPatchField<scalar>& rhop = patch().lookupPatchField<volScalarField, scalar>(rhoName_); // mass flow-rate operator==(n*U/rhop); //速度

你看下/combustion/fireFoam/LES/flameSpreadWaterSuppressionPanel的速度边界条件设置。

还没见过教材里面讨论low-Mach-number assumption，气动声学这块我还没研究过 :zoule:

Source[kth].ref() += (1 - kth)/(kth + 1)* pow(a[n], kth + 2)*g*w[n];

上面要改为

Source[kth].ref() += (1.0 - kth)/(kth + 1.0)* pow(a[n], kth + 2)*g*w[n];

否则 (1 - kth)/(kth + 1)在kth大于1的时候为0

东岳大哥您好，算几十万的网格是没有问题的，我扩展了虚拟机的内存到64g和虚拟内存到10g都不行。 。。也是奇怪

请问在OpenFOAM中，如何定义一个轴对称模型来仿真Taylor Couette flow？如何设定wedge边界？谢谢大神！

@mr-h 在 标量场壁面梯度的处理方式 中说：

在温度方程里的拉普拉斯项fvm::laplacian(ka,T) 中 存在ka的梯度的计算。

此处Ka为volScalarField，那么fvm::laplacian(ka,T)是不是等价于(fvm::grad(ka)&fvm::grad(T)) + ka*fvm::laplacian(T);这两种表达方式的意义一样吗？

@东岳 是类似与MultiRegionFoam的方法吗？存在多套网格分别对应不同的region？这种方法我没试过。我用的求解器是interfoam，现在只是想的是从一套网格上生成一个subMesh，这样绝大部分的网格节点是对应的，求解标量C方程所需要的数据可以通过原始的网格map过来。

图片
我最近在使用reactingParcelFilmFoam计算水滴与结构表面撞击并形成水膜的过程中，我想往模型中添加一个水滴的收集系数，这个系数需要知道水滴在撞击到结构时的位置坐标和初始的位置坐标，请问有人知道该怎么确定这两个坐标吗？

从代码来看，e表示回弹系数，mu表示一种损耗。e越大，回弹越完全，mu越大，速度被消耗的越大。这些值影响你颗粒撞击壁面的回弹行为。
你可以手动计算一下，看看U有什么变化

case PatchInteractionModel<CloudType>::itRebound: { keepParticle = true; active = true; vector nw; vector Up; this->owner().patchData(p, pp, nw, Up); // Calculate motion relative to patch velocity U -= Up; scalar Un = U & nw; vector Ut = U - Un*nw; if (Un > 0) { U -= (1.0 + e_)*Un*nw; } U -= mu_*Ut; // Return velocity to global space U += Up; break; }

AFAIK, you are the only one working in this area. Any update?

@东岳 好的，明白啦，我试一下。谢谢东岳老师。

@东岳 老师，您好。我目前在CFD中文网只发了一封求助帖。但因帖子中涉及到我们课题组提出的在某方面研究的新结构模型，目前还未申请专利和发表论文。为防止我们的结构创意被别人抢先用去，麻烦您把我目前唯一的这一封帖子给删了。发帖时是我有欠考虑，由此给您带来的麻烦，我十分抱歉！！谢谢您，麻烦了！！！

已删除

我在用interFoam时，笛卡尔网格可以得到合理的计算结果。但换成非结构网格（prism之类）后，计算虽没有发散，但界面有问题，请问有什么需要改进的地方吗？@东岳 @tidedrinker

无动网格，parcel量200000，k-epsilon，sprayFoam，spraycloud，solver和边界未进行任何其他编译

@东岳 好滴:xinlei: :xiezuoye:

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

@东岳 谢谢老师解答~