第九章 非线性偏微分方程

前面几章索研究的偏微分方程都是线性的，但在实际工程级数及自然科学中索遇到的方程大多都是非线性的，在有些情况下，人们为了研究方便，对问题补充了一些附加的条件或略去一些次要的项，才得到线性方程。在这一章内，我们将从一个具体问题出发引入非线性偏微分方程的概念，然后重点讨论两类重要的非线性方程。

§9.1 极小曲面问题

在第八章内已经说过，求解一个边值问题可以转化成求它所对应的一个泛函的最小值（当然，一般说来变分问题的解只是原边值问题的弱解）。其实，在数学里也已证明了相反的结论，即在一定条件下一个变分问题的解必满足一个微分方程。在这一节内，我们以极小曲面问题为例说明这个事实。

设是平面上有界区域，它的边界是充分光滑的，其方程为：

其中即是一条闭曲线。现在在上给定一条空间曲线（即作一条空间曲线，使它到所在平面的投影为）：

              （9.1）

这里。所谓极小曲面问题就是要确定一张定义在上的曲面，使得

（1）以为周界；

（2）的表面积在所有以为周界的曲面中是最小的。

假定空间曲面的方程为

则由微积分学可知，这个曲面的表面积为

             （9.2）

于是上述极小曲面问题就变成求一个函数，使得

（1）由所表示的曲面以为周界，即

，或者说，，

其中由（8.7）给出；

   （2）                                   （9.3）

这是一个变分问题。

如何求出变分问题（9.3）的解？我们先来看看假若是（9.3）的解，那么必需满足什么样的条件。为此，在任取一个元素，即任取，即。对任意，记

                （9.4）

其中由（9.2）确定，从（9.2）可知是定义在上的一个可微函数，由于是（9.3）的解，所以对任意处取得最小值，故

                 （9.5）

不难看出

代入（9.5）得

假若具有更好的光滑性，例如，由格林公式可得

由于，即，因此上式左端第二项为零，再由的任意性及被积函数的连续性可知

         （9.6）

这个方程称为变分问题（9.3）的Euler方程。

上面的推导说明，如果是（9.3）的解，且，则必满足（9.6），当然还满足边界条件

                （9.7）

因此定义在上且以空间曲线为周界的极小曲面必定在内适合方程（9.6）和在上满足边界条件（9.7）。方程（9.6）可以改写成

       （9.8）

这个方程通常称为极小曲面方程。它有什么特点？它关于二阶导数及是线性的，但它们前面的系数分别含有，及，所以对，来说它不是线性关系，特别是，如果把，，及同等看待，这个方程对它们不少线性方程，故它是一个非线性方程。

§9.2 非线性偏微分方程举例

     在上一节内，我们以极小曲面问题为例得到了一个非线性偏微分方程（9.8），其实，在力学、物理学及几何学中都有大量的非线性偏微分方程。例如，在热传导问题（第一章第一节例4）中，如果热传导系数不是常数，而是温度的函数，则三维热传导方程为

       （9.9）

这也是一个非线性方程。

在流体力学中，描述粘性气体运动的方程是著名的Navier-Stokes方程，其形式为

              （连续性方程）

            （动量方程）

（能量方程）

其中

          （9.11）

这里是流体密度，是流速，是温度，、是粘性系数，是传热系数，是压强，是定压比热，是气体常数，表示粘滞力的张量，是Kronecker记号，即

当流体为不可压缩时，是长和数；又若不计温度的变化，这（9.10）化为不可压缩流体的Navier-Stockes方程

         （9.12）

取，并利用（9.11），这上述方程组为

    （9.13）

这是关于的非线性方程组。

在热平衡问题中，如果热传导系数是常数，但物体内含有一个依赖于温度及温度梯度的热源，则可得

               （9.14）

在微分几何中，若要求出中曲率为已知的曲面时，就需要求解下列方程：

           （9.15）

其中。这个方程称为蒙日－安培尔（Monge-Ampere）方程。

