空气动力学基础（笔记自留）-第零至第一章

第零章 绪论

0.1 空气动力学的研究对象、范围和分类

1. 研究对象： 空气与物体之间有相对作用时，空气的运动规律及作用力（包括空气内部的和空气与物体之间的相互作用）所服从规律
2. 范围：地面高速交通关键技术中的减阻、稳定和降噪问题；涡轮机、轴流式或离心式压气机等叶片机中的气动力问题，鼓风机和高炉中的气流问题，大建筑物中的暖气通风问题，高大建筑物的风压问题，自然界的气象问题。
3. 分类：
   1. 飞行马赫数：低速空气动力学（Ma<0.3）和高速空气动力学，高速空气动力学又可分为亚声速（0.3≤Ma<0.8）、跨声速（0.8≤Ma≤1.2）、超声速（1.2<Ma<5）和高超声速（Ma≥5）。
4. 研究问题：外部流动（空气对物体的绕流）+内部流动（气流在通道中的流动问题）

0.2 历史概述

1. 公元前到19世纪的流体力学发展
2. 20世纪：空气动力学建立完整体系

0.3 空气动力学的研究方法

1. 理论分析方法
2. 试验方法
3. 数值方法

第一章 流体的物理属性及流体静力学

1.1 流体的连续介质模型

1.1.1 连续介质假设

流体充满一个体积时不留任何自由间隙，流体的特性连续变化

1.1.2分子平均自由程

一个气体分子在随机运动过程中两次碰撞之间所走过的距离的平均值。
当分子平均自由程与飞行器的某个特征强度远小于1时，分子间距可以忽略不计，即可以认为空气介质是无间隙的。

1.1.3 流体的宏观特性及流体物理量

流体的微观特性：流体的物理量随时间和空间的变化都是不连续的。
流体的宏观特性：物理量的统计平均值在时间上是确定的，在空间上是连续的。
流体质点：微小体积中的所有流体分子的总体。
连续介质模型：流体是由连续分布的流体质点组成的，流体流动空间中的每一点都被相应的流体质点占据。空间任意点上的流体物理量就是指位于该点的流体质点的物理量。

1.1.4 连续介质假设的适用范围

80Km以上的高空或问题的特征尺寸L变小时，连续介质假设不再成立。
无量纲参数：努森数Kn=λ/L.
Kn远小于1时，可认为连续介质假设成立。
工程上规定Kn<0.01时，可采用连续介质假设。
根据Kn数范围划分的流体状态区域并非严格和绝对的。

1.1.5 流体质点与流体微团

流体质点是大量分子的集合，是可以忽略线性尺度效应的最小单元。
流体微团是由大量流体质点所组成的，具有线性尺度效应的微小的流体团，尺度为微分dx的量级。

1.2 流体的流动性、黏性与理想无黏流体

1.2.1 流体的流动性

流体的易变形性：流体在很小的剪切力作用下将发生连续不断的变形。

1.2.2 流体黏性的表现

流体一旦运动，流体内部就具有抵抗剪切变形的特性，以内摩擦力的形式抵抗流程之间的相对运动。
流体在固体表面的无滑移条件：分子间的内聚力将流体黏附在固体表面上，随固体一起运动或静止。
流体层之间的摩擦力是黏性摩擦力。因为流体有黏性，同时相邻流体层之间有相对运动，或者说具有速度梯度时发生的。

1.2.3 牛顿内摩擦定律

“流体的两部分由于缺乏润滑而引起的阻力，同流体两部分彼此分开的速度成正比”。
黏性切应力：单位面积上的黏性摩擦力，记为τ。
黏性切应力与速度梯度间的关系可由平行平板间的黏性流体流动试验得到。
平行平板流动具有以下特点：

1. 与上板接触的流体黏在板上，以速度U随上板运动；与固定板接触的流体也黏在板上，速度为零；两板之间的流体流动速度呈线性分布。
2. 对于某些流体（如水和空气），平板和流体受到的切应力与速度成正比，与两板之间距离成反比。
3. 进一步试验表明，这个结果推广到任一层流直线运动中去：牛顿黏性切应力公式，给出了切应力与速度梯度的线性关系，满足这种关系的流体称为牛顿流体。

