CFD中文网

《A Thorough Description Of How Wall Functions Are Implemented In OpenFOAM》一文2.6节描述了nutWallFunction，nutAtmRoughWallFunction中nutw公式与之相同，其他量有改变，我再看看。

VOF+动网格是interDyMFoam，用这个求解器就行

@ibelief 谢谢，我实际计算时，把管道也建模出来的。。

@东岳 东岳老师你好，我在fluent中选的是Euler多相流模型，RNG k-e湍流模型，
5.png
6.png

aurora，可以在word里用latex

这里整理了一些 找对象 的诀窍:xinxin:

谢谢老师~~

@东岳 谢谢东岳老师的关心，后来还是有点进展，不过进度不是很大，现在所走的路就是修改求解器，文章还没有出来。等出来之后我会详细说下这个问题。

@ibelief 刚刚试了一下，还是不可以，可能是那个地方的网格不太好，差异比较大。谢谢你的回答

是用八叉树算法吗？设置part上的最大尺寸也需要是2的幂吗？对这方面有了解吗？我也有这方面的一个疑问呢，之前在论坛发帖子还没人回。帮助文档里面也没写清楚，网上查了一些也没查到

@hhj0219 和y+定义一样。其实就是：流向网格间距除以壁面第一层距离，然后再乘以y+。

应该是Sb.boundaryField() 返回const reference
Sb.boundaryField().ref() 返回 reference
根据我的经验是这样,如果错误的话把编译错误贴上来
认真看看代码很简单的,编译器的错误已经告诉你所有信息了!

痛苦，结题更痛苦

[2] [1] ##00 [38] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) at ??:? [38] #1 Foam::sigFpe::sigHandler(int) at at ??:???:? [1] [2] ##11 Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHandler(int) at ??:? [38] #2 ? at ??:? [2] #2 ? at ??:? [1] #2 ? in "/lib64/libpthread.so.0" [38] #3 ? in "/lib64/libpthread.so.0" [2] #3 ? in "/lib64/libpthread.so.0" [1] #3 ? in "/lib64/libm.so.6" [38] #4 ? in "/lib64/libm.so.6" [2] #4 ? in "/lib64/libm.so.6" [1] #4 ? in "/lib64/libm.so.6" [38] #5 ? in "/lib64/libm.so.6" [2] #5 ? in "/lib64/libm.so.6" [1] #5 ? in "/lib64/libm.so.6" in "/lib64/libm.so.6" in "/lib64/libm.so.6" _pmiu_daemon(SIGCHLD): [NID 00126] [c0-0c1s15n2] [Tue Jul 9 10:43:27 2019] PE RANK 38 exit signal Floating point exception _pmiu_daemon(SIGCHLD): [NID 00082] [c0-0c1s4n2] [Tue Jul 9 10:43:27 2019] PE RANK 1 exit signal Floating point exception [NID 00126] 2019-07-09 10:43:27 Apid 8697814: initiated application termination

I have tried it again, the same problem occurs.
Last time i have used shear gradient boundary, and this time i just use zeroGraident. However, ....

I guess. the problem is the signal floating exception???

Tks

把你的新库通过controlDict挂载一下，在controlDict里面添加

libs ( "libnew1.so" "libnew2.so" );

这个是非常常见也常用的功能
不需要独立分区
一个推荐的方式
是循环cell 比如begin_c_loop（单核模式）下c_face_loop
在c_face_loop下采用C_FACE和C_FACE_THREAD来判别网格中的循环面是否面的thread id和你所匹配的面id一致
face的面积采用F_AREA获取 注意第一形参获取的是该面积的法向矢量
实际上按照流量的公式 一般还是点乘 毕竟是垂直面积方向的流量

好长，求高手神翻译..

微信图片_20190306174704.png

微信图片_20190306174712.png

....

最近计算一个方腔流动遇到一些问题，这个流动沿流向（x方向）是周期性边界，yz方向是壁面，流动固定x方向的平均速度，自动求解压力梯度，在采用rans计算时是正常的，可以计算出一个合理的结果

pressure gradient = 0.179097

但是当采用PANS（自己编译的湍流模型）计算的时候，压力梯度会变得非常小，导致计算结果很差，

pressure gradient = 0.0903267

自己编译的模型可以在channel395中正常计算，但是在该算例中出现压力的降低，不知道是什么原因引起的

psb.jpg

春风催碧叶
夜灯盼文章
横批：诸事如意

改了一下
春风催碧叶
夜灯盼收敛
横批：哈哈哈哈

碧叶双关
收敛双关