阶梯型加速段对旋流喷嘴雾化特性的影响

旋流内芯是压力旋流式喷嘴最主要的旋流发生构件,其几何特征直接影响压力旋流式喷嘴的喷雾特性.目前采用平滑型加速段的旋流内芯导流效率较低.为减小高流量条件下的能量损失,使喷嘴旋流内芯加速段对喷雾介质产生预旋效应,增强旋流强度,本文设计喷嘴旋流内芯加速段为阶梯型,其下段阶梯相对上段阶梯旋转15°,旋向与喷嘴旋流槽方向相同.利用粒子动态分析仪(particle dynamics analysis system,PDA)和高速摄影(charge coupled device,CCD)系统实验研究了加速段结构改进前后喷嘴的喷雾流量、雾场索特尔平均直径(Sauter mean diameter,SMD)、雾滴速度以及喷雾锥角,并分析了SMD、雾滴速度的轴向和径向分布特性.结果表明,背压差0.08~0.46 MPa范围内,阶梯型加速段对喷雾介质具有较好的预旋效果.喷嘴的流量提高了48.0%~51.8%;喷雾的轴向速度提升了31.4%~32.8%,径向速度提升了1.6%~16.8%;喷雾锥角减小了4.21°~6.57°;较高背压差下喷雾下游的SMD减小了9.8%.与平滑型加速段相比,阶梯型加速段的设计有效地提高了喷嘴的雾化质量.     

应力诱发界面迁移下晶内孔洞的演化

随着微电子技术的迅猛发展, 集成电路中内连导线的失效问题引起广泛关注. 内连导线内部孔洞萌生、长大、漂移和失稳变形成狭长裂纹, 从而导致电路的开路失效. 这是内连导线失效的常见形式. 而界面迁移是导致微结构形态演化的主要机制之一. 本文基于界面迁移下微结构演化的经典理论和弱解描述, 建立了应力诱发界面迁移下微结构演化的有限单元法, 并验证了算法的可靠性. 对铜内连导线中晶内孔洞的演化进行了数值模拟, 详细分析了应力、线宽及形态比对晶内孔洞演化的影响. 研究结果表明, 椭圆形晶内孔洞存在生长和收缩两种演化分叉趋势. 通过大量数值分析得到了晶内孔洞演化的临界应力{\stackrel{\text{?}}{\sigma }}_c 、临界线宽{\stackrel{\text{?}}{h}}_c 和临界形态比{\beta }_c . 当\stackrel{\text{?}}{h}\mathrm{?}{\stackrel{\text{?}}{h}}_c, \stackrel{\text{?}}{\sigma }或\stackrel{\text{?}}{h} 时, 晶内孔洞会沿长轴长大; 反之, 晶内孔洞会收缩甚至愈合. 此外, 应力\beta 越大、线宽\stackrel{\text{?}}{\sigma } 越小或形态比\stackrel{\text{?}}{h} 越大, 晶内孔洞越易发生长大, 且孔洞面积增大速度越快; \beta 越小、\stackrel{\text{?}}{\sigma } 越大或\stackrel{\text{?}}{h} 越小, 晶内孔洞越易发生收缩, 且孔洞面积减小速度越快.      

扩展腔对方波型微混合器混合性能的影响研究

微混合器凭借节约试剂、混合强度高和易于集成等优点,在材料合成、医药制备和生化检测等领域中具有广泛的应用. 为了进一步提高混合性能,保证混合过程的安全性及生化反应结果的准确性,设计了一种带扩展腔的新型方波型微混合器. 在综合考虑混合强度和压降的前提下,通过实验研究和数值模拟分析了窄缝宽度、窄缝长度和扩展腔高度对微混合器混合性能的影响并得到了不同雷诺数\mathit{Re}条件下的最优结构参数. 与方波型微混合器的混合性能进行比较,发现\mathit{Re}=20时,带扩展腔的方波型微混合器的混合强度更高,其中\mathit{Re}时两者混合强度相差最多,可达12%. 在相同\mathit{Re}下,带扩展腔的方波型微混合器的压降要低于方波型微混合器. 对带扩展腔的方波型微混合器进行内部流场分析,发现扩展腔结构能在流体层流状态的基础上引入涡流,使通道中流体的流动状态发生改变、对流增强,进而混合性能提高.      

