2.2抛物型方程的差法分解.ppt

计 算 力 学 基 础;2.2 抛物型方程的差分解法;　　　　　在x-t平面上，取 和 分别为函数 的自变量x和t的改变量，由 （?j=0,1,…N, h= , n=0,1,…M, ）两组平行线构成的矩形网格覆盖x-t平面。h为空间步长， 为时间步长。; 定义： 网格节点上的值： 半网格节点上的值： 网格节点上的函数值 简记为u (j, n) 。 ; 在有限差分离散化时应该注意以下几点： 根据问题求解的需要，在x，t方向上离散网格时 和 可以是等分，也可以是不等分，既可以按一定规律来离散，也可以对网格进行局部加密。 对于双曲型和抛物型等发展方程在有限差分离散化时，网格的 和 不能随意选取， 需要满足一定的条件，如稳定性的CFL条件等。 为了保证边界上的计算精度，在网格边界外可设置若干虚拟网络，以保证差分格式在边界处的计算精度和内点精度保持一致。 有限差分离散网格一般选取四边形网格，对于复杂物体外形问题，也可以选择三角形网格或其他形状网格。近年来，发展了一种无结构网格和无网格的有限差分算法，它们的计算网格就更为复杂。;对于复杂外形飞行器流场的计算，一般需要通过坐标变换，可以把物理平面上的复杂的、非正交的网格转换成在计算平面上的简单、而正交的网格，这就是网格生成技术。特别要指出的是，网格生成技术在网格设计和编程中往往占有很大的工作量，网格生成技术好坏直接影响到数值计算结果的精度，网格生成技术已成为计算流体力学中的一个重要分支。;节点(j, n+1)的函数值在(j, n)点作泰勒展开：; 由（2）得： 其中 表示 一次和一次以 上的小量项．;（7）;下面引入几个概念： （1）向前差分 （forward difference） （2）向后差分 (backward difference) ;（3）一阶中心差分（central difference）;1）格式 I显示格式 (FTCS格式) 　　由（4）、（7）代入（1），有 　　式中： 　　称为截断误差（Truncation error），它不仅反映了差分算子对微分算子的逼近，也反映了差分解和方程解的误差。截断误差的阶数：就是截断误差中最低阶导数项中 或h的幂次数。;用 表示u(j ,n）的近似值；用差商近似代替式（1）中的微商后，可得相应 的差分方程 (9) 记号： 表示微分方程的解 在结点 处的准确值； 表示差分方程的解，它是 的近似值； 表示微分方程左端项在