非线性扩散型方程解的一阶导数渐近平坦性条件问询

阿华AIGC实验室

2026-4-22

嗨，这个问题挺有意思的——咱们先从你提到的经典扩散方程说起，理清这类问题的基本思路后，再延伸到非线性的情况：

经典扩散方程的证明思路

对于线性扩散方程 $\frac{\partial u}{\partial t}=\frac{\partial^2u}{\partial x^2}$，有两种常见方法能证明 $\lim_{t\to\infty}\frac{\partial u}{\partial x}=0$：

方法1：利用解的积分表达式

经典扩散方程的解是初始条件与高斯核的卷积：
$$u(x,t) = \frac{1}{\sqrt{4\pi t}} \int_{\mathbb{R}} g(y) e^{-(x-y)2/(4t)} dy$$
对 $x$ 求一阶导数：
$$\frac{\partial u}{\partial x} = -\frac{1}{\sqrt{4\pi t}} \int_{\mathbb{R}} g(y) \frac{x-y}{2t} e^{-(x-y)2/(4t)} dy$$
做变量替换 $z = \frac{x-y}{2\sqrt{t}}$，代入后化简得：
$$\frac{\partial u}{\partial x} = -\frac{1}{\sqrt{2\pi} t} \int_{\mathbb{R}} g(x-2\sqrt{t}z) z e^{-z2} dz$$
如果初始条件 $g(x)$ 有界（比如 $|g(x)| \leq M$），或在无穷远处足够快衰减（比如属于Schwartz函数类），那么积分项是有限的。当 $t\to\infty$ 时，$\frac{1}{t}$ 趋于0，整个表达式自然也趋于0。

方法2：能量估计

定义一阶导数的 $L^2$ 范数为能量泛函：$E(t) = \int_{\mathbb{R}} \left(\frac{\partial u}{\partial x}\right)^2 dx$，对 $t$ 求导：
$$\frac{dE}{dt} = 2\int_{\mathbb{R}} \frac{\partial u}{\partial x} \frac{\partial^2 u}{\partial x \partial t} dx$$
结合扩散方程替换后分部积分（假设解在无穷远处衰减足够快，边界项为0）：
$$\frac{dE}{dt} = 2\int_{\mathbb{R}} \frac{\partial u}{\partial x} \frac{\partial^3 u}{\partial x^3} dx = -2\int_{\mathbb{R}} \left(\frac{\partial^2 u}{\partial x^2}\right)2 dx \leq 0$$
这说明 $E(t)$ 是随时间递减的非负函数，必有极限 $E(\infty) \geq 0$。进一步可证明这个极限为0：若 $E(\infty) >0$，则 $\frac{\partial^2 u}{\partial x^2}$ 会几乎处处趋于0，进而推出 $\frac{\partial u}{\partial x}$ 趋于常数，但结合解的衰减性或积分守恒性，这个常数只能是0。

非线性情况：以 $F(s)=s^\alpha$ 为例

接下来聊聊你关注的幂函数形式 $\frac{\partial u}{\partial t} = \left(\frac{\partial^2 u}{\partial x^2}\right)\alpha$，这里需要分 $\alpha$ 的取值范围讨论，同时依赖初始条件的正则性和衰减性：

1. $\alpha \geq 1$ 的情况

这种情况下依然可以用能量估计的思路：
定义 $E(t) = \int_{\mathbb{R}} \left(\frac{\partial u}{\partial x}\right)^2 dx$，求导后分部积分：
$$\frac{dE}{dt} = -2\alpha \int_{\mathbb{R}} \left(\frac{\partial^2 u}{\partial x^2}\right){\alpha+1} dx$$
当 $\alpha \geq1$ 且 $\frac{\partial^2 u}{\partial x^2}$ 保持符号（比如初始条件是凸/凹函数，方程不会改变二阶导数的符号），积分项非负，因此 $E(t)$ 严格递减且有下界0，最终趋于0，即 $\frac{\partial u}{\partial x}$ 的 $L^2$ 范数趋于0。结合最大值原理控制导数的逐点界，就能证明逐点趋于0。

另外，可通过自相似解分析渐近行为：假设 $u(x,t)=t^\beta f(x/t^\gamma)$，代入方程后可推导出指数关系 $\beta-1 = \alpha(\beta-2\gamma)$。若希望解趋于平坦，需 $\beta-\gamma <0$，比如当解趋于常数时 $\beta=0$，此时 $\gamma=1/(2\alpha)$，对应的自相似解导数 $f'(z)$ 会随 $|z|$ 衰减，因此当 $t\to\infty$ 时，$\frac{\partial u}{\partial x}$ 逐点趋于0。

2. $0 < \alpha <1$ 的情况

能量估计依然成立：$\frac{dE}{dt} = -2\alpha \int_{\mathbb{R}} \left(\frac{\partial^2 u}{\partial x^2}\right){\alpha+1} dx \leq0$，$E(t)$ 递减趋于某个非负极限。但需要额外条件确保极限为0：比如初始条件 $g(x)$ 具有有限的“高阶能量”（$\int_{\mathbb{R}} |g''(x)|^{\alpha+1} dx <\infty$），或解在无穷远处的衰减速度足够快。

这种情况下方程的扩散速度比线性情况慢，但只要初始条件不是太“极端”，依然能证明 $\frac{\partial u}{\partial x}$ 趋于0。

3. $\alpha <0$ 的情况

此时 $F(s)=s^\alpha = 1/|s|^{|\alpha|}$（假设 $s\neq0$），方程行为会变得很反常：若 $\frac{\partial^2 u}{\partial x^2}$ 为正，$\frac{\partial u}{\partial t}$ 为正，解会快速增长；若 $\frac{\partial^2 u}{\partial x^2}$ 为负，$\frac{\partial u}{\partial t}$ 为负（$\alpha$ 奇数）或正（$\alpha$ 偶数），可能导致解出现奇性或越来越陡，大概率无法满足 $\lim_{t\to\infty}\frac{\partial u}{\partial x}=0$。