冲击载荷下脆性空心颗粒破碎机理

为考察脆性空心颗粒在冲击载荷作用下的应变率效应和破碎行为的细观机理,以粉煤灰漂珠为研究对象,基于低速冲击实验和有限元数值模拟,对比了典型空心颗粒材料在不同加载速率下的力学响应特性和细观压溃行为,阐释了材料宏观应变率效应产生的细观机理,获得以下结果。（1）在0.001～300 s?1应变率范围,漂珠颗粒的破碎率和Hardin破碎势平均提升了约21%和10%～30%,材料比吸能提升了50%～125%,比吸能的额外增加主要与动态颗粒滑移产生的摩擦耗能相关。颗粒平均尺寸较大的试样体现出更强的应变率效应。（2）初始压溃阶段的应力应变响应特征的数值模拟结果与实验结果较吻合,低速冲击下动态二次压溃现象产生的细观机理为动态颗粒滑移和压紧行为对加载速率的依赖性。(3) 数值模拟表明,冲击加载下产生相同应变时颗粒的损伤程度和范围大于准静态加载,这与实验所得破碎势随应变率增加的结果一致。对比低速冲击实验的相对破碎势分析和细观数值模拟结果可知,脆性颗粒堆积材料在动态冲击下表现出的宏观应变率效应主要归因于颗粒压溃行为的率敏感性和动态加载下颗粒破碎能量利用率的降低。     

颗粒破碎的微观力学模拟

基于颗粒尺度建立了考虑颗粒破碎过程中能量耗散的微观力学模型。为了模拟颗粒体的破碎过程和侧向应力的影响,把胶结在一起的三个等粒径的颗粒体当作颗粒破碎体。基于颗粒破碎准则和热力学原理,考虑了颗粒破碎过程中颗粒接触点之间的塑性滑移和颗粒转动以及侧向应力作用,通过颗粒破碎体内部颗粒之间的接触破碎来模拟颗粒体的破碎。计算结果表明:本文所提出的微观力学模型可以模拟胶结在一起的大颗粒的破碎过程;颗粒之间的微观力学参数对颗粒破碎过程中的力、塑性滑移、颗粒的转动角有不同程度的影响;力-位移关系曲线在破碎前是线性的,而颗粒破碎时竖向力迅速降低到较小值,并伴有颗粒的转动。     

微喷混合问题的初步研究

基于颗粒轨道模型,采用二元碰撞假设,在二维拉氏流体动力学程序上实现了气-粒两相混合程序的研制.使用气-粒两相混合程序对不同充气压力条件下的微喷混合问题进行了数值模拟,给出了喷射场随时间演化过程,计算结果与实验定性一致.     

岩石单颗粒压缩破碎试验研究

岩石颗粒结构复杂、形态不规则,很难用现有经典理论进行研究,分形理论可以从统计的角度研究颗粒破碎的强度。根据600颗大理岩单颗粒破碎试验,测量单颗粒破碎的应力与竖向应变关系曲线和破碎后的颗粒分布曲线,将颗粒破碎曲线按形态分为三种类型;根据颗粒破碎后的颗粒分布特性,建立了岩石颗粒破碎的分形模型,测得了大理岩颗粒破碎的分维D为2.48;根据颗粒破碎的分形模型,导出颗粒破碎强度的理论公式;将大理岩颗粒破碎强度的试验数据与理论公式预测结果进行比较,两者十分吻合,从而验证了颗粒破碎强度的分形理论。     

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.      

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.     

固体火箭冲压发动机补燃室沉积数值模拟

考虑液态颗粒碰撞和聚合过程、液态颗粒和壁面碰撞过程,建立了固体火箭冲压发动机补燃室沉积数值计算模型,对模型发动机补燃室内颗粒之间碰撞、沉积的相互作用过程进行了数值模拟,得到了颗粒沿轴向和沿出口平面径向方向的分布情况,并计算得到了补燃室壁面不同段的沉积层厚度值。将计算结果与试验结果进行对比,最大误差为0.8mm,表明该计算模型具有较高的计算精度。     

考虑破碎影响的颗粒材料亚塑性模型及应变局部化模拟

在基于2ndP-K应力率的亚塑性模型基础上,通过引入一个能够考虑颗粒破碎影响的孔隙比-平均压力临界状态方程,形成了一个能够模拟颗粒破碎影响的颗粒材料亚塑性模型,数值算例考查了颗粒破碎对应变局部化模式及位移-承载曲线的影响,结果表明,所建议模型具有模拟破碎对颗粒材料应变局部影响的良好性能。     

