MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Molecular dynamics studies on dislocations in crystallites of nanocrystalline alpha-iron

Molecular dynamics simulations have been carried our to study the atomic structure of the crystalline component of nanocrystalline alpha-iron. A two-dimensional computational block is used to simulate the consolidation process. It is found that dislocations are generated in the crystallites during consolidation when the grain size is large enough. The critical value of the grain size for dislocation generation appears to be about 9 nm. This result agrees with experiment qualitatively. AN dislocations that are preset in the original grains glide out during consolidation. It shows that dislocations in the crystallites we generated in consolidation process, but not in the original grains. Higher consolidation pressure results in more dislocations. Furthermore, new interfaces are found within crystallites. These interfaces might result from the special environment of nanomaterial. (C) 1998 Acta Metallurgica Inc.

Evolution of thermoplastic shear localization and related microstructures in Awl/SiCp composites under dynamic compression

The localized shear deformation in the 2024 and 2124 Al matrix composites reinforced with SiC particles was investigated with a split Hopkinson pressure bar (SHPB) at a strain rate of about 2.0x10(3) s(-1). The results showed that the occurrence of localized shear deformation is sensitive to the size of SiC particles. It was found that the critical strain, at which the shear localization occurs, strongly depends on the size and volume fraction of SiC particles. The smaller the particle size, the lower the critical strain required for the shear localization. TEM examinations revealed that Al/SiCp interfaces are the main sources of dislocations. The dislocation density near the interface was found to be high and it decreases with the distance from the particles. The Al matrix in shear bands was highly deformed and severely elongated at low angle boundaries. The Al/SiCp interfaces, particularly the sharp corners of SiC particles, provide the sites for microcrack initiation. Eventual fracture is caused by the growth and coalescence of microcracks along the shear bands. It is proposed that the distortion free equiaxed grains with low dislocation density observed in the center of shear band result from recrystallization during dynamic deformation.

Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a = 10(-4) to 10(-2) where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions. (C) 2002 Published by Elsevier Science Ltd.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

Singularity of dynamic stress fields around an interface crack between viscoelastic bodies

The singular nature of the dynamic stress fields around an interface crack located between two dissimilar isotropic linearly viscoelastic bodies is studied. A harmonic load is imposed on the surfaces of the interface crack. The dynamic stress fields around the crack are obtained by solving a set of simultaneous singular integral equations in terms of the normal and tangent crack dislocation densities. The singularity of the dynamic stress fields near the crack tips is embodied in the fundamental solutions of the singular integral equations. The investigation of the fundamental solutions indicates that the singularity and oscillation indices of the stress fields are both dependent upon the material constants and the frequency of the harmonic load. This observation is different from the well-known -1/2 oscillating singularity for elastic bi-materials. The explanation for the differences between viscoelastic and elastic bi-materials can be given by the additional viscosity mismatch in the case of viscoelastic bi-materials. As an example, the standard linear solid model of a viscoelastic material is used. The effects of the frequency and the material constants (short-term modulus, long-term modulus and relaxation time) on the singularity and the oscillation indices are studied numerically.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.

Molecular dynamics simulation of plastic deformation of nanotwinned copper

The plastic deformation of polycrystalline Cu with ultrathin lamella twins has been studied using molecular dynamics simulations. The results of uniaxial tensile deformation simulation show that the abundance of twin boundaries provides obstacles to dislocation motion, which in consequence leads to a high strain hardening rate in the nanotwinned Cu. We also show that the twin lamellar spacing plays a vital role in controlling the strengthening effects, i.e., the thinner the thickness of the twin lamella, the harder the material. Additionally, twin boundaries can act as dislocation nucleation sites as they gradually lose coherency at large strain. These results indicate that controlled introduction of nanosized twins into metals can be an effective way of improving strength without suppression tensile ductility. (C) 2007 American Institute of Physics.

Variation of added mass and its application to the calculation of amplitude response for a circular cylinder

