加热球体入水空泡实验研究

基于高速摄像方法,开展了17—800°C范围内不同温度球体垂直入水实验研究.呈现了随球体温度变化所产生的复杂入水现象:在1.5 m/s入水速度条件下,实验所采用的室温球体不能产生空泡,当球体温度为300°C时空泡生成,增加到400°C空泡消失,继续提高温度至700°C空泡再次形成.根据传热学与流体动力学理论,分析了温度与速度变化对空泡形成的影响机理.结果表明随着温度的升高,球体与水之间的传热效率与传热方式发生改变,汽化生成的汽泡和蒸汽膜改变了周围流体流动的湍动性和球体表面的粗糙度、疏水性,这些变化均会影响空泡的形成;在1.5—3.8 m/s入水速度范围内,当球体具有较高温度时,能否形成入水空泡主要与球体的传热性能有关,速度的提高增强了球体与水的传热效率,使高速入水条件下的较低温度球体同低速入水条件下的较高温度球体入水现象相似,速度本身仅对生成空泡的形态有所影响.     

基于高速摄像实验的开放腔体圆柱壳入水空泡流动研究

基于高速摄像方法,针对入水空泡流动特征和机理,进行了开放腔体圆柱壳垂直入水实验研究.通过对实验现象的观测,发现开放腔体圆柱壳入水运动会形成波动流动和云化流动两种流动方式,结合影像数据,分别描述了两种流动状态下的空泡形态特征,并获得了空泡波动参数的变化规律;对比不同入水速度实验,分析了入水速度对入水空泡流动方式和流动参数的影响;依据流体力学基本理论,分析了入水空泡波动和云化现象的形成机理.结果表明:随入水速度增加,入水空泡依次呈现波动和云化两种流动状态,波动频率与入水速度无关,闭合发生时间随入水速度增加而减小,与Froude数呈线性关系;入水导致开放空腔内部气体涨缩,引起开放端压力场和速度场周期性扰动,空泡截面扩展程度出现差异,形成空泡波动现象;空泡闭合后尾部形成回射流,回射流触及空泡壁面引起壁面流动转捩,形成空泡云化现象.     

基于势流理论的回转体并联入水双空泡演化动力学研究

两回转体并联入水过程中,双空泡在空间上相互干扰,使得单个入水空泡呈现非对称特性.本文基于势流理论,在现有二维轴对称入水空泡计算模型的基础上,对空泡干涉区域的三维流动进行简化,将相对流动对空泡发展的约束简化为约束势,分析了两回转体轴线内外两侧的空泡形态变化;基于非线性假设,引入影响函数,给出了空泡在三维空间演化的计算模型,并分析了同步并联入水过程空泡的三维演化特性.结果表明：回转体入水过程流场速度势可以看作由一个随回转体运动的点源和位于空泡轴线处的线源叠加产生;在并联入水过程中,双空泡演化在空间呈镜面对称,空泡间的相互扰动可以通过引入有势壁面进行分析;并联入水空泡半径随极角的变化与空泡截面所处深度有关,在靠近闭合点附近的抑制演化区空泡截面半径随极角的增大而逐渐减小,远离闭合点处的抑制演化区空泡截面半径随极角的增大而增大,空泡加权半径规律相反.     

冲击加载下铝的剪切模量

分别用Steinberg-Cochran-Guinan (SCG)模型、修正的SCG模型和有限应变理论对材料的剪切模量做了数值计算，并与一维平面应变加载下铝的实验结果进行了比较.结果表明，修正的SCG模型与实验结果较为符合.在10—80GPa的压力范围下，剪切模量随冲击压力的增加而逐渐增大，这是由于压力的影响占主要地位，发生了加工硬化.在80—125GPa的压力范围下，剪切模量随冲击压力的增大快速减小，这是因为温度的影响比较严重，发生了高温软化现象.剪切模量最终在冲击压力为125GPa处趋于零，这是由于在该压力点冲击熔化发生，剪切强度消失. 关键词： 剪切模量 Steinberg-Cochran-Guinan模型 有限应变理论 铝     

液滴低速撞击润湿球面现象观测分析

采用高速摄像仪以10000 帧/s的拍摄速度对液滴低速撞击润湿球体表面过程进行了实验观测, 分析了液滴撞击后的反弹、局部反弹和铺展等现象, 考察了黏度对撞击过程的影响; 在此基础上, 定量讨论了液滴铺展特征参数随撞击速度、球体直径和黏度的变化规律. 观测发现: 黏度较大且撞击速度较低时, 撞击后可能出现反弹和局部反弹, 黏度较小时则不发生; 铺展面积随撞击速度的增大而增大; 黏度增大时, 铺展因子减小; 在球体直径为4–20 mm范围内, 随着球体直径的增加, 铺展因子呈上升趋势. 关键词： 液滴撞击 润湿球面 铺展 黏度     

