入门不久，尝试着答一下。
我认为openFoam中的边界条件大致分两类：

强制性边界条件，直接规定数值，如fixedValue，fixedMean，这类边界条件的信息会一直存在于求解所得的结果文件，算很多步以后查看文件就会发现边界条件并没有变化； 参考性边界条件，这类边界条件大多是第二类边界条件，如出口定义inletOutlet， 对于流出的气体而言，该边界其实是zeroGradient，这时定义的value只是初始迭代的时候给的一个初值，该边界的数值由流场内部的点插值得到

1/3次幂咯

@bestucan 好的，谢谢你的指点~:ok2:

@东岳 我的边界是没有网格的，但是默认的显示就是在紧挨边界的网格上，就是你说的内网格呗，但是同样的算例，Tutorials里就能把那个数据给放到壁面上，应该是经过处理了吧。

@OF初学者周 最后没有用gdb

感谢分享 :xiexie:

@hurricane007 这个思路挺不不错。我没见过手动分区的例子，不懂了 :zoule: 加油。

最近想做一个简化的模拟，需要一种边界条件，在这个边界上，某个组分的流量等于min(最大可能的流量，一个常数）。考虑准稳态情况下，菲克扩散定律 J =-D gradient (C)*n. 在使得边界网格内这个组分不为零的情况下，从这个网格向边界扩散的流量应该小于这个网格接受其他网格扩散的量。
那么问题来了…… 这个边界上的网格，除了边界面，其他面上的扩散流怎么求……求助一下大神们……

@xpqiu 想请教一下，如果我想找使用了某个模型的算例该怎么找。例如我现在想找个使用了LangmuirHinshelwoodReactionRate 模型的算例能不能找到呢，这个模型在：
/src/thermophysicalModels/specie/reaction/reactionRate/LangmuirHinshelwood/LangmuirHinshelwoodReactionRate.H

可以算一下compressibleInterFoam下的传热:duang:

@东岳
点赞了，自己写出来的第一个OpenFOAM代码。。。

@东岳
谢谢！这下明白了！

想模拟一个二维数值波浪水池里规则波与液箱的相互作用问题。
求解器用的是waveDyMFoam ,属于interFoam里的二相流求解器。
不知道如何设置边界条件。
刚开始把水箱的外表面设为wall1,内表面为wall2，两者的边界条件如下：
P：fixedFluxPressure;
V：movingWallVelocity;
alpahwater:zeroGradient
pointdDisplacement: calculated;
Dynamicdic设置没有问题。
结果跑不起来。
我又把内边界wall2的边界条件改为：
V：
PointDisplacement：fixedValue

V：fixedValue
其余不变，结果还是跑不起来。
请问有大神知道哪里出了问题吗？感激不尽:xinxin:
!cfd2.png cfd.png

目前在收集经典的欧拉欧拉气液流动的实验数据，打算通过OpenFOAM对比一下。手头有几个，比如
Díaz, M. E., Iranzo, A., Cuadra, D., Barbero, R., Montes, F. J., & Galan, M. A. (2008). Numerical simulation of the gas–liquid flow in a laboratory scale bubble column: influence of bubble size distribution and non-drag forces. Chemical Engineering Journal, 139(2), 363-379.
打算收集20个，有人有合适的可以推荐么？
我通过OpenFOAM模拟并和实验数据对比之后，会把所有case公开供后来者研究。

@东岳 ,我有点想当然了，我看到很多算例在设置kqWallFunction时，会设置value，这个value的值一般采取与inlet的值一样，我认为既然我的入口是time-varying的，那么这个值也应该是time-varying的。paper.png
https://www.cfd-online.com/Forums/openfoam-pre-processing/70893-standard-kepsilon-kqrwallfunction-epsilonwallfunction-why-should-give-value.html
提到这个value仅仅是初始时刻的猜测值，在之后的时间步里面调用k,epsilon的WallFunction时，WallFunction表现得像zeroGradient。
我的算例初始时刻速度为0，设置k和epsilon很小应该没问题

@王大力 在 fvc::reconstruct()操作 中说：

fvc::reconstruct()的想法是通过对面心值加权平均，得到体心值。

要不然呢？类似这种操作都有加权。假如，你的四边形，如果三个面非常小，另一个面非常大，大面肯定有更大的权重影响体心值

U是体积平均速度

U = U_l\alpha_l + U_v\alpha_v

其中U_l和U_v都是filtered单相流体的速度

@东岳 好的，谢谢东岳老师~~~

以前没学过c+，这几天看了看琢磨一下，
最后我是这样设置的

INLET { type codedFixedValue; value uniform (0 0 0); name IN; code #{ const vectorField& Cf = patch().Cf(); vectorField& field = *this; forAll(Cf,faceI) { if( Cf[faceI].z()<300) { field[faceI] = vector(20*pow((Cf[faceI].z()/10), 0.15), 0, 0); }else{ field[faceI] = vector(30, 0, 0); } } #};

@浪迹天大 OpenFOAM中是先把模型分块（decomposePar），然后在进行并行计算。在某些动网格算例中，这些并行的块有的会增加网格，即网格数量变多，有的就保持不变，因此容易导致某些块中的网格量变得巨大，影响计算效率。我想是怎么能将这些网格重新分配到每一个并行进程中。
关于layerAddRemoval出错的问题，是OpenFOAM中并行的bug，这个暂时没有办法保证，避免的办法就是让移动的网格不要扫过processor的边界就好。