《飞机空气动力学.pptx》由会员分享，可在线阅读，更多相关《飞机空气动力学.pptx（122页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、高速飞行的特点高速飞行的特点v 激波阻力（波阻）激波阻力（波阻）v 声障（音障）声障（音障）v 低速、亚音速和超音速流动的区别低速、亚音速和超音速流动的区别第第7 7章章 高速可压流动基础高速可压流动基础第1页/共122页7.1 7.1 热力学基础知识热力学基础知识热力学第一定律：内能和焓热力学第一定律：内能和焓 热力学第二定律：熵热力学第二定律：熵 热力学的物系热力学的物系 气体的状态方程，完全气体和真实气体气体的状态方程，完全气体和真实气体第2页/共122页热力学的物系2 物系与外界关系:隔热体系:无物质交换,无能量交换;封闭体系:无物质交换,有能量交换;开放体系:有物质交换,有能量交换.
2、1 热力学体系:用热力学去处理的客体是和周围环境的其他物体划分开的一个任意形态的物理体系(物系).).这个体系的尺寸是宏观的.高速流中遇到的情况，绝大多数属于隔绝体系和封闭体系。第第7 7章章 高速可压流动基础高速可压流动基础第3页/共122页热力学第一定律:内能和焓1、状态方程与完全气体假设热力学指出：任何气体的压强、密度、绝对温度不是独立的，三者之间存在一定的关系。函数称为状态方程。该方程的具体表达形式与介质种类、温度、压强的不同有关。一个物系的压强、密度、温度都是点的函数，彼此之间存在一定的函数关系，但和变化过程无关，代表一个热力学状态。p,T,r r,u，h代表热力学状态参数，两个热力
3、学参数可以确定一个热力状态，如果取自变量为T,r r，其它状态变量关系为:7.1 7.1 热力学基础热力学基础第4页/共122页热力学第一定律:内能和焓2、内能、焓气体内能是指分子微观热运动（与温度有关）所包含的动能与分子之间存在作用力而形成分子相互作用的内部位能之和。对于完全气体而言，分子之间无作用力，单位质量气体的内能u仅仅是温度的函数。在热力学中，常常引入另外一个代表热含量的参数h（焓）：由于表示单位质量流体所具有的压能，故焓h表示单位质量流体所具有的内能和压能之和。焓的微分:表示气体焓的增量等于内能增量、气体膨胀功与压强差所做的功之和。7.1 7.1 热力学基础热力学基础第5页/共12
4、2页热力学第一定律:内能和焓3、热力学第一定律热力学第一定律是能量守恒定律在热力学上的具体应用。其物理意义是：外界传给一个封闭物质系统的热量等于该封闭系统内能的增量与系统对外界所做机械功之和。对于一个微小变化过程，有这是静止物系的热力学第一定律。其中，dV表示物系的体积变量，p表示物系的压强。如果用物系的质量去除上式，就变成单位质量的能量方程。单位质量流体的能量方程：其中，密度的倒数是单位质量的体积。表示外界传给单位质量流体的热量dq等于单位质量流体内能的增量与压强所做的单位质量流体的膨胀功。7.1 7.1 热力学基础热力学基础第6页/共122页热力学第一定律流动物系的能量守恒定律:(:(绝热
5、过程:dq=0):dq=0)与静止物系的能量方程相比，流动物系的能量方程多了两项，其中一项是表示流体微团在体积不变的情况下，由于压强变化引起的功（流体质点克服压差所做的功）；另一项是流体微团的宏观动能变化量。即：7.1 7.1 热力学基础热力学基础第7页/共122页4、热力学过程（1）可逆与不可逆过程在热力学中，如果将变化过程一步一步倒回去，物系的一切热力学参数都回到初始状态，且外界状态也都复旧，这样的过程则是可逆过程，否则是不可逆过程。（如高温向低温传热，机械功通过摩擦生热都是不可逆过程）可逆过程也称为准静态过程，或连续的平衡态过程。（2）绝热过程与外界完全没有热量交换，即dq=0，称为绝热
6、过程。（3）等容过程、等压过程、等温过程、绝热过程在热力学中，内能u是状态的函数，而q不是状态函数。因为其中的压力膨胀功不仅决定于过程的起点和终点，与变化过程有关。7.1 7.1 热力学基础热力学基础第8页/共122页等容过程、等压过程、等温过程、绝热过程1）等容过程如果在变化过程中，单位质量气体的容积保持不变的过程称为等容过程。此时气体的膨胀功为零。外界加入的热量全部用来增加介质的内能，即：比热定义：单位质量介质温度每升高一度所需要的热量。比热（比热容）数值的大小与具体热力学过程有关。在等容过程中，比热称为等容比热，用Cv表示。7.1 7.1 热力学基础热力学基础第9页/共122页等容过程、
7、等压过程、等温过程、绝热过程2）等压过程如果在变化过程中，气体的压强保持不变的过程称为等压过程。此时气体的膨胀功不等于零。外界加入的热量一部分用来增加介质的内能，另一部分用于气体的膨胀功。在等压过程中，单位质量介质的温度每升高一度，所需要的热量，称为定压比热，用Cp表示：定压比热与定容比热的比值，称为气体的比热比。即：在空气动力学中，在温度小于300C，压强不高的情况下，一般Cp，Cv，g g等于常数。对于水7.1 7.1 热力学基础热力学基础第10页/共122页等容过程、等压过程、等温过程、绝热过程3）等温过程在变化过程中，气体的温度保持不变的过程称为等温过程。在等温过程中，内能不变，热量与
8、膨胀功相等。单位质量气体所做的功为7.1 7.1 热力学基础热力学基础第11页/共122页等容过程、等压过程、等温过程、绝热过程4）绝热过程在热力学变化过程中，与外界完全没有热量交换。由能量方程得到：在由理想气体的状态方程，有：内能的变化为：7.1 7.1 热力学基础热力学基础第12页/共122页定容过程的比定容热容cv:内能的改变量为:du=cvdT气体作等压变化时,p=常数,dp=0:焓的变化量:比热容:物系的温度每升高1所需的热量.气体在定容变化的过程中,体积不变,1/=常数.等容过程、等压过程、等温过程、绝热过程7.1 7.1 热力学基础热力学基础第13页/共122页热力学第二定律，熵
9、通过引入熵状态参数，在不可逆过程中的变化来描述热力学第二定律。熵是一个热能可利用部分的指标。其定义如下：单位质量气体的熵定义为：其中，dqdq与dq/Tdq/T是不同的两个量。dqdq是与积分路径有关的；而dq/Tdq/T是一个与积分路径无关的量，可以表示成某一函数的全微分：在研究热力学过程中，最有意义的是熵的增量，即从状态1到状态2的熵增。即：7.1 7.1 热力学基础热力学基础第14页/共122页熵熵:热力学参数,是状态参数,和物系的具体变化过程无关.可逆过程:有(ds)(ds)inin=0;=0;不可逆过程:有(ds)(ds)inin0.0.等熵流动:流动变化过程是可逆的,则(ds)(d
10、s)exex和(ds)(ds)inin都为0,0,介质的熵没有变化的流动.一般在绕流场的绝大部分区域速度梯度和温度梯度都不大，流场可近似视为绝热可逆，熵值不变ds=0，称为等熵流动，一条流线熵值不变叫做沿流线等熵，在全流场中熵值不变，称为均熵流场。在边界层及其后的尾迹区，激波传过的流动，气体的粘性和热传导不能忽视区，流动是熵增不可逆过程ds0，等熵关系式不能用。7.1 7.1 热力学基础热力学基础第15页/共122页气体的状态方程,完全气体和真实气体质量定压热容:与比热比的关系:其中:空气 质量定容热容:气体的状态方程:p/p/=RT=RT7.1 7.1 热力学基础热力学基础第16页/共122
11、页完全气体等熵过程完全气体等熵过程的压强比对温度比的关系:在等熵流动中，有：称为等熵关系，g g为等熵指数。7.1 7.1 热力学基础热力学基础第17页/共122页7.2 7.2 声速和马赫数声速和马赫数微弱扰动传播过程与传播速度一声速微弱扰动传播过程与传播速度一声速 声速公式声速公式 现现 象象 马赫数马赫数第18页/共122页 现象Mappa尾迹在微小扰动下，介质的受绕速度也是微小的，但微小扰动的传播速度则是一定的，其值与介质的弹性和质量有关，与扰动的振幅无关。空气是一种弹性介质，在这种介质中任何一个微小扰动都会向四面传播出去，当然传播速度决定于介质的状态。不可压流中,微弱扰动的传播速度为
12、无限大。可压流中,扰动不会瞬间传遍整个流场，传播速度为一定数值。弱扰动（不可压流）:使流动参数的数值改变得非常微小的扰动强扰动（可压流）:使流动参数改变有限值的扰动第第7 7章章 高速可压流动基础高速可压流动基础第19页/共122页微弱扰动传播过程与传播速度声速V=0V=0V=00+dv+dp+dp0+dvp+dp+dpxxx声速:微弱扰动在介质中的传播速度。声速以球面波的形式传播。波：受到扰动的气体与未受到扰动的气体之间的分界面。声音以波的形式传播，声波是一种微弱的气体扰动波运动。微小扰动在弹性介质中的传递是以压力波的形式传播的，其传播速度（声速）的大小与介质的弹性存在密切的关系。假定有一根
13、十分长的管子，管子左端有一个活塞。现将活塞以微小速度dv向右推动，使管内空气产生一个压缩的微小扰动。7.27.2声速和马赫数声速和马赫数第20页/共122页这个扰动将以一定的波速a向右传播，在管道中扰动以波阵面A-A的形式向右推进。在波阵面右侧的气体未受扰动，其压强、密度、温度和速度分别为：p、r r、T、v=0；而在波阵面左侧的气体受到扰动后，其压强、密度、温度和速度分别变为：p+dp、r r+dr r、T+dT、dv。扰动的传播速度与由扰动引起介质本身的运动速度是不同的。扰动传播速度要比由扰动引起介质本身的运动速度大得多。微弱扰动传播过程7.27.2声速和马赫数声速和马赫数第21页/共12
14、2页7.2.3 声速公式由于扰动是微小的，因此有为便于分析，现采用一个相对坐标，观察者跟随波阵面一起运动，这时整个流动问题由原来非定常问题变成一个定常问题。这时波阵面不动，未扰气体以波速a向左运动，气流不断越过A-A面进入扰动区，而受扰气流以a-dv速度相对于A-A面向左流去。现围绕A-A面取一控制体，由质量守衡方程得到7.27.2声速和马赫数声速和马赫数第22页/共122页T T+dT Tdvdv p+dp p v=0v=0 +d+d A Acdtcdtdvdtdvdtc c声速公式 由动量定理得到联解可得这就是声速的微分形式公式。说明气体扰动的传播速度决定于变化过程中气体的dp和dr r的
15、比值。由于扰动变化微小、速度很快，气体既无热量交换，也无摩擦产生，可认为是一种绝热等熵过程，此时压力密度关系为：7.27.2声速和马赫数声速和马赫数第23页/共122页空气声速 空气绝热指数=1.4,=1.4,声速:c=20.1T c=20.1T1/21/2 m/s m/s 在非均匀的流场中,不同时刻,不同点上声速大小和当地的温度有关,温度越高,声速越大。声速是随着高度增大而减小。对于海平面标准大气，T=288.15K，g g=1.4，得到:对于水体而言：7.27.2声速和马赫数声速和马赫数第24页/共122页马赫数Ma数表示气流运动速度V与当地音速a之比。Ma=V/a。是一个表征流场压缩性大
16、小的无量纲参数，是高速空气动力学中的一个重要基本物理参数，反映流场压缩性大小的相似准则。衡量气体压缩程度大小的可用相对密度变化来表示，而这个相对密度变化量又与Ma数的大小存在密切的关系。说明，Ma数越大气体的压缩性越大。当Ma0，则有，说明沿着流动方向，虽然总温T0不变，但总压下降。对于一维等熵流动，在流线上任意点处的总温和总压均相等。7.3 7.3 高速一维定常流高速一维定常流第34页/共122页2.2.使用临界参考量的参数关系式n临界速度等于当地声速。V V*=c=c*=(p=(p*/*)1/21/2 =(RT=(RT*)1/21/2 =(-1)h=(-1)h*)1/21/2临界参数：临界
17、状态下的气体状态参数*、p p*、T T*、h h*。【定义】临界状态：在理想气体定常等熵流动中，流体质点速度等于当地声速(Ma=1)(Ma=1)的状态。在一维绝热流动中，沿流线某点处的流速正好等于当地声速（Ma=1），该点称为临界点或临界断面。7.3 7.3 高速一维定常流高速一维定常流第35页/共122页参数关系式n由一维绝热等熵流能量方程可得：7.3 7.3 高速一维定常流高速一维定常流第36页/共122页速度系数由一维绝热流能量方程可得：定义速度系数为：流体速度与临界速度(或临界声速)之比。由于临界点的音速a*仅是总温的函数，速度系数引入的最大好处是：在给定总温下其分母是常数，因此对速
18、度系数的各种运算只对分子就行了。7.3 7.3 高速一维定常流高速一维定常流第37页/共122页速度系数与马赫数关系n速度系数与马赫数关系：速度系数的最大值为在Ma小于1，速度系数大于Ma数；在Ma数大于1，速度系数小于Ma数。7.3 7.3 高速一维定常流高速一维定常流第38页/共122页一维等熵流总静参数比n一维等熵流总静参数比:函数随速度系数的变化曲线7.3 7.3 高速一维定常流高速一维定常流第39页/共122页例7.37.3n 飞机在h=5000m,h=5000m,以Ma=0.8Ma=0.8飞行,进气口截面A A1 1=0.5m=0.5m2 2，MaMa1 1=0.4;=0.4;出口
19、截面MaMa2 2=0.2.=0.2.求来流的总参数和进口截面处的p p1 1,1 1,T,T1 1和 质量流量 。n【解】由标准大气表,按h=5000mh=5000m查得 p ph h=54020N/m=54020N/m2 2,h h=0.73612kg/m=0.73612kg/m3 3,T,Th h=255.65K=255.65K 由Ma=0.8Ma=0.8查等熵流表或计算得：p p/p/p0 0=p=ph h/p/p0 0=0.656,=0.656,/0 0=h h/0 0=0.74,=0.74,T T/T/T0 0=T=Th h/T/T0 0=0.8865=0.8865 得:p:p0
20、0=82347.6N/m=82347.6N/m2 2,0 0=0.99476kg/m=0.99476kg/m3 3,T T0 0=288.36K=288.36K7.3 7.3 高速一维定常流高速一维定常流第40页/共122页由MaMa1 1=0.4=0.4查表或计算得：p p1 1/p/p0 0=0.8956,=0.8956,1 1/0 0=0.9243,T=0.9243,T1 1/T/T0 0=0.969=0.969 p p1 1=73750.5N/m=73750.5N/m2 2,1 1=0.91946kg/m=0.91946kg/m3 3,T,T1 1=279.44K,=279.44K,c
21、 c1 1=335.1m/s=335.1m/s V V1 1=Ma=Ma1 1c c1 1=0.4335.1=134.033m/s=0.4335.1=134.033m/s =p =p1 1V V1 1A A1 1=0.91946134.0330.5=61.62kg/s=0.91946134.0330.5=61.62kg/s例7.37.37.3 7.3 高速一维定常流高速一维定常流第41页/共122页3.等熵管流的速度与面积关系，拉瓦尔管喉道 亚声速段 超声速段 Ma1 1，Ma=1 =1 Ma1 1 超声速段 亚声速段 Ma1 1，Ma=1 =1 Ma1 1 n亚声速（M Ma a1 1）：d
22、AdA与dVdV异号，dA0dA0，dV0dV0；dA0dA0dV0。n超声速（M Ma a1 1）：dAdA与dVdV同号，dA0dA0，dV0dV0；dA0dA0，dV0dV0。n声速（Ma=1Ma=1）：dA/A=0dA/A=0，A A出现极值n连续方程的微分形式：d/+dA/A+dV/V=0d/+dA/A+dV/V=0n动量方程的微分形式：dp+VdV=0 dp+VdV=0 n得 d d/+Ma+Ma2 2dV/V=0dV/V=0n则：(Ma(Ma2 2-1)dV/V=dA/A-1)dV/V=dA/A要使气流从亚声速加速到超声速(或超声速等熵地减速到亚声速)，管道形状应该先收缩后扩张。
23、7.3 7.3 高速一维定常流高速一维定常流第42页/共122页收缩喷管与拉伐尔喷管比较收缩喷管的流道截面积是逐渐缩小的，在喷管进出口压强差的作用下，高温气体的内能转变成动能，产生很大的推力。气流速度达到音速后便不能再增大了。拉伐尔喷管即是缩放式喷管，其流道先缩小再扩大，允许气流在喉道处达到音速后进一步加速成超音速流。7.3 7.3 高速一维定常流高速一维定常流第43页/共122页低速、亚音速和超音速流动的区别低速、亚音速和超音速流动的区别低速流动（低速流动（0.3-0.4Ma）流速增加流速增加静压减小静压减小流速减小流速减小静压增加静压增加（1）对亚音速(包括低速_流动，如果管道截面收缩则流
24、速增加，面积扩大流速下降；亚音速流动（亚音速流动（0.40.85Ma）流速增加流速增加静压减小静压减小密度减小密度减小温度下降温度下降声速下降声速下降马赫数增加马赫数增加流速减小流速减小静压增加静压增加密度增加密度增加温度上升温度上升声速上升声速上升马赫数减小马赫数减小7.3 7.3 高速一维定常流高速一维定常流第44页/共122页超音速流动（超音速流动（1Ma）密度减小密度减小流速增加流速增加静压减小静压减小温度下降温度下降声速下降声速下降马赫数增加马赫数增加密度增加密度增加流速减小流速减小静压增加静压增加温度上升温度上升声速上升声速上升马赫数减小马赫数减小（2）对超音速（包括低速）流动，如
25、果管道截面收缩则流速减小，面积扩大流速增加；低速、亚音速和超音速流动的区别低速、亚音速和超音速流动的区别（3）造成超音速截面流速与截面积变化规律与亚音速相反，其原因是：密度变化对连续方程的贡献。亚音速时密度变化较速度变化为慢，而超音速时密度变化比流速变化快。要想增加流速，亚音速时截面积应缩小，超音速时截面积应放大。7.3 7.3 高速一维定常流高速一维定常流第45页/共122页拉瓦尔喷管拉瓦尔喷管(喷管)Ma 1拉瓦尔管：管道形状为先收缩后扩张,中间为最小截面(喉道)。对一维等熵管流，要想让气流沿管轴线连续地从亚声速加速到超声速，即始终保持dV0，则管道应先收缩后扩张，中间为最小截面，即喉道。
26、一个喷管在出口截面产生Ma1的超声速气流的条件如下：(1)管道形状应成为先收缩后扩张的拉瓦尔管形状；(2)在喷管上下游配合足够大的压强比。一个出口接大气的喷管，当喷管出口达到设计M数而出口压强恰等于外界大气压强时，则喷管处于设计状态，而大于1的上下游压强比（即上游总压与出口大气反压之比。则为设计压强比。如果上游压强过高或过低，喷管出口内外将出现激波或膨胀波。喉道 7.3 7.3 高速一维定常流高速一维定常流第46页/共122页质量流量喉道 亚声速段 超声速段 Ma1 1，Ma=1 =1 Ma1 1 超声速段 亚声速段 Ma1 1，Ma=1 =1 Ma1 1 A A*【质量流量】对于一维定常等熵
27、管流，流过各截面的流量是一定的，用质量流量表示。【堵塞流量】7.3 7.3 高速一维定常流高速一维定常流第47页/共122页例7.47.4n某 涡 轮 喷 气 发 动 机 喷 管 进 口 燃 气 总 压 p p*=2.310=2.3105 5P Pa a,总 温T T*=928.5K,k=1.33,A=928.5K,k=1.33,Ae e=0.1675=0.1675,大气压p pa a=0.98710=0.987105 5P Pa a,求喷管出口燃气速度和压强及通过喷管的燃气流量。n解：大气压是反压,则 p pb b/p/p*=p=pa a/p/p*=0.98710=0.987105 5/2.
28、310/2.3105 5=0.429=0.429；由k=1.33,k=1.33,得：crcr=0.54=0.54 故：p pb b/p/p*crcr，是超临界流动状态。有在出口截面：M Maeae=1=1，p pe e=crcrp p*=0.542.310=0.542.3105 5=1.24210=1.242105 5 P Pa a 气流速度:V:Ve e=c=ce e=2kRT=2kRT*/（k+1k+1）1/2 1/2 =21.33287.4928.5/(1.33+1)=21.33287.4928.5/(1.33+1)1/21/2=552m/s=552m/s 燃气流量:m:me e=K p
29、=K pe e*A Ae eq(q(e e)/(T)/(Te e*)1/21/2 =0.03972.310=0.03972.3105 50.1675/(928.5)0.1675/(928.5)1/21/2=50.3kg/s=50.3kg/s7.3 7.3 高速一维定常流高速一维定常流第48页/共122页 扰动源在静止的空气中以速度扰动源在静止的空气中以速度v作等速直线运动，根据扰动作等速直线运动，根据扰动源的不同运动速度，会出现四种可能的情况：源的不同运动速度，会出现四种可能的情况：扰动源静止不动：扰动源静止不动：M0扰动源以亚音速运动：扰动源以亚音速运动：0 M 1扰动源以等音速运动：扰动源
30、以等音速运动：M 1扰动源以超音速运动：扰动源以超音速运动：M 17.4 微弱扰动的传播区物体在静止空气中运动时，不同的运动速度其对空气的影响范围、影响方式是不同的。扰动是指引起气流发生速度、密度、压强等变化的。对于亚声速流场和超声速流场而言，扰动的传播和范围是不同的。第第7 7章章 高速可压流动基础高速可压流动基础第49页/共122页由于扰动源静止不动，所以扰动波以音速由于扰动源静止不动，所以扰动波以音速a向四周传播，形成向四周传播，形成以扰动源为中心的以扰动源为中心的同心球面波同心球面波。1、静止气体（Ma=0），V=0从某瞬间看，前i秒发出的扰动波面是以扰源O为中心、i为半径的同心球面。
31、只要时间足够长，空间任一点均会受到扰源的影响，即扰源的影响区是全流场。7.4 7.4 微弱扰动的传播区微弱扰动的传播区第50页/共122页由于扰动源以亚音速运动，所以扰动源总是落后于扰动波，由于扰动源以亚音速运动，所以扰动源总是落后于扰动波，形成偏向扰动源前进方向的形成偏向扰动源前进方向的不同心球面波不同心球面波。2、亚声速气流(Ma1)，V分界面7.4 7.4 微弱扰动的传播区微弱扰动的传播区第53页/共122页该锥面称为马赫锥，马赫该锥面称为马赫锥，马赫 锥锥 的半顶角称为马赫角的半顶角称为马赫角。显然，显然，M数越大，马赫数越大，马赫 锥就越尖锐。锥就越尖锐。超声速气流受到微小扰动后，将
32、以声速向四周传播出去，把扰动球面波包络面，称为扰动界面，也称为马赫波阵面，简称马赫波。在马赫波上游，气流未受影响，在马赫波的下游气流受到扰动影响。超声速气流(Ma1)，V分界面马赫角大小为：7.4 7.4 微弱扰动的传播区微弱扰动的传播区第54页/共122页P点的依赖域 PP点的影响域 扰动的依赖域：空间固定点P P能够接收到气流扰动信号的区域。扰动的依赖域亚声速和超声速流场微弱扰动的传播区(或影响区),),不同。超声速流场与亚声速流场主要差别：影响域和依赖域。不可压流场和亚声速流场的影响域和依赖域是全流场；超声速流场的影响域只限于扰动下游马赫锥内，依赖域在倒马赫锥内。超声速流场的影响域超声速
33、流场的依赖域7.4 7.4 微弱扰动的传播区微弱扰动的传播区第55页/共122页7.5 7.5 膨胀波膨胀波外折外折诸参数的变化趋势诸参数的变化趋势膨胀波的反射和相交膨胀波的反射和相交 超声速流绕外钝角膨胀的计算超声速流绕外钝角膨胀的计算关于微弱扰动传播区的回顾关于微弱扰动传播区的回顾 壁面外折壁面外折d d特征线法特征线法 平面无旋流的特征线法平面无旋流的特征线法第56页/共122页7.5 膨胀波关于微弱扰动传播区的回顾MaMa马赫锥MaMaoo马赫锥内的气流参数及流动方向与未受扰动气流相同。对于压强和密度存在升高的变化过程，称为压缩过程；对于压强与密度存在下降的过程，称为膨胀过程。在超声速
34、流动中，压缩和膨胀过程都是有扰界的，称为波阵面。波阵面膨胀波(或马赫线)：超声速气流因通路扩张(如壁面外折一角度),),或流动从高压区过度到低压区,气流要加速、降压,将出现膨胀波。=arcsin(1/M=arcsin(1/Ma a)马赫数：马赫角第第7 7章章 高速可压流动基础高速可压流动基础第57页/共122页壁面外折d d若在O点处壁面向外折转一个微小的角度d，使流动区域扩大。则O点是一个微小扰动源，扰动的传播范围是在O点发出的马赫波OL的下游，扰动影响的结果是，使气流也外折一个d同样大小的角度。壁面外折，相当于放宽气流的通道。对超声速气流来说，加大通道截面积必使气流速度增加，压力和密度下
35、降，气流发生膨胀。此时，马赫波线OL的作用是使超音速气流加速减压的，气流发生绝热加速膨胀过程，于是把马赫波OL称为膨胀波。n超声速气流绕凸角流动得到激波后的压强小于激波前的压强，即负转角的斜激波是膨胀过程。7.5 7.5 膨胀波膨胀波第58页/共122页壁面外折d d对于多个微小外偏角情况：在o1点，壁面外偏d1，通过膨胀波OL1在o2点，壁面外偏d2，通过膨胀波OL2第一道膨胀波与来流方向之间的夹角为：第二道膨胀波与来流方向之间的夹角为：7.5 7.5 膨胀波膨胀波第59页/共122页壁面外折d d由于气流发生膨胀，Ma2Ma1，则有：所以21，即第二道膨胀波与波前气流方向的夹角小于第一道膨
36、胀波的倾斜角。由于气流发生膨胀，Ma2Ma1，则有：所以21，即第二道膨胀波与波前气流方向的夹角小于第一道膨胀波的倾斜角。由于后产生的每一道膨胀波相对于原始气流的倾斜角都比前面的小，所以每道膨胀波不可能出现彼此相交的情况，因而形成一个连续的膨胀区域。7.5 7.5 膨胀波膨胀波第60页/共122页Ma1 O1 Ma4 Ma3 Ma2 O3 O2 O1 L1 Ln L1 L3 L2 L2 L3 1 1 1 d1 Ma2 Ma1 1 d2 d3 外折n总折角：若折点无限靠近，这些马赫波集中于一点，组成以这点为中心的扇形膨胀波束称为膨胀波。普朗特-迈耶流动：超声速气流绕外钝角的流动,在折点处产生一束
37、膨胀波。曲线可以看作是无数条微元折线的极限。超声速气流绕外凸曲壁膨胀过程情况和上面的分析完全一样，只是各道膨胀波是连成一片的连续膨胀带。7.5 7.5 膨胀波膨胀波第61页/共122页外折曲线可以看作是无数条微元折线的极限。超声速气流绕外凸曲壁膨胀过程情况和上面的分析完全一样，只是各道膨胀波是连成一片的连续膨胀带。n总折角：流线在CD段是直线，在DE段是曲线，在E之后是直线，气流完成了转折。Ma2大于Ma1。如果扰一个钝外角的流动，这时相当于壁面的外折点重合，整个马赫波形成一个扇形膨胀区，也叫膨胀波。（普朗特迈耶（Prandtl-Meyer）流动）普朗特迈耶（Prandtl-Meyer）流动（
38、绕外钝角的流动）7.5 7.5 膨胀波膨胀波第62页/共122页诸参数的变化趋势经过膨胀波以后，气流参数的变化趋势:流速V是不断增大的，dV0，因此有：由微分形式的动量方程：压强p必减小，dp0。由绝热流动的能量方程：温度T必减小，dTmax,气流 在 转 过max后 不 再贴着物面流动,而与物面分离。7.5 7.5 膨胀波膨胀波第72页/共122页2.2.流线向径r r与当地速度系数的关系式流线=1 y y r*r膨胀区中的流线方程：7.5 7.5 膨胀波膨胀波第73页/共122页3.3.数值表n角是膨胀区中的任意一道膨胀波与=1=1时流线的垂线之间的夹角。流线=1 y y r*r数值表是从
39、=1开始算起，以气流折角为自变量，给定一系列的值，算出与各个相对应的，Ma。又因膨胀过程是等熵过程，与每个相对应的，亦都列在表中。7.5 7.5 膨胀波膨胀波第74页/共122页例7.67.6n已知1 1=1=1的气流(=1.4)=1.4)绕外钝角折转10,p10,p1 1=101.325=101.325 kPa,kPa,求膨胀结束后气流的和p,p,并求通过r r*=20mm=20mm处的流线。n解：查数值表得:=10=10时,2 2=1.323,p=1.323,p2 2/p/p0 0=0.299=0.299 而p p1 1/p/p0 0=0.528=0.528得:p p2 2=(p=(p2
40、2/p/p0 0)(p)(p0 0/p/p1 1)=0.299/0.528=57.248 kPa)=0.299/0.528=57.248 kPa 查得 2 2=5550,=5550,2 2=4410=4410流线=10=1L2 r*=20mmrL12 2=55507.5 7.5 膨胀波膨胀波第75页/共122页7.6 7.6 激激 波波激波的反射与相交激波的反射与相交 圆锥激波圆锥激波 正激波正激波 斜激波斜激波 收敛收敛扩张喷管在非设计状态下的工作扩张喷管在非设计状态下的工作第76页/共122页 当飞机以等音速或超音速飞行时，在其前面也会出现由无数当飞机以等音速或超音速飞行时，在其前面也会出
41、现由无数较强的波迭聚而成的波面，这个波面就称为较强的波迭聚而成的波面，这个波面就称为激波激波。7.6 激波膨胀波是使气流发生膨胀的扰动波，而激波是以一定强度使气流发生突然压缩的波。当气体以超音速运动时，扰动来不及传到前面去，路途上的气体微团没有事先的准备，要等到物体冲到跟前，才受到压缩，因而可以造成大块气体被压缩第第7 7章章 高速可压流动基础高速可压流动基础第77页/共122页激波分类：激波分类：正激波：波面正激波：波面 与飞行速度垂直。与飞行速度垂直。斜激波：波面相对于飞行速度有倾斜角。斜激波：波面相对于飞行速度有倾斜角。激波及其分类激波及其分类7.6 7.6 激波激波第78页/共122页
42、正激波激波是很薄的一层（厚度为分子自由程的量级），物理量（速度，温度，压强）迅速地从波前值变化到波后值，梯度很大。由于激波厚度相对于流体的宏观运动非常薄，故忽略激波厚度而将激波简化成绝热的间断面。定义和气流速度垂直的物理量间断面为正激波。激波特性：激波是一层受到强烈压缩的空气层。激波是一层受到强烈压缩的空气层。气流通过激波时，压强、密度、温度突然增加，而速度却大大气流通过激波时，压强、密度、温度突然增加，而速度却大大降低。降低。激波是很薄的、具有强粘性的区域。通过激波流动是绝热的但不等熵7.6 7.6 激波激波第79页/共122页以右图活塞在一维长管中压缩为例:有一根很长的直管，管内气体原是静
43、止的热力学参数是p1,1,T1 从t=0起到t=t1为止活塞向右作急剧地加速运动，t=t1 以后匀速前进特征：居后的波比前边的波快，每道波都在追赶它前面的波过渡区AA-BB的长度随时间增长而越来越短，最后压缩到一起形成激波1正激波的形成7.6 7.6 激波激波第80页/共122页2 2 正激波的传播速度当由无数个微小压缩波叠加在一起形成激波后，其波阵面以一定的速度向右推进，现在利用积分形式的控制方程，推求激波推进速度Vs。取如图所示的控制体，设激波在初始时刻位于2-2面，在t时段，激波由2-2推进到1-1面，设控制体的长度为x=Vs t根据积分形式的控制方程来推导。设激波推进速度为Vs，激波后
44、气体的运动速度为Vg，得:7.6 7.6 激波激波第81页/共122页正激波的传播速度可得到：如果规定了激波的强度p2/p1，就可以求出激波推进速度了。由此得出，激波的推进速度总是大于微小扰动波的传播速度a1。令：这说明，激波的推进速度相对于波前气体而言必是超声速的。激波对于波后已经有Vg运动速度的气体而言，其相对速度必是亚声速的。即：7.6 7.6 激波激波第82页/共122页兰金-于戈尼奥关系式等熵过程的压强与密度关系激波的传播速度正激波的传播速度7.6 7.6 激波激波第83页/共122页激波传播速度与波前气体传播速度关系:n激波的推进速度相对于波前气体而言必是超声速的。n激波对于波后已
45、经有Vg运动速度的气体而言，其相对速度必是亚声速的。即：Vs-V气c2激波传播速度与波后气体传播速度关系:Vs-V气1时，熵增量总是正的；而当M11之下，如果机翼前缘尖劈的顶角2不太大，所形成的上下两道简单的直激波，其波面和运动方向成一定的斜角，激波依附在物体的尖端上。这种激波在形式上与正激波不同，我们把这种波阵面与来流方向斜交的激波称为斜激波。n在斜激波中，激波波阵面与来流方向之间的夹角,称之为激波斜角或简称为激波角。7.6 7.6 激波激波第93页/共122页激波面 Ma11 1 V V1 1 A AV V2 2 V V1 1 V1nV2n+SV1VtV1nV2VtV2nVtS=SSS物面
46、 (1)(1)物理图画n斜激波后的气流方向也不与激波面垂直，与波前气流方向也不平行，而是与尖劈面平行，夹角，称为气流折角,意指气流经过斜激波后所折转的角度。7.6 7.6 激波激波第94页/共122页 (2)(2)波前波后气流参数的关系n如图所示，现在斜激波波阵面上取一段12341作为控制体，其中12面、34面都平行于波阵面，且二者靠的很近。按照波阵面的方向将速度分解为与波阵面垂直和平行的分量。12面：来流速度为V1，分量为V1t、V1n；34面：合速度为V2，分量为V2t、V2n。利用积分形式的质量方程得到：然后计算切向动量关系。由切向的动量积分方程得到（在14和23面上无压差）：由此得到7
47、.6 7.6 激波激波第95页/共122页激波 Ma11 1 V V1 1 A AV2 V V1 1 斜激波的相容条件（基本方程）n去垂直于激波面的控制体。得斜激波的相容条件：n质量守恒方程：1 1V V1 n 1 n =2 2V V2 n2 n n动量方程：n法向 p p1 1+1 1V V2 21 n1 n =p=p2 2+2 2V V2 22 n2 n n 1n1n2n2n=1=1n切向 1 1V V1n1nV V1 t1 t =2 2V V2n2nV V2t2t n即 V Vt t =V=V1t1t =V=V2t2t n（切向速度无变化）n能量守恒方程：n V V1n1n2 2/2+V
48、/2+V1t1t2 2/2+h/2+h1 1=V=V2n2n2 2/2+V/2+V2t2t2 2/2+h/2+h2 2 n状态方程：p p1 1/(T/(T1 11 1)=p)=p2 2/(T/(T2 22 2)7.6 7.6 激波激波第96页/共122页1)1)与正激波相同者1)1)兰金-于戈尼奥关系式:2)2)总温(总温不变)3)3)熵增量:7.6 7.6 激波激波第97页/共122页2)2)与正激波不同者1)1)压强比与M Ma1a1数关系:2)2)密度比与M Ma1a1数关系:3)3)温度比与M Ma1a1数关系4)4)总压强比与M Ma1a1数关系7.6 7.6 激波激波第98页/共
49、122页(3)(3)激波图线及应用1）对于给定Ma1和的情况，都有两个不同的、Ma2等值。原因是：对于一定的Ma1，气流经过正激波时，方向不变，即=0：而气流经过马赫波（无限微弱的压缩波）时，仍然=0。因此，当激波斜角由马赫角增大到90时，中间必存在某个最大折角max当激波斜角由开始逐渐增大时，相应地由0逐渐增到max；而继续增大到90时气流折角却相应地由max逐渐减小7.6 7.6 激波激波第99页/共122页(3)(3)激波图线及应用2）在同一Ma1之下，一个值对应着两个大者，代表较强的激波，称为强波；小者，代表较弱的激波，称为弱波。图中的虚线表示对应于max各点的联线，这条虚线把各图分成
50、两部分，一部分是强波，一部分是弱波。实际问题中出现的究竟是强波还是弱波，由产生激波的具体边界条件来确定。根据实验观察，方向决定的斜激波，永远是只出现弱波，不出现强波。7.6 7.6 激波激波第100页/共122页在超声气流中产生激波存在三种情况 1）由气流折转所确定的激波在超声速气流中，放置一块尖劈，尖劈的斜面把气流通道挤小了，气流受到压缩，发生激波。这是的激波是被斜面的角度所确定。7.6 7.6 激波激波第101页/共122页在超声气流中产生激波存在三种情况 2）由压强条件确定的激波在自由边界上由压强条件所确定的激波。例如超声速喷管出口的压强当低于外界大气压强时，气流将会产生激波来提高压强。