Multiscale dislocation dynamics simulations of shock compression in copper single crystal

Ultra short pulse shock wave propagation, plastic deformation and evolution of dislocations in copper single crystals with (0 0 1), (0 1 1) and (1 1 1) orientations are investigated using multiscale dislocation dynamics plasticity analyses. The effects of peak pressure, pulse duration, crystal anisotropy and the nonlinear elastic properties on the interaction between shock wave and dislocations are investigated. The results of our calculations show that the dislocation density has a power law dependence on pressure with a power of 1.70 and that the dislocation density is proportional to pulse duration and sensitive to crystal orientation. These results are in very good agreement with the analytical predications of Meyers et al. [Meyers, M.A., Gregori, F., Kad, B.K., Schneider, M.S., Kalantar, D.H., Remington, B.A., Ravichandran G., Boehly, T., Wark, J., 2003. Laser-induced shock compression of monocrystalline copper: characterization and analysis. Acta Materialia 51, 1211–1228] and the experimental results of Murr [Murr, L.E., 1981. Residual microstructure-mechanical property relationships in shock loaded metals and alloys. In: Meyers, M.A., Murr, L.E. (Eds.), Shock Waves and High Strain Rate Phenomena in Metals. Plenum, New York, pp. 607–673]. It is shown also that incorporating the effect of crystal anisotropy in the elastic properties results in orientation dependent wave speed and peak pressure. The relaxed configurations of dislocation microstructures show the formation of microbands coincident with the slip planes.     

Mixed-mode interfacial adhesive strength of a thin film on an anisotropic substrate

The mixed-mode interfacial adhesion strength between a gold (Au) thin film and an anisotropic passivated silicon (Si) substrate is measured using laser-induced stress wave loading. Test specimens are prepared by bonding a fused silica (FS) prism to the back side of a 〈1 0 0〉 Si substrate with a thin silicon nitride (Si_xN_y) passivation layer deposited on the top surface. A high-amplitude stress wave is developed by pulsed laser ablation of a sacrificial absorbing layer on one of the lateral surfaces of the FS prism. Due to the negative non-linear elastic properties of the FS, the compressive stress wave evolves into a decompression shock with fast fall time. Careful selection of the incident angle between the pulse and the FS/Si interface generates a mode-converted shear wave in refraction, subjecting the Si_xN_y/Au thin film interface to dynamic mixed-mode loading, sufficient to cause interfacial fracture. A detailed analysis of the anisotropic wave propagation combined with interferometric measurements of surface displacements enables calculation of the interfacial stresses developed under mixed-mode loading. The mixed-mode interfacial strength is compared to the interfacial strength measured under purely tensile loading.     

Polychromatic microdiffraction analysis of defect self-organization in shock deformed single crystals

A spatially resolved X-ray diffraction method – with a submicron 3D resolution together with SEM and OIM analysis are applied to understand the arrangements of voids, geometrically necessary dislocations and strain gradient distributions in samples of Al (1 2 3) and Cu (0 0 1) single crystals shocked to incipient spallation fracture. We describe how geometrically necessary dislocations and the effective strain gradient alter white beam Laue patterns of the shocked materials. Several distinct structural zones are observed at different depths under the impact surface. The density of geometrically necessary dislocations (GNDs) is extremely high near the impact and back surface of the shock recovered crystals. The spall region is characterized by a large density of mesoscale voids and GNDs. The spall region is separated from the impact and back surfaces by compressed regions with high total dislocation density but lower GNDs density. Self-organization of shear bands is observed in the shock recovered Cu single crystal.     

Size dependent mechanical behaviour of tantalum

The size dependence of deformation of Ta was studied using compression tests of focused ion beam (FIB) machined microcolumns. Columns with diameters between 0.5 and 8 μm with 〈1 1 1〉 and 〈1 0 0〉 orientations along the column axis were tested. By comparing results of bcc Ta columns with results from previous experiments on fcc metals it was found that Ta shows significantly higher normalized yield stresses in combination with a weaker sample size dependence. The differences between bcc and fcc metals can be attributed to the different dislocation behaviour of bcc metals, especially to the lower mobility of screw dislocations.     

A model for plasticity kinetics and its role in simulating the dynamic behavior of Fe at high strain rates

