Relativistic Continuum Random Phase Approximation and Applications II. Applications

The fully consistent relativistic continuum random phase approximation (RCRPA) has been constructed in the momentum representation in the first part of this paper. In this part we describe the numerical details for solving the Bethe-Salpeter equation. The numerical results are checked by the inverse energy weighted sum rules in the isoscalar giant monopole resonance, which are obtained from the constraint relativistic mean field theory and also calculated with the integration of the RCRPA strengths. Good agreement between the misachieved. We study the effects of the self-consistency violation, particularly the currents and Coulomb interaction to various collective multipole excitations. Using the fully consistent RCRPA method, we investigate the properties of isoscalar and isovector collective multipole excitations for some stable and exotic from light to heavy nuclei. The properties of the resonances, such as the centroid energies and strength distributions are compared with the experimental data as well as with results calculated in other models.

Relativistic Continuum Random Phase Approximation and Applications I. Formalism

A fully consistent relativistic continuum random phase approximation (RCRPA) is constructed in terms of the Green's function technique. In this method the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single particle Green's function, which includes also the negative states in the Dirac sea in the nose aapproximation. The theoretical formalism of RCRPA and numerical details are presented. The single particle Green's function is calculated numerically by a proper product of regular and irregular solutions of the Dirac equation. The numerical details and the formalism of RCRPA in the momentum representation are presented.

Nitrogen-15 signals of leaf-litter-soil continuum as a possible indicator of ecosystem nitrogen saturation by forest succession and N loads

National Key Research and Development Program [2010CB833502]; National Natural Science Foundation of China [30600071, 40601097, 30590381]; Chinese Academy of Sciences [KZCX2-YW-432, O7V70080SZ, LENOM07LS-01]; GUCAS [O85101PM03]

Brittle-ductile transition of polymers and its percolation model

The brittle-ductile transition (BDT) of particle toughened polymers was extensively studied in terms of morphology, strain rate, and temperature. The calculation results showed that both the critical interparticle distance (IDc) and the brittle-ductile transition temperature (T-BD) of polymers were a function of strain rate. The IDc reduced nonlinearly with increasing strain rate, whereas T-BD increased considerably with increasing strain rate. The effects of temperature and plasticizer concentration on BDT were discussed using a percolation model. The results were in agreement with the experiments.

Effects of Microstructural Heterogeneity on the Spallation Behavior of Materials

It is of utmost importance to understand the spallation behaviour of heterogeneous materials. In this paper, a driven nonlinear threshold model with stress fluctuation is presented to study the effects of microstructural heterogeneity on continuum damage evolution. The spallation behavior of heterogeneity material is analyzed with this model. The heterogeniety of mesoscopic units is characterized in terms of Weibull modulus m of strength distibution and stress fluctuation parameter k. At high stress, the maximum damage increases with m; while at low stress, the maximum damage decreases. In addition, for low stress, severe stress fluctuation causes higher damage; while for high stress, causes lower damage.

Peierls-Nabarro model of interfacial misfit dislocation: An analytic solution

We propose a method to treat the interfacial misfit dislocation array following the original Peierls-Nabarro's ideas. A simple and exact analytic solution is derived in the extended Peierls-Nabarro's model, and this solution reflects the core structure and the energy of misfit dislocation, which depend on misfit and bond strength. We also find that only with beta < 0.2 the structure of interface can be represented by an array of singular Volterra dislocations, which conforms to those of atomic simulation. Interfacial energy and adhesive work can be estimated by inputting ab initio calculation data into the model, and this shows the method can provide a correlation between the ab initio calculations and elastic continuum theory.

Characteristics of Micronozzle Gas Flows

The fluid characteristics of gas flows in the micronozzle whose throat height is 20 μm were investigated by the direct simulation Monte Carlo (DSMC) method. In a series of cases, the dependence of mass flux on the pressure difference was gained, and the DSMC's results show good agreement with the experimental data. The comparison of mass flux and the Mach number contours between the DSMC and Navier-Stokes equations adding slip boundary also reveals quantitatively that the continuum model will be invalid gradually even when the average Knudsen number is smaller than 0.01. As one focus of the present paper, the phenomenon of the multiple expansion-compression waves that comes from the nozzle's divergent part was analyzed in detailed.

Molecular Dynamics Simulation of Peeling a DNA Molecule on Substrate

Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.

3D Anisotropic Elastoplastic-Damage Model and its Application in Simulating the Behavior of Rock Materials

A 3D anisotropic elastoplastic-damage model was presented based on continuum damage mechanics theory. In this model, the tensor decomposition technique is employed. Combined with the plastic yield rule and damage evolution, the stress tensor in incremental format is obtained. The derivate eigenmodes in the proposed model are assumed to be related with the uniaxial behavior of the rock material. Each eigenmode has a corresponding damage variable due to the fact that damage is a function of the magnitude of the eigenstrain. Within an eigenmodes, different damage evolution can be used for tensile and compressive loadings. This model was also developed into finite element code in explicit format, and the code was integrated into the well-known computational environment ABAQUS using the ABAQUS/Explicit Solver. Numerical simulation of an uniaxial compressive test for a rock sample is used to examine the performance of the proposed model, and the progressive failure process of the rock sample is unveiled.

The Flow Theory of Mechanism-Based Strain Gradient Plasticity

