TwoPhaseEulerFOAM的一些小问题

[yfclark]

我正在准备把在这个Foam上面植入拉格朗日模型，有一点小问题
（1）求解器里面似乎没有求解连续性方程，把连续性方程带到了动量方程里面，不知道我的理解是否正确
（2）因为要植入拉格朗日模型，Euler方法里面的曳力模型，粒子判断什么都要关掉，phaseProperties不设置任何东西，是不是就认为关掉了@东岳
@zwl 老哥似乎做过空化的，连续性方程有源项是不是需要处理后用到动量方程的求解里面

[zwl]

@yfclark 源相是需要加到压力泊松方程里面，当然由于源相的产生导致的动量源相和质量源相也要加到对应方程中的。不设置任何东西应该是不考虑曳力模型的。不清楚你要怎么加拉格朗日模型，为什么不考虑相间曳力？

[yfclark]

@zwl 谢谢老哥，我的设计是连续的两相用Euler-Euler方法，第三相是离散相，用拉格朗日方法，第三相与前面两相有质量和动量交换，Euler和Euler之间没有任何关系，所以我打算关掉这些模型

[东岳]

@yfclark 在 TwoPhaseEulerFOAM的一些小问题 中说：

Euler和Euler之间没有任何关系，

这是什么工况行为/物理行为？

[yfclark]

做雾化啊，连续的气相和液相用Euler-Euler，两相没有掺混，就没有质量，动量交换，液相破碎生成的液滴用拉格朗日追踪，更@zwl 老哥的那个空化近场雾化很相似，不过他用的是Euler追踪生成的液滴和气泡，因此Euler-Euler之间需要引入动量和质量交换

[yfclark]

这个连续性方程是怎么引入的，质量方程没有源项的话，就应该是0吧。

我对比的是OpenFOAM2.3.x的代码：

        U1Eqn =
        (
            fvm::ddt(alpha1, rho1, U1) + fvm::div(alphaRhoPhi1, U1)
          - fvm::Sp(contErr1, U1)
          + mrfZones(alpha1*rho1 + virtualMassCoeff, U1)
          + phase1.turbulence().divDevRhoReff(U1)
         ==
          - liftForce
          - wallLubricationForce
          - turbulentDispersionForce
          - virtualMassCoeff
           *(
                fvm::ddt(U1)
              + fvm::div(phi1, U1)
              - fvm::Sp(fvc::div(phi1), U1)
              - DDtU2
            )

@东岳

[东岳]

@yfclark 在 TwoPhaseEulerFOAM的一些小问题 中说：

这个连续性方程是怎么引入的

http://dyfluid.com/icoFoam.html 我建议你看看这个，双欧拉模型跟这个是一样的

[yfclark]

@东岳 非常感谢东岳老师，是我的问题没有问清楚，我只是好奇求解动量方程时引入的fvm::Sp(contErr1, U1)即
只求解一个动量方程的求解器都没有做这样的处理，因此有些好奇
正在研究OpenFOAM2.3的twoPhaseEulerFoam，有问题还要请教东岳老师，请不要见怪

[yfclark]

@东岳 @zwl 您好，我想问一下您写求解器的时候表面张力是如何处理的，我看multiphaseEulerFoam和reactingmultiphaseFoam都有surfaceTension的计算，但是在momentumtransfer里面却没有表面张力，我现在的思路主要是两种：
（1）在动量方程中左侧直接植入表面张力
（2）在动量方程中忽略表面张力，使用momentumPredict添加表面张力

[东岳]

@yfclark 处于有界性的考虑，https://openfoam.org/release/2-2-0/numerics-boundedness/

[yfclark]

@东岳 明白了，谢谢东岳老师