旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity

Nano-hardness Techniques and Its Applications in Mechanical Property Measurements for Surface Engineering

结合纳米硬度技术测量各类薄膜和块体材料表层的纳米压痕硬度、弹性模量、断裂韧性、膜厚、微结构的弯曲变形,采用纳米划痕硬度技术测量各类薄膜和块体材料的粗糙度、临界附着力、摩擦系数、划痕横剖面.纳米硬度计是检测材料表层微米乃至几十纳米力学性能的先进仪器,可广泛应用于表面工程中的质量检测.

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

Numerical verification for instrumented spherical indentation techniques in determining the plastic properties of materials

 Instrumented indentation tests have been widely adopted for elastic modulus determination. Recently, a number of indentation-based methods for plastic properties characterization have been proposed, and rigorous verification is absolutely necessary for their wide application. In view of the advantages of spherical indentation compared with conical indentation in determining plastic proper-ties, this study mainly concerns verification of spherical indentation methods. Five convenient and simple models were selected for this purpose, and numerical experiments for a wide range of materials are carried out to identify their accuracy and sensitivity characteristics. The verification results show that four of these five methods can give relatively accurate and stable results within a certain material domain, which is defined as their validity range and has been summarized for each method.

Singularity criteria for perturbation series

This paper is concerned with some mathematical aspects of the Van Dyke method inperturbation theory, i.e. the singularity criteria of perturbation series. The author suggestsa sign criterion and a Domb-syke plot for the cases with complex conjugate singularities, thussucceeding in extending the conclusions of Van Dyke's. Subsequently. effects of singularitiesof the lower order upon the criteria are taken into account. In addition, a method of locat-ing singular points is developed by analysing the new perturbation series derived by the Eulertransformation.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Improved sampling techniques for the direct simulation Monte Carlo method

The direct simulation Monte Carlo (DSMC) method is a widely used approach for flow simulations having rarefied or nonequilibrium effects. It involves heavily to sample instantaneous values from prescribed distributions using random numbers. In this note, we briefly review the sampling techniques typically employed in the DSMC method and present two techniques to speedup related sampling processes. One technique is very efficient for sampling geometric locations of new particles and the other is useful for the Larsen-Borgnakke energy distribution.

Phi 140 mm sapphire crystal growth by temperature gradient techniques and its color centers

Sapphire crystals, 140 mm in diameter and 90 turn in height, have been grown by temperature gradient techniques (TGT). The growth direction of the boule was fixed by means of Lane X-ray diffraction. A prominent 204 nm absorption band in TGT-Al2O3. which does not appear in single crystals grown by Czochralski method has been studied. Analysis further substantiates the F-center model of this band. Two relatively weaker bands absorbing at 232 nm and 254 nm were ascribed to F+ centers. F-type centers concentration was determined using Smakula's equation. (c) 2005 Elsevier B.V. All rights reserved.

Color centers and charge state change in Ce : YAG crystals grown by temperature gradient techniques

Color centers and impurity defects of Ce:YAG crystals grown in reduction atmosphere by temperature gradient techniques have been investigated by means of gamma irradiation and thermal treatments. Four absorption bands associated with color centers or impurity defects at 235, 255, 294 and 370 nm were observed in as-grown crystals. Changes in optical intensity of the 235 and 370 nm bands after gamma irradiation indicate that they are associated with F+-type color center. Charge state change processes of Fe3+ impurity and Ce3+ ions take place in the irradiation process. The variations of Ce3+ ions concentration clearly indicate that Ce4+ ions exist in Ce:YAG crystals and gamma irradiations could increase the concentration of Ce3+ ions. Annealing treatments and the changes in optical density suggest that a heterovalent impurity ion associated with the 294 nm band seems to be present in the crystals. (c) 2005 Elsevier B.V. All rights reserved.

Color centers in Yb : YAG crystals grown by temperature-gradient techniques

Yb: YAG (Yb: Y3Al5O12) crystals have been grown by temperature-gradient techniques (TGT) and their color centers and impurity defects were investigated by means of gamma irradiations and thermal treatment. Two color centers located at 255 and 290 nm were observed in the as-grown TGT-Yb: YAG. Analysis shows that the 255 nm band may be associated with Fe3+ ions. Absorption intensity changes of the 290 nm band after gamma irradiation and thermal treatment indicate that this band may be associated with oxygen-vacancy defects. (c) 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Influence of oxygen plasma treatment on properties of ZrO2 films prepared by e-beam evaporation techniques

Zirconia films were prepared by e-beam evaporation, and oxygen plasma treatment was used to modify film properties. Spectrophotometry, x-ray diffractometry (XRD), and atomic force microscopy were used to characterize refractive index, extinction coefficient, rnicrostructure, and surface roughness, respectively. The experimental results indicate that both refractive index and extinction coefficient of the films were reduced slightly after oxygen plasma treatment, with the decrease of intrinsic stress and surface roughness. From XRD spectra, the intensity decrease of the T(110) diffraction peak was clearly observed after the treatment, which was caused by the restructuring of the film atoms. (C) 2008 Elsevier Ltd. All rights reserved.

Techniques for fabrication of large-size graded reflectivity mirrors

According to the parameter requirements of a graded reflectivity mirror with a Gaussian profile, the layer structure and the mask pattern are designed using a graded-thickness middle layer. The mask and the automatic mask-switchover equipment are designed considering the actual requirement of the thin films and the specific deposit facility. The uniformity of the layer thickness is analyzed. The measurement results indicate that samples prepared with this technique are basically in accordance with the design parameter. The scattering effect between the material molecules and the mask, thickness errors, and the alignment error between the mask and the substrate are the main factors that influence the deposit result. (c) 2008 Society of Photo-Optical Instrumentation Engineers. [DOI: 10.1117/1.3027595]