差分格式学习报告.docxVIP

1 一阶线性常系数双曲型方程 先考虑线性常系数方程 TOC \o 1-5 \h \z \o Current Document —+ 白竺=0,X G R,t 0 （1.1） dt dx 其中a是常数。这是最简单的双曲型方程，一般称其为对流方程。对于方程（1.1） 给定初始条件 \o Current Document u （x,0） = o（x）, x g R （1.2） 在此之前，已经讨论过（1.1），（1.2）式的解，其解沿方程（1.1）的特征线 x 一 at = （1.3） 是常数，并且可以表示为 u （x, t） = u （g） = u （x - at） 下面讨论双曲性方程的应风格式，Lax-Friedrichs格式，Lax-wendroff格 式，Courant-Friedrichs-Lewy条件利用偏微分方程的特征线来构造有限差分格 式，蛙跳格式，数值例子。 1.1迎风格式 迎风格式的基本思想是简单的，就是在双曲型方程中关于空间偏导数用在特 征线方向一侧的单边差商来代替，（1.1）式的迎风格式为 TOC \o 1-5 \h \z \o Current Document Uj+1 - Uj + aUj - Uj-i = 0, a 0 （1.4） t h \o Current Document Uj+1 - Uj + aUj+i - Uj = 0, a 0 （1.5） t h 下面考察（1.4）式的相容性和稳定性： 1、相容性:Un+1-Un K T+ 1 竺 T 2 + 1 竺 T 3 + … j j dt 2! dt2 3! dt3 1、相容性: Un+1 - Un d U 1 d 2 U 1 d 3U —j j = + T + T 2 + t dt 2! d12 3! d13 dU ] 1 d2u 1 d3u ] Un — Un = h — h 2 + h3 + j j-1 dx 2! dx2 3! dt3 n /「j = a坐-巴竺h +竺竺h2 + 0^4) ② h dx 2! dx2 3! dx3 ① ①+②得: un un+1 — Un Un — Un 6 u t h 6 t J_ 竺 t + 上竺 t 2 + 0Y 3 a 辿 h + a 6uh2 + 0=4 2 ! 6x 2 3 ! 6x 3 =空+。竺)+ 1竺t —a竺h + a竺h2+… \6t 6x) 2! 6t2 2! 6x2 3! 6x3 所以 T(x ,t ) = 1竺t — a竺h +… j n 2! 6t2 2! 6x2 截断误差为。(+ h),故迎风格式对t 一阶精度，对h 一阶精度，当t T 0,h T 0时 T（x.,t ） T 0 ,即迎风格式相容。 2、稳定性 先把差分格式改写为 u.+i = u. - a」u.—Uj\ 令 u n = vne ikjh， j将其代入上式得到消去公因子有v 令 u n = vne ikjh， j 将其代入上式得到 消去公因子有 v n+1 e ikjh = v n e ikjh 一 a 人 一 ak + ake—ik)vn neikjh — yneik(j-1)h ) v n+1 则增长因子为 G(t ,k)= 1 — a人 + ake—ikh = 1 + a人 —ikh —]) =1 — ak —cos kh) a k i sin kh 所以有=1 — 所以有 =1 — 2ak —cos kh ) + a 2k 2 —cos kh )2 + a 2 k 2 sin 2 kh \G(t，k)2 = 1 — ak G — coskh)]2 + a2k 2 sin 2 kh =1 — 4ak—ak ) =1 — 4ak 2 即当ak 1时，差分格式（1.4）是稳定的，同样地，当|a|k 即当ak 1时， 稳定的。 类似地，用Fourier方法讨论差分格式: Un+1 — Un Un — Un J r卜 + a h_l ，a (1.6) Un+1 — Un Un — Un ^T^ + = 0，a v 0 (1.7) 对于（1.6）式，其增长因子为G（r,k）= 1 + ak-akeikh,由此有 \G(: ,k)2 = 1 + ak(1 — coskh)]2 + a2k2 sin2 kh =1 + 4ak(1 + ak) sin2 kh 2 取sin kh。0，所以G（，k） 1，从而破坏了 von Neumann条件 （1.6 ）是绝对不稳定的，同样可证，差分格式（1.7）也是绝对不稳定的。 因此差分格式 对比发现差分格式（1.4）与（1.7）在形式上是一致的，但因a的符号不同， 导致稳定性不同，从而我们可知，如果差分格式与微分方程的特征线走向一致，