Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.

Influence of Surface Stress on Frequency of Microcantilever-Based Biosensors (vol 10, pg 307, 2004)

Microcantilever-based biosensors have been found increasing applications in physical, chemical, and biological fields in recent years. When biosensors are used in those fields, surface stress and mass variations due to bio-molecular binding can cause the microcantilever deform or the shift of frequency. These simple biosensors allow biologists to study surface biochemistry on a micro or nano scale and offer new opportunities in developing microscopic biomedical analysis with unique characteristics. To compare and illustrate the influence of the surface stress on the frequency and avoid unnecessary and complicated numerical solution of the resonance frequency, some dimensionless numbers are derived in this paper by making governing equations dimensionless. Meanwhile, in order to analyze the influence of the general surface stress on the frequency, a new model is put forward, and the frequency of the microcantilever is calculated by using the subspace iteration method and the Rayleigh method. The sensitivity of microcantilever is also discussed. (19 refs.)

Piezoelectricity Of Zno Films Prepared By Sol-Gel Method

ZnO piezoelectric thin films were prepared on crystal substrate Si(111) by sol-gel technology, then characterized by scanning electron microscopy, X-ray diffraction and atomic force microscopy (AFM). The ZnO films characterized by X-ray diffraction are highly oriented in (002) direction with the growing of the film thickness. The morphologies, roughness and grain size of ZnO film investigated by AFM show that roughness and grain size of ZnO piezoelectric films decrease with the increase of the film thickness. The roughness dimension is 2.188-0.914 nm. The piezoelectric coefficient d(33) was investigated with a piezo-response force microscope (PFM). The results show that the piezoelectric coefficient increases with the increase of thickness and (002) orientation. When the force reference is close to surface roughness of the films, the piezoelectric coefficient measured is inaccurate and fluctuates in a large range, but when the force reference is big, the piezoelectric coefficient d(33) changes little and ultimately keeps constant at a low frequency.

An Improved Energy Method For Determining Young'S Modulus By Instrumented Indentation Using A Berkovich Tip

We previously proposed a method for estimating Young's modulus from instrumented nanoindentation data based on a model assuming that the indenter had a spherical-capped Berkovich geometry to take account of the bluntness effect. The method is now further improved by releasing the constraint on the tip shape, allowing it to have a much broader arbitrariness to range from a conical-tipped shape to a flat-ended shape, whereas the spherical-capped shape is just a special case in between. This method requires two parameters to specify a tip geometry, namely, a volume bluntness ratio V-r and a height bluntness ratio h(r). A set of functional relationships correlating nominal hardness/reduced elastic modulus ratio (H-n/E-r) and elastic work/total work ratio (W-e/W) were established based on dimensional analysis and finite element simulations, with each relationship specified by a set of V-r and h(r). Young's modulus of an indented material can be estimated from these relationships. The method was shown to be valid when applied to S45C carbon steel and 6061 aluminum alloy.

The Newest Developments Of Load-Unload Response Ratio (Lurr)

The Load Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It is inspiring that its predictions using LURR have been improving. Since 2004 we have made a major breakthrough in intermediate-term earthquake forecasting of the strong earthquakes on the Chinese mainland using LURR and successfully predicted the Pakistan earthquake with magnitude M 7.6 on October 8, 2005. The causes for improving the prediction in terms of LURR have been discussed in the present paper.

Statistical Interpolation Method of Turbulent Phase Screen

A relative displacement between the grid points of optical fields and those of phase screens may occur in the simulation of light propagation through the turbulent atmosphere. A statistical interpolator is proposed to solve this problem in this paper. It is evaluated by the phase structure function and numerical experiments of light propagation through atmospheric turbulence with/without adaptive optics (AO) and it is also compared with the well-known linear interpolator under the same condition. Results of the phase structure function show that the statistical interpolator is more accurate in comparison with the linear one, especially in the high frequency region. More importantly, the long-exposure results of light propagation through the turbulent atmosphere with/without AO also show that the statistical interpolator is more accurate and reliable than the linear one. (C) 2009 Optical Society of America.

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.

Self-consistent kinetic modeling of low-pressure inductively coupled radio-frequency discharges

An efficient method for solving the spatially inhomogeneous Boltzmann equation in a two-term approximation for low-pressure inductively coupled plasmas has been developed. The electron distribution function (EDF), a function of total electron energy and two spatial coordinates, is found self-consistently with the static space-charge potential which is computed from a 2D fluid model, and the rf electric field profile which is calculated from the Maxwell equations. The EDF and the spatial distributions of the electron density, potential, temperature, ionization rate, and the inductive electric field are calculated and discussed. (C) 1996 American Institute of Physics.

A BIEM optimization method for fracture dynamics inverse problem

In the present paper, based on the theory of dynamic boundary integral equation, an optimization method for crack identification is set up in the Laplace frequency space, where the direct problem is solved by the author's new type boundary integral equations and a method for choosing the high sensitive frequency region is proposed. The results show that the method proposed is successful in using the information of boundary elastic wave and overcoming the ill-posed difficulties on solution, and helpful to improve the identification precision.

A New Model for Sediment Transport

The problem of predicting sediment transportation by water waves is treated analytically with the rate of wave energy dissipation or wave damping. With resorting to the theory of shallow water waves and the basis of Yamamoto’s Coulomb-damped poroelastic model, the Boussinesq-type equation has been derived over a variation depth bed. For convenience Cnoidal wave is just discussed, The Cnoidal wave with complex wave length and wave velocity, which are as a function of wave frequency, water depth, permeability, Poisson’s ratio and complex elastic moduli of bed soil, is applied to analyse the rate of sediment transportation. Considering the sediment transportation depended on the shear stress near-bed or the horizontal velocity, the conclusion of Yamamoto’s experiment in clay bed has been extended to general situation. It could be figured out that the model should provide a method to avoid the undistinguishable factors during sediment transport processes and relate mass transport with the sediment peculiarities.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Piezoelectricity of ZnO films prepared by sol-gel method

ZnO piezoelectric thin films were prepared on crystal substrate Si(111) by sol-gel technology, then characterized by scanning electron microscopy, X-ray diffraction and atomic force microscopy (AFM). The ZnO films characterized by X-ray diffraction are highly oriented in (002) direction with the growing of the film thickness. The morphologies, roughness and grain size of ZnO film investigated by AFM show that roughness and grain size of ZnO piezoelectric films decrease with the increase of the film thickness. The roughness dimension is 2.188-0.914 nm. The piezoelectric coefficient d(33) was investigated with a piezo-response force microscope (PFM). The results show that the piezoelectric coefficient increases with the increase of thickness and (002) orientation. When the force reference is close to surface roughness of the films, the piezoelectric coefficient measured is inaccurate and fluctuates in a large range, but when the force reference is big, the piezoelectric coefficient d(33) changes little and ultimately keeps constant at a low frequency.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

An Instrumented Indentation Method for Young's Modulus Measurement with Accuracy Estimation

The relationships between indentation responses and Young's modulus of an indented material were investigated by employing dimensional analysis and finite element method. Three representative tip bluntness geometries were introduced to describe the shape of a real Berkovich indenter. It was demonstrated that for each of these bluntness geometries, a set of approximate indentation relationships correlating the ratio of nominal hardness/reduced Young's modulus H (n) /E (r) and the ratio of elastic work/total work W (e)/W can be derived. Consequently, a method for Young's modulus measurement combined with its accuracy estimation was established on basis of these relationships. The effectiveness of this approach was verified by performing nanoindentation tests on S45C carbon steel and 6061 aluminum alloy and microindentation tests on aluminum single crystal, GCr15 bearing steel and fused silica.