有限差分：1 第一章 有限差分法 基础.docVIP

第1章 有限差分法 有限差分法是数值求解微分问题的一种重要工具，很早就有人在这方面作了一些基础性的工作。到了1910年，L. F. Richardson在一篇论文中论述了Laplace方程、重调和方程等的迭代解法，为偏微分方程的数值分析奠定了基础。但是在电子计算机问世前，研究重点在于确定有限差分解的存在性和收敛性。这些工作成了后来实际应用有限差分法的指南。20世纪40年代后半期出现了电子计算机，有限差分法得到迅速的发展，在很多领域（如传热分析、流动分析、扩散分析等）取得了显著的成就，对国民经济及人类生活产生了重要影响，积极地推动了社会的进步。 有限差分法在材料成形领域的应用较为普遍，与有限元法一起成为材料成形计算机模拟技术的主要的两种数值分析方法。目前材料加工中的传热分析（如铸造成形过程的传热凝固、塑性成形中的传热、焊接成形中的热量传递等），流动分析（如铸件充型过程，焊接熔池的产生、移动，激光熔覆中的动量传递等），都可以以有限差分方式进行模拟分析。特别是在流动场分析方面，与有限元相比，有限差分法有独特的优势，因此目前进行流体力学数值分析，绝大多数都是基于有限差分法。另外，一向被认为是有限差分法的弱项——应力分析，目前也取得了长足进步。一些基于差分法的材料加工领域的应力分析软件纷纷推出，从而使得流动、传热、应力统一于差分方式下。 可以预见，随着计算机技术的飞速发展，有限差分法将得到更为广泛的应用，可以为材料成形过程提供全面的、有效的指导。 本章主要讲述有限差分的一些基本知识,包括差分原理及逼近误差，差分方程，截断误差和相容性，收敛性与稳定性以及Lax等价定理等，这些仅仅是有限差分的入门知识，为后续章节的学习奠定基础。 1差分原理及逼近误差 1.1差分原理 设有x的解析函数y=f(x)，从微分学知道函数y对x的导数为 （4.1-1） dy、dx分别是函数及自变量的微分，是函数对自变量的导数，又称微商。相应地，上式中的、分别称为函数及自变量的差分，为函数对自变量的差商。 在导数的定义中是以任意方式趋近于零的，因而是可正可负的。在差分方法中，总是取某一小的正数。这样一来，与微分对应的差分可以有3种形式： 向前差分 （4.1-2） 向后差分 （4.1-3） 中心差分 （4.1-4） 上面谈的是一阶导数，对应的称为一阶差分。对一阶差分再作一阶差分，所得到的称为二阶差分，记为。以向前差分为例，有 （4.1-5） 依此类推，任何阶差分都可由其低一阶的差分再作一阶差分得到。例如n阶前差分为 （4.1-6） n阶的向后差分、中心差分的型式类似。 函数的差分与自变量的差分之比，即为函数对自变量的差商。如一阶向前差商为 ； （4.1-7） 一阶向后差商为 ； （4.1-8） 一阶中心差商为 （4.1-9） 或 （4.1-10） 二阶差商多取中心式，即 。 （4.1-11） 当然，在某些情况下也可取向前或向后的二阶差商。 以上是一元函数的差分与差商。多元函数f(x,y,…)的差分与差商也可以类推。如一阶向前差商为 （4.1-12） （4.1-13） 1.2逼近误差 由导数（微商）和差商的定义知道，当自变量的差分（增量）趋近于零时，就可由差商得到导数。因此在数值计算中常用差商近似代替导数。差商与导数之间的误差表明差商逼近导数的程度，称为逼近误差。由函数的Taylor展开，可以得到逼近误差相对于自变量差分（增量）的量级，称为用差商代替导数的精度，简称为差商的精度。 现将函数在x的邻域作Taylor展开： （4.1-14） （4.1-15） 这里符号O（）表示与（）中的量有相同量级的量。上式表明一阶向前差商的逼近误差与自变量的增量同量级。我们把中的指数n作为精度的阶数。这里n=1，故一阶向前差商具有一阶精度。由于是个小量，因此阶数越大精度越高。 又 （4.1-16） 一阶向后差商也具有一阶精度。 将与的Taylor展开式相减可得 （4.1-17） 可见一阶中心差商具有二阶精度。 将与的Taylor展开式相加可得 （4.1-18） 这说明二阶中心差商的精度也为二阶。 在掌握了用Taylor展开分析差商精度的方法后，再回过来谈一谈函数差分和差商的定义。由于差分和差商是微分和导数的近似表达式，所以不必局限于前面的定义，而可予以扩充。设有函数f(x)，自变量x的增量为，若取 （4.1-19） 对应的函数值为，则f(x)在xi处的n阶差分可表达为 （4.1-