(四) 高精度差分格式及其数值解的逼近程度分析
差分格式
§1. 差分 1．一阶导数的差分近似（差商）导数的定义： ()()()0000limx x f x f x f x x x ®-¢=-导数的近似：()()()10010f x f x f x x x -¢»- （当 1x 与 0x 足够接近时）这样的表达式称为差商，它可作为导数的近似，称为导数的差分近似。
误差分析 － 泰勒展开：将 ()1f x 在 0x 处做泰勒展开，有()()()()()()21001001012f x f x f x x x f x x x ⅱ?=+-+-+L于是()()()()1001010f x f x f x x x x x -¢-=--各种差分近似： 取 0h >（称为步长），则可以有 向前差分近似（相当于取 100x x h x =+>）()()()000f x h f x f x h+-¢»向后差分近似（相当于取 100x x h x =-<）()()()000f x f x h f x h--¢»中心差分近似 （前差近似与后差近似的算术平均）()()()0002f x h f x h f x h+--¢»2． 差分近似的一般形式差分近似的一般形式可写成()()()()()()()()022********* m m n n f x c f x c f x c f x hc f x c f x c f x c f x ------é¢?++êë+ù++++úûL L或简写为()()01nj j j mf x c f x h =-¢»å 称为一阶导数 ()0f x ¢ 的一个 1m n ++ 点差分近似。
这里0 ( , , 2 , 1 , 0 , 1 , 2 , , )j x x jh j m n =+=---L L差分近似的精度 ： 阶 定义：若()()()01n pj jj mf x c f x h h =-¢-=å 则称表达式 ()1nj j j mc f x h =-å 是一阶导数 ()0f x ¢ 的 p 阶差分近似。
解高维抛物型方程的一个高精度显式差分格式
Ｖ０ １ ． ３ １ Ｎ ｏ ． ２
Ｍａ ｒ ． ２ Ｏ１ ３
文章编号 ： １ ０ ０ ８—１ ４ ０ ２（ ２ ０ １ ３ ） Ｄ ２— ０ ２ ７ ５一 ｏ ４
解 高维 抛 物 型 方 程 的 一个 高精 度 显 式 差分 格式①
沈高峰 ， 谷淑敏
（ １ ． 郑州轻工业学院计算 机与通 信工程学院 ． 河南 郑州 ４ ５ ０ ０ ０ ２； ２ ． 中原工学院信 息商务学院基础部 。 河南 郑州 ４ ５ ０ ０ ０ ７ ）
摘 要： 构造和研 究了五雏抛物型方程的高精度显式差分格式． 首先给 出了含参 变量的差分方 程， 并用待定 系数法适 当地选取 了这些参数的表 示式， 以使差分方程的截 断误差阶尽可能高地达 到了ｏ （ ｚ ￣ ｔ ＋△ ）； 其 次 用稳定 性分 析 的 Ｆ ｏ ｕ ｒ ｉ ｅ ｒ 方 法给 出了所得 格 式 的稳 定 性 条件 ； 接 着确 定 了高精度 显 式差分 格 式的稳 定性 条件 为 ｒ＜２ ／ ５； 最后 给 出了数 值例 子 ， 数值 结果 表 明 了本文格
件 ｒ＜１ ／ ６ 或 ｒ＜１ ／ ４ ． 本 文对 五 维热 导方 程构造 了
一
△ ［ ‰ ＋（ 叩 ＋叩 ： （ １ ◇）＋７ ７ ３口
＋７ ７ 田
＋ ｏ） ］＝（ ＋ ＋ ＋ ＋ ）
（ ７ ７ ６ ＋７ ７ ７ ｎ 胁 － １ ）＋［ （ １ ２ ◇ ＋３ （ 口
第３ １ 卷 第 ２期
２ ０ １ ３ 年 Ｏ ３月
佳 木 斯 大 学 学 报 （ 自 然 科 学 版 ） Ｊ ｏ ｕ ｎａ ｒ ｌ ｏ ｆ Ｊ ｉ ａ ｍｕ ｓ ｉ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ（ Ｎ ａ ｔ ｕ ｒ ａ ｌ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ Ｅ ｄ ｉ ｔ ｉ ｏ ｎ ）
微分方程中的数值解误差分析方法
微分方程中的数值解误差分析方法在数学领域中，微分方程是描述自然现象和物理现象的一个非常重要的工具。
然而，大多数微分方程很难用解析的方法求解，因此我们通常使用数值方法来近似求解。
然而，这些数值解不可避免地会引入误差。
本文将介绍微分方程中的数值解误差分析方法。
一、局部截断误差在使用数值方法求解微分方程时，我们通常会引入一个步长h。
在每个步长上，我们通过一系列迭代计算来逼近真实的解。
然而，由于近似计算和舍入误差等原因，我们得到的数值解与真实解之间存在误差。
这个误差被称为局部截断误差。
局部截断误差可以通过泰勒展开来近似计算。
假设我们使用的数值方法是Euler方法，那么可以得到如下的局部截断误差公式：$$LTE = \frac{y(t_{n+1}) - [y(t_n) + hf(t_n, y(t_n))]}{h}$$其中，$y(t_n)$是真实解在时间点$t_n$的值，$f(t_n, y(t_n))$是微分方程的右侧函数在$t_n$和$y(t_n)$处的取值。
二、全局截断误差除了局部截断误差之外，我们还需要考虑全局截断误差。
全局截断误差是指在整个求解过程中，数值解与真实解之间的误差累积情况。
通过对局部截断误差进行逐步累积，我们可以得到全局截断误差的估计。
例如，使用Euler方法求解微分方程，假设总共迭代了N步，步长为h，则全局截断误差的估计为：$$GTE = \frac{LTE}{h} \times N = \frac{y(T) - y(t_0)}{h} = O(h)$$其中，$y(T)$是真实解在求解区间的终点处的值，$y(t_0)$是真实解在求解区间的起点处的值。
三、稳定性分析除了局部截断误差和全局截断误差，稳定性也是数值解的一个重要性质。
在数值方法中，一个稳定的方法可以保证数值解不会因为舍入误差或者数值不稳定性而发散。
稳定性分析通常通过稳定性函数来进行判断。
对于一个给定的数值方法，我们可以将其误差传播到未来的时间点，然后观察误差是否会趋于无穷大。
五点差分格式
五点差分格式
五点差分格式是一种常用的数值计算方法，可以用于求解偏微分
方程的数值解。
它是一种前向差分格式，因为它利用了某一时刻的函
数值和其相邻的四个时刻的函数值进行计算。
五点差分格式在应用中
具有许多优点，比如计算精度高、计算效率高等。
五点差分格式的基本思想是将要求解的偏微分方程转化为差分方程，然后用迭代方法求解差分方程。
对于一个二阶偏微分方程，可以
先将其变为一组两个一阶偏微分方程，然后利用五点差分格式求解。
五点差分格式的优点在于它可以准确地估计函数在各个位置处的
导数，并且可以使用相对较大的时间步长进行计算。
而缺点则是它对
边界条件的要求比较高，需要额外的处理才能适应复杂的边界条件。
在实际应用中，五点差分格式被广泛应用于流体力学、地球物理学、材料科学等领域。
例如，它可以用于模拟流体运动过程中的压力、速度等物理量变化，以及地球内部的温度、密度等参数变化。
五点差
分格式也可以用于计算材料的力学性能，比如应力、应变等。
总之，五点差分格式是一种非常有用的数值计算方法，可以应用
于多个领域中的偏微分方程求解。
在实际应用中，需要对其进行适当
的修改和调整，以满足具体问题的要求。
【计算流体力学】第4讲-差分方法2
暗示
Re x 是某一特 征量
5
对流-扩散方程的特性：
u t
u x
1 Re
2u x2
（线性）差分方程：
u n1 j
ak
u
n j
k
(1)
k
j-2 j-1 j j+1 j+2
n+1
n
u n1 j
a un 2 j2
a un 1 j 1
......
a2
u
n j
2
某点的值是上一时刻周围几个点上值的线性组合
ui ui1 ui1 ui 则j点在间断右侧
j
（a>0 时， 右侧为“前”）
原理： 越靠近间断，振荡越剧烈
2） 根据GVC的思想构造格式
u a u 0 a 0 t x
xu j
(3u j
4u j1
u j2 ) /(2x)
(u j1 u j1) /(2x) when
when u j u j1 u j1 u j u j u j1 u j1 u j
B: 高阶人工粘性 人工粘性
系数
人工粘 性项
x2
x
N
U x
光滑区为二阶小量
分离流—— 对粘性敏感
Von Neumann MacCormack
N
u x
N
| u | c p
2 p x2
大梯度区，加大 人工粘性
转捩——对粘性敏感 8
Jameson 人工粘性法
二阶人工粘性
四阶人工粘性
u t
f (u) x
人工粘性法
TVD, 保单调限 制器 群速度控制格式
j=1
j=N n时刻： 单调函数
紧致差分格式
紧致差分格式
紧致差分格式是一种数值求解偏微分方程的方法，其主要特点是在离散化时使用了较少的节点，同时保持较高的精度。
在紧致差分格式中，我们将要求解的偏微分方程离散化为一个代数方程组，通过求解该方程组来得到数值解。
为了实现高精度，紧致差分格式通常会使用高阶的差分算子，例如二阶中心差分算子或者非中心差分算子。
常见的紧致差分格式包括：
1. 二阶中心差分格式：使用二阶中心差分算子来逼近偏微分方程中的导数项，从而得到一个二阶精度的差分格式。
2. 基于算子分裂的紧致差分格式：将整个偏微分方程分解为几个部分，在每个部分中采用不同的差分格式来逼近，然后通过交替迭代的方式求解。
3. 符号差分法：利用泰勒级数展开，将偏微分方程中的导数项用差分算子展开，然后通过合理的组合得到一个高精度的差分格式。
紧致差分格式一般适用于光滑的问题，并且需要在边界处进行一定程度的调整，以满足边界条件。
同时，紧致差分格式通常需要解一个线性方程组，因此对于大规模问题可能需要使用高效的求解算法。
四阶差分格式
四阶差分格式
u:+a4x=0
(1)其中u是以速度a传播的纯量，a是正实数我们知道，对该方程进行空间离散后，可以得到一个常微分方程组（Odes):=f(u,)Dot
再选择适当的时间离散格式，就可得到方程(1)的数值解为此，对所定义的空间区域进行一致网格剖分，设h为网格剖分步长，x=h,=u(x),为了近似计算空间变量的一阶导数，使其具有四阶精度，我们用待定系数法设
(△4)y=4h(4y+22)+2h(4+1-4-1)
其中（△w),表示对在结点x,处的近似将式中u及△u各项在w点做Taylor展开，合并后比较对应项系数，有：m1=4,m2=·1书，于是对空间变量我们得到了一个简单的四阶中心差分格式：=品1)
(2)12h(u j+2-42)再用蛙跳方法对离散得到方程()的全离散格式
(3)其中c=a k万是CFL条件数，k是时间剖分步长对问题(P):,+aux=0,(x,t)∈0,1X0,T]u(x,0)=sim(2Π)，x∈0,1
(4)period i c boundary value condition,t∈0，T]。
紧致差分格式
紧致差分格式
紧致差分格式是一种在数值计算和数值模拟中常用的数值解法。
它通过将连续的物理量分割成离散的点，并使用差分来近似导数，从而将求解微分方程的问题转化为求解代数方程的问题。
紧致差分格式的优势在于其高精度和较小的误差。
相比其他差分格式，紧致差分格式在相同离散点数的情况下能够提供更准确的解。
这是因为紧致差分格式通过使用更多的信息来近似导数，从而减小了离散误差。
紧致差分格式的核心是在相邻的离散点上使用高阶差分，以提高精度。
在一维情况下，一种常用的紧致差分格式是中心差分格式，它使用相邻的三个点来近似导数。
在二维情况下，紧致差分格式可以使用九点、五点或者七点的近似来计算二阶导数。
这些格式都可以通过解线性方程组的方式进行求解。
在应用紧致差分格式时，我们需要注意几个问题。
首先，边界条件的选择对于解的精度和稳定性至关重要。
通常，我们可以使用一阶导数的数值近似来设定边界条件。
其次，选择合适的离散点数和步长对于保证数值解的准确性也非常重要。
较小的步长会提高解的精度，但同时也会增加计算的复杂度。
总而言之，紧致差分格式是一种可靠且高精度的数值解法。
通过合理选择离散点和适当的近似方式，我们可以使用紧致差分格式对微分方程进行数值求解。
这种方法不仅可以应用于科学计算、工程仿真
等领域，还可以用于前沿科学研究中的模拟和模型验证。
因此，了解紧致差分格式的原理和应用，对于提高数值计算的准确性和效率具有重要的指导意义。