[论文解读] Analysis of the roughness regimes for micropolar fluids via homogenization

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

Advanced Mathematical Modeling in Engineering参考文献 32被引用 7

一句话总结

本文利用多尺度渐近分析方法，研究了具有周期性粗糙边界的薄域中微极流体的流动行为。基于比值 λ = ηε/ε，识别出三种不同的粗糙度 regime——Stokes、Reynolds 和高频——并推导出包含微极效应的广义 Reynolds 方程，其流动因子具有 regime 依赖性。主要贡献在于建立了一个统一的微极润滑渐近框架，其中流动因子通过求解三维微极 Stokes 问题（Stokes regime）、二维微极 Reynolds 问题（Reynolds regime）或直接恢复经典微极 Reynolds 方程（高频 regime）获得。

ABSTRACT

We study the asymptotic behavior of micropolar fluid flows in a thin domain of thickness $\eta_\varepsilon$ with a periodic oscillating boundary with wavelength $\varepsilon$. We consider the limit when $\varepsilon$ tends to zero and, depending on the limit of the ratio of $\eta_\varepsilon/\varepsilon$, we prove the existence of three different regimes. In each regime, we derive a generalized Reynolds equation taking into account the microstructure of the roughness.

研究动机与目标

• 将经典 Newton 流体的多尺度渐近框架扩展至具有周期性粗糙边界的薄域中的微极流体。
• 基于域厚度与粗糙度波长之比 λ = ηε/ε，对微极流体流动的渐近行为进行分类。
• 在每个 regime 中，推导出考虑微极微结构的广义 Reynolds 型方程。
• 建立系统化方法，用于计算极限方程中与 regime 相关的流动因子 (Aλ, bλ)。
• 证明在高频极限下，经典微极 Reynolds 方程得以恢复，原因在于振荡区域中速度与微旋转趋于零。

提出的方法

• 在厚度为 ηε、上边界为 hε(x′) = ηε h(x′/ε) 的薄域 Ωε 中建立三维微极 Stokes 方程。
• 采用尺度变换技术，将问题映射至具有周期性微结构的固定参考元胞 Y = Y′ × (0, h(y′)) 上。
• 应用改进的展开方法与薄域扩展结果，以处理边界层效应。
• 结合两尺度收敛与弱形式，令 ε → 0 时取极限，并通过 λ = lim ηε/ε 区分三种 regime。
• 求解局部单元问题：在 0 < λ < ∞ 时为三维微极 Stokes 问题，在 λ = 0 时为二维微极 Reynolds 问题，在 λ = ∞ 时无局部问题。
• 推导宏观广义 Reynolds 方程：div(−Aλ∇p + bλ) = 0，其中 Aλ 与 bλ 由局部解计算得出。

实验结果

研究问题

• RQ1 微极流体在薄且粗糙域中的渐近行为如何依赖于厚度与粗糙度波长之比 λ = ηε/ε？
• RQ2 在三种粗糙度 regime（Stokes、Reynolds、高频）中，控制微极流动的广义 Reynolds 方程具有何种形式？
• RQ3 在每个 regime 中，宏观流动因子 Aλ 与 bλ 如何计算？它们如何反映微极微结构的影响？
• RQ4 经典微极 Reynolds 方程是否在高频极限下出现？若出现，其条件为何？
• RQ5 在不同 regime 中，速度场与微旋转场之间的耦合关系是持续存在还是被破坏？

主要发现

• 微极流体存在三种截然不同的粗糙度 regime：Stokes（0 < λ < ∞）、Reynolds（λ = 0）与高频（λ = ∞），其分类与 Newton 流体类似，但已推广至微极情形。
• 在 Stokes regime（0 < λ < ∞）中，流动因子 Aλ 与 bλ 通过求解三维局部微极 Stokes 类问题获得，且保持速度与微旋转之间的耦合关系。
• 在 Reynolds regime（λ = 0）中，流动因子由二维局部微极 Reynolds 类问题求得，计算显著简化，同时仍保留微结构耦合特性。
• 在高频 regime（λ = ∞）中，振荡区域内的速度与微旋转趋于零，导致非振荡区域中恢复经典微极 Reynolds 方程，且无需求解局部问题。
• 宏观方程形式为 div(−Aλ∇p + bλ) = 0，其中 Aλ 与 bλ 依赖于 λ，表征了有效的微极行为。
• 在高频极限下，经典微极 Reynolds 方程（1.3）被恢复，验证了在极端粗糙条件下与已有结果的一致性。