CFD中文网

请问各位大佬，对于可压缩流大涡模拟的总湍动能如何定义呢？

• 对于不可压缩，大涡模型采用的亚格子模型采用的雷诺时均的形式，因此可以用雷诺时均计算大涡湍动能和亚格子湍动能，两者的和就是总湍动能。
  $K=0.5*(\overline{uu}-\bar{u}\bar{u}+\overline{vv}-\bar{v}\bar{v}+\overline{ww}-\bar{w}\bar{w})+\bar{k}_{sgs}$
• 但是对于可压缩流的大涡模拟，亚格子模型采用的为法夫尔时均形式，此时如何定义总湍动能呢？我的理解是
  $\widetilde{K}=0.5*(\overline{\rho uu}-\bar{\rho}\widetilde{u}\widetilde{u}+\overline{\rho vv}-\bar{\rho}\widetilde{v}\widetilde{v}+\overline{\rho ww}-\bar{\rho}\widetilde{w}\widetilde{w})+\widetilde{k}_{sgs}$
  请问各位大佬，我的理解对么

这比较复杂了。一般系统会平衡CPU的负载，避免一核有难，多核观望的情况出现，所以算例实际调用的CPU资源和设置的核心数可能是不匹配的。另外，多核情况下，并行算法和分区规则对计算也是有影响的。

@李东岳李老师 您回复实在是太及时了
这是原文出处，Title:Large eddy simulations and analyses of hydrocarbon fuel heat transfer in vertical upward flows at supercritical pressures;
https://www.sciencedirect.com/science/article/pii/S0017931021000910

看了很多做圆管整管LES的论文，做1/4的确实是少见，不知道实际结果能差多少

今天看了一篇大涡模拟的论文，用的1/4 圆周做的大涡模拟，径向采用了周期性边界条件，感觉这种做法应该是欠妥当的。这样做相当于在周向上扰动同步，这样不就与湍流随机性相违背了么，而且靠近轴线附近的湍动能是不是也就变成0了呀。
有没有老铁做过对比测试，比较过1/4管和整管计算的差别？

@yingdong 遇“数”不决找原文，书和原文不一样，以期刊原文为主。 。P.S.好多系数都是经（cou）验(chu)性(lai)的，有些系数可能需要自己根据实际情况调整。

LienLeschziner低雷诺数湍流模型勘误。

Openfoam *.c文件中源码和出处

tmp<volScalarField> LienLeschziner::fMu() const
{
    const volScalarField yStar(sqrt(k_)*y_/nu());

    return
        (scalar(1) - exp(-Anu_*yStar))
       /((scalar(1) + SMALL) - exp(-Aeps_*yStar));
}

以下是原论文出处
Lien, F., and Leschziner, M. A. (December 1, 1993). "A Pressure-Velocity Solution Strategy for Compressible Flow and Its Application to Shock/Boundary-Layer Interaction Using Second-Moment Turbulence Closure." ASME. J. Fluids Eng. December 1993; 115(4): 717–725.

@桑葚熟了 第二种挺常见的，板式换热器的板片间的线接触和点接触一般都是做悬空

两种处理方法吧，一种是类似“黑滴”现象那种，取一个小段，线接触改成面接触。再一种就是不让他们接触了，悬空一段小缝隙。

@maud 你用的k-omega模型啊，那是不能直接得到C_D的，你可以用C_O去算C_D

@关伯兰在笑 楼主的代码是没问题的，我在自己电脑编译通过了。应该是编译环境的问题，试试重装个编译器吧，我一直用VS2010，比较稳定

我没做过这个方向，不过看这里这么多thread，是不是你的C_D里面的thread用错了？要不检查一下tm、tp和thread这三个指针？

在植入湍流模型的时候，需要调用热膨胀系数$\beta$, 其中beta是温度的函数$\beta=f(T)$。但是看了一下热物理模型类的所有equationOfState都不包含beta这一项。老版本boussinesqBuoyantPimpleFoam里面也只是定义了一个常数beta。请问应该怎么植入一个靠多项式计算的$\beta$呢？

@simthere 可以看看流量计压降信号转流量信号公式怎么算的，如果不修正悬臂压降的影响，用二维算是最经济的。如果里面有对悬臂压降的修正，那用二维是不合适的。说到底就是这个悬臂的压降能不能被忽略，如果不能，必须用三维

@ejector 可以用. ansys官方文档，"best practice:scale_resolving simulation in ansys cfd"里面有说明

@李东岳 谢谢，我再认真读一下

这是FLUENT中k方程浮升力生成项
$$G_{b}=-g_{i}\frac{\mu_t}{\rho Pr_t}\frac{\partial\rho}{\partial x_i}$$
把$\mu_t =\rho C_\mu k^2 /\varepsilon$ 带入上式以后应该是
$$G_{b}=-C_\mu g_{i}\frac{1}{ Pr_t}\frac{k}{\varepsilon}\frac{\partial\rho}{\partial x_i} k$$
但是openfoam buoyantKEpsilon.C中的形式却是
$$G_{b}=-C_\mu g_{i}\alpha\frac{k}{\varepsilon}\frac{\partial\rho}{\partial x_i} k$$
openfoam这里$\alpha=kappa/cp$么？感觉这里这个源项怪怪的。

另外。
下面是k方程源项的源码

buoyantKEpsilon<BasicTurbulenceModel>::kSource() const
{
    const uniformDimensionedVectorField& g =
        this->mesh_.objectRegistry::template
        lookupObject<uniformDimensionedVectorField>("g");

    if (mag(g.value()) > small)
    {
        return -fvm::SuSp(Gcoef(), this->k_);
    }
    else
    {
        return kEpsilon<BasicTurbulenceModel>::kSource();
    }
}

这里k方程的源项不应该是

    if (mag(g.value()) > small)
    {
        return kEpsilon<BasicTurbulenceModel>::kSource()-fvm::SuSp(Gcoef(), this->k_)
    }

@发芽的土豆 LES会过滤一部分小扰动转化为湍流粘性，想要自发生成扰动，就得保证自发的小扰动不被粘性耗散掉，否则流动会向层流发展。

做静电吸附还是范德华吸附？如果是范德华吸附，可以写一个概率模型。

@袁宝强 windows为了平衡负载，开启的线程数和实际调用的线程并不一致。

@李子橙 机器内存256G实际占用只有不到15%