CFD中文网

@尚善若水 那请问您做的是单相问题吗，看着像是气体射流(个人瞎猜)

@尚善若水 对了老师，我想请问一下，您这里面的射流是气体吗，还是液体，我最近在用hybridCentralSolvers里面的vofTwoPhase算超声速液膜问题，但是算出来的液膜并没有被超声速来流剪切形变破碎，一直找不到原因

@尚善若水 谢谢老师，问题已经解决啦

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

下面为一个名为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改掉