Wallner lines, crack velocity and mechanisms of crack nucleation and growth in a brittle bulk metallic glass

Mode I fracture experiments were conducted on brittle bulk metallic glass (BMG) samples and the fracture surface features were analyzed in detail to understand the underlying physical processes. Wollner lines, which result from the interaction between the propagating crack front and shear waves emanating from a secondary source, were observed on the fracture surface and geometric analysis of them indicates that the maximum crack velocity is similar to 800 m s(-1), which corresponds to similar to 0.32 times the shear wave speed. Fractography reveals that the sharp crack nucleation at the notch tip occurs at the mid-section of the specimens with the observation of flat and half-penny-shaped cracks. On this basis, we conclude that the crack initiation in brittle BMGs is stress-controlled and occurs through hydrostatic stress-assisted cavity nucleation ahead of the notch tip. High magnification scanning electron and atomic force microscopies of the dynamic crack growth regions reveal highly organized, nanoscale periodic patterns with a spacing of similar to 79 nm. Juxtaposition of the crack velocity with this spacing suggests that the crack takes similar to 10(-10) s for peak-to-peak propagation. This, and the estimated adiabatic temperature rise ahead of the propagating crack tip that suggests local softening, is utilized to critically discuss possible causes for the nanocorrugation formation. Taylor's fluid meniscus instability is unequivocally ruled out. Then, two other possible mechanisms, viz. (a) crack tip blunting and resharpening through nanovoid nucleation and growth ahead of the crack tip and eventual coalescence, and (b) dynamic oscillation of the crack in a thin slab of softened zone ahead of the crack-tip, are critically discussed. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Wave propagation in a 2D nonlinear structural acoustic waveguide using asymptotic expansions of wavenumbers

Nonlinear acoustic wave propagation in an infinite rectangular waveguide is investigated. The upper boundary of this waveguide is a nonlinear elastic plate, whereas the lower boundary is rigid. The fluid is assumed to be inviscid with zero mean flow. The focus is restricted to non-planar modes having finite amplitudes. The approximate solution to the acoustic velocity potential of an amplitude modulated pulse is found using the method of multiple scales (MMS) involving both space and time. The calculations are presented up to the third order of the small parameter. It is found that at some frequencies the amplitude modulation is governed by the Nonlinear Schrodinger equation (NLSE). The first objective here is to study the nonlinear term in the NLSE. The sign of the nonlinear term in the NLSE plays a role in determining the stability of the amplitude modulation. Secondly, at other frequencies, the primary pulse interacts with its higher harmonics, as do two or more primary pulses with their resultant higher harmonics. This happens when the phase speeds of the waves match and the objective is to identify the frequencies of such interactions. For both the objectives, asymptotic coupled wavenumber expansions for the linear dispersion relation are required for an intermediate fluid loading. The novelty of this work lies in obtaining the asymptotic expansions and using them for predicting the sign change of the nonlinear term at various frequencies. It is found that when the coupled wavenumbers approach the uncoupled pressure-release wavenumbers, the amplitude modulation is stable. On the other hand, near the rigid-duct wavenumbers, the amplitude modulation is unstable. Also, as a further contribution, these wavenumber expansions are used to identify the frequencies of the higher harmonic interactions. And lastly, the solution for the amplitude modulation derived through the MMS is validated using these asymptotic expansions. (C) 2015 Elsevier Ltd. All rights reserved.

Weakly nonlinear acoustic wave propagation in a nonlinear orthotropic circular cylindrical waveguide

