multiphaseEulerFoam--blending功能理解

ryjiao

各位老师好，关于 multiphaseEulerFoam--blending的实现过程，关于BlendedInterfacialModel.C中提及的两个量{model_.valid() 与subtract}不知该如何理解，想烦请老师与大神们指点，万分感谢。
另附OpenFOAM-v8中BlendedInterfacialModel.C的68-93行代码，再次感谢~

 if (model_.valid() && subtract)
 {
     FatalErrorInFunction
         << "Cannot treat an interfacial model with no distinction between "
         << "continuous and dispersed phases as signed"
         << exit(FatalError);
 }

 if (model_.valid() || model1In2_.valid())
 {
     f1 =
         blendedInterfacialModel::interpolate<scalarGeoField>
         (
             blending_.f1(phase1_, phase2_)
         );
 }

 if (model_.valid() || model2In1_.valid())
 {
     f2 =
         (subtract ? -1 : +1)
        *blendedInterfacialModel::interpolate<scalarGeoField>
         (
             blending_.f2(phase1_, phase2_)
         );
 }

ryjiao

各位老师与大神们好，我想再补充一下我的问题。
首先对于model_.valid() ，BlendedInterfacialModel.H中的注释为“//- Model for region with no obvious dispersed phase”，我浅薄的理解为alpha_1=alpha_2=0.5；
而对于subtract，我初步理解为subtract=f1-f2，subtract=0时为false，非0时为true；
但按照以上两种理解，我不明白对于计算f2时subtract为何要在原f2上进行变号，也无法推导出类似M=f1M1+f2M2的blending公式形式。
我认为一定是我的理解存在问题，想请各位指点，非常感谢~