应用空气动力学-2

合集下载

飞行空气动力学--固定翼飞机结构详解---2

飞行空气动力学--固定翼飞机结构详解---2

⼒。

外壳摩擦⼒是最难降低的寄⽣阻⼒类型。

没有完全光滑的表⾯。

甚⾄是机械加⼯的表⾯,通过放⼤来检测的话,仍然可以看到粗糙的不平坦的外观。

这种粗糙的表⾯会使表⾯的空⽓流线型弯曲,对平滑⽓流产⽣阻⼒。

通过使⽤光滑的磨平的表⾯,和去掉突出的铆钉头,粗糙和其他的不规则物来最⼩化外壳摩擦⼒。

设计飞机时必须要增加另⼀个对寄⽣阻⼒的考虑。

这个阻⼒复合了形阻⼒效应和外壳摩擦,称为所谓的⼲涉阻⼒。

如果两个物体靠近放置,产⽣的合成紊乱会⽐单个测试时⼤50%到200%。

形阻⼒,外壳摩擦⼒和⼲涉阻⼒这三个阻⼒都要被计算以确定⼀个飞机的寄⽣阻⼒。

寄⽣阻⼒中⼀个物体的外形是⼀个很⼤的因素。

然⽽,说道寄⽣阻⼒时指⽰空速也是⼀个同样重要的因素。

⼀个物体的外形阻⼒保持在⼀个相对⽓流固定的位置,⼤约以速度的平⽅成正⽐增加;这样,空速增加为原来的两倍,那么阻⼒就会变成原来的四倍,空速增加为三倍的话阻⼒也就增加为九倍。

但是,这个关系只在相当的低⾳速时维持很好。

在某些更⾼速度,外形阻⼒的增加会随速度⽽变的突然很快。

第⼆个基本的阻⼒类型是诱导阻⼒。

以机械运动⽅式⼯作的系统没有⼀个可以达到100%的效率,这是⼀个确定的物理事实。

这就意味着⽆论什么特性的系统,总是以系统中消耗某些额外的功来获得需要的功。

系统越⾼效,损失就越⼩。

在平飞过程中,机翼的空⽓动⼒学特性产⽣要求的升⼒,但是这只能通过某种代价才能获得。

这种代价的名字就叫诱导阻⼒。

诱导阻⼒是内在的,在机翼产⽣升⼒的任何时刻,⽽事实上,这种阻⼒是升⼒的产物中不可分离的。

继⽽,只要有升⼒就会有这种⼒。

机翼通过利⽤三种⽓流的能量产⽣升⼒。

⽆论什么时候机翼产⽣升⼒,机翼下表⾯的压⼒总是⼤于机翼上表⾯的压⼒。

结果,机翼下⽅的⾼压区空⽓有向机翼上⽅的低压去流动的趋势。

在机翼的翼尖附近,这些压⼒有区域相等的趋势,产⽣⼀个从下表⾯到机翼上表⾯的向外的侧⾯⽓流。

这个侧向⽓流给予翼尖的空⽓和机翼后⾯的尾流⼀个旋转速度。

第二章 空气动力学

第二章 空气动力学
➢ 流体微团在宏观上无限小,在微观上无限大。
2.1.3 流场、 定常流和非定常流
➢ 流场:流体流动所占据的空间称为流场。 ➢ 流场的选取可根据研究的需要进行确定。可大可小。
非定常流与定常流
➢ 非定常流与非定常流场:
在流场中的任何一点处,如果流体微团流过时的流动参 数——速度、压力、温度、密度等随时间变化,这种流动 就称为非定常流,这种流场被称为非定常流场。
的地方, 却流得比较快。 夏天乘凉时, 我们总喜欢坐在两座房屋之间的过道中, 因
为那里常有“ 穿堂风”。 在山区你可以看到山谷中的风经常比平原开阔的地方来得
大。
连续方程
质量守恒定律
➢ 质量守恒定律是自然界基本的定律之一, 它说明物质既不 会消失, 也不会凭空增加。
➢ 应用在流体的流动上: 在定常流动中,当流体低速、稳 定、连续不断地流动时, 流进任何一个截面的流体质量
➢ 只要相对气流速度相同 , 产生的空气动力也就相等。
(非定常流动转换为定常流动)
风洞实验
➢ 将飞机的飞行转换为空气的流动 ,使空气动力问题的研 究大大简化。
➢ 风洞实验就是根据这个原理建立起来的。
风洞应用
相对气流方向的判定
➢ 相对气流的方向与飞机运动的方向相反 。
平飞时:
相对气流方向 飞行速度方向
➢ 对于不可压缩的、理想的流体( 没有粘性) 表示为:
p1v2
2
p0
常数
静压
动压
总压
➢ 静压:单位体积流体具有的压力能。在静止的空气中, 静压等于大气压力。 ➢ 动压:单位体积流体具有的功能。 ➢ 总压:静压和动压之和。
p1v2
2
p0
常数
➢ 上式即为:不可压缩的、理想的流体( 没有粘性) 的伯努利 方程。

第二章_空气动力学(民航大学)

第二章_空气动力学(民航大学)
早期低速飞机可以通过调整外撑杆的长度来调整 机翼的安装角:加大安装角叫内洗,减小安装角 叫外洗。
2.3 机体几何外形和参数
上反角和下反角:机翼底 面与垂直机体立轴平面之 间的夹角,ψ 。
纵向上反角:机翼安装角 与水平尾翼安装角之差。
机身的几何形状和参数
机身长度Lsh、最大当量 直径Dsh及其所在轴向相 对位置和长细比 λ sh=Lsh/Dsh。
附面层转变的原因
气流流过机体表面的距离越长,附面层越厚。 机体表面过于粗糙、凹凸不平。
层流附面层和紊流附面层
紊流附面层VS层流附面层
紊流附面层比层流附面层厚,底部的横向速度 梯度也比层流的大。紊流附面层对气流的阻滞 作用比层流附面层大。
层流附面层和紊流附面层的速度型
气流在机体表面的流动状态
机翼的空气动力
α 小迎角下作用在机翼上的空气动力
伯努利定理的应用
阻力
阻力是与飞机运动轨迹平行,与飞行速度方向相反的 力。阻力阻碍飞机的飞行,但没有阻力飞机又无法稳 定飞行。
阻力的分类
对于低速飞机,根据阻力的形成原因,可将阻力 分为:
•摩擦阻力(Skin Friction Drag) •压差阻力(Form Drag) •干扰阻力(Interference Drag)
影响因素:
空气的粘性 附面层内气流的流动状态(紊流大于层流)。 机体与气流的接触面积越大,机体表面越粗糙,
摩擦阻力越大。
层流附面层和紊流附面层的速度型
摩擦阻力在飞机总阻力构成中占的比例较大
超音速战斗机 大型运输机 小型公务机 水下物体 船舶
摩擦阻力占总阻力的比例 25-30% 40% 50% 70% 90%
废阻力

空气动力学与热工基础 2—3三角翼的空气动力特性

空气动力学与热工基础 2—3三角翼的空气动力特性
在气流尚未分离的引角下, 升力系数随迎角的增大而直线增长,升力系数与迎角表 现为线性变化关系。例如歼7飞机就是这样,参见图3— 1—15所示。

而细长三角翼具有不寻常的升力特性,其不同点为: 升力系数曲线的斜率比大展弦比机翼小得多;其随迎角 的变化呈现非线性,升力系数的增长比迎角更快一些, 如图3—2—34所示。其所以如此,是因为升力由两部 分组成。一部分是翼面的附着流(整个下表面和部分上 表面)所产生的升力,叫做“位流升力”,其变化与迎 角成线性关系。另一部分是上表面脱体涡所产生的升力, 叫“涡升力”,其变化与迎角成非线性关系。
• 在很小的迎角下,上式可写成 • 2 C y K p K N
Kp 式中第一项是位流升力,第二项是涡升力; 与 Kn 均为常值,其大小取决于展弦比。图3—2—36表明了 按上式计算的结果与实验结果的比较。 • 当迎角增大到一定程度,脱体涡在机翼上表面后缘 发生破碎,变得不规则,这会使流谱发生变化。迎角进 步增大,破碎点向前移动,能量进一步耗散,涡升力减 小。再后,出现失速,升力相应下降。临界迎角可高 达 35 ~ 40 。 • 三角翼虽然有这样大的临界迎角,但起飞、着陆, 还很难得到充分利用。因为起飞、着陆,增大迎角或迎 角过大,势必影响飞行员的视界,还会造成机身尾部擦 地。例如歼7飞机起飞的着陆迎角,不超过 9 ~ 10,远 远小于临界迎角。为此,有的超音速飞机将机头向下折 转,改善视界。同时,加高起落架,防止机尾擦地。


