Correlative reference model and molecular dynamics simulation of dislocation emission process

A correlative reference model for computer molecular dynamics simulations is proposed. Based on this model, a flexible displacement boundary scheme is introduced and the dislocations emitted from a crack tip can continuously pass through the border of the inner discrete atomic region and pile up at the outer continuum region. The effect of the emitted dislocations within the plastic zone on the inner atomistic region can be clearly demonstrated. The simulations for a molybdinum crystal show that a full dislocation in a bcc crystal is dissociated into three partial dislocations and interaction between the crack and the emitted dislocations results in gradual decrease of the local stress intensity factor.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Two-way coupling model for shock-induced laminar boundary-layer flows of a dusty gas

The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute two-phase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersed-phase flow fields (including the number density of particles) in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.

Validation of component assembly model and extension to plasticity

Material potential energy is well approximated by '' pair-functional '' potentials. During calculating potential energy, the orientational and volumetric components have been derived from pair potentials and embedding energy, respectively. Slip results in plastic deformation, and slip component has been proposed accordingly. Material is treated as a component assembly, and its elastic, plastic and damage properties are reflected by different components respectively. Material constitutive relations are formed by means of assembling these three kinds of components. Anisotropy has been incorporated intrinsically via the concept of component. Theoretical and numerical results indicate that this method has the capacity of reproducing some results satisfactorily, with the advantages of physical explicitness, etc. (c) 2007 Elsevier Ltd. All rights reserved.

Failure process analysis of cement under explosive loading with a generalized driven nonlinear threshold model

The microstructural heterogeneity and stress fluctuation play important roles in the failure process of brittle materials. In this paper, a generalized driven nonlinear threshold model with stress fluctuation is presented to study the effects of microstructural heterogeneity on continuum damage evolution. As an illustration, the failure process of cement material under explosive loading is analyzed using the model. The result agrees well with the experimental one, which proves the efficiency of the model.

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

An Atomistic-Continuum Hybrid Simulation Of Fluid Flows Over Superhydrophobic Surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper, an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces, in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths, which cannot be obtained by molecular dynamics simulation alone.

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.

An atomistic-continuum hybrid simulation of fluid flows over superhydrophobic surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths which cannot be obtained by molecular dynamics simulation alone.

Two-component model for the interstellar scattering in direction to the Vela pulsar

This paper discusses a rigorous treatment of the refractive scintillation of pulsar PSR B0833-45 caused by a two-component interstellar scattering medium. It is assumed that the interstellar scattering medium is composed of a thin screen ISM and an extended interstellar medium. We consider that the scattering of the thin screen concentrates in a thin layer presented by a delta function distribution and that the scattering density of the extended irregular medium satisfies the Gaussian distribution. We investigate and develop equations for the flux density structure function corresponding to this two-component ISM geometry in the scattering density distribution and compare our result with that of the Vela pulsar observations. We conclude that the refractive scintillation caused by this two-component ISM scattering gives a more satisfactory explanation for the observed flux density variation of the Vela pulsar than does the single extended medium model. The level of refractive scintillation is strongly sensitive to the distribution of scattering material along the line of sight. The logarithmic slope of the structure function is sensitive to thin screen location and is relatively insensitive to the scattering strength of the thin screen medium. Therefore, the proposed model can be applied to interpret the structure function of flux density observed in pulsar PSR B0833-45. The result suggests that the medium consists of a discontinuous distribution of plasma turbulence embedded in the Vela supernova remnant. Thus our work provides some insight into the distribution of the scattering along the line of sight to the Vela pulsar.

CONTINUUM-THEORIES OF OPTICAL PHONONS AND POLARITONS IN SUPERLATTICES - A BRIEF CRITIQUE - COMMENT

The so-called hydrodynamic (HD) model on optical-phonon modes in superlattices is critically examined. Contrary to the HD model, a comparison between TM polaritons and the Fuchs-Kliewer-type interface modes has shown that the Fuchs-Kliewer interface modes do possess Frohlich potentials.

Isoscalar Giant Monopole Resonance in Relativistic Continuum Random Phase Approximation

The isoscalar giant monopole resonance (ISGMR) in nuclei is studied in the framework of a fully consistent relativistic continuum random phase approximation (RCRPA). In this method the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single particle Green's function technique. The negative energy states in the Dirac sea are also included in the single particle Green's function in the no-sea approximation. The single particle Green's function is calculated numerically by a proper product of the regular and irregular solutions of the Dirac equation. The strength distributions in the RCRPA calculations, the inverse energy-weighted sum rule m(-1) and the centroid energy of the ISGMR in Sn-120 and Pb-208 are analysed. Numerical results of the RCRPA are checked with the constrained relativistic mean field model and relativistic random phase approximation with a discretized spectrum in the continuum. Good agreement between them is achieved.

Importance of self-consistency in relativistic continuum random-phase approximation calculations

A fully consistent relativistic continuum random phase approximation (RCRPA) is constructed, where the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single-particle Green's function technique. The full consistency of the calculations is achieved that the same effective Lagrangian is adopted for the ground state and the excited states. The negative energy states in the Dirac sea are also included in the single-particle Green's function in the no-sea approximation. The currents from the vector meson and photon exchanges and the Coulomb interaction in RCRPA are treated exactly. The spin-orbit interaction is included naturally in the relativistic frame. Numerical results of the RCRPA are checked with the constrained relativistic mean-field theory. We study the effects of the inconsistency, particularly the currents and Coulomb interaction in various collective multipole excitations.