[论文解读] Elements of mathematical foundations for a numerical approach for weakly random homogenization problems

Arnaud Anantharaman, Claude Le Bris|arXiv (Cornell University)|May 21, 2010

Advanced Mathematical Modeling in Engineering参考文献 1被引用 10

一句话总结

本文通过推导同胚张量在扰动参数 η 的二阶展开，为弱随机均质化问题构建了一个严格的微扰框架。利用泛函分析与形式随机展开，证明了一阶和二阶校正项的收敛性，并通过数值测试展示了其计算效率与准确性。

ABSTRACT

This work is a follow-up to our previous work "A numerical approach related to defect-type theories for some weakly random problems in homogenization" (preprint available on this archive). It extends and complements, both theoretically and experimentally, the results presented there. Under consideration is the homogenization of a model of a weakly random heterogeneous material. The material consists of a reference periodic material randomly perturbed by another periodic material, so that its homogenized behavior is close to that of the reference material. We consider laws for the random perturbations more general than in our previous work cited above. We prove the validity of an asymptotic expansion in a certain class of settings. We also extend the formal approach introduced in our former work. Our perturbative approach shares common features with a defect-type theory of solid state physics. The computational efficiency of the approach is demonstrated.

研究动机与目标

• 将弱随机非均质材料的微扰方法从伯努利分布扰动推广至更一般的随机分布。
• 通过勒贝格空间中的弱收敛，严格证明同胚张量的二阶渐近展开。
• 基于扰动场的分布律，形式化一种随机展开方法，实现校正项的确定性计算。
• 与直接随机均质化相比，展示微扰方法在计算效率与鲁棒性方面的优势。
• 在特定情形（特别是在一维情况下）建立形式推导与严格推导之间的等价性。

提出的方法

• 在遍历性与平稳性假设下，利用 $ L^\nu(\text{D}, L^p(\text{D})) $ 中的弱收敛（$ p > 1 $），推导出同胚张量的一阶渐近展开 $ A_{\text{per}}^* + \theta \bar{A}_1^* + o(\theta) $。
• 通过推论 3 和 4，在附加假设下建立二阶收敛性，依赖于泛函分析技术。
• 提出一种基于随机场 $ b_\theta $ 分布律的正式方法，假设其在 $ \theta $ 上存在分布展开，从而能够通过超胞元上的确定性单元问题极限计算 $ \bar{A}_1^* $ 与 $ \bar{A}_2^* $。
• 通过显式构造与弱极限，证明当 $ N \to \bigotimes $ 时，一阶校正 $ A_{1}^{*,N} $ 的收敛性。
• 将二阶校正 $ A_{2}^{*,N} $ 定义为提取子列的极限，并通过技术性分析建立其有界性（详见 [1]）。
• 将该正式方法与固态物理中的缺陷理论相联系，强调罕见但具有影响力的扰动在改变有效材料性能中的作用。

实验结果

研究问题

• RQ1 在一般随机分布下，且具备适当的收敛性质时，弱随机均质化问题的微扰方法是否可严格证明，超越伯努利分布的情形？
• RQ2 在同胚张量渐近展开中，一阶与二阶校正项 $ \bar{A}_1^* $ 与 $ \bar{A}_2^* $ 的收敛性需满足何种条件？
• RQ3 形式随机展开方法与严格的泛函分析方法在推导校正项时有何异同？
• RQ4 通过基于确定性单元问题的微扰展开，能否显著降低随机均质化的计算成本？
• RQ5 在一维情形下，能否通过显式计算使形式展开完全严格化，且其结果是否与严格推导一致？

主要发现

• 一阶校正 $ \bar{A}_1^* $ 被严格证明存在，并可通过 $ L^\nu(\text{D}, L^p(\text{D})) $ 中的弱极限计算，将先前结果推广至一般随机分布。
• 二阶校正 $ \bar{A}_2^* $ 在 $ N \to \infty $ 时通过提取子列收敛，且在一维情形下通过显式计算得到验证。
• 在一维设置下，基于 $ b_\theta $ 分布律的形式展开与泛函分析所得的严格结果完全一致，验证了形式方法的有效性。
• 二阶校正 $ \bar{A}_2^* $ 显式表达为三项之和，分别涉及 $ \mathbb{E}(\bar{B}_0) $、$ \text{Var}(\bar{B}_0) $ 与 $ \mathbb{E}(\bar{B}_0)^2 $，其积分涉及周期胞元中的校正函数与辅助函数。
• 数值测试表明，即使在中等强度扰动下，该微扰方法仍具有计算高效与高精度的特性，优于直接随机均质化方法。
• 该形式方法产生仅依赖于周期结构与扰动二阶矩的确定性、可计算校正项，而无需完整路径实现。