[论文解读] A Multiscale Method for Porous Microstructures

Donald L. Brown, Daniel Peterseim|arXiv (Cornell University)|Nov 7, 2014

Advanced Mathematical Modeling in Engineering参考文献 16被引用 11

一句话总结

本文提出了一种基于截断校正函数的局部正交分解（LOD）多尺度有限元方法，用于求解多孔微观结构中的Neumann问题。通过严格分析穿孔区域中的Poincaré常数，建立了最优收敛率，表明纤维状结构的Poincaré常数依赖于微观结构，而孤立颗粒结构则具有统一的有界性。

ABSTRACT

In this paper we develop a multiscale method to solve problems in complicated porous microstructures with Neumann boundary conditions. By using a coarse-grid quasi-interpolation operator to define a fine detail space and local orthogonal decomposition, we construct multiscale corrections to coarse-grid basis functions with microstructure. By truncating the corrector functions we are able to make a computationally efficient scheme. Error results and analysis are presented. A key component of this analysis is the investigation of the Poincaré constants in perforated domains as they may contain micro-structural information. Using a constructive method originally developed for weighted Poincaré inequalities, we are able to obtain estimates on Poincaré constants with respect to scale and separation length of the pores. Finally, two numerical examples are presented to verify our estimates.

研究动机与目标

• 解决具有精细尺度特征的复杂、高度非均质多孔微观结构中Neumann问题的挑战。
• 开发一种计算高效的多尺度格式，能够捕捉微观结构效应，而无需在全局尺度上求解。
• 分析Poincaré常数对孔隙尺寸和孔间距等微观结构参数的依赖关系。
• 在穿孔区域中，为全局支持和局部化的多尺度基函数建立误差估计。
• 提出一种构造性方法，利用加权Poincaré不等式技术估计多孔介质中的Poincaré常数。

提出的方法

• 构建一个粗粒度穿孔网格，并通过基于局部L^2投影的拟插值算子定义细尺度空间。
• 应用局部正交分解（LOD）将函数空间分解为粗尺度和细尺度分量，从而实现基函数的多尺度校正。
• 通过将校正函数截断至局部子域，确保计算效率的同时保持精度。
• 利用Galerkin正交性和衰减估计，控制校正函数截断引入的误差。
• 推导出依赖于粗网格尺寸H和截断参数m的误差界，收敛率为O(H)，且在m上呈指数衰减。
• 提出一种构造性程序，利用加权Poincaré不等式方法估计穿孔区域中的Poincaré常数。

实验结果

研究问题

• RQ1 多尺度方法如何适应处理孔隙等几何特征而非振荡系数的微观结构？
• RQ2 在穿孔区域中，Poincaré常数对孔隙尺寸和孔间距等微观结构参数的依赖关系如何？
• RQ3 局部正交分解框架能否扩展至穿孔区域中的Neumann问题，并实现严格的误差控制？
• RQ4 微观结构如何影响校正函数的衰减速率，以及由此带来的局部化格式的收敛性？
• RQ5 在何种条件下Poincaré常数相对于微观结构保持一致有界，而在何种情况下对几何特征敏感？

主要发现

• 该方法在粗网格逼近中实现了O(H)的最优收敛率，且在截断参数m上呈指数衰减。
• 穿孔区域中的Poincaré常数依赖于微观结构几何：对于极细纤维状结构可能恶化，但对于孤立颗粒则保持一致有界。
• 用于估计Poincaré常数的构造性方法可得到显式边界，其与孔间距和孔隙尺寸成比例，尤其在网状纤维微观结构中表现显著。
• 局部化多尺度格式保持了与全局版本相同的收敛率，误差受截断参数m和网格尺寸H的控制。
• 数值示例验证了理论收敛率在H、m和微观结构参数方面的有效性。
• 该方法中使用的拟插值算子具有鲁棒性和高效性，其性能在不同微观结构配置下均得到验证。