Piston Velocities Research Articles

Atomic-level deformation of CuxZr100-x metallic glasses under shock loading

Plastic deformation mechanisms in CuxZr100-x bulk metallic glasses (MGs) subjected to shock are investigated using molecular dynamics simulations. MGs with Cu compositions between 30 and 70 at. % subjected to shock waves generated via piston velocities that range from 0.125 to 2.0 km/s are considered. In agreement with prior studies, plastic deformation is initiated via formation of localized regions of high von Mises shear strain, known as shear transformation zones (STZs). At low impact velocities, but above the Hugoniot elastic limit, STZ nucleation is dispersed behind the shock front. As impact velocity is increased, STZ nucleation becomes more homogeneous, eventually leading to shock-induced melting, which is identified in this work via high atomic diffusivity. The shear stress necessary to initiate plastic deformation within the shock front is independent of composition at shock intensities near the elastic limit but increases with increasing Cu content at high shock intensities. By contrast, both the flow stress in the plastically deformed MG and the critical shock pressure associated with melting behind the shock front are found to increase with increasing Cu content over the entire range of impact velocities. The evolution of the short-range order in the MG samples during shock wave propagation is analyzed using a polydisperse Voronoi tessellation method. Cu-centered polyhedra with full icosahedral symmetry are found to be most resistant to change under shock loading independent of the MG composition. A saturation is observed in the involvement of select Cu-centered polyhedra in the plastic deformation processes at a piston velocity around 0.75 km/s.

Numerical analysis for irreversible processes in a piston-cylinder system

A numerical analysis for the irreversible process in an adiabatic piston-cylinder system has been conducted. Axisymmetric compressible momentum and energy equations were solved numerically to obtain the state quantities of the system using the laminar flow model. The numerical method is based on the combined Implicit Continuous-fluid Eulerian technique and the Arbitrary-Lagrangian-Eulerian method. The computations were performed for a single compression process and a single expansion process with the piston velocities of ±1 m/s, ±2 m/s, ±4 m/s, ±6 m/s, ±8 m/s and ± 10 m/s and for cyclic compression and expansion processes with sinusoidal velocity variation. It is found that the piston velocity has effects on the state quantities of the piston-cylinder system and it experienced an irreversible process when the piston moved with an infinite velocity. However, the process can be treated as a polytropic process and the polytropic exponent is approximately equal to the adiabatic exponent, n ≈ γ when the piston velocity is less than ±10 m/s. In the cyclic process of 10,000 rpm, the internal energy increases 0.037% of the compression work in each cycle.

Shock compression of Cu x Zr100−x metallic glasses from molecular dynamics simulations

The shock response of Cu x Zr100−x (x = 30, 50 and 70) metallic glasses (MGs) is characterized using large-scale molecular dynamics simulations. A wide range of piston velocities U p = 0.125–2.5 km/s are simulated corresponding to shock pressures from 3 to 130 GPa. Independent of composition, the metallic glasses exhibit the following shock wave propagation regimes: (1) single elastic shock wave for U p < 0.25 km/s, (2) split elastic and plastic shock waves for 0.25 < U p < 0.75 km/s and (3) overdriven plastic shock wave with a narrow elastic precursor for U p > 0.75 km/s. Within the split wave and overdriven regimes, the amplitude of the elastic precursor increases with increasing shock intensity, thereby indicating a pressure-dependent yield criterion. Hugoniot states are strongly dependent on the Cu content of the MG with Cu70Zr30 exhibiting a much higher resistance to plastic deformation than either Cu50Zr50 or Cu30Zr70.

Shock induced deformation response of single crystal copper: Effect of crystallographic orientation

We have carried out multimillion atom non-equilibrium molecular dynamics simulations for investigating the effect of crystallographic orientation over the evolution of deformation pathway of single crystal copper under shock compression. Based on symmetry, three different crystallographic directions, 〈100〉,〈110〉 and 〈111〉 are selected and taken as shock directions. Shock Hugoniot points has been calculated and compared among these different directions up to ∼450GPa of shock pressure i.e. piston velocity of 3.0km/s. Orientational anisotropy has been observed for the bulk Cu single crystals shock loaded along these three different directions. Even though this feature may not show up explicitly in experimental investigations which typically measures shock-velocity and density Hugoniot curve, it is apparent from large scale atomistic simulations which measures the temperature Hugoniot curve quite accurately. Differences are observed in the von-Mises strain and stress plot distributions for shock loading of different intensities along the three directions. Large directional dependency is also evident in the evolution mechanism of deformation. Temperature profiles at different piston velocities for the shock front and the shock equilibrated regions shows significantly different and interesting patterns along the three orientations. Maxwell-Boltzmann distribution is observed in the atomic velocities (thereby the temperature profiles also) for both the shock front region as well as the shock equilibrated region for shock loading along all the three directions.

