block的划分对质量分数是否有影响

[liujunCFD]

各位前辈好，有一个问题困扰我很久了。计算中途总是会报错 “质量分数总和大于1”，无论把来流设成空气（氧气，氮气，水）还是单组分气体（纯氧），都会出现这个报错，并且报错的都是三角形的上下两个block， 同样的情况下如果我关闭粘性和湍流模型则能完成计算。

我已经重新划分了block开始计算，目前已经算完了一半，暂时还没有报错，不知道大家有没有遇到过类似的问题？

然后还要谢谢东岳老师，观摩了这么久，终于可以来这里发帖了

模型如图所示

报错的block

计算停止的地方

报错信息

[liujunCFD]

@liujunCFD 对了，可能有的前辈会说三角形上下两个block的节点和相邻的点没有对上，这是因为靠近三角形的地方是边界层，另外就算对上了，也是会报同样的错，这个我已经试过了

[李东岳]

没见过你这个报错，你用rhoCentralFoam或者rhoSimpleFoam试试报错么

[liujunCFD]

@李东岳 老师好！我重新划分了block 就能完成计算了

话说我的计算结果还是不理想（下图均为密度）

刚开始看着还挺正常的

但是随着计算进行，激波不断向前移动

而且边界层好像被激波拉起来了？

壁面条件用的都是绝热壁面（考虑粘性），现在正在尝试壁面设成固定温度或者所有壁面不考虑粘性，不知道能不能改善结果

这个是原始文献的图

话说第一次见到李老师真人，激动

[李东岳]

你的激波看起来被耗散掉了。感觉是数值格式的问题。你用的什么高阶格式。也可以尝试去掉粘度看看会不会更尖锐。

另外你这个为什么不做结构网格，为什么要用交界面？

[liujunCFD]

@李东岳 老师好 我用的这个求解器虽然基于OF 不过还是有些差别 我用的是该求解器自带的网格处理工具，做的是结构网格，由“点→线→面”生成不同的block，然后block组合成整个几何，您说的交界面我不太理解是什么意思。网格代码比较长，如果老师对该求解器划分网格的方式感兴趣的话我可以把网格代码发上来。

这是我的一些参数选择 双破折号后面的为注释

另外像您说的，壁面去掉粘性后，计算结果正常多了

忽略左下角哈哈哈 那个block的壁面我忘记去掉粘性了