开放空腔壳体入水扰动流场结构及空泡失稳特征

采用高速摄像实验和数值计算相结合的方法,对开放空腔壳体入水过程中空腔内气体涨缩对入水空泡的扰动机理和扰动作用下空泡局部失稳特征开展研究.基于实验观测结果,对比开放空腔和封闭空腔两种壳体入水空泡形态差异性,获得开放空腔壳体入水空泡波动特征,并根据能量守恒定律和能量转化关系,定性分析空腔自激扰动机理和扰动引起的空泡波动机理.基于数值计算结果,结合实验观测到的空泡阶段性局部失稳现象,定量分析空泡局部的结构性失稳特征和流动性失稳特征,并参照边界层理论和漩涡理论,揭示了空泡局部失稳机理.结果表明:入水冲击压缩空腔气体形成扰动源,对流场结构形成周期性扰动,导致空泡波动;撞水前空腔气体经冲击压缩密度升高,导致入水后气体首次膨胀阶段部分气体外泄,改变分离点流动,空泡局部结构失稳;空泡壁面流动具有类边界层流动特性,波动形态空泡形成多级回流现象,并逐级作用在空泡凹陷位置,使局部掺混区厚度增加,产生涡旋、转捩流动,空泡局部云化失稳;空泡逐级在波谷位置闭合、脱落,伴随大尺度漩涡生成,脱落过程形成拟序结构流场,漩涡导致脱落空泡迅速溃灭,但不会对附着空泡的流动产生影响.     

冲击波作用下铝的等效剪切模量

提出了冲击波作用下通过测量沿着准弹性卸载过程纵波波速与体波波速得到等效剪切模量的方法,应用该方法对LY12铝在20—70 GPa冲击压力下的等效剪切模量进行了测量.实验结果表明,等效剪切模量随应力线性减少直至反向屈服时为零,而初始卸载时的等效剪切模量可以由修正的Steinberg-Cochran-Guinan模型进行描述.应用得到的等效剪切模量进行了初步数值模拟,计算结果与实验测量的粒子速度剖面符合很好,重现了准弹性卸载到塑性卸载的光滑过渡. 关键词： 等效剪切模量 波速 准弹性卸载 冲击波 铝     

不同头型弹体模型入水现象的实验研究

弹体入水问题作为理论研究方面的复杂问题,涉及到气-液-固三者之间相互影响,对其很难建立出一个能够广泛应用的、能精准地描述入水现象的模型.本文通过实验手段,利用高速摄像技术,记录了不同头型弹体入水过程.通过不同视角拍摄弹体垂直入水和垂直落入加有表面活性剂的水体实验,研究弹体头型对反弹水花速度、水面附近水花演变以及空泡的演变等现象的影响规律.实验验结果表明:同一下落高度下,顶角越大的弹体,水下空泡越大,空泡断裂时其底端的位置越靠下,水压加速距离越长,反弹水花速度越快,水面现象越剧烈,半球形头型弹体的排水现象最轻微.水体表面张力的改变对头部顶角为180°弹体水花形态有影响,而对60°,90°,120°和半球型弹体水花形态无影响.     

石英玻璃球撞击刚性壁的破碎过程

利用高速枪对石英玻璃球撞击刚性靶板进行了实验研究,分析了不同速度下球体的破碎过程和失效模式。当冲击速度低于临界破坏速度时,石英玻璃球以略低于原速从靶板回弹;当超过临界破坏速度时,球体呈现“压缩破碎区-表面剥落区-剪切破坏区”的破坏结构;进一步提高碰撞速度,剪切破坏区的扩展导致球体碎裂为若干“月牙状”的碎块;更高撞击速度下,石英玻璃球发生坍塌式破碎,在远离撞击端处产生层裂现象。利用离散元软件对球体的撞击破坏过程进行了模拟研究,球体在高速碰撞下的破碎可以分为弹性压缩、整体破碎和二次撞击3个阶段。球体碎裂前Hertz接触理论可以较好描述其撞击力,而破碎后的撞击力由于碎裂卸载远小于理论值,且偏差随冲击速度逐渐增加。     

单液滴冲击超疏水壁面的压力特性研究

