Fast calculation technique for scattering in T-matrix method

In scattering calculations using the T-matrix method, the calculation of the T-matrix involves multiplication and inversion of matrices. These two types of matrix operations are time-consuming, especially for the matrices with large size. Petrov et al. [D. Petrov, Y. Shkuratov, G. Videen, Opt. Lett. 32 (2007) 1168] proposed an optimized matrix inversion technique, which suggests the inversion of two matrices, each of which contains half the number of rows. This technique reduces time-consumption significantly. On the basis of this approach, we propose another fast calculation technique for scattering in the T-matrix method, which obtains the scattered fields through carrying out only the operations between matrices and the incident field coefficient. Numerical results show that this technique can decrease time-consumption by more than half that of the optimized matrix inversion technique by Petrov et al. (c) 2008 Elsevier B.V. All rights reserved.

Method for calculating ecological water storage and ecological water requirement of marsh

As one of the most typical wetlands, marsh plays an important role in hydrological and economic aspects, especially in keeping biological diversity. In this study, the definition and connotation of the ecological water storage of marsh is discussed for the first time, and its distinction and relationship with ecological water requirement are also analyzed. Furthermore, the gist and method of calculating ecological water storage and ecological water requirement have been provided, and Momoge wetland has been given as an example of calculation of the two variables. Ecological water use of marsh can be ascertained according to ecological water storage and ecological water requirement. For reasonably spatial and temporal variation of water storage and rational water resources planning, the suitable quantity of water supply to marsh can be calculated according to the hydrological conditions, ecological demand and actual water resources.

Calculation of the thermal expansion coefficient for Bi2Sr2CaCu2O8

An estimation method of thermal expansion coefficient in term of lattice energy which was developed earlier for simple materials is extended to a complex material of Bi2Sr2CaCu2O8 (Bi-2212). The calculation of the chemical bond property and thermal expansion coefficient of Bi-2212 has been carried out and the theoretical values were in good agreement with the corresponding experimental results. The dependence of the thermal expansion coefficient on the different structures and on the flexible oxidation states of Bi and Cu are investigated. The results indicate that the thermal expansion coefficients of Bi-2212 are insensitive to the low lattice distortion of the average structure and the changes of formal valences of Bi and Cu ions.

Numerical approximation improves well survey calculation

It is now possible to improve the precision of well survey calculations by order of magnitude with numerical approximation.

Although the most precise method of simulating and calculating a wellbore trajectory generally requires more calculation than other, less-accurate methods, the wider use of computers in oil fields now eliminates this as an obstacle.

The results of various calculations show that there is a deviation of more than 10 m among the different methods of calculation for a directional well of 3,000 m.¹ Consequently, it is important to improve the precision and reliability of survey calculation-the fundamental, necessary work of quantitatively monitoring and controlling wellbore trajectories.

Thermal Strain Analysis of a Flip-Chip Package by a Modified Hybrid Method

介绍一种可用于微电子封装局部应变场分析的实验/计算混合方法,该方法结合了有限元的整体/局部模型和实时的激光云纹干涉技术,利用激光云纹干涉技术所测得的应变场来校核有限元整体模型的计算结果,并用整体模型的结果作为局部模型的边界条件,对实验难以确定的封装结构局部位置的应力、应变场进行分析.用这种方法对可控坍塌倒装封装结构在热载荷作用下焊球内的应变场分布进行了分析,结果表明该方法能够提供封装结构内应力-应变场分布的准确和可靠的结果,为微电子封装的可靠性分析提供重要的依据. For the reliability analysis of electronic packages, strains in very localized areas, such as an interconnection or a corner, need to be determined. In this paper, a modified hybrid method of global/local modeling and real time moire interferometry is presented. In this method, a simplified, coarsely meshed global model is developed to get rough information about the deformation of the microelectronic package. In order to make sure the global model has been reasonably simplified and the material properties ...

