圆柱突然起动的快速涡解法

 文中将Chorin的涡方法、格子涡(VIC)方法和快速Fourier变换(FFT)相结合,发展了一套快速有效的求解高Re数下非定常N-S方程的算法(快速涡方法)。通过对二维圆柱突然起动问题的数值计算,得到流动分离及其演化过程(包括二次涡及二次涡对的生成);此外还得到回流区对称轴线上速度分布。这些结果与实验数据及其它数值解相比,是令人满意的。该程序对其它高Re数非定常流动问题也是适用的。     

宽速域翼尖涡及其与斜激波相互作用

针对宽速域飞行器产生的翼尖涡及其与斜激波相互作用问题,采用数值模拟方法,在来流马赫数0.2～6.0范围内,研究了翼尖涡的演化特性,并探讨了斜激波诱导翼尖涡破碎的机理。结果表明,在宽速域范围,翼尖涡的切向速度和环量分布仍然符合低速流动中的自相似关系式;在超声速至高超声速流动中,翼尖涡强度沿流向更快地减弱;当翼尖涡与斜激波相互作用时,涡核处压力越低的翼尖涡更容易破碎,而经典的涡破碎理论难以在宽速域范围内准确预测涡破碎现象;通过引入涡核处的压力亏损效应修正经典理论,大幅提升了预测斜激波诱导翼尖涡破碎的能力。     

模拟有限长涡片旋转的一种新数值模型

本文把涡格法中的超收敛性推广到自由涡片的演化问题，涡片的演化由ＢｉｒＫ－ｈｏｆｆ－Ｒｏｔｔ方程所控制，该方程具有强烈的非线性，难以用解析法求解。在数值模拟中通常离散涡方法。本研究表明对自由涡片的旋转问题在精确解和离散涡方法结果之间的误差主要来自点涡近似。使用基于超收格式的新数值模型可以显著地提高数值解的精度。     

Multhopp网格在弯曲涡层计算中的超收敛性问题

本文从理论上分析了Ｍｕｌｔｈｏｐｐ网桥弯曲涡面积分方程求解中的数值方法和超收敛性问题;并通过算例验证,证实Ｍｕｌｔｈｏｐｐ网格应用于弯曲涡层时,的确同样具有远大于其它任何网格计算的超收敛速度。     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。     

随机涡法及对钝体后气流结构的模拟

一、问题与方法 本文应用随机涡法,对三棱柱在无界均匀不可压流场中的绕流问题进行了计算。图1为计算简图。      

等离子体涡电磁散射特性及隐身性能

 针对非均匀等离子体在飞行器隐身中的应用,采用分段线性电流密度递归卷积时域有限差分（PLJERC-FDTD）方法计算等离子体涡及涡串电磁散射特性,分析等离子体涡对飞行器隐身性能影响。计算表明,等离子体涡在很大频率区间对电磁波吸收效果显著,RCS降低很大,具有明显的隐身效果。等离子体涡表现出一定规律性的极化特性,对L,S和C波段电磁波具有不同的吸收、反射特性。     

固体火箭发动机转角涡脱落分类

应用三维大涡模拟(LES)方法对C1x模拟试验器进行了数值模拟,以研究转角涡脱落引起的声涡耦合现象.结果表明:与二维和轴对称LES相比,三维LES方法可以准确捕获转角涡脱落运动的三维特性及锁频现象.根据转角涡脱落是否在腔内形成声反馈循环过程可将其分为两类问题,第1类形成声反馈循环涡;第2类形成自耗散涡.燃烧室的台阶扩张比、台阶前后长径比是影响两类转角涡脱落现象的重要结构参数.     

旋流预混燃烧室流动混合的大涡模拟

采用大涡模拟的方法研究旋流预混燃烧室中的流动及燃料空气混合过程,选取亚格子湍动能模式为基准的涡黏模型,通过模拟涡核的运动过程定性的表征了拟序结构的发展;引入混合不均匀度的概率密度函数,研究入口边界条件尤其是入口速度对于涡团破碎及混合程度的影响规律;通过研究速度分量与燃料质量分数脉动的相关性,分析了不同方向的速度脉动对于混合过程的影响因子.这些均为进一步研究旋流预混燃烧室的燃烧不稳定性及污染物排放提供了理论上的指导.      

