抱歉，找到原因了，是我忘记wclean了

RT，我在做管流的模拟，现在想在一个横截面上做出雷诺剪切力（u'v' VS r/D）的曲线图象。现在用openfoam的sample功能提取了大量的速度-位置数据（100个时间步长, 每个时间步长在该横截面上有360个径向直线，在一共36000组，每组直线取了200个点的速度（u，v，w）的数据，所以现在在MATLAB上得到了一个200×36000的矩阵。现在想根据这个矩阵得到雷诺剪切力的图像，请问怎么操作呀？

这个是我画了一下午的wedge网格，先给所需，在这个基础上改吧........

@东岳 谢谢大神回答！

https://www.ubuntu.com/download/desktop

琢磨要不要换一个

@浪迹天大 非常感谢您的帮助，我已经完美解决这个问题。

http://www.cfd-china.com/topic/1548/

好像跟这个帖子有关。气动噪声算法上我还不太理解。目前只知道是计算压力小波的传递？

你用的是非结构网格？很有意思 哈哈

在你画出来流线图之后，在你原始的.foam基础上，再做contour plot，然后同时显示流线和contour。

你试试运行这个命令：

sudo chmod u+x /home/lnq/OpenFOAM/

有时候相间的物质转换会造成局部高压，然后碰到wall就会反弹。。。
一般non reflective boundary都是用于出口入口的，那么它
可以用在wall上吗。

他是不是想截这个图 :mihu:

替代文字

而且用的是压力基 两相流

谢谢！经过测试，完美地解决了我的问题！

谢谢大神，我看到你原来的代码里有相关的。再次感谢。
http://www.cfd-china.com/topic/1228/owner-neighbour-网格寻址/2

在外挂硬盘上计算非常慢，不写log可以，写log的话，在机械硬盘上写log严重耽误了计算的时间。

另外，log文件有时候变的特别大，我这个算例时间步长比较小，写文件间隔比较大，导致结果没有多大，但是log上G。可以通过 > /dev/null来禁止写log，当然如果你确定不需要看log的话。

0_1525840011230_2018-05-09 12-25-51屏幕截图.png

约2G左右

https://www.sciencedirect.com/science/article/pii/S0010465517303375

See this dsmcFoam+ published on Computer Physics Communications.

5.5.1 fluentMeshToFoam
......
There is currently no support for embedded interfaces and refinement trees.
貌似对这个是不支持的。我也用过snappyHexMesh， 这种refinement的过渡处是会生成非结构网格的，也就是没有interface。现在snappyhexmesh是做了一些改进吗？

我的理解是fvc::div(phiHbyA)已经将原控制方程中的散度表示为高斯积分的求和形式了

是的。正确。