空气动力学基础

我把Introduction to flight的第四章Basic aerodynamics略读了一遍，提炼了其中的重点要点，将其总结在一起分享给同学们，希望对大家空气动力学的学习有所帮助。

这个文档内容涉及的气流都是无黏的（书134—228页），没有包含黏性研究的部分。

因为领域导论书对黏性没怎么研究，基本都是只给结论，所以就不总结了。

本文档包括两部分，一是一些基本方程，二是这些方程的一些应用。

我读书只是蜻蜓点水，对一些公式的理解可能有错误；写的只是大致的推导过程，难免有不细致严谨之处；对一些英文的翻译可能不标准，同时可能输入有误。

希望大家批评指正、私下交流。

真心希望我们共同为之润色添彩，使其更加准确无误。

同时，大家有什么学习资料都记得共享啊，让我们共同进步！大家可以再看看领域导论书，看了这个总结，再看书就比较简单了。

看书最好也看看例题，例题不仅是对公式的简单应用，而且有些还包含新的知识，能增进我们对公式的理解。

这些内容只能算是一些变来变去的简单代数问题，大家不要有压力。

不过有几条注意事项：1、注意公式的限定条件，避免错误地加以应用。

2、大物书上的理想气体方程是Pv=RT，其中的R是普适气体常量（universal gas constant），领域导论书上的P=ρRT是经过变换的等价形式，其中的R是个别气体常量(specific gas constant)，等于普适气体常量R普适/M，大家变一下马上就懂了。

2、谈谈我的一个理解：本书中的研究好像不太强调质量和体积，可能是因为空气动力学研究没必要也不方便强调。

在一、基本方程——7、能量方程的推导中，v=1/ρ，这里的1应理解为单位质量，后面的能量方程中的V2也包含单位质量1，不然与h的量纲就不统一了；在二、公式应用——3、空速测定——C、高速亚声速流中，我们可以看出在本书中，Pv=RT，同样把大物书上的状态方程Pv=R普适T中的m当成单位质量1，并利用普适气体常量和个别气体常量的关系R个别=R普适/M，即可推出Pv=RT。

3、本书中涉及到比热（specific heat），用c v（对于等体过程）和c p（对于等压过程）在表示。

我们在大物中也学有c v和c p，不过它们不一样，不要混淆。

大物中那两个是摩尔热容（molar heatcapacity），分别为定体摩尔热容（molar heatcapacity at constant volume）和定压摩尔热容（molar heat capacity at constant pressure）。

对比起来有（下式中R个指个别气体常量，R普指普适气体常量，i指分子自由度，γ指热容比）：比热摩尔热容c v=R个，c p=R个c v=R普，c p=R普c p- c v= R个c p- c v= R普γ==γ==4、小写v代表体积，大写V代表速度，注意区分，其他字母符号的意义大家应该都能弄懂。

一、基本方程1、连续方程dm1=ρ1 dv1=ρ1A1V1 dt =ρ2A2V2 dt=dm2 则ρ1A1V1 =ρ2A2V2即ρAV=const对于不可压缩流，ρ1=ρ2，则A1V1= A2V22、欧拉方程（忽略了黏性和重力）在一个边长分别为dx dy dz的长方体流体元的x方向进行研究，忽略重力和黏性，朝向x正方向的力为P dy dz压强的变化率为则朝向x负方向的力为(P+dx) dy dz则合力F=P dy dz - (P+dx) dy dz=-(dx dy dz) 又m=ρdv=ρ(dx dy dz)，a===V 由F=ma化简得dP=-ρV dV3、伯努利方程（忽略了黏性和重力，适用于不可压缩流）对于不可压缩流，ρ不变，对欧拉方程进行积分，易得P1+ρV12= P2+ρV22即P+ρV2在一条流线上是常量，其中ρV2就是传说中的动压，用q表示，对于不可压缩流，P+ρV2等于总压，我们在方程的应用中会再提及。

4、关于热力学第一定律系统的内能增量=外界传热+外界做功，即de=δq +δw其中δw=-P dv(压缩，所以v减小，dv是负值，所以有负号)则δq=de+P dv定义焓h=e+Pv做微分得dh=de+v dP + P dv与上式一起消去de得δq=dh- v dP5、内能与焓定义比热（specific heat）c=，即系统增加单位温度所吸收的热量等体过程的比热写作c v，等压过程的比热写作c p 对于等体过程dv=0代入δq=de+P dv可得de=δq=c v dT从e=0和T=0积分得e= c v T我们在大物中学的是e=R普T，m还是要当做单位质量1，推出e=R个T=c v T。

因此，它们是等价的。

对于等压过程dP=0代入δq=dh- v dP则dh=δq=c p dT从h=0和T=0积分得h= c p Tde=c v dT，e= c v T，dh=c p dT，h= c p T四式虽然是从等体过程和等压过程推出的，但对于理想气体是普遍适用的。

6、等熵过程（适用于等熵过程）对于等熵流（绝热可逆）δq=0代入δq=de+P dv和δq=dh- v dP则-P dv=de = c v dT，v dP=dh= c p dT两式相除得=-γ其中定义了热容比γ=c p/c v对于空气，γ=，应该是因为空气的绝大部分是氮气和氧气，都是双原子分子，分子自由度i=5，根据大物中学的热容比γ=，可得γ=。

