偏微分方程-浙大

《偏微分方程-浙大》由会员分享，可在线阅读，更多相关《偏微分方程-浙大（33页珍藏版）》请在金锄头文库上搜索。

1、偏微分方程 浙江大学 2008 9 2008 12 主要内容 一 一阶方程 二 二阶方程 椭圆型方程 典型的例子 Laplace方程 抛物型方程 典型的例子 热传导方程 Li Yau s Harnack inequality 双曲型方程 典型的例子 波动方程 三 特征流形的Cauchy问题 Cauchy Kowalevski定理 四 无解的线性方程 H Lewy例子 习题 必做 思考题 Open problems 参考书 F John Partial Diff erential Equations Springer Verlag 1982 谷超豪 李大潜等人 数学物理方程 高教出版社 2002

2、年 姜礼尚 孔德兴等人 应用偏微分方程 高教出版社 2008年 1 绪 言 1 基本概念 偏微分方程 PDE 关于函数u x y 的PDE是形如 F x y u ux uy uxx uxy 0 1 的关系式 其中F是自变量x y 未知函数u以及u的 有 限多个偏导数的已知函数 解 solution 称u是 1 的 解 如果把u x y 及其相应的偏导数代入 1 式后 在x y 空间的 某个区域 中 1 式关于这些变量恒等地成立 除非有相反的说明 在本课程中我们总是要求x y 是实的 u以及在方程 1 中 出现的u的偏导数在实空间的区域 中都是关于x y 的连续函数 为了简单起见 我 们有时也常

3、常省略区域 的明确描述 而把所述的命题也 局部地 适用于x y 空 间中一点的某一适当领域 偏微分方程组 PDEs 涉及一个或几个未知函数及其偏导数的多个偏微分方程组成一个方程组 记n为未知函数的个数 m为PDE的个数 当n m时 此时方程组称为欠定的 under determined 当n 0常数 5 其中u u t x1 xn n 1 弦的振动 波在管中的传播波 c表示传播速度 n 2 浅水面上的水波 n 3 声波或光波 例3 Maxwell 方程 Maxwell equations 在真空中且无自由电荷和电流的情况下 关于电场强度向量E E E1 E2 E3 及磁 场强度向量H H1 H

4、2 H3 的Maxwell方 程 实质上是由六个一阶方程所组成的线性方 程组 Et curlH Ht curlE divE divH 0 6 其中 是描述电磁介质的常数 分别称为真空介电常数和导磁系数 特别地 作为关 系式 Et curlH Ht curlE 的推论 如果t 0时 关系式 divE divH 0 成立 则上式对所有的t均成立 不难验证 这里的每个分量Ei Hk均满足具有c2 1 的波动方程 5 事实上 在方程组中消去磁场强度便得到电场强度的偏微分方 程 这只需对 6 式中的第二式求旋度 再用其第四式可得 curl curlE curlH t Ett 4 又因为 curl cur

5、lE divE E 再利用 6 式中的第三式得到 Ett 1 E 类似地 我们可以得到磁场强度向量H所满足的偏微分方程 例4 在经典弹性理论中 弹 性 波 可由线性方程组 2ui t2 ui xi divu i 1 2 3 7 描述 其中ui t x1 x2 x3 是位移向量u的分量 是密度 而 是弹性材料的Lame常 数 可以证明 每一个分量ui都满足由两个不同的波动算子所组成的四阶方程 2 t2 2 2 t2 ui 0 8 当弹 性 平衡 即ut 0 时 我们便得到重 调 和 方 程 2u 0 9 例5 当密度和比热都是常数时 导热体中的温度分布满足热 传 导 方 程 ut k u 10