上面我们已经从不同的问题引入了一些非线性方程或方程组，现在再对它们作一些比较。方程（9.14）中的最高阶导数（即二阶导数）部分纯粹是线性得，它的非线性只出现在函数及其一阶导数项，这样的方程称为半线性方程，方程组（9.13）也是半线性的；方程（9.8）对最高阶导数（二阶导数项）来说是线性的，但它们的系数依赖于未知函数的非最高阶导数（那里是一阶导数），这样的方程称为拟线性的；方程（9.15）的特点是对最高阶导数（二阶导数）也是非线性的，这样的方程称为完全非线性（或真正非线性）方程。显而易见，完全非线性方程的非线性程度最高，半线性方程的非线性程度最低，拟线性方程的非线性程度介于两者之间。

对于非线性偏微分方程，一般说来是无法求出解的表达式，只能求其近似解。但对一些很特殊的情形，通过适当的未知函数的变换将方程化成线性方程，或者经过适当的数学处理化成可以求解的方程，下面举例说明。

例1 在流体力学中有一个很重要的比尔吉斯（Burgers）方程

                 （9.16）

这是一个二阶偏微分方程，为了解这个方程，令，再对积分一次可得

再令

则得

这是一维的线性热传导方程，对它的各种定解问题可以用第二、三中的方法求出它的解，有了之后可以求出。

例2 在微分几何中遇到如下Liouville方程

                    （9.17）

这是一个半线性的二阶方程。若令是

                   （9.18）

的解，再构造一个偏微分方程组

             （9.19）

其中是常数。通过计算可以验证：若是（9.19）的解，则必是（9.17）的解。有第三章已知，（9.18）的解总可以写成

其中是任意可微函数。有了再解（9.19），最后可得（9.17）的通解

与线性方程相比，非线性方程还有一个特点，即它的解即使存在，也不一定对所有的时间都存在（当然假定方程中含有时间变量），而只是再某个有限时间内存在，见下例

例3 考虑Riccati方程的初值问题

            （8.20）

容易求出它的解

显然，若，则（9.20）的解对所有都存在，简称整体存在；若，则当时，，这时解在时刻产生破裂，所以（9.20）只在内有解，简称解是局部存在的。

§9.3 单个守恒律 激波

在这一节内，我们将研究形如

                    （9.21）

的一阶非线性双曲型方程初值问题。由于方程（9.21）左端是散度形式，通常将它称为守恒律。我们将要指出这个方程的解可能产生间断，并着重介绍激波的概念。

先看一个特例，即考虑Burgers方程的初值问题

在这个问题古典解（类解）存在的范围内，可由

                    （9.24）

定义其特征线。显然，沿着特征线有

即在每一条特征线上取常数值。由（9.24）知特征线是直线，通过点的特征线为

                 （9.25）

在其上

                （9.26）

设的模有界（即），且较小，则有

故由（9.25）可得

代入（9.26）即得问题（9.22），（9.23）的解为

这说明，这个问题总存在惟一的局部古典解。

但是，只要不是一个单调不减函数，在轴上必存在两点，使得：

但

这时，过此二点的特征线

必在有限时刻相交，在交点处解就不能惟一地确定，即初值问题不可能在内存在整体的古典解。这种现象在力学里对应于形成激波。

上述讨论可以用于一般的单个守恒律（9.21），考虑这个方程具有初始条件（9.23）的初值问题。和前面类似，如果存在点，使得：

则从点，出发的特征线与将在内某点相交，沿。因此，在点解必须间断，即使初值与函数充分光滑，甚至是解析的，仍然如此，即这个现象完全是由于方程的非线性所致。

为了把这个问题说得更清楚一点，设满足凸性条件（图9.1）

                     （9.27）

在内任取一点，以表示通过的特征线与轴的交点的横坐标。因为沿着特征线是常数，且

可见必须由下列隐式关系来确定

             （9.28）

若是可微的，利用隐函数存在定理，对充分小的，由（9.28）解得，且

因此，若，则对所有，和保持有界，且解是整体存在的；如果在某点，则当时，与变成无界的，即古典解不可能整体存在。

上面的分析表明，为了研究初值问题（9.21），（9.23），必须推广解的概念。

为了推广解的概念，令。在内任取一个函数，使它当及充分大时恒等于零，用乘方程（9.21）的两端后再在内积分，利用分部积分法及的条件可得

         （9.29）

我们把满足（9.29）的函数称为问题（9.21），（9.23）的广（或弱解）。显然，若是古典解，则它必是广。

