[论文解读] Localized bases for finite dimensional homogenization approximations with non-separated scales and high-contrast

Houman Owhadi, Lei Zhang|arXiv (Cornell University)|Nov 3, 2010

Advanced Mathematical Modeling in Engineering参考文献 35被引用 25

一句话总结

该论文提出了一种局部化、与网格无关的方法，用于构造具有非分离尺度和高对比度系数的椭圆型、抛物型及双曲型PDE解的有限维近似空间。通过利用当源项属于 $ L^2 $ 时解空间在 $ H^1 $ 中的紧嵌入性质，该方法通过求解局部PDE来构造局部化基函数，实现了 $ \tilde{O}(h^{2-2\beta}) $ 的精度，且近似空间维数为 $ O(h^{-d}) $，其中基函数在大小为 $ O(h^{\beta} \ln(1/h)) $ 的子域上具有紧支集，且 $ \beta \ne 1 $。当包含缓冲区时，该方法在高对比度介质中仍能保持精度。

ABSTRACT

We construct finite-dimensional approximations of solution spaces of divergence form operators with $L^\infty$-coefficients. Our method does not rely on concepts of ergodicity or scale-separation, but on the property that the solution space of these operators is compactly embedded in $H^1$ if source terms are in the unit ball of $L^2$ instead of the unit ball of $H^{-1}$. Approximation spaces are generated by solving elliptic PDEs on localized sub-domains with source terms corresponding to approximation bases for $H^2$. The $H^1$-error estimates show that $\mathcal{O}(h^{-d})$-dimensional spaces with basis elements localized to sub-domains of diameter $\mathcal{O}(h^α\ln \frac{1}{h})$ (with $α\in [1/2,1)$) result in an $\mathcal{O}(h^{2-2α})$ accuracy for elliptic, parabolic and hyperbolic problems. For high-contrast media, the accuracy of the method is preserved provided that localized sub-domains contain buffer zones of width $\mathcal{O}(h^α\ln \frac{1}{h})$ where the contrast of the medium remains bounded. The proposed method can naturally be generalized to vectorial equations (such as elasto-dynamics).

研究动机与目标

• 开发一种适用于具有 $ L^∞ $ 系数的散度型PDE的有限维近似方法，且不依赖于尺度分离或遍历性。
• 解决在经典均质化方法失效的高对比度与多尺度介质中逼近解空间的挑战。
• 构造局部化、计算高效的基函数，以实现精确的 $ H^1 $-范数近似。
• 通过相同的局部化基框架，将该方法推广至抛物型与双曲型问题。
• 建立与系数对比度或系数矩阵最小特征值无关的严格 $ H^1 $-误差估计。

提出的方法

• 利用当源项 $ g ∈ L^2(Ω) $ 时解空间在 $ H^1_0(Ω) $ 中的紧嵌入性质，而非 $ H^{-1}(Ω) $，从而实现有限维近似。
• 通过在子域上求解具有 $ g $-类源项的局部椭圆型PDE来构造近似基函数，确保基函数具有局部支集。
• 使用通量范数 $ \|\cdot\|_{a\text{-flux}} $，其与能量范数等价且与 $ a $ 无关，以推导误差估计。
• 应用文献[10]中的转移性质，以控制局部解在其支集外的衰减，从而实现具有指数衰减的局部化。
• 在高对比度介质中引入宽度为 $ O(h^\beta \ln(1/h)) $ 的缓冲区，以保持精度。
• 通过相同的局部化基构造方法，将该框架推广至矢量型方程（如弹性动力学）。

实验结果

研究问题

• RQ1 能否在不假设尺度分离或遍历性的前提下，为多尺度PDE构造有限维近似空间？
• RQ2 如何构造局部化基函数，使其 $ H^1 $-误差随网格尺寸 $ h $ 衰减，即使在高对比度介质中亦成立？
• RQ3 为实现给定的 $ H^1 $-范数精度，局部化子域的最优尺寸是什么？
• RQ4 相同的局部化基框架能否推广至如抛物型与双曲型方程等时间依赖问题？
• RQ5 系数 $ a $ 中的高对比度如何影响局部近似的精度，又该如何缓解？

主要发现

• 当源项属于 $ L^2(Ω) $ 时，PDE的解空间在 $ H^1_0(Ω) $ 中是紧嵌入的，从而可在 $ H^1 $-范数下实现任意精度的有限维近似。
• 维数为 $ O(h^{-d}) $ 的近似空间，其基函数局部化于直径为 $ O(h^\beta \ln(1/h)) $ 的子域上，且 $ \beta \in [1/2, 1) $，可在 $ H^1 $-范数下实现 $ O(h^{2-2\beta}) $ 的精度。
• 在高对比度介质中，只要局部子域包含宽度为 $ O(h^\beta \ln(1/h)) $ 的缓冲区（其中对比度保持有界），该方法仍能保持 $ O(h^{2-2\beta}) $ 的精度。
• 通量范数误差估计与 $ \lambda_{\min}(a) $ 和 $ \lambda_{\max}(a) $ 无关，从而在广泛系数范围内具有鲁棒性。
• 该方法可自然推广至矢量型方程（如弹性动力学），并保持相同的收敛速率。
• 在 $ g \in H^{-\nu} $，$ \nu < 1 $ 的强紧性条件下，为经典均质化、数值均质化与模型降阶方法提供了统一的理论基础。