Selective Nanopores in Graphene Sheet for Separation I-129 Isotope from Air

Molecular dynamic simulation was used to investigate the ability of nanoporous graphene membrane in gas separation process. Three di-atomic gases (I2, N2 and O2) were considered, in which different pore sizes were modeled on graphene. The structure contains an impermeable movable wall (piston) to push the mixture gases toward the nanoporous graphene membrane. Two different simulations were carried out, with two different piston velocities. Two key factors in gases separation process are the pore size of graphene and the velocity of movable wall. The results revealed that I-129 separation was improved by using proper size of pore and by decreasing the velocity of movable wall. It was also found that the I-129 gas radionuclides could be completely separated from nitrogen and oxygen molecules in the pore-12 graphene configuration. It was also found that nitrogen was more strongly adsorbed onto the membrane than oxygen, while I-129 was not adsorbed.

Evolution of dislocation mechanisms in single-crystal Cu under shock loading in different directions

Even though there are numerous experiments and molecular dynamic simulations of Cu under shock loading, there appears to be no literature on the evolution of different types of dislocation mechanisms and their mutual interactions during the process of shock loading, which this article addresses through molecular dynamic simulations using the Mishin EAM potential for Cu. Three different directions , , and that have been considered in this article are subjected to shock compression with piston velocities ranging between 0.3–3 km s−1. The evolution of Hirth locks, Lomer–Cottrell locks, cross-slips, jogs, and dislocation-originated stacking-fault tetrahedra are demonstrated in this article for different direction shock loading of single-crystal Cu.

릴리프밸브 쿠션기구 내장형 공기압 실린더의 구동 특성

This paper presents the meter-out and meter-in speed control characteristics of a pneumatic cylinder with relief valve type cushion device. The piston displacement and velocity are measured to investigate high speed driving performance with variation of the pressure setting in relief valve, air supply pressure, load mass, the supply and exhaust flow rate from the cylinder. Also, the internal pressures and temperatures driving pressure and cushion chamber are measured. The piston displacements and velocities of meter-out and meter-in control are compared experimentally determined data. A comparison experimental data meter-out and meter-in control show that a relief valve type cushion device is suitable for high speed pneumatic cylinders. The desired response characteristics of piston displacement and velocity are satisfactory adjust the pressure setting of a relief valve with varying system parameters such as air supply pressure, load mass and controlled flow rate.

Finite-time output feedback control for a pneumatic servo system

In this paper, the position tracking control problem of pneumatic servo systems is investigated. These systems usually have high nonlinearities and unmeasurable piston velocities. Firstly, by using adding a power integrator technique, a global finite-time state feedback controller is proposed. Secondly, based on homogeneous theory, a nonlinear observer is developed to estimate the piston velocity. Finally, the corresponding output feedback controller is derived, which local finite-time stabilizes the position tracking error system. Compared with the conventional backstepping output feedback control scheme, the developed nonsmooth output feedback control scheme offers a faster convergence rate and a better disturbance rejection property. Numerical simulations illustrate the effectiveness of the proposed control scheme.

Dimensionality effects in void-induced explosive sensitivity

The dimensionality of defects in high explosives controls their heat generation and the expansion of deflagrations from them. We compare the behaviour of spherical voids in three dimensions to that of circular voids in two dimensions. The behaviour is qualitatively similar, but the additional focusing along the extra transverse dimension significantly reduces the piston velocity needed to initiate reactions. However, the reactions do not grow as well in three dimensions, so detonations require larger piston velocities. Pressure exponents are seen to be similar to those for the two-dimensional system.

Shockwave Energy Dissipation in Metal–Organic Framework MOF-5

We investigate the response of the metal–organic framework MOF-5 to shock loading along various crystallographic orientations using molecular dynamics simulations with the reactive force field ReaxFF. The dynamical compressive load leads to the volumetric collapse of the open structure of MOF-5. A two-wave structure forms with a pore-collapse wave following a leading elastic precursor that propagates at a higher speed. Interestingly, the propagating pore-collapse wave weakens the elastic precursor indicating possible applications of these materials for shockwave dissipation. Within the two-wave regime, an increase in piston velocity leads to an increase in the velocity of the pore-collapse wave, but the pressure and velocity of the elastic precursor remain unchanged. For piston velocities between 2 and 3 km/s (depending on orientation) the pore-collapse wave catches up with the elastic precursor leading to an overdriven regime. The simulations yield insight into the molecular process of pore collapse and ...

Shock wave propagation and spall failure in single crystal Mg at atomic scales

