单次海恩斯跳跃中前缘回缩的线性响应理论

《Journal of Fluid Mechanics》：Linear response theory of front back-contraction during a single Haines jump

　　

研究人员研究了局部瞬态扰动下多孔介质中流体-流体界面的线性响应。一个单孔因微小位移而偏离力学平衡，引发界面的回缩（back-contraction）。在假设总流量可忽略的前提下，界面以随时间与空间变化的剖面 h(x,t) 松弛，针对平面几何建模。该设置反映了利用孔喉尺寸局部非均质性在实验上可实现的条件下的情况。研究人员通过解析计算与孔隙网络模拟（pore-network simulations），确定了弛豫过程中演化的界面形状，以及相关的压力场与速度场。

研究背景与问题提出

多孔介质内流体运动研究长期受关注，不混溶流体驱替可产生从致密到分形多种形态，主要由黏性、毛细（capillary）与重力力平衡主导，润湿特性及局部孔隙几何亦起关键作用。慢速排液（drainage）中以间歇性Haines跳跃为特征，由界面毛细力欠压驱动：前缘某孔因暴露于更低毛细压力阈值的邻孔而失稳，流体自失稳区被抽向前缘直至重新平衡。已有改进入侵渗流、孔隙网络模型（PNMs）、密度泛函水动力学与格子玻尔兹曼等数值研究刻画了Haines跳跃尺寸分布与动态，但针对单次跳跃详细动态的解析理论仍缺失。刻画前缘回缩（front back-contraction）的空间与时间尺度，对建立两相流连续介质描述及CO₂封存等应用至关重要。本文目标即给出单次孔失稳后界面弛豫剖面 h(x,t) 及压力、速度场的解析理论，并在二维与三维下与数值对照。

主要技术方法

研究人员采用线性响应框架：将两流体Darcy定律（忽略重力与惯性）与不可压缩条件结合，对压力扰动 p_i(x,z,t) 取Fourier变换得 Laplace 方程解；由Young–Laplace关系将界面位移 h 与毛细压力变化 Δp_c 线性耦合，定义 Γ=γ/(r₀²α)；通过界面法向速度连续条件引入阻力项 R=(R₁+R₂)/2，其中 R_i=μ_i?(0)/κ，导出傅里叶空间一阶常微分方程 2R?(q,t)+|q|Γh(q,t)=0，以脉冲初值 h(x,0)=?ΔAδ(x) 解得指数衰减核，反演得 h(x,t) 闭式。数值侧采用孔隙网络模拟（pore-network simulations）在二维与三维平面几何中复现局部大孔喉设定，校验解析剖面与压力场。

研究结果

1. 2.
   理论（Theory）
   2.1 表面动力学（Surface dynamic）：研究人员推导界面扰动远场压力衰减 p₁∝e^|q|z、p₂∝e^?|q|z；由曲率半径变化 Δr 得 h=?αΔr 及 Δp_c=γh/(r₀²α)=Γh；在 |h|?a 下略去高阶项，由 z 向Darcy通量 ?(0)?=?(κ/μ_i)?_zp_i|_z=0 得傅里叶空间式，相加消去背景压得弛豫方程（2.31），解为 h(q,t)=?ΔAe^?|q|Γt/(2R)，反演积分（2.34）给出时空剖面。结论：单次扰动下界面以特征时间 τ(q)=2R/(|q|Γ) 按波数依赖指数弛豫，总流量近似为零。
2. 3.
   数值与对比（§3 概述）
   研究人员将解析 h(x,t) 与孔隙网络模拟在二维/三维、等孔喉加局部大孔设定下比对，界面形状、压力场及速度场均高度吻合，验证了线性响应假定与 Γ、R 参数化在跨孔尺度后的适用性。

讨论与结论翻译

研究人员指出，该线性响应理论首次为单个Haines跳跃的前缘回缩提供了可解析追踪的连续介质描述，明确了弛豫长度与时间尺度由 Γ（界面刚度）与 R（两相黏度-孔隙度阻力）共同控制；该框架可嵌入更大规模排液过程，为两相流宏观方程提供微观基准，并对X射线CT观测的3D单次跳跃分析、CO₂封存界面动力学建模具支撑意义。结论：在总流可忽略、扰动小于孔尺寸 a 的限定下，局部孔失稳引发的界面回缩可由 h(x,t) 指数傅里叶谱完全描述，解析与孔隙网络模拟一致。