Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

Refractive index and thickness analysis of natural silicon dioxide film growing on silicon with variable-angle spectroscopic ellipsometry

The refractive index and thickness of SiO2 thin films naturally grown on Si substrates were determined simultaneously within the wavelength range of 220-1100 nm with variable-angle spectroscopic ellipsometry. Different angles of incidence and wavelength ranges were chosen to enhance the analysis sensitivity for more accurate results. Several optical models describing the practical SiO2-Si system were investigated, and best results were obtained with the optical model, including an interface layer between SiO2 and Si, which proved the existence of the interface layer in this work as described in other publications.

Effects of Molecular Length, Ortientation and Amino Acid Mutation on Kinetics of Selectin/Ligand interactions

Receptor/ligand interactions are basic issues to cell adhesion, which are important to many physiological and pathological processes such as lymphocyte-mediated cytotoxicity, tumor metastasis and inflammatory reactionl. Selectin/carbohydrate ligand bindings have been found to mediate the fast rolling of leukocytes on activated endothelial monolayer. Kinetic rate and binding affinity constants are essential determinants of cell adhesion...

Effects Of The Length Of A Cylindrical Solid Shield On The Entrainment Of Ambient Air Into Turbulent And Laminar Impinging Argon Plasma Jets

When materials processing is conducted in air surroundings by use of an impinging plasma jet, the ambient air will be entrained into the materials processing region, resulting in unfavorable oxidation of the feedstock metal particles injected into the plasma jet and of metallic substrate material. Using a cylindrical solid shield may avoid the air entrainment if the shield length is suitably selected and this approach has the merit that expensive vacuum chamber and its pumping system are not needed. Modeling study is thus conducted to reveal how the length of the cylindrical solid shield affects the ambient air entrainment when materials processing (spraying, remelting, hardening, etc.) is conducted by use of a turbulent or laminar argon plasma jet impinging normally upon a flat substrate in atmospheric air. It is shown that the mass flow rate of the ambient air entrained into the impinging plasma jet cannot be appreciably reduced unless the cylindrical shield is long enough. In order to completely avoid the air entrainment, the gap between the downstream-end section of the cylindrical solid shield and the substrate surface must be carefully selected, and the suitable size of the gap for the turbulent plasma jet is appreciably larger than that for the laminar one. The overheating of the solid shield or the substrate could become a problem for the turbulent case, and thus additional cooling measure may be needed when the entrainment of ambient air into the turbulent impinging plasma jet is to be completely avoided.

Effects of physical parameter range on dimensionless variable sensitivity in water flooding reservoirs

The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.

Effects of physical parameter range on dimensionless variable sensitivity in water flooding reservoirs

The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.

Variable power singular interface elements for a crack normal to the bimaterial interface

Zero thickness crack tip interface elements for a crack normal to the interface between two materials are presented. The elements are shown to have the desired r(lambda-1) (0 < lambda < 1) singularity in the stress field at the crack tip and are compatible with other singular elements. The stiffness matrices of the quadratic and cubic interface element are derived. Numerical examples are given to demonstrate the applicability of the proposed interface elements for a crack perpendicular to the bimaterial interface.

3-dimensional transient wave response in a cracked elastic solid

The three-dimensional transient wave response problem is presented for an infinite elastic medium weakened by a plane crack of infinite length and finite width. Tractions are applied suddenly to the crack, which simulates the case of impact loading. The integral transforms are utilized to reduce the problem to a standard Fredholm integral equation in the Laplace transform variable and sequentially invert the Laplace transforms of the stress components by numerical inversion method. The dynamic mode I stress intensity factors at the crack tip are obtained and some numerical results are presented in graphical form.

The effect of variable currents on internal solitary waves

The effect of variable currents on internal solitary waves is described within the context of a variable coefficient Korteweg-de Vries (KdV) equation, and the approximate slowly varying, solitary-wave solution of this equation. The general theory which leads to the variable coefficient KdV equation is described; a derivation for the special case when the solitary wave and the current are aligned in the same direction is given in the Appendix. Using further simplifications and approximations, a number of analytical expressions are obtained for the variation in the solitary wave amplitude resulting from variable shear in the basic current or from when the basic current is a depth-independent flow which is a simple representation of a geostrophic current, tidal flow or inertial wave.

Statistical modeling of damage evolution in spallation

In order to understand the mechanism of the incipient spallation in rolled metals, a one dimensional statistical mode1 on evolution of microcracks in spallation was proposed. The crack length appears to be the fundamental variable in the statistical description. Two dynamic processes, crack nucleation and growth, were involved in the model of damage evolution. A simplified case was examined and preliminary correlation to experimental observations of spallation was made.

Distributions of temperature, velocity and viscosity at symmetrical axis of cylindrical plume in earth's mantle with variable viscosity

The parameters at the symmetrical axis of a cylindrical plume characterize the strength of this plume and provide a boundary condition which must be given to investigate the structure of a plume. For Newtonian fluid with a temperature-and pressure-dependence viscosity, an asymptotical solution of hydrodynamic equations at the symmetrical axis of the plume is found in the present paper. The temperature, upward velocity and viscosity at the symmetrical axis have been obtained as functions of depth, The calculated results have been given for two typical sets of Newtonian rheological parameters. The results obtained show that the temperature distribution along the symmetrical axis is nearly independent of the theological parameters. The upward velocity at the symmetrical axis, however, is strongly dependent on the rheological parameters.

Discussions on slip length measurements by microPIV/PTV in microchannels

Fluid transportation in microfluidic system could be benefit from the slip on solid-liquid interface. Slip length on many kinds of hydrophilic/hydrophobic surfaces have been measured recently. The two common-used experimental methods for boundary slip measurement include: (1) surface force measurement, such as surface force apparatus (SFA), atom force microscope (AFM), and (2) velocity measurement, like microPIV/PTV (Particle image velocimetry / Particle tracking velocimetry), total internal reflection velocimetry (TIRV). However, the measured results are rather scattered, larger measured slip lengths were reported by microPIV/PTV experiments. In this paper, we will investigate the deviations of the measured slip length on smooth hydrophilic surface. After measuring detailed velocity profiles very close to hydrophilic glass wall, we give a discussion on the effects influencing the slip measurements.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

A block SPP algorithm for multidimensional tridiagonal equations with optimal message vector length

A parallel strategy for solving multidimensional tridiagonal equations is investigated in this paper. We present in detail an improved version of single parallel partition (SPP) algorithm in conjunction with message vectorization, which aggregates several communication messages into one to reduce the communication cost. We show the resulting block SPP can achieve good speedup for a wide range of message vector length (MVL), especially when the number of grid points in the divided direction is large. Instead of only using the largest possible MVL, we adopt numerical tests and modeling analysis to determine an optimal MVL so that significant improvement in speedup can be obtained.

On optimal message vector length for block single parallel partition algorithm in a three-dimensional ADI solver

It has long been recognized that many direct parallel tridiagonal solvers are only efficient for solving a single tridiagonal equation of large sizes, and they become inefficient when naively used in a three-dimensional ADI solver. In order to improve the parallel efficiency of an ADI solver using a direct parallel solver, we implement the single parallel partition (SPP) algorithm in conjunction with message vectorization, which aggregates several communication messages into one to reduce the communication costs. The measured performances show that the longest allowable message vector length (MVL) is not necessarily the best choice. To understand this observation and optimize the performance, we propose an improved model that takes the cache effect into consideration. The optimal MVL for achieving the best performance is shown to depend on number of processors and grid sizes. Similar dependence of the optimal MVL is also found for the popular block pipelined method.