[论文解读] Travelling wave solutions for gravity fingering in porous media flows

Koondanibha Mitra, Andreas Rätz+1 more|arXiv (Cornell University)|Jan 1, 2020

Fluid Dynamics and Thin Films参考文献 29被引用 1

一句话总结

本文通过采用具有动态毛细压力的滞后模型，研究了非饱和多孔介质中重力指流的行波解。通过将时变问题转化为移动参考系，作者建立了单个指流以恒定速度传播的自由边界问题，利用变分法证明了解的存在性，并通过数值结果展示了与时变模拟和实验一致的指流形状与速度。

ABSTRACT

We study an imbibition problem for porous media. When a wetted layer is above a dry medium, gravity leads to the propagation of the water downwards into the medium. In experiments, the occurrence of fingers was observed, a phenomenon that can be described with models that include hysteresis. In the present paper, we describe a single finger in a moving frame and set up a free boundary problem to describe the shape and the motion of one finger that propagates with a constant speed. We show the existence of solutions to the travelling wave problem and investigate the system numerically.

研究动机与目标

• 建立重力驱动的非饱和多孔介质中非标准 Richards 方程无法捕捉指流流动模式的模型。
• 通过引入滞后与动态毛细压力效应（特别是 p = pc(s) + τ∂ts 的 τ 修正模型）来克服标准模型的局限性。
• 分析系统在移动参考系中的行波形式，以描述单个指流的形状与传播速度。
• 在具有预设通量与边界条件的半无限域上，建立相应自由边界问题解的存在性。
• 通过数值方法计算指流剖面与传播速度，并与时变模拟及实验观测结果进行验证。
• 采用符号函数的正则化方法，推导出 sH 与 ∇pH 的 Lipschitz 有界性，确保在 H → ∞ 极限下的稳定性与收敛性。

提出的方法

• 通过坐标变换 (y, z + ct) 将时变的入渗问题转化为行波框架，将原始 PDE 系统简化为移动参考系中的稳态问题。
• 将系统表述为耦合 PDE 系统：c∂zs = ∇⋅(k(s)[∇p + gez]) 与 cτ∂zs = [p − pc(s)]+，分别表示质量守恒与动态毛细滞后效应。
• 利用变分原理，在截断域 ΩH = (0,L)×(0,H) 上证明解的存在性，通过最小化包含 k(s) 与压力梯度的能量泛函。
• 施加边界条件：在 z = 0 处 s = s0 与 p = p0（底部），在 z → ∞ 处通量为 F∞，以及侧向无通量条件。
• 应用数值延拓与有限元方法求解截断问题，提取指流剖面与传播速度。
• 通过符号函数的正则化方法，推导出 sH 与 ∇pH 的 Lipschitz 有界性，确保在 H → ∞ 极限下的稳定性与收敛性。

实验结果

研究问题

• RQ1 能否为具有滞后与动态毛细压力的多孔介质中单个重力指流构造出行波解？
• RQ2 在移动参考系中，此类解的存在性需满足何种条件？其与通量、速度及边界数据的关系如何？
• RQ3 截断问题的数值解在多大程度上逼近完整的无限域问题？在 H → ∞ 极限下，压力与饱和度的行为如何？
• RQ4 计算得到的指流形状与传播速度是否与时间依赖模拟及实验观测结果一致？
• RQ5 动态毛细参数 τ 在决定指流剖面与稳定性方面起到何种作用？

主要发现

• 通过变分公式在截断域 ΩH 上证明了行波问题解的存在性，最小化相应能量泛函。
• 仅当在固定高度 z0 上对 k(sH) 的积分远离零时，压力在 H → ∞ 极限下才保持有界；否则压力将无界。
• 数值模拟生成的饱和度剖面呈指流形状，与时间依赖模拟及实验观测结果一致，验证了行波方法的有效性。
• 传播速度 c 与时间依赖计算所得结果高度一致，尽管不同数值格式间存在微小差异。
• 解在空间与时间上表现出 Lipschitz 连续性，|∂zsH| 与 |∂ysH| 的有界性依赖于 c、τ 与初始数据 s0、p0。
• 分析表明，解的非唯一性可能是导致计算 c 值出现差异的原因，提示在相同边界条件下可能存在多个可能的指流速度。