值得指出的是，并不是每个间断函数都可以作为问题（9.21），（9.23）的广，表达式（9.29）对解的间断性强加可苛刻的，为了说明这一点，设是一光滑曲线，在其上由阶跃间断，但除以外，是光滑的。在上任取一点，以表示以为圆心的小球（图9.2），假设在内由方程表示。任给，由（9.29）得：

（9.30）

由于在内是的，利用散度定理得

     （9.31）

因在上，这些线积分仅在上是非零的。记，则

因此

          （9.32）

其中

      （9.33）

分别是与沿着的阶跃。因为任意，故

             （9.34）

其中为间断线的速度。

公式（9.34）反映了与在上的阶跃与的速度之间的关系，这个关系称为阶跃关系，在气体动力学总就是熟知的Rankine-Hugoniot条件。

对于广而言，有可能丧失惟一性，例如Burgers方程（9.22），如果初值取为：

容易验证上述初值问题有两个解

这两个解在扇形内完全不同，有趣的是，是连续函数。这说明连续解可能具有部连续的初值，这一点与线性方程也完全不同。

既然广不一定是惟一的，那么如何区分物理上说需要的解？可以证明单个守恒律方程（9.21）存在惟一的满足下列所谓熵条件

            （9.35）

的解，其中不依赖于及。

从条件（9.27）及（9.35）可知，当增大时解的阶跃是下降的，及，如果我们将阶跃条件（9.34）写成

其中，则在凸性条件（9.27）之下，得到下列熵不等式

                （9.36）

这表明间断线的速度介于介于间断线两侧的特征速度（注意，守恒律（9.21）的特征线由确定）之间。

满足条件（9.34）与（9.36）的间断线称为激波，也称为激波的速度。

§9.4 KdV方程 孤立子

15年脱维克（Korleweg）和德伏莱斯（de Vries）在研究水波时导出了如下形式的半线性方程（简称KdVfangcheng）：

                （9.37）

不失一般性，不妨设（事实上，通过变换

总可以将情形化成的情形）。

现在来寻求方程（9.37）的平面前进波（或简称行波）解，令

，是常数，       （9.38）

将（9.38）代入（9.37）得

对积分一次得

，为任意常数         （9.39）

用乘（9.39）并对积分，得

         （9.40）

或

           （9.41）

如果将与分别理解为时间坐标与空间坐标，则（9.41）能视为举单位质量的质点在势的作用下运动的能量关系。

从方程（9.41）可以看出，如果，则是KdV方程的解，这个解称为定常解；只有当时，KdV方程才可能有实的行波解。由于是的三次方程，不外有两种可能性：

（ⅰ）仅有一个实的零点（如图9.3所示），这时为了获得实的行波解，只能在内考虑，在这个区域内

                 （9.42）

假设（否则不可能是的单重实零点），则（9.42）存在满足

的解。由于当时，，故这个解是无界的。

（ⅱ）有三个实的零点

不妨设，这时可以表示成

            （9.43）

我们比较（9.40）与（9.43）可得

再分三种情况：

（A）；（B）；（C）.

对应于这三种情况的草图如图9.4，其中曲线A与轴有三个交点，曲线B与轴在点相切，曲线C与轴在点相切。为了使解是实的且有界。只能在下列区域内考虑（9.42）

对于曲线A：；

对于曲线B：；

（对于曲线C：，这时与（i）相同，解是无界的）

在情形A，容易求出

         （9.44）

或

  （9.45）

其中是Jacobi椭圆函数，

             （9.46）

由（9.45）可以看出关于是周期的，其周期为

    （9.47）

其中是第一类完全椭圆积分，故在情形A，KdV方程的有界行波解是一个周期性的波。

在情形B

              （9.48）

这时有

              （9.49）

这里因时，，故记，且

故

由（9.49）或（9.50）所表示的行波有什么特点呢？它在传播种波形不变，基本上集中在的周围，在时的定常状态与时的定常状态相同，都是，其形状如图9.5所示。

这种波还有一个特征：如果有两个这样的波，一个波速快，一个波速满，开始时波速快的波在后面，波速慢的波在前面，后来被赶上了，经过相互作用后，过一段时间后又分开为原来形式的两个波，波速快的走在前面了。这种相互作用受干扰后仍恢复原状运动的现象引起了人们很大的兴趣，称这种波为孤立波或孤立子。