数值计算法.docVIP

第十五章 工程流体力学中的数值计算方法基础 流体力学数值计算是利用计算机和数值计算方法求解流体力学具体问题的近似解，并对各种可能出现的条件进行数值模拟。随着计算机日益普及和计算方法不断发展，它已成为解决流体力学问题的重要手段。本章介绍工程流体力学数值计算方法的一些基础知识，包括代数方程的牛顿迭代法，偏微分方程的有限差分法和有限元法，以及实验数据拟合方法。 代数方程的牛顿迭代法 牛顿迭代法用于求解非线性方程 f(x) =0 的根。 基本思想 设xn是方程 f(x) =0 的1个近似根，把非线性函数f(x)在xn处作1阶泰勒级数展开，即 f(x)︽f(xn) + f’(xn) (x—xn) 则有如下近似方程 f(xn) + f’(xn) (x—xn) = 0 （15—1） 设f’(xn)=\0, 则其解为 Xn+1= xn—f(xn)/ f’(xn) (n=0,1,2, …) （15—2） 称式（15—2）为牛顿法的迭代公式。 牛顿迭代法具有明显的几何意义。方程f(x)=0的根就是曲线 y=f(x)与x轴交点的横坐标x＊，如图15—1， 设xn是x＊的第n次近似值，，过（xn ，f(xn)） 作y=f(x)切线，其切线为与x轴交点的横坐标为 Xn+1= xn—f(xn)/ f’(xn) 即用切线与x轴交点的横坐标近似代替曲线与x轴交点的横坐标。 牛顿迭代法比较简单，使用计算器就能完成计算。如果函数f(x)比较复杂。迭代次数比较多，可以用计算机完成。 求解步骤 写出f(x)的表达式，并求出f’(x)； （2）选取初始值x0，利用牛顿迭代公式（15—2）计算，迭代次数视精度要求而定。设ｅ为给定精度，当 ︱Xn+1- xn︱＜ ｅ时， 迭代结束，xn就是满足精度的近似解。 算例 例15-1，某圆形有压涵管如图15-2所示，上游水深H1.4d(管径)，此为涵管形成有压流的条件，涵管长度l=20m，上、下游水头差z=1m，通过流量Q=2m3/s，沿程阻力系数，进口、出口的局部阻力系数分别为，试确定涵管管径d。 解： 涵管为淹没出流，略去上、下游过流断面的流速水头可得 代入各项数据得： 整理后得 用牛顿迭代法求d。上式可改写为的函数形式，其一阶导数为；其n+1次迭代与第n次迭代的关系为： 式中dn+1为第n+1次迭代的管径；dn为第n 次迭代的管径。设第一次迭代以代入，则 d4=0.9182 m, d5=0.9181 m 迭代的次数视所要求的精度而定，工程实践中往往选与迭代值接近的或稍大一些的标准管径，本题可选d=1m或d=900mm的管径。现取d=1.0m的管径。 复核：当d=1m,z=1m时的流量满足所需通过的流量。 ＄15-2有限差分法 有限差分法是最早建立的一种数值解法，也是目前采用的最主要的一种数值计算法。有限差分法由于它具有许多优点，例如列式简便，求解容易，计算容量少，常被人们所采用。 1. 有限差分的基本思想 基本思想是以差商代替微商，对微分方程问题离散化，用包含有限多个未知量的差分方程近似代替微分方程，将微分方程化为代数方程，最后求得微分方程问题离散点的数值解。 有限差商的格式 一维问题： 解析函数可以在点邻域展开成泰勒级数， 如图15-3所示，设有，，三个差分节点，其坐标为，，。令函数在这三个节点的值为=，=，=，设节点间距为，则有泰勒展开式 = （15-4a） = (15-4b) 记号表示截断误差是的三次方的量级。由上两式得 （15-5） （15-6） 式（15-5）表示一阶导数的向后差分式，而式（15-6）则是向前差分。它们都具有一阶精度。式（15-4a,b）相减则得 （15-7） 上式是一阶导数的中心差分式，具有二阶精度。 如果将式（15-4a,b）相加，则有 （15-8） 上式是二阶导数的差分式，具有二阶精度。 二维问题 首先将求解域划分成许多矩形网格如图15-4所示。在x,y方向上的网络节点间距分别取=和=,这里，称为步长。将横坐标，纵坐标的空间点简写为，如空间