0 引 言

高超声速扁平形升力体的气动力和气动热问题是当前航空航天领域遇到的重要问题，其短轴处的转捩位置和热流明显与两侧不同[1-2]。这是由于短轴处基本流存在流向涡[3-6]，与两侧的基本流有本质的区别，导致短轴处稳定性与两侧有本质区别，最终导致转捩位置和热流不同。为了预测短轴处的转捩和热流，可以椭圆锥作为模型，研究其短轴处的稳定性。该处基本流沿展向和法向的变化同量级，而沿流向的变化要小得多，所以是二维全局稳定性问题。

人们认识到短轴处是二维全局稳定性问题，经历了一个过程。20世纪90年代中期，人们开始注意到短轴处的不稳定性[7]，但一直用一维剪切层不稳定性来考虑问题[8-9]，所得计算结果始终无法得到实验数据[10]的支持。直至2009年，Choudhari等[11]才明确指出该处的稳定性与两侧完全不同，需要用全局稳定性分析研究。直到最近，Paredes和Theofilis[12-13]才在长短轴比为2∶1椭圆锥，来流马赫数7，单位雷诺数分别为8×105/m和1.89×106/m的短轴处找到了不稳定的外模态，其特征是特征函数峰值位置远离壁面，靠近边界层外缘，相速度较高。除了外模态，Zhang和Luo[14]在不可压展向局部型自由湍流带条纹边界层中还找到过内模态，其特征是特征函数峰值位置靠近壁面，相速度较低。

在Paredes和Theofilis[12-13]研究的基础之上，有以下几个更深入的具体问题需要回答：(1) 是否有不稳定的内模态；(2)有哪些不稳定的外模态；(3)同一流向位置，不稳定模态随波数(频率)的变化规律；(4)不稳定模态沿流向的演化规律；(5)不稳定模态的物理本质，是否都由无粘的拐点不稳定性造成的。本文正是针对这几个问题进行研究。

1 数值计算方法

本文使用Arnoldi方法和Rayleigh商迭代法计算椭圆锥短轴处二维全局稳定性问题，Arnoldi方法能够求解指定范围的全部特征值，能够很方便的用来寻找特征模态的种类[14]。Rayleigh商迭代法，对于追踪某一特定模态随着某一参量的演化更为方便。

1.1 基本流求解方法

椭圆锥基本流是通过直接数值模拟求解三维守恒型可压缩N-S方程得到的，控制方程如下：

(1)

第二阶段：通过标准修订及数据更新，将军有民无中的部分非军事特有要素纳入到第一部分中，扩充军民共同表达要素。

本文采用Sutherland公式计算粘性系数及热传导系数，即：

采取动态心电图分析仪予以检测，采用数字式的动态记录盒，应用mV5、mV1以及mV3三通道心电图对患者的24小时心电图情况予以详细记录，且配有起搏脉冲相关标记，采取计算机对记录予以回放，进行误差纠正后，分析和评估患者的起搏器起搏功能异常及起搏器间歇性感知功能字符异常情况。

(3)

其中，Cκ、Cμ为参考温度，本文均取为110.4 K。在实际情况中，因为椭圆锥模型较为复杂，使得计算网格具有不规则性和复杂性，在一般曲线坐标系计算更为方便。通过坐标变换，得到一般曲线坐标系下的控制方程：

改变总是困难的，“数据分析”可以帮助企业做出显著的改变：为了在全新的协作式、数字化、以“数据分析为主导”的世界生存，人们被迫改变思维方式和技能。

(4)

式中：

(5)

1.2 稳定性分析方法

二维全局稳定性分析方法(Bi-Global)，是一维线性稳定性理论的扩展，是在当地求解特征值问题。无量纲化后的非守恒型N-S方程如下：

·V)=0

[λ(·V)]+

式中，(x,y,z)表示空间坐标，t表示时间变量，U=(ρ,ρu,ρv,ρw,ρe)T为守恒变量，ρ为密度，e=cvT+(u2+v2+w2)/2，cv为比定容热容，T为温度；u、v、w分别表示x、y、z方向的速度分量，E、F、G分别为x、y、z方向的对流项通量，Ev、Fv、Gv分别为x、y、z方向的粘性项通量。

该计划根据社区业余俱乐部人数、硬件设施、训练及比赛数量、竞赛成绩等因素将社区足球俱乐部分为3个级别，在提供专业教练员及场地支持外，还会依据级别赞助不同数量及质量的服装与训练装备。

·(κT)

(6)

式中：

