甲烷平面射流扩散火焰的大涡模拟

本文对甲烷-空气平面自由射流扩散火焰进行了大涡模拟,采用分步投影法求解动量方程,湍流亚格子项采用动态模式模拟,化学反应速率亚格子项采用动态相似模式模拟,压力泊松方程采用修正的循环消去法快速求解,空间方向采用二阶精度的差分格式,在时间方向上采用二阶精度的显式差分格式。模拟结果给出了湍流扩散火焰的瞬态发展变化过程,表明射流扩散火焰的发展过程存在着“湍流控制”和“化学反应控制”两个不同阶段。 “湍流控制”阶段仅存在于火焰发展初期的极短时间内。     

超音速氢气燃烧火焰结构特性的直接数值模拟

在课题组前期工作的基础上,对一马赫数为1.2的三维超音速氢气射流抬升火焰进行了直接数值模拟研究,其中空间离散采用波带优化的四阶WENO格式,时间积分采用带有TVD性质的三步三阶龙格库塔格式,边界条件采用了无反射特征边界条件,总的计算网格数达到9.75亿。结果表明:超音速射流氢气燃烧火焰可分为根部层流状的高温高热量释放率稳燃区、高度褶皱的湍流剧烈混合区和远场燃烧区。火焰自燃稳燃点出现在喷口附近的x/D=0.86处,对应着最易反应混合分数。在此下游,预混燃烧和扩散燃烧两种模式同时存在,其中在剧烈混合区和远场区火焰以扩散燃烧为主,但在火焰根部的局部区域预混燃烧热量释放率达到35%左右。     

脉冲激励下超音速混合层涡结构的演化机理

采用大涡模拟方法对脉冲激励作用下的超音速混合层流场进行数值模拟,所得结果清晰展示了流场中涡结构的独特生长机理.基于涡核位置提取方法,对超音速混合层流场中涡结构的空间尺寸和瞬时对流速度等动态特性进行了定量计算.通过分析流场中涡结构的动态特性在不同频率脉冲激励下的变化,揭示出受脉冲激励超音速混合层流场中涡结构的演化机理:涡结构的生长不再是依靠相邻涡-涡结构之间的配对与融合,而是通过涡核外围的一串小涡旋结构被依次吸进涡核来实现,且受激励流场中各个涡结构的空间尺寸变化较小;流场中的涡结构数量与脉冲频率成正比例关系,而涡结构的空间尺寸与脉冲频率成反比例关系;涡结构的平均对流速度随脉冲频率的增大而减小.针对受脉冲激励超音速混合层,给出了能够表征涡结构特性与脉冲激励参数之间关系的方程式,即受激励流场中涡结构的平均对流速度与脉冲周期的乘积近似等于流场中涡结构的空间尺寸(涡结构平均直径).     

基于FGM的湍流射流扩散火焰超大涡模拟

为了研究超大涡模拟(VLES)预测湍流燃烧问题的能力,本文采用VLES结合基于假定概率密度函数的火焰面生成流型(FGM)建表湍流燃烧模型对值班甲烷/空气湍流射流扩散火焰(Sandia Flame D)开展了高精度数值研究,并与实验结果进行了详细比较.结果表明本文发展的VLES-FGM方法可以较准确地预测出湍流射流扩散火...     

三维超声速射流火焰的直接模拟研究

本文使用了直接模拟（DNS）程序计算了三维超声速湍流射流火焰。首先介绍了本程序使用到的计算方法,包括热物性参数和输运参数的计算、反应机理、边界条件、离散方法。然后,本文使用该方法计算了马赫数为1.2的三维超声速H₂-AIR湍流射流火焰,描述了火焰的整作结构。最后,本文对火焰底部的稳燃的机理—自燃进行了分析和论证。     

湍流射流与扩散火焰大涡拟序结构的波动特性研究

１前言湍流射流扩散燃烧方式提供了射流火焰与涡团相互作用的最基本形式,对研究在湍流射流剪切边界层内的反应物的卷吸混合、热量及动量的输运及湍流射流结构方面具有其特殊的意义。Ｋａｔｔａ［１］利用浮力与单步反应机理模型对Ｎ２－Ｈ２射流扩散火焰进行了直接数值模...     

射流火焰微团特征反应时间

