本文研究如下带有扩散项的非线性抛物方程的初边值问题
其中, $ (\int_\Omega|\nabla u|^2dx)^\gamma\Delta u $为非局部扩散项, $ \Omega\subset R^n $为边界充分光滑的有界区域, 常数$ \gamma\in R $, $ p>1 $.
(1.1)可用来描述热传导现象、半导体中的电子与空穴流、流体在多孔介质中的运动规律等问题. 当$ \gamma = 0 $时, (1.1)是经典的抛物型方程, 与其相关的抛物方程整体解的存在和爆破问题已被许多学者进行了广泛的研究, 详见[1-6]. 当$ \gamma\neq0 $时, (1.1)即为非局部反应扩散方程, 与之相关的抛物方程的研究也有一些成果. 例如, 在文献[7]中, Sismen等考虑了如下非线性抛物方程的初边值问题
其中, $ \Omega\subset R^n $为边界充分光滑的有界区域, 常数$ p\geq 2 $. $ f:R\rightarrow R $是利普希茨连续函数. $ l:L^2(\Omega)\rightarrow R $是线性连续函数, 即存在一个函数$ g\in L^2(\Omega) $, 使得对所有的$ u\in L^2(\Omega) $, 都有$ l(u) = l_g(u) = \int_\Omega g(x)u(x)dx $成立. $ 0< m\leq a(s)\leq M(s\in R) $. 采用Galerkin方法和Aubin紧性定理, 研究了弱解的整体存在性、唯一性以及关于初值的连续性.
事实上, 在以往大多数文献中, 均假设扩散项系数是有界的, 即存在正常数$ m, M $, 使得$ 0< m\leq a(s)\leq M<\infty(s\in R) $, 此时问题总是非退化的, 如文献[8-10]. 但在刻画扰动传播的有限性等问题的反应扩散方程时总会出现退化情形, 因此有必要研究具有退化性的反应扩散方程. 如Almeida等在文献[11]中研究了如下非局部退化抛物问题
此时, 扩散项系数$ (\int_\Omega u^2dx)^\gamma(\gamma\in R) $可能是零或无穷大. 文[11]采用扰动方法, 将退化问题转化成非退化问题, 利用文献[12]中非退化问题弱解的存在性, 通过建立解的相应先验估计, 得到了退化问题弱解整体唯一存在和局部唯一存在的条件, 并研究了解的渐近行为.
本文旨在研究另一种扩散项为$ (\int_\Omega |\nabla u|^2dx)^\gamma\Delta u(\gamma\in R) $的情形, 此时, 扩散项系数$ (\int_\Omega |\nabla u|^2dx)^\gamma(\gamma\in R) $也可能是零或无穷大, 且进一步将文献[11]中非齐次项$ f(x, t) $拓展为非线性项$ |u|^{p-1}u $, 受文献[13]的启发, 结合Gagliardo-Nirenberg不等式、Sobolev嵌入定理和势井理论, 利用Galerkin方法, 建立了弱解整体唯一存在和爆破的充分条件.
首先给出弱解的定义.
定义1 若$ u\in L^\infty(0, T; H^1(\Omega)) $, $ \frac{\partial u}{\partial t}\in L^2(0, T; L^2(\Omega)) $, 对任意的$ w\in H_0^1(\Omega) $, $ t\in[0, T] $, 都有等式
成立, 且
则称函数$ u $是问题(1.1)的一个弱解.
现给出本文主要结论如下.
定理1 假设$ u_0(x)\in H_0^{1}(\Omega) $, $ \gamma>\frac{p-1}{2} $, 当$ n> 2 $时, $ 1< p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 1< p< \infty $. 则问题(1.1) 存在一个整体弱解$ u(x, t) $, 满足
若$ n< 2 $, 则由初值函数$ u_0(x) $决定的整体弱解$ u(x) $是唯一的.
定理2 假设$ 0<\gamma< \frac{p-1}{2} $, 当$ n> 2 $时, $ 1< p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 1< p< \infty $. 则对每一个包含在势井$ W $中的初值函数$ u_0(x) $, 初边值问题(1.1)都存在一个包含在$ \overline W $中的整体弱解$ u(x, t) $, 使得
进一步的, 有$ t> 0 $, 因此$ u(x, t) $满足
其中, $ \overline W\rm{是}W\rm{在}H_0^1(\Omega)\rm{中的闭包} $, $ W $的定义见(2.6).
定理3 假设$ u_0(x)\in H_0^1(\Omega) $, $ 0<\gamma<\frac{p-1}{2} $, $ 1< p<\frac{n+4}{n} $. 则存在一个正常数
使得当$ 0\leq t\leq T_0 $时, 问题(1.1)存在一个局部弱解$ u(x, t) $, 满足
并且有
和
成立. 其中$ J(u) $的定义见(2.5).
定理4 假设定理3的条件全部满足. 进一步的假设
且
则问题(1.1)满足初值函数的弱解不是全局的, 即在有限时间
内爆破.
定理5 假设定理1的条件成立, 对任意固定的$ T>0 $, 令$ u $, $ v $是问题(1.1)满足$ u(0) = u_0 $, $ v(0) = v_0 $的两个整体弱解, 若$ n<2 $, 则有
全文将$ \|\cdot\|_{L^p(\Omega)} $和$ \|\cdot\|_{H^s(\Omega)} $分别简记为$ \|\cdot\|_{L^p} $和$ \|\cdot\|_{H^s} $, 记$ (u, v) = \int_\Omega uvdx $, $ c $表示不依赖于未知函数的正常数, 在不同地方不一定相同.
首先给出结论证明过程中需要的两个引理.
引理1 [12] (Gagliardo-Nirenberg不等式) 对任意函数$ u(x)\in W_0^{1, p}(\Omega) $, $ p\geq1 $且$ r\geq1 $, 有如下不等式成立,
其中
(1) 当$ p\geq n = 1 $时, $ r\leq q\leq\infty; $
(2) 当$ n> 1 $且$ p< n $时, 若$ r\leq\frac{np}{n-p} $, 则$ r\leq q\leq\frac{np}{n-p}; $若$ r\geq\frac{np}{n-p}, $则$ \frac{np}{n-p}\leq q\leq r; $
(3) 当$ p = n> 1 $时, $ r\leq q\leq\infty; $
(4) 当$ p> n> 1 $时, $ r\leq q\leq\infty. $
常数$ C_1 $仅依赖于$ n, p, q, r. $
引理2 [12] (Sobolev嵌入定理) 对任意函数$ u(x)\in H_0^{1}(\Omega) $, 有
其中, 当$ n> 2 $时, $ 0\leq p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 0\leq p< \infty $, 常数$ C_2 $仅依赖于$ n, p, \Omega. $当$ n< 2 $时, $ H_0^1(\Omega) $中的函数是连续的, 且
常数$ C_3 $仅依赖于$ n, \Omega. $
下面引入能量泛函
定义势井
其中势井深度
接下来, 给出势井深度和势井的相关性质.
引理3 假设$ u\in H_0^{1}(\Omega) $, $ -1< \gamma< \frac{p-1}{2} $, 当$ n> 2 $时, $ 0\leq p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 0\leq p< \infty $. 则$ d>0. $
证 显然
考虑对任意$ \lambda> 0 $, 有
令$ h(\lambda) = \lambda^{2\gamma+1-p}\|\nabla u\|_{L^2}^{2\gamma+1}-\|u\|_{L^{p+1}}^{p+1}. $因为$ -1< \gamma< \frac{p-1}{2} $, 所以
并且
由(2.9), (2.10)可以得到, 存在唯一一个$ \lambda^* = \left(\frac{\|\nabla u\|_{L^2}^{2\gamma+2}}{\|u\|_{L^{p+1}}^{p+1}}\right)^{\frac{1}{p-2\gamma-1}}> 0 $使得$ h(\lambda^*) = 0 $. 进一步的, 由(2.8)有, $ \frac{d}{d\lambda}J(\lambda u)|_{\lambda = \lambda^*} = \lambda^ph(\lambda). $因为当$ \lambda\in(0, \lambda^*) $时, $ h(\lambda)> 0 $; 当$ \lambda\in(\lambda^*, +\infty) $时, $ h(\lambda)< 0 $. 因此可以得到, 当$ \lambda\in(0, \lambda^*) $时, $ J(\lambda u) $单调递增; 当$ \lambda\in(\lambda^*, +\infty) $时, $ J(\lambda u) $单调递减. 由于$ J(\lambda u)|_{\lambda = 0} = 0, $$ -1< \gamma< \frac{p-1}{2} $及当$ n> 2 $时, $ 0\leq p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 0\leq p< \infty $, 根据引理2可以得到
引理4 令$ W_* = \left\{u|u\in H_0^{1}(\Omega), \|\nabla u\|_{L^2}^{2\gamma+2}-\|u\|_{L^{p+1}}^{p+1}>0, J(u)< d\right\}, $则
证 $ (1) $假设$ u\in W\rm{且}u\not = 0 $, 则
因此$ \left(\frac{\|\nabla u\|_{L^2}^{2\gamma+2}}{\|u\|_{L^{p+1}}^{p+1}}\right)^{\frac{1}{p-2\gamma-1}}> 1, $故$ u\in W_*. $
(2) 反过来, 假设$ u\in W_* $, 则有
而当$ 0<\lambda\leq 1 $时, $ \frac{\partial}{\partial\lambda}J(\lambda u)> 0 $且$ J(\lambda u)|_{\lambda = 0} = 0 $, 因此$ \sup_{\lambda\in[0, 1]}J(\lambda u) = J(u)< d. $
综上$ W = W_*\cup\left\{0\right\}. $
下面的定理与$ d $的上界有关.
引理5 考虑以下非线性正特征值问题
其中, $ -1< \gamma< \frac{p-1}{2} $, 当$ n> 2 $时, $ 0\leq p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 0\leq p<\infty $. 若对$ \lambda> 0 $, $ (2.11) $有一特征函数$ u_\lambda $, 则$ d $是有上界的, 且势井$ W $在$ H_0^1(\Omega) $中是有界的.
证 一方面,
另一方面, 由(2.11)可以得到
因此
由以上可知$ d $是有上界的.
令$ u\in W $且$ u\not = 0 $, 则由引理4得到
由于$ -1< \gamma< \frac{p-1}{2} $, 即
所以, 势井$ W $包含在球$ \left\{v|v\in H_0^1(\Omega), \|\nabla v\|_{L^2}^{2\gamma+2}\leq\frac{(2\gamma+2)(p+1)}{p-2\gamma-1}d\right\} $中.
证明定理1 利用Galerkin方法. 取空间$ H_0^{1}(\Omega) $的一组基$ \{\omega_i\}_{i = 1}^\infty $, 且$ \omega_i $是$ -\Delta $算子的狄利克雷问题的特征函数:
将(3.1)两边同时与$ \omega_i $做内积, 可知特征值$ \lambda_i\geq 0, i\in N^*, $将$ \omega_i $标准化, 使得$ \|\omega_i\|_{L^2} = 1 $. 接下来, 寻找(1.1)具有以下形式的近似解
其中未知函数$ g_{im} $是由以下常微分方程决定,
且初值条件满足
因此, 由常微分方程中Picard迭代法, (3.3), (3.4)在某一区间$ [0, t_m], 0< t_m< T $中存在一个局部解$ u_m(t) $. 下面证明, 对任意的$ T> 0 $, 这样的解可以通过先验估计延拓到整个区间$ [0, T] $上.
在(3.3)的两边同时乘以$ g'_{im} $, 对$ i $从1到$ m $求和, 并关于$ t $积分, 得到
令$ f(t) = \|\nabla u_m\|_{L^2}^2 $, 则(3.5)左边第二项
换回原函数名, 即
将(3.6)代入(3.5)整理可得
由于$ \gamma>\frac{p-1}{2} $, 当$ n> 2 $时, $ 1< p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 1< p< \infty $. 利用引理2和Young不等式, 则有
由此可以得到
利用(3.9)和Poincar$ \acute{e} $不等式,
因此, 由先验估计(3.8)–(3.10)可知, 存在一个函数$ u $和$ \left\{u_m\right\} $的子序列, 不妨仍记为$ \left\{u_m\right\} $, 使得
一方面, 由Aubin紧性定理, $ W = \{v|v\in L^\infty(0, T; H_0^1(\Omega)), v_t\in L^2(0, T; L^2(\Omega))\} $紧嵌入到$ L^\infty(0, T; H_0^1(\Omega)) $中, 根据(3.11)–(3.13)可以得到
于是, 存在$ \left\{u_m\right\} $的子序列, 不妨仍记为$ \left\{u_m\right\} $, 使得$ u_m\stackrel{a.e}{\longrightarrow}u $\ 在$ \Omega\times\left[0, T\right] $中. 因为$ s\rightarrow |s|^{p-1}s $是一个连续函数, 所以
由$ u_m $在$ L^\infty(0, T; H_0^1(\Omega)) $中的一致有界性和引理2, $ |u_m|^{p-1}u_m $在$ L^\infty(0, T; L^{\frac{p+1}{p}}(\Omega)) $上一致有界. 则
另一方面, 由(3.14)可知
为了讨论方便, 记$ a(u) = (\int_\Omega|\nabla u|^2dx)^\gamma $, 由于$ a $是连续的, 可以得到
固定(3.3)中的$ j $, 由于$ \frac{p+1}{p}< 2 $, 且每一项均在$ L^{\frac{p+1}{p}}[0, T] $中弱收敛, 令$ m\longrightarrow\infty $, 可以得到
根据基$ \{\omega_i\}_{i = 1}^\infty $的稠密性,
由(3.4), (3.13), (3.14)知, 在$ L^2(\Omega) $中, $ u_m(0)\rightharpoonup u(0) $, $ u_m(0)\rightharpoonup u_0. $根据极限的唯一性, $ u(0) = u_0. $因此$ u $是问题(1.1)的整体弱解.
唯一性的证明如下.
令$ u_1 $和$ u_2 $是(1.1)的两个满足初值条件的整体弱解, 则$ w = u_1-u_2 $满足
接下来, 将(3.15)与$ w $做内积可得
利用拉格朗日中值定理和Hölder不等式, (3.16)右边
而左边第二项和第三项利用Hölder不等式可化简为
因此, 由引理2可得
所以, 根据Gronwall不等式, $ w = 0. $
证明定理2 再次利用Galerkin方法. $ \{\omega_i\}_{i = 1}^\infty $, $ \{u_m\} $与定理1中所述相同, 令$ \left\{u_{0m}\right\} $是一个满足如下条件的序列:
则在区间$ \left[0, t_m\right](t_m> 0) $上存在解$ u_m(t) $, 且根据(3.7), 在该区间上$ u_m(t) $满足
下面证明
假设$ (4.3) $不成立, 设$ t^* $是最小的时间, 使得$ u_m(t^*)\not\in W $, 则由$ u_m(t) $的连续性可知$ u_m(t^*)\in \partial W. $因此, 由引理4可以得到
或者
这与$ (4.1), (4.2) $矛盾. 事实上, 当$ (4.4) $成立时, 显然矛盾. 当$ (4.5) $成立时, 有
这也意味着矛盾. 因此, 由$ (4.2) $和引理4得到
因为$ 0<\gamma<\frac{p-1}{2} $, 所以$ \frac{1}{2\gamma+2}-\frac{1}{p+1}> 0 $, 由此可知
同定理1的证明可以得到
由(5.7)–(5.9)和Aubin紧性定理, 重复定理1的证明过程, 可得到$ u $是问题(1.1)的整体弱解.
下面证明(1.3).
根据引理4, $ u\in \overline W $, 所以
将(1.1)与$ u $做内积并关于$ t $积分, 整理化简可得
因此, 由(4.10)和(4.11)立即可以得到(1.3).
唯一性的证明重复定理1的证明过程可以得到.
证明定理3 依旧采用Galerkin方法. $ \{\omega_i\}_{i = 1}^\infty $, $ \{u_m\} $, $ \{u_{0m}\} $与定理1中所述相同. 在(3.2)的两边同时乘以$ g_{im} $并对$ i $从1到$ m $求和, 得到
因为当$ n> 2 $时, $ 1< p\leq\frac{n+2}{n-2} $; 当$ n\leq2 $时, $ 1< p< \infty $, 所以$ \frac{1}{2}-\frac{1}{n}\leq\frac{1}{p+1}<\frac{1}{2}, $于是根据引理1可以得到
由于$ \gamma>\frac{np-n-4}{4} $, 所以
又因为$ \gamma<\frac{p-1}{2} $, 故有
其中, $ \theta $是由$ \frac{1}{p+1} = (\frac{1}{2}-\frac{1}{n})\theta+\frac{1}{2}(1-\theta) $唯一确定的[0, 1]中的数. 于是根据Young不等式,
其中$ \delta^{''} = \frac{\delta^{'}}{2(1-\delta)}>1. $因此, 由(5.1), (5.3)得到
于是, 微分不等式(5.4)的解$ \|u_m\|_{L^2}^{2} $可由以下初值问题的解来表示:
如果
则(5.5)的解$ y $是有限的. 因此, 在区间$ 0\leq t\leq T_0 = \frac{t_\infty}{2} $内, 可以得到
由(3.7), (5.3), (5.7), 立即得到
由(5.8)–(5.10)和Aubin紧性定理, 重复定理1的证明过程, 可以得到$ u $是问题(1.1)在区间$ 0\leq t\leq T_0 $上的弱解.
下满证明$ u $满足(1.4)和(1.5).
令$ V(0, T; H_0^1(\Omega)) $是由满足以下条件的所有函数$ v(t) $构成的空间, 即
在(3.2)的两边同时乘以任意函数$ f(t)\in C^1([0, T]) $, 并关于$ t $积分, 得到
固定$ j $, 令$ m\rightarrow \infty $则有
这意味着对任意$ \psi(\tau)\in V(0, T; H_0^1(\Omega)) $都有
在(5.11)中令$ \psi(\tau) = u(\tau) $可立即得到(1.4).
从(3.7)可以得到
令$ \theta\in C([0, T_0]) $是非负函数, 则有
其中, 当$ m\longrightarrow \infty $时, (5.13)第二项趋于$ \int_{0}^{T_0}J(u_0(t))\theta(t)dt. $而第一项关于$ L^2(0, T_0; H_0^1(\Omega)) $的弱拓扑是下半连续的. 因此
因为$ \theta $是任意的, 所以
定理3的证明完成.
注1 当$ \delta^{''}>1 $时, $ T_0> 0 $. 下面证明$ \delta^{''} = \frac{\delta^{'}}{2(1-\delta)} = \frac{(n-np+2p+2)(\gamma+1)}{n-np+4\gamma+4}> 1. $整理化简即证$ 2\gamma-n\gamma+2> 0 $. 当$ n\leq 2 $时, 显然成立. 当$ n> 2 $时, 即证$ \gamma<\frac{2}{n-2} $. 因为当$ n> 2 $时, $ 1< p\leq\frac{n+2}{n-2} $, 所以$ \frac{2}{n-2}-\frac{p-1}{2} = \frac{n-np+2p+2}{2n-4}> 0 $. 于是$ \gamma<\frac{p-1}{2}<\frac{2}{n-2} $. 综上, 结论成立.
证明定理4 假设定理4的结论不成立, 令$ u(x, t) $是由初值函数$ u_0(x) $决定的满足定理条件的全局解, 则对$ \forall t> 0 $, $ u(x, t) $满足(1.4)和(1.5). 从(1.5)和(1.6)可以得到
将(5.14)代入(1.4)中, 则有
这里, 用到不等式$ \|u(t)\|_{L^2}^2\leq(\rm{mes}(\Omega))^{\frac{p-1}{p+1}}\|u(t)\|_{L^{p+1}}^2. $由(5.15)立即得到
因此, 当$ t\longrightarrow\frac{1}{L}\|u_0\|_{L^2}^{1-p}> T $时, $ \|u(t)\|_{L^2}^{p-1}\longrightarrow\infty $. 这与$ u $是全局解矛盾. 故定理成立.
证明定理5 令$ u, v $是问题(1.1)满足$ u(0) = u_0, v(0) = v_0 $的两个弱解, 则
同定理1中弱解的唯一性证明过程,
将(6.2), (6.3)代入(6.1)中, 且由于$ u, v\in L^\infty(0, T; H_0^1(\Omega)), n< 2 $, 根据引理2可以得到
由于$ v $是问题(1.1)的解, 将$ v $代入方程(1.1), 两边同时与$ v_t $做内积并关于$ t $积分可以得到
故
同理
将(6.5), (6.6)代入(6.4), 整理化简可得
在(6.7)的两边同时乘以$ e^{-\eta t} $, 计算得
其中$ C(T, \|\nabla u_0\|_{L^2}, \|u_0\|_{L^{p+1}}, \|\nabla v_0\|_{L^2}, \|v_0\|_{L^{p+1}}) = e^{\eta T}. $
注2 假设定理2的条件成立, 对任意固定的$ T>0 $, 令$ u $, $ v $是问题(1.1)满足$ u(0) = u_0 $, $ v(0) = v_0 $的两个整体弱解, 若$ n<2 $, 通过类似的证明过程, 可得到相似的结果.