有限差分法基础课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《有限差分法基础课件.ppt》由用户(三亚风情)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限 差分法 基础 课件
- 资源描述:
-
1、第二章 有限差分法主讲人:胡才博主讲人:胡才博中国科学院大学地球科学学院中国科学院大学地球科学学院中国科学院计算地球动力学重点实验室中国科学院计算地球动力学重点实验室第二章 有限差分法 2.1 有限差分法基础 2.2 网格剖分 2.3 差分格式 2.4 差分方程 2.5 应用实例1. 地球内部介质,不仅存在纵向非均匀结构(一维地球模型),也存在横向非均匀结构(不同块体、断层系统);2. 几何模型也呈现出相当的复杂性;3. 另外,边界条件和初始条件对于不同问题具有特殊性。解析方法的局限性解析方法的局限性Hu, C., Y. Cai, and Z. Wang (2012), Effects of
2、large historical earthquakes, viscous relaxation, and tectonic loading on the 2008 Wenchuan earthquake, Journal of Geophysical Research, 117, B06410, doi:10.1029/2011JB009046. (SCI, IF: 3.303)汶川大地震的动力学成因对于存在复杂介质和几何、特殊边界条件和初始条件的实际地质问题,一般不存在解析解,需要近似的数值求解方法。有限差分方法是地球物理方法中最常见的一种。有限差分方法(Finite Difference
3、 Method, FDM)是计算机数值模拟最早采用的方法,至今仍被广泛使用。有限差分方法的基本特点该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。2.1 有限差分法基础 有限差分法以变量离散取值后对应的函数值来近似微分方程中独立变量的连续取值。 我们放弃了微分方程中独立变量可以取连续值的特征,而关注独立变量离散取值后对应的函数值。 有限差分法的具体操作分为两个部分: (1)用差分代替微分方程中的微分,将连续变化的变量离散化,从而得到差分方程组的数学形式; (2)求解差分方程组。 有限差分方法的基本原理该方法将求解域划分为差分网格,用
4、有限个网格节点代替连续的求解域。有限差分方法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的函数值为未知数的代数方程组。2.1 有限差分法基础有限差分法的主要内容1. 建立地球物理问题的离散有限差分模型(1)如何根据问题的特点将定解区域做网格划分;(2)如何在所有网格节点上用有限差分格式对导数求近似,对函数、初始条件和边界条件求近似;(3)如何把原方程离散化为代数方程组,即有限差分方程组。2.从理论上研究有限差分模型的形态,以保证计算过程的可行性和计算结果的正确性(1)解的相容性;(2)解的稳定性;(3)解的收敛性。3. 如何数值求
5、解差分方程组 网格剖分就是研究区域和边界的离散化 1.矩形分割 2.三角形分割 3.极网格分割2.2 网格剖分对地球物理问题的连续求解区域通过网格划分离散为空间上得一系列网格点,接下来需要利用一定的差分格式对偏微分方程组中的导数用差商进行近似,从而将偏微分方程组离散化为差分方程组。对于函数f(x),通常意义下的导数(微商)定义为:()()lim()()lim()()lim2xxxfxdxfxdxfxfxdxdxfxdxfxdxdx当dx0时,以上三种形式都是微商的正确定义。如果dx是有限的,如何给出微商的近似定义?2.3 差分格式用Taylor级数级数展开可以给出微商的近似形式。对于连续函数f
6、(x),它在相邻点上的值f(x+x)和f(x- x)可以用Taylor 级数展开为dx变为x223344234223344234()()()()()()()()()2 !3!4 !()()()()()()()()()2 !3!4 !dxdxdxdfxxfxxfxfxfxfxdxdxdxdxdxdxdxdfxxfxxfxfxfxfxdxdxdxdx 如果x很小,f(x)可微,则以上级数收敛。次数越高,收敛级数的项的绝对值越小。由(1)得到,223344234()()()()()()()()()2 !3!4 !()()()()dxdxdxdfxxfxxfxfxfxfxdxdxdxdxdfxxfxf
7、xOxdxx (1)(2)(3)(4)()()()()dfxxfxfxOxdxx 式中的O(x)项表示忽略掉的所有项中的最大项的量级是x,也就是说,忽略掉这些项带来的误差中的最大项和x成正比。()()()dfxxfxfxdxx 由(4)给出导数的一阶精度(first order accurate)近似为:(4)(5)(5)式称为向前差分格式向前差分格式(forward-difference formula)由(2)式得到223344234()()()()(-)()-()()-()2 !3!4 !()()()+()dxdxdxdfxfxxxfxfxfxfxdxdxdxdxdfxfxxfxOxdx
8、x 由(7)式得到导数的另一个一阶精度近似:()(-)()dfxfxxfxdxx(6)(7)(8)(8)式称为向后差分形式向后差分形式(backward-difference formula)。(1)式减去(2)式,得到:2()(-)()=+()2dfxxfxxfxOxdxx (9)式中的O(x2)项表示忽略掉这些项带来的误差中的最大项和x2成正比。由(9)式得到导数的二阶精度(second order accurate)近似为:()(-)()2dfxxfxxfxdxx (10)式称为中心差分形式中心差分形式(central-difference formula)。(9)(10)2233442
9、34223344234()()()()()()()()()2 !3!4 !()()()()()()()()()2 !3!4 !dxdxdxdfxxfxxfxfxfxfxdxdxdxdxdxdxdxdfxxfxxfxfxfxfxdxdxdxdx (1)(2)223344234223344234()()()()()()()()()2 !3!4 !()()()()()()()()()2 !3!4 !dxdxdxdfxxfxxfxfxfxfxdxdxdxdxdxdxdxdfxxfxxfxfxfxfxdxdxdxdx (1)(2)2222()2()(-)()=+()()dfxxfxfxxfxOxdxx
10、(1)式和(2)式相加,得到:222()2()(-)()()dfxxfxfxxfxdxx (12)式称为二阶导数的二阶精度中心差分形式为二阶导数的二阶精度中心差分形式。忽略x的四次方及更高阶项(11)(12) f(xi+h)-f(xi): 节点xi的一阶向前差分 f(xi)-f(xi-h): 节点xi的一阶向后差分 f(xi+h)-f(xi-h): 节点xi的一阶中心差分 前后是相对x轴正方向而言()()()dfxxfxfxdxx ()(-)()dfxfxxfxdxx222()2()(-)()()dfxxfxfxxfxdxx ()(-)()2dfxxfxxfxdxx 总结:总结:1、向前差分形
11、式:2、向后差分形式:3、中心差分形式:单侧,一阶精度单侧,一阶精度对称,二阶精度对于二阶导数二阶精度对一阶导数 定解问题的有限差分解法1离散 x = ih, y= jh, i= 0, 1, 2,. n, h: 步长(正方形的边长)2根据泰勒级数建立差商格式:对于一维情况:在x处的一阶导数可以用3. 建立和求解差商方程组差分格式差分格式的另一种推导的另一种推导为了寻求更精确的差分格式,我们引入两个待定常数,由泰勒展开,构造如下关系式(i,j)(i+1,j)(i-1,j)01234(i,j)(i+1,j-1)(i-1,j-1)(i,j+1)(i+1,j+1)(i-1,j+1)i-1ii+1j-1
展开阅读全文