脱体涡具有增大上表面吸力,使升力增大的作 用。因为脱体涡从前缘连续产生,形成稳定的低压区, 上表面正处于脱体涡低压之下,所以吸力很大。迎角 大,低压区吸力也大,所以升力增大更多。图3—2— 35是一个展弦比为1的三角翼,在20°迎角下的各个 横断面上压强分布图。它说明了上表面在脱体涡所复 盖的区域,吸力很大。 • 据理论分析结果:细长三角翼的升力系数( C y ) 与迎角( )之间的关系,如下式所示: C y K p sin a cos 2 K N cos a sin 2

《空气动力学》课件

《空气动力学》课件

未来挑战与机遇
环境保护需求
新能源利用
随着环境保护意识的提高,对空气污 染和气候变化的研究需求增加,这为 空气动力学带来了新的挑战和机遇。
新能源的利用涉及到流动、传热和燃 烧等多个方面,需要空气动力学与其 他学科合作,共同解决相关问题。
航空航天发展
航空航天领域的发展对空气动力学提 出了更高的要求,需要不断改进和完 善现有技术,以满足更高性能和安全 性的需求。
04
翼型与机翼空气动力学
翼型空气动力学
翼型概述
翼型分类
翼型是机翼的基本截面形状,具有特定的 弯度和厚度。
根据弯度和厚度的不同,翼型可分为超临 界、亚音速和超音速翼型等。
翼型设计
翼型与升力
翼型设计需考虑气动性能、结构强度和稳 定性等多个因素。
翼型通过产生升力使飞机得以升空。
机翼空气动力学
01
机翼结构
课程目标
掌握空气动力学的基本概 念和原理。
提高分析和解决实际问题 的能力。
了解空气动力学在各领域 的应用和发展趋势。
培养学生对空气动力学的 兴趣和热爱。
02
空气动力学基础
流体特性
01
02
03
04
连续性
流体被视为连续介质,由无数 微小粒子组成,彼此之间存在
相对运动。
可压缩性
流体的密度会随着压力和温度 的变化而变化。
《空气动力学》PPT课件
目 录
• 引言 • 空气动力学基础 • 流体动力学 • 翼型与机翼空气动力学 • 空气动力学应用 • 未来发展与挑战
01
引言
主题介绍
空气动力学:一门研 究空气运动规律和空 气与物体相互作用的 科学。
课件内容涵盖了基础 理论、应用实例和实 验演示等方面。

空气动力学第二章第二部分讲解

空气动力学第二章第二部分讲解
(Cp ) (Cpn )n cos2 (Cy ) (Cyn )n cos2
Cy (Cy )n cos
(Cx ) (Cxn )n cos3
四、后掠翼低速气动特性-后掠翼流动特点
翼根前段:流管粗,扩张,V ,Cp ; 翼根后段:流管变细,V ,Cp ,Cpmin 后移;
亚音速前(后)缘
或 m 1
超音速前(后)缘
或 m 1
一、主要概念回顾(续)
V
V
d d
(z) 1 l 2 d
4V l 2 z
改变了实际迎角,有效迎角为
e(z) a(z) (z)
二、升力线理论 — 升力
剖面假设:各剖面展向速度分量 以及流动参数沿展向的变化比其 他方向小得多,剖面流动为二维。
库塔-儒可夫斯基定理
R(z) Ve (z)(z)
机翼单位展长翼段升力可表示为:
Y
(z)

1 2
V2C

y
(
z)
b( z)
1

1 2
V2b( z )

C
y
a
(z)


(
z)
Y(z) V(z)
(z)

1 2
Vb( z )C y
a
(z)


(z)
d d
(z) 1 l 2 d
翼根效应:翼根剖面最小压强点后移, 升力贡献下降; 翼尖前段:流管变细,V ,Cp ,Cpmin 前移。 翼尖后段:流管变粗,V ,Cp 。
翼尖效应:翼尖剖面最小压强点前移升 力增加。
翼尖先失速
四、后掠翼低速气动特性-后掠翼气动特性

2024年空气动力学总结(2篇)

2024年空气动力学总结(2篇)

2024年空气动力学总结____年是一个科技进步迅速的年代,空气动力学领域也在不断发展和创新。

在这篇总结中,我将为您介绍____年空气动力学的最新进展,包括技术发展、应用领域和研究成果。

1. 技术发展1.1 高效翼型设计:____年,在空气动力学领域,一项重要的技术进展是高效翼型的设计和优化。

通过使用先进的计算流体力学模拟和优化算法,工程师们能够设计出更加优化和高效的翼型,以减小飞机的阻力并提高飞机的升力系数。

1.2 翼尖涡减阻技术:翼尖涡是飞机在飞行过程中产生的一种涡旋,会增加飞机的阻力。

____年,工程师们开发出了一项翼尖涡减阻技术,通过在翼尖上安装一种新型的尖状装置,能够有效减小翼尖涡的产生,从而降低飞机的阻力和燃油消耗。

1.3 多孔翼面技术:多孔翼面是一种新型的翼面结构,____年,科研人员取得了一系列突破性进展。

多孔翼面的优点是能够有效减小风阻,提高飞机的升力系数,并且在一定程度上能够吸收和减小噪音。

这项技术在未来有望得到广泛应用。

2. 应用领域2.1 商用航空:商用航空是空气动力学的一个重要应用领域。

____年,随着技术的不断进步,商用航空公司将能够开发出更加高效和环保的飞机,减少对化石燃料的依赖,并降低排放。

2.2 无人机:无人机是近年来迅速发展起来的一种飞行器,广泛应用于农业、测绘、物流等领域。

随着空气动力学技术的不断提高,____年的无人机将具备更长的续航能力、更高的飞行速度和更稳定的飞行性能。

2.3 超音速运输:超音速运输是一个具有巨大潜力的领域,通过超音速运输,人们能够更快速、更高效地到达目的地。

____年,随着技术的进步,科研人员将为超音速飞行器研发出更加高效和稳定的空气动力学设计,推动超音速运输的发展。

3. 研究成果3.1 超轻复合材料:在____年,科研人员取得了重要的突破,开发出了一种超轻复合材料。

这种材料具有高强度、高刚度和低密度的特点,能够极大地降低飞机的重量,提高飞机的燃油效率和航程。

空气动力学课后答案

空气动力学课后答案

空气动力学课后答案【篇一:空气动力学复习题】txt>第一章低速气流特性1.何谓连续介质?为什么要作这样的假设?2.何谓流场?举例说明定常流动与非定常流动有什么区别。

流场——流体所占居的空间。

定常流动——流体状态参数不随时间变化;非定常流动——流体状态参数随时间变化;3.何谓流管、流谱、流线谱?低速气流中,二维流谱有些什么特点?流线谱——由许多流线及涡流组成的反映流体流动全貌的图形。

流线——某一瞬间,凡处于该曲线上的流体微团的速度方向都与该曲线相应点的切线相重合。

流管——通过流场中任一闭合曲线上各点作流线,由这些流线所围成的管子。

二维流谱——1.在低速气流中,流谱形状由两个因素决定:物体剖面形状,物体在气流中的位置关系。

2.流线的间距小,流管细,气流受阻的地方流管变粗。

3.涡流大小决定于剖面形状和物体在气流中的关系位置。

4.写出不可压缩流体和可压缩流体一维定常流动的连续方程,这两个方程有什么不同?有什么联系?方程可变为:va=c(常数)气流速度与流管切面积成反比例。

方程可变为:适用于理想流体和粘性流体5.说明气体伯努利方程的物理意义和使用条件。

方程表达式: p?1?v2??gh?常量 21?v2?p0?常量2高度变化不大时,可略去重力影响,上式变为:p?即:静压+动压=全压 (p0相当于v=0时的静压)方程物理意义:空气在低速一维定常流动中,同一流管的各个截面上,静压与动压之和(全压)都相等。

由此可知,在同一流管中,流速快的地方,压力(p)小;流速慢的地方,压力(p)大。

方程应用条件1.气流是连续的、稳定的气流(一维定常流);2.在流动中空气与外界没有能量交换;3.空气在流动中与接触物体没有摩擦或摩擦很小,可以忽略不计(理想流体);4.空气密度随流速的变化可忽略不计(不可压流)。

图1-7 一翼剖面流谱p1+?v12=p2+?v22=p3+?v32v1a1=v2a2=v3a3v2=200米/秒p2=-3273675帕斯卡v3=83米/秒p3=445075帕斯卡7.何谓空气的粘性?空气为什么具有粘性?空气粘性——空气内部发生相对运动时,相邻两个运动速度不同的空气层相互牵扯的特性。