1.2.4 黏度（粘性系数）

黏度/动力黏度：黏性应力公式中的比例系数μ，国际单位为N·s/m^2，反映了流体抵抗剪切应变率的能力。

1. 黏性系数随温度变化特点及原因
   不同介质的黏度μ不同，液体的黏度比气体的大得多。同一种介质的μ值则随温度变化。
   黏度随温度的变化规律是由流体黏性的微观机制决定的。
   液体的黏性起因于分子内聚力，当温度升高时，液体分子运动幅度增大，使分子平均距离增大，而分子间吸引力是随时间增大而减小的，因此内聚力减小，黏度相应减小。
   气体的黏性起因于分子热运动，当温度升高时，分子热运动加剧，气体通过热运动进行宏观动量交换的能力增强，黏性相应增大。
2. 黏度计算公式
   苏士南公式/幂次公式
   运动黏度（运动黏性系数）v=μ/ρ

1.2.5 理想无黏性流体

一般流体的黏度都很小，在流动的速度梯度不太大时，黏性力也就不大。如果黏性力比作用在流体上的其他力小得多，可以认为黏性力为零。

1.2.6 气体的导热性

导热性：气体中沿某个方向n存在温度梯度时，热量会由温度高的区域传向温度低的区域。
单位时间内通过垂直于n方向的单位面积所传递的热量q按傅里叶导热定律确定。
同黏度一样，气体的热传导系数随温度的增高而增大，同时气体热传导系数的数值也非常小。

1.3 完全气体的热力学特性

1.3.1 完全气体状态方程

气体流动时，各点的热力学参数可以是不同的，处于热力学平衡态的气体，其压强、密度和温度之间存在着某个函数关系。
气体的状态方程：p=p(ρ，T)。
对于远离液态的气体，也就是温度不是特别低 、压力不是特别高时，可以采用完全气体的模型给出状态方程。
完全气体状态方程：p=ρRT。R=普适气体常数/气体摩尔质量。
极高压或低温是，完全气体假设不成立。
称未发生化学反应的完全气体为热完全气体。

1.3.2 气体的内能

1. 定义：气体分子热运动的能量、分子间相互作用的能量和分子内部的能量的总和。
   对于热完全气体，假定只有分子热运动，不计分子间相互作用，内能为分子热运动的能量，包括分子的平动能、转动能、振动能和电子激发能。
2. 结果：
   单位质量的双原子分子气体
   1 . 平动能：etr=3RT/2
   2 . 转动能：erot=RT
   3 . 振动能：evib=Rθv/(e^(θv/t)-1),θv为振动特征温度
   4 . 总内能：e双原子= 3Rt/2+RT+Rθv/(e^(θv/t)-1)
   单原子分子气体
   1 .内能=平动能 ，e单原子=3RT/2
   热完全气体的内能e仅为温度T的函数。
   温度不高时，视振动能未激发，为零。此时气体称为量热完全气体，双原子分子的内能为e双原子=5RT/2，单原子分子的内能仍为3RT/2。
   对于空气，在800k以下视为量热完全气体。

1.3.3 热力学第一定律、焓和比热容

