气流作用下同轴带电射流的不稳定性研究

通过对气体驱动同轴电流动聚焦的实验模型进行简化,开展了电场力和惯性力共同作用下同轴带电射流的不稳定性理论研究.在流动为无黏、不可压缩、无旋的假设下,建立了三层流体带电射流物理模型并得到了扰动在时间域内发展演化的解析形式色散关系,利用正则模方法求解色散方程发现了流动的不稳定模态,进而分析了主要控制参数对不稳定模态的影响.结果表明,在参考状态下轴对称模态的最不稳定增长率最大,因此轴对称扰动控制整个流场.外层气流速度越高,气体惯性力越大,射流的界面越容易失稳.内外层液-液同轴射流之间的速度差越大,射流越不稳定.表面张力对射流不稳定性起到促进作用.轴向电场对射流不稳定性具有双重影响:当加载电场强度较小时,射流不稳定性被抑制;当施加电压大于某一临界值时,轴向电场会促进射流失稳.临界电压的大小与界面上自由电荷密度和射流表面扰动发展关系密切.这些结果与已有的实验现象吻合,能够对实验的过程控制提供理论指导.     

轴对称射流场涡结构的离散涡段方法研究

本文用三维离散涡方法,模拟了轴对称圆射流涡结构的发展。     

同轴气流式液体射流分裂液滴粒径研究

针对同轴气流式液体射流分裂液滴粒径预测模型缺乏的现状,结合射流线性稳定性理论,建立了基于临界模数的同轴气流式黏性液体射流分裂液滴粒径表达式,在此基础上,分别研究了气流旋拧(气流同时存在轴向和周向运动)及流体物性(气体可压缩性、液体黏性、气液密度比和表面张力)对液滴粒径的影响规律.研究发现:周围气流轴向引射作用和同轴旋转...     

与温度相关的非均质圆筒的非线性耦合广义热弹性分析

基于能量守恒方程和L–S广义热弹性理论,借助状态空间技术和Newmark法求解了材料性质沿径向任意梯度分布同时又与温度相关的非均质圆筒非线性耦合广义热弹性问题。通过对材料性质与温度无关和相关功能梯度圆筒的算例分析,给出了在线性和非线性耦合下圆筒内温度和应力沿径向和随时间的变化关系,验证了本文解的正确性和有效性。数值结果表明,考虑材料性质是否与温度相关,能量守恒方程中耦合项是线性还是非线性,得到的温度与应力均存在不同程度的差异。本文解可方便地应用于不同边界条件和初始条件下圆筒的广义热弹性分析。

     

不可压黏性流问题的无网格法分析

不可压缩黏性流问题一般采用Navier-Stokes方程来描述，基于加权残值法，推导了问题的无网格伽辽金法（EFGM）离散Navier-Stokes方程，在时间域上采用分步方法计算，速度和压力由相互独立的方程以解耦的形式求解，并采用同阶移动最小二乘近似，在每一时间步中，对压力解和速度解采用了Newton-Raphson迭代法进行修正，最后将所得到的方法应用到剪切驱动空腔流问题中，验证了方法的有效性，且解的精度高、稳定性好。     

NUMERICAL SIMULATION OF UNSTEADY VISCOUS FLOWS USING AN IMPLICIT PROJECTION METHOD

In this paper an implicit projection method for the solution of the two-dimensional, time-dependent, incompressible Navier– Stokes equations is presented. The basic principle of this method is that the evaluation of the time evolution is split into intermediate steps. The computational method is based on the approximate factorization technique. The coupled approach is used to link the equations of motion and the turbulence model equations. The standard k-ϵ turbulence model is used. The current methodology, which has been tested extensively for steady problems, is now applied for the numerical simulation of unsteady flows. Several cases were tested, such as plane or axisymmetric channels, a backward-facing step, a square cavity and an axisymmetric stenosis.     

同轴旋转可压缩流动中液体射流稳定性

基于线性稳定性理论,建立了描述同轴旋转可压缩流动中超空化条件下液体射流稳定性的数学模型,并对数学模型及其求解方法进行了验证;在此基础上,对模型中考虑的射流及气体可压缩性、气体同轴旋转以及超空化等因素对射流稳定性的影响进行了分析. 分析结果表明,模型中考虑射流及气体的可压缩性后,与不考虑可压缩性相比,计算得到的射流稳定性明显变差,最小液滴直径减小,分裂液滴直径变化范围变宽,且小液滴数量增多. 气体的同轴旋转在轴对称与非轴对称扰动下对射流稳定性的影响完全相反;轴对称扰动时,气体旋转使射流稳定性增强,而非轴对称扰动时则正好相反;气体旋转有可能导致影响射流稳定性的扰动模式发生根本性变化. 超空化使射流稳定性变差;超空化程度较弱时,超空化使分裂液滴最小直径减小,分裂液滴直径变化范围增大;而超空化达到一定程度后,进一步提高超空化程度,分裂液滴最小直径几乎保持不变.     

圆环旋转黏性液体射流空间不稳定性研究