Kármán涡街中旋涡三维变形的初步研究

本文用涡动力学及LIA方法作为基本理论模型,数值研究了尾迹中孤立涡和细涡丝的三维演化规律。结果表明,圆柱分离尾迹中的旋涡存在三维不稳定性。卡门涡在尾迹平均流场中演化产生三维的类似于马蹄形-勺子形的流场结构。细涡丝在涡辫区的三维演化形成趋向流场最大拉伸变形方向的流向涡结构。     

七叶大侧斜螺旋桨空泡特性及梢涡演化数值模拟研究

为研究七叶大侧斜螺旋桨尾流场及梢涡特性,本文基于DDES（延迟分离涡方法）建立了螺旋桨空化流场数值预报模型。为验证所建立数值模型的准确性,进行了多套不同尺度网格的不确定性分析,同时将非定常流动中螺旋桨（VP1304）螺旋桨空泡及梢涡特性计算结果与试验结果进行了对比。随后基于该模型对七叶大侧斜螺旋桨的尾流场及梢涡特性进行了数值分析。计算结果表明：结本文所建立数值模型精度较高,可以准确地捕捉到梢涡空泡及梢涡尾流场特性;同时DDES方法相比RANS方法在对复杂湍流流动的捕捉能力更强,更适用于螺旋桨梢涡捕捉;尾流场网格加密对尾流场模拟及准确捕捉梢涡十分重要,但对螺旋桨水动力性能影响不大;尾流区域的轴向速度场可以分为加速流动区和自由流动区,进速系数越小,自由流动区与加速流动区之间的界限向外扩张越明显;E1619桨在重载工况有明显的梢涡空泡产生,而轻载工况空泡面积较小且无梢涡空泡发生;梢涡在向下游发展过程中会发生相互融合,进速系数越大,融合发生的越晚,梢涡强度也越小;七叶大侧斜螺旋桨在尾流区域会产生一个分支涡,分支涡起始于吸力弯曲面下缘,且与梢涡呈一定的夹角,进速系数越大,夹角越小。     

湍流的耗散及弥散相互作用理论

 湍流具有涡团散裂、耗散和弥散的特性。根据Prigogine倡导的耗散结构理论,推导了表征耗散与弥散相互作用的新的湍流控制方程组。其特点是:用稳定性分析得到湍流动能产生项,再根据广义熵增原理推出并列存在的,分别适用于强弱涡量的二个湍流动量方程。运用该理论已成功地计算了一些典型的湍流问题:湍流边界层中的马蹄涡拟序结构、钝体涡尾区的湍流能量逆转、湍流涡团散裂弛豫及各向异性分布。文中给出了部分算例。     

湍流边界层中马蹄涡的实验研究

采用~种新的流场显示方法(激光片光运动法)和不同实验方案对湍流边界层中的马蹄涡(或称发卡涡)进行了实验研究,其中包括马蹄涡的生成和发展过程。实验结果表明,湍流边界层中马蹄涡的~种主要生成方式是二次不稳定。其生成和发展过程既不同于层流中人工生成的马蹄涡;也不同于目前存在于人们观念中的那种完美的马蹄涡。这种马蹄涡有以下特征:(1)涡头部分展向尺度不断线性增长,最后固定在某个尺度;(2)涡头部沿流向作匀速运动;(3)流向涡管存在拉伸和变形。     

计算超声速战术导弹大攻角气动力的数值方法

本文给出了超声速大攻角情况下战术导弹气动力特性及压力分布的计算方法。该方法以Woodward,F.A.等人的有限基本解方法为基础。考虑了弹身头部涡、翼面前缘涡、翼面侧缘涡、后缘涡及诸运动涡系的影响,使其应用范围超出了小攻角的限制。 本方法适用于有、滚转有任意舵面偏转情况下正常式或前控导弹外形的气动力特性及压力分布计算。攻角范围直到25°。     

确定性涡方法的精确度研究

