廖敦明《有限差分法基础》1 第一章 有限差分法 基础.pdfVIP

第 1 章 有限差分法 有限差分法是数值求解微分问题的一种重要工具，很早就有人在这方面作了 一些基础性的工作。到了 1910 年，L. F. Richardson 在一篇论文中论述了 Laplace 方程、重调和方程等的迭代解法，为偏微分方程的数值分析奠定了基础。但是在 电子计算机问世前，研究重点在于确定有限差分解的存在性和收敛性。这些工作 成了后来实际应用有限差分法的指南。20 世纪 40 年代后半期出现了电子计算机， 有限差分法得到迅速的发展，在很多领域（如传热分析、流动分析、扩散分析等） 取得了显著的成就，对国民经济及人类生活产生了重要影响，积极地推动了社会 的进步。 有限差分法在材料成形领域的应用较为普遍，与有限元法一起成为材料成形 计算机模拟技术的主要的两种数值分析方法。目前材料加工中的传热分析（如铸 造成形过程的传热凝固、塑性成形中的传热、焊接成形中的热量传递等），流动 分析（如铸件充型过程，焊接熔池的产生、移动，激光熔覆中的动量传递等）， 都可以以有限差分方式进行模拟分析。特别是在流动场分析方面，与有限元相比， 有限差分法有独特的优势，因此目前进行流体力学数值分析，绝大多数都是基于 有限差分法。另外，一向被认为是有限差分法的弱项——应力分析，目前也取得 了长足进步。一些基于差分法的材料加工领域的应力分析软件纷纷推出，从而使 得流动、传热、应力统一于差分方式下。 可以预见，随着计算机技术的飞速发展，有限差分法将得到更为广泛的应用， 可以为材料成形过程提供全面的、有效的指导。 本章主要讲述有限差分的一些基本知识,包括差分原理及逼近误差，差分方 程，截断误差和相容性，收敛性与稳定性以及 Lax 等价定理等，这些仅仅是有限 差分的入门知识，为后续章节的学习奠定基础。 1 差分原理及逼近误差 1.1 差分原理 设有 x 的解析函数 y=f(x)，从微分学知道函数 y 对 x 的导数为 x xfxxf x y dx dy xx Δ ?Δ+ = Δ Δ = →Δ→Δ )()(limlim 00 （4.1-1） dy、dx 分别是函数及自变量的微分， dx dy 是函数对自变量的导数，又称微商。相 应地，上式中的 yΔ 、 xΔ 分别称为函数及自变量的差分， x y Δ Δ 为函数对自变量的 差商。 在导数的定义中 xΔ 是以任意方式趋近于零的，因而 xΔ 是可正可负的。在差 分方法中， xΔ 总是取某一小的正数。这样一来，与微分对应的差分可以有 3 种 形式： 向前差分 )()( xfxxfy ?Δ+=Δ （4.1-2） 向后差分 )()( xxfxfy Δ??=Δ （4.1-3） 中心差分 ) 2 1() 2 1( xxfxxfy Δ??Δ+=Δ （4.1-4） 上面谈的是一阶导数，对应的称为一阶差分。对一阶差分再作一阶差分，所 得到的称为二阶差分，记为 y2Δ 。以向前差分为例，有 )()(2)2( )]()([)]()2([ )()( )]()([ )(2 xfxxfxxf xfxxfxxfxxf xfxxf xfxxf yy +Δ+?Δ+= ?Δ+?Δ+?Δ+= Δ?Δ+Δ= ?Δ+Δ= ΔΔ=Δ （4.1-5） 依此类推，任何阶差分都可由其低一阶的差分再作一阶差分得到。例如 n 阶前差 分为 )]}()(([{ )]}([{ )]([ )( 2 1 xfxxf y y yy n nn ?Δ+ΔΔΔ= ΔΔΔΔ= ΔΔΔ= ΔΔ=Δ ? ? L L LL （4.1-6） n 阶的向后差分、中心差分的型式类似。 函数的差分与自变量的差分之比，即为函数对自变量的差商。如一阶向前差 商为 x xfxxf x y Δ ?Δ+ = Δ Δ )()( ； （4.1-7） 一阶向后差商为 x xxfxf x y Δ Δ?? = Δ Δ )()( ； （4.1-8） 一阶中心差商为 xxxfxxf x y Δ Δ??Δ+ = Δ Δ )2 1() 2 1( （4.1-9） 或 x xxfxxf x y Δ Δ??Δ+ = Δ Δ 2 )()( （4.1-10） 二阶差商多取中心式，即 22 2 )( )()(2)( x xxfxfxxf x y Δ Δ?+?Δ+ = Δ Δ 。 （4.1-11）