利用线性稳定性理论, 进行了液体黏性对不同旋转强度下圆环旋转液体射流 空间不稳定性影响的研究. 在推导出的三维扰动下具有固体涡核型旋转速度分布的圆环旋转 黏性液体射流色散方程的基础上, 针对中低速射流, 进行了类反对称模式与类对称模式下圆 环旋转黏性液体射流的空间不稳定性分析. 研究结果表明, 对于旋转强度较大的圆环旋转液 体射流, 液体黏性的增加, 不利于射流的破碎; 随着液体黏性的增加, 射流的特征频率和最 不稳定波数减小. 然而, 对于旋转强度较小的圆环旋转液体射流, 液体黏性的增加, 有利于 射流的破碎; 随着液体黏性的增加, 类反对称模式下射流特征频率先减小后增大, 类对称模 式下射流特征频率增大; 随着液体黏性的增加, 类反对称模式下射流最不稳定波数先减小后 增大, 类对称模式下射流最不稳定波数增大.     

TRANSIENT SIMULATION OF INTERNAL SEPARATED FLOWS USING AN INTELLIGENT HIGHER-ORDER SPATIAL DISCRETIZATION SCHEME

This paper summarizes the method-of-lines (MOL) solution of the Navier–Stokes equations for an impulsively started incompressible laminar flow in a circular pipe with a sudden expansion. An intelligent higher-order spatial discretization scheme, which chooses upwind or downwind discretization in a zone-of-dependence manner when flow reversal occurs, was developed for separated flows. Stability characteristics of a linear advective–diffusive equation were examined to depict the necessity of such a scheme in the case of flow reversals. The proposed code was applied to predict the time development of an impulsively started flow in a pipe with a sudden expansion. Predictions were found to show the expected trends for both unsteady and steady states. This paper demonstrates the ease with which the Navier–Stokes equations can be solved in an accurate manner using sophisticated numerical algorithms for the solution of ordinary differential equations (ODEs). Solutions of the Navier–Stokes equations in primitive variables formulation by using the MOL and intelligent higher-order spatial discretization scheme are not available to date. © 1997 by John Wiley & Sons, Ltd.     

静气动弹性计算方法研究

对基于结构网格的Euler方程及N-S方程求解器和基于非结构网格的Euler方程求解器,采用结构模态分析方法和柔度矩阵方法,对无人机大展弦比机翼在Ma=0.6,α=2?, 飞行高度20km的巡航状态下的静气动弹性特性进行了数值模拟. 验证了两种求解器对静气动弹性模拟的准确性. 同时,对模态分析方法和柔度矩阵方法进行了对比研究,发现柔度矩阵方法更适用于静气动弹性数值模拟. 另外,对应用物面法向偏转方法替代网格变形技术模拟静气动弹性进行了研究,计算表明物面法向偏转方法可以大大提高静气动弹性计算效率和克服机翼结构变形过大时动网格技术无法处理的不足.      

DYNAMIC STRESS ANALYSIS OF GAS TURBINE ROTOR AIRFOILS USING THERMOELASTIC TECHNIQUES

The SPATE system provided a means for experimental ther-moelastic analysis of high frequency modes of vibration of a complicated turbine impeller which did not facilitate analysis with conventional methods due to low stress response and high modal density.
Testing also demonstrated a unique experimental stress analysis technique which was applied to a blade and vane airfoil and exhibited excellent stress distribution correlation with an established experimental technique (brittle lacquer) and with analytical stress distributions predicted by NASTRAN.
Full field dynamic stress analysis of gas turbine compressor airfoils accomplished using thermoelastic techniques provided an efficient method of determining optimum strain gage locations for monitoring rig and engine vibratory stress.
Significant cost savings could be realized using this technique from: (1) reduction in the amount of instrumentation required to characterize airfoil dynamic behavior and (2) increased vibratory response mode coverage thereby preventing the re-running of engine test programs to measure vibratory stress response levels of unexpected resonances encountered during initial engine testing.
Care is required when using SPATE results for vibratory modes resulting in high bi-axial stress distributions since the system output indicated stress magnitudes based on the sum of the principal stresses which could be significantly higher than either of the individual principal stresses.     

非牛顿流体雾化力学机理研究

本文用稳定性分析和数值计算方法研究了粘弹性液体射流雾化的力学机理。结果表明,雾化液体的 Ohnesorge 数和 Deborah 数的增大均使波频 ω_i 降低,但前者使空间增长率 k_1上升,而后者的作用相反。本文发现,对于雾化液体应使其表面扰动波的波数大于某一确定波数,并尽量降低液体的 Deborah 数。     

AN ARTIFICIAL BOUNDARY CONDITION FOR TWO-DIMENSIONAL INCOMPRESSIBLE VISCOUS FLOWS USING THE METHOD OF LINES

We design an artificial boundary condition for the steady incompressible Navier–Stokes equations in streamfunction–vorticity formulation in a flat channel with slip boundary conditions on the wall. The new boundary condition is derived from the Oseen equations and the method of lines. A numerical experiment for the non-linear Navier–Stokes equations is presented. The artificial boundary condition is compared with Dirichlet and Neumann boundary conditions for the flow past a rectangular cylinder in a flat channel. The numerical results show that our boundary condition is more accurate.     

