超声速凹槽大涡数值模拟研究

在超声速燃烧室内常采用凹槽增强燃料与空气的混合,实现火焰稳定,本文采用大涡模拟方法对开式凹槽流场进行了数值模拟研究,计算采用Smagorinsky 亚网格模型模拟小尺度涡的作用.数值模拟结果给出了凹槽流动的自激振荡的发展过程,以及凹槽流动中拟序结构的演变过程,如涡的卷起、增长,涡与涡之间的合并、破碎等,计算得到的凹槽波动的峰值频率和理论结果一致.     

SLAU格式在大涡模拟方法中的应用

采用格心格式的有限体积法求解三维大涡模拟控制方程,成功地将最新改进的SLAU格式推广到大涡数值模拟方法中。SLAU格式的空间重构方式采用五阶WENO格式,时间推进采用LU-SGS隐式时间推进。湍流模式为最新提出的基于随机分析的动态亚格子模型。计算了圆柱绕流、NACA4412翼型绕流、NACA0012翼型绕流三个算例,分别从可压缩较大迎角翼型绕流、可压缩跨音速翼型绕流以及不可压缩低马赫数圆柱绕流三个方面验证了将SLAU格式应用到湍流大涡数值模拟计算中的可行性。结果表明,结合SLAU格式的大涡模拟方法可以达到较好的数值模拟效果,同时间接检验了基于随机分析的动态亚格子模型的数值模拟效果。     

洁净室流场大涡模拟

应用湍流大涡模拟方法对矢流洁净室和全顶棚送风下部两侧回风式洁净室进行了数值模拟。计算中使用了隐含流线抑风耗散作用的具有三阶部分展开的Taylor-Galerkin离散格式，导出了经Gauss滤波后的基本方程组的有限元列式，应用涡粘性亚格子模式的湍流大涡模拟方法进行数值模拟。数值计算结果与实验结果符合较好，表明湍流大涡模拟能够较准确地捕捉到洁净室流场中旋涡的位置和尺度，在洁净室的数值仿真方面，与标准的K－ε湍流模式相比，湍流大涡模拟结果具有较高的可靠性。     

绕翼型分离流结构的数值研究

利用大涡模拟技术,对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟,详细地给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数。在数值模拟中,采用了弱压缩流的控制方程,用贴体坐标技术进行了网格生成。     

超声速喷流混合流场大涡模拟

以光学窗口外冷喷流为研究背景，采用大涡模拟方法对后台阶外形切向喷流混合流场进行了研究。数值方法基于隐式亚格子模型，采用高精度WENO格式进行空间离散，并通过超声速平面混合层流动对数值方法进行了考核验证。喷流混合流场计算模型与试验一致，来流和喷流马赫数分别为3.4和2.5。数值模拟清晰地捕捉到了流场波系以及混合剪切层、壁面边界层等典型流场结构，并精细预测了混合层发生失稳、转捩及发展为充分发展湍流的时空发展过程。数值模拟得到的湍流大尺度结构的位置和形态与实验图像一致。通过对瞬时流场、统计平均流场和脉动参数的分析，揭示了流场结构特征及其时空演化规律，并获得了流场密度脉动特性。      

分区LES/DES混合方法及缝翼三维流场模拟

提出一种分区大涡/脱体涡模拟混合方法,采用高阶空间有限差分方法和二阶隐式LU-SGS时间推进方法,对有限展长缝翼的三维流场进行了数值计算。该混合方法吸收了脱体涡模拟方法和大涡模拟方法各自的优点,在非核心区域使用脱体涡方法进行计算,相对于大涡模拟方法具有较高的计算效率;相对于分区脱体涡模拟中采用的雷诺平均方法,在计算量没有明显增加的条件下提高了对流场各尺度流动的模拟能力。计算结果与风洞实验数据吻合良好,同时缝翼流场主要观测点的压力脉动数据为后续缝翼噪声分析和低噪声优化提供了基础。     

细长旋成体大迎角俯仰振荡的数值研究

