内容简介
《边界层转捩》注重于基础性和实用性研究的结合,总结了在边界层转捩方面的研究工作。主要内容分为三个部分:一是转捩过程,研究各种边界层稳定性问题及边界层转捩的预测;二是转捩机理,研究转捩流场的涡系结构,分析相关物理现象及其演化过程;三是转捩控制,以航空上的实际应用为背景,讨论转捩控制的方法和控制技术的新发展。
目录
目录
前言
第1章转捩引论1
1.1转捩现象1
1.2转捩类型1
1.3稳定性理论3
1.4转捩机理8
1.5转捩控制8
参考文献10
第2章线性稳定性理论10
2.1引言13
2.2Orr-Sommerfeld方程14
2.3特征值问题及其数值解16
2.3.1OSE的全局特征值差分解17
2.3.2局部特征值迭代法19
2.3.3二维扰动问题23
2.4非平行流边界层稳定性26
2.4.1非平行流边界层稳定性方程组26
2.4.2多重尺度法数值解27
2.4.3非平行性影响分析29
2.5线性稳定性理论与转捩预测的eN方法31
参考文献34
第2章抛物化稳定性方程方法36
3.1引言36
3.2稳定性方程的抛物化处理36
3.3扰动增长率表达式38
3.4正规化条件39
3.5空间推进数值解法40
3.5.1差分方法40
3.5.2初始条件和边界条件42
3.5.3线性PSE空间推进算法44
3.5.4二维PSE数值解45
3.6稳定性分析49
3.7PSE方法的优势51
3.8PSE的弱椭圆特性分析54
3.8.1弱椭圆特性55
3.8.2消除方法56
3.8.3解的一致性57
参考文献58
第4章非平行流非线性边界层稳定性61
4.1引言61
4.2二维非线性边界层稳定性62
4.2.1二维非线性抛物化稳定性方程63
4.2.2初始条件与局部法64
4.2.3均流变形与模态分析66
4.2.4算法与算例分析69
4.3空间演化的二次稳定性73
4.4三维非线性边界层稳定性77
4.4.1三维非线性抛物化稳定性方程77
4.4.2三维模态分析80
4.4.3数值方法81
4.5型三维扰动的非线性稳定性82
4.6型稳定性分析85
4.7附录A 非线性项的完整表达式87
参考文献94
第5章高速边界层稳定性97
5.1引言97
5.2可压缩流边界层方程及基本流参数98
5.3可压缩流边界层稳定性101
5.3.1扰动方程101
5.3.2线性抛物化稳定性方程102
5.3.3数值方法104
5.3.4非平行流稳定性分析105
5.4超/高超声速流的线性边界层稳定性108
5.4.1无黏稳定性方程及数值解108
5.4.2无黏稳定性与多重模态111
5.4.3黏性多重不稳定模态113
5.5超/高超声速流的非线性边界层稳定性117
5.5.1非线性抛物化稳定性方程117
5.5.2NPSE数值解118
5.5.3非线性稳定性分析119
附录B 方程05,10)与式(5.47)展开式124
参考文献131
第6章三维气动体可压缩流边界层稳定性133
6.1引言133
6.2三维可压缩流边界层135
6.2.1边界层方程135
6.2.2横向分段推进法137
6.2.3基本流速度分布138
6.3曲线坐标系的三维抛物化稳定性方程140
6.3.1Navier-stokes方程140
6.3.2扰动方程142
6.3.3抛物化稳定性方程143
6.4数值方法144
6.5三维边界层稳定性分析145
6.6三维边界层流动的转捩预测150
6.6.1转捩位置的预测问题150
6.6.2转捩预测的方法150
6.6.3转捩预测的DNS和PSE方法152
附录C 方程(6.19)展开式153
参考文献159
第7章边界层感受性问题162
7.1引言162
7.2感受性理论163
7.2.1自然感受性与强迫感受性163
7.2.2前缘感受性理论164
7.2.3局部感受性理论165
7.3渐近分析法165
7.3.1前缘感受性分析166
7.3.2渐近Orr-Sommerfeld分析167
7.3.3结合PSE方法167
7.4感受性问题的数值模拟170
7.5进展与分析174
参考文献175
第8章边界层转捩的直接数值模拟方法178
8.1引言178
8.2守恒形式N-S方程179
8.3初始条件179
8.4方程离散181
8.4.1紧致差分格式181
8.4.2空间离散185
8.4.3滤波函数185
8.4.4时间推进186
8.4.5离散方法比较187
8.5边界条件188
8.5.1壁面边界条件'188
8.5.2入口边界条件189
8.5.3展向边界条件189
8.5.4特征无反射边界条件189
8.6计算网格与并行算法193
参考文献197
第9章转捩边界层的典型涡结构200
9.1引言200
9.2主流向涡和次生流向涡202
9.3多种涡结构的形成207
9.4环状涡结构分析211
9.5 K型和H型转捩的后阶段涡结构214
参考文献216
第10章转捩过程的物理现象219
10.1引言219
10.2上喷和下扫220
10.2.1涡与一次上喷和下扫220
10.2.2环状涡与二次上喷和下扫222
10.3负尖峰和正尖峰224
10.4高低速条带结构227
10.4.1条带形成和演化227
10.4.2条带特性分析229
10.5高剪切层230
参考文献235
第11章边界层转捩的后期流场238
11.1引言238
11.2 U形涡和桶形涡239
11.3湍流斑的演化242
11.4小涡结构与无序化过程245
11.5转捩后期湍流形成的新认识247
11.5.1湍流形成的经典理论247
11.5.2刘超群对湍流形成的新认识247
参考文献251
第12章边界层转捩控制255
12.1引言255
12.2影响转捩因素256
12.2.1压力梯度与物体外形256
12.2.2表面粗糙度与壁面温度258
12.2.3湍流度与压缩性260
12.3转捩控制的基本途径262
12.3.1改变边界层基本流262
12.3.2改变扰动波265
12.4转捩控制技术的应用267
12.4.1转捩预测问题267
12.4.2 NLF机翼268
12.4. 3 LFG和HLFC机翼269
12.4.4 LFG机身272
12.5转捩、分离与激波272
12.5.1分离流转捩272
12.5.2转捩与激波边界层干扰254
12.5.3有层流分离的转捩点预计253
12.6转捩控制技术的新发展277
参考文献279
彩图
精彩书摘
第1章 转捩引论
1.1 转捩现象
1883年,雷诺第一次在实验中观察到在圆管水流动中存在层流和湍流两种不同的流态!当组合参数后来定义为雷诺数为平均流速度为管直径,为运动黏性系数,大于一定值时,流动会从层流转变为湍流。流体运动中的这种流态转变的物理现象称为转捩现象,广泛存在于许多流动问题中。转捩的发生与层流的不稳定性有关!是由于层流失稳导致了流动的转捩。从层流到湍流的极为复杂的转捩过程及研究工作的特别困难,使得这个重大的基础科学问题及其在工程技术上的应用,长期以来一直吸引着人们的高度关注。
层流向湍流的转捩是两种完全不同的流动形态的转变,图1.1展示了通过烟流法实验所看到的层流经过格栅后转变为湍流的过程"对应的格栅雷诺数。由图可以清楚地看到两种流态的显著差别,与规则的层流不同,湍流的流态十分紊乱,故又称为“紊流”(日文又名“ 乱流”)。古人曰:“堆出于岸,流必湍之”(三国+李康“ 运命论”),说明了湍流是自然界中一种客观存在的流动现象。
图1.1层流经过格栅后向湍流的转变(烟流法显示)
1940年,PRANDT提出了边界层(boundary layer,又称“附面层”)概念以后,人们开始研究边界层的流动特性,分析不同流动形态及其转捩现象'所谓边界层就是流动在物面边界附近所形成的一个黏性薄层(图1.2),开始于物面的前缘(或者是绕物体流动的前驻点),层内与层外的流动是一种渐近关系而没有确定的分界线,一般是在层内速度与外部速度相差1%的地方,作为边界层的外边界。边界层的厚度随流动不断增长,平板层流边界层的厚度可表为,其中,为来流速度,为雷诺数'由于空气,水等流体的运动黏性系数,很小,所以在一般流动中边界层的厚度是很小的(当雷诺数时,厚度)。在边界层内的流动黏性作用十分显著,黏性力项是与惯性力项同量级的,这是因为在很薄的边界层内,速度从物面的过渡到外边界的,导致速度梯度及其黏性力很大,特别是在物面附近'在边界层之外,黏性影响可以忽略而看成是理想流体流动。
图1.2平板边界层流动示意图
在通常研究的剪切流动中,包括边界层流,平面Poiseuille流及各种自由剪切流等,它们的流动转捩有许多共同特性,本书着重讨论相对较复杂的边界层的流动转捩,也可以供其他流动的研究作参考。
对于一般的转捩边界层来说,通常有三种不同的区域:层流区、转捩过渡区及湍流区。图1.3是采用阴影法显示的实验结果(超声速来流马赫数为),清晰地展示了三维旋成体边界层流动的三个区域。如图所示,随着流动向下游的发展,边界层厚度在不断增加,边界层内的流动也由规则逐步变为紊乱,尤其是在湍流区。通过这个边界层侧视图可以看到,与物体长度相比,边界层厚度的确是很薄的-若与图中的前缘激波很薄的厚度相比,它同样也是非常小的,这就进一步佐证了边界层厚度是很薄的真实情况。
图1.3三维边界层从层流到湍流的过渡(阴影法显示)
边界层转捩现象的另一个重要标识,就是在边界层转捩区中,流向速度沿法向分布(常称为速度型,或者速度剖面)的显著改变。图1.4是在平板边界层转捩区中测得的不同流向位置的速度型(这里的横向和纵向坐标分别为无量纲速度以及离壁面的高度#。在图中,“1层流”是Blasius层流边界层速度型,“ 2湍流”是湍流边界层(1/7幂次律)速度型,层流边界层对应于较小雷诺数(位置靠前),湍流边界层对应于较大雷诺数(位置靠后)。边界层的速度型曲线显示了两种流态的不同特性,尤其是在壁面处速度的法向梯度,湍流要比层流大得多。由于壁面摩擦阻力直接与该梯度成正比,这就意味着湍流的摩擦阻力更大。研究表明,在转捩区不同位置处速度型的变化,是与流动的稳定性特性及其转捩过程密切相关的。
图1.4平板边界层转捩区的速度型
1.2 转捩类型
从层流到湍流的边界层转捩,通常是由扰动引起的,是扰动随时间和空间演化的结果。转捩的类型主要与初始扰动有关,其过程也有不同。一般来说,对于初始扰动较小的边界层流动,不断增长的扰动波经历了线性和非线性阶段的发展和演化,从层流到湍流的转捩过程(图1.5)大致如下。
开始是外界的扰动进入边界层,产生不稳定扰动波。然后是扰动波在向下游传播过程中不断增长,起始阶段的扰动振幅很小,且各自独立演化,常称为“线性稳定性阶段”,可以忽略扰动的高阶小量,用线性稳定性理论来描述。这样的线性近似使得该理论仅能用于转捩的前期阶段。随着扰动的进一步增长,当扰动振幅达到一定值时(如自由流速度的1%的量级),需要考虑扰动波之间的相互影响,非线性作用已经不能再忽略不计了,此时流动进入“非线性阶段”。初期的非线性阶段常称为弱非线性阶段,在这个阶段中,各种不同频率和波数的特征扰动模态的相互作用(尤其是那些共振模态#已经较强,并能迅速放大,形成广泛频率谱的三维不稳定波。经过弱非线性阶段后的扰动继续增长,进入了边界层转捩的后阶段(又称为强非线性阶段,或实质非线性阶段),主要特征是各种复杂涡结构的形成和发展,并逐步向湍流过渡。
研究表明,环境扰动和基本流的特性能够影响转捩的过程,根据目前的理论,转捩的类型大致分为:
自然转捩(natural transition,或称正常转捩,常规转捩#,一般出现在背景扰动较低的情况。通过初始小扰动在层流边界层中激发形成Tollmien-Schlich-Ting波(T-S波),而后经过线性放大和非线性演化,形成如图1.6(a)所示的三维扰动波,不同涡系结构,强剪切层和湍流斑等,最终演变成湍流。转捩的扰动引入形式可以是自然扰动,也可以是人工扰动。在实验研究中,常用人工小扰动形式,如用振动带,声波或狭缝射流等产生扰动波。绝大多数的转捩属于自然转捩,也是书中关注的重点和主要研究对象。
旁路转捩(bypass transition,或称强迫转捩,“逾捩”型转捩。klebanoff等在边界层转捩实验中,观察到只要初始扰动足够大,自然转捩中的线性阶段被跳过,边界层扰动出现突变式的增长,而没有出现特征模态的增长方式,这样的转捩称为旁路转捩。这种转捩的初始阶段与自然转捩的不稳定性增长有显著的区别,一般是背景湍流度(如自由流扰动、糙表面等)情况下的平板边界层的转捩过程。显然,高湍流度与通常的低湍流度情况有很大差别,它跳过了线性阶段,经过(I)自由流局部涡扰动形成条纹、(II)初期湍流斑的出现及(III)湍流斑的聚合等阶段,最后完成层流向湍流的转捩[图1.6(b)]
图1.6不同自由流湍流度的平板边界层转捩过程比较
此外,在转捩中存在的斜波转捩(oblique wave transition)现象,也有将其单独归结为另一类转捩。Schmid等通过直接数值模拟方法首次研究在槽道流动中由一对斜波模态作为初始激励的转捩过程,在开始时没有加入、在过程中也没有观测到二维ABC波扰动。Berlin等进行了边界层的斜波转捩研究,展示了斜波转捩中的一些涡系结构的形成和演化。Wu等详细研究了一对斜波从线性到非线性的各个阶段,并分析引起转捩的演化过程。
需要指出的是,不同于从光滑壁面观察到的结果,若表面十分粗糙,也可能是直接转捩,而不出现湍流斑;若流动遇到较强的逆压梯度,或是存在某种分离的情况下,也可能会缺少自然转捩过程中的某些阶段。
1.3稳定性理论
一个多世纪以来,流动稳定性的研究经历了很长的发展过程,转捩开始于不稳定扰动波幅值的增长,从线性稳定性(包括平行流与非平行流)到非线性稳定性,形成了比较系统的理论。
1.线性平行流稳定性
从层流向湍流转捩的发生是由基本流场中早期的不稳定性引起的,与流动中很小的、有时也是不确定的扰动有关,这就产生了稳定性理论的基本思想:从一种了波扰动的放大和最终的层流溃变。历史上关于稳定性和转捩研究的成功和不成功的理论模型是很多的,现在用来研究稳定性的大多数模型是基于Prandtl假设:转捩是由小扰动放大引发的。假如流动使小扰动逐渐减弱并最终消失,流动恢复到原状态,那么该流动就是稳定的;反之如果该扰动逐渐增长,并不能恢复到原状态,则流动是不稳定的。因此,“稳定性”被定义为对抗小扰动的性质。最初,orr和sommerfeld分别用平行平板间的流动(流线相互平行#模拟二维波的放大过程,并假设波的振幅很小而忽略非线性项,控制方程简化为常微分方程,称为Orr-Sommerfeld方程(OSE)。Tollmien和schlichting在边界层稳定性计算方面做了开创性工作,先后通过求解OSE及其在稳定性分析中取得了突破,计算了平行流(忽略边界层缓慢增长)的中性稳定性曲线,以及包括中性曲线点之间的增长率。因此,上述分析中的二维波就被称为T-S波。但是仍然有很多人对稳定性理论持怀疑态度,直到Schubauer和Skramstad进行的风洞实验,揭示了转捩中不稳定波的测定规律,这是第一次对TOLLMIEN的理论预言进行的实验验证,得到基本一致的结果(鉴于他们所做的贡献,Morkovin建议T-S符号代表Tollmien、Schlichting、Schubauer和Skramstad)。
线性稳定性理论从不可压缩流扩展到可压缩流问题,在20世纪四五十年代就已取得了成功,随着计算机的广泛应用,能够得到可压缩线性稳定性方程的精确解。其中,Mack作出很大贡献,他探讨了线性稳定性的许多未知领域,特别是高马赫数时的各种模态,所发现的高阶模态,在不可压缩流中是没有相对应的,现在称其为mack模态。
2.线性非平行流稳定性
Gaster通过直接求解OSE以获得空间增长率,并给出空间与时间增长率的关系。然而,实验与理论的结果在较低雷诺数等情况并不吻合,使得人们开始怀疑OSE理论的平行性近似假设的局限性。将边界层流动当做平行流处理是基于这样的前提之下,即认为边界层的流向尺度远大于T-S波的波长。实际上,这对于靠近前缘或者边界层厚度变化激烈的区域来说是有疑问的,此时认为在边界层中所有位置上都具有相同增长率的平行流结果,显然是不准确的。因此,非平行性对稳定性的影响是需要考虑的。Bouthier采用多重尺度(摄动)技术,在分析线性边界层稳定性时,考虑了主流的非平行性,随后又有许多的研究和探讨,得到了非平行作用的可靠修正结果。
3.非线性稳定性
由于线性稳定性理论还远不能解决转捩问题,所以需要进一步研究非线性稳
前言/序言
边界层转捩 下载 mobi epub pdf txt 电子书 格式