CFD中文网

@calf-z-dns 你解决了吗?

我也在查 Reynolds Shear Stress, 和 Reynolds stress 有点分不太清楚.
http://brennen.caltech.edu/fluidbook/basicfluiddynamics/turbulence/reynoldsstresses.pdf (http://brennen.caltech.edu/fluidbook 公式的索引)
这篇上还有个 Reynolds normal stress ...

我的理解，对于可压缩流来说，求解的变量都是Favre平均量（密度除外，密度是时均密度）。所以，UPrime2Mean得到的也是基于Favre平均的脉动量，因为对于求解器来说它并不知道时均或者Favre平均，你让它算什么它就算什么。我看你贴出来的图的纵轴用的是两撇，这个一般就是Favre平均之后的脉动，因为一般时均后的脉动用一撇表示。

从这个角度讲，如果实验提供了Favre平均的相关统计量，我认为直接用模拟的结果去比对是没有问题的。如果实验只提供了时均的相关统计量，直接用模拟的结果去比对严格来讲不太正确。这里我也还没完全搞明白，看看有没有人来讨论了。

另外，我不清楚你这个画的是什么，只有一个类似的阶梯状的图。类似这个图我在今年8月份复现过，不过不是你这个算力，很多高阶格式都会把smooth的型线压缩成square的，例如Superbee，即使Minmod也有这种情况存在。

Figure 6.2(b) shows the same test problem as before but with the superbee method. The discontinuity stays considerably sharper.On the other hand, we see a tendency of the smoothhump to become steeper and squared off. This is sometimes a problem with superbee – by choosing the larger of the neighboring slopes it tends to steepen smooth transitions near inflection points.

Leveque 2002 p. 112

不过，我觉得你的问题可能不是这个...只不过图跟我遇到的有点像，给你点启发

我想了解下OpenFOAM是否可以批量截取截面？比如轴向有1000个网格点，就截取1000个面，然后计算每个面上的某物理量平均值。我看cfdpost可以联合Python进行批处理。我尝试了foamDataToFluent，但导入fluent流场都错了，就做不了批处理。

@星星星星晴 然而也没其他办法了吧:zoule:

@Calf-Z-DNS 好久没联系 你的DNS算的怎么样了

@calf-z-dns 你好，我也在做DNS的管道传热，壁面是恒定热流密度，可是我的温度云图结果几乎和层流一样，可以请教下你这种情况下该怎么解决么

@Calf-Z-DNS delta Y+的计算公式你已经写出来了，和delta X+是一样的。平板边界层的估算算法：https://www.cfd-online.com/Wiki/Y_plus_wall_distance_estimation

这个应该是超算那边给的测试内存不够引起的，我用少点的网格就可以正常跑。

另外，不知是否了解国家超算中心？具体哪家资源会更好一些呢？谢谢

我正在做加热管道空气流动例子的DNS，有新旧两个算例，其中旧算例的残差情况是：

初始.png

然后我把中间结果文件copy到新算例里面，接着继续跑，然后现在残差一直都比较水平：

中间.png

但是量级还是有点高，看了平均壁温，也没有达到均匀和期望效果。

我目前用的是基于buoyantPimpleFoam改编的求解器，设置成piso模式来跑，内迭代两次。U跟h的离散格式是QUICK.

不知道出现这种情况暗示着什么？有什么办法可以改变这种残差一直水平的情况呢？谢谢！

之前在安装openfoam5，以为安装ok了，但是在hpc提交脚本文件时，出现如下错误，rmall-.png

在cfdonline上面找到如下相关帖子：
https://www.cfd-online.com/Forums/openfoam-solving/216750-error-when-running-through-slurm.html
不知道是编译还未完成，或是需要打个补丁？还望解答，谢谢！

原实验论文里面把这个case称作‘subturbulent’，雷诺数是6020，T0 = 297.15 K,
P0 = 0.0927 MPa, D = 0.0274 m，V0=3.73 m/s，文中好像有讲到流动有趋于层流的趋势。

@东岳 请问东岳老师，可以用速度基求解器比如rhoCentralFoam求解这个算例么？还得加个温度去解压力？
没算过DNS，望指点！谢谢！

@Calf-Z-DNS 你好，看到你提的这个问题，请问你解决了吗？
我也想计算整个流动区域近壁面好几层网格的yPlus，而不是仅仅贴壁面的一层网格的yPlus，暂时我看到的wallShearStress和yPlus的计算都是贴壁面一层的值；
或者像东岳老师的示例代码那样直接给定了不变的uTau，而其实2维流动中，沿着流线方向uTau是变化的，所以我觉得比较复杂的是，怎么用某一壁面网格上的wallShearStress（或uTau）确定 对应法线方向上的多层网格yPlus，这个对应关系不知道如何寻找。
希望我表达清楚了，非常感谢！

@Dingcy 哈哈我是在用低雷诺数湍流模型算超临界CO2流动传热的例子，物性变化很剧烈，不好收敛。

@队长别开枪 谢谢指导。
这么看来，OpenFOAM对可压缩跟不可压缩的问题实现了代码的复用。
如果是不可压的，收敛后∇⋅U=0，2/3这项自动为0，也就省略了，不知道这么理解对不对。

那么对于动量方程2.png
这种形式的动量方程，没有体现2/3μ这一项，是因为不够严谨，还是因为这一项因为∇⋅U=0而省略掉？

如果是后者的话，那么这就变成不可压缩的动量方程？那为何还有密度这一项？

还是说这是一个统称的写法，具体还要根据问题判断。谢谢。

可以用吗？如何解决，求分享

externalWallHeatFluxTemperature，将q设为0.

不过我用这个BC加恒定热流密度，后处理用wallheatflux计算，发现有个min/max区间，integral的值比应该注入的能量少0.3W（48.8W VS 48.5W），不知道这是否正常。

问题已解决。
在.C文件中少了下面这段代码：

Foam::tmp<Foam::scalarField> Foam::rhotabularThermo::rho(const label patchi) const { return rho_.boundaryField()[patchi]; }

看看 这个 满足你的需求吗？