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.

Effect of plasma temperature on electrostatic shock generation and ion acceleration by laser

The effect of plasma temperature on electrostatic shock generated by a circularly polarized laser pulse in overdense plasma is studied by particle-in-cell simulation. Ion reflection and transmission in the collisionless electrostatic shock (CES) are investigated analytically. As the initial ion temperature is varied, a distinct transition from the laser-driven piston scenario with all ions being reflected to the CES scenario with partial ion reflection is found. The results show that at low but finite temperatures the ions are much more accelerated than if they were cold.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Micro-ramp control for oblique shock wave / boundary layer interactions

Supersonic engine intakes operating supercritically feature shock wave / boundary layer interactions (SBLIs), which are conventionally controlled using boundary layer bleed. The momentum loss of bleed flow causes high drag, compromising intake performance. Micro-ramp sub-boundary layer vortex generators (SBVGs) have been proposed as an alternative form of flow control for oblique SBLIs in order to reduce the bleed requirement. Experiments have been conducted at Mach 2.5 to characterise the flow details on such devices and investigate their ability to control the interaction between an oblique shock wave and the naturally grown turbulent boundary layer on the tunnel floor. Micro-ramps of four sizes with heights ranging from 25% to 75% of the uncontrolled boundary layer thickness were tested. The flow over all sizes of microramp was found to be similar, featuring streamwise counter-rotating vortices which entrain high momentum fluid, locally reducing the boundary layer displacement thickness. When installed ahead of the shock interaction it was found that the positioning of the micro-ramps is of limited importance. Micro-ramps did not eliminate flow separation. However, the previously two-dimensional separation was broken up into periodic three-dimensional separation zones. The interaction length was reduced and the pressure gradient across the interaction was increased.

Contrasting high-resolution characteristic shock-capturing methods in aeroacoustic test problems

Two shock-capturing methods are considered. One is based on a standard conservative Roe scheme with van Leer's MUSCL variable extrapolation method applied to characteristic variables and a Runge-Kutta time stepping scheme. The other is based on the novel CABARET space-time scheme, which uses two sets of staggered variables, one for the conservation step and the other for characteristic splitting into local Riemann invariants. The methods are compared in a range of 2-D inviscid compressible flow test cases. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Transient response of structures with uncertain properties to nonlinear shock loading

A method is presented to predict the transient response of a structure at the driving point following an impact or a shock loading. The displacement and the contact force are calculated solving the discrete convolution between the impulse response and the contact force itself, expressed in terms of a nonlinear Hertzian contact stiffness. Application of random point process theory allows the calculation of the impulse response function from knowledge of the modal density and the geometric characteristics of the structure only. The theory is applied to a wide range of structures and results are experimentally verified for the case of a rigid object hitting a beam, a plate, a thin and a thick cylinder and for the impact between two cylinders. The modal density of the flexural modes for a thick slender cylinder is derived analytically. Good agreement is found between experimental, simulated and published results, showing the reliability of the method for a wide range of situations including impacts and pyroshock applications. © 2013 Elsevier Ltd. All rights reserved.

The effect of wind tunnel size on incident shock boundary layer in-teraction experiments

We examine the effects of varying the tunnel width to height ratio on the shock boundary layer interac-tion of an incident oblique shock with a turbulent boundary layer. The computational domain is a simpli-fied representation of typical wind tunnel experiments; the top wall of the tunnel is not modeled; only the flow conditions imposed by the shock are modeled on the top of the computational domain. A hy-pothesis of the expected effect of width to height ratio is presented and tested computationally. All flows are found to be three dimensional for the single shock strength range of width to height ratios considered. The effect of tunnel width is a function of the boundary layer thickness which decreases the effective width.

Canonical normal shock wave/boundary-layer interaction flows relevant to external compression inlets

