CFD中文网

@李东岳 好的 谢谢东岳老师建议，我去研究一下这个文章

下面为一个名为hybridCentralFoam的可压缩两相流求解器的动量方程UEqn代码，我仔细阅读了一下，发现在方程等号右边有重力和表面张力项，但没有压力梯度项，即类似 fvc::grad(p) 这样的形式。个人觉得奇怪，这样不就无法描述相界面压差对流动的影响了吗？ 恳求各位大佬指点一二。
求解器方程这块我是刚开始研究的小白，如果我想要加入压力梯度项的话，具体要怎么加入呢，别的地方还需要修改什么吗？

void Foam::vofTwoPhaseCentralFoam::UEqn()
{
    const auto& rho = mixture_model_.rho();

    surfaceScalarField phiU_own
    (
        vF1face_*phi1_own_ +
        vF2face_*phi2_own_
    );
    surfaceScalarField phiU_nei
    (
        vF1face_*phi1_nei_ +
        vF2face_*phi2_nei_
    );
    phiU_own.rename("phiU_own");
    phiU_nei.rename("phiU_nei");

    E_ = fvc::ddt(rho) + fvc::div(phiU_own) + fvc::div(phiU_nei);
    fvVectorMatrix UEqn
    (
          fvm::ddt(rho,U_)
        + fvm::div(phiU_own,U_) + fvm::div(phiU_nei,U_)
        - fvm::Sp(E_,U_)
        + turbulence_->divDevRhoReff(rho,U_)
    );
    oneByA_  = 1.0/UEqn.A();
    HbyA_ = UEqn.H()*oneByA_;
    HbyA_.boundaryFieldRef() == U_.boundaryField();

    B_ = fvc::reconstruct
    (
        (
           interface_.surfaceTensionForce()
            - ghf_*fvc::snGrad(rho)
        ) * B_.mesh().magSf()
    );
}

@瑞水南 您好，请问层流里如何制定的这个问题解决了吗，能不能分享一下

@Shihang-Chen 十分感谢，我去查一下

@Zhujh @Shihang-Chen 两位大佬好，请问你们用到的CLSVOF 方法有相关的代码资源吗，我最近也打算结合CLS和VOF方法，但是没有找到合适的切入点，openfoam这块属于刚开始学

@mark 请问楼主是如何植入相变的呢，能向您请教一下吗，或者有什么相关文献分享一下

@Joker 好的，感谢大佬建议，我去试一下这个

明白了，谢谢东岳老师，我接下来试试3D模型。老师那您觉得是用层流还是湍流比较好呢

@李东岳 东岳老师，这个case是2D的，液相是VOF模型，这个求解器之前用来求解激波下的液滴演化破碎，默认湍流模型是层流，我现在换成k-epsilon，但是一直算的不合理，调了各种参数，但是没有动过fvScheme和fcSolution这两个文件，因为了解的不够，还请您指点

@李东岳

感谢东岳老师回复，这是湍流粘度云图，液膜射流入口给的是calculated条件，然后射流下游的壁面给的是壁面函数，这个分布我看的不是很明白，您能帮忙分析一下问题吗，十分感谢东岳老师！
下面是湍流粘度的边界条件

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0.0001;

boundaryField
{
    UPPER-WALL
    {
type            fixedValue;
      value           uniform 0.0001;
        
}

    F-INLET
    {
     type            calculated;
  value           uniform 0.0001;
         
    }

    OUTLET
    {
       type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

   BOTTOM-WALL
    {
      type            nutkWallFunction;
        value           uniform 0.0001;
    }
    C-INLET
    {
       type            calculated;
        value           uniform 0;

    }
    frontAndBackPlanes
    {
        type            empty;
    }
}

@Joker 我按照超音速来流算的y+，特征长度取的是计算域的x方向长度的10%，跟目前网格的第一层厚度基本上一致，

在网上找到了一个叫hybridCentralFoam的可压缩两相流求解器，用于模拟激波下液滴破碎，默湍流模型为laminar，我打算用该求解器来计算超音速气流下液膜演化的问题。第一张是用求解器默认的laminar湍流模型计算的液相体积分数结果，并没有出现液膜在超音速来流下扭转破碎的现象

因此我将湍流模型换成了k-epsilon，得到了下面的结果

液膜出现了一定程度的破碎，但奇怪的是变薄了很多，感觉跟层流差异特别大。
我是按照公式算的k,nut,的大概初值，求助各位大佬

@Joker 现在的问题就是文献对不上，只是有类似的样子，但是具体的破碎特点、温度、速度分布都有差异，这个问题我自己再研究研究。您说的湍流边界条件是这样的：
入口的k,omega都是用公式算的fixed value，nut 给的是calculated条件
壁面都用的是壁面函数
出口的k,omega是zero gradient，nut 是calculated。
目前我只根据湍流模型修改了0文件，fvScheme和fvSolution我没有改动，因为这块还是不太懂

@Joker 目前我在用hybridCentralFoam这个求解器算超音速下液膜演化问题，就是分析超音速来流对下壁面流出的液膜的影响，这个求解器本身默认的湍流模型是laminar，下面是算出来的液相体积分数，红色就是液膜，很平滑，没有体现出超音速气流对液膜的剪切破碎作用。

所以我打算把laminar换成看k-omega模型，算出来的液膜是下面这样的
感觉液相突然间少了许多，就像凭空消失了一样，而且温度场也跟层流下的差异很大
我是按照公式大概估算了入口的k,nut,omega这些值，所以不知道出了什么问题

@Joker 可以讨论一下吗，我在做液膜射流问题，也打算换k-omega做，但是初始值边界条件的确定遇到了点问题

主流速度较低的时候还是有液膜破碎分离的，上图是300m/s的工况，但是速度提高之后破碎现象就渐渐消失了，变得十分平滑

@李东岳 是的，李老师，用的是vof方法，抱歉我是初学小白，很多问题确实不太明白，现在在试着把湍流模型的laminar改掉

@李东岳

李老师，我已经把phaseproperties里面的液相粘度增加到10-3量级了，还是没啥效果，需要把湍流模型从laminar改成k-epsilon或者别的吗？

@李东岳 感谢李老师建议，我试下网格加密。 另外您说的没有考虑粘性，也就是不会出现剪切吧，那如何将粘性考虑进去呢，但是我读这个求解器来源的文献里讲的是可以模拟可压缩两相粘性流体，DOI: 10.20944/preprints202205.0412.v1

我看到phaseProperties文件里，液相的粘度是1e-6，比气相还要小，这点也很让我纳闷，还希望李老师能给些建议

目前在用 hybridCentralFoam 这个可压缩两相流求解器计算超声速平板液膜射流问题，正常来讲液膜会在超声速主流的剪切作用下发生波动直至破碎，形成液丝、液滴等，但是我算出来的液膜出流后十分平滑地覆盖在下表面上，几乎没有任何不稳定波动，下面是液相分数场和温度场。求助各位大佬，能否给予指点。


另外，我看了下该求解器的湍流模型选的是laminar，难道是这个原因吗，要怎么改进呢