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

想模拟一个二维数值波浪水池里规则波与液箱的相互作用问题。
求解器用的是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的边界就好。

@东岳 哈哈，是OpenFOAM.

Ubuntu出错.PNG
2.PNG
打开Ubyntu后，无法进入任何文件夹，点击就会出现上面操作失败的提示。

哦，对我来说，可能就是编程上比较麻烦，好像我问过的人都说难，:zoule:

/OpenFOAM/OpenFOAM-2.4.0/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/0/T.water
这是我想修改的温度文件的位置。
初始蒜粒里面wall给的边界条件是zeroGradient.我想做一个壁面散热的效果（壁面以热通量W/m2的形式向外散热）很明显zeroGradient不再适合。
我想求教：这种热通量向外散热的边界条件该如何设置？是采用fixedGradient还是其他形式的壁面函数？
如果是fixedGradient，那么它的具体值的单位是什么？如果采用其他方式，那么有没有能使twoPhaseEulerFoam运行起来的成功的壁面函数？如果有，还望各路大神不吝赐教！！！
最后还是麻烦东岳老师帮我看一下，或者号召一下大家伙，看看我的问题。拜谢拜谢！
@东岳

solidTraction边界条件是HronTurekFsi算例中柔性固体的边界条件，tractionDisplacement是flappingConsoleSmall算例中的柔性固体边界条件，请问有人知道两者有什么区别吗？我给定同样的固体参数，分别用icoFisFoam和icoFsiElasticNonLinULSolidFoam算，但是结果差别很大。用第二个求解器求解的时候会产生特别大的变形，固体不会弹回来一直到网格变形太大算死。有人了解吗？求解答

https://github.com/OpenFOAM/OpenFOAM-3.0.x/blob/master/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C

去这里加个源项不就好了？

// Turbulent kinetic energy equation tmp<fvScalarMatrix> kEqn ( fvm::ddt(alpha, rho, k_) + fvm::div(alphaRhoPhi, k_) - fvm::laplacian(alpha*rho*DkEff(), k_) == alpha*rho*G - fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_) - fvm::Sp(alpha*rho*epsilon_/k_, k_) + kSource() // ? + K*k_ );

@东岳 嗯嗯 谢谢东岳老师！

@buaa_hgj 非常感谢

最近在看dynamicRefineFvMesh类，发现他和普通的八叉树网格细化有一些区别，有没有人知道这个算法最早是谁写的？有没有相关的参考文献？

@bestucan 时间太短了，烧不起来，只能强算，往后算

@史浩 ANSYS用的是单元生死技术，帝国理工好像是移动标架法。