icoReactingMultiphaseInterFoam中没有组分浓度梯度的扩散项

[xr]

openfoam中的icoReactingMultiphaseInterFoam求解器应该是用来求解多相多组分的问题的，组分浓度Y的求解使用的是其中链接的库函数src/phaseSystemModels/multiphaseInter/phasesSystem/phaseModel/MultiComponentPhaseModel/MultiComponentPhaseModel.C
其中YEqn的求解如下所示

想请教一下各位前辈这个里面为什么不包含组分扩散项呢，就是比如使用最简单的Fick定律计算得到的用于浓度梯度的扩散项，并且也没有相关的输入质量扩散系数D的地方，请教一下各位前辈

[李东岳]

我没有太详细看这个求解器。从这个方程来看，看起来没有植入。同时，这是个迎风格式的对流。精度也比较低。我不确定是否在其他地方有进行相关处理。

[xr]

谢谢李老师，我看起来也觉得确实没有处理，不知道是不是在其他的库函数中。看看其他的前辈有没有了解的，非常感谢

[lrl3512]

组分扩散在最下面啊
在constant里面有是否打开扩散的开关

[xr]

@lrl3512 谢谢前辈， 这个里面的计算用的是湍流粘性除以湍流Sc数，这部分难道不是只包含了湍流引起的扩散吗，但是分子之间的扩散是不是没有包含呢？

[lrl3512]

@xr 你说的好像对 我说算出来怎么扩散不明显

[lrl3512]

@xr 你对 https://cfd-china.com/topic/6398/关于icoreactingmultiphaseinterfoam中的组分方程 怎么看

[xr]

@lrl3512 谢谢您的回复，对于您的问题，具体我也不太清楚但是看这两个方程的话，如果用带相的方程扩散项怎么写呢，就是分子扩散相，乘D的话在界面处会不会存在一些误差呢

[lrl3512]

@xr 我昨天又重新看了一下代码，发现Yi并不是质量分数，而是体积分数，所以如果你直接在下面扩散方程里面加扩散系数，应该是有问题的

[xr]

@lrl3512 谢谢您的回复，但是总的来说就是他确实没有分子扩散项还是得加这一项的对吧，感觉是个大工程

[lrl3512]

@xr 如果你想继续研究的话，我们可以加一下联系方式，一起研究一下

[xr]

@lrl3512 好的

[xr]

@lrl3512 这个怎么加联系方式呢 我有点没太用过这个 我不知道怎么私聊

[xr]

@lrl3512 要不您加一下我qq 2559727310