受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)中间

@dzw05 好的，谢谢

你这一看就是几何体有缝隙..

@东岳 多谢！找到了！！

@dingcy PIMPLE确实可以进行内迭代，它是SIMPLE和PISO的耦合，增加内迭代的目的是为了可以采用更大的时间步长。我自己写过的程序，用LU-SGS时间推进，采用内迭代的目的是为了提高时间方向上的精度。
如你所说，不知道CFX采用了什么时间推进方法，所以也并不知道它内迭代的目的是什么。

变成32位重新编译解决了这个问题

@izumi 你可以贴一下你的fvSolution文件，因为不管使用哪种方法，只要收敛，结果大体上是应该差不多的

A是系数矩阵的对角部分，H是H操作符，是用于速度压力解耦的。

多谢大佬。大体了解了OpenFOAM类的调用结构。在声明指针变量的时候，往往是基类指针，运行的时候，通过读取我们fvScheme文件中的设置，这些指针指向不同子类中的函数进行计算

源项是一个什么样的表达式？是一个常数项还是与方程变量相关的源项？

各位大神，我有两个问题想向大家请教。如下图所示，为一固体热扩散问题，固体中具有热源。整体图形是一个圆柱，为了

简化，取了1/12作为对称单元进行模拟。六边形是一个个的组件，（整个固体）物性都是一样的。为了显示美观（或者说为

了把组件边界显示出来），这些组件边界都是在画网格的时候设定成interior，以便于在后处理软件中能够将其选择显示出来

该模拟是别人用COMSOL做的，我想用OpenFOAM来做一遍，遇到两个问题：

我发现用paraview不能显示这些组件的边界，paraview是一款很强大的后处理软件，我觉得应该是我没有找到显示的方法

六边形组件中的圆管流的是液体，用来冷却具有热源的固体，在该模拟中，只考虑固体热扩散。因此将圆管壁面设置成非

均匀的温度梯度边界条件即可。我在网上看了一些非均匀边界条件设置方法，都需要知道这些边界网格的位置信息。网上的

一些例子是简单的二维问题，而且是结构化网格，很容易得到边界网格的位置信息，然而该例边界处为非结构化网格，图形

也不规则，我不知道该如何提取边界处的位置信息，请求大家帮我解答，或者有其他不需要知道边界处网格的信息就能设置

非均匀边界条件的方法，求告知，不胜感激！！

0_1537691185627_a.PNG

我自己编写了一个小型数据库fin，里面有cfdTools finiteVolume fvMesh lnInclude surfaceMesh fields fvMatrices interpolation Make volMesh。其中wmake文件最后一句写成LIB = $(FOAM_USUER_LIBBIN)/libmyfiniteVolume，现在我想把库链接到icoFoam
中，请问他的options我该怎么修改呢
icoFoam的options文件
EXE_INC =
-I$(LIB_SRC)/finiteVolume/lnInclude
-I$(LIB_SRC)/meshTools/lnInclude

EXE_LIBS =
-lfiniteVolume
-lmeshTools

foamer新手，求指教如何往风场中加颗粒？

感谢东岳老师！我发现先读我里面已经写了的，是我自己没仔细！再次感谢