[论文解读] A homogenized limit for the 2D Euler equations in a perforated domain

Matthieu Hillairet, Christophe Lacave+1 more|arXiv (Cornell University)|Jul 9, 2019

Advanced Mathematical Modeling in Engineering参考文献 28被引用 2

一句话总结

本文在临界孔隙率 regime 下首次建立了二维欧拉方程在穿孔区域中的均匀化极限，其中孔径大小与孔间距之比趋于一个正常数，导致多孔介质具有非零但较小的体积分数。通过反射法，推导了散度-旋度问题的精细估计，得到一个带有有效矩阵项的均匀化欧拉型系统，该系统描述了涡度与速度之间的耦合，适用于先前稀疏或密集极限未覆盖的小但非零孔隙率情形。

ABSTRACT

We study the motion of an ideal incompressible fluid in a perforated domain. The porous medium is composed of inclusions of size $a$ separated by distances $ ilde d$ and the fluid fills the exterior. We analyse the asymptotic behavior of the fluid when $(a, ilde d) o (0,0)$. If the inclusions are distributed on the unit square, this issue is studied recently when $\frac{ ilde d}a$ tends to zero or infinity, leaving aside the critical case where the volume fraction of the porous medium is below its possible maximal value but non-zero. In this paper, we provide the first result in this regime. In contrast with former results, we obtain an Euler type equation where a homogenized term appears in the elliptic problem relating the velocity and the vorticity. Our analysis is based on the so-called method of reflections whose convergence provides novel estimates on the solutions to the div-curl problem which is involved in the 2D-Euler equations.

研究动机与目标

• 分析当孔径大小 a 和孔间距 ˜d 同时趋于零且满足 ˜d/a → ¯k > 0 时，二维不可压缩欧拉流在穿孔区域中的渐近行为。
• 填补先前研究仅考虑稀疏（˜d/a → ∞）和密集（˜d/a → 0）情形所留下的理论空白。
• 在多孔介质具有小但非零体积分数的临界 regime 下，推导欧拉方程的均匀化系统。
• 利用反射法建立穿孔区域中散度-旋度问题的新估计，从而实现向均匀化极限的收敛。

提出的方法

• 分析分为两步：首先，为穿孔区域 FN 中的散度-旋度问题导出关于速度与通过流函数 ψN 关联的涡度的精细一致估计。
• 采用反射法构造近似解，并基于均匀化解 ψc 推导流函数 ψN 的误差界。
• 均匀化问题表述为：在 R² 中，有 div[(I₂ + kMK)∇ψc] = f，其中 MK 是依赖于孔洞形状 K 的矩阵，k 为极限体积分数。
• 关键技术创新在于：利用加权范数和涡度梯度的对数型索伯列夫型界，推导出真实速度 uN 与均匀化速度 uc 之间差值的各向异性 L∞ 估计。
• 通过将能量估计和 Gronwall 型论证应用于涡度输运方程，建立涡度与速度场的收敛性。
• 采用Bootstrap方法统一延长解的存在时间，确保收敛性在固定时间 T 内成立。

实验结果

研究问题

• RQ1 当多孔介质具有小但非零体积分数时，即在 ˜d/a → ¯k > 0 的临界 regime 下，二维欧拉方程的均匀化极限是什么？
• RQ2 在欧拉框架下，非零多孔介质的存在如何影响宏观流体动力学，特别是对有效粘性或各向异性的影响？
• RQ3 反射法是否能提供足够估计，以控制穿孔区域中真实解与均匀化解在临界 regime 下的误差？
• RQ4 控制流体运动的均匀化方程的结构是什么？它与经典欧拉方程或达西型定律有何不同？
• RQ5 穿孔区域中的涡度与速度场如何收敛到其均匀化对应量？该收敛过程的统一有界性如何？

主要发现

• 在临界 regime 下，二维欧拉方程的均匀化极限由一个修正的椭圆问题控制：在 R² 中，有 div[(I₂ + kMK)∇ψc] = f，其中 MK 是依赖于孔洞形状 K 的矩阵。
• 对于单位圆盘，MK = 2I₂，因此有效算子变为 div[3I₂∇ψc] = f，表明对经典拉普拉斯算子存在显著的各向异性修正。
• 速度场 uN 到均匀化场 uc 的收敛性在 L∞(OT) 中建立，收敛速率由 F(N,k) = (a/d)^{3−η} + ||µ−k||_{W^{−1,p}}^{(1−η)(p+2)/p} + ||µ−k||_{W^{−1,p}}^{1/2} + ||k||_{L∞}^2 控制。
• 通过对数型 Gronwall 估计，涡度梯度在时间上保持一致有界，这对控制欧拉方程中的非线性项至关重要。
• 流体粒子轨迹 XN(t,x) 在穿孔区域中以误差 O(F(N,k)) 均匀收敛于均匀化轨迹 Xc(t,x)，在区间 [0,T] 上成立。
• 该结果对具有紧支集且 L∞ 范数有界的初始涡度成立，且在固定时间 T 内收敛性一致成立，前提是体积分数 k 较小且属于 FV(ε₀)。