北航空气动力学课件.pptVIP

从“野渡无人舟自横”到“香蕉球”技术 浅谈“香蕉球”的力学原理 2010年版本 3．4 二维对称物体绕流的数值解 　 把直匀流和分布的偶极子（或总强度为零的分布点源和点汇）叠加起来，所得到的组合流动为对称封闭物体绕流。设直匀流沿 x 轴正向流来，其速度为V∞，在 x 轴上 x =a和 x =b 范围内连续分布一系列的偶极子，单位长度内偶极子的强度设为m（偶极子密度）。 　 如果偶极子密度的分布形式已知，则离原点距离为 ξ 的小区间内由偶极子产生的流函数为： 总流函数为： 。物体的外形可以用零流线来表示。改变不同的偶极子密度分布，可以获得不同形状的封闭物体，由流函数和速度以及速度与压强的关系确定流场中各点及物体表面的速度分布和压强分布。 奇点叠加数值解法 2010年版本 北京航空航天大学《空气动力学》国家精品课 3．4 二维对称物体绕流的数值解 　 2010年版本 北京航空航天大学《空气动力学》国家精品课 3．4 二维对称物体绕流的数值解 　 对于实际问题，往往是给定物体的外形来确定其流动的特性。在这种情况下，偶极子密度分布函数的确定需要由流函数求解。对偶极子密度来说，流函数是一个积分方程，求它的解是比较困难的。但是随着计算机技术的发展，可以用数值方法比较迅速地获得这种方程的有一定准确度的数值解。 下面简单地叙述用数值方法求解已知物体形状确定绕物体流动特性的过程。　 2010年版本 北京航空航天大学《空气动力学》国家精品课 3．4 二维对称物体绕流的数值解 1、数值解法步骤 　 首先，我们把偶极子分布区域分成等宽度的 n 段，设每段的宽度为△ξ，段数 n 可根据计算机容量及结果的准确度要求而确定。流场中某一定点P处的流函数为： 　 式中 为第 j 段中点离原点的距离； 为第 j 段内偶极子密度的平均值； 　 表示第 j 段内偶极子的强度。 2010年版本 北京航空航天大学《空气动力学》国家精品课 3．4 二维对称物体绕流的数值解 用物面边界条件来确定待求的偶极子密度 对于给定物体外形上的n个已知点（xi，yi），就可以得到一个对未知函数的 n 元一次联立代数方程组： 其中 为影响系数，表示 处的单位偶极子 密度对物体表面某点Pi (xi, yi) 处的流函数贡献。 2010年版本 北京航空航天大学《空气动力学》国家精品课 3．4 二维对称物体绕流的数值解 展开上式，即 利用解一次方程组的各种计算方法，求解上面方程组，确定偶极子密度。 2010年版本 北京航空航天大学《空气动力学》国家精品课 3．4 二维对称物体绕流的数值解 一旦所给定物体外形的偶极子密度分布已经解得，则可以确定流场内任意点处的流函数。此后即可由流函数与速度的关系式及伯努利方程，确定流场内各点处的速度及压强值。 在上述过程中，我们实际上是把第j段中分布的偶极子用集中在该段中点处的等强度偶极子来代替了。显然，如果分段数量较多，这种近似表示才有一定的准确性。理论上，当段数 n 趋于无限大时，偶极子密度分布的数值结果趋近于精确解。在实际应用时，由于计算机容量和计算机机时的限制，以及多元一次联立方程组解的不稳定性。分段的数目不宜太多。 也可以由位函数出发，用位函数对应的物面条件来解决实际流动问题。这两种方法是等价的。在实际应用中，用位函数叠加法比用流函数法更广泛。 本章基本要求 掌握平面不可压位流中位函数与流函数的性质与关系。 掌握平面不可压位流的基本方程即拉普拉斯方程的特点、叠加原理和边界条件； 掌握四种基本而重要的位流流动即：直匀流，点源（点汇）、偶极子和点涡的表达； 重点掌握直匀流与偶极子和点涡的叠加； 掌握儒可夫斯基升力定律； 了解二维对称物体绕流数值解法步骤 本章讨论怎样求解不可压理想流体无旋运动的规律 在理想不可压条件下欧拉方程和连续方程包括四个方程和四个未知函数（u,v,w,p），理论上是可解的 由于飞行器的外形都比较复杂，要在满足如此复杂的边界条件下求该偏微分方程组的解析解是非常困难的，原因在于方程包含非线性项，而且方程中速度与压强相互耦合，需要一并求出 人们发现在无旋条件下问题可以得到大大简化，尤其是可以将速度和压