CFD中文网

好像去掉平均流没啥变化。。。现在就是不知道怎么确定求出来的能谱是否正确？
导入的这个速度场已经是稳态的湍流了，无论是对比平均流还是budget项，都和别人的数据能对的上

@coolhhh 对，是三维能谱。我刚才发现，能谱计算应该用脉动速度场，我再去试试

搞了几天总算是把代码调通了，但是看起来不是很对。。。。。输入的是瞬时速度场，计算出来的能谱。。。

@尚善若水 哦哦，一般计算能谱都是OF里先输出速度场，然后再去python或matlab里计算的

@coolhhh 谢谢大佬，我去仔细研究一下。前面看一些书和网页，写的太抽象了根本看不懂

各位大佬 。。。。虽然问题很简单。。。。@李东岳 @尚善若水 @coolhhh

各位大佬，最近在学能谱，有一些问题实在没搞懂

1、一维能谱和三维能谱区别就是一维输入的速度场只有u一个方向，三维是输入uvw三个方向吗？它们计算得出的曲线也只能是各自方向的能量随波数k变化的吧？

2、能谱分析输入的是一个任意大小的速度场吗？就是直接输入我网格对应的整个速度场？输入的速度场是脉动速度场吗？

@RQLI_CFD 在 请教：第三代涡识别方法Liutex的代码在linux中运行的问题 中说：

@lizhisongsjtu 感觉这个liutex没有那么好啊，我计算的是槽道流，liutex识别的结果很混乱（liutex的等值面），感觉还没Q准则显示的涡结构清晰，但是看刘老师文章里的效果挺好的。。。

您好，我也遇到这个问题了，Q准则识别的槽道流涡结构比Liutex好太多了，Liutex识别出来的涡都是断断续续的，怎么调阈值都不行。请问您解决了吗？

参考：https://www.cnblogs.com/Xiwang-Sun/p/17463243.html#:~:text=通过bloc

各位大神，帮帮我吧@李东岳 @尚善若水 @coolhhh @学流体的小明

各位大佬，有一个问题困扰我好久了，就是已知离心力Fcentrifugal=- fvc::grad(CenPot)，也就是离心力可以写成势函数CenPot=0.5omega^2(x^2+y^2)的梯度格式，目前我已经验证过CenPot这个量求出来是没问题的，关键就在于OF里面求一次梯度直接就求错了，求得非常离谱。
理论上求梯度，x方向梯度应该是omega^2*x，y方向梯度应该是omega^2*y吧，omega只有z分量，并且已经排除了网格、边界条件问题，网格用的结构化质量没毛病。求出来y方向梯度是正确的，x方向就不对，结果如下：

上图是-fvc::grad(CenPot)求出来的，居然是对称的？？？

这是标准的离心力，非常正确越往右边越大

这是标准的离心势，也是没毛病
那为啥就是OF算梯度算出来不对呢？？
这是我计算时梯度离散代码：无论是 Gauss linear还是leastSquares算出来都不对

gradSchemes
{
    default         Gauss linear;
}

这是我求解器源码：

Info<< "Reading field CenPot\n" << endl;//读取离心力势
volScalarField CenPot
(
    IOobject
    (
        "CenPot",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    1 * SRF->CenPot()
);

volVectorField GradCenPot
(
    IOobject
    (
        "GradCenPot",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    - fvc::grad(CenPot)
);

@Amadeus 是的，一个是Cp一个是Cv，但是目前看的各种版本推导的温度方程，有的用Cp有的用Cv不太一致。。。。

@Amadeus @aurora16 感谢各位回复，最后发现是由于对流-扩散方程中，离心力项太大了(高达几百量级)，而其它项还在10的量级，离散求解方程时，量级差异太大了，导致发散。计算还是要保证方程中各项量级差异不能太大

各位大佬，小弟我最近研究了一下CFD里面的温度能量方程，网上搜到了很多推导，但是我目前发现有量个版本，就是温度T的拉普拉斯算子项，有时候前面是k/Cp有时候是k/Cv，我看了好几个版本的推导，发现都没毛病呀？？？
1、从焓方程推导

把这个方程h=CpT代入，并且认为压力全导数Dp/Dt=0，就得到温度方程如下：

2、从内能e方程推导

把这个里面的e=CvT代入，温度拉普拉斯前面就变成了k/Cv了

到底该是哪个呀？？？？

各位大佬，小弟我最近总是出现这个问题，就是wmake报下面错误

mxy@LAPTOP-CODMI2H5:~/OpenFOAM/OpenFOAM-v2312/src/finiteVolume/cfdTools/general$ wmake
wmake: 'Make' directory does not exist in /home/mxy/OpenFOAM/OpenFOAM-v2312/src/finiteVolume/cfdTools/general
    Searching up directories tree for Make directory

有时候它就能自动找到上级Make文件

受不了了。。。。现在一个OpenFOAM上的问题就能卡我一个星期，进度太慢了，有没有大佬专家能帮我指导指导啊，可有偿。。。。。0.0

@尚善若水 这个程序我也下了，也看了一下，但也没有实现只对某一方向求平均。这个程序新增的功能就是
1、对时间目录中的所有字段进行平均
2、对边界面数据进行平均
3、Does not average across the channel centerline, so you get the full profile across the channel.这条没看懂啥意思
但是我也下载调试了，还是只能一下两个方向平均

@李东岳 在 OpenFOAM PostChannel只对一个方向平均该怎么改呀？ 中说：

那我不明白了，只对一个方向进行平均，这该怎么改呀？

变量的某一个分量么？

对于一个槽道流，就是只进行流向上的平均，或者只进行展向/法向平均。现有的PostChannel是直接对两个方向平均，比如说对流向、展向同时平均得到法向分布。

@李东岳 东岳老师救一下TuT

各位大佬，小弟我最近花时间终于看懂了Postchannel的原理，他首先对计算域分层，然后对每一层Layer进行类似于面积平均，最终得到某个方向的分布。那我不明白了，只对一个方向进行平均，这该怎么改呀？这原理都不一样了？？？