水射流辅助破岩机理研究(1)——气泡空蚀

利用MSC.Dytram软件中的具有强度的多材料欧拉网格模拟了高速运动气泡对靶板(岩石)的冲击过程,研究了靶板破坏深度、破坏宽度与气泡初始压力、速度以及入射方向等的关系,试图指出这些关系的物理意义并分析射流破岩机理.

Polymers with aligned carbon nanotubes: Active composite materials

We review the current state of the polymer-carbon nanotube composites field. The article first covers key points in dispersion and stabilization of nanotubes in a polymer matrix, with particular attention paid to ultrasonic cavitation and shear mixing. We then focus on the emerging trends in nanocomposite actuators, in particular, photo-stimulated mechanical response. The magnitude and even the direction of this actuation critically depend on the degree of tube alignment in the matrix; in this context, we discuss the affine model predicting the upper bound of orientational order of nanotubes, induced by an imposed strain. We review how photo-actuation in nanocomposites depend on nanotube concentration, alignment and entanglement, and examine possible mechanisms that could lead to this effect. Finally, we discuss properties of pure carbon nanotube networks, in form of mats or fibers. These systems have no polymer matrix, yet demonstrate pronounced viscoelasticity and also the same photomechanical actuation as seen in polymer-based composites. © 2008 Elsevier Ltd. All rights reserved.

Strength of nanotubes, filaments, and nanowires from sonication-onduced scission

A model to describe the cavitation-induced breakage of nanofilaments during their sonication in solution is proposed. The model predicts a limiting length below which scission no longer occurs, and accurately describes experimental results for materials ranging from carbon nanotubes to protein fibrils. Sonication-induced breakage can now be used as a probe for the strength of nanostructures. © 2009 WILEY-VCH Verlag GmbH & Co. KGaA.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Investigation on mechanism of critical cavitating flow in liquid jet pumps under operating limits

The critical cavitating flow in liquid jet pumps under operating limits is investigated in this paper. Measurements on the axial pressure distribution along the wall of jet pumps indicate that two-phase critical flow occurs in the throat pipe under operating limits. The entrained flow rate and the distribution of the wall pressure upstream lowest pressure section does not change when the outlet pressure is lower than a critical value. A liquid-vapor mixing shockwave is also observed under operating limits. The wave front moves back and forth in low frequency around the position of the lowest pressure. With the measured axial wall pressures, the Mach number of the two-phase cavitating flow is calculated. It's found that the maximum Mach number is very close to I under operating limits. Further analysis infers a cross-section where Mach number approaches to I near the wave front. Thus, the liquid-vapor mixture velocity should reach the local sound velocity and resulting in the occurrence of operating limits.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

Study of the origins of toughness in amorphous metals

Amorphous metals that form fully glassy parts over a few millimeters in thickness are still relatively new materials. Their glassy structure gives them particularly high strengths, high yield strains, high hardness values, high resilience, and low damping losses, but this can also result in an extremely low tolerance to the presence of flaws in the material. Since this glassy structure lacks the ordered crystal structure, it also lacks the crystalline defect (dislocations) that provides the micromechanism of toughening and flaw insensitivity in conventional metals. Without a sufficient and reliable toughness that results in a large tolerance of damage in the material, metallic glasses will struggle to be adopted commercially. Here, we identify the origin of toughness in metallic glass as the competition between the intrinsic toughening mechanism of shear banding ahead of a crack and crack propagation by the cavitation of the liquid inside the shear bands. We present a detailed study over the first three chapters mainly focusing on the process of shear banding; its crucial role in giving rise to one of the most damage-tolerant materials known, its extreme sensitivity to the configurational state of a glass with moderate toughness, and how the configurational state can be changed with the addition of minor elements. The last chapter is a novel investigation into the cavitation barrier in glass-forming liquids, the competing process to shear banding. The combination of our results represents an increased understanding of the major influences on the fracture toughness of metallic glasses and thus provides a path for the improvement and development of tougher metallic glasses.

Generation of strong quasistatic magnetic fields in interactions of ultraintense and short laser pulses with overdense plasma targets

An analytical fluid model is proposed for the generation of strong quasistatic magnetic fields during normal incidence of a short ultraintense Gaussian laser pulse with a finite spot size on an overdense plasma. The steepening of the electron density profile in the originally homogeneous overdense plasma and the formation of electron cavitation as the electrons are pushed inward by the laser are included self-consistently. It is shown that the appearance of the cavitation plays an important role in the generation of quasistatic magnetic fields: the strong plasma inhomogeneities caused by the formation of the electron cavitation lead to the generation of a strong axial quasistatic magnetic field B-z. In the overdense regime, the generated quasistatic magnetic field increases with increasing laser intensity, while it decreases with increasing plasma density. It is also found that, in a moderately overdense plasma, highly intense laser pulses can generate magnetic fields similar to 100 MG and greater due to the transverse linear mode conversion process.

Simulation of shock-induced bubble collapse with application to vascular injury in shockwave lithotripsy