雨滴撞击索类结构表面可能会激发振动及表面积冰等问题,现有研究多关注于疏水壁面及液滴撞击壁面铺展、回缩特性,鲜见涉及雨滴冲击超疏水壁面压力特性的研究.为此,采用CLSVOF方法对单液滴冲击超疏水固壁面这一过程进行数值计算,分析了液滴速度、初始直径等因素对液滴冲击超疏水壁面的压力特性、液滴动态行为特性及液滴与壁面接触时间的影响.结果 表明:单液滴撞击超疏水壁面的过程中,接触瞬间在接触点附近产生局部高压区,而在液滴铺展过程中,对壁面几乎没有压力冲击;在回缩反弹阶段,壁面受到持续较长时间的压力波动,且压力波动区域不局限于初始接触点附近较小范围.撞击速度或液滴初始直径的增大使壁面受到的冲击更为剧烈,且初始速度对壁面受压的影响更为明显.一定范围内液滴初始直径的增大则会导致接触时间延长.     

单液滴正碰球面动态行为特性实验研究

在考虑空气阻力影响,确定液滴撞击球面速度的基础上,对较高韦伯数液滴撞击干燥球面动态行为过程进行了实验研究,分析了球面曲率与韦伯数对液滴撞击行为和铺展因子的影响,并与前人撞击平面结果进行了对比.实验表明,靠近撞击球面时,液滴降落速度出现明显波动;球面曲率对液滴撞击后行为影响明显,曲率较大时,液滴撞击后铺展液膜会超出球面直径并滑落,曲率较小时,液滴撞击后在球面上呈现明显的铺展、回缩、震荡、着附动态变化行为,此时最大铺展因子受曲率影响小,随曲率减小,逐渐趋向于撞击平面时的最大铺展因子;韦伯数对液膜铺展速率影响较小,但对液膜回缩时间影响明显,最大铺展因子随韦伯数增加逐渐增大,获得的关联式呈指数变化.     

Dynamic strength behavior of a Zr-based bulk metallic glass under shock loading

Dynamic strength behavior of Zr51Ti5Ni10Cu25Al9 bulk metallic glass(BMG) up to 66 GPa was investigated in a series of plate impact shock-release and shock-reload experiments.Particle velocity profiles measured at the sample/Li F window interface were used to estimate the shear stress,shear modulus,and yield stress in shocked BMG.Beyond confirming the previously reported strain-softening of shear stress during the shock loading process for BMGs,it is also shown that the softened Zr-BMG still has a high shear modulus and can support large yield stress when released or reloaded from the shocked state,and both the shear modulus and the yield stress appear as strain-hardening behaviors.The work provides a much clearer picture of the strength behavior of BMGs under shock loading,which is useful to comprehensively understand the plastic deformation mechanisms of BMGs.     

液滴碰撞液膜润湿壁面空气夹带数值分析

采用复合水平集-流体体积法并综合考虑传热及接触热阻的作用, 对液滴碰撞液膜润湿壁面空气夹带现象进行了数值分析. 揭示了夹带空气形成机理, 探索了夹带空气特性参数随碰撞速度和液膜厚度的变化规律, 获得了夹带空气作用下液滴碰撞润湿壁面的传热机理. 研究结果表明: 撞壁前气液两相压力差是引起气液相界面拓扑结构变化以及夹带空气形成的主要原因; 液滴碰撞速度与压缩空气层内压力以及相界面形变高度密切相关; 液滴接触液膜时, 碰撞轴上液滴底部和液膜表面速度相等, 大约是碰撞速度的1/2; 碰撞速度对夹带空气层底部到破碎点的无量纲弧长和最大无量纲夹带空气直径均存在较大的影响; 液滴和液膜的无量纲形变高度与斯托克斯数密切相关; 液膜初始厚度对液滴和液膜的无量纲形变高度和最大无量纲夹带空气直径影响较大; 撞壁初始阶段, 碰撞中心区域夹带空气对壁面热流密度分布存在较大的影响.     

Water entry of a superhydrophobic low-density sphere

Abstract  

This paper virtually demonstrates an air cavity formed by a superhydrophobic low-density sphere entering into water with the use of a high-speed camera. Unlike the previous results for a high-density sphere, the low-density sphere pinches off the resulting air cavity at the shorter depth before the cavity has fully grown, and attracts minimal air on the rear surface due to rapid deceleration of the sphere in water. The analytical model previously proposed on the time evolution of the cavity shape is thus shown to be no longer valid for the low-density sphere.     

Toughening through multilayering in TiN–AlTiN films

The damage response of columnar multilayers of TiN and AlTiN to Vickers indentation is studied through focused ion beam machining and elastic modelling. Multilayers display an enhanced resistance, which increases with layer refinement, to the multiple fracture modes that appear at high loads in these materials, including edge (nested) cracks and inclined shear cracks. Measurements of layer thickness reveal that multilayers display additional modes of plastic deformation that lead to permanent compression and bending of the film. An elastic model of contact deformation in a bilayer where plasticity is mimicked by greatly enhanced elastic compliance of the film is used to rationalize the trends in crack resistance. It is shown that the enhanced toughness is not due to any increase in the strain capacity (hardness/modulus) of the film material, brought about by multilayering.     

