CFD中文网

看了各位大佬的评论，我觉得无论是什么样的矩形通道，应该都是用的最小那一边的半槽高度。例如进出口内外均为周期性的平板槽道，它就等价于一个AR=宽/高=无穷的矩形通道，这种情况不可能用水力直径的，应该都是用的最小半槽高度

各位大佬，小弟我最近遇到一个问题，已知摩擦雷诺数为：

目前的研究基本上都是对于一个平板槽道或者长=宽的方形管道，那这样特征长度取半槽高度h没问题，那对于一个矩形通道，比如说高y=0.2m，宽z=0.8m，它的摩擦雷诺数里面的特征长度该取半槽高0.1，还是水力直径呢？因为这个影响到了我压力梯度的取值
@尚善若水 @coolhhh @学流体的小明

哦哦 好的，感谢大佬分享这么多资料！

@coolhhh 大佬，我这几天又研究了一下，理解的差不多了，但是还有一个小问题没搞明白。按各种教材书上写的能谱计算，都是通过在傅里叶空间中按球壳（即波数的模值）积分来计算能谱，那这不就是都基于各向同性湍流了嘛？对于各向异性湍流也是用这套方法计算的吗？

好像去掉平均流没啥变化。。。现在就是不知道怎么确定求出来的能谱是否正确？
导入的这个速度场已经是稳态的湍流了，无论是对比平均流还是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