The normal shock wave/boundary-layer interaction is important to the operation and performance of a supersonic inlet, and the normal shock wave/boundary-layer interaction is particularly prominent in external compression inlets. To improve understanding of such interactions, it is helpful to make use of fundamental flows that capture the main elements of inlets, without resorting to the level of complexity and system integration associated with full-geometry inlets. In this paper, several fundamental flowfield configurations have been considered as possible test cases to represent the normal shock wave/boundary-layer interaction aspects found in typical external compression inlets, and it was found that the spillage diffuser more closely retains the basic flow features of an external compression inlet than the other configurations. In particular, this flowfield allows the normal shock Mach number as well as the amount and rate of subsonic diffusion to all be held approximately constant and independent of the application of flow control. In addition, a survey of several external compression inlets was conducted to quantify the flow and geometric parameters of the spillage diffuser relevant to actual inlets. The results indicated that such a flow may be especially relevant if the terminal Mach number is about 1.3 to 1.4, the confinement parameter is around 10%, and the width is around twice or three times the height. In addition, the area expansion downstream of the shock should be limited to the conservative side of incipient stall based on incompressible diffusers. Copyright © 2013 by the authors.

The first study on the effects of microcystin-RR on gene expression profiles of antioxidant enzymes and heat shock protein-70 in Synechocystis sp PCC6803

Microcystins are heptapeptide toxins produced by cyanobacteria. Microcystin-RR(MC-RR) is a common variant among the 80 variants identified so far. There have been many investigations documenting the toxic effects of microcystins on animals and higher plants, but little is known on the toxic effects of microcystins on algae, especially at molecular level. We studied the effects of MC-RR on gene expression profile of a few antioxidant enzymes and heat shock protein-70 (Hsp70) in Synechocystis sp. PCC6803. After two days post-exposure, a high dose toxin (5 mg/l, about 4.8 x 10(-3) mM) significantly increased expression levels of the genes gpx1, sodB, katG, acnB, gamma-TMTand dnaK2, while a relatively low dose toxin (1 mg/l, about 9.63 x 10(-4) mM) induced a moderate and slow increase of gene expression. Our results indicate that MC-RR could induce the oxidative stress in Synechocystis sp. PCC6803 and the increase in gene expression of antioxidant enzymes and Hsp70 might protect the organism from the oxidative damage. in addition, cell aggregation was observed during the early period of exposure, which might be a specific oxidative stress reaction to MC-RR. (C) 2008 Elsevier Ltd. All rights reserved.

Effects of pure microcystin-LR on the transcription of immune related genes and heat shock proteins in larval stage of zebrafish (Danio rerio)

Microcystins (MCs) are cyanobacterial toxins in water blooms that have received increasing attention as a public biohazard for human and animal health. Previous studies were mainly focused on the toxic effects on adult fish, rather than juvenile or larvae, and the response of fish immune system were usually neglected. This paper presents the first data of the effects of microcystin-LR (MC-LR) on transcription of several genes essential for early lymphoid development (Rag1, Rag2, Ikaros, GATA1, Lck and TCR alpha) and heat shock proteins (HSP90, HSP70, HSP60, HSP27) in zebrafish larvae. Relative changes of mRNA transcription were analyzed by real time PCR. The transcription of Rag1, Rag2, Ikaros, GATA1, Lck and TCR alpha were up-regulated when following exposure to 800 mu g/L MC-LR, which may indicate that specific lymphocytes differentiation and TCR/lg arrangement are induced to counteract the toxic effects of MC-LR. It was also interesting to note the dramatically increased transcription of HSP90. HSP70, HSP60 and HSP27, which may indicate their important roles as molecular chaperones under oxidative stress. (C) 2009 Elsevier B.V. All rights reserved.

Feeding rates affect heat shock protein levels in liver of larval white sturgeon (Acipenser transmontanus)

