- 1、本文档共35页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
计算方法8.1-8.2
练习: 用类似方法分析一阶前差的空间近似性如何? * * * * 2. 二阶空间差商 * * 分析 想办法将上式分子分母凑成这种形式 则 * * 通过本节分析,得到下列结论: (1)空间差分是空间微商的一种近似; (2)近似好坏与几种因素有关 -空间网格大小(似乎越小越好?) -函数本身性质(函数越平滑,波长越长好) -使用空间差商,对于那些短波,更明确的说, 波长小于2△x的波是无法描述的。 * * 第八章 差分方程概论 * * * * 8.1 有限差分离散化 第八章 差分方程概论 8.2 离散近似 8.3 初值问题差分格式的有效性 * 大气运动方程组是一组非线性偏微分方程组,用来描述系统的状态及其运动和变化规律,这类方程通常无法求出解析解,常常借助计算机采用数值方法求它们的近似解,在这一过程中,首先要将微分方程离散化。数值计算方法中,离散化的方法一般有三种: * * (1)用差商代替微商,使偏微分方程变成差分方程,而差分方程是一个代数方程。这样,微分方程组变成代数方程组,然后求解线性代数方程组,得到微分方程的数值解。这种方法称为“差分法”。应用最为广泛。 (2)谱方法:利用一些基函数(如球谐函数),把解展开成其有限项的线性组合,可以在一定意义下近似满足方程和定解条件,再利用基函数的特性(如正交性等),使方程化为展开系数和其对时间微商的常微分方程组,然后用差商代微商求数值解。 (3)把偏微分方程问题变成相应的泛函数极小问题求解的方法,如利兹法,有限元法。 * * 8.1 有限差分离散化 * * 8.1.1 差商的概念 设f(x)在x0的某领域内具有直到n+1阶的连续导数,由Taylor公式,有 (截断误差) (1) (2) 注:h为一较小数 * * (3) (4) (5) (6) f(x)在x0处向前差商 f(x)在x0处向后差商 f(x)在x0处中心差商 * * (7) f(x)在x0处二阶差商 * * 同理,若设二元函数u(x,y)在(x,y)的某领域内有连续的n+1阶偏导数,则由二阶Taylor展开式,有 (8) (10) (9) 前差 后差 中心差 * * 关于x的二阶中心差商 关于y的二阶中心差商 以上差商与网格系统密切相关,h和k分别为x和y方向的步长 * * 8.1.2 解域离散化和差分网格的建立 采用数值方法求解控制方程时,都是想办法将控制方程在空间区域上进行离散,然后求解得到的离散方程组。要想在空间域上离散控制方程,必须使用网格将计算域进行有规则的分割,分割的交点称为网格点,网格点间的距离称为格距。如果给定边界条件及初始时刻的气象要素值,就可计算出在这些格点上以后时刻的气象要素值。由于大气运动是近似水平的,以下将主要说明在二维平面上的网络分割。 一、解域离散化 正三角形网络 正方网络 正六角形 * * 常见的有规则的网格分割主要有三种:正方形网格或矩形网格,正三角形网格和正六角形网格。 * * 在相同的区域范围内和相等的格距条件下,取正方形、正三角形和正六角形网格,所能得到的格点数目之比依次近似为5:8:4。 正三角形网格:密度最高,具有较高精度,但是由于程序设计比较复杂,在实际应用中较少使用。 正六角形网格:密度最低,精度也较差,程序设计也复杂,实际应用中也很少采用。 正方形(矩形)网格:精度较高,程序设计简单,是应用最为广泛的网格。 以x-y平面的二维矩形网格为例,具体说明网格的建立,最简单常用的是“等步长”差分分割。取适当的间隔△x, △y的平行直线群进行分割。 二、网格的建立 * * (i,j) ﹡ ﹡ ﹡ ﹡ ﹡ (i,j+1) (i+1,j) (i,j-1) (i-1,j) 用 两族直线将x-y平面的第一象限剖分成矩形网格(如图) 其中的(xi, yj)(简写为(i, j))称为网格点, △x, △y为x, y方向的网格距(相当于前面的h, k) 函数u(x,y)在(i,j)处的值可表示为 * * 以上是空间网格的建立。 同理,对时间也可建立分割,取适当的分割间隔△t,称为时间步长,其分割点记为 如此建立了时间的网格系统。 注: △x ,△y, △t不能随意取,要根据实际问题和差分方法内在性质的要求适当选取。上述是最简单的等格距正交差分网格的分割,但在实际问题中,有时要求非等距分割,即 (8) 前差 * * 则 (8)′ 同理 (9) 后差 (9)′ (10) 中心差 * * (10)′ (11) 二阶中心差 (11)′ 同理可写出关于t求导的差商 前差 后差 * * 中心差 二阶中心差 (同学们试着自己写出
文档评论(0)