在过去的几十年间, 非线性抛物方程解的爆破问题已经得到了很广泛的研究, 文献[1, 2]有较详细的介绍.其中, 关于解爆破时间$ T $的研究, 仍然有许多文章仅针对$ T $的上界, 但下界的估计却更加符合实际情况.随着Payne[3]开创性的工作以后, 这种讨论爆破时间下界估计的研究得到了许多学者的关注.需要指出的是, Li[4]在齐次Dirichlet边界条件下讨论了方程
的爆破问题, 获得了方程的解发生爆破时爆破时间的下界估计.这里拓展文[4]的研究范围, 将在Robin边界条件
下讨论相同的问题.在这种边界条件下, 文献[4]使用的Sobolev不等式不再可以使用.另外, 由于方程中梯度项的存在, 需要克服一些困难.因此, 当方程的解发生爆破时, 通过建立合适的Sobolev型微分不等式, 获得了爆破时间的下界估计.关于爆破时间下界估计的更多结论见文献[5-11].
设$ \Omega \subset \mathbb{R}^3 $是具有足够光滑边界$ \partial \Omega $与凸性的有界区域.考虑如下具有非线性梯度项的抛物问题
其中$ p, k > 0, q > 2 $, $ \Delta $和$ \nabla $表示拉普拉斯和梯度算子, $ T $是可能的爆破时间, $ (\partial u/\partial \nu) $表示在边界向外的法向单位导数, $ u_0(x) $满足适当紧条件的连续函数.
根据最大值原理有$ u\geq 0 $.另外, 除了在某些时刻会发生爆破之外, 还假设问题存在古典正解.当解发生爆破时, 本文的目的是要确定爆破时间的下界估计.为此, 定义下面的辅助函数
利用分部积分与边界条件, 简单的计算就可以得到
另外, 由于$ | \nabla u^{\frac{8p+1}{2}}|^2 = \frac{(8p+1)^2}{4} u^{8p-1} |\nabla u|^2, \ \quad | \nabla u^{\frac{8p+q}{q}}|^q = (\frac{8p+q}{q})^q u^{8p} |\nabla u|^q. $因此(2.3)式可以写成
下面考虑$ J_1 $.由于$ |\nabla u^{\frac{8p+q}{q}}| $的指数是$ q $而不是$ 2 $, 因此需要一点技巧去克服这个困难使得其指数变成$ 2 $.为了方便, 令$ v = u^{\frac{8p+q}{q}} $并通过Hölder不等式就有
再从Poincare不等式, $ \lambda_1 \int_\Omega v^q \mathrm{d}x \leq \int_\Omega |\nabla v^{\frac{q}{2}}|^2 \mathrm{d}x - \int_{\partial \Omega} v^{\frac{q}{2}} \frac{\partial v^{\frac{q}{2}} }{\partial \nu} \mathrm{d}S, $这里$ \lambda_1 $是下列弹性膜问题
的第一正特征值.利用Robin边界条件, 又有
再根据散度定理
利用$ \Omega $的凸性可以令
以及Young不等式得到
这里$ \theta $是任意正的常数.这样从(2.8)式和(2.6)式就可以获得
如果$ \Omega $, 也就是$ \lambda_1 $满足
然后选择充分小的$ \theta $, 使得
这样(2.9)式就可以写成:
其中$ C_0 = 1+\frac{kd(8p+q)}{2 l_0 \theta} $.把(2.11)式代入(2.5)式得到
因此根据$ J_1 $的定义就有
下面考虑(2.13)式右边的第一项.如果$ p > 1 $, 利用Hölder不等式和$ a^r b^q \leq ra + qb, r+q = 1, a, b > 0 $得到
这里
借助文献[3]的思想, 下面寻找一个合适的Sobolev型不等式来计算$ u^{9(8p+q)/8} $的积分.首先, 令$ x_{im} $和$ x_{iM} $分别表示$ \Omega $在坐标$ x_i $上的最小值与最大值, $ D_z $表示$ \Omega $与平面$ x_3 = z $相交的截面.其次, 为了方面计算, 设$ \omega: = u^{(8p+q)/4} $.根据Schwarz不等式得到
令$ P = (\bar{x}_1, \bar{x}_2, z) $是截面$ D_z $上的一个点, $ P_1 $与$ P_2 $表示截面$ D_z $的边界与直线$ x_2 = \bar{x}_2 $的两个交点.同样, $ Q_1 $与$ Q_2 $表示$ D_z $的边界与直线$ x_1 = \bar{x}_1 $的两个交点, 那么
这样可以获得
类似的,
把上面两式左右相乘并在$ D_z $上积分, 即有
再代回(2.15)式得到
其中
最后, 为了计算$ \omega^3 $在$ D_z $上的积分, 用$ \Omega^+ $表示$ \Omega $在$ D_z $上的部分, 而相应的边界是$ \partial \Omega^+ $, 而下面部分与边界则分别用$ \Omega^- $和$ \partial \Omega^- $表示.这样根据散度定理,
由此得到
即可以扩展(2.19)式到$ i = 3 $的情况.又由于
因此从(2.18)式得
其中$ A_i = \frac{1}{2} \int_{\partial\Omega} \omega^3 |\nu_i| \mathrm{d}S, B_i = \frac{3}{2}\int_\Omega \omega^2 |\omega_{x_i}| \mathrm{d}x, i = 1, 2, 3. $
另外根据Schwarz不等式, 有
类似地,
因此(2.20)式可以变成
类似于(2.7)与(2.8)式, 可以得到
把(2.22)式代入(2.21)式,
再使用不等式$ (a+b)^{\frac{3}{2}} \leq 2^{\frac{1}{2}}(a^{\frac{3}{2}} + b^{\frac{3}{2}}), a, b > 0, $ (2.23)式就可以变成
再依次使用Young和Hölder不等式, 还有
这里$ \varepsilon_1 $是任意正的常数, 并且还假设$ p > \frac{3q-4}{8} $.
结合(2.13), (2.14), (2.24)–(2.27)式, 并选择$ \varepsilon_1 $满足
由此就可以获得
下面估算$ J_2 $.类似于(2.14)式, 首先有
然后从$ J_2 $的定义得到
类似上面的讨论, 也可以估算$ u^{\frac{9(8p+1)}{8}} $的积分,
这里$ \varepsilon_2 $也是一个任意正的常数.把(2.30)式代进(2.29)式, 并且选择$ \varepsilon_2 $满足
(2.29)式就变成
把(2.28)和(2.31)式代入(2.4)式, 就可以获得关于函数$ \varphi $的一阶微分不等式
再积分得到
由此得到爆破时间$ T $的一个下界估计
其中$ \varphi(0) = \int_\Omega [u_0(x)]^{8p+1} \mathrm{d}x. $
把上面的分析总结成如下定理.
定理1 如果$ p > \max\{ 1, \ \frac{3q-4}{8}\} $, 以及(2.10)式成立, 那么当问题(2.1)的非负古典解$ u $在(2.2)式测度$ \varphi $的意义下有限时刻$ T $发生爆破时, 则$ T $的下界$ T_0 $由(2.35)式给出.
这里考虑$ p < 1 $的情况, 此时问题(2.1)的解不会在有限时刻内发生爆破.令$ \varphi(t) $依然是(2.2)式定义的辅助函数, 从(2.4)式得到
类似(2.11)式的讨论, 有
因此从(2.13)式得到
另外, 根据Hölder不等式, 也有
和
把(3.4), (3.5)式代入(3.3)式得
下面估算$ J_2 $.类似(2.11)式的讨论也得到
利用(3.4)式与(3.7)式有
最后, 结合(3.1), (3.6)与(3.8)式就可以得到下面的一阶微分不等式
显然, 如果问题(2.1)的解$ u $在测度$ \varphi $的意义下发生爆破, 那么从(3.9)式导出$ \varphi^\prime(t) \leq 0 $, 这是一个矛盾.这样就得到了下面的定理.
定理2 假设$ u $是问题(2.1)的一个非负古典解, 如果$ p < 1 $和第一特征值$ \lambda_1 $满足$ 2 l_0 \lambda_1 - 3k(16p+q) > 0 $, 那么$ u $不能在测度$ \varphi $的意义下发生爆破.
注 为了给出定理2的证明, 可以使用另外的辅助函数来代替(2.2)式定义的函数$ \varphi(t) $.例如, 类似文献[6], 可以定义函数$ \Phi(t) = \int_\Omega u^2 \mathrm{d}x $来完成定理的证明.