为研究俯仰振荡对细长旋成体流场与气动特性的影响,利用雷诺平均方法与大涡模拟的混合方法相结合的脱体涡模拟方法对大迎角下细长旋成体俯仰运动进行了数值模拟研究。通过对背风面分离涡的强度和位置、细长体绕流形态、截面压强系数和侧向力系数的观察和分析,发现施加小振幅、特定频率的俯仰振荡对大迎角细长旋成体背风面流场有明显的控制能力,...     

超声速横向射流三维流场结构特征

为了研究超声速燃烧室内燃料与空气快速掺混过程的流场特性，基于可压缩Navier-Stokes方程，采用大涡模拟方法和高精度WENO-TCD混合格式对来流马赫数为2.68，喷压比为36的超声速横向射流流场结构进行数值研究。数值结果清晰描述了超声速主流与横向射流相互作用过程的流场结构特征，得到了三维激波形态的演变规律以及它们在强化混合过程中的作用。另外，因桶形激波背风面低压区处的斜压效应，射流气体在桶形激波背风面形成一对螺旋向上的反向涡对，反向涡对的卷吸作用诱导进入壁面边界层的主流向上运动，形成冲击射流。冲击射流以v=557m/s的法向速度向上冲击桶形激波背风面，因而在桶形激波背风面留下类三角锥面凹痕。     

直方管内充分发展湍流的大涡模拟 第二部分

本文第一部分用大涡模拟的方法对直方管内Ｒｅ＝４．９×１０￣４的充分发展湍流运动进行了数值模拟。所得结果证实大涡模拟是高雷诺数复杂湍流运动数值模拟的一个有效手段。在本文第二部分用所得的数据库对湍流运动的结构和形态进行了分析，并且用流体质点运动图形显示了湍流运动的二次流。     

翼身组合体复杂涡系相互干扰的数值模拟

本文应用有限体积法求解Ｅｕｌｅｒ方程，数值模拟翼身组合大攻角、跨声速绕流。成功地模拟了机身体涡、机翼脱体涡、翼尖侧缘涡的形成、发展及其相互间的干扰。     

混合流场控制的大涡模拟

采用大涡模拟方法研究了三维小扰动下,混合流场大尺度拟序结构的产生和演化过程,捕捉了展向涡的卷起、配对、合并,以及二次流向涡的出现等大尺度的三维拟序结构,分析了拟序结构与入口扰动方式之间的内在联系,再现了涡卷自身撕裂而引发转捩的现象。     

基于Kirchhoff方法的亚声速平面混合层主涡对并声场分析

为拓展主涡对并气动声场计算分析方法研究,针对亚声速平面混合层主涡对并流动,采用直接数值模拟(DNS)方法求解获取近场声源,利用二维Kirchhoff频域法和时域法外推声辐射远场,并应用Kirchhoff频域法对主涡二次对并声辐射特征进行分析。研究结果表明:Kirchhoff频域法和时域法具有一定差异性,但与主涡对并直接声计算结果对比,显示二者计算精度相当且均具有较高的精确性。通过Kirchhoff频域法应用分析,清晰揭示出主涡二次对并各声模态与主涡卷起、对并过程的对应相关性,利用频域法这种良好的计算分析性能,本文进一步研究了混合层入口扰动相位差对涡并声场的影响,展现出主导声波在声场中的重要作用。     

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

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

用随机涡法研究混合层中的颗粒扩散

由于 Stuart的涡街模型不能模拟大尺度旋涡的发展与合并过程 ,因而影响了其所得结论的一般性。本文应用随机涡法对混合层中大尺度旋涡的产生、发展与合并过程进行直接模拟 ,定性研究了颗粒在混合层中的扩散 ,并在更一般的意义下获得了与 Crowe相同的结论。     

三维超声速混合层内标量混合的大涡模拟