将瞬时物理量写成基本流和扰动量之和，消去基本流满足的方程，忽略二阶及以上扰动量，得到非守恒形式的线性扰动方程。通过链式求导法则，将直角坐标系的导数转化到一般曲线坐标系下，得到一般曲线坐标系下的线性扰动方程:

(7)

椭圆锥短轴附近的流向涡，在展向变化剧烈，不能像LST一样设成行进波形式，但流向变化依然较为缓慢，因此流向引入局部平行假设，设成行进波：

注意均衡营养，维生素与矿物质的补充，养成良好的生活习惯与规律的起居作息，戒烟戒酒，不熬夜；缓和情绪反应与调适工作压力，在充分认知之下，避免太多无谓的错觉妄想，保持平常心；改善工作环境，避开所有可能有危害的污染物质。

(8)

代入上式，即可以得到二维线性稳定性方程：

(9)

其中:

(10)

2 计算模型及程序验证

2.1 计算模型

计算采用的是HIFiRE-5飞行试验的38.1%缩比模型，模型为长短轴比2∶1的椭圆锥，短轴所在平面头半径1 mm，半锥角7°。考虑出口对基本流的影响，模型计算域长度延长到380 mm。模型的示意图如图1所示，基本流计算的来流条件，参考了Juliano实验[2]的工况,如表1所示。

  

图1 HIFiRE-5模型Fig.1 Sketch of HIFiRE-5 model

模型零攻角工况下流场具有对称性，为节省计算资源，只取1/4模型进行计算，网格数601×251×201，为了得到一个更好的流场，计算网格在模型头部、边界层内和激波处适当加密了网格。特征长度取1 mm，数值计算中，时间推进采用具有TVD性质的三阶Runge-Kutta格式，对流项采用AUSM+格式进行通量分裂。流场中含有激波，采用2阶NND格式捕捉激波。展向采用对称边界条件，极轴处存在奇点，采用极轴边界条件，出口采用线性外推边界条件。外边界采用远场边界条件，壁面采用无滑移等温条件。具体格式以及各控制方程具体系数矩阵可以参见张绍龙博士论文[15]。

 

表1 基本流计算的来流条件Table 1 Flow condition for the basic flow calculation

  

Ma∞Reunit/mT∞/KTw/KAOA/(°)6.010.2×10652.0300.00.0

HIFiRE-5模型由于在长轴子午面和短轴子午面之间存在压力梯度，导致边界层内流线的弯曲，基本流在短轴附近存在一个“蘑菇”状的流向涡，该涡结构沿流向变化缓慢，且缓慢抬升，沿展向和法向变化较为剧烈。图2给出了短轴附近流向涡的稳定性分析计算域选取范围示意图以及流向涡沿流向演化数值计算结果。

无论是哪种风格的影视动画作品，都包含着特定的故事内容和情节，后期制作要对其加以丰富和深化，从而加强作品的戏剧性效果。后期处理要尊重剧本的初始设计，避免对作品的精神内核造成破坏。为了引发观众产生共鸣，后期制作者要在遵循故事主线的前提下，对故事情节进行适当的剪辑和增加，并通过网络技术强化视觉效果，通过特效等使作品更加惊险刺激、撼动人心。这一过程中，后期制作不能打乱故事发展顺序，也不能另起炉灶，要确保故事合理，抓住情绪的起伏点，深化作品的内涵，优化动画影视作品的逻辑结构，使观众走进情境，从而增加认同感。

  

(a)

  

(b)

 

图2 (a)全局稳定性分析选取计算域和

 

(b)椭圆锥短轴流向涡Fig.2 (a)The computational domain of global stabilityanalysis; (b) Streamwise vortices near minor-axismeridian of hypersonic elliptic cone.

2.2 线性稳定性分析程序验证

为了验证二维全局稳定性分析的程序，与一维线性稳定性方法Muller法(LST)和直接数值模拟计算结果进行了对比。

2.2.1 与Muller法对比

二维全局稳定性分析与一维线性稳定性分析的最主要区别体现在考虑了基本流沿展向变化，因此扰动在展向不能设成波动形式。但当基本流展向均匀的时候，二维全局稳定性分析仍有效，且与一维线性稳定性计算结果应该是一致的。选取平板边界层，具体参数见表2，全局稳定性分析与用Muller法求解一维线性稳定性分析的特征值结果对比如表3所示，所得特征函数如图3所示，可见两种方法所得结果一致，程序是可靠的。需要指出，这里我们计算了时间模式和空间模式两种情况。

 

