[论文解读] Numerical approach to $L_1$-problems with the second order elliptic operators

N. Limić, Mladen Rogina|ArXiv.org|Dec 21, 2007

Differential Equations and Boundary Problems参考文献 14被引用 3

一句话总结

本文提出了一种用于求解具有狄利克雷边界条件的散度形式二阶椭圆算子的 $L^1$-问题的数值方法。该方法通过具有分 compartment（M-矩阵类似）结构的有限差分格式构造网格解，将其嵌入到帽函数中，并证明了在 $p \in [1, d/(d-1))$ 条件下强 $L^p$ 收敛和弱 $\dot{W}^1_p$ 收敛。关键贡献在于：在任意维数下，该方法具有可证明收敛性、单调性，并保持物理结构，同时具备严格的误差控制。

ABSTRACT

For a second order differential operator $A(\msx) =- abla a(\msx) abla + b'(\msx) abla+ abla \big(\msb''(\msx) \cdot\big)$ on a bounded domain $D$ with the Dirichlet boundary conditions on $\partial D$ there exists the inverse $T(λ, A)= (λI+A)^{-1}$ in $L_1(D)$. If $μ$ is a Radon (probability) measure on Borel algebra of subsets of $D$, then $T(λ, A)μ\in L_p(D), p \in [1, d/(d-1))$. We construct the numerical approximations to $u =T(λ, A)μ$ in two steps. In the first one we construct grid-solutions ${\bf u}_n$ and in the second step we embed grid-solutions into the linear space of hat functions $u(n) \in \dot{W}_p^1(D)$. The strong convergence to the original solutions $u$ is established in $L_p(D)$ and the weak convergence in $\dot{W}_p^1(D)$.

研究动机与目标

• 开发一种稳定、单调的数值格式，用于求解涉及散度形式二阶椭圆算子的 $L^1$-问题。
• 确保离散问题的系统矩阵具有分 compartment 结构（主对角线为正，非对角线为非正，列和非负），以模仿 M-矩阵。
• 在 $p \in [1, d/(d-1))$ 条件下，建立 $L^p(D)$ 中的强收敛与 $\dot{W}^1_p(D)$ 中的弱收敛。
• 将现有数值方法扩展至 $\mathbb{R}^d$ 中具有 Lipschitz 边界的任意区域，并处理奇异数据（Radon 测度）。
• 为 $d=2$ 提供构造性算法，并为 $d \geq 3$ 提供一般构造方法，同时给出扩散张量的充分条件。

提出的方法

• 在结构化网格上使用有限差分法对散度形式的二阶椭圆算子 $A(x)$ 进行离散化，确保系统矩阵 $A_n$ 的一致性和分 compartment 结构。
• 将网格解 $u_n$ 构造为离散线性系统 $A_n u_n = \mu_n$ 的解，其中 $\mu_n$ 为离散 Radon 测度。
• 通过插值将网格解 $u_n$ 嵌入到连续分片线性（帽）函数空间 $u^{(n)} \in \dot{W}^1_p(D)$ 中。
• 利用紧致性与一致性论证，证明 $u^{(n)}$ 在 $L^p(D)$ 中强收敛于真实解 $u = T(\lambda, A)\mu$，且在 $\dot{W}^1_p(D)$ 中弱收敛。
• 使用双线性型 $a(v,u)$ 和格林公式，将连续与离散的变分形式联系起来。
• 对于 $d \geq 3$，对扩散张量 $\{a_{ij}(x)\}$ 施加假设 4.1，以保证 $A_n$ 的分 compartment 结构。

实验结果

研究问题

• RQ1 能否为任意维数下具有一般扩散张量的散度形式二阶椭圆算子构造出单调的有限差分格式？
• RQ2 当离散解被嵌入帽函数后，是否在 $p \in [1, d/(d-1))$ 条件下于 $L^p(D)$ 中强收敛、于 $\dot{W}^1_p(D)$ 中弱收敛？
• RQ3 扩散张量的何种条件可保证系统矩阵 $A_n$ 的分 compartment 结构？
• RQ4 该方法在具有奇异数据（如狄拉克测度）的问题中数值表现如何？
• RQ5 该格式能否在 $d=2$ 中以直接算法结构高效实现？

主要发现

• 对于 $d=2$，该格式具有显式算法结构，可直接应用，且在 $W^1_2$ 与 $W^1_p$ 空间中已证明收敛。
• 在 $d \geq 3$ 时，若张量值函数 $\hat{a}$ 严格正定，则根据引理 8.1，$A_n$ 的分 compartment 结构可保证成立。
• 数值示例 7.1 显示，对于含狄拉克测度的问题，$\varepsilon_{1,\text{rel}} = 0.003371$ 且 $\varepsilon_{\infty,\text{rel}} = 1.6254$，最大误差出现在奇点附近。
• 示例 7.2 经过 210 次迭代后，$\varepsilon_{1,\text{rel}} = 8 \times 10^{-6}$ 且 $\varepsilon_{\infty,\text{rel}} = 4 \times 10^{-3}$，确认了高精度。
• 该方法通过使用扩展 stencils 避免了在 $d=2$ 中的旋转操作，保持了单调性，并避免了解决方案曲面上的锯齿状伪影。
• 该格式对奇异数据（Radon 测度）具有鲁棒性，其逆算子 $T(\lambda, A)$ 将 $\mu \in R(D)$ 映射至 $L^p(D)$ 与 $\dot{W}^1_p(D)$，其中 $p < d/(d-1)$。