再积分=-γ得=()-γ把体积换成密度得=()γ同时借助状态方程ρ=P/(RT)在有ρ的那个式子中消去ρ或借助我们熟悉的形式（大物书上的）Pv=RT在有v的那个式子中消去v可得=（）γ/(γ-1)总结：=()γ=（）γ/(γ-1)，即Pρ-γ=常量，Pγ-1T-γ=常量把大物书上的式子中的体积换为密度，就跟这个完全一样了7、能量方程（适用于无黏）对于绝热过程δq=dh- v dP=0代入欧拉方程dP=-ρV dV得dh+vρV dV=0v=1/ρ（这里的v应理解为单位质量的体积）则dh+V dV=0做积分得h1+V21= h2+V22，即h+V2=常量代入h= c p T得c p T1+V21= c p T2+V22，即c p T+V2=常量对于非绝热过程δq≠0可得δq=dh+ V dV做积分=+得h1+V21+Q12= h2+V22也可写为c p T1+V21+ Q12= c p T2+V228、一个重要结论对于等熵流，总温T0，总压P0，总密度ρ0是定值总温（total temperature），总压（total pressure），总密度（total density）定义：Total temperature/pressure/density at a given point ina flow is the temperature/pressure/density thatwould exist if the flow were slowed downisentropically（等熵地）to zero velocity二、公式应用1、声速公式的推导由于声波穿过气流与气流以声速穿过声波等价，因此可用后者来研究在声波两侧，设气流压强分别为P和P+dP密度ρ和ρ+dρ温度T和T+dT速度a和a+da应用连续方程有ρA1a=（ρ+dρ）A2（a+da）A1=A2，则ρa=（ρ+dρ）（a+da）展开，再忽略无穷小量dρda,可得a=-ρ代入欧拉方程dP=-ρa da，即da=-可得a2=通过声波的气流是等熵流，则P/ργ=const=c因此=cργ=cγργ-1代入c= P/ργ,得a2=γ对于理想气体，可以再代入状态方程P=ρRT最终得出a=可以看出，理想气体中的声速仅与温度有关2、低速亚声速风洞设Settling chamber（reservoir）和Test section的气流速度分别为V1，V2压强分别为P1，P2面积分别为A1，A2通过低速亚声速风洞的气流可以看作不可压缩流，由连续方程和伯努利方程可得V1=V2，P1+ρV12= P2+ρV22联立两式消去V1，可得V2=A2/A1对于给定风洞是定值，要想调节Test section的速度大小，可以调节P1-P2。

以前人们用U型管分别连接Settling chamber（reservoir）和Test section来测P1-P2，现在我们工艺先进，通过压力传感器实现3、空速测定A、设备：总压管（Pitot tube），空速管（Pitot-static tube）B、对于低速亚声速流(M<在上图中空速管上的A点压强为静压P，速度为V1在B点压强为总压P0，速度为0应用伯努利方程得P0=P+ρV12可得V1=定义动压q=ρV2（此定义式对所有气流都成立）可得P0= P+q（注意此式和P0= P+ρV12只对不可压缩流成立）可见：只要设法获得P0- P和ρ的值，就能求出速度，P0- P的测定通过空速管或总压管可以实现。

对于ρ，若使用真值（true value，即设法测的飞机周围的ρ），则获得真实空速（true airspeed）V true=但是测定飞机周围的ρ比较难，所以低速飞机计算时都是用的标准海平面密度ρs，获得当量空速（equivalent airspeed）V e=其实当量空速有更深层次的意义：Consider an airplaneflying at some true airspeed at some altitude .Itsequivalent airspeed at this condition is defined as the velocity at which it would have to fly at standard sea level to experience the same dynamicpressure.给定了当量空速，就相当于给定了动压。

当量空速的概念十分重要，在研究飞行表现时很有用。

C、对于高速亚声速流（<M<1）由h=e+Pv=e+RT,即c p T=c v T+RT可得c p-c v=R根据γ= c p/ c v可得c p=气流在空速管或总压管前的探针前的停滞点（stagnation point）处从最开始温度T1和速度V1的状态等熵静止,速度变为0，因此温度为总温（total temperature）T0，压强为总压（total pressure）P0，在此过程应用能量方程得c p T1+V12= c p T0变换等式可得=1+代入c p=可得=1+又声速a12=γRT1则=1+=1+M12(只要求是绝热过程)结合等熵过程方程=()γ=（）γ/(γ-1)可得=（1+M12）γ/(γ-1)，=（1+M12）1/(γ-1)（要求是等熵过程）由上式可得M12=[()(γ-1)/γ-1]因此，通过总压和静压的比值可以直接求出马赫数代入M1=V1/a1则V12=[()(γ-1)/γ-1]也可写为V12=[(+1)(γ-1)/γ-1]因为实践中一般获得P0-P1,所以上式用得较多而且，由于T难以测量，即a难以获得，静压P1也难测，所以高速亚声速飞机一般用标准海平面的声速和压强a s和P s代入上式。