流动传热与燃烧的数值计算2差分课程.ppt

第一章：离散化方法 离散化就是将原微分方程在整个求解域内所考虑的物理量求解问题，转化为对求解域内有限个位置（网格结点）上物理量的求解问题。?用各网格结点上定义的离散状态量来代替在求解域内连续分布的状态量。 对于给定的微分方程，对应有不同的离散化方程，它们之间的差别是由于所取的分布以及推导离散化方程的方法不同造成的。 §1.1有限差分法 差分法的基本思想是用差商代替微商。 用差分法离散微分方程的具体方法有很多种，对同一微分方程可以建立不同的差分方程。 §1.1.1差分方程的建立 以一维问题为例，设x3-x2=x2-x1=Δx 由Taylor级数展开可得： 则 可以表示为： 则 可以表示为： §1.1.2网格划分 离散过程包括两个方面：①微分方程的离散化； ②求解域的离散化。 设微分方程有两个自变量x，t 。 微分方程的解应是x和t的函数u(x,t)，而差分法（数值方法）求的是若干个离散点（x=j Δx，t=nΔt ）上u(x,t)的近似值 因此，首先要在求解域上指定一系列离散点。 一般取时间步长Δt =const.，空间步长Δx =const. 对于二维和三维问题，不要求Δx = Δy= Δz。 tn=nΔt 和xj=jΔx 的直线为网格线，其交点( xj , tn)称为网格结点。 网格设计好之后，就可以把连续求解域上的微分方程离散化为有限个网格结点上的差分方程。 例：用差分法求解抛物型方程 设： tn=nΔt ，xj=jΔx ，n=0,1,2,….,N，j= 0,1,2,….,J 点( xj , tn)可以用(j ,n)表示 u( xj , tn)用 表示 令微分方程中的一阶导数用向前差分，二阶导数用中心差分，则与原微分方程对应的差分方程为： 显式格式：解可以由明显的公式计算出来，不需要求解方程组。FTCS 隐式格式：解不能简单地求出，需要求解代数方程组。全隐式格式 §1.1.3 差分方程的相容性 用差分方程代替微分方程之后产生了一个问题，即所采用的差分方程是否逼近原来的微分方程？ 相容性是差分方程的一个基本特性，它既反映了差分方程是否收敛到相应的微分方程，也反映了截断误差。 截断误差：用ujn代替 u(xj, tn) 的差分方程与在点 (xj, tn)处的微分方程之差。 E=O(Δtq+Δxp) 差分方程的精度：如果E=O(Δtq+Δxp)，则称此差分方程对t是q阶精度，对x是p阶精度。 相容性：如果当网格步长Δx→0, Δt→0时，差分方程的截断误差也E→0 ，此时在每一网格结点上的差分方程与原微分方程是等同的，则该差分方程与相应的微分方程是相容的。 §1.1.4 差分方程的收敛性 差分方程能否应用，要看差分方程的近似解能否任意地逼近微分方程的准确解。 首先要考虑差分方程理论上的精确解能否任意地逼近微分方程的准确解，其次还要考虑求解差分方程过程中产生的误差。 离散误差：差分方程的精确解ujn与微分方程的准确解 u(xj, tn) 之差。 ejn =ujn - u(xj, tn) 收敛性：如果当网格步长Δx→0, Δt→0时，任何网格结点上的离散误差ejn →0 ，即 ujn →u(xj, tn)，则该差分方程是收敛的。 离散误差的大小与方程的截断误差有关。 在网格步长相同的情况下，截断误差的阶数提高会使离散误差ejn 减小。 对同一离散格式，网格加密也会使离散误差ejn 减小。 收敛性是讨论当 Δx→0, Δt→0 时，差分方程的精确解是否收敛于原微分方程的准确解。 相容性是讨论当Δx→0, Δt→0 时，差分方程是否逼近原微分方程，相容并不能保证收敛。 举例：对于FTCS格式，当S=1, 1/2, 1/6时，考察结点xj=0.4，tn=8处的数值解，并与微分方程在该处的准确解 u(0.4, 8)=45.03963比较。 其中α=0.01，l=1.0 §1.1.5 差分方程的稳定性 求解差分方程的过程中不可避免地会引入误差（舍入误差、初始误差等），稳定性讨论差分方程解的误差在计算过程中的发展问题。 稳定性：如果求解差分方程过程中引入的数值误差在以后计算步逐步消失或保持有界，则称此差分方程是稳定的；如果引入的数值误差在以后计算步被逐渐放大，以致物理问题的解被完全破坏，则称此差分方程是不稳定的。 1.直观法（离散摄动法） 在某计算步引入误差（离散摄动）后，直接考察在以后各计算步中差分方程对误差发展所起的作用。 以FTCS格式为例，讨论其稳定性 ： ： ： ： ： ： ： : : 2.Von Neumann方法（傅里叶级数法、分离变量法