火烧炉内气相燃烧过程的三维数值模拟

本文对某工业火烧炉内的三维湍流气相燃烧场特性进行了数值模拟,湍流模型采用可实现k-ε双方程模型,湍流燃烧采用有限速率/涡耗散模型,辐射换热采用P1模型,数值方法采用SIMPLE算法,实现了炉内火烧试验的数值再现,特别考察了入口压力、燃料/空气比、多烧嘴不同组合工况等因素对炉内温度场和流场的影响,并对该火烧炉的运行提出了一些合理化建议,为其实际运行提供了有益的参考。     

高强驻波声场中的湍流

反射端声压级高于140dB以后驻波管中的声流结构会发生明显变化。利用流动显示和LDV（Laser Doppler Velocimetry）技术对该问题进行了实验研究。流动显示结果表明，在160dB左右的高强驻波声场中随着旋涡结构的合并、破裂等现象的出现，声致流动有明显的湍流特性。LDV测量结果显示，随着驻波管中声压级的提高，一些测点的速度频谱明显加宽，出现高频速度脉动，而另一些点的速度频谱则没有显著变化。径向湍流强度测量结果显示，此时管中流动已是复杂的三维结构,而非轴对称的二维结构。     

Using gamma-ray burst prompt emission to probe relativistic shock acceleration

It is widely accepted that the prompt transient signal in the 10 keV–10 GeV band from gamma-ray bursts (GRBs) arises from multiple shocks internal to the ultra-relativistic expansion. The detailed understanding of the dissipation and accompanying acceleration at these shocks is a currently topical subject. This paper explores the relationship between GRB prompt emission spectra and the electron (or ion) acceleration properties at the relativistic shocks that pertain to GRB models. The focus is on the array of possible high-energy power-law indices in accelerated populations, highlighting how spectra above 1 MeV can probe the field obliquity in GRB internal shocks, and the character of hydromagnetic turbulence in their environs. It is emphasized that diffusive shock acceleration theory generates no canonical spectrum at relativistic MHD discontinuities. This diversity is commensurate with the significant range of spectral indices discerned in prompt burst emission. Such system diagnostics are now being enhanced by the broad-band spectral coverage of bursts by the Fermi Gamma-Ray Space Telescope; while the Gamma-Ray Burst Monitor (GBM) provides key diagnostics on the lower energy portions of the particle population, the focus here is on constraints in the non-thermal, power-law regime of the particle distribution that are provided by the Large Area Telescope (LAT).     

数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

高主流湍流度下密度比对涡轮导叶全气膜冷却特性的影响

为了研究高主流湍流度下二次流密度比对涡轮导叶全气膜冷却特性的影响,使用热色液晶测量了在主流湍流度为15%,二次流密度比为1.0和1.5下三维涡轮导叶的气膜冷却效率和换热系数。二次流与主流质量流量比为7.0%和12.5%。结果表明:二次流密度比增大可以降低冷气射流的动量,小流量比工况下,在叶片前缘和压力面前半段,动量较低的二次流在高主流湍流度的影响下更易耗散,增大二次流密度比使冷却效率明显降低;大流量比工况下,二次流动量降低使气膜孔后区域冷气贴附性增强,气膜冷却效率和冷气覆盖效果均得到提升。小流量比工况下,二次流密度比增大对叶片表面换热的影响较小;大流量比工况下,二次流密度比增大使吸力面中弦区域和压力面后半段的平均换热系数比分别降低15%和25%。     

均匀各向同性大气湍流对声爆传播特性的影响

基于离散Fourier模态有限和生成的随机大气湍流场，采用修正波形参数方法，开展了均匀各向同性大气湍流对典型超声速客机声爆传播特性的影响分析。计算采用的超声速客机模型为自行生成的简化超声速公务机模型。首先，应用航空工业空气动力研究院自主研制的CFD软件ARI_Overset在三维空间求解Navier-Stokes方程，得到作为声爆远场传播初始值的近场空间压力分布；其次，基于Von Karman能量谱，采用离散Fourier模态有限和形式生成随机均匀各向同性大气湍流场；最后，采用修正波形参数方法模拟了声爆信号在随机速度湍流场中的传播过程。数值结果表明：各向同性大气湍流对于地面声爆特征有重要影响。与无湍流状态相比，均匀各向同性大气湍流使得地面声爆特征增强的概率约为55%，使得地面声爆特征减弱的概率约为45%，故总体而言大气湍流效应更倾向于增强地面声爆特征；均匀各向同性大气湍流对于声爆传播路径影响相对较小，但是这一变化仍会导致地面信号接收点的不确定性。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

不同钝体宽度下的驻涡燃烧室排放试验

为了研究空气流量分配对驻涡燃烧室对排放特性的影响,了解对驻涡燃烧室内污染物生成的过程及其影响因素,设计了一个能够改变中心钝体宽度、仅凹腔供油的驻涡燃烧室.在常压下对该驻涡燃烧室进行了排放特性试验,进口温度保持200℃.试验中,燃烧室进口马赫数为0.15~0.3.影响排放的因素主要包括雾化质量、凹腔当量比以及与进口马赫数相关的驻留时间等.总体来说雾化质量、凹腔当量比的提高对降低CO和HC的排放是有利的,但是这会使NOx排放增加.在低凹腔当量比时,CO排放曲线变化下降比较平缓,甚至出现上升趋势,而HC排放曲线比较陡峭.这是由于HC的消耗速度比CO消耗速度快,随着凹腔当量比的增加,供油压力提高,燃油雾化粒径变小,燃油蒸发时间缩短,使HC排放快速减少,中间产物CO大量产生而来不及消耗.凹腔当量比进一步上升时,由于燃烧温度的提高,使得CO排放快速减少.在燃烧室内燃烧过程中,NOx的形成和消耗是非常复杂的过程,目前只能作一些定性的分析,而CO和HC的反应过程相对简单.通过对不同钝体槽宽下,具有相似凹腔前壁流量的工况的比较,发现CO和HC的形成主要受凹腔内工作状况影响,而NOx的形成过程更复杂,主流也对其产生着重要的影响.      

New constitutive relation for eddy viscosity models

A new turbulent constitutive relation was directly derived from Boussinesq's hypothesis and mixing length theory, and then implemented in the standard k-ε model. The performance of this constitutive relation was validated in zero pressure gradient flat-plate boundary layer flow, fully-developed turbulent channel flow and separated flow in a plane asymmetric diffuser. The investigation demonstrated that, this new constitutive relation gave very accurate results in the former two basic cases and provided significant improvement in prediction of separated and reattachment points in the plane asymmetric diffuser. Separation and reattachment points at x/H=7.5 and 29 were calculated accurately in comparison to experimental results, and the static pressure coefficient of 0.82 was very close to large eddy simulation calculation. These results are very encouraging but further verification and extensive application of the new constitutive relation to other two-equation eddy viscosity model are needed.      

基于自然激励技术的颤振边界预测

为了预测紊流激励条件下机翼的颤振边界,基于自然激励技术提取紊流响应的自由衰减信号,采用矩阵束方法识别模态参数,最后通过Z-W(Zimmerman-Weissenburger)方法计算稳定性判据,拟合判据变化曲线并外推颤振边界.对平板机翼模型进行了数值仿真分析,对单独机翼模型风洞颤振试验数据进行了计算.结果表明:采用自然激励技术与矩阵束方法能够较准确地识别紊流激励响应的模态参数,频率识别误差小于6％,阻尼比识别误差小于30％,结合Z-W方法能够在较低风速较早地预测颤振边界,有助于提高试验的安全性.