An improvement and proof of OGY method

OGY method is the most important method of controlling chaos. It stabilizes a hyperbolic periodic orbit by making small perturbations for a system parameter. This paper improves the method of choosing parameter, and gives a mathematics proof of it.

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.

High-order asymptotic field of tensile plane-strain nonlinear crack problems

In this paper, the governing equations and the analytical method of the secondorder asymptotic field for the plane-straln crack problems of mode I have been presented. The numerical calculation has been carried out. The amplitude coefficients k2 of the second term of the asymptotic field have been determined after comparing the present results with some fine results of the finite element calculation. The variation of coefficients k2 with changes of specimen geometry and developments of plastic zone have been analysed. It is shown that the second term of the asymptotic field has significant influence on the near-crack-tip field. Therefore, we may reasonably argue that both the J-integral and the coefficient k2 can beeome two characterizing parameters for crack initiation.

Design of a Self-Adaptive Down-Hole Layered Steam Injector and Its Application

Layered steam injection, widely used in Liaohe Oilfield at Present, is an effective recovery technique to heavy oil reserves. Which makes the steam front-peak push forward uniformly, the amount of steam injection be assigned rationally, and the effect of injection steam be obtained as expected. To maintain a fixed ratio of layered steam injection and solve the problem of nonadjustable hole diameter with the change of layer pressure in the existing injectors, a new method is proposed in this paper to design layered steam injectors based on the dynamic balance theory. According to gas-liquid two-phase flow theory and beat transfer theory, the energy equation and the heat conduction equation in boreholes are developed. By analyzing the energy equilibrium of water-steam passing through the injector hole, we find an expression to describe the relation between the cross-sectional area of injector hole and the layer pressure. With this expression, we provide a new set of calculation methods and write the corresponding computer program to design and calculate the main parameters of a steam injector. The actual measurement data show that the theoretically calculated results are accurate, the software runs reliably, and they provide the design of self-adjustable layered steam injectors with the theoretical foundation.

Design of a Self-Adaptive Down-Hole Layered Steam Injector and Its Application

Layered steam injection, widely used in Liaohe Oilfield at present, is an effective recovery technique to heavy oil reserves. Which makes the steam front-peak push forward uniformly, the amount of steam injection be assigned rationally, and the effect of injection steam be obtained as expected. To maintain a fixed ratio of layered steam injection and solve the problem of nonadjustable hole diameter with the change of layer pressure in the existing injectors, a new method is proposed in this paper to design layered steam injectors based on the dynamic balance theory According to gas-liquid two-phase flow theory and heat transfer theory, the energy equation and the heat conduction equation in boreholes are developed. By analyzing the energy equilibrium of water-steam passing through the injector hole, we find an expression to describe the relation between the cross-sectional area of injector hole and the layer pressure. With this expression, we provide a new set of calculation methods and write the corresponding computer program to design and calculate the main parameters of a steam injector. The actual measurement data show that the theoretically calculated results are accurate, the software runs reliably, and they provide the design of self-adjustable layered steam injectors with the theoretical foundation.

Theoretical analysis of optical properties of dielectric coatings dependence on substrate subsurface defects

A model for refractive index of stratified dielectric substrate was put forward according to theories of inhomogeneous coatings. The substrate was divided into surface layer, subsurface layer and bulk layer along the normal direction of its surface. Both the surface layer (separated into N-1 sublayers of uniform thickness) and subsurface layer (separated into N-2 sublayers of uniform thickness), whose refractive indices have different statistical distributions, are equivalent to inhomogeneous coatings, respectively. And theoretical deduction was carried Out by employing characteristic matrix method of optical coatings. An example of mathematical calculation for optical properties of dielectric coatings had been presented. The computing results indicate that substrate subsurface defects can bring about additional bulk scattering and change propagation characteristic in thin film and Substrate. Therefore, reflectance, reflective phase shift and phase difference of an assembly of coatings and substrate deviate from ideal conditions. The model will provide some beneficial theory directions for improving optical properties of dielectric coatings via substrate surface modification. (c) 2005 Elsevier B.V. All rights reserved.

