Abstract Multiphase flow equations for two components (for instance water and hydrogen) and two phases (liquid and gas), with equilibrium phase exchange, have been used to simulate the process of miscible displacement in porous media. Laboratory and field studies have shown that this assumption fails under certain circumstances especially for accurate description of the pollution and a finite transfer velocity, called here the kinetics. We propose a numerical scheme based on a two-step convection/diffusion-relaxation strategy to simulate the non-equilibrium model. In the first step we solve the intra-phase transfer (convection/diffusion) working with liquid saturation, liquid pressure and dissolved hydrogen concentration as primary variables. In the second step, we address the interphase transfer using the solubility relation that is solved by projection on the equilibrium state. This technique also ensures the positivity for the liquid saturation and produces energy estimates. One important advantage of this approach is the fact that the simulation can be easily adapted to different linear and non-linear equilibrium laws between phases. Another advantage is that our proposed method includes the kinetics in a projection step and keeps the displacement of the components. We implemented this new model in in-house simulation code and present numerical results comparing the model with equilibrium and non-equilibrium one.

中文翻译：

---

具有动力学传质的部分混相两相流的建模和模拟

摘要 具有平衡相交换的两组分（例如水和氢）和两相（液体和气体）的多相流动方程已被用于模拟多孔介质中的混相置换过程。实验室和现场研究表明，这种假设在某些情况下是不成立的，特别是对于污染和有限传输速度的准确描述，这里称为动力学。我们提出了一种基于两步对流/扩散-弛豫策略的数值方案来模拟非平衡模型。在第一步中，我们以液体饱和度、液体压力和溶解氢浓度作为主要变量来求解相内转移（对流/扩散）。第二步，我们使用通过平衡态投影求解的溶解度关系来解决相间转移。该技术还确保液体饱和度的正性并产生能量估计。这种方法的一个重要优点是模拟可以很容易地适应相之间不同的线性和非线性平衡定律。另一个优点是我们提出的方法包括投影步骤中的动力学并保持组件的位移。我们在内部模拟代码中实现了这个新模型，并提供了将模型与平衡和非平衡模型进行比较的数值结果。这种方法的一个重要优点是模拟可以很容易地适应相之间不同的线性和非线性平衡定律。另一个优点是我们提出的方法包括投影步骤中的动力学并保持组件的位移。我们在内部模拟代码中实现了这个新模型，并提供了将模型与平衡和非平衡模型进行比较的数值结果。这种方法的一个重要优点是模拟可以很容易地适应相之间不同的线性和非线性平衡定律。另一个优点是我们提出的方法包括投影步骤中的动力学并保持组件的位移。我们在内部模拟代码中实现了这个新模型，并提供了将模型与平衡和非平衡模型进行比较的数值结果。