本文以无界流场中圆形涡扩散问题为模型问题，对几种为典型的确定性涡方法进行了讨论，分析了方法的参数条件，数值结果表明，确定性涡方法比随机涡方法有更高的精度。     

基于DMD方法的翼型大迎角失速流动稳定性研究

基于剪切应力输运湍流模型的SST-DDES混合方法对NACA0012翼型大迎角分离流动进行非定常数值模拟,采用动力学模态分解(Dynamic Mode Decomposition,DMD)数学工具对失速初始状态、浅失速状态以及深失速状态的流场进行稳定性分析。结果表明:DMD方法准确地提取了翼型大迎角流动中的主频和高阶倍频及对应的流场模态结构;与FFT分析结果相比,频率最大差异小于0.16%;且发现两者提取的频率在流动中的主导作用顺序也一致。通过特征值对相应的模态进行稳定性分析,所有模态的放大率均非常小,所有模态处于弱发散、弱收敛或稳定极限环状态。DMD提取的一阶模态主要表现为分离涡演化过程中最主要的静止分离涡结构,前三阶低频对应的模态涡结构与流动中以此频率进行演化的涡结构比较一致,更高阶的倍频主要表现为尾涡和尾迹区的涡结构。且发现不同模态系数之间存在相位差,说明分离涡流动中不同频率对应的涡结构运动不同步。     

一个扩压叶栅角区流向涡的物理模式

 本文采用烟显示和三维流场测量的方法,研究扩压叶栅角涡的物理模式。试验结果揭示了角涡发生、发展程中的流动形态和结构特点,以及角涡与叶面附面层相互作用的物理过程。在这个物理模式的基础上讨论了角涡引起高流动损失的物理机制。     

低雷诺数下横流-射流中剪切涡的试验研究

为深入分析横流-射流(JICF)的流动特性及其中的复杂涡系结构,从流动机理上研究燃机叶片气膜冷却,揭示高温燃气流与冷却流的掺混机理,本文对横向流中单孔射流所形成的剪切涡开展了试验研究。主要研究了速度比、雷诺数及射流角对JICF所形成剪切涡的影响。结果表明:速度比、雷诺数以及射流角会改变主流与射流之间的掺混程度,从而改变射流轨迹的曲率、高度及垂向渗透能力,最终改变剪切涡的特性;迎风涡与背风涡分别是由射流边界层涡与主流边界层涡形成的,当主流边界层涡强度大于射流边界层涡时,背风涡是流场的主导结构,反之,迎风涡将成为流场的主要涡系结构。     

本征正交分解法对正十二烷喷雾射流的解析

应用大涡模拟方法对航空煤油单组分替代燃料正十二烷冷态喷雾现象进行模拟计算,探究喷雾射流中的多物理场耦合现象。从宏观层面分析燃油破碎蒸发及油气混合过程,发现流场涡团在其中起着重要的作用。运用本征正交分解方法对正十二烷射流质量分数场和涡量场数据库进行加工,从“微观”层面深入探究两者的联系。结果表明:正十二烷质量分数场喷雾结构与射流涡量场涡团结构存在显著的空间分布相关性,反应出两者结构演变具有重要关联特征。湍流脉动致使涡团尺度发生变化,从而影响燃油破碎蒸发,最终使得喷雾不同尺度结构分布形态发生改变。      

自由剪切层中的三维不稳定性

本文是在文[1]的基础上研究自由剪切层中由Kelvin-Helmholtz不稳定波发展而形成的展向大涡结构的三维不稳定性。以大涡结构为基本流动,将稳定性分析归结为二维特征值问题,用pseudo spectral(伪谱)方法数值求解。研究发现:在没有亚谐波存在的情况下,大涡结构的最不稳定的扰动波是流向波长和其相同。展向波数较高,有对流特性的三维扰动波。它在剪切层中的发展与展向涡量的分布有关,大涡结构的涡核不稳定性和辫子不稳定性是流向涡形成的主要力学机制。本文还给出了不同雷诺数下三维扰动波增长率与展向波数的关系,这些结论与实验及数值模拟结果基本一致。