Splitting the Fast and Slow Motions in Molecular Dynamics Simulations Based on the Change of Cold Potential Well Bottom

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.

Thermodynamic model of helium and hydrogen co-implanted silicon surface layer splitting

A thermodynamic model of the evolution of microcracks in silicon caused by helium and hydrogen co-implantation during annealing was studied. The crack growth rate relies on the amount of helium atoms and hydrogen molecules present. Here, the crack radius was studied as a function of annealing time and temperature, and compared with experimental results. The mean crack radius was found to be proportional to the annealing temperature and the helium and hydrogen implanted fluence. The gas desorption should be considered during annealing process. (C) 2009 Elsevier B.V. All rights reserved.

Autler-townes splitting in the multiphoton resonance ionization spectrum of molecules produced by ultrashort laser pulses

We report for the first time the proper conditions to observe Autler-Townes splitting (ac-Stark splitting) from vibrationally coherent states belonging to the different electronic terms of a diatomic molecule. Wave packet dynamics simulations demonstrate that such a process is feasible by multiphoton resonance ionization of the molecule Na-2 with a single ultrashort intense laser pulse. With the ultrahigh time resolution of a femtosecond laser pulse, one can directly measure the absolute value of the transition dipole moment between any kinds of molecular states by this kind of Autler-Townes splitting, which is a function of the internuclear distance R.

Effect of N fertilizer on the productivity and nutritive values of perennial grass mixtures in the alpine region of Qinghai-Tibetan Plateau, China

A study was conducted on grass mixtures that included smooth bromegrass (SB) + drooping wild ryegrass (DW), smooth bromegrass + Siberian wild ryegrass (SW) + crested wheatgrass (CW) and smooth bromegrass + Siberian wild ryegrass + drooping wild ryegrass + crested wheatgrass in the alpine region of Qinghai-Tibetan Plateau. The study was conducted from 1998 to 2000 to investigate the effects of N application rates and growing year on herbage dry matter (DM) yield and nutritive values. Herbage DM production increased linearly with N application rates. The effect of N application on DM yields was greater (P < 0.05) in the 2nd and 3rd production years than in the establishment year. Dry matter yields of SB + SW + CW and SB + SW + DW + CW can reach as high as 15 000 kg ha(-1) at 345 kg ha(-1) N rate in the 3rd growing year. With increased N application rates, crude protein (CP) contents and 48 h in sacco DM degradability of grasses increased (P < 0.05). No effect (P > 0.05) of N application was detected on organic matter (OM) and acid detergent fibre (ADF) concentration. It can be concluded that for increased biomass production in the alpine region of the Qinghai-Tibetan Plateau, a minimum of 345 kg N ha(-1) should be applied to grass stands in three split application of 115 kg N ha(-1), in early June, early July and late July

Nanocrystalline Formation During Mechanical Attrition of Cobalt

The nanocrystalline (nc) formation was studied in cobalt (a mixture of c (hexagonal close packed) and gamma (face-centered cubic) phases) subjected to surface mechanical attrition treatment. Electron microscopy revealed the operation of {10(1) over bar 0}< 11(2) over bar 0 > prismatic and {0001}< 11(2) over bar 0 > basal slip in the E phase, leading to the successive subdivision of grains to nanoscale. In particular, the dislocation splitting into the stacking faults was observed to occur in ultrafine and nc grains. By contrast, the planar dislocation arrays, twins and martensites were evidenced in the gamma phase. The strain-induced gamma ->epsilon martensitic transformation was found to progress continuously in ultrafine and nc grains as the strain increased. The nc formation in the gamma phase was interpreted in terms of the martensitic transformation and twinning.

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed. to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number.

Studies on the phase behavior of lyotropic liquid crystal in DAD/C4OH/n-C8H18/D2O system

The lytropic liquid crystals in dodecanic acid diethanolamine (DAD)/n-butanol (C4OH)/octane (n-C8H18)/deuteron (D2O) system were studied to determine the phase regions and were investigated by H-2-NMR spectroscopy,optical polarizing microscope and small-angle X-ray diffraction (SAXD) methods. The results indicate that the lamellar, hexagonal and cubic liquid crystals all exist in the above system. Keeping the weight ratio of DAD and C4OH constant,the microphase structure, H-2 quadruple splitting and the interlayer spacing are all changed with the addition of deuteron.

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

Numerical simulation of hypersonic viscous flow around space shuttle with method of diffusion analogy

Hypersonic viscous flow around a space shuttle with M(infinity) = 7, Re = 148000 and angle of attack alpha = 5-degrees is simulated numerically with the special Jacobian matrix splitting technique and simplified diffusion analogy method. With the simplified diffusion analogy method the efficiency of computation and resolution of the shock can be improved.

Viscosity-induced mode splitting and potential energy criterion for mode stability

This paper points out that viscosity can induce mode splitting in a uniform infinite cylinder of an incompressible fluid with self-gravitation, and that the potential energy criterion cannot be appropriate to all normal modes obtained, i.e., there will be stable modes with negative potential energy (<0). Therefore the condition >0 is not necessary, although sufficient, for the stability of a mode in an incompressible static fluid or magnetohydrodynamics (MHD) system, which is a correction of both Hare's [Philos. Mag. 8, 1305 (1959)] and Chandrasekhar's [Hydrodynamic and Hydromagnetic Stability (Oxford U.P., Oxford, 1961), p. 604] stability criterion for a mode. These results can also be extended to compressible systems with a polytropic exponent.

Nanocrystalline Formation during Mechanical Attrition of Cobalt

The nanocrystalline (nc) formation was studied in cobalt (a mixture of c (hexagonal close packed) and gamma (face-centered cubic) phases) subjected to surface mechanical attrition treatment. Electron microscopy revealed the operation of {10(1) over bar 0}< 11(2) over bar 0 > prismatic and {0001}< 11(2) over bar 0 > basal slip in the E phase, leading to the successive subdivision of grains to nanoscale. In particular, the dislocation splitting into the stacking faults was observed to occur in ultrafine and nc grains. By contrast, the planar dislocation arrays, twins and martensites were evidenced in the gamma phase. The strain-induced gamma ->epsilon martensitic transformation was found to progress continuously in ultrafine and nc grains as the strain increased. The nc formation in the gamma phase was interpreted in terms of the martensitic transformation and twinning.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.