第二章理论基础(全章文本文件).pdfVIP

第二章 有限差分法理论基础 有限差分方法是计算流体力学中应用最多的离散化数值方法； 作为计算技术它是历史最悠久，理论上相对成熟的数值方法。 本章主要参考书： 1.Anderson Dale A. Computational Fluid Mechanics and Heat Transfer （前4 章） 2.张涵信、沈孟育，计算流体力学－差分方法的原理和应用，国防工业出版社 3.忻孝康等 计算流体动力学， 国防科技大学出版社 4. 4.John D. Anderson, Computational Fluid Dynamics,The Basics with Applications , 计算流体力学入门，清华大学出版社（影印版）（第4 章） $ 2.1 有限差分离散化方法 一个定解的流体动力学问题的数学描述； ⎧∂u ⎪∂t L(u ) (x,t )=∈Ω×[0,T ] ⎪ ⎨ u(x.o) φ(x) ⎪ ⎪ l(u) |δΩ 0 ⎩ 其中L(u) 为对空间坐标的微分算子。 问题的解法: ⑴ 理论（解析）解 ⑵ 差分数值解； 方程的离散， 求解域（时+空）的离散， 代数方程的求解。 数值解的讨论。 （一）求解域的离散化——差分分割。 分割尺寸（空间网格步长，时间步长） （网格）结（节）点，（网格）单元 （边界外）虚网格点——网格延拓。 （二）微商（偏导数）的差商近似 1）。差商近似；一阶，二阶导数的偏心差分、中心差分格式 一阶偏微商的定义； ∂u u(x +Δx, y ) −u(x , y ) Lim 0 0 0 0 ∂x Δ→x 0 Δx (x0 ,y0 ) 1 u x x y u x y ( 0 +Δ , 0 ) − ( 0 , 0 ) ∂u 若取消取极限过程，用 ，代替 |(x0 ,y 0 ) 就是一种差商近 Δx ∂x 似。称为差分格式， 2 ）。差分格式的导出方法 a). Taylor’s 公式展开； 2 2 n n ∂u ∂ u Δx ∂ u Δx u(xo +Δx , y 0 ) u(x0 , y 0 ) + 0 Δx + 2 0 +⋅⋅⋅⋅⋅⋅+ n ζ (x 0 ≤ξ ≤x0 +Δx) ∂x ∂x 2! ∂x n! ( , ) ( , ) ∂u u x 0 +Δx y 0 −u x 0 y 0 ⇒ +T .E . T.E=Truncation Error ∂x (x0