基于CV法的空间薄膜充气管的展开速度分析

为了有效控制充气薄膜结构模拟过程中的展开速度,本文首先从CV法(控制体积法)的基本理论入手,通过改变材料和结构的相关参数,对模型的设置进行了改进,以研究影响充气结构展开速度的因素;其次将两种不同充气速率对充气管展开速度影响的模拟结果与实验结果进行了对比,发现数值结果与实验结果趋势基本相同,从而验证了CV法的可行性;最后研究了材料参数、充气温度、排气面积等参数对充气管展开速度的影响。研究发现:充气速率越大,充气管的展开速度越大;材料的弹性模量影响较小;质量密度越大、温度越低、排气孔的面积越大,则充气管的展开速度越低。因此,从充气结构的材料本身和结构模型考虑,选用高密度材料、采用较低的充气温度、加大气囊的排气面积均能有效降低充气结构的展开速度,从而为充气结构的设计和应用提供参考。     

装置研制中爆炸力学的发展

介绍爆炸力学结合其他学科知识, 开展工程应用的理论和实验工作,进行爆炸装置的研究和设计.包括细观结构对材料力学性能的影响、金属自由面微喷以及多尺度数值模拟的应用.描述了作为应用基础学科的爆炸力学,在某些现代爆炸装置研制中多学科融合的爆炸力学将发挥重要作用.      

Simulation of the primary breakup of a high-viscosity liquid jet by a coaxial annular gas flow

The conversion of low-grade fossil and biogenic energy resources (petcoke, biomass) to a synthesis gas in a high pressure entrained flow gasification process opens a wide spectrum for high efficient energy conversion processes. The synthesis gas can be used for production of methane (SNG), liquid fuels (BtL, CtL) or as fuel for operation of a gas turbine in a combined cycle power plant (IGCC). The production of a tar free high quality syngas is a challenging objective especially due to the fact that typical liquid or suspension fuels for entrained flow gasifiers feature viscosities up to 1000 mPas. Fuel droplet conversion at typical entrained flow gasification conditions is characterized by heat up, evaporation and subsequent degradation of the vapour phase. To guarantee a high fuel conversion rate in the gasifier an efficient atomization of the fuel is required. Mainly twin-fluid burner nozzles are used for atomization of those typically high viscous fuels. The present study is focused on the assessment of the accuracy of CFD computations for the primary breakup of high-viscosity liquids using an external mixing twin fluid nozzle. In a first step experiments were performed with a Newtonian glycerol-water-mixture featuring a liquid viscosity of 400 mPas. Jet breakup was investigated using a high speed camera as well as PIV and LDA-System for a detailed investigation of the flow field. In a second step the experimental results serve as reference data to assess the accuracy of CFD computations. Compressible large eddy simulations (LES) were performed to capture the morphology of the primary breakup as well as the important flow field characteristics. A Volume of Fluid (VOF) approach was used to track the unsteady evolution and breakup of the liquid jet. Comparison of experimental and numerical results showed good agreement with respect to breakup frequency, velocity fields and morphology.     

Numerical investigation of turbulent-jet primary breakup using one-dimensional turbulence

Primary breakup to form droplets at liquid surfaces is an important fundamental process to study as it determines the initial properties of the dispersed phase, which affect mixing rates, secondary breakup, droplet collisions, and flow separation within the dispersed flow region. Primary breakup can be regarded as one of the least developed model components for simulating and predicting liquid jet breakup. However, it is of paramount importance in many technical applications, e.g. fuel injection in engines and spray painting. This paper presents a numerical investigation of primary breakup of a turbulent liquid jet in still air at standard conditions using the one-dimensional turbulence (ODT) modeling framework. ODT is a stochastic model that simulates turbulent flow evolution along a notional 1D line of sight by applying instantaneous maps to represent the effect of individual turbulent eddies on property profiles. An important feature of ODT is the resolution of all relevant scales, both temporal and spatial. The restriction to one spatial dimension in ODT permits affordable high resolution of interfacial and single-phase property gradients, which is key to capturing the local behavior of the breakup process and allows simulations at high Reynolds and Weber numbers that are currently not accessible to direct numerical simulations (DNS).This paper summarizes our extensions of the ODT model to simulate geometrically simple jet breakup problems, including representations of Rayleigh wave breakup, turbulent breakup, and shear-driven breakup. Each jet breakup simulation consists of a short temporal channel section to initialize a turbulent velocity profile at the nozzle exit followed by an adjacent jet section. The simulations are carried out for jet exit Reynolds number of 11,500, 23,000, 46,000 and 92,000 while the Weber number is varied within the range 10²–10⁷. We present results on breakup statistics including spatial locations of droplet release, droplet sizes and liquid core length. The results on primary breakup are compared to experimental results and models.     

