我手头有一个 CFD 求解器，它可以处理两个阶段，到目前为止，它能够在两个阶段相遇的界面上交换动量和能量。这是如何完成的，在 Hu、Khoo 和 Adams 的以下论文“可压缩流的保守接口方法”中有详细记录，该论文可公开获得。

在本文中，等式 (10) 给出了如何更新所谓的“切割单元”中的保守变量的一般概念，这些单元由界面切割（通过水平集方法跟踪）。

在这个更新程序持有平均能量和动量交换。现在我想扩展代码，使其可以处理反应前沿（如燃烧或空化）。$\hat{X}$

出现的黎曼问题现在由反应黎曼求解器解决，该求解器提供接触不连续处的压力和速度，以及通过反应前沿的质量流量。

到目前为止一切都很好，但现在我很难调整以使其将质量（密度）传输到前面。$\hat{X}$

如果有任何想法或提示，我将不胜感激。