CFD中文网

@东方白杨 波谱是速度-空间位置序列的自相关函数的傅里叶变换，数学上等价于序列傅里叶变换的模的平方(Pope的书里有证明)。能谱是对波谱在不同半径的球面上的积分。各项均匀同性湍流的时候，能谱就乘上对应半径的球面面积就行。

Zahtila, T., Lu, W., Chan, L., & Ooi, A. (2023). A systematic study of the grid requirements for a spectral element method solver. Computers & Fluids, 251, 105745.可以参考这个文章嗷

钝体绕流来说，升阻力的主频和流场监测点的速度主频应该是相同的。

@李东岳 好像不能是patch来着，必须要是wall类型的边界。

@黒洛克 wallShearStress存在对应“wall”的面上，不在体上。可以在paraview里直接输出对应面的数据再处理。

wallShearStress可以计算tau_w。然后用来算y+就好啦。

@李东岳 李老师，对于周期边界下的压力场，piso循环中有pEqn.setReference(pRefCell,pRefValue);用来确定压力场。那在计算周期边界下的被动标量时，应该用什么来对温度场进行修正呢？

@李东岳 !李老师，就是这样：

求得脉动速度后，直接做双内积，然后把这个标量场存下来，再用fieldAverage做时间平均。

@李东岳 是对压力边界采用fixedJump边界吗？
我先用层流稳态算例尝试一下看看~

@李东岳 东岳老师，jump边界上的问题解决了，可能是因为在并行计算的时候，没有将关联的周期边界放在一起引起的，加了：

这几句以后的计算结果种jump关联边界上的问题基本就没有了。

@李东岳 李老师您好~
求得时均速度之后，就可以用时均速度和当前时间步的脉动速度求瞬时的脉动速度应变率张量S’，可以直接对S’求双内积（&&）嗷，为什么要用这一步呢？

@李东岳 没用rhosimpleFoam是因为它解出来的好像是内能e，还涉及到热物性相关内容，不敢乱来hhh

@李东岳 但是就计算结果来说，确实实现了温度边界的阶跃，而且楼上有老师对二维算例的压力边界计算结果，效果也都还好。
就是不太明白为什么会出现这样部分网格计算结果不好

@李东岳 就是在simpleFoam求解器中添加了一个温度的标量输运方程

@李东岳 李老师，我看很多文献里，槽道湍流DNS计算时，除了时间平均，还可以做展向和流向的空间平均，在OpenFOAM自带LES Ret=395的算例中，有一个postChannel脚本可以改来用，个人觉得用了这个效果挺好的hhh~

@气固两相流 您好，我在simpleFoam里加了能量方程想测试一下fixedJump边界，是三维槽道算例，边界上会出现如下图的情况，在槽道出入口处有些网格上的结果感觉比较异常。

该算例中采用fvoptions添加动量源项实现动量方程的周期边界，因此速度的边界就是cyclic，在出入口处云图就很平滑：

除此以外，我在对一个横掠管束的算例中对压力采用fixedJump边界的时候，也出现了这种现象(压力场)：

但是在采用fvoptions添加动量源项时边界处结果就比较正常(虽然对管束这么实现周期边界是不正确的)：

还有就是采用这种边界，温度场好像会因为数值误差的累计(存疑)使得整体计算区域温度越算越高，下图分别是采用simpleFoam刚开始计算和很久之后的温度场分布，明显发现计算很久之后的温度场温度整体偏高，按理说稳态计算应该是不变的才对。(另：这个网格原来是用来算DNS的，所以我觉得算很久之后出现的这种温度场波动应该是算出湍流脉动了)：

这是我最近调试fixedJump边界遇到的一些问题，烦请大佬有时间帮忙瞅瞅~
如有描述不清的地方，还请多多指教！