激波诱导火焰变形、混合和燃烧的数值研究

为了深入研究激波诱导的火焰变形以及由此带来的混合和燃烧变化特性,采用带单步化学反应的Navier-Stokes方程和高网格分辨率,对平面入射激波及其反射激波诱导球形火焰变形的现象进行了二维数值研究,计算结果与实验结果较好吻合。研究结果显示,在反射激波作用前,火焰的变形和皱褶主要受入射激波诱导等物理过程影响;而在反射激波与失稳火焰作用后,燃烧放热率、火焰有效面积和界面长度均迅速增加,控制火焰变形的机制逐渐向化学反应(燃烧)过程过渡;在失稳火焰发展的后期,增强的燃烧过程能够削弱火焰界面的皱褶,进而抑制未燃气和可燃气的混合。由此可以得出结论,激波诱导的火焰界面通过变形可促进界面两侧未燃气与可燃气的混合,进而强化燃烧过程,但燃烧的增强却反过来能抑制混合;认识两者之间的关系有助于利用或控制激波 火焰相互作用过程。     

激波冲击火焰的涡量特性研究

激波冲击火焰的现象涉及一系列复杂的物理化学过程,其中涡量的生成与演化对控制火焰发展起重要作用。为系统分析激波冲击火焰过程中的涡量特性,采用二维带化学反应的Navier-Stokes方程对平面入射激波及其反射激波与球形火焰作用的现象进行了数值研究,通过引入并行计算达到高网格分辨率的要求。计算结果表明,斜压项对火焰区内涡量生成起主导作用,压缩项和耗散项在火焰膨胀阶段抑制涡量生成,此外,火焰在激波压缩阶段主要受物理过程而非化学反应过程影响。     

聚心火焰与激波相互作用的数值研究

基于带化学反应的2维轴对称Euler方程,利用带有monotonized centered(MC)限制器的波传播算法,在两端敞开的圆桶中对惰性介质的聚心激波和氢气-空气混合物的聚心火焰与激波的相互作用进行了数值模拟。数值结果表明,在惰性介质中激波在轴心的每次汇聚均可成长为马赫干,马赫干的追赶使激波得到一定程度的增强,但整体还呈下降趋势。在氢气-空气混合物中,燃烧诱导的激波,由于与火焰的反复作用,使激波在轴心处产生马赫干的频率和强度皆高于惰性介质中的情形。同时,火焰在与激波的相互作用过程中发生失稳变形,使其形状呈扁平头部的蘑菇云。     

激波诱导火焰失稳与爆轰的条件研究

采用九阶WENO和十阶中心差分格式数值求解激波与火焰作用过程,考察了激波强度、火焰尺寸对激波与球形火焰作用过程的影响。结果表明,增大激波强度或火焰尺寸均可在流场中引发爆轰,但激波强度的影响更大,并且其引发的爆轰可使火焰迅速膨胀,放热率提高,从而影响燃烧特性;此外,爆轰波传播过程中会迅速消耗可燃预混气,合并原有的反射激波,并在流场中形成局部高压区,极大地改变流场结构。     

INTERFACE EVOLUTIONS AND GROWTH PREDICTIONS OF MIXING ZONE ON PREMIXED FLAME INTERFACE DURING RM INSTABILITY 1)

预混火焰界面的RM (Richtmyer-Meshkov)不稳定导致的界面混合区增长过程在自然界和工程实践中十分常见,但化学反应对其增长的影响机理仍不明确,反应性界面混合区增长速率的预测也未见报道, 因此,开展RM不稳定过程中火焰界面演化和混合区预测的研究十分必要.本文采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了正弦形预混火焰界面在平面入射激波及其反射激波作用下的RM不稳定过程.结果表明, 在入射激波作用后的阶段,除RM不稳定本身导致的界面演化为"钉-帽"和"泡"形结构外,化学反应一方面以预混火焰传播的方式促进了界面中"泡"结构的增长,另一方面通过与涡结构的复杂相互作用促进了"钉-帽"结构的增长.化学反应活性越强, 火焰界面的"泡" 结构和"钉-帽"结构的增长越快.在第一次反射激波作用后的阶段,化学反应以相同的火焰传播方式对"泡"和"钉-帽"结构产生影响, 两者效应相抵,因而导致反射激波作用后的阶段中界面混合区增长不受化学反应活性的影响.根据以上分析,分别针对入射激波和第一次反射激波作用后的火焰界面混合区增长速率提出了相应的预测模型,为探索反应性RM不稳定过程的理论预测方法提供了有益参考.     

RM不稳定过程中预混火焰界面演化及混合区增长预测