本文试图从非稳定射流火焰的燃烧噪声角度来研究火焰中的微团特性。定义火焰特征反应时间来反映火焰的化学反应动力学和火焰流场流动动力学特性。通过自相关系数测量火焰涡团的特征反应时间，结果表明：相同的燃气流量下，特征反应时间与射流速度成反比；特征反应时间随射流速度的变化斜率与燃气流量相关。     

甲烷湍流射流火焰锋面结构的激光PLIF测量

利用OH-PLIF技术,研究了甲烷空气预混射流火焰,在不同出口雷诺数下,不同氮气稀释比例下火焰OH自由基的二维分布以及火焰的前锋面结构特性,以此实现对柔和燃烧的基础特性研究。实验结果显示,在获得大量氮气稀释后,预混火焰的OH浓度峰值减少了30%以上,并且整个燃烧反应明显延迟,但火焰面面密度没有发生明显变化。同时发现,改变出口雷诺数对火焰OH浓度影响较小,但是增加出口雷诺数可以扩大燃烧反应区。     

煤粉射流火焰的DNS研究

本文用直接数值模拟(DNS)的方法对煤粉气固两相射流燃烧进行了研究。文章介绍了DNS的数值计算方法,包括气相、颗粒相的控制方程以及煤粉燃烧模型。分析了煤粉射流燃烧火焰的结构形态,发现火焰区域中存在三种火焰形态,通过火焰特性的统计研究,发现虽然火焰形态不同,但是其燃烧特性相似,非预混火焰占主导地位。分析火焰释热量,发现上游处单颗粒燃烧释热量较小,条带状释热量占24%,对下游稳定火焰的形成至关重要。     

变密度随机涡模型在湍流射流扩散火焰中的应用

采用变密度随机涡模型,对H₂/O₂/N₂湍流射流扩散火焰进行数值模拟,湍流过程通过涡的采样、涡的抑制和涡的翻转实现.其中,针对变密度反应流问题,提出一种大涡抑制的新机制,并详细讨论各种参数对模型预测效果的影响.计算结果表明,修改后的模型可以合理预测H₂/O₂/N₂射流火焰结构,能够反映湍流的涡特性;模型中与涡采样和涡抑制有关的参数对预测结果有一定影响.     

Topological Structure of Chern-Simons Vortex

Based on the gauge potential decomposition theory and the φ-mapping method, the topological inner structure of the Chern-Sirnons-Higgs vortex has been studied strictly. It is shown that there exits a multi-charged vortex at every zero point of the Higgs scalar field φ. The multivortex solutions in the Chern-Simons-Higgs model are obtained strictly.     

Topological Structure of Vortex Solution in Jackiw-Pi Model

By using o-mapping method, we discuss the topological structure of the self-duality solution in Jackiw-Pi model in terms of gauge potential decomposition. We set up relationship between Chern-Simons vortex solution and topological number, which is determined by Hopf index and Brouwer degree. We also give the quantization of flux in this case. Then, we study the angular momentum of the vortex, which can be expressed in terms of the flux.     

间断有限元方法在弹尾超音速喷流计算中的应用

采用间断有限元方法对超音速无粘喷流流动进行数值模拟.将二维双曲守恒方程的间断有限元方法发展到轴对称Euler方程,并就某导弹尾部超音速伴随射流进行数值计算.计算结果与实验照片反映的流动特征吻合较好,与高精度、高分辨率TVD格式的计算结果相比,间断有限元方法的计算结果在轴线反射点附近具有较高的分辨率,表明该方法对激波具有较强的捕捉能力,在激波阵面上不会产生振荡或抹平间断现象.     

Topological Structure of Vortex Solution in Jackiw-Pi Model

By using φ-mapping method, we discuss the topological structure of the self-duality solution in Jackiw-Pi model in terms of gauge potential decomposition. We set up relationship between Chern-Simons vortex solution and topological number, which is determined by Hopf index and Brouwer degree. We also give the quantization of flux in this case. Then, we study the angular momentum of the vortex, which can be expressed in terms of the flux.     

H2/air旋流预混燃烧的流场特性和火焰结构分析