表2 验证Bi-Global程序所用基本流的工况Table 2 Basic flow freestream condition forverifying of Bi-Global program

  

Ma∞Tw/KT∞/KRe6.0300.079.01.0×105

  

(a) 时间模式

  

(b) 空间模式

 

图3 Bi-Global方法和Muller方法计算的平板边界层特征函数的比较Fig.3 Flat plate boundary layer eigenfunction, Bi-Global and Muller comparison

 

表3 Bi-Global方法和Muller法计算的平板边界层特征值的比较Table 3 Flat plate boundary layer eigenvalues comparisonbetween Bi-Global and Muller

  

αrαiωrωiMuller(temporal)1.601.46570.0466Bi-Global(temporal)1.601.46650.0462Muller(spatial)1.5917-0.06311.460Bi-Global(spatial)1.5907-0.06241.460

2.2.2 与直接数值模拟对比

为了验证基本流沿展向变化较大时,全局稳定性程序的计算是否可靠，与直接数值模拟做了对比。基于全局稳定性分析基本流，以X*=192mm处的流场剖面，沿着流向均匀延拓，构造了计算域为20倍波长的平行流。将全局稳定性分析结果作为直接数值模拟的入口条件，向下游进行计算。图4是用密度、速度、温度和压力不同参量得到的增长率，可见，经过入口附近的调整，增长率最终都与全局稳定性分析(空间模式)吻合。DNS所得特征值见表4，特征函数见图5，都与全局稳定性分析的结果吻合。

·[μ((V)+(V)T)]}

  

图4 全局稳定性与直接数值模拟得到的扰动增长率Fig.4 Perturbation growth rate obtained by globalstability and direct numerical simulation

 

表4 全局稳定性分析与直接数值模拟对比：特征值Table 4 Comparison between global stability analysisand direct numerical simulation: eigenvalues

  

ωαr-αiDNS0.30.50500.0222Bi-Global0.30.50510.0219

 

 

 

 

 

 

 

图5 全局稳定性分析与直接数值模拟对比： 特征函数(等值线为全局稳定性分析结果， 彩色云图为直接数值模拟结果)Fig.5 Comparison between global stability analysisand direct numerical simulation:eigenfunction

3 计算结果及讨论

3.1 是否有不稳定的内模态

Parades和Theofilis[12-13]在短轴处找到了不稳定的外模态。而在不可压二维全局稳定性分析中，既发现了外模态，也发现了有不稳定的内模态，且内模态的增长率较大[14,17]。因此，本文认为有必要在短轴处观察是否有内模态。

这一年图科回归，总经理叶鹏又选来了一个瘦高灵活，能里能外的外援拉米扎纳。这两人特点互补性格开朗，每场都能给对手带来强大的压力，人称“黄金双枪”。

在跨境电商旅游英语教学过程中教师应结合学生们的实际情况、学习需求及教学计划等通过实施微课视频的制作教学及翻转课堂的教学模式促进学生们在跨境电商旅游英语学习过程中自主探究及学习，在自主探究与学习理论与实践技能的过程中提升学生的专业知识水平、实践技能，在微课与翻转课堂中提升学生们的学习兴趣及自主学习能力，转变传统教学模式，营造积极、自由及愉悦的学习氛围，促进学生们专业素养的提升，为社会的发展及跨境电商旅游行业提供高素质的全能型人才。

在X*=192mm处做全局稳定分析，所得特征值分布如图6所示。

  

图6 全局稳定性分析结果：特征值分布Fig.6 Globalstability analysis result: the distrbution of eigenfuction

找到了不稳定的外模态，尚未发现不稳定的内模态。外模态有两个分支，即Y模态和Z模态，二者都存在对称和反对称两种模态。在同一流向位置，Y模态最大增长率高于Z模态，但在低波数段，Z模态增长率略高于Y模态。观察不稳定模态沿流向的演化，发现几个模态一直在相互竞争，没有一个模态始终是最不稳定的。所找到的不稳定模态都出现在基本流沿其梯度方向的Fjortoft拐点处，这表明它们都具有相同的产生机理，即由基本流剪切的无粘拐点不稳定性造成。因为没有符合内模态特性的Fjortoft拐点，所以没有内模态。

3.2 不稳定的外模态种类