2空气动力学基础-第2章 流体动力学

2空气动力学基础-第2章 流体动力学

§2.1.3 流线、流管、流面与流量
人们希望用一些曲线将流场上的流动情况表现出来。在某一 瞬间看流场的话,从某点出发,顺着这一点的速度指向画一 个微小的距离到达邻点,再按邻点在同一瞬间的速度指向再 画一个微小距离,一直画下去便得一条曲线。这条某瞬时的 空间曲线,其切线都和该点的微团速度指向相一致。这样的 空间曲线称为流线,这样的线可以画无数条。
§2.1.1 拉格朗日方法与Euler方法
流体质点的其它物理量也都是 a,b,c,t 的函数。例如流体 质点(a,b,c)的温度可表为T(a,b,c,t) 2、Euler方法(Euler方法,空间点法,流场法) •Euler方法的着眼点不是流体质点而是空间点。考察不同流 体质点通过空间固定点的流动行为,通过记录不同空间点流 体质点经过的运动情况,从而获得整个流场的运动规律。 •在固定空间点看到的是不同流体质点的运动变化,无法像 拉格朗日方法那样直接记录同一质点的时间历程。
p p( x, y, z, t ),
( x, y , z , t , )
T T ( x, y , z , t )
§2.1.1 拉格朗日方法与Euler方法
如果场只是空间坐标的函数而与时间无关则称为定常场, 否则为非定常场,例如,定常速度场的表达为:
u u ( x, y, z ), v v( x, y, z ), w w( x, y, z )
注意上式并非全导数的表达(在微积分中当复合函数 只是一个自变量 t 的函数时才有全导数),因为在Euler观 点下 x、y、z 等与时间 t 无关,不能写出 dx/dt 的表达。
§ 2.1.2 Euler法的加速度表达式
算子: u v w t x y z Material derivative:往往用D/Dt这样一个符号来表示。 这个导数称为随流体运动的导数,或称随体导数、实质导 数或物质导数。

空气动力学习题答案

空气动力学习题答案

