[论文解读] Modelling and Numerical Simulation of Gas Migration in a Nuclear Waste Repository
本文提出了一种新颖的、无退化问题的组合型两相流模型，用于模拟深地质核废料处置库中氢气的迁移，采用总氢质量浓度和液相压力作为主要变量。该公式成功模拟了饱和与非饱和条件下相态的出现与消失，数值测试证实其在处理尖锐饱和度前缘和动态平衡过渡时具有鲁棒性。
We present a compositional compressible two-phase, liquid and gas, flow model for numerical simulations of hydrogen migration in deep geological radioactive waste repository. This model includes capillary effects and the gas diffusivity. The choice of the main variables in this model, Total or Dissolved Hydrogen Mass Concentration and Liquid Pressure, leads to a unique and consistent formulation of the gas phase appearance and disappearance. After introducing this model, we show computational evidences of its adequacy to simulate gas phase appearance and disappearance in different situations typical of underground radioactive waste repository.
研究动机与目标
- 解决在核废料处置库模拟中，多孔介质内气相出现与消失的建模挑战。
- 通过采用非退化变量重新表述系统，克服相饱和度方程中的数值退化问题。
- 确保在饱和与非饱和流态下，氢气迁移模拟的一致性与稳定性。
- 验证模型处理毛细效应、气体扩散和压力不连续性等复杂物理现象的能力。
- 展示该模型在典型地下放射性废料处置库中复杂、具有挑战性的实际场景中的适用性。
提出的方法
- 构建一个考虑液相与气相的组合型可压缩两相流模型，重点关注水和氢气（H₂）作为组分。
- 采用总氢质量浓度和液相压力作为主要变量，以避免饱和度方程中的退化问题。
- 引入毛细压力关系和理想气体定律描述气相行为，同时考虑液相中溶解氢的扩散。
- 应用达西-穆斯卡特定律计算相流动速率，其中相对渗透率函数依赖于相饱和度。
- 采用控制方程的弱形式以处理饱和度前缘处的不连续性和陡峭梯度。
- 通过自适应时间步长和空间离散化进行数值模拟，以解析饱和与非饱和状态之间的动态转换。
实验结果
研究问题
- RQ1如何构建一个两相流模型，使其在气相出现与消失过程中保持数值稳定性？
- RQ2何种变量选择可确保在饱和与非饱和条件下基于饱和度的方程具有一致性并避免退化？
- RQ3毛细效应与溶解氢扩散如何影响气相饱和度和压力分布的演化？
- RQ4该模型能否准确模拟从饱和状态向非饱和状态转变、伴随移动饱和度前缘的过程？
- RQ5对流与扩散通量在饱和度前缘附近氢浓度梯度的形成中起到何种作用？
主要发现
- 模型成功模拟了气相饱和度前缘在 1.92 × 10⁵ s 内从 x = 0.5 m 向左边界（Γ_in）移动的过程，前缘高度从约 0.16 降低至约 0.03。
- 在饱和区域，液相压力迅速上升（从约 1 MPa 上升至约 1.6 MPa，t = 10³ s 时），随后缓慢下降，表明上方气相对流体的动态压缩作用。
- 在饱和度前缘处形成了强烈的溶解氢浓度梯度，持续至 t = 1.92 × 10⁵ s，这是由于从左向右的对流向运与相反方向的扩散通量相互竞争所致。
- 在 t = 1.92 × 10⁵ s 之后，整个计算域变为非饱和（S_g > 0），所有变量（压力、氢气密度、饱和度）单调演化至空间均匀的稳定状态。
- 系统达到平衡状态，压力与浓度分布连续，与封闭域中的质量守恒一致。
- 模型准确捕捉了压力连续性的即时恢复以及浓度梯度的渐进平滑过程，验证了其在处理相变动力学方面的鲁棒性。
更好的研究，从现在开始
从阅读论文到最终审阅，大幅缩短您的研究时间。
无需绑定信用卡
本解读由 AI 生成，并经人工编辑审核。