Nonlinear acoustic wave propagation is considered in an infinite orthotropic thin circular cylindrical waveguide. The modes are non-planar having small but finite amplitude. The fluid is assumed to be ideal and inviscid with no mean flow. The cylindrical waveguide is modeled using the Donnell's nonlinear theory for thin cylindrical shells. The approximate solutions for the acoustic velocity potential are found using the method of multiple scales (MMS) in space and time. The calculations are presented up to the third order of the small parameter. It is found that at some frequencies the amplitude modulation is governed by the Nonlinear Schrodinger Equation (NLSE). The first objective is to study the nonlinear term in the NLSE, as the sign of the nonlinear term determines the stability of the amplitude modulation. On the other hand, at other specific frequencies, interactions occur between the primary wave and its higher harmonics. Here, the objective is to identify the frequencies of the higher harmonic interactions. Lastly, the linear terms in the NLSE obtained using the MMS calculations are validated. All three objectives are met using an asymptotic analysis of the dispersion equation. (C) 2015 Acoustical Society of America.

Microcracks Propagation in One Al/SiCp under Impact Loading

Numerous microcracks propagation in one metal matrix composite, Al/SiCp under impact loading was investigated. The test data was got with a specially designed impact experimental approach. The analysis to the density, nucleating locations and distributions of the microcracks as well as microstructure effects of the original composite was received particular emphasis. The types of microcracks or debonding nucleated in the tested composite were dependent on the stress level and its duration. Distributions of the microcracks were depended on that of microstructures of the tested composite while total number of microcracks in unit area and unit duration, was controlled by the stress levels. Also, why the velocity was much lower than theoretical estimations for elastic solids and why the microcracks propagating velocities increased with the stress levels' increasing in current experiments were analysed and explained.

Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media

In this paper, we examine the characteristics of elastic wave propagation in viscoelastic porous media, which contain simultaneously both the Biot-flow and the squirt-flow mechanisms (BISQ). The frequency-domain Green's functions for viscoelastic BISQ media are then derived based on the classic potential function methods. Our numerical results show that S-waves are only affected by viscoelasticity, but not by squirt-flows. However, the phase velocity and attenuation of fast P-waves are seriously influenced by both viscoelasticity and squirt-flows; and there exist two peaks in the attenuation-frequency variations of fast P-waves. In the low-frequency range, the squirt-flow characteristic length, not viscoelasticity, affects the phase velocity of slow P-waves, whereas it is opposite in the high-frequency range. As to the contribution of potential functions of two types of compressional waves to the Green's function, the squirt-flow length has a small effect, and the effects of viscoelastic parameter are mainly in the higher frequency range. Crown Copyright (C) 2006 Published by Elsevier Ltd. All rights reserved.

Microcracks Propagation in One Al/SiCp under Impact Loading

Numerous microcracks propagation in one metal matrix composite, Al/SiCp under impact loading was investigated. The test data was got with a specially designed impact experimental approach. The analysis to the density, nucleating locations and distributions of the microcracks as well as microstructure effects of the original composite was received particular emphasis. The types of microcracks or debonding nucleated in the tested composite were dependent on the stress level and its duration. Distributions of the microcracks were depended on that of microstructures of the tested composite while total number of microcracks in unit area and unit duration, was controlled by the stress levels. Also, why the velocity was much lower than theoretical estimations for elastic solids and why the microcracks propagating velocities increased with the stress levels' increasing in current experiments were analysed and explained.

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.

Group velocity of probe light pulse in an open lambda system

The group velocity of the probe light pulse (GVPLP) propagating through an open Lambda-type atomic system with a spontaneously generated coherence is investigated when the weak probe and strong driving light fields have different frequencies. It is found that adjusting the detuning or Rabi frequency of the probe light field can realize switching of the GVPLP from subluminal to superluminal. Changing the relative phase between the probe and driving light. elds or atomic exit and injection rates can lead to GVPLP varying in a wider range, but cannot induce transformation of the property of the GVPLP. The absolute value of the GVPLP always increases with Rabi frequency of the driving light field increasing. For subluminal and superluminal propagation, the system always exhibits the probe absorption, and GVPLP is mainly determined by the slope of the steep dispersion.

Group velocity of the light pulse in an open V-type system

