反应流中Navier-Stokes方程组分源项技术问询

嘿，这个问题问到点子上了——在多组分反应流的数值模拟里，这是个很基础但绝不能搞错的关键点。我给你掰扯清楚：

首先，当你给煤油、氧气这类单独组分分别建立NS方程时，它们的连续性方程必须包含反应源项，原因如下：

• 反应本质是组分质量的转化：煤油和氧气发生燃烧反应时，煤油的质量会被消耗（源项为负），氧气的质量也会被消耗（源项为负），同时生成CO₂、H₂O等产物（这些产物组分的源项为正）。如果不在各自的连续性方程里加入这个源项ω̇_k，你的模型就完全忽略了反应对组分质量的影响，相当于假设这两种物质只是混在一起，根本没发生反应——这显然和实际情况完全背离。
• 多组分质量守恒的核心要求：单组分的连续性方程没有源项，因为不存在组分转化；但多组分体系中，每个组分的质量守恒必须考虑“反应带来的质量增减”。对应到方程上，就是每个组分的连续性方程形式为：

  ∂ρ_k/∂t + ∇·(ρ_k **v**_k) = ω̇_k
  

  其中ρ_k是组分k的密度，**v**_k是组分k的速度，ω̇_k就是反应源项（单位体积内组分k的质量生成/消耗速率）。

补充个小细节：如果把所有组分的连续性方程加起来，由于总质量守恒，所有组分的源项之和为0，所以混合物整体的连续性方程是没有源项的：

∂ρ/∂t + ∇·(ρ**v**) = 0

这里ρ是混合物的总密度，**v**是质量平均速度。

总结一下：对每个独立的反应物组分单独建模时，反应源项是连续性方程里必不可少的一部分，少了它，你的反应流模型就失去了最核心的物理基础。

内容的提问来源于stack exchange，提问作者Abdoulaye ndiongue