The recent diagnostic capability of the Omega laser to study solid-solid phase transitions at pressures greater than 10 GPa and at strain rates exceeding 10⁷ s⁻¹ has also provided valuable information on the dynamic elastic-plastic behavior of materials. We have found, for example, that plasticity kinetics modifies the effective loading and thermodynamic paths of the material. In this paper we derive a kinetics equation for the time-dependent plastic response of the material to dynamic loading, and describe the model’s implementation in a radiation-hydrodynamics computer code. This model for plasticity kinetics incorporates the Gilman model for dislocation multiplication and saturation. We discuss the application of this model to the simulation of experimental velocity interferometry data for experiments on Omega in which Fe was shock compressed to pressures beyond the α-to-ε phase transition pressure. The kinetics model is shown to fit the data reasonably well in this high strain rate regime and further allows quantification of the relative contributions of dislocation multiplication and drag. The sensitivity of the observed signatures to the kinetics model parameters is presented.     

A pulse-expansion wave tube for nucleation studies at high pressures

 The design and performance of a new pulse-expansion wave tube for nucleation studies at high pressures are described. The pulse-expansion wave tube is a special shock tube in which a nucleation pulse is formed at the endwall of the high pressure section. The nucleation pulse is due to reflections of the initial shock wave at a local widening situated in the low pressure section at a short distance from the diaphragm. The nucleation pulse has a duration of the order of 200 μs, while nucleation pressures that can be achieved range from 1 to 50 bar total pressure. Droplet size and droplet number density can accurately be determined by a 90°-Mie light scattering method and a light extinction method. The range of nucleation rates that can be measured is 10⁸ cm^-3 s^-1<J<10¹¹ cm^-3 s^-1. We will illustrate the functioning and possibilities of the new pulse-expansion wave tube by nucleation rate measurements in the gas-vapour mixture nitrogen/water in the temperature range 200–260 K, and in the mixture methane/n-nonane as a function of supersaturation S at various total pressures up to 40 bar and temperatures around 240 K. Received: 5 June 1996/Accepted: 9 December 1996     

Rate dependent critical strain energy density factor of Huanglong limestone

Critical strain energy density of rock can be defined as a fundamental parameter in rock fracture mechanics, an intrinsic material property related to resistance to crack initiation and propagation. By means of the three-point bending experiments, the critical strain energy density factor of Huanglong limestone was measured over a wide range of loading rates from 8.97 × 10⁻⁴ MPam^1/2 s⁻¹ to 1.545 MPam^1/2 s⁻¹. According to the approximate relationship between static and dynamic critical strain energy density factor of Huanglong limestone, relationship between the growth velocity of crack and magnitude of load is obtained. The main conclusions are summarized as follows: (1) when the loading rate is higher than 0.0279 MPam^1/2 s⁻¹, the critical strain energy density factor of rock increased markedly with increasing loading rate. However, when loading rate is lower than 0.0279 MPam^1/2 s⁻¹, the critical strain energy density factor slightly increased with an increase in loading rate. It is found from experimental results that the critical strain energy density factor is linear proportional to the exponential expression of loading rate, (2) for Huanglong limestone, when the growth velocity of crack is lower than 100 m/s, value of the maximum load was nearly a constant. However, when the growth velocity of crack is higher than 1000 m/s, value of the maximum load dramatically increases with increasing the crack growth velocity, and (3) the critical SED of Huanglong limestone is higher as the loading rate is higher.     

Indentation Schmid factor and orientation dependence of nanoindentation pop-in behavior of NiAl single crystals

