[论文解读] Homogenization of the Neumann problem for higher-order elliptic equations with periodic coefficients

T. A. Suslina|arXiv (Cornell University)|May 21, 2017

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

一句话总结

本文为具有周期系数的高阶椭圆诺伊曼问题的均匀化建立了精确的算子误差估计。它推导了在 $ L_2 $ 和 $ H^p $ 范数下对预解式 $ (A_{N, ho}- heta I)^{-1} $ 的逼近，显示出误差界为 $ \varepsilon^{1/2}\rho_\flat(\zeta)^{1/2} + \varepsilon|\zeta+1|\rho_\flat(\zeta) $，并引入了校正项 $ \varepsilon^p K_N^0(\zeta;\varepsilon) $，在微分算子和系数矩阵满足最优正则性和秩条件的前提下成立。

ABSTRACT

Let $\mathcal{O}\subset\mathbb{R}^d$ be a bounded domain of class $C^{2p}$. In $L_2(\mathcal{O};\mathbb{C}^n)$, we study a selfadjoint strongly elliptic operator $A_{N,\varepsilon}$ of order $2p$ given by the expression $b({\mathbf D})^* g({\mathbf x}/\varepsilon) b({\mathbf D})$, $\varepsilon >0$, with the Neumann boundary conditions. Here $g({\mathbf x})$ is a bounded and positive definite $(m imes m)$-matrix-valued function in ${\mathbb R}^d$, periodic with respect to some lattice; $b({\mathbf D})=\sum_{|α|=p} b_α{\mathbf D}^α$ is a differential operator of order $p$ with constant coefficients; $b_α$ are constant $(m imes n)$-matrices. It is assumed that $m\geqslant n$ and that the symbol $b({\boldsymbol ξ})$ has maximal rank for any $0 e {\boldsymbol ξ}\in {\mathbb C}^d$. We find approximations for the resolvent $\left(A_{N,\varepsilon}-ζI ight)^{-1}$ in the $L_2(\mathcal{O};\mathbb{C}^n)$-operator norm and in the norm of operators acting from $L_2(\mathcal{O};\mathbb{C}^n)$ to the Sobolev space $H^p(\mathcal{O};\mathbb{C}^n)$, with error estimates depending on $\varepsilon$ and $ζ$.

研究动机与目标

• 分析在有界 $ C^{2p} $ 域中，具有周期系数的自伴、强椭圆、高阶 $ 2p $ 阶微分算子的均匀化问题。
• 在 $ \varepsilon \to 0 $ 时，推导预解式 $ (A_{N,\rho} - \zeta I)^{-1} $ 在 $ L_2 $ 和 $ H^p $ 空间中的算子范数误差估计，其中 $ \varepsilon $ 为周期尺度。
• 构造一个校正项 $ \varepsilon^p K_N^0(\zeta;\varepsilon) $，以改善在 $ H^p $ 范数下对预解式的逼近。
• 建立依赖于谱参数 $ \zeta \in \mathbb{C} \setminus [c_\flat, \infty) $ 的误差界，同时捕捉低频与高频行为。

提出的方法

• 分析基于缩放变换、弗洛凯-布洛赫理论和解析摄动理论的算子理论方法。
• 将微分算子分解为 $ A_{N,\varepsilon} = b(\mathbf{D})^* g(\mathbf{x}/\varepsilon) b(\mathbf{D}) $，其中 $ b(\mathbf{D}) $ 是一个具有最大秩符号的 $ p $ 阶常系数算子。
• 通过 $ g(\mathbf{x}) $ 的周期结构导出的均匀化矩阵 $ g^0 $ 构造有效算子 $ A_N^0 $。
• 通过傅里叶域中的谱分析，利用预解式恒等式以及对谱投影算子和校正算子的估计，推导误差估计。
• 显式定义了校正项 $ K_N^0(\zeta;\varepsilon) $，并证明其能捕捉解在尺度 $ \varepsilon $ 下的振荡行为，且满足 $ \|K_N^0(\zeta;\varepsilon)\|_{L_2 \to H^p} = O(\varepsilon^{-p}) $。
• 在两种算子范数下证明了估计：$ L_2 \to L_2 $ 和 $ L_2 \to H^p $，其界依赖于 $ \varepsilon $、$ \zeta $ 和谱函数 $ \rho_\flat(\zeta) $。

实验结果

研究问题

• RQ1 当 $ \varepsilon \to 0 $ 时，预解式 $ (A_{N,\varepsilon} - \zeta I)^{-1} $ 在 $ L_2 $-算子范数下收敛到均匀化预解式 $ (A_N^0 - \zeta I)^{-1} $ 的收敛速率是多少？
• RQ2 如何通过在 $ H^p $-范数下引入校正项，将收敛性改进至超越一阶逼近？
• RQ3 误差对谱参数 $ \zeta \in \mathbb{C} \setminus [c_\flat, \infty) $ 的依赖关系如何？
• RQ4 在何种条件下校正项消失或简化，这又如何影响误差估计？
• RQ5 能否为具有周期系数的高阶椭圆诺伊曼问题推导出精确的算子误差估计，从而推广此前针对二阶算子的结果？

主要发现

• 预解式 $ (A_{N,\varepsilon} - \zeta I)^{-1} $ 在 $ L_2(\mathcal{O};\mathbb{C}^n) \to L_2(\mathcal{O};\mathbb{C}^n) $ 范数下收敛至 $ (A_N^0 - \zeta I)^{-1} $，误差满足 $ \leq \mathfrak{C}_1 \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) $，其中 $ 0 < \varepsilon \leq \varepsilon_1 $。
• 在 $ L_2(\mathcal{O};\mathbb{C}^n) \to H^p(\mathcal{O};\mathbb{C}^n) $ 范数下，误差被控制在 $ \widetilde{\mathfrak{C}}_{10} \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) $ 之内，且包含校正项 $ \varepsilon^p K_N^0(\zeta;\varepsilon) $。
• 当 $ g^0 = \overline{g} $ 时，校正项消失，误差估计简化为 $ \| \mathbf{u}_\varepsilon - \mathbf{u}_0 \|_{H^p} \leq \mathfrak{C}_{10} \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) \| \mathbf{F} \|_{L_2} $。
• 对于通量 $ \mathbf{p}_\varepsilon = g^\varepsilon b(\mathbf{D}) \mathbf{u}_\varepsilon $，其在 $ L_2 $-范数下的误差被控制在 $ \widetilde{\mathfrak{C}}_{11} \left( \varepsilon^{1/2} \rho_\flat(\zeta)^{1/2} + \varepsilon |\zeta+1| \rho_\flat(\zeta) \right) \| \mathbf{F} \|_{L_2} $ 之内。
• 当 $ g^0 = \underline{g} $ 时，误差估计变为 $ \| \mathbf{u}_\varepsilon - \mathbf{v}_\varepsilon^0 \|_{H^p} \leq \widehat{\mathfrak{C}}_{10} \varepsilon \left( \rho_\flat(\zeta)^{1/2} + |\zeta+1| \rho_\flat(\zeta) \right) \| \mathbf{F} \|_{L_2} $，显示出更优的 $ \varepsilon $-依赖性。
• 这些估计是精确的，且显式依赖于谱参数 $ \zeta $，其中 $ \rho_\flat(\zeta) $ 编码了预解式的谱间隙和衰减性质。