本文采用完全可压缩的N-S方程,对当量比为1.0的H₂/air旋流预混火焰进行了直接数值模拟研究。氢气和空气的化学反应采用9种组分19步的详细机理。模拟结果表明,强旋流流场中存在回流区,碗形旋流火焰稳定在回流区的外围。在火焰面上沿火焰法向提取了局部火焰结构,将局部湍流火焰结构与层流预混火焰的火焰结构进行了比较,发现局部湍流火焰比层流预混火焰更薄,燃烧强度更高。     

DNS assessment of relation between mean reaction and scalar dissipation rates in the flamelet regime of premixed turbulent combustion

The linear relation between the mean rate of product creation and the mean scalar dissipation rate, derived in the seminal paper by K.N.C. Bray [‘The interaction between turbulence and combustion’, Proceedings of the Combustion Institute, Vol. 17 (1979), pp. 223–233], is the cornerstone for models of premixed turbulent combustion that deal with the dissipation rate in order to close the reaction rate. In the present work, this linear relation is straightforwardly validated by analysing data computed earlier in the 3D Direct Numerical Simulation (DNS) of three statistically stationary, 1D, planar turbulent flames associated with the flamelet regime of premixed combustion. Although the linear relation does not hold at the leading and trailing edges of the mean flame brush, such a result is expected within the framework of Bray's theory. However, the present DNS yields substantially larger (smaller) values of an input parameter c_m (or K₂ = 1/(2c_m ? 1)), involved by the studied linear relation, when compared to the commonly used value of c_m = 0.7 (or K₂ = 2.5). To gain further insight into the issue and into the eventual dependence of c_m on mixture composition, the DNS data are combined with the results of numerical simulations of stationary, 1D, planar laminar methane–air flames with complex chemistry, with the results being reported in terms of differently defined combustion progress variables c, i.e. the normalised temperature, density, or mole fraction of CH₄, O₂, CO₂ or H₂O. Such a study indicates the dependence of c_m both on the definition of c and on the equivalence ratio. Nevertheless, K₂ and c_m can be estimated by processing the results of simulations of counterpart laminar premixed flames. Similar conclusions were also drawn by skipping the DNS data, but invoking a presumed beta probability density function in order to evaluate c_m for the differently defined c's and various equivalence ratios.     

Interaction Mechanism and Loss Analysis of Mixing between Film Cooling Jet and Passage Vortex

The interaction between the film-cooling jet and vortex structures in the turbine passage plays an important role in the endwall cooling design. In this study, a simplified topology of a blunt body with a half-cylinder is introduced to simulate the formation of the leading-edge horseshoe vortex, where similarity compared with that in the turbine cascade is satisfied. The shaped cooling hole is located in the passage. With this specially designed model, the interaction mechanism between the cooling jet and the passage vortex can therefore be separated from the crossflow and the pressure gradient, which also affect the cooling jet. The loss-analysis method based on the entropy generation rate is introduced, which locates where losses of the cooling capacity occur and reveals the underlying mechanism during the mixing process. Results show that the cooling performance is sensitive to the hole location. The injection/passage vortex interaction can help enhance the coolant lateral coverage, thus improving the cooling performance when the hole is located at the downwash region. The coolant is able to conserve its structure in that, during the interaction process, the kidney vortex with the positive rotating direction can survive with the negative-rotating passage vortex, and the mixture is suppressed. However, the larger-scale passage vortex eats the negative leg of the kidney vortices when the cooling hole is at the upwash region. As a result, the coolant is fully entrained into the main flow. Changes in the blowing ratio alter the overall cooling effectiveness but have a negligible effect on the interaction mechanism. The optimum blowing ratio increases when the hole is located at the downwash region.     

DNS study on the formation of Lambda rotational core and the role of TS wave in boundary layer transition

Λ-vortex plays a pivotal role in the boundary layer transition process. In this paper, the generation of Lambda rotational core (Λ-vortex) in boundary layer transition and the role of TS (Tollmien–Schlichting) waves performed are studied by direct numerical simulation (DNS) in detail. It shows that the growth of TS wave under non-linear NS equations differs from the one under linearised NS (Navier–Stokes) equations at very early stage. With the evolution of TS waves, a newly found evolution procedure under non-linear NS equations helps to form the perturbation structures and lead them to evolve to Λ-vortices. The transition process undergoing the non-linear stage is also found to be much earlier than the previous studies had suggested.