[论文解读] Weak solution of the Hele-Shaw problem: shocks and viscous fingering

Seung-Yeop Lee, Razvan Teodorescu+1 more|arXiv (Cornell University)|Dec 2, 2008

Cold Atom Physics and Bose-Einstein Condensates参考文献 1被引用 3

一句话总结

本文提出了一种弱解来解决Hele-Shaw问题，通过黏性激波（viscous shocks）消除有限时间内的尖点奇点——即压力和涡度中的不连续性，形成分支且不断增长的树状结构。该解在宏观尺度上保持可积性与无旋流动，同时允许微观尺度的弛豫，其中自相似的(2,3)-尖点奇点触发一个普适的、尺度不变的激波分支事件，由超越方程和唯一的过渡比η ≈ 0.915决定。

ABSTRACT

In Hele-Shaw flows, boundaries between fluids develop unstable viscous fingers. At vanishing surface tension, the fingers further evolve to cusp-like singularities. We show that the problem admits a {\it weak solution} where shock fronts triggered by a singularity propagate together with a fluid. Shocks form a growing, branching tree of a mass deficit, and a line distribution of vorticity where pressure and velocity of the fluid have finite discontinuities. Imposing that the flow remain curl-free at macroscale determines the shock graph structure. We present a self-similar solution describing shocks emerging from a generic (2,3)-cusp singularity -- an elementary branching event.

研究动机与目标

• 解决在表面张力趋近于零时，黏性指进（viscous fingering）中出现的有限时间尖点奇点问题。
• 建立一个弱解框架，保持宏观尺度上的不可压缩性和无旋流动，同时允许物理约束在微观尺度上的弛豫。
• 表征黏性激波作为压力和涡度不连续性的形成机制，其结构为不断增长的分支树。
• 识别从(2,3)-尖点奇点出发的基元分支事件的普适自相似解。
• 证明正则化过程保持Hele-Shaw流动的可积性，从而将其与流体动力学中其他类型的激波区分开来。

提出的方法

• 使用复势理论和解析高度函数Y(X,t)来表述Hele-Shaw问题，将Darcy定律表示为复势φ = ψ + ip的形式。
• 通过泛函Ω(X) = −i∫ρ₀Y dX引入Darcy定律的积分形式，将时间演化与复势通过ẌΩ = iφ联系起来。
• 将激波识别为Ω的反Stokes线，即满足Re Ω = 0的线，并通过可接受性条件Ω > 0在激波边界上侧来选择。
• 利用由椭圆函数（Weierstrass ℘-函数）导出的超越方程来描述激波的几何形状和端点运动。
• 应用Rankine-Hugoniot条件计算激波速度，表明激波传播速度超过流体吸收速率（t⁵/⁴ vs. t）。
• 计算普适的过渡比η = limₜ→₀ ẇCₜ>₀ / ẇCₜ<₀ ≈ 0.915，表征奇点处流量容量的突变。

实验结果

研究问题

• RQ1 在缺乏表面张力的条件下，如何对Hele-Shaw问题中的有限时间尖点奇点进行正则化？
• RQ2 何种物理机制使得流动在越过奇点后仍能保持宏观尺度上的不可压缩性和无旋条件？
• RQ3 由(2,3)-尖点奇点引发的黏性激波网络的普适结构是什么？
• RQ4 激波结构如何与复势及积分Darcy定律的反Stokes线相关联？
• RQ5 容量过渡比η在表征激波特形成过程中的作用是什么？

主要发现

• Hele-Shaw问题存在一种弱解，其中尖点奇点会触发黏性激波的形成——即有限压力不连续性和涡度分布的线状结构。
• 激波形成一个不断增长的分支树结构，其传播速度超过流体排出速度，激波速度满足V⊥ = X⊥ / t的标度关系。
• 激波图由宏观尺度上保持无旋流动的条件唯一确定，从而确保物理一致性。
• 从(2,3)-尖点出发的基元分支事件具有自相似性且与流动细节无关，由涉及Weierstrass椭圆函数的超越方程描述。
• 流量容量的转变是突变的，由一个普适比值η ≈ 0.91522表征，该值由临界时刻的g₃/g₂比值导出。
• 积聚在激波上的质量亏损可量化为∫ₑ₃ᵉ₁₂ σ dℓ = (4/5)|ψ(e₁,₂)| ≈ (4/5)(6.34513) t⁵/⁴，表明激波导致指端流体后退。