Instrumented nanoindentation techniques have been widely used to characterize the small-scale mechanical behavior of materials. The elastic-plastic transition during nanoindentation is often indicated by a sudden displacement burst (pop-in) in the measured load-displacement curve. In defect-free single crystals, the pop-in is believed to be the result of homogeneous dislocation nucleation because the maximum shear stress corresponding to the pop-in load approaches the theoretical strength of the materials and because the statistical distribution of pop-in stresses is consistent with what is expected for a thermally activated process of homogeneous dislocation nucleation. This paper investigates whether this process is affected by crystallography and stress components other than the resolved shear stress. A Stroh formalism coupled with the two-dimensional Fourier transformation is used to derive the analytical stress fields in elastically anisotropic solids under Hertzian contact, which allows the determination of an indentation Schmid factor, namely, the ratio of maximum resolved shear stress to the maximum contact pressure. Nanoindentation tests were conducted on B2-structured NiAl single crystals with different surface normal directions. This material was chosen because it deforms at room temperature by {1 1 0}〈0 0 1〉 slip and thus avoids the complexity of partial dislocation nucleation. Good agreement is obtained between the experimental data and the theoretically predicted orientation dependence of pop-in loads based on the indentation Schmid factor. Pop-in load is lowest for indentation directions close to 〈1 1 1〉 and highest for those close to 〈0 0 1〉. In nanoindentation, since the stress component normal to the slip plane is typically comparable in magnitude to the resolved shear stress, we find that the pressure sensitivity of homogeneous dislocation nucleation cannot be determined from pop-in tests. Our statistical measurements generally confirm the thermal activation model of homogeneous dislocation nucleation. That is, the extracted dependence of activation energy on resolved shear stress is almost the same for all the indentation directions considered in this study, except for those close to 〈0 0 1〉. Because very high pop-in loads are measured for orientations close to 〈0 0 1〉, which implies a large contact area at pop-in, there is a higher probability of activating pre-existing dislocations in these orientations, which may explain the discrepancy near 〈0 0 1〉.     

Effects of heat flux,mass flux,vapor quality,and saturation temperature on flow boiling heat transfer in microchannels

Flow boiling heat transfer with the refrigerants R-134a and R-245fa in copper microchannel cold plate evaporators is investigated. Arrays of microchannels of hydraulic diameter 1.09 and 0.54 mm are considered. The aspect ratio of the rectangular cross section of the channels in both test sections is 2.5. The heat transfer coefficient is measured as a function of local thermodynamic vapor quality in the range −0.2 to 0.9, at saturation temperatures ranging from 8 to 30 °C, mass flux from 20 to 350 kg m⁻² s⁻¹, and heat flux from 0 to 22 W cm⁻². The heat transfer coefficient is found to vary significantly with heat flux and vapor quality, but only slightly with saturation pressure and mass flux for the range of values investigated. It was found that nucleate boiling dominates the heat transfer. In addition to discussing measurement results, several flow boiling heat transfer correlations are also assessed for applicability to the present experiments.     

Enhanced nucleate boiling on copper micro-porous surfaces

Saturation boiling of PF-5060 dielectric liquid on Cu micro-porous surface layers (95, 139, 171, 197 and 220-μm thick) is investigated. These layers are deposited on 10 × 10 mm Cu substrates using two-stage electrochemical process. The basic micro-structure, obtained in the first stage using current density of 3 A/cm² for 15–44 s, depending on thickness, is strengthened by continuing electrochemical deposition using much lower current density for 10’s of minutes. For conditioned surface layers, after a few successive boiling tests, the pool boiling curves are reproducible and the temperature excursion prior to boiling incipience is either eliminated or reduced <7 K. Present nucleate boiling results are markedly better than those reported for dielectric liquids on micro- and macro-structured surfaces. Present values of CHF (22.7–27.8 W/cm²) and h_MNB (2.05–13.5 W/cm² K) are ∼40–70% higher than and >17 times those reported on plane surfaces (<16 W/cm² and ∼0.8 W/cm² K). Best results are those of the 171-μm thick layer: CHF of 27.8 W/cm² occurs at ΔT_sat of only 2.1 K and h_MNB of 13.5 W/cm² K occurs at ΔT_sat = 2.0 K.     

Non-coaxiality of strain increment and stress directions in cross-anisotropic sand

An experimental program was carried out in a recently developed torsion shear apparatus to study the non-coaxiality of strain increment and stress directions in cross-anisotropic deposits of Fine Nevada sand. Forty-four drained torsion shear tests were performed at constant mean confining stress, σ_m, constant intermediate principal stress ratios, as indicated by b = (σ₂ − σ₃)/(σ₁ − σ₃), and constant principal stress directions, α. The experiments were performed on large hollow cylinder specimens deposited by dry pluviation and tested in an automated torsion shear apparatus. The specimens had height of 40 cm, and average diameter of 20 cm, and wall thickness of 2 cm. The stress–strain behavior of Fine Nevada sand is presented for discrete combinations of constant principal stress direction, α, and intermediate principal stress. The effects of these two variables on the non-coaxiality are presented. The experiments show that the directions of the strain increments do not in general coincide with the directions of stresses, and there is a switch from one to the other side between the two quantities.     

U–6Nb shear stress relaxation in compression waves (IJP 585-AV)