In the present study, analyzed are the variation of added mass for a circular cylinder in the lock-in ( synchronization) range of vortex-induced vibration (VIV) and the relationship between added mass and natural frequency. A theoretical minimum value of the added mass coefficient for a circular cylinder at lock-in is given. Developed are semi-empirical formulas for the added mass of a circular cylinder at lock-in as a function of flow speed and mass ratio. A comparison between experiments and numerical simulations shows that the semi-empirical formulas describing the variation of the added mass for a circular cylinder at lock-in are better than the ideal added mass. In addition, computation models such as the wake oscillator model using the present formulas can predict the amplitude response of a circular cylinder at lock-in more accurately than those using the ideal added mass.

Accommodation of large plastic strains and defect accumulation in nanocrystalline Ni grains

A transmission electron microscopy (TEM) study has been carried out to uncover how dislocations and twins accommodate large plastic strains and accumulate in very small nanocrystalline Ni grains during low-temperature deformation. We illustrate dislocation patterns that suggest preferential deformation and nonuniform defect storage inside the nanocrystalline grain. Dislocations are present in individual and dipole configurations. Most dislocations are of the 60 degrees type and pile up on (111) slip planes. Various deformation responses, in the forms of dislocations and twinning, may simultaneously occur inside a nanocrystalline grain. Evidence for twin boundary migration has been obtained. The rearrangement and organization of dislocations, sometimes interacting with the twins, lead to the formation of subgrain boundaries, subdividing the nanograin into mosaic domain structures. The observation of strain (deformation)-induced refinement contrasts with the recently reported stress-assisted grain growth in nanocrystalline metals and has implications for understanding the stability and deformation behavior of these highly nonequilibrium materials.

Atomistic simulations of crack nucleation and intergranular fracture in bulk nanocrystalline nickel

We report large scale molecular dynamics simulations of dynamic cyclic uniaxial tensile deformation of pure, fully dense nanocrystalline Ni, to reveal the crack initiation, and consequently intergranular fracture is the result of coalescence of nanovoids by breaking atomic bonds at grain boundaries and triple junctions. The results indicate that the brittle fracture behavior accounts for the transition from plastic deformation governed by dislocation to one that is grain-boundary dominant when the grain size reduces to the nanoscale. The grain-boundary mediated plasticity is also manifested by the new grain formation and growth induced by stress-assisted grain-boundary diffusion observed in this work. This work illustrates that grain-boundary decohesion is one of the fundamental deformation mechanisms in nanocrystalline Ni.

Deformation mechanisms of face-centered-cubic metal nanowires with twin boundaries

This letter addresses the issue of deformation mechanisms and mechanical tensile behavior of the twinned metal nanowires using atomistic simulations. Free surfaces are always the preferential dislocation nucleation sites in the initial inelastic deformation stage, while with further plastic deformation, twin boundary interfaces will act as sources of dislocations with the assistance of the newly formed defects. The smaller the twin boundary spacing, the higher the yielding stresses of the twinned nanowires. Twin boundaries, which serve both as obstacles to dislocation motion and dislocation sources, can lead to hardening effects and contribute to the tensile ductility. This work illustrates that the mechanical properties of metal nanowires could be controlled by tailoring internal growth twin structures. (c) 2007 American Institute of Physics.

Attractor crisis and bursting in a fluid flow with two no-slip directions

Characteristic burtsing behavior is observed in a driven, two-dimensional viscous flow, confined to a square domain and subject to no-slip boundaries. Passing a critical parameter value, an existing chaotic attractor undergoes a crisis, after which the flow initially enters a transient bursting regime. Bursting is caused by ejections from and return to a limited subdomain of the phase space, whereas the precrisis chaotic set forms the asymptotic attractor of the flow. For increasing values of the control parameter the length of the bursting regime increases progressively. Passing another critical parameter value, a second crisis leads to the appearance of a secondary type of bursting, of very large dynamical range. Within the bursting regime the flow then switches in irregular intervals from the primary to the secondary type of bursting. Peak enstrophy levels for both types of bursting are associated to the collapse of a primary vortex into a quadrupolar state.

On the topological bifurcation of flows around a rotating circular cylinder

Flow fields around a rotating circular cylinder in a uniform stream are computed using a low dimensional Galerkin method. Results show that the formation of a Fopple vortex pair behind a stationary circular cylinder is caused by the structural instability in the vicinity of the saddle located at the rear of the cylinder. For rotating cylinder a bifurcation diagram with the consideration of two parameters, Reynolds number Re and rotation parameter a, is built by a kinematic analysis of the steady flow fields.