6、其中k 0是常数 表示介质的热传导系数 例6 在势能为V x y z 的场中 运动的质量为m的单个质点所满足的Schr odinger方 程是 i t 2 2m V 11 其中h 2 是Planck常数 例7 Tricomi方程 uxx xuyy 12 另一类重要的方程是 uxx yuyy 5 这两类方程在平面跨音速流的研究中具有十分重要的作用 上述例子中的方程都是 线 性 的 非 线 性方程也是常见的 但是求解非线性方程实际 上更为困难 因此在实际中常用线性方程近似的表示它们 下面是几个非线性方程的 例子 例8 3维Euclid空间中的极 小 曲 面 z u x y 即通过给定周线而具有最小

7、面积的曲 面 满足下述二阶拟线性方程 1 u2 y uxx 2uxuyuxy 1 u 2 x uyy 0 13 例9 1 n维Minkowski空间中的极 值 曲面x x t Rn满足下述二阶拟线性方程 x 2xtt 2hxt x ixt xt 2 1 x 0 14 例10 密度为 的二维定态绝热无旋等熵流的速 度 势 x y 其速度分量为 x y 满 足下述二阶拟线性方程 1 c 2 2 x xx 2c 2 x y xy 1 c 2 2 y yy 0 15 其中c是速率q p 2 x 2y的已知函数 例如 对于状态方程为 p A 16 的多方气体 或称为 气体 c2 1 1 2 q2 17

8、例11 关于不可压缩液体的粘性流的Navier Stokes方 程 是速度分量uk和压力p之间的 一组偏微分方程 ui t 3 X k 1 ui xk uk 1 p xi ui i 1 2 3 3 X k 1 uk xk 0 或写成 divu 0 18 其中 是常密度而 是运动的粘性系数 6 例12 空气动力学方程 t 3 X j 1 xj vj 0 t vi 3 X j 1 xj vivj ijp 0 t E 3 X j 1 xj vjE pvj 0 19 其中 t x 表示气体的密度 v v1 t x v2 t x v3 t x 为速度 p为压力 E E t x 为内能 由热力学知识可知

9、所有的热力学量只有两个是相互独立的 因此密 度 压力p 温度T以及内能E中 它们之间有一个确定的关系式 p p E 或 p p T 20 对不同的气体 上式具有不同的表达式 20 式通常称为气体的 状 态 方 程 注意到状态 方程 20 式 19 式构成一个封闭的一阶拟线性偏微分方程组 例13 函数u t x 的三阶非线性方程的一个典型例子是Korteweg de Vries方 程 ut cuux uxxx 0 21 它是在水波的研究中被首先提出的 可用来描述浅水波中的孤立子的传播 例14 双曲Monge Ap ere方程 S S2 1 S S 22 上述方程是在研究平均曲率流时提出的 通常我

10、们试图描述或理解所考虑的偏微分方程的解 流 形 不同类型的方程 其结果 是十分不同的 偏微分方程的有意义的 适 定 问 题常常受具体的物理背景及物理意义所 启示 7 第一章一阶方程 一阶方程是一类最基本的偏微分方程 它在数学 物理学以及工程技术中具 有广泛的应用背景 本章我们通过具体的例子着重介绍一阶方程的一些基本概 念 方法和结果 1 一个简单线性方程 在本节中我们考虑关于函数u u t x 的一个最简单的方程 ut cux 0 1 1 并用它来说明将在后面起重要作用的某些概念 其中c 0是常数 一 解析求解 特征线方法 在 t x 平面上 定义 特 征 线 族 dx dt c 1 2 沿着

11、特征线族中的任意一条直线 x ct const 1 3 方程 1 1 的解u满足 du dt d dtu t ct ut cux 0 1 4 因此 沿着这样的一条直线 u保持为常数 它仅与区分线族中不同直线的参数 有关 于是方程 1 1 的 通 解具有下面的形式 u t x u 0 f f x ct 1 5 其中f 是一任意给定的函数 它表示u的初 始 值 上式表明通解u由初始值 u 0 x f x 1 6 唯一确定 反之 如果f属于C1 R 类 则任何形如 1 5 的函数必是 1 1 的具 有初始值f的解 我们注意到u在任意点 t x 处的值仅与初始函数f在单个变 1 元 x ct时的值有关

