Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

A characteristics study on the performance of a 2-stage light gas gun

In order to obtain an overall and systematic understanding of the performance of a two-stage light gas gun (TLGG), a numerical code to simulate the process occurring in a gun shot is advanced based on the quasi-one-dimensional unsteady equations of motion with the real gas effect,;friction and heat transfer taken into account in a characteristic formulation for both driver and propellant gas. Comparisons of projectile velocities and projectile pressures along the barrel with experimental results from JET (Joint European Tons) and with computational data got by the Lagrangian method indicate that this code can provide results with good accuracy over a wide range of gun geometry and loading conditions.

Two-dimensional simulation of inductive plasma sources with self-consistent power deposition

A time averaged two-dimensional fluid model including an electromagnetic module with self-consistent power deposition was developed to simulate the transport of a low pressure radio frequency inductively coupled plasma source. Comparsions with experiment and previous simulation results show, that the fluid model is feasible in a certain range of gas pressure. In addition, the effects of gas pressure and power input have been discussed.

Fractal features of oscillatory convection in the half-floating zone

Unsteady and two-dimensional numerical simulation is applied to study the transition process from steady convection to turbulence via subharmonic bifurcation in thermocapillary convection of a liquid bridge in the half-floating zone. The results of numerical tests show clearly the fractal structure of period-doubling bifurcations, and frequency-locking at f/4, f/8, f/16 with basic frequency f is observed with increasing temperature difference. The Feigenbaum universal constant is given by the present paper as delta(4) = 4.853, which can be compared with the theoretical value 4.6642016.

Molecular dynamics simulation of crack-tip processes in copper

The crack tip processes in copper under mode II loading have been simulated by a molecular dynamics method. The nucleation, emission, dislocation free zone (DFZ) and pile-up of the dislocations are analyzed by using a suitable atom lattice configuration and Finnis & Sinclair potential. The simulated results show that the dislocation emitted always exhibits a dissociated fashion. The stress intensity factor for dislocation nucleation, DFZ and dissociated width of partial dislocations are strongly dependent on the loading rate. The stress distributions are in agreement with the elasticity solution before the dislocation emission, but are not in agreement after the emission. The dislocation can move at subsonic wave speed (less than the shear wave speed) or at transonic speed (greater than the shear wave speed but less than the longitudinal wave speed), but at the longitudinal wave speed the atom lattice breaks down.

Numerical simulation of evolution-induced catastrophe

A numerical simulation of damage evolution in a two-dimensional system of micocracks is presented. It reveals that the failure is induced by a cascade of coalescences of microcracks, and the fracture surface appears fractal. A model of evolution-induced catastrophe is introduced. The fractal dimension is found to be a function of evolution rule only. This result could qualitatively explain the correlation of fractal dimension and fracture toughness discovered in experiments.

Numerical-simulation of ac plasma-arc thermodynamics

A mathematical model and approximate analysis for the energy distribution of an ac plasma arc with a moving boundary is developed. A simplified electrical conductivity function is assumed so that the dynamic behavior of the arc may be determined, independent of the gas type. The model leads to a reduced set of non-linear partial differential equations which governs the quasi-steady ac arc. This system is solved numerically and it is found that convection plays an important role, not only in the temperature distribution, but also in arc disruptions. Moreover, disruptions are found to be influenced by convection only for a limited frequency range. The results of the present studies are applicable to the frequency range of 10-10(2) Hz which includes most industry ac arc frequencies. (C) 1994 Academic Press, Inc.

Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

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.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

Energy conversions and storage caused by an unsteady poloidal flow in active solar regions

In this paper we discuss coupling processes between a magnetic field and an unsteady plasma motion, and analyze the features of energy storage and conversions in active region. It is pointed out that the static force-free field is insufficient for a discussion of storage processes, and also the pure unsteady plasma rotation is not a perfect approach. In order to analyze the energy storage, we must consider the addition of poloidal plasma motion. The paper shows that because the unsteady poloidal flow is added and coupling occurs between the magnetic field and both the toroidal and the poloidal plasma flows, an unsteady process is maintained which changes the force-free factor with time. Hence, the energy in the lower levels can be transferred to the upper levels, and a considerable energy can be stored in the active region. Finally, another storage process is given which is due to the pure poloidal flow. The article shows that even if there is no twisted magnetic line of force, the energy in the lower levels may still be transferred to the upper levels and stored there.