预混火焰界面的RM (Richtmyer-Meshkov)不稳定导致的界面混合区增长过程在自然界和工程实践中十分常见,但化学反应对其增长的影响机理仍不明确,反应性界面混合区增长速率的预测也未见报道, 因此,开展RM不稳定过程中火焰界面演化和混合区预测的研究十分必要.本文采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了正弦形预混火焰界面在平面入射激波及其反射激波作用下的RM不稳定过程.结果表明, 在入射激波作用后的阶段,除RM不稳定本身导致的界面演化为"钉-帽"和"泡"形结构外,化学反应一方面以预混火焰传播的方式促进了界面中"泡"结构的增长,另一方面通过与涡结构的复杂相互作用促进了"钉-帽"结构的增长.化学反应活性越强, 火焰界面的"泡" 结构和"钉-帽"结构的增长越快.在第一次反射激波作用后的阶段,化学反应以相同的火焰传播方式对"泡"和"钉-帽"结构产生影响, 两者效应相抵,因而导致反射激波作用后的阶段中界面混合区增长不受化学反应活性的影响.根据以上分析,分别针对入射激波和第一次反射激波作用后的火焰界面混合区增长速率提出了相应的预测模型,为探索反应性RM不稳定过程的理论预测方法提供了有益参考.      

多次激波诱导正弦扰动预混火焰界面失稳的数值研究

激波诱导火焰失稳是实际中常见的现象,为深入研究火焰失稳特性,采用三维单步化学反应的Navier-Stokes方程和9阶weighted essentially non-oscillatory (WENO)的高精度格式,对不同马赫数的入射激波及其反射激波多次诱导正弦型预混火焰界面失稳的现象进行了三维数值模拟,并对计算结果的可靠性进行了验证。研究结果显示,在激波的多次作用下,火焰界面的演变过程主要受Richtmyer-Meshkov (RM)不稳定特性和化学反应特性的双重影响,且随着入射激波强度的增强,上述2种特性均得到进一步强化。为构造体现反应性RM不稳定特性的参数,根据火焰界面混合区平均涡量和化学反应速率,提出了表征界面受不稳定性和化学反应影响的量纲一参数。通过分析发现,在同一入射激波强度下,该参数的对数形式随入射激波和反射激波的多次作用呈基本相同的线性增长趋势;对不同马赫数的入射激波,该参数对数形式的线性增长率也基本一致。这样的变化表明该量纲一参数能够反映反应性RM不稳定过程中火焰界面发展的内在规律。     

铝粉湍流火焰诱导激波现象的实验研究

粉尘湍流火焰诱导激波问题是工业灾害研究中的重要课题．本文在自行设计的大型卧式燃烧管内，对铝粉火焰诱导激波现象进行了实验研究，测试了湍流火焰阵面前压缩波到激波的转捩过程，并将实验结果与数值模拟结果进行了比较．     

两相介质中火焰诱导激波现象的研究

基于双流体模型，采用跟踪法处理火焰内边界，ＴＶＤ格式计算气相激波，ＭａｃＣｏ－ｍａｃｋ格式计算颗粒相流场，对两相可燃介质中火焰诱导激波现象进行了理论分析与实验模拟，研究了颗粒相参数对火焰诱导激波现象的影响。     

氢燃料缓燃向爆震转捩过程中波与火焰的匹配特性

为了解氢燃料爆震过程中压力波与火焰之间相互匹配的特性,在60mm60mm2000mm 方 爆震管内,用氢气和空气混合物进行了单爆震性能研究。在爆震转捩区内布置压力传感器与离子探针,用来 监控压力波和火焰的信号,同时利用高速摄影仪集中拍摄转捩区域。根据压力波和火焰面在爆震管不同时刻 的强度特性、速度特性及位置特性来分析爆震过程中波与火焰匹配的规律。结果表明:压力波和火焰的强度 呈现为相互正反馈匹配性质;缓燃向爆震转捩(DDT)过程中,压力波和火焰的速度表现为相互交替的变化过 程,且缓燃阶段中火焰速度的增幅大于压力波速度的增幅;当火焰面追赶上激波时,产生过爆,火焰面会临时 位于激波前面;在过爆衰减为正常爆震波的过程中,激波在火焰前面。     

Simulation of spatial multibody systems with friction

A formulation for modeling and simulation of friction effects in spatial multibody systems is presented. Constraint reaction forces on rigid bodies that are connected by joints that support friction are derived as functions of Lagrange multipliers, using D’Alembert’s principle. Friction forces acting on bodies are calculated as a function of joint geometry, constraint reaction forces that are functions of Lagrange multipliers, and relative velocities at constraint contact points that are determined by system kinematics. Friction forces are implemented in index 0 differential-algebraic equations of motion that are solved numerically using explicit and implicit numerical integration methods. Spatial examples are presented, yielding accurate results and demonstrating that the systems are not stiff, even in the presence of friction and stiction.     

正交各向异性厚板的边界元解法

