[论文解读] Stability of solutions to nonlinear diffusion equations

Teemu Lukkari|arXiv (Cornell University)|Jun 12, 2012

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

一句话总结

本文通过基于弱解和紧致性的初等估计，建立了非线性扩散方程（特别是多孔介质方程和快扩散方程）在非线性指数 $ m $ 的扰动下解的稳定性。证明了当 $ m_i \to m $ 时，解在 $ L^{1+m} $ 范数下收敛，即使跨越退化（$ m>1 $）和奇异（$ m<1 $）情形，且无需依赖正则性理论或半群方法。

ABSTRACT

We prove stability results for nonlinear diffusion equations of the porous medium and fast diffusion types with respect to the nonlinearity power $m$: solutions with fixed data converge in a suitable sense to the solution of the limit problem with the same data as $m$ varies. Our arguments are elementary and based on a general principle. We use neither regularity theory nor nonlinear semigroups, and our approach applies to e.g. Dirichlet problems in bounded domains and Cauchy problems on the whole space.

研究动机与目标

• 建立非线性扩散方程解对非线性幂 $ m $ 变化的稳定性，特别是临界范围 $ m > m_c = (n-2)_+/n $ 内的稳定性。
• 提供一个适用于退化（$ m>1 $）和奇异（$ m<1 $）情形的稳定性结果，包括极限情形 $ m=1 $（热方程），且无需假设 $ m \geq 1 $ 或 $ m \leq 1 $。
• 避免依赖非线性半群理论、Hölder 连续性或 Harnack 不等式等高级工具，仅使用基本的 $ L^2 $-型估计和紧致性。
• 将稳定性结果扩展至具有有限质量的正测度初值（如 Barenblatt 解）的情形，其具有 Dirac 初值迹。
• 为有界区域上的 Dirichlet 问题和 $ \mathbb{R}^n $ 上的柯西问题，提供 $ L^{1+m} $ 范数下收敛性的直接、初等证明。

提出的方法

• 使用方程 $ \partial_t u - \Delta u^m = 0 $ 的弱形式，其中 $ u^m $ 和 $ \nabla u^m $ 属于 $ L^2 $，从而可进行基于 $ L^2 $ 的能量估计。
• 应用紧致性原理：局部一致有界的弱解序列存在几乎处处收敛的子序列。
• 采用测试函数 $ \eta(x,t) = \int_t^T (u^m - u_i^{m_i}) \, ds $ 比较不同 $ m $ 下的解 $ u $ 和 $ u_i $，从而导出能量型估计。
• 利用不等式 $ |a-b|^{1+m} \leq c (a^m - b^m)(a-b) $（对 $ a,b > 0 $ 成立，且当 $ m > \frac{n-2}{n+2} $ 时成立），以控制解之间的差异。
• 结合 Hölder 不等式和点态估计 $ |u_i^{m_i} - u_i^m| \leq c(1 + u_i^{m_i+\varepsilon} + u_i^{m+\varepsilon})|m - m_i| $，以控制非线性项之间的差异。
• 结合 Sobolev 嵌入与能量估计，获得 $ u_i^{m_i(1+\varepsilon/m_i)(1+1/m)} $ 的一致 $ L^1 $ 有界性，确保对小的 $ \varepsilon $ 具有可积性。

实验结果

研究问题

• RQ1 非线性扩散方程 $ \partial_t u - \Delta u^m = 0 $ 的解是否对指数 $ m $ 连续依赖，尤其是在退化与奇异情形之间？
• RQ2 是否可以在不依赖非线性半群理论或 Harnack 不等式等正则性估计的前提下，建立解的稳定性？
• RQ3 当初始数据为正测度（如 Dirac delta）时，如 Barenblatt 解所示，解的收敛性是否仍然成立？
• RQ4 收敛性能否以 $ |m - m_i| $ 的形式量化，且在 $ L^{1+m} $ 范数下的收敛速率如何？
• RQ5 稳定性结果是否可扩展至 $ m \to 1 $ 的极限情形，包括极限方程为经典热方程的情况？

主要发现

• 当 $ m_i \to m $ 时，多孔介质方程和快扩散方程的解在 $ L^{1+m} $ 范数下收敛，且收敛速率满足 $ \|u - u_i\|_{L^{1+m}(\Omega_T)} \leq c |m - m_i|^{1/m} $（当 $ i $ 足够大时），前提是 $ m > \frac{n-2}{n+2} $。
• 稳定性结果对 $ m > 1 $（退化）和 $ m < 1 $（奇异）情形均成立，包括极限情形 $ m = 1 $，后者对应于热方程。
• 收敛性在有界区域上的 Dirichlet 问题和 $ \mathbb{R}^n $ 上的柯西问题中均成立，初始数据为有限质量的正测度，如 Barenblatt 解。
• Barenblatt 解 $ \mathcal{B}_m $ 满足 $ \mathcal{B}_{m_i} \to \mathcal{B}_m $ 在 $ L^p $ 中（当 $ m_i \to m $ 时），此结果此前在文献中尚未建立。
• 证明过程避免使用非线性半群理论、Hölder 正则性或反向 Hölder 不等式，仅依赖于弱解和紧致性。
• 对于 $ m \geq 1 $，存在一种替代证明方法，可给出解差值关于 $ |m - m_i| $ 的估计，但该方法仅适用于 $ m \geq 1 $，且依赖于强单调性。