When uranium alloyed with 6-wt% niobium (U–6Nb) is rapidly compressed in uniaxial strain experiments, shear stress is observed to relax with a characteristic time of 30 ± 7 ns. In shock wave experiments, this relaxation inhibits the development of an elastic precursor commonly seen in other materials. When U–6Nb is cold-rolled to pre-twin and significantly increase the density of dislocations in the material, stress relaxation effects are diminished suggesting that twinning causes relaxation in the un-worked material. Separate ramp wave compression experiments produce effects that agree with those observed in shock-loading experiments. A phenomenological model is introduced that allows accurate simulation of all experiments. Estimates of residual shear stress after relaxation are obtained.     

The effects of couple stresses on dislocation strain energy

The correspondence theorem which relates the solutions of displacement boundary value problems in classical and couple stress elasticity is formulated and applied to derive the solutions for edge and screw dislocations in an infinite medium. The effects of couple stresses on the dislocation strain energy is evaluated for both types of dislocations. It is shown that within a radius of influence of each dislocation in a metallic crystal with dislocation density of 10¹⁰ cm⁻², the strain energy contribution from couple stresses (excluding the core energy) is about 15% in the case of an edge dislocation, and about 11% in the case of a screw dislocation. It is then shown that couple stresses make large effect on the total work of tractions acting on the dislocation core surface.     

Towards laboratory produced relativistic electron–positron pair plasmas

We review recent experimental results on the path to producing electron–positron pair plasmas using lasers. Relativistic pair-plasmas and jets are believed to exist in many astrophysical objects and are often invoked to explain energetic phenomena related to Gamma Ray Bursts and Black Holes. On earth, positrons from radioactive isotopes or accelerators are used extensively at low energies (sub-MeV) in areas related to surface science positron emission tomography and basic antimatter science. Experimental platforms capable of producing the high-temperature pair-plasma and high-flux jets required to simulate astrophysical positron conditions have so far been absent. In the past few years, we performed extensive experiments generating positrons with intense lasers where we found that relativistic electron and positron jets are produced by irradiating a solid gold target with an intense picosecond laser pulse. The positron temperatures in directions parallel and transverse to the beam both exceeded 0.5 MeV, and the density of electrons and positrons in these jets are of order 10¹⁶ cm⁻³ and 10¹³ cm⁻³, respectively. With the increasing performance of high-energy ultra-short laser pulses, we expect that a high-density, up to 10¹⁸ cm⁻³, relativistic pair-plasma is achievable, a novel regime of laboratory-produced hot dense matter.     

Discrete element simulation of shock wave propagation in polycrystalline copper

The implementation of the characteristic of compressive plasticity into the Discrete Element Code, DM2, while maintaining its quasi-molecular scheme, is described. The code is used to simulate the shock compression of polycrystalline copper at 3.35 and 11.0 GPa. The model polycrystal has a normal distribution of grain sizes, with mean diameter 14 μm, and three distinct grain orientations are permitted with respect to the shock direction; 〈1 0 0〉, 〈1 1 0〉, and 〈1 1 1〉. Particle velocity dispersion (PVD) is present in the shock-induced flow, attaining its maximum magnitude at the plastic wave rise. PVD normalised to the average particle velocity of and are yielded for the 3.35 and 11.0 GPa shocks, respectively, and are of the same order as those seen in the experiment. Non-planar elastic and plastic wave fronts are present, the distribution in shock front position increasing with propagation distance. The rate of increase of the spread in shock front positions is found to be significantly smaller than that seen in probabilistic calculations on nickel polycrystals, and this difference is attributed, in the main, to grain interaction. Reflections at free surfaces yield a region of tension near to the target free surface. Due to the dispersive nature of the shock particle velocity and the non-planarity of the shock front, the tensile pressure is distributed. This may have implications for the spall strength, which are discussed. Simulations reveal a transient shear stress distribution behind the shock front. Such a distribution agrees with that put forward by Lipkin and Asay to explain the quasi-elastic reloading phenomenon. Simulation of reloading shocks show that the shear stress distribution can give rise to quasi-elastic reloading on the grain scale.     

Grain size,strain rate,and temperature dependence of flow stress in ultra-fine grained and nanocrystalline Cu and Al: Synthesis,experiment, and constitutive modeling