The flow theory of mechanism-based strain gradient (MSG) plasticity is established in this paper following the same multiscale, hierarchical framework for the deformation theory of MSG plasticity in order to connect with the Taylor model in dislocation mechanics. We have used the flow theory of MSG plasticity to study micro-indentation hardness experiments. The difference between deformation and flow theories is vanishingly small, and both agree well with experimental hardness data. We have also used the flow theory of MSG plasticity to investigate stress fields around a stationary mode-I crack tip as well as around a steady state, quasi-statically growing crack tip. At a distance to crack tip much larger than dislocation spacings such that continuum plasticity still applies, the stress level around a stationary crack tip in MSG plasticity is significantly higher than that in classical plasticity. The same conclusion is also established for a steady state, quasi-statically growing crack tip, though only the flow theory can be used because of unloading during crack propagation. This significant stress increase due to strain gradient effect provides a means to explain the experimentally observed cleavage fracture in ductile materials [J. Mater. Res. 9 (1994) 1734, Scripta Metall. Mater. 31 (1994) 1037; Interface Sci. 3(1996) 169].

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.

Mathematical Modeling of Near-Wall Flows of Two-Phase Mixture with Evaporating Droplets

In the framework of the two-continuum approach, using the matched asymptotic expansion method, the equations of a laminar boundary layer in mist flows with evaporating droplets were derived and solved. The similarity criteria controlling the mist flows were determined. For the flow along a curvilinear surface, the forms of the boundary layer equations differ from the regimes of presence and absence of the droplet inertia deposition. The numerical results were presented for the vapor-droplet boundary layer in the neighborhood of a stagnation point of a hot blunt body. It is demonstrated that, due to evaporation, a droplet-free region develops near the wall inside the boundary layer. On the upper edge of this region, the droplet radius tends to zero and the droplet number density becomes much higher than that in the free stream. The combined effect of the droplet evaporation and accumulation results in a significant enhancement of the heat transfer on the surface even for small mass concentration of the droplets in the free stream. 在双连续介质理论框架下,采用匹配渐进展开方法导出并求解了具有蒸发液滴的汽雾流中层流边界层方程,给出了控制汽雾流的相似判据。对于沿曲面的流动,边界层方程的形式取决于是否存在液滴的惯性沉积。给出了热钝体验点附近蒸汽。液滴边界层的数值计算结果。它们表明：由于蒸发,在边界层内近壁处形成了一个无液滴区域；在该区上边界处,液滴半径趋于零而液滴数密度急剧增高。液滴蒸发及聚集的联合效应造成了表面热流的显著增加,甚至在自由来流中液滴质量浓度很低时此效应依然存在。

Implication of scaling hierarchy associated with nonequilibrium: field and particulate

The advent of nanotechnology has necessitated a better understanding of how material microstructure changes at the atomic level would affect the macroscopic properties that control the performance. Such a challenge has uncovered many phenomena that were not previously understood and taken for granted. Among them are the basic foundation of dislocation theories which are now known to be inadequate. Simplifying assumptions invoked at the macroscale may not be applicable at the micro- and/or nanoscale. There are implications of scaling hierrachy associated with in-homegeneity and nonequilibrium. of physical systems. What is taken to be homogeneous and equilibrium at the macroscale may not be so when the physical size of the material is reduced to microns. These fundamental issues cannot be dispensed at will for the sake of convenience because they could alter the outcome of predictions. Even more unsatisfying is the lack of consistency in modeling physical systems. This could translate to the inability for identifying the relevant manufacturing parameters and rendering the end product unpractical because of high cost. Advanced composite and ceramic materials are cases in point. Discussed are potential pitfalls for applying models at both the atomic and continuum levels. No encouragement is made to unravel the truth of nature. Let it be partiuclates, a smooth continuum or a combination of both. The present trend of development in scaling tends to seek for different characteristic lengths of material microstructures with or without the influence of time effects. Much will be learned from atomistic simulation models to show how results could differ as boundary conditions and scales are changed. Quantum mechanics, continuum and cosmological models provide evidence that no general approach is in sight. Of immediate interest is perhaps the establishment of greater precision in terminology so as to better communicate results involving multiscale physical events.

Contingent sensitivity to configuration in a nonlinear evolution system

A trans-scopic sensitivity of macroscopic failure to slight differentiation in the meso-scopic structure of a system with nonlinear evolution is reported. A periodical chain following a non-local load-sharing evolution was applied as a propotype in failure study. The results demonstrate that there is a transition region composed of globally stable (GS) and evolution induced catastrophic (EIC) modes. That is different from a critical threshold as predicted by percolation and renormalization group theories. Moreover, the EIC mode shows a distinctive sample specific behaviour. For instance, some neighbouring initial states may evolve into completely different final states, though different initial states can evolve into the same final states. As an example, a marginal configuration of EIC mode, a quasi-Fibonacci skeleton, is constructed.

Plastic constitutive behavior of short-fiber/particle reinforced composites

A cylindrical cell model based on continuum theory for plastic constitutive behavior of short-fiber/particle reinforced composites is proposed. The composite is idealized as uniformly distributed periodic arrays of aligned cells, and each cell consists of a cylindrical inclusion surrounded by a plastically deforming matrix. In the analysis, the non-uniform deformation field of the cell is decomposed into the sum of the first order approximate field and the trial additional deformation field. The precise deformation field are determined based on the minimum strain energy principle. Systematic calculation results are presented for the influence of reinforcement volume fraction and shape on the overall mechanical behavior of the composites. The results are in good agreement with the existing finite element analyses and the experimental results. This paper attempts to stimulate the work to get the analytical constitutive relation of short-fiber/particle reinforced composites.