Shockwave lithotripsy is a noninvasive medical procedure wherein shockwaves are repeatedly focused at the location of kidney stones in order to pulverize them. Stone comminution is thought to be the product of two mechanisms: the propagation of stress waves within the stone and cavitation erosion. However, the latter mechanism has also been implicated in vascular injury. In the present work, shock-induced bubble collapse is studied in order to understand the role that it might play in inducing vascular injury. A high-order accurate, shock- and interface-capturing numerical scheme is developed to simulate the three-dimensional collapse of the bubble in both the free-field and inside a vessel phantom. The primary contributions of the numerical study are the characterization of the shock-bubble and shock-bubble-vessel interactions across a large parameter space that includes clinical shockwave lithotripsy pressure amplitudes, problem geometry and tissue viscoelasticity, and the subsequent correlation of these interactions to vascular injury. Specifically, measurements of the vessel wall pressures and displacements, as well as the finite strains in the fluid surrounding the bubble, are utilized with available experiments in tissue to evaluate damage potential. Estimates are made of the smallest injurious bubbles in the microvasculature during both the collapse and jetting phases of the bubble's life cycle. The present results suggest that bubbles larger than 1 μm in diameter could rupture blood vessels under clinical SWL conditions.

Simulação numérica da cavitação em turbomáquinas usando uma formulação Euler-Lagrange.

Turbomáquinas são máquinas operacionais que transferem energia mecânica entre um rotor e um fluido. Estas máquinas têm muitas aplicações industriais. Um dos componentes de uma turbomáquina responsável pela transferência da energia, ou receber a rotação do eixo e transformar em energia de fluido em caso de bomba ou transferir a energia do fluido para o eixo em caso de uma turbina, é o impelidor ou rotor. O fenómeno da cavitação envolve escoamento bifásico: o líquido a ser bombeado e as bolhas de vapor que são formadas durante o processo de bombeamento. O processo de formação dessas bolhas é complexo, mas ocorre principalmente devido a presença de regiões de pressões muito baixas. O colapso dessas bolhas pode muitas vezes levar a deterioração do material, dependendo da intensidade ou da velocidade de colapso das bolhas. O principal objetivo deste trabalho foi estudar o comportamento hidrodinâmico do escoamento nos canais do impelidor de uma turbomáquina do tipo radial usando recursos de fluidodinâmica computacional (CFD). Uma abordagem Euler-Lagrange acoplada com o modelo da equação de Langevin foi empregada para estimar a trajetória das bolhas. Resultados das simulações mostram as particularidades de um escoamento líquido-bolha de vapor passando em um canal de geometria curva, fornecendo assim informações que podem nos ajudar na prevenção da cavitação nessas máquinas.

EFFECT OF GRAIN-BOUNDARY MICROSTRUCTURE ON CAVITY NUCLEATION IN A NICKEL-BASED SUPERALLOY.

A study is presented of grain-boundary cavitation produced in Nimonic 80A by cold-deformation and stress-free annealing. The cavities were found to originate either from transverse cracking of carbide particles, or from decohesion of the particle-grain boundary interfaces. This decohesion could occur either during deformation, or during annealing. The cavities were invariably located at or close to the point of impingement of a matrix slip band on the grain boundary, but not all slip bands at a particular boundary were associated with cavitation. Quantitative evidence is presented showing that the mean number of dislocations associated with each slip band increases with macroscopic strain, but there is considerable variation between slip bands. This accounts for the differential ability of slip bands to result in cavity nucleation.

动网格技术在非稳态空化流计算中的应用

该文基于对钝锥头航行体的空化流动的实验观测, 利用商业软件FLUENT~[1]考核了分块网格模型对降低动网格方法误差影响、提高空化流计算精度的可行性和有效性. 通过与实验的对比表明, 采用分块网格模型能够使基于动网格方法的非稳态空化流计算很好的捕捉到空泡生长, 回射, 翻卷, 脱落, 溃灭这一过程, 并获得较为准确的航行体的运动轨迹和所受的水动力. 进而分析了流场涡结构与空泡演化的机理

Propulsion characteristics of wing-in-ground effect dual-foil propulsors

Propulsion characteristics of wing-in-ground effect propulsors were investigated using a comparative analysis of thrust and powering characteristics between wing-in-ground (WIG) effect thrusters and traditional screw propellers. WIG thrusters were found to have constant thrust production and efficiency, nearly independent of speed of advance, as contrary to screw propellers, whose optimum efficiency occurs at only one speed point. To produce the same amount of thrust as equivalent screw propellers, WIG thrusters have to work under heavily loaded operating conditions. WIG thrusters were also found to produce a relatively lower but nearly constant efficiency and thrust, independent of speed. Another distinguishing propulsion characteristic revealed for WIG thrusters is that they are capable of operating at much higher speeds, in a range of three to six times that of screw propellers of the same size. While the speed range of screw propellers is mainly limited by their geometric pitch, the speed range of WIG thrusters has no speed limit in ideal fluid. In reality, the speed range is only limited by viscous drag and cavitation, or compressibility, in water or air, respectively. This suggests a potential for WIG thrusters of higher speed application than screw propellers. An experimental investigation and validation of the propulsion system is warranted. Crown Copyright (C) 2010 Published by Elsevier Ltd. All rights reserved.

Design Study of a Heavily Loaded Ice Class Propeller Using an Advance Panel Method

A design and optimization procedure developed and used for a propeller installed on a twin-semitunnel-hull ship navigating in very shallow and icy water under heavy load conditions is presented. The base propeller for this vessel was first determined using classic design routines under open-water condition with existing model test data. In the optimization process, a panel method code (PROPELLA) was used to vary the pitch values and distributions and take into account the inflow wake distribution, tunnel gap, and cavitation effects. The optimized propeller was able to improve a ship speed of 0.02 knots higher than the desired speed and 0.06 knots higher than the classic B-series propeller. The analysis of the effect of inflow wake, hull tunnel, cavitation, and blade rake angle on propulsive performance is the focus of this paper.