在X*=192 mm处做全局稳定分析，图6中有两个分支，都是外模态，分别为Y模态和Z模态(Y、Z模态见文献[15-16])，其特征函数如图7所示。Y模态特征函数在蘑菇状慢速区的“两肩”，沿展向伸展成条带状，表明该模态主要由法向剪切造成。Z模态特征函数在蘑菇状慢速区的“腰部”，为沿法向伸展成条带状，表明该模态主要由展向剪切造成。Y模态又分为两种，分别为对称和反对称模态，特征函数的实部和虚部如图8所示。Z模态也分为对称和反对称模态两种，特征函数如图9所示。

图3为热解不同质量浓度(简称浓度)重镁水溶液所得产物的XRD图，图4为热解不同浓度Mg(HCO3)2溶液所得产物的SEM图。

3.3 在流向同一个位置，不稳定模态随波数的变化规律

在X*=192 mm处做全局稳定分析，改变波数α，观察不稳定模态随其变化，结果如图10所示，图中画出了三个最不稳定的Y模态和三个最不稳定的Z模态。可见，Y模态的最大增长率高于Z模态，但是在低波数(低频率)段，Z模态增长率略高于Y模态。我们在X*=175mm处也做了计算，结果与之类似，如图11所示。

  

图7 全局稳定性分析结果：特征函数(黑实线表示边界层外缘)Fig.7 Globalstability analysis result:eigenfunction (the black solid lines refer to the outer edge of the boundary layer)

  

图8 对称和反对称的Y模态特征函数：左中(右)Fig.8 Eigenfunction of symmetric and anti-symmetric

  

图9 对称和反对称的Z模态特征函数：左中(右)Fig.9 Eigenfunction of symmetric and anti-symmetric

 

  

图10 X*=192mm处增长率随波数α和频率ωr的变化Fig.10 Growth rate varies with wave-number and frequency at X*=192mm

3.4 在固定波数下，不稳定模态沿流向的演化规律

为研究模态沿流向的演化，取固定波数α=0.5，观察不稳定模态沿流向的变化，结果如图12所示，图中画出了三个最不稳定的Y模态和三个最不稳定的Z模态。可见，几个模态一直在相互竞争，没有一个模态始终是最不稳定的。从相速度随流向演化的曲线来看，Y模态相速度较高，主要分布在0.9左右，Z模态相速度与Y模态有明显差别，比其小很多。值得注意的是，在流向X*=210 mm附近，有很多不稳定的模态在此站位交汇，很可能会引起模态间的转换，具体转换机理还得做深入的研究。我们又取波数α=0.7做了计算，结果与之类似，如图13所示。

 

  

图11 X*=175 mm处增长率随波数α和频率ωr的变化Fig.11 Growth rate varies with wave-numberand frequency at X*=175 mm

 

  

图12 不稳定模态沿流向的变化(α=0.5): 增长率ωi和相速度CFig.12 Variation of unstable modes along streamwise(α=0.5): the growth rateωi and phase velocity C

(The contour is the global stability analysis result, the colorcloud image is the direct numerical simulation result)

 

  

图13 不稳定模态沿流向的变化(α=0.7)： 增长率ωi和相速度CFig.13 Variation of unstable modes along streamwise (α=0.7): the growth rate ωi and phase velocity C

3.5 不稳定模态的物理本质

Y模态和Z模态的特征函数在形式上有所区别，但都是由于基本流的剪切造成。由一维基本流的稳定性理论可知，无粘的拐点不稳定性需要满足Fjortoft拐点条件。椭圆锥短轴处基本流流向速度分布是二维的，在Y*-Z*平面内既沿Y*方向有变化，又沿Z*方向有变化。而变化最剧烈的方向，即为基本流流向速度在Y*-Z*平面内的梯度方向，基本流沿这个方向变化最剧烈，即剪切最强。所用Fjortoft拐点条件为：

Calculation Method Study for the Horizontal Bearing Capacity of Offshore Wind Turbine Steel-Pile BI Mingjun,LIU Donghua,WANG Hongqing(77)

(12)

其中∂/∂l为当地基本流在Y*-Z*平面内沿流向速度梯度方向的导数。对X*=192mm处的基本流，将Fjortoft拐点位置用绿色圆圈标记在图14中，相应的Y模态和Z模态的流向速度特征函数也画在图中。可见，特征函数集中的位置与Fjortoft拐点吻合。这表明它们都具有相同的产生机理，即由基本流剪切的无粘拐点不稳定性造成。

由图14还可以看出，“蘑菇”所在区域没有符合内模态特性的Fjortoft拐点，这就解释了为什么我们找不到内模态。图4中红色云图表示梯度方向二阶导数大于零；蓝色表示小于零；白色表示等于零；绿色圆圈表示Fjortoft拐点位置；黑色等值线为特征函数。

  