Theoretical analysis of influence of substrate microdefects on polarized light properties of dielectric coatings

A model for refractive index of stratified dielectric substrate was presented according to inhomogeneous coatings theories. The substrate was divided into surface layer, subsurface layer and bulk layer along the normal direction of its surface. The former two layers were equivalent to inhomogeneous coatings. Theoretical deduction was executed by employing the characteristic matrix method of optical coatings, and one mathematical calculation example was presented. The results indicate that reflectance, reflective phase shift and phase difference of polarized light deviate from ideal conditions. It shows that substrate microdefects can induce volume scattering and change propagation characteristic of light both in coatings and substrate. (c) 2005 Elsevier GmbH. All rights reserved.

Fabrication of photonic crystals on several kinds of semiconductor materials by using focused-ion beam method

In this paper, we introduced the fabrication of photonic crystals on several kinds of semiconductor materials by using focused-ion beam machine, it shows that the method of focused-ion beam can fabricate two-dimensional photonic crystal and photonic crystal device efficiently, and the quality of the fabricated photonic crystal is high. Using the focused-ion beam method, we fabricate photonic crystal wavelength division multiplexer, and its characteristics are analyzed. (c) 2007 Elsevier B.V. All rights reserved.

Determination of the values of hole-mixing coefficients due to interface and electric field in GaAs/AlxGa1-xAs superlattices

Numerical calculations within the envelope function framework have been performed to analyze the relations between the magnitude of in-plane optical anisotropy and the values of the additional hole-mixing coefficients due to interface and electric field in (001) symmetric GaAs/AlxGa1-xAs superlattices for light propagating along the [001] direction. It is found that the heavy- and light-hole states are mixed independently by interface and electric field. The numeric results demonstrate that the line shape of the in-plane anisotropic spectrum is determined by the ratio of the two hole-mixing coefficients. Theoretical analysis shows that with the help of simple calculation of the anisotropy at k=0, reliable values of the hole-mixing coefficients can be determined by reflectance-difference spectroscopy (IDS) technique, demanding no tedious fitting of experimental curves. The in-plane optical anisotropy measured by RDS provides a new method of getting the information on buried interfaces through the Value of the hole-mixing coefficient due to interface.

Fast Discrimination of Commerical Corn Varieties Based on Near Infrared Spectra

The existing methods for the discrimination of varieties of commodity corn seed are unable to process batch data and speed up identification, and very time consuming and costly. The present paper developed a new approach to the fast discrimination of varieties of commodity corn by means of near infrared spectral data. Firstly, the experiment obtained spectral data of 37 varieties of commodity corn seed with the Fourier transform near infrared spectrometer in the wavenurnber range from 4 000 to 12 000 cm (1). Secondly, the original data were pretreated using statistics method of normalization in order to eliminate noise and improve the efficiency of models. Thirdly, a new way based on sample standard deviation was used to select the characteristic spectral regions, and it can search very different wavenumbers among all wavenumbers and reduce the amount of data in part. Fourthly, principal component analysis (PCA) was used to compress spectral data into several variables, and the cumulate reliabilities of the first ten components were more than 99.98%. Finally, according to the first ten components, recognition models were established based on BPR. For every 25 samples in each variety, 15 samples were randomly selected as the training set. The remaining 10 samples of the same variety were used as the first testing set, and all the 900 samples of the other varieties were used as the second testing set. Calculation results showed that the average correctness recognition rate of the 37 varieties of corn seed was 94.3%. Testing results indicate that the discrimination method had higher precision than the discrimination of various kinds of commodity corn seed. In short, it is feasible to discriminate various varieties of commodity corn seed based on near infrared spectroscopy and BPR.