[论文解读] Homogenization of elliptic problems: error estimates in dependence of the spectral parameter

T. A. Suslina|arXiv (Cornell University)|Jun 29, 2014

Advanced Mathematical Modeling in Engineering参考文献 12被引用 5

一句话总结

本文为具有周期系数的强椭圆微分算子的同伦化建立了精确的算子误差估计，重点关注有界域中狄利克雷与诺伊曼问题的预解。推导了在算子范数下依赖于缩放参数 $\varepsilon$ 和谱参数 $\zeta$ 的收敛速率，实现了对预解逼近在 $H^1$-范数下的 $O(\varepsilon)$ 误差界，常数由谱和几何参数控制。

ABSTRACT

We consider a strongly elliptic differential expression of the form $b(D)^* g(x/\varepsilon) b(D)$, $\varepsilon >0$, where $g(x)$ is a matrix-valued function in ${\mathbb R}^d$ assumed to be bounded, positive definite and periodic with respect to some lattice; $b(D)=\sum_{l=1}^d b_l D_l$ is the first order differential operator with constant coefficients. The symbol $b(ξ)$ is subject to some condition ensuring strong ellipticity. The operator given by $b(D)^* g(x/\varepsilon) b(D)$ in $L_2({\mathbb R}^d;{\mathbb C}^n)$ is denoted by $A_\varepsilon$. Let ${\mathcal O} \subset {\mathbb R}^d$ be a bounded domain of class $C^{1,1}$. In $L_2({\mathcal O};{\mathbb C}^n)$, we consider the operators $A_{D,\varepsilon}$ and $A_{N,\varepsilon}$ given by $b(D)^* g(x/\varepsilon) b(D)$ with the Dirichlet or Neumann boundary conditions, respectively. For the resolvents of the operators $A_\varepsilon$, $A_{D,\varepsilon}$, and $A_{N,\varepsilon}$ in a regular point $ζ$ we find approximations in different operator norms with error estimates depending on $\varepsilon$ and the spectral parameter $ζ$.

研究动机与目标

• 推导具有狄利克雷与诺伊曼边界条件的周期椭圆算子预解的算子误差估计。
• 量化同伦化逼近误差对小参数 $\varepsilon$ 和谱参数 $\zeta$ 的依赖关系。
• 在 $H^1$-范数下建立预解逼近的精确 $O(\varepsilon)$ 收敛速率，优于先前的 $O(\varepsilon^{1/2})$ 估计。

提出的方法

• 利用校正理论构造 $\mathcal{A}_{D,\varepsilon}$ 和 $\mathcal{A}_{N,\varepsilon}$ 预解的改进逼近。
• 通过预解恒等式与泛函演算，应用依赖谱参数的估计。
• 采用缩放变换与弗洛quet-布洛赫理论分析周期结构与谱行为。
• 在内部子区域 $\mathcal{O}'$ 上推导 $H^1(\mathcal{O}')$-范数下的误差界，显式依赖于到边界距离 $\delta$。
• 在预解中引入谱参数 $\zeta$，研究收敛速率对谱位置的依赖性。
• 通过辅助算子与投影技术控制 $L^2$ 与 $H^1$-范数下的误差，估计涉及 $c(\vartheta)$ 与 $\rho_\flat(\zeta)$。

实验结果

研究问题

• RQ1 同伦化椭圆算子的误差估计如何依赖于预解中的谱参数 $\zeta$？
• RQ2 尽管存在边界效应，是否可在诺伊曼问题中实现 $H^1$-范数下的精确 $O(\varepsilon)$ 收敛速率？
• RQ3 谱参数 $\zeta$ 在决定预解逼近误差界之精确性与结构中起何作用？
• RQ4 到边界的距离 $\delta$ 如何影响内部子区域中的收敛速率？
• RQ5 能否将诺伊曼型问题在 $H^1$-范数下的算子误差估计从 $O(\varepsilon^{1/2})$ 提升至 $O(\varepsilon)$？

主要发现

• 在 $H^1(\mathcal{O}')$ 中，$\mathcal{A}_{N,\varepsilon}$ 预解逼近的算子范数误差被控制为 $C \varepsilon \left( \delta^{-1} \rho_\flat(\zeta) + \delta^{-1} c(\vartheta)^{5/2} \rho_\flat(\zeta)^{3/4} \right)$，其中 $\delta = \mathrm{dist}(\mathcal{O}^\prime, \partial\mathcal{O})$。
• 在 $L^2(\mathcal{O}')$ 中，通量 $\mathbf{p}_\varepsilon = g^\varepsilon b(\mathbf{D}) \mathbf{u}_\varepsilon$ 的误差估计被控制为 $C \varepsilon \left( \delta^{-1} \rho_\flat(\zeta) + c(\vartheta)^{5/2} \rho_\flat(\zeta)^{3/4} \right)$，表明其为 $O(\varepsilon)$ 收敛。
• $\mathcal{A}_{N,\varepsilon}^{-1}$ 的 $L^2 \to H^1$ 算子范数误差被控制为 $C \varepsilon \left( \delta^{-1} \rho_\flat(\zeta) + c(\vartheta)^{1/2} \rho_\flat(\zeta)^{5/4} \right)$，常数与 $\varepsilon$ 和 $\zeta$ 无关。
• 校正项 $\varepsilon \widehat{K}_N(\varepsilon; \zeta)$ 在诺伊曼问题中实现了 $H^1(\mathcal{O}')$ 中的 $O(\varepsilon)$ 收敛，优于先前的 $O(\varepsilon^{1/2})$ 结果。
• 谱参数 $\zeta$ 的依赖性通过 $\rho_\flat(\zeta)$ 与 $c(\vartheta)$ 体现，二者控制谱间隙中误差的增长。
• 结果在阶数上是精确的，并可推广至狄利克雷与诺伊曼问题，显式依赖于谱参数与区域几何。