For the first time, high quality bulk nanocrystalline (nc) fcc metals, with least amounts of imperfections, exhibiting high strength and ductility at room and different temperatures, under quasi-static and dynamic types of loading, were prepared and a comprehensive study on their post-yield mechanical properties was performed. This investigation included study of the effect of temperature on stress–strain responses of mechanically milled bulk nc Cu and Al. The samples after preparation through mechanical milling and consolidation processes were subjected to uniaxial compressive loading at quasi-static and dynamic strain rates of 10⁻² s⁻¹ and 1840–3105 s⁻¹, respectively, at temperatures ranging from 223 to 523 K. In both materials strong dependency of flow stress to temperature was observed; this dependency was rather more pronounced when the materials were tested at the quasi-static strain rate. Further, a new grain size and temperature dependent viscoplastic phenomenological constitutive equation, Khan–Liang–Farrokh (KLF) model was developed based on the Khan–Huang–Liang (KHL) constitutive equation. The model was featured to correlate different characteristic behaviors of polycrystalline materials in the plastic regime, as the result of grain refinement. In addition, the viscoplastic responses of bulk Cu and Al of different grain sizes (from sub-micron to nanometer range), and those from bulk nc Cu and Al at different strain rates (quasi-static to dynamic), recently published (21 and 22), were simulated using the newly developed equation. The results confirmed reasonable capability of the developed model to correlate a wide spectrum of the viscoplastic responses of these fcc metals.     

Studying astrophysical collisionless shocks with counterstreaming plasmas from high power lasers

Collisions of high Mach number flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large. Therefore, most shock waves in astrophysics are “collisionless”, since they form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at the laser facilities can achieve the conditions necessary for the formation of collisionless shocks, and will provide a unique avenue for studying the nonlinear physics of collisionless shock waves. We are performing a series of experiments at the Omega and Omega-EP lasers, in Rochester, NY, with the goal of generating collisionless shock conditions by the collision of two high-speed plasma flows resulting from laser ablation of solid targets using ∼10¹⁶ W/cm² laser irradiation. The experiments will aim to answer several questions of relevance to collisionless shock physics: the importance of the electromagnetic filamentation (Weibel) instabilities in shock formation, the self-generation of magnetic fields in shocks, the influence of external magnetic fields on shock formation, and the signatures of particle acceleration in shocks. Our first experiments using Thomson scattering diagnostics studied the plasma state from a single foil and from double foils whose flows collide “head-on”. Our data showed that the flow velocity and electron density were 10⁸ cm/s and 10¹⁹ cm⁻³, respectively, where the Coulomb mean free path is much larger than the size of the interaction region. Simulations of our experimental conditions show that weak Weibel mediated current filamentation and magnetic field generation were likely starting to occur. This paper presents the results from these first Omega experiments.     

Pool boiling heat transfer of ultra-light copper foam with open cells

High speed visualizations and thermal performance studies of pool boiling heat transfer on copper foam covers were performed at atmospheric pressure, with the heating surface area of 12.0 mm by 12.0 mm, using acetone as the working fluid. The foam covers have ppi (pores per inch) from 30 to 90, cover thickness from 2.0 to 5.0 mm, and porosity of 0.88 and 0.95. The surface superheats are from −20 to 190 K, and the heat fluxes reach 140 W/cm². The 30 and 60 ppi foam covers show the periodic single bubble generation and departure pattern at low surface superheats. With continuous increases in surface superheats, they show the periodic bubble coalescence and/or re-coalescence pattern. Cage bubbles were observed to be those with liquid filled inside and vented to the pool liquid. For the 90 ppi foam covers, the bubble coalescence takes place at low surface superheats. At moderate or large surface superheats, vapor fragments continuously escape to the pool liquid.     

Three phase liquid–liquid–gas flows in 5.6 mm and 7 mm inner diameter pipes

Three phase liquid–liquid–gas flow maps in pipes of medium inner diameters (5.6 mm and 7 mm), are presented. A low viscosity paraffin oil (4.5 × 10⁻³ Pa s viscosity and 818.5 kg m⁻³ density at 20 °C), deionised water and air are flowing concurrently in Schott Duran^® glass pipes. A decreasing pipe diameter changes the flow pattern maps and also the behavior of the transition boundaries. Flow patterns are determined by high speed photography. To illuminate the pipe, laser induced fluorescence (LIF) is applied. The laser sheet is cutting through the axial vertical plane of the pipe. The laser light excites a fluorescent dye (uranine) in the water phase to separate the phases optically. The resulting flow maps are compared with literature data and a theoretical model.