12、 而这个 值正是过 x t 点的特征线与初始曲线x 轴交 点的横坐标 u t x 关于初始值的依 赖 区 域 由单个点 组成 在点 处的初始值 只影 响 特征线 1 3 上解u t x 的值 见图1 1 6 t t x x 0 x ct 图 1 1 特征线 如果对每一个固定的时间t 用 x u 平面上的图象表示函数u 我们会发 现t T时的图象是由t 0时的图象沿平行于x 轴的方向平移cT的距离得到 u x 0 u x cT T f x 1 7 解的图象描述了一个以速度c向右传播而不改变形状的波 见图1 2 6 u c xxx cT u 0 x u T x 图 1 2 波的平移 波 方程的解

13、线性波 线性方程的解 非线性波 非线性方程的解 二 近似求解 有限差分方法 2 下面我们用这个具有显示解的例子引出用有 限 差 分 方 法 求偏微分方程数值解 的某些概念 用x方向步长为h t方向步长为k的矩形网格覆盖 t x 平面 换句话说 即只 考虑x是h的整数倍 t是k的整数倍的点 t x 为了求数值近似 用差分方程 v t k x v t x k cv t x h v t x h 0 1 8 代替偏微分方程 1 1 是自然的 当h k 0时 这个方程在形式上就变成 vt cvx 0 我们自然要问 当h k很小时 在网格点上 1 8 的具有初始值 v 0 x f x 1 9 的解v与初值

14、问题 1 1 1 6 的解相差多少 为了回答这个问题 我们引入 k h 这样 1 8 式可写成递推公式 v t k x 1 c v t x cv t x h 1 10 上式表明我们可以用t时刻v的值表示t k时刻v的值 引入平 移 算 子 E Ef x f x h 1 11 于是 1 10 式就变成 v t k x 1 c cE v t x 1 12 其中t nk 由此 通过迭代就得到关于 1 8 的初值问题的解 v t x v nk x 1 c cE nv 0 x n X m 0 Cm n 1 c m cE n mf x n X m 0 Cm n 1 c m c n mf x n m h 1

15、 13 3 显然v t x v nk x 的依赖区域是由x 轴上的点 x x h x 2h x nh x t 1 14 组成 所有这些点均在x和x nh之间 而偏微分方程解的依赖区域是由单 点 x ct x c nh组成 它在区间 x x nh 之外 当h k 0时 我们显 然不能指望v收敛到偏微分方程的精确解 因为在构造v t x 时我们没有用到对确 定u t x 至关重要的值f 反而用到了与u t x 毫无关系的f在区间 x x t 1 上的值 这个差分格式违背了Courant Friedrichs Lewy准 则 差 分 方 程 的 依 赖 区 域 的 极 限 必 须 包 含 偏 微 分

16、 方 程 的 依 赖 区 域 差分格式 1 8 的严重不 稳 定 性 也说明它是不适当的 在应用问题中数据f不 是完全精确的 而且在数值计算的每一步我们也不易 甚至不可能 用精确值 而总会有小的误差 从 1 13 式显然可以看出 如果f有适当交错符号的误差并 且f的绝对值不超过 那么其结果将可能导致v t x v nk x 有大小为 n X m 0 Cm n 1 c m c n m 1 2 c n 1 15 的误差 因此 当网格比 固定时 v的可能产生的误差将按t方向的步数n的指数 函数增长 一个更为适当的差分格式是用向后差商写出的 v t k x v t x k cv t x v t x h h 0 1 16 或者用记号写成 v t k x 1 c cE 1 v t x 1 17 1 17 的初值问题的解变成 v t x v nk x n X m 0 Cm n 1 c m c n mf x n m h 1 18 在这个格式中 v t x 关于f的依赖区域由点 x x h x 2h x nh x t 1 19 4 所组成 所有这些点均在x t 1和x之间 令h k 0 而网格比 保持不