Adiabatic shear localization: occurrence, theories, and applications

On the origin of shear banding instability in metallic glasses

To uncover the physical origin of shear-banding instability in metallic glass (MG), a theoretical description of thermo-mechanical deformation of MG undergoing one-dimensional simple shearing is presented. The coupled thermo-mechanical model takes into account the momentum balance, the energy balance and the dynamics of free volume. The interplay between free-volume production and temperature increase being two potential causes for shear-banding instability is examined on the basis of the homogeneous solution. It is found that the free-volume production facilitates the sudden increase in the temperature before instability and vice versa. A rigorous linear perturbation analysis is used to examine the inhomogeneous deformation, during which the onset criteria and the internal length and time scales for three types of instabilities, namely free-volume softening, thermal softening and coupling softening, are clearly revealed. The shear-banding instability originating from sole free-volume softening takes place easier and faster than that due to sole thermal softening, and dominates in the coupling softening. Furthermore, the coupled thermo-mechanical shear-band analysis does show that an initial slight distribution of local free volume can incur significant strain localization, producing a shear band. During such a localization process, the local free-volume creation occurs indeed prior to the increase in local temperature, indicating that the former is the cause of shear localization, whereas the latter is its consequence. Finally, extension of the above model to include the shear-induced dilatation shows that such dilatation facilitates the shear instability in metallic glasses.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Few-cycle relativistic solitons in plasmas

A set of exact one-dimensional solutions to coupled nonlinear equations describing the propagation of a relativistic ultrashort circularly polarized laser pulse in a cold collisionless and bounded plasma where electrons have an initial velocity in the laser propagating direction is presented. The solutions investigated here are in the form of quickly moving envelop solitons at a propagation velocity comparable to the light speed. The features of solitons in both underdense and overdense plasmas with electrons having different given initial velocities in the laser propagating direction are described. It is found that the amplitude of solitons is larger and soliton width shorter in plasmas where electrons have a larger initial velocity. In overdense plasmas, soliton duration is shorter, the amplitude higher than that in underdense plasmas where electrons have the same initial velocity.

μ^-子催化核聚变中强脉冲激光对介原子μ^3He的电离

数值求解了一维含时的Schroedinger方程，研究了μ子催化核聚变反应中激光强度和波长对介原子μ^3He电离的影响．发现当激光强度为10^19-10^23W/cm^2量级时，介原子μ^3He有2．7％左右的电离率；当激光强度达到6．0×10^24W/cm^2时，对介原子μ^3He有显著的电离，并且电离率随着激光的强度、波长而递增，进而会有效提高μ子的催化效率．

双色激光场中激光强度和脉宽对线性多原子分子离子电离的影响

利用强激光场电离和离解分子来研究分子激发态的波包结构是强场物理的重要研究方向。利用短时指数传播子对称分割法和快速傅里叶变换技术。数值求解了一维含时Schr（oe）dinger方程，探讨了双色激光场中激光的基波和谐波强度之间的不同配比以及脉宽对线性多原子分子离子电离的影响。理论计算结果表明：基波和谐波的相对相位为π时，尽管随着激光的基波和谐波强度之间配比的变化，电离几率随原子间距变化的趋势基本保持不变，但在一定的激光基波强度下(1．2×10^13～1．2×10^15W／cm^2)，激光基波强度的变化可以明显

Generation of plasma intrinsic oscillation at the front surface of a target irradiated by a circularly polarized laser pulse

In laser-target interaction, the effects of laser intensity on plasma oscillation at the front surface of targets have been investigated by one-dimensional particle in cell simulations. The periodical oscillations of the ion density and electrostatic field at the front surface of the targets are reported for the first time, which is considered as an intrinsic property of the target excited by the laser. The oscillation period depends only on initial plasma density and is irrelevant with laser intensity. Flattop structures with curves in ion phase space are found with a more intense laser pulse due to the larger amplitude variation of the electrostatic field. A simple but valid model is proposed to interpret the curves.

Standing ultrashort relativistic solitons in overdense plasmas

By using a one-dimensional self-consistent relativistic fluid model, an investigation is made numerically on relativistic electromagnetic solitons with a high intensity in cold overdense plasmas with an electrons' initial velocity opposite to the laser propagating direction. Two types of standing solitons with zero group velocity are found at the given electrons' initial velocities. One is single-humped with a weakly relativistic intensity; the another is multi-humped with a strong relativistic amplitude. The properties of these two types of solitons are presented in detail.

A pseudo-spin model for the wall of cytoskeletal microtubule

A pseudo-spin model is intended to describe the physical dynamics of unbound electrons in the wall of cytoskeletal microtubule (MT). Due to the inherent symmetry of the structure and the electric properties in the MT, one may treat it as a one-dimensional ferroelectric system, and describe the nonlinear dynamics of dimer electric dipoles in one protofilament of the MT by virtue of the double-well potential. Consequently, the physical problem has been mapped onto the pseudo-spin system, and the mean-field approximation has been taken to get some physical results.

Steady state ion acceleration by a circularly polarized laser pulse

The steady state ion acceleration at the front of a cold solid target by a circularly polarized flat-top laser pulse is studied with one-dimensional particle-in-cell (PIC) simulation. A model that ions are reflected by a steady laser-driven piston is used by comparing with the electrostatic shock acceleration. A stable profile with a double-flat-top structure in phase space forms after ions enter the undisturbed region of the target with a constant velocity. (C) 2007 Elsevier B.V. All rights reserved.

Multistaged acceleration of ions by circularly polarized laser pulse: Monoenergetic ion beam generation

A multiple-staged ion acceleration mechanism in the interaction of a circularly polarized laser pulse with a solid target is studied by one-dimensional particle-in-cell simulation. The ions are accelerated from rest to several MeV monoenergetically at the front surface of the target. After all the plasma ions are accelerated, the acceleration process is repeated on the resulting monoenergetic ions. Under suitable conditions multiple repetitions can be realized and a high-energy quasi-monoenergetic ion beam can be obtained.

Electron acceleration by a propagating laser pulse in vacuum

Electrons accelerated by a propagating laser pulse of linear or circular polarization in vacuum have been investigated by one-dimensional particle-in-cell simulations and analytical modeling. A stopping target is used to stop the laser pulse and extract the energetic electrons from the laser field. The effect of the reflected light is taken into account. The maximum electron energy depends on the laser intensity and initial electron energy. There is an optimal acceleration length for electrons to gain maximum energy where electrons meet the peak of the laser pulse. The optimal acceleration length depends strongly on the laser pulse duration and amplitude. (C) 2007 American Institute of Physics.