本文利用 Hormander 算子法和平面波分解法导出了计入剪切变形的正交各向异性厚板的基本解。建立了计入剪切变形的正交各向异性厚板的边界积分方程。文中详细地讨论了基本解的数值计算,并用边界元法分析了一些算例。     

张拉膜结构力密度法混合找形分析

论述了张拉膜结构力密度法混合找形的基本理论;所谓力密度法混合找形,即部分单元力密度控制,部分单元弹性控制;力密度控制采用线性求解,弹性控制采用非线性求解,通过迭代计算混合找形求出各结点的坐标值。据此编制了相应的计算软件;对工程实例进行了验算,结果表明,本文给出的计算结果与德国著名软件easy的计算结果相吻合。     

Effects of initial stress on the reflection and transmission waves at the interface between two piezoelectric half spaces

The effects of initial stress on the reflection and transmission waves at the interface between two piezoelectric half spaces are studied in this paper. First, the secular equations in the traverse isotropic piezoelectric half space are derived from the general dynamic equation with initial stress taken into consideration. Then, the interface conditions that displacement, stress, electric potential, and electric displacement are continuous across interface are required to be satisfied by three sets of coupled waves, namely, quasi-longitudinal wave, quasi-transverse wave and the electric–acoustic wave. The algebraic equations resulting from the interface conditions are solved to obtain the amplitude ratio of various waves and furthermore the energy reflection and transmission coefficients of various waves. The numerical results are shown graphically and the effects of initial stress are discussed.     

根部柔性梁的不确定性建模与确认

首先导出了包含根部挠性参数的悬臂梁动力学方程近似解,将平动和转动柔性参数作为不确定性源,根据若干预紧力矩下的模态实验结果完成了柔性参数识别,并假定其服从正态随机分布,识别了均值与标准差,进而通过回归分析建立了不确定参数随预紧力矩变化的数学模型.最后通过新的模态实验结果对所建模型进行了确认,通过确认结果,明确了该模型的使用范围.     

Supersonic Inviscid Corner Flows with Regular and Irregular Shock Interaction

The results of a numerical and analytical investigation of steady-state supersonic inviscid flows in corners formed by intersecting compression wedges are presented. The flows considered are symmetric about the corner bisector. The distinctive features of flow pattern formation related with the reflection of wedge-generated shocks from the bisector plane are studied. The wedge angles at which transition from regular to irregular shock reflection occurs are determined both numerically and analytically using the criteria available for plane flows; the data thus obtained are found to be in agreement. Flow patterns with irregular shock reflection, namely, single, transitional, and double Mach, as well as von Neumann reflection, are identified; they are similar to the known types of reflection for plane quasi-steady-state flows. Varieties of these types not observed in the plane flows are found to exist. The effects of the angle of inclination of the plane surfaces of the corner to the freestream direction, the sweep angle of the leading edges, and the dihedral angle are investigated. Some previously unknown parameters of corner configurations for which transition may occur in accordance with the von Neumann criterion are determined.     

Applications and developments of aeroelasticity of flexible structure in flow controls

柔性结构与空气动力耦合形成的系统呈现出丰富的非定常、非线性流动和结构动力学行为,对其气动弹性效应合理地控制和利用,能够大幅度提高飞机机翼、风力机叶片等结构的气动性能,并使其具有气动自适应能力.本文总结了近年来与气弹效应应用相关的研究进展及存在的问题,具体介绍了薄膜翼型的流动控制特性、柔性壁面减阻技术以及Sinha扰流装置的发展过程、主要成果以及未来发展趋势,着重对相关试验、流固耦合数值分析、Lagrangian拟序结构动力学等理论分析方法进行总结,展示了气弹效应在流动控制方面的巨大潜力和深远的学术意义,以便更多的研究人员开展该领域的研究工作.     

A Hierarchical Finite Element for Geometrically Non-linear Vibration of Thick Plates

A p-version, hierarchical finite element for moderately thick isotropic plates is derived and free vibrations are studied. The effects of the rotatory inertia, transverse shear and geometrical non-linearity, due to moderately large displacements, are taken into account. The time domain free vibration equations of motion are obtained by applying the principle of virtual work, and are mapped into the frequency domain by the harmonic balance method. The ensuing frequency domain equations are solved by a predictor–corrector method. The convergence properties of the element, the influence of the plate's thickness and of the width to length ratio on the backbone curves and on the non-linear mode shapes are investigated. The first and higher order modes are analysed and results are compared with published results.     

论工程地质模型——涵义、意义、建模与应用

本文重点论述了工程地质模型。阐明了工程地质模型的概念与涵义、阐述了建模的依据与方法；探讨了其应用的范畴与要点；尤其充分论述了工程地质模型的学科意义与深远影响。工程地质模型是条件研究的归宿, 是分析、计算与评论的基础。