有老师在多组分的条件下使用过质量源项吗
-
最近在使用OpenFOAM中的质量源项,发现了一些奇怪的情况,想跟各位老师交流一下
massSource1 { type scalarSemiImplicitSource; timeStart 0; duration 10.0; selectionMode points; points ((0.3 0.3 0.3)); volumeMode absolute; injectionRateSuSp { Air (1.0 0); rho (1.0 0); } }
以上述源项为例,我流场中初始包含两种气体,一个是Air,另一个是水蒸气。流场完全封闭,我在流场的中心位置放置了一个空气的源项,发现:
①质量不守恒,在我添加了空气源项后,水蒸气的质量会减少,但是不应当减小(我没有启用什么冷凝蒸发模型,这里就当是两组分气体就行)
②如果Air的输入值大于rho,那么流场中会出现质量分数大于1的情况,且最大值等于输入的Air和rho的比值,在Air有定义,而rho无定义时,该值会趋于无穷(1e15)(paraFoam导出的流场中,Air场最大值大于1);
③如果不输入Air,仅输入rho,那么还是空气的质量会增加,我猜测这与我定义的背景气体是Air相关。(热力学模型为multiComponentMixture,需要定义inertSpecie)
我的问题是,OpenFOAM的这个质量源项的定义,在多组分下应该如何正确的定义(求解器为chtMultRegionFoam,是有多组分方程求解的)。现在我尝试了很多,但是没有一个是满足质量守恒的,实在难受 -
质量增加了,体积不变,密度是应该跟着变的。
Ar 和 rho 之间应该有个关系。以单一网格为例,密度的变化 就是 质量变化 除以 网格体积。
- 添加了空气源项,水蒸气质量减少。因为密度没增加,总体质量没增加,强制增加其中一个组分的质量,另一个组分的质量可不就下降了。总体质量守恒。
2、3. air 和 rho 的比例导致的质量分数问题,也是总体质量守恒造成的。质量分数大于1的时候。分母是总体质量(原本质量 + rho 源项),分子是组分质量(原本 Air质量 + Air 源项)。上下对不上,可不就偏离 1 了。大概这个意思,代码里兴许算的时候会有些系数什么的
我猜的