(a) Y-mode

  

(b) Z-mode

 

图14 流向速度在梯度方向的Fjortoft拐点与特征函数Fig.14 Fjortoft inflection point at the gradient direction ofstreamwise velocity and

4 结 论

图中有两个分支，都是外模态，没有内模态的结果。图6中两个分支上的特征函数如图7所示，都是典型的外模态。其它流向位置结果与之类似，不再赘述。

致谢:感谢国家重点研发计划“大科学装置前沿研究”重点专项“高超声速边界层转捩机理、预测及控制方法研究”(2016YFA0401200)提供资助。

参 考 文 献：

[1]Wheaton M B, Juliano J T, Berridge C D, et al.Instability and transition measurements in the Mach-6 Quiet Tunnel[R].AIAA 2009-3559.

[2]Juliano J Thomas.Instability and transition on the HIFiRE-5 in a Mach-6 quiet tunnel[D].Purdue University.Indiana, US, 2010.

[3]Schmisseur J, Schneider S, Collicott S.Receptivity of the Mach-4 boundary-layer on an elliptic cone to laser-generated localized free stream perturbations[R].AIAA-98-0532, 1998.

[4]Schmisseur J, Schneider S, Collicott S.Response of the Mach-4 boundary layer on an elliptic cone to laser-generated free stream perturbations[R].AIAA-99-0410, 1999.

[5]Poggie J, Kimmel R.Traveling instabilities in elliptic cone boundary-layer transition at Ma=8[R].AIAA-98-0435, 1998.

[6]Huntley M, Smits A.Transition studies on an elliptic cone in Mach 8 flow using filtered Rayleigh scattering[J].European Journal of Mechanics B-Fluids, 2000, 19(5): 695-706.

[7]Lyttle I, Reed H.Use of transition correlations for three-dimensional boundary layers within hypersonic flows[R].AIAA-95-2293, 1995.

周所长仔细一看肯定地说是，喊来站前所民警小张认一认，小张跑进来一看就说：“这是老赵和老常，他们在车站做搬运工已经好多年了的。”

[8]Kimmel R, Klein M, Schwoerke S.Three-dimensional hypersonic laminar boundary-layer computations for transition experiment design[J].Spacecraft Rockets, 1997, 34(4): 409-415.

[9]Gosse R, Kimmel R.CFD study of three-dimensional hypersonic laminar boundary layer transition on a Mach 8 ellipticcone[R].AIAA 2009-4053, 2009.

[10]Poggie J, Kimmel R, Schwoerke S.Traveling instability waves in a Mach 8 ow over an elliptic cone[J].AIAA J, 2000, 38(2): 251-258.

[11]Choudhari M, Chang C, Jentink T, et al.Transition analysis for the HIFiRE-5 vehicle[R].AIAA 2009-4056, 2009.

[12]Paredes P, Theofilis V.Spatial linear global instability analysis of the HIFiRE-5 elliptic cone model flow[R].AIAA 2013-2880.

[13]Paredes P, Theofilis V.Traveling global instabilities on the HIFiRE-5 elliptic cone model flow[R].AIAA 2014-0075, 2014.

相关系数是反映变量之间关系密切程度的统计指标，按积差方法计算，以两变量与各自平均值的离差为基础，以两个离差相乘来体现两变量之间的相关程度。

[14]Zhang Y M, Luo J S.Application of Arnoldi method to boundary layer instability[J].Chinese Physics B, 2015, 24(12): 124701.

[15]Zhang S L.The instability and wave propagation in the hypersonic 2∶1 elliptic cone boundary layer[D].Tianjin University, 2016.(in Chinese)张绍龙.高超声速2∶1椭圆锥边界层的稳定性特征及扰动演化[D].天津大学, 2016.

[16]Li F, Choudhari M M, Duan L, et al.Nonlinear development and secondary instability of traveling crossflow vortices[J].Phys Fluids, 2014, 26: 064104.

双星行星机构的速比验证条件与单行星机构相同,所不同的是由于配齿条件的差异导致行星机构特征参数α不同,致使双星与单星行星机构的实际速比序列会有一定的偏差,这也使整个的离散速比序列分布更密实,使实际速比有更多机会接近理论速比.

[17]Zhang Yongming, Zaki Tamer, Sherwin Spencer, et al.Nonlinear response of a laminar boundary layer to isotropic and spanwise localized free-stream turbulence.AIAA 2011-32[R].Reston: AIAA, 2011.