A two-week trial was conducted to study the effect of feeding rates on heat shock protein levels in larval white sturgeon. The larvae (30 day post hatch, 230 mg initial body weight) were fed a commercial feed (12.6% moisture, 49.5% crude protein. 20.7% Crude fat, and 8.6% ash) at 5, 15. or 25% body weight per clay (BW d(-1)). Liver heat shock proteins (Hsp) were measured before and after the larvae were subjected to a heat shock from 18 to 26 degrees C at 1 degrees C/15 min and maintained at 26 degrees C for 4 h thereafter. Before heat shock, larvae fed 5% BW d(-1) had significantly (P<0.05) lower final body weight, RNA/DNA ratio, whole body lipid and protein content, and Hsp60 and Hsp70 levels but higher protein efficiency ratio, and whole body moisture content than larvae fed the two higher feeding rates. Heat shock significantly induced Hsp60 and Hsp70 levels in the liver of all fish but they were lower in larvae fed the 5% than those fed 15 and 25% BW d(-1). Hsp70 level increased much more than Hsp60 after the heat shock Suggesting that Hsp70 is a more sensitive biomarker under our experimental conditions. (c) 2008 Elsevier B.V. All rights reserved.

Active control of normal shock wave/turbulent boundary layer interaction using "Smart" piezoelectric flap actuators

This paper looks at active control of the normal shock wave/turbulent boundary layer interaction (SBLI) using smart flap actuators. The actuators are manufactured by bonding piezoelectric material to an inert substrate to control the bleed/suction rate through a plenum chamber. The cavity provides communication of signals across the shock, allowing rapid thickening of the boundary layer approaching the shock, which splits into a series of weaker shocks forming a lambda shock foot, reducing wave drag. Active control allows optimum control of the interaction, as it would be capable of positioning the control region around the original shock position and control the rate of mass transfer. © 2004 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Differential expression of three Paralichthys olivaceus Hsp40 genes in responses to virus infection and heat shock

Heat shock proteins (Hsps) are a family of highly conserved cellular proteins present in all organisms, mediating a range of essential housekeeping and cytoprotective functions as well-known molecular chaperons and recently as regulators of the immune response. By subtractive suppression hybridization, three Hsp40 homologues have been identified in the flounder (Paralichthys olivaceus) embryonic cells (FEC) after treatment with UV-inactivated turbot (Scophthalmus maximus L.) rhabdovirus (SMRV), termed PoHsp40A4, PoHsp40B6 and PoHsp40B11, whose encoded proteins all possess the conserved DnaJ domain, a signature motif of the Hsp40 family. Based on different protein structure and phylogenetic analysis, they can be categorized into two subfamilies, PoHsp40A4 for Type I Hsp40, PoHsp40B6 and PoHsp40B11 for Type 11 Hsp40. Further expression analysis revealed two very different types of kinetics in response either to heat shock or to virus infection, with a marked induction for PoHsp4OA4 and a weak one for both PoHsp40B6 and PoHsp40B11. A very distinct tissue distribution of mRNA was also revealed among the three genes, even between PoHsp40B6 and PoHsp40B11. This is the first report on the transcriptional induction of Hsp40 in virally stimulated fish cells, and the differential expressions might reflect their different roles in unstressed and stressed cells. (c) 2005 Elsevier Ltd. All rights reserved.

All-fiber Yb-doped laser mode-locked by nanotubes

Single-wall carbon nanotubes (SWNTs) and graphene have emerged as promising saturable absorbers (SAs), due to their broad operation bandwidth and fast recovery times [1-3]. However, Yb-doped fiber lasers mode-locked using CNT and graphene SAs have generated relatively long pulses. All-fiber cavity designs are highly favored for their environmental robustness. Here, we demonstrate an all-fiber Yb-doped laser based on a SWNT saturable absorber, which allows generation of 8.7 ps-long pulses, externally compressed to 118 fs. To the best of our knowledge, these are the shortest pulses obtained with SWNT SAs from a Yb-doped fiber laser. © 2013 IEEE.

Self-starting passively mode-locking all-solid-state laser with GaAs absorber grown at low temperature

We realize a stable self-starting passively mode-locking all-solid-state laser by using novel GaAs mirrors as the absorber and output coupler. The GaAs mirror is grown by the technology of metal organic chemical vapour deposition at low temperature. With such an absorber as the output coupler in the laser resonator, laser pulses with duration of 42ps were generated at a repetition rate of 400MHz, corresponding to the average power of 590mW.