We investigate the group velocity of the probe light pulse in an open V-type system with spontaneously generated coherence. We find that, not only varying the relative phase between the probe and driving pulses can but varying the atomic exit rate or incoherent pumping rate also can manipulate dramatically the group velocity, even make the pulse propagation switching from subluminal to superluminal; the subliminal propagation can be companied with gain or absorption, but the superluminal propagation is always companied with absorption. (c) 2006 Elsevier GmbH. All rights reserved.

Upper level and phase controlling group velocity in V-type system

With the external field coupling the two upper levels, we investigate the light pulse propagation properties with weak probe field in a V-type system. Due to the external upper level (UL) coupling field, the dispersion of the system has been influenced by the relative phase. It is shown that the UL field and the relative phase can be regarded as switches to manipulate light propagation between subluminal and superluminal. (c) 2007 Elsevier B.V. All rights reserved.

A study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas

This thesis is a study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas externally biased with a magnetostatic field. The present study is restricted to the nonlinear phenomena rising from the interaction of electromagnetic waves in the ionized gas. The important effects of nonlinearity are wave-form distortion leads to cross modulation of one wave by a second amplitude-modulated wave.

The nonlinear effects are assumed to be small so that a perturbation method can be used. Boltzmann’s kinetic equation with an appropriate expression for the collision term is solved by expanding the electron distribution function into spherical harmonics in velocity space. In turn, the electron convection current density and the conductivity tensors of the nonlinear ionized gas are found from the distribution function. Finally, the expression for the current density and Maxwell’s equations are employed to investigate the effects of nonlinearity on the propagation of electromagnetic waves in the ionized gas, and also on the reflection of waves from an ionized gas of semi-infinite extent.

High frequency wave propagation in the earth: theory and observation

The wave-theoretical analysis of acoustic and elastic waves refracted by a spherical boundary across which both velocity and density increase abruptly and thence either increase or decrease continuously with depth is formulated in terms of the general problem of waves generated at a steady point source and scattered by a radially heterogeneous spherical body. A displacement potential representation is used for the elastic problem that results in high frequency decoupling of P-SV motion in a spherically symmetric, radially heterogeneous medium. Through the application of an earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as a Weyl integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle:

1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), the critically refracted wave reduces to the classical head wave for flat, homogeneous layers.

2) For gradients more negative than critical, the amplitude of the critically refracted wave decays more rapidly with distance than the classical head wave.

3) For positive, null, and gradients less negative than critical, the amplitude of the critically refracted wave decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

These theoretical results applied to published amplitude data for P-waves refracted by the major crustal and upper mantle horizons (the Pg, P*, and Pn travel-time branches) suggest that the 'granitic' upper crust, the 'basaltic' lower crust, and the mantle lid all have negative or near-critical velocity gradients in the tectonically active western United States. On the other hand, the corresponding horizons in the stable eastern United States appear to have null or slightly positive velocity gradients. The distribution of negative and positive velocity gradients correlates closely with high heat flow in tectonic regions and normal heat flow in stable regions. The velocity gradients inferred from the amplitude data are generally consistent with those inferred from ultrasonic measurements of the effects of temperature and pressure on crustal and mantle rocks and probable geothermal gradients. A notable exception is the strong positive velocity gradient in the mantle lid beneath the eastern United States (2 x 10^-3 sec^-1), which appears to require a compositional gradient to counter the effect of even a small geothermal gradient.

New seismic-refraction data were recorded along a 800 km profile extending due south from the Canadian border across the Columbia Plateau into eastern Oregon. The source for the seismic waves was a series of 20 high-energy chemical explosions detonated by the Canadian government in Greenbush Lake, British Columbia. The first arrivals recorded along this profile are on the Pn travel-time branch. In northern Washington and central Oregon their travel time is described by T = Δ/8.0 + 7.7 sec, but in the Columbia Plateau the Pn arrivals are as much as 0.9 sec early with respect to this line. An interpretation of these Pn arrivals together with later crustal arrivals suggest that the crust under the Columbia Plateau is thinner by about 10 km and has a higher average P-wave velocity than the 35-km-thick, 62-km/sec crust under the granitic-metamorphic terrain of northern Washington. A tentative interpretation of later arrivals recorded beyond 500 km from the shots suggests that a thin 8.4-km/sec horizon may be present in the upper mantle beneath the Columbia Plateau and that this horizon may form the lid to a pronounced low-velocity zone extending to a depth of about 140 km.