Low-pressure twin-fluid atomization: Effect of mixing process on spray formation

The present work comparably examines four different twin-fluid atomizers operated under the same operating conditions. Spray formation was examined by several approaches. The internal flow pattern was estimated using a simplified analytical approach, and the results were supported by the observation of the liquid discharge in the near-nozzle region. A high-speed back illumination was used for visualisation of the primary breakup. In the region of fully developed spray, the dynamics of droplets was studied using a phase-Doppler analyser (PDA). The information obtained from all methods was then correlated. Results show that the spray formation process depends mainly on the internal design of twin-fluid atomizer at low gas to liquid ratios (GLR). The amount of gas influences the character of the internal two-phase flow, a mechanism of the liquid breakup, droplet dynamics and a resulting drop size distribution. Differences among the atomizers are reduced with the increase in GLR. Moreover, it was shown that a certain mixing process can inherently create the annular internal flow which generates a stable spray characterized by relatively low mean droplet size.     

Dynamics of two-phase downwards flows in submerged entry nozzles and its influence on the two-phase flow in the mold

Gas–liquid flows inside the submerged entry nozzle (SEN) of a slab mold and its influence in the flow field in the mold were studied using video recording, mathematical simulations and particle image velocimetry (PIV) approaches. Bubbly and annular flows in the SEN yield structurally-uncoupled and structurally coupled flows in mold, respectively. High gas loads (ratio of mass flux of gas and mass flux of liquid) at high casting rates lead to increases of bubble population and bubbles sizes due to coalescence processes whose rate exceeds that of their breakup. The presence of gas bubbles or gas layers inside the SEN lead to periodical twisting of the liquid flow that induces biased flows through both ports yielding uneven flows in the mold. A multiphase mathematical model predicts acceptably well the flow dynamics of two-phase flows inside the SEN.     

A unified spray model for engine spray simulation using dynamic mesh refinement

This study is based on dynamic mesh refinement and uses spray breakup models to simulate engine spray dynamics. It is known that the Lagrangian discrete particle technique for spray modeling is sensitive to gird resolution. An adequate spatial resolution in the spray region is necessary to account for the momentum and energy coupling between the gas and liquid phases. This study uses a dynamic mesh refinement algorithm that is adaptive to spray particles to increase the accuracy of spray modeling. On the other hand, the accurate prediction of the spray structure and drop vaporization requires accurate physical models to simulate fuel injection and spray breakup. The present primary jet breakup model predicts the initial breakup of the liquid jet due to the surface instability to generate droplets. A secondary breakup model is then responsible for further breakup of these droplets. The secondary breakup model considers the growth of the unstable waves that are formed on the droplet surface due to the aerodynamic force. The simulation results are compared with experimental data in gasoline spray structure and liquid penetration length. Validations are also performed by comparing the liquid length of a vaporizing diesel spray and its variations with different parameters including the orifice diameter, injection pressure, and ambient gas temperature and density. The model is also applied to simulate a direct-injection gasoline engine with a realistic geometry. The present spray model with dynamic mesh refinement algorithm is shown to predict the spray structure and liquid penetration accurately with reasonable computational cost.     

Shock wave interactions with particles and liquid fuel droplets

Shock waves traveling through a multiphase flow environment are studied numerically using the Flux Corrected Transport (FCT) algorithm. Both solid particles and liquid droplets are used as the dispersed phase with their trajectories being computed using a Lagrangian tracking scheme. The phases are coupled by including source terms which account for mass transfer, momentum, and energy exchange from the dispersed phase in the governing equations of motion for the gas phase. For solid particles, droplet size effects are examined at constant mass loading. Deceleration of the shock wave is observed with effects increasing with decreasing particle size. The equilibrium velocity attained is found to agree with analytical results for an equivalent dense gas with a modified specific heat ratio. For liquid droplets, a droplet breakup model is introduced and the results show a faster attenuation rate than with the solid particle model. The inclusion of vaporization to the breakup model is seen to increase the attenuation rate but does not alter the final equilibrium velocity. When an energy release model is used in the simulations, behavior resembling a detonation is observed under certain conditions, with energy release coupling with and accelerating the shock front. Received 17 July 2000 / Accepted 20 August 2002 / Published online 4 December 2002 Correspondence to: Dr. K. Kailasanath (e-mail: kailas@lcp.nrl.navy.mil)     

