OpenFOAM中，默认的多组分方程全部有误？

[四季之中的隐星]

一般而言，组分方程中的拉普拉斯项，其广义扩散系数为D，即一种流体在另外一种流体中的扩散系数。然而，在OpenFOAM中的组分方程中，没有一个使用的是D，均使用的是alpha（热扩散系数）和mu（粘度，或者叫动量扩散系数）。

源代码中：
在chtMultiRegionFoam的组分方程中

            fvScalarMatrix YiEqn
            (
                fvm::ddt(rho, Yi)
              + mvConvection->fvmDiv(phi, Yi)
              - fvm::laplacian(turbulence.muEff(), Yi)//注意，这里调用了粘度
             ==
                reaction.R(Yi)
              + fvOptions(rho, Yi)
            );

在fireFoam的组分方程中

fvScalarMatrix YiEqn
            (
                fvm::ddt(rho, Yi)
              + mvConvection->fvmDiv(phi, Yi)
              - fvm::laplacian(turbulence->alphaEff(), Yi)//注意，此处为热扩散系数
              ==
                parcels.SYi(i, Yi)
              + surfaceFilm.Srho(i)
              + combustion->R(Yi)
              + fvOptions(rho, Yi)
            );

那么，为什么好像从来没有人提到或者询问这个问题（甚至可以说是很根本的问题），还是说我对代码的理解有误。望各位老师和同僚赐教

[wangfei9088]

@四季之中的隐星
$D_k$ = $\alpha$/$Le_k$，$\alpha$是热扩散系数，$k$是组分。
OpenFOAM中默认各组分Lewis number都等于1。
同样，各组分Schmidt number也都默认等于1。

[四季之中的隐星]

@wangfei9088 好的，谢谢 那看来要用OpenFOAM做一些对流扩散问题的研究还有很多需要修正的地方

[wangfei9088]

@四季之中的隐星
是的，根据你的需要进行修改。因为多数时候用不上，或者仅用作演示，所以都精简了。