Large scale molecular dynamics (MD) simulations are carried out to investigate the wave propagation and failure behavior of single crystal Mg under shock loading conditions. The embedded atom method interatomic potential, used to model the Mg systems, is first validated by comparing the predicted Hugoniot behavior with that observed using experiments. The first simulations are carried out to investigate the effect of loading orientation on the wave propagation and failure behavior by shock loading the system along the [0001] direction (c-axis) and the [101¯0] direction using a piston velocity of 1500 m/s. The spall strength (peak tensile pressure prior to failure) is predicted to be higher for loading along the [101¯0] direction than that predicted for loading along the [0001] direction. To investigate the effect of shock pressure on the failure behavior and spall strength of the metal, the MD simulations are carried out using piston velocities of 500 m/s, 1000 m/s, 1500 m/s, and 2000 m/s for loading along the c-axis. The results indicate that the higher piston velocities result in higher shock pressures, and the predicted values for the spall strength decrease with an increase in the shock pressure. In addition, the simulations reveal that the various piston velocities result in variations in the interactions between the reflected waves and the tail of the pressure waves and, hence, variations in the failure behavior. In addition, MD simulations are also carried out to investigate the effect of temperature on the wave propagation behavior and spall strength by equilibrating the initial system at temperatures of 300 K, 600 K, and 800 K prior to shock loading simulations using a piston velocity of 1000 m/s. The results suggest a decrease in spall strength of the single crystal metal with an increase in the initial temperature of the system. The strain rates generated, the evolution of temperature, the variations in the wave interactions, and the spall strength computed using MD simulations are discussed and compared to experimental results in the literature.

Shock Loading of Granular Ni/Al Composites. Part 2: Shock-Induced Chemistry

We use molecular dynamics simulations to characterize the chemical processes resulting from the shock compaction of a loosely packed granular reactive composite of Ni and Al. For all of the impact strengths studied (with piston velocities up in the range 0.5–2.5 km/s), we find that reactions initiate in the vicinity of the collapsed pores. For the lowest impact velocities (up ≤ 0.75 km/s), the reactions that initiate at the collapsed pores subsequently slow down as thermal transport dissipates the initial temperature excursion and outpaces the exothermic energy release rate. At intermediate impact velocities (up ≈ 1.0 km/s), the localization of thermal kinetic energy is sufficient to establish a reaction rate that is self-sustaining in exothermic energy release, and the sample reacts within a few nanoseconds. At the highest impact velocities (up ≥ 1.5 km/s), the localization of translational kinetic energy as well as thermal energy following pore collapse drives the rapid propagation of the reaction from the collapsed pores.

Shock Wave-Induced Damage of a Protein by Void Collapse

In this study, we report on a series of molecular dynamics simulations that were used to examine the effects of shock waves on a membrane-bound ion channel. A planar shock wave was found to compress the ion channel upon impact, but the protein geometry resembles the crystal structure as soon as the solvent density begins to dissipate. When a void was placed in close proximity to the membrane, the shock wave proved to be more destructive to the protein due to formation of a nanojet that results from the asymmetric collapse of the void. The nanojet was able to cause significant structural changes to the protein even at low piston velocities that are not able to directly cause poration of the membrane.

Mesoscale numerical modeling of plastic bonded explosives under shock loading

Mesoscale responses of plastic bonded explosives under shock loading are investigated using material point method as implemented in the Uintah Computational Framework. The two-dimensional geometrical model which can approximately reflect the mesoscopic structure of plastic bonded explosives was created based on the Voronoi tessellation. Shock loading for the explosive was performed by a piston moving at a constant velocity. For the purpose of investigating the influence of shock strength on the responses of explosives, two different velocities for the piston were used, 200 m/s and 400 m/s, respectively. The simulation results indicate that under shock loading there forms some stress localizations on the grain boundary of explosive. These stress localizations lead to large plastic deformations, and the plastic strain energy transforms to thermal energy immediately, causing temperature to rise rapidly and form some hot spots on grain boundary areas. The comparison between two different piston velocities shows that with increasing shock strength, the distribution of plastic strain and temperature does not have significant change, but their values increase obviously. Namely, the higher the shock strength is, the higher the hot spot temperature will be.

Shock Loading of Granular Ni/Al Composites. Part 1: Mechanics of Loading

We present molecular dynamics simulations of the thermomechanical response under shock loading of a granular material consisting of laminated Ni/Al grains. We observe two regimes: At low piston velocities (up ≲ 1km/s), the shock wave is diffuse, and the width of the shock front decreases with increasing piston velocity. Beyond a critical shock strength, however, the width remains relatively constant at approximately the mean grain radius. This change in behavior follows from an evolution of the mechanism of compaction with increasing insult strength. Furthermore, the mechanism evolves from plastic deformation-mediated pore collapse for relatively weak shocks, to solid extrusion and fluid ejecta filling pores ahead of the shock front at intermediate strengths, and finally to atomic jetting into the pore for very strong shocks (up ≳ 2 km/s). High-energy fluid ejecta into pores leads to the formation of flow vorticity and can result in a large fraction of the input energy localizing into translational kinetic energy components including the formation of hot spots. This has implications for the mechanical mixing of Ni and Al in these reactive composites.