Computer simulation of the matrix-inclusion interphase in bulk metallic glass based nanocomposites

Atomistic models for matrix-inclusion systems are generated. Analyses of the systems show that interphase layers of finite thickness appear interlinking the surface of the nanocrystalline inclusion and the embedding amorphous matrix. In a first approximation, the interphase is characterized as an amorphous structure with a density slightly reduced compared to that of the matrix. This result holds for both monatomic hard sphere systems and a Cu(47.5)Zr(47.5)Al(5) alloy simulated by molecular dynamics (MD). The elastic shear and bulk modulus of the interphase are calculated by simulated deformation of the MD systems. Both moduli diminish with decreasing density but the shear modulus is more sensitive against density reduction by one order of magnitude. This result explains recent observations of shear band initiation at the amorphous-crystalline interface during plastic deformation.     

Material property estimates from ultrasound attenuation in fibre suspensions

An investigation of a new method for measuring fibre material properties from ultrasonic attenuation in a dilute suspension of synthetic fibres of uniform geometry is presented. The method is based on inversely solving an ultrasound scattering and absorption model of suspended fibres in water for the material properties of the fibres. Experimental results were obtained from three suspensions of nylon 66 fibres each with different fibre diameters. A forward solution to the model with reference material values is compared to experimental data to verify the model’s behaviour. Estimates of the shear and Young’s modulus, the compressional wave velocity, Poisson’s ratio and loss tangent from nylon 66 fibres are compared to data available from other sources. Experimental data confirms that the model successfully predicts that the resonance features in the frequency response of the attenuation are a function of diameter. Consistent estimated values for the compressional wave velocity and the Poisson’s ratio were found to be difficult to obtain but in combination gave values of shear modulus within previously reported values and with low sensitivity to noise. Young’s modulus was underestimated by 54% but was consistent and had low sensitivity to noise. The underestimation is believed to be caused by the assumption of isotropic material used in the model. Additional tests on isotropic fibre would confirm this. Further analysis of the model sensitivity and the reasons for the resonance features are required.     

Regional stretch method to measure the elastic and hyperelastic properties of soft materials

Characterizing the mechanical properties of soft materials and biological tissues is of great significance for understanding their deformation behaviors. In this paper, a regional stretching method is proposed to measure the elastic and hyperelastic properties of a soft material with an adhesive surface or with the aid of glue. Theoretical and dimensional analyses are performed to investigate the regional stretch problem for soft materials that obey the neo-Hookean model, the Mooney-Rivlin model, or the Arruda-Boyce model. Finite element simulations are made to determine the expressions of the dimensionless functions that correlate the stretch response with the constitutive parameters. Thereby, an inverse approach is established to determine the elastic and hyperelastic properties of the tested materials. The regional stretch method is also compared to the indentation technique. Finally, experiments are performed to demonstrate the effectiveness of the proposed method.     

Identification of material properties of orthotropic elastic cylinders immersed in fluid using vibroacoustic techniques

A numerical study is presented to show the potential for using vibroacoustic-based experiments to identify elastic material properties of orthotropic cylindrical vessels immersed in fluids. Sensitivity analyses and a simulated inverse problem are shown to quantify the potential for material characterization through the use of acoustic emissions. For comparison purposes, the analyses are also shown with the normal component of the velocity at the surface of the cylinder as the measured response in place of the acoustic pressure. The simulated experiment consisted of an orthotropic cylinder immersed in water with an impact force applied to the surface of the cylinder. The material parameters of the cylinder considered in the analyses were the circumferential and longitudinal elastic moduli, and the in-plane shear modulus. The velocity response is shown to provide sufficient information for characterizing all three moduli from a single experiment. Alternatively, the acoustic pressure response is shown to provide sufficient information for characterizing only the two elastic moduli from a single experiment. The analyses show that the acoustic pressure response does not have sufficient sensitivity to the in-plane shear modulus for characterization purposes.     

Cavity formation during water entry of heated spheres

We experimentally study the cavity formation when heated spheres impact onto water at low and high subcooling.The observations present that the formation and appearance of the cavity are affected by the boiling modes and the heat transfer intensity. In the nucleate-boiling regime, a rough cavity can be formed at a rather low impact velocity, while at the same velocity, the cavity formed in the film-boiling regime may have a very smooth interface with a stable vapor layer around the sphere. We discuss the effects of the impact speed, water and sphere temperatures on the stability of the vapor layer. For low subcooled water, the stable vapor layer will be disturbed when increasing the impact velocity, leading to a disturbed cavity. For high subcooled water, the film boiling has a particular boiling model in which the vapor layer around the sphere cannot keep its stability. In this particular film-boiling regime, no cavities can be formed at low impact velocities and only broken cavities can be formed at high impact velocities.