错列角度对双圆柱涡激振动影响的数值模拟研究

为研究错列角度\alpha对双圆柱涡激振动问题的影响,采用自主研发的基于CIP (constrained interpolation profile)方法的数值模型,对雷诺数\mathit{Re}=100、错列角度\alpha =0^{°}\mathrm{~}90^{°} (间隔15 #x00B0;)的等直径双圆柱涡激振动问题进行数值模拟. 模型在笛卡尔网格系统下建立,采用具有三阶精度的 CIP 方法求解 N-S (Navier--Stokes)方程,采用浸入边界法处理流--固耦合问题,避免了任意拉欧方法下的网格畸变和重叠动网格技术中的大量信息交换问题,保证了模型的计算效率. 重点分析不同错列角度\alpha上下游圆柱的升阻力系数、位移响应、涡脱频率和尾涡模态等. 结果表明:折合速度U_r=2.0\mathrm{~}3.0时,上下游圆柱升阻力随错列角度的增大基本呈单调增大的趋势;U_{\mathrm{r}}=5.0\mathrm{~}8.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“上凸”趋势,在\alpha =15^{°}\mathrm{~}30^{°}取得最大值;U_r=10.0\mathrm{~}13.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“下凹”趋势,在\alpha =30^{°}\mathrm{~}45^{°}取得最小值,且柱体横流向振幅和升力没有明显的对应关系. 最后,结合尾涡模态对以上规律的成因进行分析. 研究结果可为相关海洋工程设计提供参考.      

高雷诺数下多柱绕流特性研究

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性. 计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在范围内未出现类似于单圆柱绕流的阻力危机现象,其平均阻力系数\overline{C_d}、升力系数均方根{C\mathrm{\text{'}}}_l及斯特劳哈尔数\mathit{St}维持在一定范围内波动. 串列双圆柱绕流与串列双方柱绕流中,均选取L/D=2.0, 2.5, 3.0, 3.5和4.0这五中间距比进行计算. 串列双圆柱绕流中,当\mathit{Re}=2.2\times 10^4时,在内存在一临界间距比(L_c/D)使得L_c/D前后上下游圆柱的升阻力系数发生跳跃性变化,且当时,下游圆柱的阻力系数为负数. 而当\mathit{Re}=3.0\times 10^6时,则不存在临界间距比,且下游圆柱的阻力系数始终为正数. 串列双方柱绕流在\mathit{Re}=1.6\times 10^4和\mathit{Re}=1.0\times 10^6两种工况下的临界间距比分别处于和区间内,且当时,两个雷诺数下的下游方柱阻力系数均为负数.      

复杂加载条件下的砂土本构模型

试验表明,饱和砂土的应力应变关系具有显著的密度以及压力依存性,上述两点构成了描述砂土静力加载下变形特性无法忽视的因素. 此外,在循环加载等复杂加载作用下,砂土还会表现出明显的应力诱导各向异性以及相变转换特性. 基于在e--p空间中存在唯一的临界状态线这一基本假定,通过在e--p空间中引入当前状态点与临界状态线的距离R来作为反映密度与压力依存特性的状态参量, 将变相应力比以及峰值应力比表达为状态参量的指数函数,将上述应力比参量引入到统一硬化参量中可准确地反映初始状态下围压、密度 对于单调加载下应力应变关系的影响规律,能描述砂土剪缩、剪胀,应变软化、硬化等特性. 采用非相关联流动法则,p--q空间中采用水滴型屈服面,塑性势面为椭圆面,松砂在单调加载下的静态液化现象也可描述. 为反映循环加载下塑性体积应变的累积特性以及塑形偏应变的滞回特性,在循环加载下将状态参量R表达为应力比参量,并在硬化参数中引入描述应力诱导各向异性特性的旋转硬化部分,所提模型可有效地描述循环加载下剪切模量的衰减特性、刚度衰化性质、强度减小特性,在不排水约束作用下,则会产生往返活动性现象. 通过一系列的模型模拟与试验结果对比,验证了本构模型的有效性及适用性.      

功能梯度材料动力学问题的POD模型降阶分析

为了快速分析非均质材料结构在复杂载荷作用下的动态响应, 提出一种模型降阶方法, 只需计算结构在简单均质材料情况下的动力学问题, 进而用其计算结果对非均质材料结构进行分析. 首先, 采用结构内部任意一点处的材料参数值作为整个结构的材料参数, 利用有限元分析软件计算该均质材料结构在动态载荷作用下的位移场建立数据库, 该数据库包含计算模型各个节点(自由度为L)在某时间段内L个时刻的位移; 其次, 对数据库中的信息按照时间离散的特定方式组集成瞬像矩阵, 并利用特征正交分解方法对其进行分解, 得到该模型的H个特征正交基底, 选取其中能反应模型主要特征的个(其中\mathrm{~})作为一组最优基底, 通过这组基底建立模型的低阶离散控制方程; 最后, 求解低阶离散微分方程组, 得到功能梯度材料结构在复杂载荷作用下的位移场. 文中分别给出二维和三维算例, 比较了降阶模型和全阶模型计算结果, 验证了该方法的有效性, 并且计算效率能提高1{\mathrm{}}^{1)}2个数量级.      

PMMA膨胀环动态拉伸碎裂实验研究

在强动载作用下, 脆性材料的碎裂问题是一个重要的研究课题, 而脆性材料在冲击拉伸载荷下的力学行为的实验研究相对较匮乏. 提出了一种动态拉伸断(碎)裂的液压膨胀环实验技术, 可用于准脆性/脆性材料的动态拉伸. 利用该技术对有机玻璃(PMMA)圆环试件进行了不同膨胀速度下的动态碎裂实验研究. 从回收碎片的断口形貌和碎片内部残余裂纹观察可知试件的破碎由环向拉伸应力造成, 碎片断口处发出的稀疏波会将周围的拉伸应力卸载, 从而抑制其他裂纹的进一步发展. 利用超高速相机记录了试件的膨胀碎裂过程, 利用DISAR激光速度干涉仪获得了试件外表面粒子的径向膨胀速度历史, 通过试件上的应变片获得了试件的应变历史和断裂应变. 实验结果表明: 在拉伸应变率150\mathrm{~}500s^{-1}范围, 材料的动态断裂应变低于准静态加载下的断裂应变, 体现出“动脆”现象; 随着加载应变率的提高, PMMA 材料的碎片尺寸减小; 无量纲化的PMMA圆环的平均碎片尺寸介于韧性碎裂模型和脆性碎裂模型的预测数值之间, 反映出材料的准脆性特性.      

一种考虑液相惯性的一维准静力固结模型

经典Terzaghi一维固结理论不考虑孔隙流体惯性影响,且该理论在不同时期模型推导和表述结果差别较大,导致当前仍存在诸多困惑甚至认识混乱的现象. 在笔者前期研究大变形动力固结理论框架内,忽略固相惯性而重点考虑液相惯性影响,经过合理简化建立反映孔隙流体惯性的一维小变形固结波动模型. 该固结波模型具有频散和耗散特性. 采用分离变量法,可得到单面排水和瞬时加载条件下无量纲形式固结波解析解答. 算例分析结果表明:固结波发展规律受无量纲数D_{\mathrm{c}}变化影响而呈现不同性态;D_{\mathrm{c}}数值较大时固结波响应会出现阶跃和正负波动现象;当D_{\mathrm{c}}值较小时,可能出现Mandel-Cryer效应等特殊现象. 通过对早期和后期Terzaghi固结模型的分析和对比,初步探明Terzaghi固结理论模型内部的矛盾性,在普通土体坐标和固相体积坐标两种不同解读条件下,早期Terzaghi (1923,1925)模型可以分别诠释为具有小变形和大变形属性的不同固结模型. 在经典一维固结理论模型的不同诠释背景下,固结波模型也可以据此作出相应拓展和表述. 固结波理论揭示缩尺固结试验中土体物理力学参数与固结波响应两种因素之间存在一种不确定性矛盾,据此建议微观土力学研究重视尺度效应. 固结波模型的意义还在于,可为Terzaghi经典固结模型理论精度分析提供新的依据.      

超临界CO2磨料射流流场影响因素的模拟分析

为了揭示超临界CO₂磨料射流流场特性,利用计算流体动力学模拟软件,对超临界CO₂磨料射流结构及不同因素对射流流场的影响规律进行了研究。结果表明:超临界CO₂磨料射流轴向速度和冲击力随着喷距的增大,先增大后减小,即存在最优喷距,喷射压差为10~30 MPa时最优喷距为3~6倍喷嘴直径;喷射压差一定时,围压由10 MPa增至30 MPa对射流速度场及液相冲击力会造成较小的负面影响。通过超临界CO₂射流破岩实验对上述2因素进行了辅助对比验证;流体温度由333 K增至413 K,固液两相轴向速度增大,而流体密度降低,导致液相冲击力减弱;磨料浓度由3.0%连续增至11.0%,射流固液两相轴向速度逐渐降低,降幅逐渐减小。     

Controlled atomization using a twin-fluid swirl atomizer

This paper presents the results of an experimental study of a twin-fluid internally mixed swirl atomizer. In this type of injectors, atomization is attained by injecting a small amount of air (i.e. of the order of less than 16% of the mass flow rate of liquid) into a liquid stream within the injector and the two-phase air liquid mixture is passed through a swirling passage to impart a swirling motion to the flow. Since most of the energy for atomization is supplied to the liquid by the atomizing air, a significantly small pressure drop can produce very fine spray with a small amount of atomizing air. At low values of air–liquid mass ratio (ALR), the appreciable tangential component of velocity with respect to the axial velocity provides a hollow cone spray structure, which turns into a solid cone spray with the increase in axial momentum, through either an increase in ALR or the liquid supply pressure. The results presented in this paper suggest that the investigated injector could be used to control the flow rate and spray characteristics (e.g. spray cone angle, spray solidity, breaking distance, and the droplet diameter) independent of each other by simultaneously varying the supply pressure of the liquid and the atomizing air flow rate. The controlled atomization studied in this paper for a twin-fluid internally mixed swirl atomizer makes it attractive to be used for various commercial applications as the atomizer is capable of providing various spray characteristics depending upon the application requirement.     

An experimental investigation of swirl atomizer sprays

In our previous studies (Chu et al. in Heat Mass Transf 43(11):1213–1224, 2007), a theoretical model of swirl atomizers was successfully established. From the analysis, the equations for the droplet size, velocity components, the boundary layer thickness and the spray cone angle were deduced based on the fundamental governing equations. The purpose of this study is to further compare the experimental result with the theoretical one already gained by a satisfactory embodiment of series of experiments. The aim is to corroborate the analytical results of the influence of atomizer construction and controlled pressure difference on typical swirl chambers. The results provide the droplet diameter as a function of pressure difference, swirl atomizer geometry, flow rate, spray cone angle. The experimental outputs also show a good confirmation of theoretical results and can also be used for further optimization on existing swirl chambers. Based on the results obtained, an optimization methodology on characteristics of swirl atomizers is proposed with the adjustment of individual design parameter and the matching flow number.     

Imaging and PDA analysis of a GDI spray in the near-nozzle region

This paper describes a combined LDA, PDA and imaging analysis of the pressure swirl spray in the near-nozzle region of a GDI injector. This innovative approach in the use of multiple, complementary diagnostics facilitates the interpretation of a complex spray flow field.The LDA and PDA data were ensemble-averaged into time bins to produce comprehensive time-history and spatial profiles of liquid velocity, droplet velocity and size and the sample number. They indicated times at which the spray exhibited seven different characteristics. These were identified as: (a) pre-swirl spray, (b) spray cone develops, (c) spray cone relaxes, (d) maximum velocity in spray cone, (e) fully developed steady state, (f) spray cone collapses and (g) the spray detaches from the nozzle. The most effective method to present the spatial and temporal development of the spray was to superimpose the velocity vector and drop size field plots onto the spray images.This article is part of the special issue 11^th International Symposium on Applications of Laser Techniques to Fluid Mechanics, Lisboa, Portugal, July 2002, January 2004, Vol. 36, Issue no. 1     

Theoretical analysis of heat and mass transfer in swirl atomizers

In this study, theoretical analyses have been performed to investigate the effects of atomizer construction and controlled pressure difference of swirl atomizers. The analysis of fluid field in the swirl chamber is governed by mass/energy conservation rules; in the region outside the nozzle, the analysis of oscillation of liquid sheet is based on Squire’s expression for the amplitude growth rate. With some physical assumptions of control volume, initial values and model correlation, analytical results make it possible to predict film thickness, velocity distribution, spray cone angle and droplet size directly. The distribution of velocity profile and boundary layer thickness in the swirl chamber have been established with the aid of MATLAB. Based on the results we obtained, we here propose the change of individual design parameter and its corresponding flow number to optimize the performance of swirl atomizers.     

The spray characteristics of a pressure-swirl injector with various exit plane tilts

The gasoline spray characteristics of a pressure-swirl injector were investigated with various exit plane tilts. The analysis focused on the correlation between tilt angle and flow angle. Mie-scattering technique and phase Doppler anemometry were employed to analyze the macroscopic spray development and droplet size distribution of the spray. An analytical method for mass flux estimation was applied to understand the velocity distribution at the nozzle exit. The results showed that the spray shape and velocity distribution of the spray were more asymmetrical at high tilt angles. In particular, an opened hollow cone spray was formed when the tilt angle is greater than the complementary flow angle. The pressure drop inside the spray, one of the crucial factors for the swirl spray collapse at various surrounding conditions, was attenuated in this opened hollow cone spray since the pressure inside the spray was assimilated to the surrounding air pressure. The spray collapse at high fuel temperature and back pressure conditions did not appear when the tilt angle is larger than the complementary flow angle due to the reduced pressure drop inside the spray. However, tilt angle should be optimized to fulfill the requirements of spray robustness and avoid the locally rich area. The droplet size of 70° tilted nozzle spray shows a value similar to that of the original swirl spray in the plane that includes nozzle axis and the major axis of exit surface ellipse (Major axis plane) while it shows an increased value in the plane that includes nozzle axis and the minor axis of exit surface ellipse (Minor axis plane).     

圆环旋转黏性液体射流破碎液滴粒径与速度数量密度分布相关性研究

压力旋流喷嘴被广泛应用于航空发动机、船用发动机、车用汽油缸内直喷发动机、燃气轮机等动力机械的燃油喷射系统中.以压力旋流喷嘴射流为研究对象,开展了圆环旋转黏性液体射流破碎液滴粒径与速度数量密度分布相关性问题研究.对于液体射流,以往的研究往往对破碎液滴粒径数量密度分布或速度数量密度分布进行单独研究,对于这两种数量密度分布之间关系的研究较少;从相关性的角度对圆环旋转黏性液体射流破碎液滴粒径与速度数量密度分布之间的关系进行研究.采用最大熵原理方法建立了圆环旋转黏性液体射流破碎液滴粒径与速度联合概率密度函数.对圆环旋转黏性液体射流破碎液滴粒径与速度联合概率密度函数进行了讨论,对圆环旋转黏性液体射流破碎液滴粒径数量密度分布与速度数量密度分布的相关性问题进行了研究.研究结果表明,为了给出正确的圆环旋转黏性液体射流破碎液滴粒径与速度联合概率密度函数,射流守恒约束条件中必须同时包括质量守恒定律、动量守恒定律以及能量守恒定律;破碎液滴粒径的数量密度分布与速度数量密度分布密切相关;射流旋转强度对破碎液滴粒径数量密度与速度数量密度分布结构影响不大,对破碎液滴粒径数量密度和速度数量密度的分布区域影响较大.      

双柱单锥型液-液旋流管内流场的激光诊断

应用激光测速仪，对一种双柱单锥型液 液旋流管内的流动结构，进行了全场范围内的多工况流动诊断研究．揭示出其切向速度由内旋流区和外旋流区构成，其中内旋流区中的速度分布符合准强制涡关系，外旋流区中的速度分布符合准自由涡关系；轴向速度由上行流动区和下行流动区构成，两者之间在直管段以零速点作分界，在锥体段则以零速区作过渡并伴随有一定的回流出现，且该过渡区或回流区的大小随锥体截面的收缩而减小，直到进入直管段后消失；各湍流量的分布以管芯处最大向外逐渐减小，沿轴向是直管段中的湍流度大于锥体段中的湍流度，而且湍流度在旋流管内具有各向异性的特性．     

Flow and atomization characteristics of a twin-fluid nozzle with internal swirling and self-priming effects

A twin-fluid nozzle was proposed for low-pressure atomization. The nozzle is featured by swirling air flows in the mixing chamber. Liquid medium is thereby inhaled due to the pressure difference. An experimental work was performed to investigate the atomization performance of the nozzle and the hydrogen peroxide solution served as the liquid medium. Droplet size and droplet velocity were measured. Effects of the diameter of the air-injection orifice and the air-injection pressure were investigated. The results show that small droplet size is achieved with the proposed nozzle. As the spray develops, Sauter mean diameter (SMD) of the droplets decreases first and then increases, irrespective of the variation of the air-injection orifice diameter and the air-injection pressure. Overall SMD varies inversely with the air-injection orifice diameter and air-injection pressure. Near the nozzle, cross-sectional velocity distribution exhibits a peak-valley pattern, which is replaced with uniformized velocity distributions away from the nozzle. Similarity of cross-sectional radial velocity distribution at different air pressures is evidenced. Furthermore, the correlation between droplet size and droplet velocity is established.     

Regimes of spray formation in gas-centered swirl coaxial atomizers

Spray formation in ambient atmosphere from gas-centered swirl coaxial atomizers is described by carrying out experiments in a spray test facility. The atomizer discharges a circular air jet and an axisymmetric swirling water sheet from its coaxially arranged inner and outer orifices. A high-speed digital imaging system along with a backlight illumination arrangement is employed to record the details of liquid sheet breakup and spray development. Spray regimes exhibiting different sheet breakup mechanisms are identified and their characteristic features presented. The identified spray regimes are wave-assisted sheet breakup, perforated sheet breakup, segmented sheet breakup, and pulsation spray regime. In the regime of wave-assisted sheet breakup, the sheet breakup shows features similar to the breakup of two-dimensional planar air-blasted liquid sheets. At high air-to-liquid momentum ratios, the interaction process between the axisymmetric swirling liquid sheet and the circular air jet develops spray processes which are more specific to the atomizer studied here. The spray exhibits a periodic ejection of liquid masses whose features are dominantly controlled by the central air jet.     

Viscous flow through the swirl chamber of a pressure-swirl atomizer

We study theoretically the viscous flow through the swirl chamber of a pressure-swirl atomizer. The aim is to predict quantitatively the flow rate of a given liquid driven by a given pressure difference across the atomizer in order to explain the counter-intuitive experimental finding that, for moderate liquid dynamic viscosities, at a given driving pressure difference and atomizer geometry, a higher liquid viscosity results in a higher flow rate. The concept for the flow analysis is to subdivide the flow field in the swirl chamber into zones allowing for neglect of velocity components or boundary-layer simplification. The result is a quantitative prediction of the liquid flow rate for a given driving pressure difference and atomizer geometry, and for given liquid properties relevant for the discharge from the atomizer. Flow rates are compared to experimental data from various sources and show good agreement. Another part of the results is the diameter of the air core formed around the symmetry axis of the swirl chamber, which is of sub-millimetric order here. This result is compared to different experimental correlations and also shows very good agreement. For small values of the swirl velocity and/or very high liquid viscosities, the air core breaks down. The phenomenon of air core break down has been analysed, and it is shown that for high viscosities the air core breaks down due to weak swirl velocity.