Experimental validation of a magnetorheological energy absorber design optimized for shock and impact loads

A linear stroke adaptive magnetorheological energy absorber (MREA) was designed, fabricated and tested for intense impact conditions with piston velocities up to 8 m s−1. The performance of the MREA was characterized using dynamic range, which is defined as the ratio of maximum on-state MREA force to the off-state MREA force. Design optimization techniques were employed in order to maximize the dynamic range at high impact velocities such that MREA maintained good control authority. Geometrical parameters of the MREA were optimized by evaluating MREA performance on the basis of a Bingham-plastic analysis incorporating minor losses (BPM analysis). Computational fluid dynamics and magnetic FE analysis were conducted to verify the performance of passive and controllable MREA force, respectively. Subsequently, high-speed drop testing (0–4.5 m s−1 at 0 A) was conducted for quantitative comparison with the numerical simulations. Refinements to the nonlinear BPM analysis were carried out to improve prediction of MREA performance.

Elastic-plastic collapse of super-elastic shock waves in face-centered-cubic solids

Shock waves in the [110] and [111] directions of single-crystal Al samples were studied using molecular dynamics (MD) simulations. Piston-driven simulations were performed to investigate the split shock-wave regime. At low piston velocities, the material is compressed initially to a metastable over-compressed elastic state leading to a super-elastic single shock wave. This metastable elastic state later collapses to a plastic state resulting in the formation of a two-wave structure consisting of an elastic precursor followed by a slower plastic wave. The single two-zone elastic-plastic shock-wave regime appearing at higher piston velocities was studied using moving window MD. The plastic wave attains the same average speed as the elastic precursor to form a single two-zone shock wave. In this case, repeated collapse of the highly over-compressed elastic state near the plastic shock front produces ultrashort triangle pulses that provide the pressure support for the leading elastic precursor.

On shock response of nano-void closed/open cell copper material: Non-equilibrium molecular dynamic simulations

Non-equilibrium molecular-dynamic simulations were carried out on model three-dimensional nano-void copper material with different idealised pore structure and porosity to highlight differences in response behaviour between them when subjected to various piston velocities simulating planar shock loading of different intensities. This article demonstrates and explains from a mechanistic perspective the differences in response observed with respect to Hugoniot elastic limits, dislocation line and jet formation, void collapse mechanism and hot spot generation, specific volume, partial recrystallisation and temperature evolution in void collapsed regions, shock and particle velocity curves.

Experimental validation of a magnetorheological energy absorber design analysis

A key challenge when designing linear stroke magnetorheological energy absorbers for high-speed impact is that high piston speeds in linear stroke magnetorheological energy absorbers induce high Reynolds number flows in the magnetic valve of the magnetorheological energy absorber, so that achieving high controllable dynamic range can be a design challenge. So far, the research on magnetorheological energy absorbers has typically assumed that the off-state force increases linearly with piston velocity. But at the higher piston velocities occurring in impact events, the off-state damping exhibits nonlinear velocity squared damping effects. This problem was recognized in our prior work, where it was shown that minor losses are important contributing factors to off-state damping. In this study, a nonlinear analytical magnetorheological energy absorber model is developed based on a Bingham-plastic nonlinear flow model combined with velocity squared dependent minor loss factors. This refined model is denoted as the Bingham-plastic nonlinear flow model with minor losses. From this Bingham-plastic nonlinear flow model with minor losses, an effective design strategy is presented for conventional magnetorheological energy absorbers. The Bingham-plastic nonlinear flow model with minor losses is validated via computational fluid dynamics simulation, so that magnetorheological energy absorber performance can be analytically verified before being manufactured. The magnetorheological energy absorber is fabricated and tested up to an effective piston velocity of 5 m/s by using the high-speed drop tower facility at the GM R&amp;D Center. Comparison of our analysis with measured data is conducted, and the effective design of the magnetorheological energy absorber using the Bingham-plastic nonlinear flow model with minor losses is validated.

Atomistic investigation of phase transition in solid argon induced by shock wave transmission

A molecular dynamics (MD) simulation is employed to study the phase transition process in argon induced by shock wave transmission. Deriving the relation between the shock and piston velocities, the theoretical equation of state for argon is presented. Also, argon equation of state is obtained by measuring the quantities directly from simulations to be able to detect the phase transitions. The phase transition is also detected by using argon phase diagram and free energy calculations. A comparison shows good agreement between the theoretical and MD results for the phase transitions. Based on these simulations, it is concluded that under a shock wave transmission with suitable energy, the solid argon experiences a phase transition from solid to liquid and another from liquid to supercritical fluid. By reflecting the shock wave back at the end of its passage, the whole argon may reach the supercritical state.