Molecular dynamics simulation of the role of dislocations in microcrack healing

The molecular dynamics method is used to simulate microcrack healing during heating or/and under compressive stress. A centre microcrack in Cu crystal would be sealed under compressive stress or by heating. The role of compressive stress and heating in crack healing was additive. During microcrack healing, dislocation generation and motion occurred. When there were pre-existing dislocations around the microcrack, the critical temperature or compressive stress necessary for microcrack healing would decrease, and, the higher the number of dislocations, the lower the critical temperature or compressive stress. The critical temperature necessary for microcrack healing depended upon the orientation of the crack plane. For example, the critical temperature for the crack along the (001) plane was the lowest, i.e. 770K.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Structural Failure Analysis and Numerical Simulation of Micro-Accelerometers under Impulsive Loading

Micromachined accelerometer is a kind of inertial MEMS devices, which usually operate under intensive impact loading. The reliability of micromachined accelerometers is one of the most important performance indices for their design, manufacture and commer

Burst Detection in Turbulent Channel Flows Based on Large Eddy Simulation Databases

Reliable turbulent channel flow databases at several Reynolds numbers have been established by large eddy simulation (LES), with two of them validated by comparing with typical direct numerical simulation (DNS) results. Furthermore, the statistics, such as velocity profile, turbulent intensities and shear stress, were obtained as well as the temporal and spatial structure of turbulent bursts. Based on the LES databases available, the conditional sampling methods are used to detect the structures of burst events. A method to deterimine the grouping parameter from the probability distribution function (pdf) curve of the time separation between ejection events is proposed to avoid the errors in detected results. And thus, the dependence of average burst period on thresholds is considerably weakened. Meanwhile, the average burst-to-bed area ratios are detected. It is found that the Reynolds number exhibits little effect on the burst period and burst-to-bed area ratio.

Monte-Carlo simulation of surface crack growth rate for offshore structural steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.

Two-dimensional simulation of an electron cyclotron resonance plasma source with self-consistent power deposition

Using a refined two-dimensional hybrid-model with self-consistent microwave absorption, we have investigated the change of plasma parameters such as plasma density and ionization rate with the operating conditions. The dependence of the ion current density and ion energy and angle distribution function at the substrate surface vs. the radial position, pressure and microwave power were discussed. Results of our simulation can be compared qualitatively with many experimental measurements.

Intensity Distribution Design for Laser-Induced Thermal Loading Based on Numerical Simulation

To accomplish laser-induced thermal loading simulation tests for pistons,the Gaussian beam was modulated into multi-circular beam with specific intensity distribution.A reverse method was proposed to design the intensity distribution for the laser-induced thermal loading based on finite element(FE) analysis.Firstly,the FE model is improved by alternating parameters of boundary conditions and thermal-physical properties of piston material in a reasonable range,therefore it can simulate the experimental resul...

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.

Numerical simulation of an adaptive optics system with laser propagation in the atmosphere

A comprehensive model of laser propagation in the atmosphere with a complete adaptive optics (AO) system for phase compensation is presented, and a corresponding computer program is compiled. A direct wave-front gradient control method is used to reconstruct the wave-front phase. With the long-exposure Strehl ratio as the evaluation parameter, a numerical simulation of an AO system in a stationary state with the atmospheric propagation of a laser beam was conducted. It was found that for certain conditions the phase screen that describes turbulence in the atmosphere might not be isotropic. Numerical experiments show that the computational results in imaging of lenses by means of the fast Fourier transform (FFT) method agree well with those computed by means of an integration method. However, the computer time required for the FFT method is 1 order of magnitude less than that of the integration method. Phase tailoring of the calculated phase is presented as a means to solve the problem that variance of the calculated residual phase does not correspond to the correction effectiveness of an AO system. It is found for the first time to our knowledge that for a constant delay time of an AO system, when the lateral wind speed exceeds a threshold, the compensation effectiveness of an AO system is better than that of complete phase conjugation. This finding indicates that the better compensation capability of an AO system does not mean better correction effectiveness. (C) 2000 Optical Society of America.

Experimental study and simulation principles of an oil-gas multiphase transportation system

Presented is an experimental study on the performance of an oil-gas multiphase transportation system, especially on the multiphase flow patterns, multiphase pumping and multiphase metering of the system. A dynamic simulation analysis is conducted to deduce simulation parameters of the system and similarity criteria under simplified conditions are obtained. The reliability and feasibility of two-phase flow experiment with oil and natural gas simulated by water and air are discussed by using the similarity criteria.

A Generalized Soft-Sphere for Monte Carlo Simulation

A new collision model, called the generalized soft-sphere (GSS) model, is introduced. It has the same total cross section as the generalized hard-sphere model [Phys. Fluids A 5, 738 (1993)], whereas the deflection angle is calculated by the soft-sphere scattering model [Phys. Fluids A 3, 2459 (1991)]. In virtue of a two-term formula given to fit the numerical solutions of the collision integrals for the Lennard-Jones (6-12) potential and for the Stockmayer potential, the parameters involved in the GSS model are determined explicitly that may fully reproduce the transport coefficients under these potentials. Coefficients of viscosity, self-diffusion and diffusion for both polar and nonpolar molecules given by the GSS model and experiment are in excellent agreement over a wide range of temperature from low to high.

Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Numerical Simulation of Nox Formation in Coal Combustion with Inlet Natural Gas Burning

A full two-fluid model of reacting gas-particle flows and coal combustion is used to simulate coal combustion with and without inlet natural gas added in the inlet. The simulation results for the case without natural gas burning is in fair agreement with the experimental results reported in references. The simulation results of different natural gas adding positions indicate that the natural gas burning can form lean oxygen combustion enviroment at the combustor inlet region and the NOz concentration is reduced. The same result can be obtained from chemical equilibrium analysis.