1.1解:)(k s m 84.259mk R 22328315∙===-RT p ρ=36m kg 63.5063032.5984105RT P =⨯⨯==ρ气瓶中氧气的重量为354.938.915.0506.63G =⨯⨯==vg ρ 1.2解:建立坐标系根据两圆盘之间的液体速度分布量呈线性分布 则离圆盘中心r ,距底面为h 处的速度为0u kn u +=当n=0时 u=0推出0u 0= 当n=h 时 u=wr 推出hwr k = 则摩擦应力τ为hwr u dn du u==τ 上圆盘半径为r 处的微元对中心的转矩为θθτdrd hwr u r rdrd h wr u r dA d 3=⋅=⋅=T则⎰⎰==T 2D 0332032D u drd hr uωπθωπ1.4 解:在高为10000米处T=288.15-0.0065⨯10000=288.15-65=223.15压强为⎪⎭⎫ ⎝⎛=T a T Pa P 5.2588MKN43.26Ta T pa p 2588.5=⎪⎭⎫ ⎝⎛=密度为2588.5T a T a ⎪⎭⎫⎝⎛=ρρmkg4127.0Ta T a 2588.5=⎪⎭⎫⎝⎛=∴ρρ1-7解:2M KG 24.464RTPRT p ==∴=ρρ空气的质量为kg 98.662v m ==ρ 2-3解:将y 2+2xy=常数两边微分2ydy+2xdx+2ydx=0整理得ydx+(x+y )dy=0 (1) 将曲线的微分方程yx V dyV dy =代入上式得 yVx+(x+y )V y =0 由22y 2xy 2x V ++=得V x 2+V y 2=x 2+2xy+y 2 ((2)由(1)(2)得()y v y x v y x =+±=, 习题二2-2解:流线的微分方程为yx v dyv dx =将v x 和v y 的表达式代入得ydy xdx yx 2dyxy 2dx 22==, 将上式积分得y 2-x 2=c ,将(1,7)点代入得c=7 因此过点(1,7)的流线方程为y 2-x 2=482-5解:直角坐标系与柱坐标系的转换关系如图所示 速度之间的转换关系为{θθθθθθcos v sin v v sin v cos v v r y r x +=-=由θθθθθθcos r1y v sin yrsin r 1xvcos x rrsin y rcos x =∂∂=∂∂⎪⎩⎪⎨⎧-=∂∂=∂∂⇒⎭⎬⎫==()()⎪⎭⎫⎝⎛--∂∂+-∂∂=∂∂∂∂+∂∂⋅∂∂=∂∂θθθθθθθθθsin r 1sin V cos V cos sin V cos V r x v v x r r v x v r r x x xθθθθθθθθθθθθθs i n c o s V s i n V s i n V c o s V r 1c o s s i n r V c o s r V r r r ⎪⎭⎫ ⎝⎛-∂∂--∂∂-⎪⎭⎫ ⎝⎛∂∂-∂∂=θθθθθθθθθθθθθθcos sin V r1sin V r 1sin V r 1cos sin V r 1cos sin r V cos r V 22r r 2r +∂∂++∂∂-∂∂-∂∂=()()θθθθθθθθθcos r1cos V sin V sin cos V sin V r y v v V y r V V V V r r y x y xy +∂∂++∂∂=∂∂⋅∂∂+∂∂⋅∂∂=∂∂θθθθθθθθθθθθθcos r1sin V cos V cos V sin V sin cos r V sin r V r r r ⎪⎭⎫ ⎝⎛-∂∂++∂∂+⎪⎭⎫ ⎝⎛∂∂+∂∂=θθθθθθθθθθθθθcos sin V r1cos V r 1cos V r 1cos sin v V r 1cos sin r V sin r V 22r r 2r -∂∂++∂∂+∂∂+∂∂=zV V V r 1r V z V y V x V div z r r z y x ∂∂+⎪⎭⎫ ⎝⎛∂∂++∂∂=∂∂+∂∂+∂∂=∴θυθ 2-6解:(1)siny x 3x V 2x -=∂∂ s i n y x 3y V 2y =∂∂ 0yV x V y x =∂∂+∂∂ ∴此流动满足质量守恒定律(2)siny x 3x V 2x =∂∂ s i n y x 3y V 2y =∂∂ 0siny x 6yV x V 2y x ≠=∂∂+∂∂ ∴此流动不满足质量守恒定律(3)V x =2rsin rxy2=θ V y =-2rsin 2ry 22-=θ33ry 2x Vx =∂∂332yr 2y y x 4y V +-=∂∂0ryx 4y V x V 32y x ≠-=∂∂+∂∂∴此流动不满足质量守恒方程(4)对方程x 2+y 2=常数取微分,得xdy dy dx -= 由流线方程yx v dy v dx =(1) 由)(得2r k v v r k v 422y 2x =+= 由(1)(2)得方程3x r ky v ±= 3y rkx v = 25x r kxy3x V =∂∂∴ 25y rkxy 3yV ±∂∂ 0yV x V yx =∂∂+∂∂∴此流动满足质量守恒方程2—7解:0xVz V 0r yz 23r yz 23z V y V z x 2727y z =∂∂-∂∂=⋅+⋅-=∂∂-∂∂同样 0y V x V x y =∂∂-∂∂∴该流场无旋()()()2322222223222z y x z y x z y x d 21zy x z d z y d y x d x dz v dy v dx v d ++++⋅=++++=++=Φ c zy x 1222+++-=Φ∴2—8解:(1)a x V x x =∂∂=θ a yV y y =∂∂=θ a z Vz z -=∂∂=θ 021v ;021v ;021v z y x =⎪⎪⎭⎫⎝⎛∂∂+∂∂==⎪⎭⎫ ⎝⎛∂∂+∂∂==⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=y V x V x V z V z V x V x x z x y z (2)0y V x V 210x V z V 210z V y V 21x y z z x y y z x =⎪⎪⎭⎫⎝⎛∂∂-∂∂==⎪⎭⎫ ⎝⎛∂∂-∂∂==⎪⎪⎭⎫ ⎝⎛∂∂-∂∂=ωωω;; 位该流线无旋,存在速度∴ (3)azdz 2aydy axdx dz v dy v dx v d z y x -+=++=ϕc az ay 21ax 21222+-+=∴ϕ 2—9解:曲线x 2y=-4,()04y x y x f 2=+=, 切向单位向量22422422y2x 2y2x yx 4x xy 2i yx 4x x j f f fx i f f fy t +-+=+-+=t t v v v t ⋅∇=⋅=∇=ϕϕ切向速度分量 把x=2,y=-1代入得()()j x 2x i y x 2x j yi x v 2+-+--=∂∂+∂∂=∇=ϕϕϕ j 21i 21j y x 4x 2xy i y x 4x x t 2242242+=⎪⎪⎭⎫ ⎝⎛+-+= 23t v v t -=⋅= j 23i 23j 21i 2123t v v t t --=⎪⎭⎫⎝⎛+-== 2—14解:v=180hkm =50sm根据伯努利方程22V 21V 21p ρρρ+=+∞∞ pa p =∞ 驻点处v=0,表示为1531.25pa 501.22521V 21pa p 22=⨯⨯==-∞ρ相对流速为60sm 处得表示为75.63760225.12125.1531V 21V 21pa p 222-=⨯⨯-=-=-∞ρρ 习题三3—1解:根据叠加原理,流动的流函数为()xyarctg 2Q y V y x πϕ+=∞, 速度分量是22y 22x yx y2Q x V y x x 2Q V y V +⋅=∂∂-=+⋅+=∂∂=∞πϕπϕ; 驻点A 的位置由V AX =0 V Ay =0求得 0y V 2Qx A A =-=∞;π 过驻点的流线方程为2x y arctg 2y x y arctg 2y y Q V Q V A A A =+=+∞πθπ θθππθππsin 2r x y arctg 2y -⋅=⎪⎭⎫ ⎝⎛-=∞∞V V Q 或即 在半无限体上,垂直方向的速度为θπθθππ-sin v r sin 2y x y 2v 222y ∞==+=Q Q 线面求极值()0-sin v -cos sin v 2d dv 22y=+=∞∞θπθθπθθθ 当0sin =θ 0v v miny y ==2-tg -=θπθmaxyy v v =用迭代法求解2-tg -=θπθ得 取最小值时,y 1v 2183.1139760315.1==θ 取最大值时,y 2v 7817.2463071538.4 ==θ由θπθθππ-sin v r sin 2y x y 2v 222y ∞==+=Q Q θπθθθππ-cos sin v r cos 2v y x x 2v v 22x +=+=++=∞∞∞Q Q 可计算出当∞∞===v 6891574.0v v 724611.0v x y 1,时,θθ 6891514.0v v 724611.0v x y 2=-==∞,时,θθ合速度∞=+=v v v 2y 2x V3—3解:设点源强度为Q ,根据叠加原理,流动的函数为xa3-y a r c t g2a x y a r c t g 2a x y a r c t g 2πθπθπθϕ+++-=两个速度分量为()()()⎥⎥⎦⎤⎢⎢⎣⎡+++++++--=222222a 3-y x xy a x a x y a x a x 2x πθ()()()⎥⎥⎦⎤⎢⎢⎣⎡++++++-=222222y a 3-y x a3-y y a x y y a x y 2v πθ 对于驻点,0v v y x ==,解得a 33y 0x ==A A , 3—4解:设点源的强度为Q ,点涡的强度为T ,根据叠加原理得合成流动的位函数为Q ππθϕ2l n r 2Γ+=πθϕπθϕθ2r 1r 12r 1r r Γ=∂∂==∂∂=V V ; 速度与极半径的夹角为Qarctg arctgr Γ==V V θθ 3—5根据叠加原理得合成流动的流函数为⎪⎪⎭⎫ ⎝⎛+--+=∞y a y yaarctg a y y aarctgV ϕ 两个速度分量为()()()()⎥⎦⎤⎢⎣⎡++---+++=∂∂=∞1y v 2222x y a x a x a y a x a x a V ϕ ()()⎥⎦⎤⎢⎣⎡+--++=∂∂-=∞2222y y v y a x yy a x y a V ϕ 由驻点()0a 30,得驻点位置为±==y x v v 零流线方程为0ay yaarctg a y y xaarctgy =--++∞∞V V 对上式进行改变,得⎪⎭⎫ ⎝⎛-=-+a y tan ay2a y x 222当0x =时,数值求解得a 03065.1y ±= 3—9解:根据叠加原理,得合成流动的流函数为a y y a r c t g 2a y y a r c t g 2y v -++-=∞ππϕQ Q速度分量为()()2222x y a x ax 2y a x a x 2y v v +-+++++-=∞ππQ Q ()()2222y ya x ax 2y a x a x 2v +-+++++-=ππQ Q 由0v v y x ==得驻点位置为⎪⎪⎭⎫ ⎝⎛+±∞0v a a 2,πQ 过驻点的流线方程为ay y arctg 2a y y arctg 2y v =-++--∞ππQ Q 上面的流线方程可改写为ay y arctga y y arctg y v 2--+=∞Q π 222a y x ay2a y y arctg a y y arctg tan y v 2tan -+=⎪⎪⎭⎫ ⎝⎛--+=⎪⎪⎭⎫ ⎝⎛∴∞Qπ 容易看出y=0满足上面方程当0y ≠时,包含驻点的流线方程可写为⎪⎭⎫ ⎝⎛-=-+∞Q y v 2tan ay2a y x 222π当12v a ===∞πQ 时,包含驻点的流线方程为tanyy21y x 22--=-+ 3—10解:偶极子位于原点,正指向和负x 轴夹角为α,其流函数为 22yx x s i n y c o s 2+--=ααπϕM 当45=α时22yx xy 222+--=πϕM 3—11解:圆柱表面上的速度为a2sin v 2v πθΓ--=∞ 222222a4a 2s i n v 4v ππθΓ+Γ=∞ 222222v a 4av 2sin 4sin 4v v ∞∞∞Γ+Γ+=⎪⎪⎭⎫ ⎝⎛ππθθ压强分布函数为222p v asin 41sin 41v v 1⎪⎪⎭⎫ ⎝⎛Γ+-=⎪⎪⎭⎫ ⎝⎛-=∞∞θπθC习题四4—1解:查表得标准大气的粘性系数为nkg 1078.1u 5-⨯= 65el 1023876.11078.16.030225.1u⨯=⨯⨯⨯==-∞LV R ρ 平板上下两面所受的总得摩擦阻力为N S V L R F 789.021e 664.0222=⨯⨯=∞ρ 4—2解:沿边阶层的外边界,伯努利方程成立代表逆压梯度代表顺压梯度,时;当时当0m 0m 00m 00m m v v v 21p 12201002〈〉∴〉∂∂〈〈∂∂〉-=-=∂∂-=∂∂=+--xpx p x v x v x v xx p c m m m ρρρρδδδ 4—4解:(a )将2x y 21y 23v v ⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛=δδδ带入(4—90)中的第二式得δδδδδ28039dy vv 1v v 0x x =⎪⎪⎭⎫ ⎝⎛-=⎰** 由牛顿粘性定律δτδuu 23y v u 0y xw =⎪⎪⎭⎫⎝⎛∂∂==下面求动量积分关系式,因为是平板附面层 0dx dv =∴δ积分关系式可表示为dxd v 2w **=δρτδ 将上述关系式代入积分关系式,得δρδδv dxud 14013=边界条件为x=0时,0=δ 积分上式,得平板边界层的厚度沿板长的变化规律()64.428039646.0x x x64.4ll ⨯==∴=**R R δδ(b )()74.164.483x x 83dy v v 1lx =⨯=∴=⎪⎪⎭⎫ ⎝⎛-=*∞*⎰R δδδδ(c )由(a )知()64.4x x l =R δ(d )646.0x x646.0v 21324xx 64.4u23l f l 2wf l w =∴====R C R C R δρτδδδτ)得—由(; (e )单面平板的摩擦阻力为()292.1x x 292.1s v 21b bdx v 21l f l 2f l02f=∴===⎰R C R X C C X F F δδρρ摩阻系数为假设版宽为4—6解:全部为层流时的附面层流厚度由式(4—92)得 ()01918.048.5L e ==LR Lδ 全部为湍流时的附面层流厚度由式(4—10)得()0817.037.0L 51e ==-LLR δ第五章5—3证明(1)将r (θ)表示为下列三角级数()⎪⎭⎫⎝⎛+=∑∞=∞1n 0n s i n n s i n c o s v 2r θθθθA A 将其代入(5—35)得()∑∞==+-1n f10dxdy n ncos θαA A 可得⎰⎰=-=ππθθπθπα011fn 01f 0d cosn dxdy 2d dx dy 1A A ;对于平板,0dx dy f =,故有α=0A ,()θθαθsin cos v 2r 0n 21∞=∴===A A A 当πθ→时,()0r ≠π,不满足后缘条件(2)将()⎪⎭⎫⎝⎛++=∑∞=∞1n 0nsins sin cos 1v 2r θθθθA A 将其带入(5—35)积分得()αθθθθθθθθθπππ-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-+-+-⎰⎰∑∞=∞dxdy d cos cos sin nsinn cos cos d cos 1v 1f 0021n 2A()∑∞==+-1n 1f10s i n dy n ncos θθαA A⎰-=1f 0d dx dy 1θπαA ⎰=πθθπ011fn d c o s n dx dy 2A 对于平板0dxdyf =,0n 210====∴A A A A ;α()θθαθsin cos 1v 2r +=∴∞当πθ→时,()0r =θ,满足后缘条件5—2解:设在41弦线处布涡的强度为Γ,则该涡在43弦线处产生的诱导速度为c2c 2v yi ππΓ=Γ=若取43弦点为控制点,在改点满足边界条件⎪⎭⎫ ⎝⎛-=Γ∴⎪⎭⎫ ⎝⎛-=Γ∞∞απαπdx dy cv dx dy v c f f 因此开力为⎪⎭⎫⎝⎛-=Γ-=∞∞dx dy cv v f 2αρπρL 开力系数为⎪⎭⎫ ⎝⎛-==∞dx dy 2c v 21f 2απρLC L 对于平板0dx dy f =ππαα22==∴L L C C ;5—4解对于薄翼型,πα2=LC 对于2412翼型,()()1x 4.0x 28.00555.0dxdy 4.0x 0x 28.081dx dy ff ≤≤-=≤≤-=;; 令()1cos 121x θ-=,则当x=0.4时,2.0arccos 1=θ ()()π≤≤-=≤≤-=x 2.0a r c c o s 0.28.00555.0dxdy 2.0arccos x 00.28.081dx dy ff ;;()()()112.0a r c c o s1101f 0d c o s 12.0c o s 811d c o s 1dx dy 1θθθπθθπαπ--=-=∴⎰⎰()()112.0a r c c o s1d c o s 12.0c o s 0555.01θθθππ--+⎰101fn d c o s n dxdy 2θθππ⎰=A()()()()⎥⎦⎤⎢⎣⎡--+--=⎰⎰12.0arccos 1112.0arccos 011cos 12.0cos 0555.02d cos 12.0cos 812θθπθθθππA ()()⎥⎦⎤⎢⎣⎡-+-=⎰⎰πθθθθθθπ2.0arccos 111112.0arccos 012d cos22.0cos 0555.0d cos22.0cos 812A ()214mp 4A A C -=π5—5解:根据余弦定理9924.0c 9849.0abcosc 2b a c 222=∴=-+=9962.0cb c o s ca ac 2b abcosc 2b a a 2ac b c a cos 2222222=-=--++=-+=B 059878.4==∠∴B折算后的迎角为010,()()1x 32170tan dx dy 32x 05tan dx dy d cos 1dxdy 120f 0f 101f00≤≤=≤≤=-=-=⎰;;;θθπαααππL C令()弧度时当9106.131arccos 32x cos 121x 11=⎪⎭⎫⎝⎛-==-=θθ ()()119106.1019106.10100d cos 1tan1701d cos 15tan 1θθπθθπαπ-+-=∴⎰⎰()()⎰⎰-=-+-=9106.10119106.101101253.0d cos 1tan170d cos 15tan θθπθθπ()8837.11253.018010220=⎥⎦⎤⎢⎣⎡+⨯=-=∴ππααπL C 5—7解:()()()x 2x 3x k 2x 1-x kx y 23f +-=-=()2x 6x 3k dx dy 2f +-= 令0dx dy f =得()正号舍去331±=x ()6x 6k dx y d 2f 2-=将331-=x 代入,得0dx y d 2f2〈 因此f y 在331-=x 处取得极大值,2f =%将331-=x 代入f y 得k=0.052 令()1cos 121x θ-=代入(1)得k 41cos 23cos 43dx dy 112f ⎪⎭⎫ ⎝⎛-+=θθ ()110f0d cos 1dxdy 1θθπαπ-=∴⎰()()0235.11105.00524.0220=-=-=∴πααπL C 07794.0d cos dx dy 2110f1==⎰θθππA 04587.0d dxdy 110f0=-=⎰θπαπA0186.0d cos2dx dy 2110f 2=⎪⎭⎫⎝⎛=⎰θθππA ()533.0210=+=πA A C L ()1798.041412-=--=L L C A A C π6—5解:根据开力线理论()()ζζδζπδd d d 41v 22yi Γ-=⎰-LL已知()2122021202112d d 21⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛-Γ-=Γ⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛-Γ=ΓL L L ζζζζδ; ()11122220yi d sin 2d cos 2cos 2d 213v 21θθζθζθζζζδζζπδL L L L L L L =-=-=-⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫⎝⎛-Γ=∴⎰-;;;令 则⎪⎭⎫ ⎝⎛-Γ-=-Γ-=⎰θθθθθθθππsin 3sin 183d cos cos cos sin 3v 01011122yi L L当LLL L 43v 283v 3240yi 0yi Γ-===Γ-===,时,时πθζπθζ6—6解(1)有叠加原理可知,a 处的下洗速度为⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎪⎭⎫ ⎝⎛+Γ-=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎪⎭⎫ ⎝⎛+Γ-⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛Γ-=a a 21a 2a 1242a 22a 22a 4v 22222222yi L L L L L L L L πππa 处的下洗角α为L V V L C L LV V L ∞∞∞∞Γ==⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎪⎭⎫ ⎝⎛+Γ=-=λρπα221a a 21v 222yi ; 因此a 2L V C L ∞=Γ代入下洗角中得⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎪⎭⎫ ⎝⎛+=a a 21222L C L πλα (2)对于椭圆翼()()00222121ααλπλπλππααπλαα-+=+=-+=∞∞L L L C C C()02222i 1aa 2211a a 22d ααλπλ-⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡++⎪⎭⎫ ⎝⎛+=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡++⎪⎭⎫ ⎝⎛=L L C L ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡++⎪⎭⎫ ⎝⎛+=∴1aa 221d dd 22i L λα当4.0a 8==,λ时26.0d dd i=α6—9解:1268.41;274.0s 21-∞∞∞=+===rad C C C V L C L L LL αααρ00013.22.1354.3;354.3=-===-ααααLLC C00385.02==πλLDi C C 6—11解:()09985.01;846.0s 2122=+===∞δπλρLDi L C C V L C71.41017N;s 212===∞Lx V C x i Di i ρ% 第七章7—1解状态方程RT ρ=p3212312123121321300v v w v v 21a 25.1019a 62.506a 62.506T T K T KP P KP P KP P ;;;;;;;;========ρρρρρ(1)由状态1等压膨胀到2的过程中,根据质量守恒方程12v 2v =所以1221ρρ=等压变化K T T T T T T 600221221122211====∴=;ρρρρ 由32→等容变化,根据质量方程23ρρ= 等容变化2323223322T T T T T P T P ==∴=; (2)介质只在21→过程中膨胀做功KJ 53.21v p w =∇= (3)()996.182m v p =+=T C T C Q δ(4)161.466KJ pdv -q du pdv du q ==∴+=δδ(5)k kj 298.0ln s r 2112v =∆∴⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛=δρρδP P C 7—3解根据质量守恒小截面与2A 截面的流量相等即()()()()25.0388.0q q q c q c2211220201010=∴==∴=λλλλλA A T A P T A P7—4解:气流从Ma=1加速到Ma1=1.5需要的外折角度为091.11='δ总的外折角度091.2615=+'=δδ 查表得Ma2=2.02456.010********=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫⎝⎛=⋅=P P P P P P P P P P 7—5解:经过正激波时绝热,总温度0T 不变根据总静温之比1r 2a 21r 1020+=*∴-+=T T M T T 1r r 2r 1r 200+=*=+=*∴*RT RT C T T ;波后的速度系数为1r r 2v v 0222+==*RT C λ根据波前波后的速度关系121=λλ 1r r 2v 1021+=∴RT λ 根据马赫数与速度系数的关系,得得波德马赫数2121211r 1r 11r 2a λλ+--+=M 总压损失系数δ为()()1r 121211r 1212a 1r a 1r 1r 1r a 1r r 2---⎥⎦⎤⎢⎣⎡+-+⎪⎭⎫⎝⎛+--+=M M M δ。

空气动力学期末复习题 (2)

空气动力学期末复习题 (2)

第一章一:绪论;1.1大气的重要物理参数 1、最早的飞行器是什么?——风筝2、绝对温度、摄氏温度和华氏温度之间的关系。

——95)32(⨯-T =T F C15.273+T =T C K6、摄氏温度、华氏温度和绝对温度的单位分别是什么?——C ο F ο K ο 二:1.1大气的重要物理参数1、海平面温度为15C ο时的大气压力为多少?——29.92inHg 、760mmHg 、1013.25hPa 。

3、下列不是影响空气粘性的因素是(A)A 、空气的流动位置B 、气流的流速C 、空气的粘性系数D 、与空气的接触面积4、假设其他条件不变,空气湿度大(B)A 、空气密度大,起飞滑跑距离长B 、空气密度小,起飞滑跑距离长C 、空气密度大,起飞滑跑距离短D 、空气密度小,起飞滑跑距离短 5、对于音速.如下说法正确的是: (C)A 、只要空气密度大,音速就大B 、只要空气压力大,音速就大C、只要空气温度高.音速就大D、只要空气密度小.音速就大6、大气相对湿度达到(100%)时的温度称为露点温度。

三:1.2 大气层的构造;1.3 国际标准大气1、大气层由内向外依次分为哪几层?——对流层、平流层、中间层、电离层和散逸层。

2、对流层的高度.在地球中纬度地区约为(D)A、8公里。

B、16公里。

C、10公里。

D、11公里3、现代民航客机一般巡航的大气层是(对流层顶层和平流层底层)。

4、云、雨、雪、霜等天气现象集中出现于(对流层)。

5、国际标准大气指定的依据是什么?——国际民航组织以北半球中纬度地区大气物理性质的平均值修正建立的。

6、国际标准大气规定海平面的大气参数是(B)A、P=1013 psi T=15℃ρ=1、225kg/m3B、P=1013 hPA、T=15℃ρ=1、225 kg/m3C、P=1013 psi T=25℃ρ=1、225 kg/m3D、P=1013 hPA、T=25℃ρ=0、6601 kg/m37. 马赫数-飞机飞行速度与当地音速之比。

空气动力学试题 (2)

空气动力学试题 (2)

编号:西北工业大学考试试题(卷)2009 -2010 学年第二学期开课学院航空学院课程流体力学基础学时与自由来流叠加形成该圆柱绕流的偶极子的强度应如何变化?(10)答:1.由Cp=1-4 sin*sin=0得 =30 150 210 3302.加倍2.注:1. 命题纸上一般不留答题位置,试题请用小四、宋体打印且不出框。

2. 命题教师和审题教师姓名应在试卷存档时填写。

共2页 第1页西北工业大学命题专用纸三、计算题(共45分)注:标准海平面大气压强:,密度:1. 一无限薄的平板(即忽略其厚度),其弦长为m,超音速来流攻角为,平板上下表面的压强和剪应力分布分别为:,,,,其中x为平板上的点据平板前缘的距离,压强和剪应力的单位为:N/m2,求平板所受到的法向力N和轴向力A,升力L和阻力D,以及绕前缘和处的力矩。

(10分)2. 翼型在海平面的飞行速度为60m/s,气流沿翼型表面某点的速度为90m/s,求该点的压强。

若翼型放置在入口连接海平面大气的直流风洞(风扇在实验段后)实验段内,测得该点有相同的压强,则该点处的速度为多少?(10分)3. 薄翼型的零升攻角为,弦长2米,以60 m/s的时速在海平面飞行,若其单位展长的机翼产生的升力是1695牛,则该翼型的升力系数是多少?其飞行攻角是多少?(10分)4. 椭圆无扭转机翼的展弦比为10,展长10米,重2400kg。

构成该机翼的薄翼型零升攻角为-10,剖面型阻系数。

若机翼在海平面以100m/s匀速平飞,求此时的升力、升力系数、诱导阻力系数、推力、飞行时的几何攻角。

()(10分)教务处印制 共2页 第2页。

空气动力学课后答案(北航)

空气动力学课后答案(北航)

钱第一章 1.1解:)(k s m 84.259m k R 22328315•===-RT p ρ=36m kg 63.5063032.5984105RT P =⨯⨯==ρ 气瓶中氧气的重量为354.938.915.0506.63G =⨯⨯==vg ρ1.2解:建立坐标系根据两圆盘之间的液体速度分布量呈线性分布 则离圆盘中心r ,距底面为h 处的速度为0u kn u +=当n=0时 u=0推出0u 0= 当n=h 时 u=wr 推出hwr k =则摩擦应力τ为hwr u dn du u ==τ上圆盘半径为r 处的微元对中心的转矩为θθτdrd hwr u r rdrd h wr u r dA d 3=⋅=⋅=T则⎰⎰==T 2D 0332032D u drd hr uωπθωπ1.4解:在高为10000米处T=288.15-0.0065⨯10000=288.15-65=223.15压强为⎪⎭⎫ ⎝⎛=T a T Pa P 5.2588MKN43.26Ta T pa p 2588.5=⎪⎭⎫ ⎝⎛=密度为2588.5T a T a ⎪⎭⎫⎝⎛=ρρmkg4127.0Ta T a 2588.5=⎪⎭⎫⎝⎛=∴ρρ1-7解:2M KG 24.464RTPRT p ==∴=ρρ空气的质量为kg 98.662v m ==ρ第二章2-2解流线的微分方程为yx v dyv dx =将v x 和v y 的表达式代入得ydy xdx yx 2dyxy 2dx 22==, 将上式积分得y 2-x 2=c ,将(1,7)点代入得c=7 因此过点(1,7)的流线方程为y 2-x 2=482-3解:将y 2+2xy=常数两边微分 2ydy+2xdx+2ydx=0整理得ydx+(x+y )dy=0 (1) 将曲线的微分方程yx V dy V dy =代入上式得 yVx+(x+y )V y =0 由22y 2xy 2x V ++=得 V x 2+V y 2=x 2+2xy+y 2 ((2)由(1)(2)得()y v y x v y x =+±=,2-5解:直角坐标系与柱坐标系的转换关系如图所示 速度之间的转换关系为{θθθθθθcos v sin v v sin v cos v v r y r x +=-=由θθθθθθcos r1y v sin yrsin r 1xvcos x rrsin y rcos x =∂∂=∂∂⎪⎩⎪⎨⎧-=∂∂=∂∂⇒⎭⎬⎫==()()⎪⎭⎫⎝⎛--∂∂+-∂∂=∂∂∂∂+∂∂⋅∂∂=∂∂θθθθθθθθθsin r 1sin V cos V cos sin V cos V r x v v x r r v x v r r x x xθθθθθθθθθθθθθsin cos V sin V sin V cos V r 1cos sin r V cos r V r r r ⎪⎭⎫ ⎝⎛-∂∂--∂∂-⎪⎭⎫ ⎝⎛∂∂-∂∂=θθθθθθθθθθθθθθcos sin V r1sin V r 1sin V r 1cos sin V r 1cos sin r V cos r V 22r r 2r +∂∂++∂∂-∂∂-∂∂=()()θθθθθθθθθcos r1cos V sin V sin cos V sin V r y v v V y r V V V V r r y x y xy +∂∂++∂∂=∂∂⋅∂∂+∂∂⋅∂∂=∂∂θθθθθθθθθθθθθcos r1sin V cos V cos V sin V sin cos r V sin r V r r r ⎪⎭⎫ ⎝⎛-∂∂++∂∂+⎪⎭⎫ ⎝⎛∂∂+∂∂=θθθθθθθθθθθθθcos sin V r1cos V r 1cos V r 1cos sin v V r 1cos sin r V sin r V 22r r 2r -∂∂++∂∂+∂∂+∂∂=zV V V r 1r V z V y V x V div z r r z y x ∂∂+⎪⎭⎫ ⎝⎛∂∂++∂∂=∂∂+∂∂+∂∂=∴θυθ2-6解:(1)siny x 3x V 2x -=∂∂ siny x 3y V 2y =∂∂ 0yV x V yx =∂∂+∂∂ ∴此流动满足质量守恒定律(2)siny x 3x V 2x =∂∂ siny x 3y V 2y =∂∂ 0siny x 6y V x V 2y x ≠=∂∂+∂∂ ∴此流动不满足质量守恒定律(3)V x =2rsin rxy2=θ V y =-2rsin 2ry 22-=θ33ry 2x Vx =∂∂332yr 2y y x 4y V +-=∂∂0ryx 4y V x V 32y x ≠-=∂∂+∂∂ ∴此流动不满足质量守恒方程(4)对方程x 2+y 2=常数取微分,得xdy dy dx -= 由流线方程yx v dy v dx =(1) 由)(得2r k v v r k v 422y 2x =+=由(1)(2)得方程3x r ky v ±= 3yr kxv = 25x r kxy3x V =∂∂∴25y rkxy 3yV ±∂∂0yV x V yx =∂∂+∂∂ ∴此流动满足质量守恒方程2—7解:0xVz V 0r yz 23r yz 23z V y V z x 2727y z =∂∂-∂∂=⋅+⋅-=∂∂-∂∂同样 0y V x V x y =∂∂-∂∂ ∴该流场无旋()()()2322222223222z y x z y x z y x d 21zy xzdzydy xdx dz v dy v dx v d ++++⋅=++++=++=Φ c zy x 1222+++-=Φ∴2—8解:(1)a x V x x =∂∂=θ a yV y y =∂∂=θ a z Vz z -=∂∂=θ 021v ;021v ;021v z y x =⎪⎪⎭⎫ ⎝⎛∂∂+∂∂==⎪⎭⎫ ⎝⎛∂∂+∂∂==⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=y V x V x V z V z V x V x x z x y z (2)0y V x V 210x V z V 210z V y V 21x y z z x y y z x =⎪⎪⎭⎫⎝⎛∂∂-∂∂==⎪⎭⎫ ⎝⎛∂∂-∂∂==⎪⎪⎭⎫ ⎝⎛∂∂-∂∂=ωωω;; 位该流线无旋,存在速度∴ (3)azdz 2aydy axdx dz v dy v dx v d z y x -+=++=ϕc az ay 21ax 21222+-+=∴ϕ2—9解:曲线x 2y=-4,()04y x y x f 2=+=, 切向单位向量22422422y2x 2y2x yx 4x xy 2i yx 4x x j f f fx i f f fy t +-+=+-+=t t v v v t ⋅∇=⋅=∇=ϕϕ切向速度分量 把x=2,y=-1代入得()()j x 2x i y x 2x j yi x v 2+-+--=∂∂+∂∂=∇=ϕϕϕj 21i 21j y x 4x 2xy i y x 4x x t 2242242+=⎪⎪⎭⎫ ⎝⎛+-+= 23t v v t -=⋅= j 23i 23j 21i 2123t v v t t --=⎪⎭⎫⎝⎛+-==2—14解:v=180h km =50s m根据伯努利方程22V 21V 21p ρρρ+=+∞∞ pa p =∞ 驻点处v=0,表示为1531.25pa 501.22521V 21pa p 22=⨯⨯==-∞ρ相对流速为60s m 处得表 示为75.63760225.12125.1531V 21V 21pa p 222-=⨯⨯-=-=-∞ρρ第三章3—1解:根据叠加原理,流动的流函数为()xy arctg 2Q y V y x πϕ+=∞, 速度分量是22y 22x yx y2Q x V y x x 2Q V y V +⋅=∂∂-=+⋅+=∂∂=∞πϕπϕ; 驻点A 的位置由V AX =0 V Ay =0求得 0y V 2Qx A A =-=∞;π 过驻点的流线方程为2x y arctg 2y x y arctg 2y y Q V Q V A A A =+=+∞πθπ θθππθππsin 2r x y arctg 2y -⋅=⎪⎭⎫ ⎝⎛-=∞∞V V Q 或即 在半无限体上,垂直方向的速度为θπθθππ-sin v r sin 2y x y 2v 222y ∞==+=Q Q 线面求极值()0-sin v -cos sin v 2d dv 22y=+=∞∞θπθθπθθθ 当0sin =θ 0v v min y y == 2-tg -=θπθmaxyy v v =用迭代法求解2-tg -=θπθ得 取最小值时,y 1v 2183.1139760315.1==θ取最大值时,y 2v 7817.2463071538.4==θ由θπθθππ-sin v r sin 2y x y 2v 222y ∞==+=Q Q θπθθθππ-cos sin v r cos 2v y x x 2v v 22x +=+=++=∞∞∞Q Q 可计算出当∞∞===v 6891574.0v v 724611.0v x y 1,时,θθ 6891514.0v v 724611.0v x y 2=-==∞,时,θθ 合速度∞=+=v v v 2y 2x V3—3解:设点源强度为Q ,根据叠加原理,流动的函数为 xa3-y arctg2a x y arctg 2a x y arctg 2πθπθπθϕ+++-=两个速度分量为()()()⎥⎥⎦⎤⎢⎢⎣⎡+++++++--=222222a 3-y x xy a x a x y a x a x 2x πθ()()()⎥⎥⎦⎤⎢⎢⎣⎡++++++-=222222y a 3-y x a3-y y a x y y a x y 2v πθ 对于驻点,0v v y x ==,解得a 33y 0x ==A A ,3—4解:设点源的强度为Q ,点涡的强度为T ,根据叠加原理得合成流动的位函数为Q ππθϕ2lnr 2Γ+= πθϕπθϕθ2r 1r 12r 1r r Γ=∂∂==∂∂=V V ; 速度与极半径的夹角为Qarctg arctg r Γ==V V θθ3—5根据叠加原理得合成流动的流函数为⎪⎪⎭⎫ ⎝⎛+--+=∞y a y yaarctg a y y aarctg V ϕ 两个速度分量为()()()()⎥⎦⎤⎢⎣⎡++---+++=∂∂=∞1y v 2222x y a x a x a y a x a x a V ϕ()()⎥⎦⎤⎢⎣⎡+--++=∂∂-=∞2222y y v y a x yy a x y a V ϕ 由驻点()0a 30,得驻点位置为±==y x v v零流线方程为0ay yaarctg a y y xaarctg y =--++∞∞V V对上式进行改变,得⎪⎭⎫ ⎝⎛-=-+a y tan ay2a y x 222当0x =时,数值求解得a 03065.1y ±=3—9解:根据叠加原理,得合成流动的流函数为a y y arctg 2a y y arctg 2y v -++-=∞ππϕQ Q速度分量为()()2222x ya x ax 2y a x a x 2y v v +-+++++-=∞ππQ Q ()()2222y y a x ax 2y a x a x 2v +-+++++-=ππQ Q 由0v v y x ==得驻点位置为⎪⎪⎭⎫ ⎝⎛+±∞0v a a 2,πQ 过驻点的流线方程为ay yarctg 2a y y arctg 2y v =-++--∞ππQ Q 上面的流线方程可改写为ay yarctg a y y arctg y v 2--+=∞Q π 222a y x ay2a y y arctg a y y arctg tan y v 2tan -+=⎪⎪⎭⎫ ⎝⎛--+=⎪⎪⎭⎫ ⎝⎛∴∞Qπ 容易看出y=0满足上面方程当0y ≠时,包含驻点的流线方程可写为⎪⎭⎫ ⎝⎛-=-+∞Q y v 2tan ay2a y x 222π当12v a ===∞πQ 时,包含驻点的流线方程为tanyy21y x 22--=-+3—10解:偶极子位于原点,正指向和负x 轴夹角为α,其流函数为 22yx xsin ycos 2+--=ααπϕM 当45=α时 22yx xy 222+--=πϕM3—11解:圆柱表面上的速度为a2sin v 2v πθΓ--=∞ 222222a4a 2sin v 4v ππθΓ+Γ=∞ 222222v a 4av 2sin 4sin 4v v ∞∞∞Γ+Γ+=⎪⎪⎭⎫ ⎝⎛ππθθ压强分布函数为222p v asin 41sin 41v v 1⎪⎪⎭⎫ ⎝⎛Γ+-=⎪⎪⎭⎫ ⎝⎛-=∞∞θπθC第四章4—1解:查表得标准大气的粘性系数为nkg 1078.1u 5-⨯=65el 1023876.11078.16.030225.1u ⨯=⨯⨯⨯==-∞L V R ρ 平板上下两面所受的总得摩擦阻力为N S V L R F 789.021e 664.0222=⨯⨯=∞ρ 4—2解:沿边阶层的外边界,伯努利方程成立代表逆压梯度代表顺压梯度,时;当时当0m 0m 00m 00m m v v v 21p 12201002〈〉∴〉∂∂〈〈∂∂〉-=-=∂∂-=∂∂=+--xpx p x v x v x v xx p c m m m ρρρρδδδ4—4解:(a )将2x y 21y 23v v ⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛=δδδ带入(4—90)中的第二式得δδδδδ28039dy vv 1v v 0x x =⎪⎪⎭⎫ ⎝⎛-=⎰**由牛顿粘性定律δτδu u 23y v u 0y x w =⎪⎪⎭⎫ ⎝⎛∂∂==下面求动量积分关系式,因为是平板附面层0dx dv =∴δ积分关系式可表示为dxd v 2w **=δρτδ 将上述关系式代入积分关系式,得δρδδv dxud 14013=边界条件为x=0时,0=δ 积分上式,得平板边界层的厚度沿板长的变化规律()64.428039646.0x x x64.4ll ⨯==∴=**R R δδ(b )()74.164.483x x 83dy v v 1lx =⨯=∴=⎪⎪⎭⎫ ⎝⎛-=*∞*⎰R δδδδ(c )由(a )知()64.4x x l =R δ(d )646.0x x646.0v 21324xx 64.4u23l f l 2wf l w =∴====R C R C R δρτδδδτ)得—由(; (e )单面平板的摩擦阻力为()292.1x x 292.1s v 21b bdx v 21l f l 2f l02f=∴===⎰R C R X C C X F F δδρρ摩阻系数为假设版宽为4—6解:全部为层流时的附面层流厚度由式(4—92)得()01918.048.5L e ==LR Lδ 全部为湍流时的附面层流厚度由式(4—10)得()0817.037.0L 51e ==-LLR δ第五章5-1 一架低速飞机的平直机翼采用NACA2415翼型,问此翼型的f ,f x 和c 各是多少?解:此翼型的最大弯度f =2% 最大弯度位置f x =40% 最大厚度c =15%5-2 有一个小α下的平板翼型,作为近似,将其上的涡集中在41弦点上,见图。

列车空气动力学问题(2)

列车空气动力学问题(2)
• 压力脉冲的绝对值约为车速的二次方函数,同时也 是车头车尾形状的函数;
• 压力脉冲对车体侧墙特别是侧窗玻璃有较大的破坏 作用,广深线曾发生大面积侧窗玻璃破损事故;
• 隧道内的会车导致更为复杂的气流脉动。
2020/8/19
机械工程学院
8
影响会车压力脉冲的因素
• 车头车尾形状(特别是车头尾变截面长度); • 线间距(线间距越大,会车压力脉冲的影响越弱。但
线间距越大,线路占用土地越多,线路造价越昂贵) • 列车速度(绝对速度和相对速度); • 环境风速和风向。
2020/8/19
机械工程学院
9
6. 列车与路旁建筑、设备和人员间气动效应问题
• 列车高速通过时对路旁建筑、设备和人员造成两个 方面的气动作用:
1. 列车诱导滑流引起的列车风(速度较大)
列车风的速度与列车运行速度大致呈线性规律变化,与设 备(人员)距列车距离大致呈二次方函数规律变化。
4. 列车过隧道的气动效应问题
• 列车通过隧道时引起的空气动力学问题体现在两方面
1. 列车运行气动阻力增加
由于隧道中列车周围气流运动受到更多阻碍,车体的压差阻 力和摩擦阻力均有较大上升。
2. 隧道出入口产生压力波并以声速来回传递
车头进入隧道时会产生一个压缩波,以声速向隧道出口传递, 并在出口处产生一个膨胀波,并以声速向反方向传递;当车尾 进入隧道入口时会产生一个膨胀波,以声速向隧道出口方向传 递,膨胀波在隧道出口处产生一压缩波并向反方向传递。因此 压力波的运动使得隧道内空气流动状况非常复杂。
2020/8/19
机械工程学院
1
压力波的产生过程示意
列车进入隧道时产生压缩波
2020/8/19
在隧道出口处形成膨胀波

汽车空气动力学-第二章

汽车空气动力学-第二章

1 2 M R v AaC RM 2
绕y轴的纵倾力矩MP (Pitching Moment) 纵倾力矩
1 2 M P v AaC PM 2
1 2 M Y v AaCYM 2
绕z轴的横摆力矩MY(Yawing Moment) 横摆力矩
力和力矩
横摆角 0 时车身纵向作用的 气动阻力D 阻力系数
Re
vd

局部损失——弯管
对于弯管,在拐角处将发 生流动分离现象。拐角外 壁压力高而内壁压力低, 因此,在流动接近拐角外 壁时和离开拐角内壁时发 生分离现象。
分离区将随转弯半径r的减小 和转角的增加而愈加明显。
局部损失——进气口
在进气口处,流动将产生总压损失。特别是对 于尖角入口,将发生分离现象,沿程损耗系数 较大。 为得到小的损耗系数,进气口必须很好地圆化。
飞机:达到0.08
目前雨滴的风阻系数最小 :0.05左右
下面是一些物体的风阻
一般轿车风阻系数: 0.28-0.4 好些的跑车在:0.25左右 赛车可以达到:0.15左右
载货汽车 公共汽车
二轮车
0.40~0.60
0.50~0.80
0.60~0.90
新甲克虫
CD = 0.38
车进行比较,以每小时88km的时速行驶了
100km,燃油消耗后者比前者节约了1.7L。
气动阻力
压差阻力 摩擦阻力
气动阻力
诱导阻力 干扰阻力 内流阻力
压差阻力 压差阻力有时又称为形状阻力,因为它的大小直
接取决于物体的形状
压差阻力是由于运动空气的粘性导致汽车前后产
生压力差而形成的阻力。约占汽车总气动阻力的
10万以下热点车型
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

三.面元法和涡格法是基于什么方程求解,有什么异同,简述涡格法求解步骤答:面元法和涡格法都是拉普拉斯方程
2,同:a基本求解都是基于一个面上
b边界条件在控制点上是不可穿透的
c求解高维线性方程组得知每个基本解的强度
异:a涡格法强调升力,不能模拟厚度
b边界条件不一样,涡格法布置在中性面上,不在实际的面上
c基本解布置位置不一样,涡格法不是布置在整个面上
d涡格法考虑的是薄面,面元法对厚度没限制
3,涡格法求解步骤 a 对某个近似平面用四边形划分涡格 b 在每个涡格上布置马蹄涡(1/4c) c 每个涡格控制点满足不可穿透条件 d 根据每个马蹄涡的环量强度,计算每个马蹄涡的升力,然后据此计算全机的升力 4.面元法求解步骤:
将翼型上下表面打断城直线段,假定在每一线段或者每一块面内点源强度是一个常量,每个快之间是不同的值,而涡格强度对每个面内都是常量。

四,简述CFD 求解方程
答:1建立控制方程——2确立初始条件及边界条件——3划分网格,生成计算节点——4建立离散方程——5离散初始条件和边界条件——6给定求解控制参数——7求解离散方程——8判断解是否收敛(不收敛则返回4)——9显示和输出计算结果
五.什么是离散化?常用离散化的方法,各自的特点。

答:1,离散化:在对制定问题进行CFD 计算之前,首先要将计算区域离散化,即对空间上
连续计算的区域进行划分,把他划分成多个子区域,并确定每个区域中得节点,从而生成网格。

然后将控制方程在网格上离散,即将偏微分格式的控制方程转化为各个节点上的代数方程组。

对于瞬态问题,还要进行时间域的离散。

即对计算区域进行空间和时间方向的离散。

2,常用的离散化有:
a 有限差分法:直接将微分问题变成代数问题的近似数值解法,这种方法发展较早,比较成
熟,适用于求解双曲型和抛物型问题,但求解边界条件复杂、尤其是椭圆问题不如有限元或有限体积法方便。

b 有限元法:具有广泛的适应性,特别适用于集合无力条件比较复杂的问题(尤其是对椭圆
问题有更好的适用性)。

求解速度比有限差分法和有限体积法慢,故在CFD 软件里应用并不普遍。

c 有限体积法:简单地说,子域法加离散,就是有限体积法的基本方法。

特点是计算效率高,
在CFD 领域得到了广泛应用。

九.当地时间步长,多重网格,预处理的加速收敛机理? 答:当地时间步长:max max t λx ∆≈∆,ij NF
j i s c n v CFL t ∆+∙∑Ω=∆→→-)(1
多重网格:
多重网格是一种非常有效的加速收敛技术,即可用于显式格式,又可用于隐式格式。

其思想是,为了是密网格上的流场计算尽可能快地收敛到最终的定常解,同时在另外几套依次变稀的网格上做计算,稀网格的计算结果再反馈给密网格。

其加速收敛的机理是:
大多数的计算是在稀网格上进行的,可取较大的时间步长,而且计算量较小。

(收敛快,计算机时少)
绝大部分显式或隐式时间推进和迭代求解方法降低高频误差有效,对降低低频误差效果很差。

计算域内所有频段的误差都得到降低才能达到最终的定常解(一般地,对一给定的网格,经过若干迭代步,可以很快消除掉高频误差,而低频误差则需要更多的迭代步数)。

多重网
格正是在快速消除高低频误差这一点上有很大帮助:密网格上的低频误差相当于稀网格上的高频误差,所以在各自不同密度的网格上快速降低各自的高频误差,相当于同时降低了密网格上从高到低的所有频率的误差。

相关文档
最新文档