@东岳 谢谢东岳老师，我还想问一下，OpenFOAM自带算例中的topoSet一般都是指定一个方形区域（boxToCell），如果我想指定一个圆柱区域该如何使用呢？关于topoSet，老师有没有什么较为系统和详细的资料推荐呀？

在计算粒子的Eo的时候，OpenFOAM定义

template<class ParcelType> inline Foam::scalar Foam::KinematicParcel<ParcelType>::Eo ( const vector& a, const scalar d, const scalar sigma ) const { vector dir = U_/(mag(U_) + ROOTVSMALL); return mag(a & dir)*(rho_ - rhoc_)*sqr(d)/(sigma + ROOTVSMALL); }

注意上面的(rho_ - rhoc_)，如果粒子密度小于连续相密度，在计算升力的时候：

scalar Eo = p.Eo(g, p.d(), sigma_); scalar dH = p.d()*cbrt(1.0 + 0.163*pow(Eo, 0.757));

对于Eo的负值，pow(Eo, 0.757)为虚数

最近在用interFoam算板上液膜流动，用固定步长计算时，液膜表面观察不到任何变化，用变步长计算到同样的时间时，又可以观察到液膜表面的波动，请问有人遇到过吗？
固定步长计算结果：
0_1539420931486_ae58dccc-52ef-4d14-be19-14a6aeb0090e-图片.png
变步长计算结果：
0_1539421027533_e4c454a3-88ec-4842-a05b-480f07a0f73c-图片.png

解决了，我太粗心了，谢谢老师！！！！

@东岳
另一个alphatJayatillekeWallFunction和Fluent中的壁面函数一样，根据$y^+$使用不同的公式计算，打算试一下这个

@maoyanjun_dut 好的，谢谢小毛哥！:quwan:

不知道你测试过没有，我之前的回复是可行的，所有

对‘yyFlexLexer::yywrap()’未定义的引用

的错误提示，都可以这样处理。近期做了个OpenFOAM-2.2.x的现成的Ubuntu虚拟机，可在这里直接下载

@东岳 嗯，谢谢您。
这样的确有点复杂，不过因为我要模拟街道峡谷污染物扩散，如果用循环的话污染物一直在计算域内出不去，然后结果误差较大，所以是不得已而为之。

没用过tecplot，有什么功能paraview不能实现么？

@东岳 thanks

@wwzhao 太感谢了，新手上路，昨天遇到好多首字母大小写报错的问题，最后看得眼花了没看到这个标点错误
感谢指正

这个是我自用的，不知道tutorials里面有没有现成的：

functions { minMax { // Type of functionObject type fieldMinMax; // Where to load it from (if not already in solver) libs ("libfieldFunctionObjects.so"); // Function object enabled flag enabled true; // Log to output (default: false) log true; // Write information to file (default: true) write false; location false; // Fields to be monitored - runTime modifiable fields ( particleU0 //particleU1 //particleU2 //alphaSequence0 //alphaSequence1 //alphaSequence2 ); } }

@piteqiu 我觉得没有

受OpenFOAM中wall function实现方式的影响，贴近壁面处的du/dy并不是真实的速度梯度，因此用du/dy求出来的并不是真实的wall shear stress。

也就是说du/dy的算法只适用于laminar flow，而不适用于turbulent flow。

@regina #3 Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) at ??:?应该是GAMG代数方程求解器的DIC预处理器遇到值为0的主对角元素了，除0错误。需要检查一下算例。

@random_ran 万分感谢:chouchou: :chitang: :chitang:

看错了看错了，符号在pEqn.H里面处理了

@东岳 谢谢东岳老师

http://dyfluid.com/docs/book/_book/chapter2/section2.3.html

我再后续要更新你这些东西，要是看这个资料看不明白，欢迎反馈我来更新

那个-L后面不需要<>号

@东岳

在rhoPimpleFoam/LES/pitzDaily 算例里找到了这个应用，

---Properties Second order Unbounded Attempts to limit the unboundedness of the linear scheme ---Usage The scheme is specified using: divSchemes { default none; div(phi,U) Gauss filteredLinear2 <k> <l>; } ---Where the coefficients --k: Scales the rate at which the correction is applied 0 : linear 1 : fully limited --l: Maximum allowed overshoot/undershoot relative to the difference across the face. 0 : no overshoot/undershoot 1 : overshoot/undershoot equal to the difference across the face

对于矢量时候 用filteredLinear2V， 求解出来流场的unboundedness 介于 Gauss Linear (unbounded) 和 Gauss GammaV (bounded)中间