1. 热力学第一定律
   能量守恒定律在热能和机械能相互转换的情况中的应用。
   微分形式：δq+δw=de
   对于准静态无摩擦这样的可逆过程，外界对系统做工为压力做工：
   δw=-pd(1/ρ)
   δq=de+p(1/ρ）
2. 气体的焓
   h=e+p/ρ
   代入完全气体状态方程式，得到：h=e+RT
   对于800k以下的空气，有e空气=5R空气T/2，h空气=7R空气T/2‘
   引入焓后，在只有压力做工的情况下热力学第一定律可以改写为：
   δq=dh-dp/ρ
3. 比热容与比热容比
   1 . 比热容是热力学中衡量物质吸收热量和其温度升高之间关系的量。
   气体的比热容与具体的热力学过程有关，吸热引起气体温度变化与吸热过程中的气体的压力和体积变化有关。
   比定容热容（定容比热）指保持气体体积（比容）不变的情况下，单位质量气体升高一摄氏度所需吸收的热量，表示为cv:
   cv=δq|v=const/dT
   比定压热容（定压比热）指保持气体压强不变的情况下，单位质量气体升高一摄氏度所需吸收的热量，表示为cp:
   cp=δq||p=const/dT
   比热容比（比热比）定义为γ=cp/cv
   2 . 内能和焓与比热容的关系
   由焓与内能的关系，cp和cv之间存在以下关系：
   cp-cv=R
   由gamma=cp/cv，得cp=γR/(γ-1)，cv=R/(γ-1)
   3 。 量热完全气体的比热容
   对于量热完全气体，内能和焓都与温度成正比，因此比定容热容和比定压热容都为常数：
   cv,双原子=5R/2，cv,单原子=3R/2
   cp,双原子=7R/2，cp,单原子=5R/2
   内能和焓与温度的关系可写为：
   e=cvT，h=cpT
   此时比热容比也为常数：
   γ双原子=cp/cv=1.4，γ单原子=5/3
   空气是双原子分子的混合物，88K以下可视为量热完全气体，比热容和比热容比分别为：
   cv,空气=5R空气/2，cp,空气=7R空气/2，γ空气=1.4

1.3.4 热力学第二定律和熵

决定过程向哪个方向进行。

1. 熵的定义与热力学第二定律
   经典热力学中定义了一个状态函数——熵：ds=δqrev/T
   上式中，s为单位质量系统的熵，δqrev为可逆过程中外界传给该系统的热量，T为热源的温度（可逆过程中，系统和热源温度相等）。
   若从状态1到状态2经历的是一个不可逆过程，熵的变化依旧是ds，满足ds=δq/T+dsirrev.
   上式中，δq/T称为熵流，是过程量，不是全微分。dsirrev是不可逆增熵，也是与过程有关的量，是不可逆过程中由于黏性摩擦耗散、不可逆热传导、质量扩散产生的熵增。
   热力学第二定律通常的表达式为：ds≥δq/T。对于绝热过程，有ds≥0.
   热力学第二定律指出，任意热力学过程中系统的熵增大于或等于系统吸收的热量除以热量温度，绝热过程中系统（孤立系统）的熵永不减小。
2. 熵的计算式
   根据熵的定义，可逆过程中熵增与系统吸收热量之间的关系为：
   δqrev=Tds
   分别代入热力学第一定律表达式，得到：
   Tds=de+pd(1/ρ)和Tds=dh-dp/ρ
   对于热完全气体，有de=cvdT和p=ρRt，代入上式，有：
   ds=cvdT/T-Rdρ/ρ
   对于某一热力学过程从初态1到终态2积分，对于量热完全气体，cv=R/(γ-1)为常数，
   s2-s1=cvInT2ρ1^(γ-1)/T1ρ2^(γ-1)
   代入dh=cpdT和p=ρRT，积分，对于量热完全气体，cp=γR/(γ-1)，得
   s2-s1=cvp2ρ1^γ/p1ρ2^γ
3. 熵含义的理解
   熵其实是系统无序性的度量。
   1 . 升高温度T，熵s增大。
   2 . 增大体积，熵s增大。
   3 . 升高压强p，熵s减小。
   熵的定义表明，可逆吸热过程会使系统熵增加。以热完全气体的两个热力学过程为例，说明具体过程。
   1 . 等温吸热过程中de=0，d(1/ρ)>0，密度减小，压强减小，熵增加。
   2 . 等容吸热过程d (1/ρ)=0，de>0，温度升高，熵增加。

1.3.5 等熵关系式

等熵过程中，满足以下公式：
1 . T/ρ^γ-1=const或T2/T1=（ρ2/ρ1）^γ-1
2 . p/ρ^γ=const或p2/p1=(ρ2/ρ1)^γ
3 . p2/p1=(ρ2/ρ1)^γ=（T2/T1）^γ/γ-1

1.3.6 声速

声波定义：对弹性介质世家一个任意的小扰动，就会在介质中引起微小的压力增量，并以波的形式向四周传播。小扰动波的传播速度为声速。
物理学中声速关系式：a=sqrt(Es/ρ)
上式中，Es为等熵过程的体积弹性模量。（通常认为声波的传播为等熵过程，来不及与环境进行热交换）
对于热完全气体，a=sqrt(γRT)，声速只与温度有关。
声速是热力学参数的导数，也是介质自身的属性，与宏观流速无关，与分子热运动状态有关。
分子运动论给出的分子平均热运动速率为：C平均=sqrt(8RT/π)
可见声速的平方也可以作为气体分子热运动内能的度量。

1.4 流体的可压缩性与不可压流动

1.4.1 流体的可压缩性与可压缩系数

流体的可压缩性指流体的体积/密度在压力作用下发生改变的性质。
可压缩系数β，表示在一定温度下升高一个单位压强时，液体体积的相对减小量
β=-dv/vdp
上式中，用负号后可压缩系数为正，v为比容，表示单位质量流体的体积，v=1/ρ。可压缩系数也可以写为β=dρ/ρdp。
可压缩系数常要指明压缩过程的特点。

1.4.2 流体的弹性和弹性模量

1. 定义：流体抵抗压缩变形的能力和特性称为弹性，可用体积弹性模量来表征。体积弹性模量定义为产生单位的体积相对变化（减小）所需的压强增量:E=1/β。
   流体的体积弹性模量就是可压缩系数的倒数。

1.4.3 不可压流体/流动与可压缩流体/流动

一般液体很难压缩，将其当做不可压流体。
在研究的具体问题中，压强变化很小，虽然气体的弹性模量较小，但密度变化也很小，将其作为不可压流体处理，认为其密度为常数。
低速流动又称为不可压流动。
高速流动又称为可压缩流动。

1.4.4 采用马赫数划分流动范围的物理内涵

马赫数为飞行速度/气流速度与声速之比：Ma=V/a
Ma<0.3时流动为低速（不可压）流动。
Ma>0.3时流动为高速（可压缩）流动。
声速的平方可以作为分子随机热运动内能的度量，马赫数的平方可以反映气体的宏观动能与内能的比较。
对于量热完全气体，单位质量气体的内能为：e=cvT=RT/γ-1
气流宏观动能与内能之比为v²/2e=γ(γ-1）Ma²/2

1.5 作用在流体微团上的力

作用在流体微团上的力按照作用力的性质和作用方式可以分为彻体力和表面力两类，彻体力是由外力场施加的，表面力是由与微团接触的其他流体或物体施加的。

1.5.1 彻体力

1. 定义：由外力场作用于物体微团的质量中心、大小与微团质量成正比的非接触力，也称体积力或质量力。包含重力、惯性力和电磁场运动的带电气体所受的电磁力等。
   分量形式：f=fxi+fyj+fzk，f在数值上等于加速度，重力场中，彻体力为f=-gj

1.5.2 表面力

1. 定义：由物体或相邻流体作用在流体微团的外表面上的、大小与微团表面积成正比的接触力称为表面力。
   单位面积上的1表面力称为应力，单位为Pa，N/m²。

1.5.3 流体内一点处的压强

1. 定义：对于无黏流体或静止流体，只有垂直于表面的法向力。
   单位：1atm=760torr=1.013x10⁵Pa，1torr=1mmHg=133.3Pa
2. 理想流体压强的方向无关性
   对于理想无黏流体，同一点的各个不同方向上的压强相等。

1.6 流体静力学

研究流体处于平衡状态时的力学规律，包括绝对平衡和相对平衡（研究流体处于平衡状态时，内部压力与彻体力的关系，由此确定流体内部压力的分布规律及流体对于周围固体壁面的作用力）

1.6.1 流体静力学的基本方程——欧拉静平衡方程

1. 推导：∂p/∂x=ρfx，∂p/∂y=ρfy，∂p/∂z=ρfz
   矢量形式：▽p=ρf
2. 静平衡方程其他形式
   1 . 静压力差公式：dp=ρ（fxdx+dydy+fzdz）
   2 . 彻体力有势，记势函数∏（x,y,z），有dp=-ρd∏
   密度不变时，p=-ρ∏+c
   3 . 非惯性系中的静平衡方程
   流体整体具有加速的的情况同样适用，建立以流体整体运动为基准的非惯性系即可。
   欧拉静平衡方程为：▽p=ρf=ρ（fg+fa）
   上式中，fg为流体真实受到的彻体力，惯性系中的彻体力；fa为单位质量流体在加速度为a的非惯性系中受到的惯性力，fa=-a

1.7 国际标准大气

略