CFD中文网

这里如果使用SMALL，当接近零的时候，Y_s的值会接近零，滑移速度为0，也是符合原理的。但是我不确定OpenFOAM在计算时是否Y_s会真的趋向于0，请大佬们看看是否有更好的处理方法

在单流体mixture模型中，植入曳力模型中除以0的问题，想问下该怎么处理：
整个框架中使用omega-Y方程计算液滴粒径

再使用计算得到的粒径计算曳力，求解滑移速度


但是在纯液体以及纯气体中会出现液滴粒径为0的情况，这就导致其中nd(网格中的粒子数目)会出现为零的情况，想请教下这种情况下直接使用SMALL有用么，还是需要别的特殊的处理方法，非常感谢~

就是在存在Y=1的区域的时候，如果Ur不合理，会导致发散

方程19处理成固定值是指Ur的大小方程都恒定么？这个时候不会出现震荡，但是遇到了一个问题，就是Y在计算时会超过1，应该算是过度堆积？这种问题应该怎么处理呢~ 这里我也使用了MULES进行限制。
MULES的调用是这样的

void Foam::mulesWithDiffusionImplicitLimiter
(
    const volScalarField& rho,
    volScalarField& Y,
    const surfaceScalarField& phi,
    scalarField& lambdaFace,
    surfaceScalarField& rhoPhif,
    surfaceScalarField& diffFlux,
    const surfaceScalarField& Dmi,
    const fvScalarMatrix& SuSp
)
{
    const fvMesh& mesh = rho.mesh();
    const word Yname (Y.name());
    Y.rename("Yi");

    upwind<scalar> UDs(mesh, phi);

    fvScalarMatrix YConvection
    (
        fv::gaussConvectionScheme<scalar>(mesh, phi, UDs).fvmDiv(phi, Y)
    );

    surfaceScalarField rhoPhifBD (YConvection.flux());

    surfaceScalarField& rhoPhifCorr = rhoPhif;
    rhoPhifCorr -= rhoPhifBD;

    volScalarField Su
    (
        "Su",
        SuSp & Y
    );

    MULES::limiter
    (
        lambdaFace,
        1.0/mesh.time().deltaTValue(),
        rho,
        Y,
        rhoPhifBD,
        rhoPhifCorr,
        zeroField(),
        Su,
        oneField(), //psiMax,
        zeroField() //psiMin,
    );
    Y.rename(Yname);
}

我参考的是这篇文章~
1-s2.0-S0301932221001397-main.pdf

现在植入的Udm的方程是这样的：
Udm =rhoc * alphac / rho * UdcMag * A /mag(A+epsilon)*
( exp(-a1_*max(alphad_ - residualAlpha_, scalar(0)))

• exp(-a_*max(alphad_ - residualAlpha_, scalar(0))) );
  其中A是加速度
  出现这种类似于棋盘一样的状态，可能是因为这里使用的参数的问题么~

请教个问题，就是再滑移速度模型植入时，出现了类似于棋盘一样的相对速度Udm，想问下这种问题该怎么处理呢~

求解器更换成了twoPhaseMixingCentralFoam（Kraposhin Matvey的求解器https://github.com/unicfdlab/hybridCentralSolvers），在其中的YEqn中添加了滑移速度Udm的通量，但是出现了Udm的棋盘一样的变化，不清楚问题出在了哪里，请大佬指条明路~
Udm植入时使用了固定大小的Udc，Udc的方向由加速度方向决定（认为气液两相惯性的差异导致速度的差异）。
这个算例在孔内是50%质量含气率的气液两相，入口压力10MPa，温度383K，（这是个测试，并非真实场景）

我在计算的课题是欠膨胀气液两相流内的激波，想看看滑移速度的影响，有没有可能计算的场景的瞬态特征导致了计算的不稳定啊~

这个模型是不是存在一个问题，就是滑移速度在计算时除了和加速度相关，与雷诺数相关，似乎与别的参数没有什么关系。这就会导致随着加速度的增加，滑移速度急剧增加，这里是否还缺少能量上的约束以避免出现不合理的结果呢

把源代码放上来，我用的是v2312版本的OpenFOAM，tutorials里面有参考算例
compressibleDriftFoam.zip

在将这个方程植入到OpenFOAM中之后，出现了明显的相对速度的震荡，惯性加速度最大能到10^5m/s2。植入时直接使用- fvc::ddt(U) - (U & fvc::grad(U))进行计算，使用的时1D的气液激波管问题做的验证，此时无滑移时的速度为20m/s左右，

给定滑移速度时计算的速度场没有震荡，使用的是OpenFOAM中的general方法。

想了解下出现这种问题的可能原因是什么呢~

谢谢李老师，我先试试~

fluent中提供了不同的滑移速度的计算方法，我想尝试将其移植到OpenFOAM中去，在处理滑移速度的时候遇到了问题：相对速度计算时需要用到相对雷诺数，相对雷诺数计算时需要用相对速度，有点套娃的感觉。这里在计算时需要使用迭代计算的方法么，有没有资料可以参考呢~ 先行谢过各位大佬~

关联的问题，线没有关联好，放到了wedge类中了

我先用一维气液激波管试试看看会不会有显著的震荡，谢谢李老师详细的解答~

了解了，密度梯度那里是个点乘，v1m·grad(rho1)，在这种情况下是不是也有可能出现有激波时的震荡啊，因为直接使用了密度的梯度~

了解了，谢谢李老师~为什么右边的密度梯度是错的呢~

我大概理解李老师的意思了，我对OpenFOAM中的MULES的原理理解的不透彻，李老师的意思应该是直接基于第一个方程： 直接求解此方程，漂移对流项放到前面去。
我这里处理的时候是想同时使用这两个方程：
最终构造成
，把漂移的影响分开到了源项和对流项中，需要同时修改alphaEqn和alphaSuSp。您觉得这样处理有必要么~
73348540-9b27-449d-9432-26ec0c6fe88b-alphaEqn.H d6ed1b9c-4078-46e3-b212-f7ed1662c6c0-alphaSuSp.H

是的，所以我参照着compressibleInterFoam中的方法对方程做了调整以使最终形式编程\partial alpha/ \partial t+ alpha\div vm的形式

采用这样处理的原因而不是直接使用$\partial(\alpha*\rho)/\partial t$是想直接套用compressibleInterFoam中的MULES方法，保证alpha的有界性