BOUNDARY TREATMENT AND AN EFFICIENT PRESSURE ALGORITHM FOR INTERNAL TURBULENT FLOWS USING THE PDF METHOD

The generalized Langevin model, which is used to model the motion of stochastic particles in the velocity–composition joint probability density function (PDF) method for reacting turbulent flows, has been extended to incorporate solid wall effects. Anisotropy of Reynolds stresses in the near-wall region has been addressed. Numerical experiments have been performed to demonstrate that the forces in the near-wall region of a turbulent flow cause the stochastic particles approachi ng a solid wall to reverse their direction of motion normal to the wall and thereby, leave the near-wall layer. This new boundary treatment has subsequently been implemented in a full-scale problem to prove its validity. The test problem considered here is that of an isothermal, non-reacting turbulent flow in a two-dimensional channel with plug inflow and a fixed back-pressure. An efficient pressure correction method, developed in the spirit of the PISO algorithm, has been implemented. The pressure correction strategy is easy to implement and is completely consistent with the time- marching scheme used for the solution of the Lagrangian momentum equations. The results show remarkable agreement with both k–ϵ and algebraic Reynolds stress model calculations for the primary velocity. The secondary flow velocity and the turbulent moments are in better agreement with the algebraic Reynolds stress model predictions than the k– ϵ predictions. © 1997 by John Wiley & Sons, Ltd.     

离轴受限圆射流的三维湍流流场分析

利用原始变量,交叉网格的方法,以及采用κ-ε湍流模型,计算了离轴受限射流三维湍流流场,计算结果与用LDV测量结果作了比较,结果是满意的。计算结果提供了这样一个复杂的三维流场的清晰图象,为研究具有离轴射流的预燃室的稳燃机理提供了很多有用的信息,对计算结果作了详细的分析和讨论后,指出这样一个流场的重要三维特征是在θ方向梳线有大角度的扭曲。     

含能气体射流与液体相互作用的实验研究

为了研究整装式液体发射药的燃烧稳定性的控制方法,设计了点火喷射模拟装置及4种多级渐扩型观察室,利用数字高速摄像系统,观察含能气体射流在液体模拟工质中的扩展过程,并对实验中出现的喷孔壅塞现象进行了分析. 结果表明: 射流在渐扩型结构中扩展稳定,喷射压力、喷孔直径和渐扩结构对射流扩展形态和气液掺混过程有显著影响,通过合理调整这些参数,可以实现对射流扩展过程的有效控制;喷孔壅塞时射流扩展形态非对称,影响气液掺混,不利于控制射流的稳定性.     

基于模态分解的轴对称超声速射流啸声产生位置数值分析

不完全膨胀超声速射流的势核中会产生准周期的激波栅格结构, 其与剪切层内拟序结构的相互作用会产生激波噪声. 啸声是主要向上游方向传播的、具有离散频率的高强度激波噪声, 其产生是受一种非线性的声反馈环机制驱动. 精确定位啸声的声源位置是定量理解啸声反馈环机制和发展准确的啸声预测模型的一个关键所在. 为了分析近场啸声, 本文采用高精度数值方法直接求解轴对称可压缩Navier-Stokes方程, 数值模拟了完全膨胀射流马赫数为1.10和1.15的圆形声速喷管欠膨胀超声速冷射流, 得到了A₁和A₂两种轴对称模态啸声. 通过傅里叶模态分解、本征模态分解和动态模态分解, 分析了射流时序压力场和速度场, 研究了啸声关联拟序流动结构的空间演化, 精确定位了轴对称模态啸声的声源位置. 研究表明: 啸声关联拟序流动结构存在饱和态区域, 啸声声波是在其饱和态区域产生并向外传播; 在本文所涉及的射流马赫数范围内, A₁和A₂两种轴对称模态啸声的有效声源位置分别是在第4和第3个激波栅格结构的尾缘.      

CARBON NANOTUBES VIA METHANE DECOMPOSITION ON AN ALUMINA SUPPORTED COBALT AEROGEL CATALYST

An alumina-supported cobalt aerogel catalyst prepared from a sol-gel and a supercritical drying method was used in the catalytic decomposition of methane.The physical-chemical properties of the catalyst were characterized and its activity for methane decomposition was investigated.The effects of calcination and reaction temperatures on the activity of the catalyst and the morphology of the carbon nanotubes produced were discussed.A COAl2O4 spinel structure formed in the calcined catalyst.The quantity of the nanotubes produced in the reaction increases with the amount of cobalt in the reduced catalyst.A higher reaction temperature leads to a higher reaction rate,though faster deactivation of the catalyst occurs with the change.The carbon nanotubes grown on the catalyst have smooth walls and uniform diameter distribution.     

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

RNG湍流模型在冲击射流数值计算中的应用

将3种版本的RNG,k-ε湍流模型应用于轴对称冲击射流场的数值模拟. 通过与实验结果的比较,考察了模型系数的改变和非线性附加项对数值结果的影响,评价了三种模型对冲击射流的数值预测能力.