对三维对流马赫数0.62的超声速混合层流动的标量混合进行大涡模拟, 控制方程对流项采用五阶精度的WENO格式求解, 小尺度涡的作用采用一方程LDKM亚格子模型处理, 过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了混合层流场大尺度拟序结构以及标量场的演化过程, 研究表明标量混合过程受混合层内涡系演化所控制, 标量场具有明显的三维特征.模拟得到的速度、组分及其脉动的统计时均结果和实验结果相符较好.      

Characteristics of reattached boundary layer in shock wave and turbulent boundary layer interaction

The reattached boundary layer in the interaction of an oblique shock wave with a flatplate turbulent boundary layer at Mach number 2.25 is studied by means of Direct Numerical Simulation(DNS). The numerical results are carefully compared with available experimental and DNS data in terms of turbulence statistics, wall pressure and skin friction. The coherent vortex structures are significantly enhanced due to the shock interaction, and the reattached boundary layer is characterized by large-scale...     

亚声速等膨胀平面射流的初始流场结构

基于大涡模拟方法与5阶精度混合TCD/WENO格式,对等膨胀射流在马赫数为Ma=0.6条件下的初始流场结构进行了数值模拟。数值结果描述了射流等膨胀过程中,射流初始流场的形成、变化及发展过程。并发现了射流初始流场中涡核为射流失稳的起始点,研究了射流混合层与涡核剪切层上涡的产生、发展与破碎以及它们卷吸环境气体并与射流气体混合的整个过程,揭示了射流初始过程中涡核与涡环直径的变化趋势,得到了射流轴向速度与压力的分布。数值结果可为相关发动机的设计提供重要参考。     

Planar laser scattering visualization of streamwise vortex pairs in a Mach 6 flow

A series of cross-sectional flow fields of Counterrotating Vortex Pairs(CVPs) generated by a large-scale ramp vortex generator is observed using an ice-cluster-based Planar Laser Scattering(PLS) method in a shock tunnel with a nominal flow Mach number of 6. Combined with a numerical simulation, two streamwise CVPs with opposite rotating directions are identified in the wake flow of the vortex generator with an absence of a boundary layer, namely, a Primary CVP(PCVP) and a Secondary CVP(SCVP). Th...     

交叉剪切混合层中拟序结构的数值研究

采用拟谱方法对时间模式的交叉剪切混合层进行了直接数值模拟。计算结果表明:与平面混合层一样,展向涡的拉伸作用是交叉剪切混合层中流向涡形成的主要机制。当展向剪切强度较大时(如两主流交叉角为40°),与初期展向KelvinHelmholtz相关的单向旋转流向涡在拉伸作用下很快增长起来,并“坍缩”成“肋状”涡。当交叉角为40°时,涡核区存在类似平面混合层中“方块状”涡的流向涡结构,展向涡辫区还存在一组符号相反的流向涡,不过与“肋状”涡对应的涡结构呈扁平状,始终没有“坍缩”。当交叉角为60°时,“肋状”涡非常强,以致完全抑制了平面混合层“对称模式”的发展。当交叉角小到20°时,流向涡结构更接近于对称分布,然而“肋状”涡却没有形成。另外,计算结果还证实:与二维混合层相比,大强度展向剪切的引入能够加强流场的混合,同时,适当增加展向扰动波初始强度和波数也是提高混合效率的有效手段。     

DBD涡流发生器及其在角区流动控制中的数值研究

为明确介质阻挡放电(DBD)涡流发生器对马蹄涡的影响,在采用唯象模型的基础上,通过数值方法研究了DBD涡流发生器诱导产生流向涡的结构特性及其对马蹄涡的控制特性。研究结果表明,流动在激励器上游边缘处形成羊角涡,自由剪切层卷入流向涡的涡核并为其提供持续涡量;在柱体根部角区流动中,当对称面两侧激励器诱导流动指向对称面,诱导涡与马蹄涡环绕方向相反时,马蹄涡可以得到有效抑制,反之,则控制效果不佳。最后得出,诱导涡对下游马蹄涡的控制机制体现在其黏性扩散作用、掺混作用以及低压效应3个方面。