CFD中文网

@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进行类似于面积平均，最终得到某个方向的分布。那我不明白了，只对一个方向进行平均，这该怎么改呀？这原理都不一样了？？？

@aurora16 这些因素已经排除过了，况且你说的这几个点也不会导致我上面提到的问题

最近一直在尝试用SRFpimpleFoam求解器算一个旋转槽道流，如图1，可是发现计算一直发散，经过再三尝试，目前我测出来的结论如下：

计算参数，进口主流0.1335m/s，旋转轴为Z轴（0 0 1），旋转中心为原点（0 0 0），转速rpm=120，槽道几何比例为4:2:2，我尝试了不同scale下旋转槽道的计算，发现有问题。

首先，计算域为0.04m-0.02m-0.02m时，旋转状态下是可以算的，监测通道进口、中心点相对速度，都可以稳定。但是当计算域为4m-2m-2m时，通道进口、中心点相对速度就无限增大，然后发散了。

经过各种测试，排除原因，我得到一个结论，就是这种情况是离心力导致的，SRF求解的动量方程如下：

离心力已经写成势函数形式，我发现上面的计算发散就是由于离心力太大导致的，经过测试，有三种情况会导致发散，
1转速omga很大
2旋转半径很大，就如前面的0.02m变为了2m的旋转半径，导致离心力很大
3旋转中心里通道进口很远，如果我把旋转中心从0点调为(-5 0 0)，它也会发散
综上原因，就是当离心力很大的时候，难道计算就无法稳定了吗？？？

@学流体的小明 说错了，我的意思是v=w×r，这个算的是圆周速度，而我们看的是相对坐标系下的相对速度Urel，从NS方程来说它多受到了离心力和哥氏力，哥氏力和主流不在一个方向，而离心力是有势力，如果我进口主流速度都一样的话，不管旋转半径多大，与进口差距一定距离的点，其受到的离心力做功都是一样的，所以理论上改变旋转半径大小，那个点的相对速度应该基本是一样的吧。。。。

@学流体的小明 你说的这个是绝对速度，看相对速度的话，和转轴没多大关系
我突然有个问题，目前算槽道流都是按6.4h-2h-3.2h比例算的，那我是不是可以计算域改成0.64m-0.2m-0.32m，和计算域为6.4m-2m-3.2m是一样的吧，就是只要比例一致就行了嘛？？

调了半天还是不行。。。。网格不管加密多大，都算不了。
channel395 网格尺度为0.04m×0.02m×0.02m，旋转中心为x=-0.2m，转速120rpm，设置进口来流0.1335m/s。probe监测x=0.02m，y=z=0.01m处，计算可以稳定，稳定后速度为0.160271m/s
当保持上述设置，只把旋转中心为x=-2m，计算没法稳定，速度在1~2m/s波动，显然不正确。
按理说离心力是一种有势力的话，他只和起点终点有关系，我改了旋转中心应该不影响离心力加功的多少啊？？？上述速度都是相对速度Urel哈

进出口还都设置的是普通的patch条件，不是周期性条件，已经排除了离心力做功这个因素
@李东岳 老师

如题，最近我用SRF求解器算DNS，同样的进出口条件、网格大小，以channel395为例，原始的网格尺度是4m×2m×2m，主流速度0.1335m/s，设置旋转速度rpm=120，这样计算出来使用probe监测任意一点速度，发现相对速度很离谱，随着时间一直升高，最后都10的6、7次方相对速度了。
反而，我在blockMesh里把网格scale设置为0.01，也就是计算域变成了0.04m×0.02m×0.02m，其他的设置都一样，然后计算就能收敛稳定了，如下图是通道中心相对速度随时间变化：

这到底是为什么呀？？？其它条件都一样，就改了网格尺度就算不了了，求大神解答~