2.2抛物型方程的差法探究.pptVIP

定义： 网格节点上的值： 半网格节点上的值： 网格节点上的函数值 简记为u (j, n) 。 在有限差分离散化时应该注意以下几点： 根据问题求解的需要，在x，t方向上离散网格时 和 可以是等分，也可以是不等分，既可以按一定规律来离散，也可以对网格进行局部加密。 对于双曲型和抛物型等发展方程在有限差分离散化时，网格的 和 不能随意选取， 需要满足一定的条件，如稳定性的CFL条件等。 为了保证边界上的计算精度，在网格边界外可设置若干虚拟网络，以保证差分格式在边界处的计算精度和内点精度保持一致。 有限差分离散网格一般选取四边形网格，对于复杂物体外形问题，也可以选择三角形网格或其他形状网格。近年来，发展了一种无结构网格和无网格的有限差分算法，它们的计算网格就更为复杂。 对于复杂外形飞行器流场的计算，一般需要通过坐标变换，可以把物理平面上的复杂的、非正交的网格转换成在计算平面上的简单、而正交的网格，这就是网格生成技术。特别要指出的是，网格生成技术在网格设计和编程中往往占有很大的工作量，网格生成技术好坏直接影响到数值计算结果的精度，网格生成技术已成为计算流体力学中的一个重要分支。 (2) 直接差分逼近法 其中 为 之间的常数。 　把这些表达式带入微分方程后，在一阶近似条件下略去小量 ，得到相应的差分方程。例如，采用上述方法可得到： 控制体积法不是从连续的微分方程出发，而是从物理量守恒规律出发建立离散的差分方程。考虑一维平面流动问题，设流体以速度 沿 轴正方向流动。流体中某一物理量，例如某物质浓度 ，应满足守恒定律。 在空间位置 附近画出一个控制体积单元 ，如果控制体积内无源项，且不考虑扩散作用，在 内所含物质总量M应满足守恒定律，即：在控制体积 中物质总增量=由流进 内物质M的净通量。把物质浓度 近似看作是控制体积 内点j上的平均值，在时间 内控制体积 内总增量为： 从左界面流入控制体积 内物质的平均通量速率为： 在 时间内通过左边界面流入控制体积 内总通量为： 同理，在 时间内通过右边界面流出控制体积 总通量是： 流进控制体积 物质M净通量为流入通量与流出通量差： 由守恒定律得到： 在界面上函数值可以取相邻两个节点值的算术平均值，即： 把它代入守恒方程，整理后得到 若用 标记，则可以得到差分方程： 在上述推导中如果我们在界面上函数值C采用不同的近似方式，则可得到不同的差分方程。一般而言，可取： 其中 。当 时为中心差分方程，当 时为前差方程，当 时为后差方程。用控制体积法构造差分方程总是守恒型差分方程。 考虑对流方程，在矩形网格中控制体积 内，对时间和空间都取前差，即对流方程从 到 ，从 到 进行积分： 整理后得到： 这就是对流方程FTFS差分方程。 如果采用不同的积分方式，就形成不同的差分方程。但是，如果积分区域不是矩形，则采用积分方法构造差分方程就不会这么简单，但基本原理是一样的。 或 通过积分运算： 则得到： 采用积分方法构造差分方程基本思想是把微分方程在一定的控制体积内进行积分得到相应的差分方程。 (4) 积分方法 把上述方程用数值积分近似表达式来表示，即得到： 计 算 力 学 基 础 第二章 有限差分方法 2.2 抛物型方程的差分解法 2.2 抛物型方程的差分解法 一维热传导方程为： 或 对这样一个问题的求解，分为以下三个步骤来离散。 （1） 　　　　　在x-t平面上，取 和 分别为函数 的自变量x和t的改变量，由 （?j=0,1,…N, h= , n=0,1,…M, ）两组平行线构成的矩形网格覆盖x-t平面。h为空间步长， 为时间步长。 t x n j （j-1,n） （j ,n） （j+1,n） （j,n+1） （j,n-1） 2.2.1定解区域的离散 节点(j