大速差射流燃烧室中烟煤与贫煤燃烧的数值模拟

本文基于颗粒相的轨道模型,对大速差射流燃烧室中烟煤粉与贫煤粉的二维流动,混合及燃烧进行了数值模拟,模拟结果从两相耦合的角度,阐明了煤粉颗粒在燃烧室中运动的规律,煤粉与大速差射流诱导的中心气体逆流之间的混合及其对煤粉火焰稳定的影响,指出此种燃烧室中煤粉火焰稳定的回流区燃烧机理,气相流场及回流区的预报结果与实验符合良好。     

Impact Failure of Planetary Materials:

The dynamic fragmentation of a fine grained granitoid material has been examined. Target thicknesses ranged from 7 to 40 mm and impact energies from 12 to 2,500 J. Combined particle image velocimetry and image enhancement techniques are introduced and have been used to measure the size and velocity of material ejected laterally from the rear of the target during impact testing. Non-dimensional groups were formed and fitted with coefficients to predict median values of the distribution of mass and kinetic energy among radial distance, R, from the impact centre, ejecta velocities, v, and ejecta lengths, L. The statistics are well correlated with increasing non-dimensional impact energy (positive correlation for radial distance and velocity, and negative correlation for ejecta length). Median values were used to collapse cumulative distributions and non-centred Gaussian fits were used to describe these curves. Approximately 85 % of the total mass and kinetic energy is captured between R/R_{50 %} = 0.3 to 2, v/v_{50 %} = 0.3 to 2, and L/L_{50 %} = 0.2 and 3. This data facilitates a better comparison among a wide range of test conditions, especially when attempting to extrapolate principal features of impacts into brittle materials at higher velocities. The ejecta tracking techniques and methodologies can be used to improve current impact testing experiments and computer modelling validation.     

爆炸驱动液体分层现象和液滴尺寸模拟

爆炸驱动液体介质外界面的分散和破碎是气溶胶云团形成的重要过程。采用基于维数分裂的欧拉程序和Youngs混合界面处理方法,对中心药爆炸驱动甘油和水介质流场的液体分层现象进行了数值模拟。结合试验结果推断提出了液滴形成过程的三种并存机制:外层射流破碎、内层R-T失稳和中间液层"空化"破碎,分别建立了不同液层破碎液滴的尺寸模拟方法。对比给出抛撒甘油和水装置初级液滴的尺寸分布及最外层理论射流量。     

Nonlinear instability of an annular liquid sheet exposed to gas flow

Nonlinear instability and breakup of an annular liquid sheet has been modeled in this paper. The liquid sheet is considered to move axially and is exposed to co-flowing inner and outer gas streams. Also, the effect of outer gas swirl on sheet breakup has been studied. In the developed model a perturbation expansion method has been used with the initial magnitude of the disturbance as the perturbation parameter. This is a comprehensive model in that other geometries of planar sheet and a coaxial jet can be obtained as limiting cases of very large inner radius and inner radius equal to zero, respectively. In this temporal analysis, the effect of liquid Weber number, initial disturbance amplitude, inner gas-to-liquid velocity ratio, outer gas-to-liquid velocity ratio and outer gas swirl strength on the breakup time is investigated. The model is validated by comparison with earlier analytical studies for the limiting case of a planar sheet as well as with experimental data of sheet breakup length available in literature. It is shown that the linear theory cannot predict breakup of an annular sheet and the developed nonlinear model is necessary to accurately determine the breakup length. In the limiting case of a coaxial jet, results show that gas swirl destabilizes the jet, makes helical modes dominant compared to the axisymmetric mode and decreases jet breakup length. These results contradict earlier linear analyses and agree with experimental observations. For an annular sheet, it is found that gas flow hastens the sheet breakup process and shorter breakup lengths are obtained by increasing the inner and the outer gas velocity. Axially moving inner gas stream is more effective in disintegrating the annular sheet compared to axially moving outer gas stream. When both gas streams are moving axially, the liquid sheet breakup is quicker compared to that with any one gas stream. In the absence of outer gas swirl, the axisymmetric mode is the dominant instability mode. However, when outer gas flow has a swirl component higher helical modes become dominant. With increasing outer gas swirl strength, the maximum disturbance growth rate increases and the most unstable circumferential wave number increases resulting in a highly asymmetric sheet breakup with shorter breakup lengths and thinner ligaments.