Posts made by vbcwl

@vbcwl 替代文字](图片地址)
这个是用的of7 decomposePar之后得出的，左边是process0-4右边是5-19，我感觉这个有可能是正确的，process5-19应该是没有分到外壁面的部分，所以是nonuniform 0（）

@李东岳 图没上传成功重新发下，
process0-4都是这种，5块的值都是-0.06788236085，

@李东岳感谢东岳老师回复，这个图是我设置的初始条件
用的scotch拆分成20块，没有设置任何参数

这个是拆分后的processor0-4都是下图这种
![3)HI(1U[8I96}(H@)D%(~X.png

processor5-19都是下图这种形式

有前辈遇到过类似的情况么？就是decomposePar分块后壁面边界条件fixedGradient uniform -100会变成别的值

请问一下各位前辈，我在设定温度T的边界条件的时候，壁面设置为
type fixedGradient
gradient uniform 10;
并行分块后就变成了 type fixedGradient
gradient nonuniform List<scalar>；
值的话就是壁面的初始温度值。
这是哪里出错了呢？

@小白好海牛 https://zhuanlan.zhihu.com/p/315686722

@李东岳 明白了，谢谢东岳老师

@vbcwl 我后面用grad求解出du/dx dv/dy dw/dz这些，相加确实等于这个div求解结果，但是du/dx 与(u2-u1)/(x2-x1)相差巨大，有好几倍，就算这个手算精度差也不至于这么大差距吧，请问有前辈了解么

请教一下各位前辈，我做的槽道流，体积平均速度为15.43，然后使用postProcess -func 'div(U)'求出来散度，值都在0.03 0.04这附近，甚至能到0.13，求解格式是div(U) Gauss linear，用的不可压缩pimpleFoam，这个散度这么大是什么情况？

@bestucan 感谢前辈回复，已经解决了

如题，望各位前辈指点

@bestucan 谢谢前辈，我去试试去

@yfleung 自带的有meanVelocityForce，可以维持恒定的体积(面)流速，但是这种方法压力不是定值。想要定值 https://www.cfd-china.com/topic/280/动量方程如何加入加压力梯度源项?_=1656051870046

@bestucan 前辈指的是这样么？切换为6，foam确认确实切换成功，然后回到算例文件夹，在并行？还是不行啊

@李东岳 这是我师兄之前安装好的，也是用的东岳流体上的方法。我是最近又用东岳流体上的方法安装了of9。安装完之后切换为6就无法并行了
![KZ~U2DKGQS$8AMDJC@}Q15.png 这是报错信息。

@bestucan 我是终端输入of6回车 输入 mpirun -np 4 pimpleFoam -parallel运行of6自带的槽道流算例，直接报错，如果不并行是of多个版本都能正常使用，但是现在只有最后安装的9可以使用并行

@bestucan
前辈您好，我这个第一种写法就是东岳流体上多版本共存的写法，第二种那里是我的bashrc文件的路径，可还是没有效果咋办啊？我有用of6切换为op6版本

请教一下各位前辈，我在有openFoam6的情况下，安装了of9，结果只要运行mpirun -np -parallel这个命令就会报错，我感觉应该是of9给我安装了更新版本的mpirun导致的，这个应该怎么解决啊？

@李东岳 谢谢东岳老师，现在正常了

请教一下各位前辈，我用of6的自带算例channel395计算Uprime2Mean，可是of9计算出来的结果，第1 4 6列为什么是负数？of6的结果貌似是正确的，是同一个算例文件，取得同时刻的结果，除了openFOAM版本没有任何区别

如图，后面的是of6，前面的是of9，of9数据大小与6差距不小，但为负数着实无法理解。

@vbcwl 明白了，是因为of9没有这个函数了，这个是之前的功能

@李东岳 东岳老师您好，我把这个加入到controlDict里后，终端和每个时间步的结果文件都找不到最大最小值，请问这个值应该在哪里呢？

@李东岳 谢谢东岳老师，我去研究下这个boxTurb去

@vbcwl

@李东岳 这个虽然不符合物理，但是达到湍流状态后就不会再加这个随机数了，只是为了尽快达到湍流状态。请问东岳老师从算法上应该怎么做呢？

@李东岳 没有。我现在是想，对每一步的速度结果Ux都随机加上一个值，这个值在-0.5到0.5之间，然后在继续运算，来使流动尽快到达湍流。我通过groovyBC添加的边界扰动效果不理想，请问下东岳老师这个OPenFoam自身有办法实现么？

@zhouxu 都可以，最好内部场

请问一下，我想给速度场加一个扰动来使流动更快的到达稳态，需要在前一千秒的速度场中逐步添加随机数，有前辈可以指点下应该怎么做么？

@李东岳 东岳老师可以重发一下湍流粘度的链接么？这个链接有点问题

@xpqiu 感谢前辈的回复，我用的laminar做qDNS，没有用到nut这些，只用到了p p_rgh U T。请问您了解过没有nut文件的情况么？

@李东岳 能请教下这个具体怎么使用么？我没能输出utau

@benqing 请问下这个是计算utau=sqrt(tauw/密度)里的tauw么？我求出来怎么全是负数

@xpqiu 您好，请问下我想利用utau来对距壁面距离和流体速度做无量纲化的，utau=sqrt(nuu0/y0) u0是壁面上方第一个网格的速度，y0是writeCellCentre命令写出的第一个网格的中心坐标，然后由于yPlus小于1，直接用mu来求解壁面剪应力，yPlus=utauy/nu,这样子算出来的yPlus与postProcess写出的yPlus有较大差距是什么原因？

请问下有没有前辈研究过meanVelocityForce或者知道哪里能找到它的详细解析的？我一直没搞明白他具体是什么原理，参考了东岳流体上的http://www.dyfluid.com/boundaryFoam.html boundaryFoam的解析之后大概了解他和这个类似，但是具体的公式推导是怎么样的还是不明白

meanVelocityForce只用于不可压缩，但是涉及传热又不区分可压缩与不可压缩，那么像buoyantPimpleFoam这种求解器可以用meanVelocityForce来驱动么？有了解的前辈解释下么？

@李东岳 东岳老师我在做和这个前辈类似的东西，请问您有他的论文或者算例文件可以给我学习下么？

@shrine 你多试几个，我也是新手，之前遇到过类似的，换几个总有能用的

@shrine 你去看看别的算例的fvSolution,有定义cellMot*****的复制过来就行了

@李东岳 东岳老师，我了解了以下结构化网格的排序，然后刚才通过postProcess -func writeCellCentres 写出了channel395的网格，这里网格是按照先x在y后z的顺序排列的，时间步中的结果应该是和这个对应的请问下东岳老师，有办法把结果输出成三维列表来么？

@李东岳 对的

@vbcwl 想把数据导出去求平均，本来以为是先x后y在z排列的，可是发现好像不是这么排列的

openfoam的速度和压力这些都是输出的一维列表，这是具体是怎么排序的？或者可以怎么转化成存储为三维列表呢？

@浪迹天大 请问下这个怎么使用啊

请问如何同时对流场做流向和展向的平均呢？我现在只会在paraview里通过切出平行于流向展向的平面来做平均，但是这样的话想要画出速度沿法向的分布图需要切很多个面，太过于费时间了，有没有什么简单点的方法呢？

@yfleung 请问下用paraview该怎么做流向和展向的平均呢？

请教一下各位，我做的圆环柱的传热，流向为周期性边界条件，内壁面恒定热流，外壁面绝热，内外壁面无滑移，fvOptions加了一个meanVelocityForce，（0.7 0 0），0.7就是流向，可是最终的速度云图，产生了一块速度较大的区域，并且该区域还会沿环向流动

如图，可以明显看到红色高速区在缓慢的逆时针移动，请问下为什么会产生这个高速区，以及为什么高速区会移动？

@bestucan 试出来效果差不多，但是我水平不足以判断正确与否不敢使用哈