Solar flare particle propagation--comparison of a new analytic solution with spacecraft measurements

A new analytic solution has been obtained to the complete Fokker-Planck equation for solar flare particle propagation including the effects of convection, energy-change, corotation, and diffusion with ĸ_r = constant and ĸ_Ɵ ∝ r². It is assumed that the particles are injected impulsively at a single point in space, and that a boundary exists beyond which the particles are free to escape. Several solar flare particle events have been observed with the Caltech Solar and Galactic Cosmic Ray Experiment aboard OGO-6. Detailed comparisons of the predictions of the new solution with these observations of 1-70 MeV protons show that the model adequately describes both the rise and decay times, indicating that ĸ_r = constant is a better description of conditions inside 1 AU than is ĸ_r ∝ r. With an outer boundary at 2.7 AU, a solar wind velocity of 400 km/sec, and a radial diffusion coefficient ĸ_r ≈ 2-8 x 10²⁰ cm²/sec, the model gives reasonable fits to the time-profile of 1-10 MeV protons from "classical" flare-associated events. It is not necessary to invoke a scatter-free region near the sun in order to reproduce the fast rise times observed for directly-connected events. The new solution also yields a time-evolution for the vector anisotropy which agrees well with previously reported observations.

In addition, the new solution predicts that, during the decay phase, a typical convex spectral feature initially at energy T_o will move to lower energies at an exponential rate given by T_KINK = T_oexp(-t/Ƭ_KINK). Assuming adiabatic deceleration and a boundary at 2.7 AU, the solution yields Ƭ_KINK ≈ 100h, which is faster than the measured ~200h time constant and slower than the adiabatic rate of ~78h at 1 AU. Two possible explanations are that the boundary is at ~5 AU or that some other energy-change process is operative.

Guided propagation of surface acoustic waves and piezoelectric field enhancement in ZnO/GaAs systems

The characteristics and dispersion of the distinct surface acoustic waves (SAWs) propagating in ZnO/GaAs heterostructures have been studied experimentally and theoretically. Besides the Rayleigh mode, strong Sezawa modes, which propagate confined in the overlayer, arise due to the smaller sound velocity in ZnO than in the substrate. The design parameters of the structure providing the strongest piezoelectric field at a given depth within the layered system for the different modes have been determined. The piezoelectric field of the Rayleigh mode is shown to be more than 10 times stronger at the interface region of the tailored ZnO/GaAs structure than at the surface region of the bulk GaAs, whereas the same comparison for the first Sezawa mode yields a factor of 2. This enhancement, together with the capacity of selecting waves with different piezoelectric and strain field depth profiles, will facilitate the development of SAW-modulated optoelectronic applications in GaAs-based systems. © 2011 American Institute of Physics.

Wave propagation in nonlinear one-dimensional soil model

The objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of shear modulus as a function of strain. In this paper the meaning of the material stiffness obtained from impact and harmonic excitation tests on a surface slab is discussed. A one-dimensional discrete model with the nonlinear material stiffness is used for this purpose. When a static load is applied followed by an impact excitation, if the amplitude of the impact is very small, the measured wave velocity using the cross-correlation indicates the wave velocity calculated from the tangent modulus corresponding to the state of stress caused by the applied static load. The duration of the impact affects the magnitude of the displacement and the particle velocity but has very little effect on the estimation of the wave velocity for the magnitudes considered herein. When a harmonic excitation is applied, the cross-correlation of the time histories at different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loop under steady-state condition. Copyright © 2008 John Wiley & Sons, Ltd.