[论文解读] Theoretical derivation of Darcy's law for fluid flow in thin porous media

Francisco J. Suárez‐Grau|arXiv (Cornell University)|Jan 1, 2020

Advanced Mathematical Modeling in Engineering参考文献 22被引用 9

一句话总结

本论文在 aε ≪ ε 条件下，针对具有两个小参数 ε（膜厚）和 aε（微结构尺寸）的薄多孔介质中不可压缩牛顿流体的流动，严格推导了达西定律。通过一种新颖的双参数展开方法，推导出基于三维局部单元问题求解的二维达西定律，其渗透率张量与以往在薄介质区域中的均质化模型不同。

ABSTRACT

In this paper we study stationary incompressible Newtonian fluid flow in a thin porous media. The media under consideration is a bounded perforated $3D$ domain confined between two parallel plates. The description of the domain includes two small parameters: $\\varepsilon$ representing the distance between pates and $a_\\varepsilon$ connected to the microstructure of the domain such that $a_\\varepsilon\\ll \\varepsilon$. We consider the classical setting of perforated media, i.e. $a_\\varepsilon$-periodically distributed solid (not connected) obstacles of size $a_\\varepsilon$. The goal of this paper is to introduce a version of the unfolding method, depending on both parameters $\\varepsilon$ and $a_\\varepsilon$, and then to derive the corresponding $2D$ Darcy's law.

研究动机与目标

• 将经典的穿孔介质均质化理论扩展至具有两个小参数 ε（膜厚）和 aε（微结构尺寸）的薄多孔介质。
• 构建一种依赖于 ε 和 aε 的双参数展开方法，用于分析薄多孔域中的流体流动。
• 在 aε ≪ ε 的情形下，严格推导出二维达西定律，此时微结构远小于膜厚。
• 建立当 ε 和 aε 同时趋于零时的均质化过程收敛性，确保所推导达西定律的数学有效性。

提出的方法

• 域被建模为一个薄的三维穿孔圆柱 Qε = ω × (0, ε)，其中周期性分布着尺寸为 aε 的固体障碍物，且满足 limε→0 aε/ε = 0。
• 引入一种双参数展开算子，将周期性展开方法推广至同时处理 ε（膜厚）和 aε（微结构尺寸）的情形。
• 在穿孔域 Ωε = Qε ∩ Tε 中提出斯托克斯方程，障碍物上施加无滑移边界条件，上下板面施加无通量边界条件。
• 通过展开方法推导弱形式并取极限，得到在乘积空间 ω × Yf 上的均质化问题，其中 Yf 是单位胞腔 Y 的流体部分。
• 在三维流体域 Yf 中求解局部单元问题，以确定渗透率张量 K，其中 Kij = ∫Yf wⁱⱼ(y) dy，wⁱ 为 Yf 中斯托克斯型问题的解。
• 压力 p 满足二维达西方程：−∇x′·(K∇x′p) = f′ 在 ω 中成立，通过积分与极限识别过程推导得出。

实验结果

研究问题

• RQ1 在薄多孔介质中，流体流动的渐近行为如何依赖于膜厚 ε 与微结构尺寸 aε 的相对尺度？
• RQ2 当 aε ≪ ε 时，特别是在微结构不与膜厚方向对齐的情况下，应采用何种数学框架来推导薄多孔介质中的达西定律？
• RQ3 能否构造一种双参数展开方法，以同时处理横向（ε）与微尺度（aε）方向的均质化？
• RQ4 在薄介质区域中，有效达西定律的形式是什么？它与在均质薄介质中推导出的经典二维达西定律有何不同？
• RQ5 在 aε ≪ ε 情形下，渗透率张量 K 如何从三维局部问题中导出？其数学性质是什么？

主要发现

• 论文在 aε ≪ ε 的情形下，严格推导出薄多孔介质中的二维达西定律，此时微结构远小于膜厚。
• 渗透率张量 K 由单位胞腔 Y 的流体部分 Yf 中的三维局部单元问题确定，其中 Kij = ∫Yf wⁱⱼ(y) dy，wⁱ 为 Yf 中斯托克斯型问题的解。
• 压力 p 满足达西方程 −∇x′·(K∇x′p) = f′ 在 ω 中成立，且 K 对称正定，确保解的存在性与唯一性。
• 速度场表示为 u(x′) = K(f′ − ∇x′p)，由展开方法与弱收敛性论证得到的均质化极限推导得出。
• 证明了序列 (aε⁻²ũε, pε) 对均质化解 (u, p) 的收敛性，确立了所推导达西定律的数学有效性。
• 所推导的达西定律与以往在均质薄介质中建立的模型（如 Fabricius 等 [9]）不同，由于 aε ≪ ε 的尺度关系，其保留了三维局部问题，而其他薄介质区域模型通常采用二维局部问题。