Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

Iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact

We present in this paper an iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact. To retain the advantage of integration on an element-by-element basis which is at the heart of modern hydrocodes, we suggest that the first step should be to solve for accelerations at an advanced time step by using the lumped mass approach, then iterate using a consistent mass matrix to improve the estimate. Examples are given to show the improved resolution with the new method.

Surface waves in a vertically excited circular cylindrical container

The nonlinear free surface amplitude equation, which has been derived from the inviscid fluid by solving the potential equation of water waves with a singular perturbation theory in a vertically oscillating rigid circular cylinder, is investigated successively in the fourth-order Runge-Kutta approach with an equivalent time-step. Computational results include the evolution of the amplitude with time, the characteristics of phase plane determined by the real and imaginary parts of the amplitude, the single-mode selection rules of the surface waves in different forced frequencies, contours of free surface displacement and corresponding three-dimensional evolution of surface waves, etc. In addition, the comparison of the surface wave modes is made between theoretical calculations and experimental measurements, and the results are reasonable although there are some differences in the forced frequency.

A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

区域尺度降雨径流估算方法研究Ⅰ-算法设计

径流是区域水土流失的基础。基于已有试验研究和观测数据,结合GIS空间分析功能,初步提出了区域径流计算的基本思路,并对降雨径流产生的各个环节做出了算法设计,考虑的过程包括降水、植被截留、入渗、微地形存储、地表径流等。基于DEM将流域划分为规则网格并以此为基本计算单元,将月降水过程划分为若干时段作为计算迭代的基本单元,本算法可以计算出区域内每一网格单元任一时段末的地表径流量,为进一步建立区域水土流失模型奠定了基础。

地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

地下水位上升造成的黄土远程滑坡形成机理与数值模拟研究

Based on field survey, laboratory testing and numerical modeling, engineering characteristics of undisturbed loess and the mechanism of long-runout loess landslides caused by underground water level rise, as well as the formation conditions and spatial distribution of landslides, are systematically studied and analyzed. Loess landslides at south Plateau of Jingyang County are mainly classified as flowslide, slide and fall. Flowslide is the main type characteristic of high velocity, long runout and multi-stages. The steep relief composed of loose structured loess-old aged soil serials and the rise of groundwater table are the predominant conditions for landslides in the study area. To study loess mechanic poperties and loess landslides mechanisims, isotropically and anisotropically consolidated undrained compression(ICU and ACU) tests and constant-deviator-drained compression (CQD) tests were carried out on undisturbed samples. The results of undrained compression tests performed at the in-situ stress level show that the soils are of consistently strain-softening in the stress-strain relations and cause high excess pore pressure. The steady-state line and the potential region of instability are obtained from ICU and ACU test results. A necessary condition for liquefaction is that the soil state initially lies in or is brought into the potential instability region. In addition, a strong strain-softening model is also formed. CQD tests demonstrate that the mobilized friction angle is far less than the steady-state angle and that the soil experiences undrained contractive failure suddenly at very small strains when its stress path during drained loading tries to cross the potential instability region，thus validates the proposed instability region. Based on the location of the region of potential instability and the stress state of slope soil, a method of static liquefaction analysis is proposed for loess landslides caused by rise in groundwater table. Compared with other liquefaction analysis methods, this method overcomes the limitations inherent in conventional slope stability method and undrained brittleness index method. Triaxial tests composed of constant water content (CW) and wetting tests at constant deviator stress are performed on undisturbed unsaturated samples. The stress-strain relation of CW tests takes on strain-hardening behavior; The results of wetting tests at constant deviator stress designed to study the mechanics of failure of unsaturated loess caused by an increase in the degree of saturation (wetting) shows that a contractive failure occurs in the undisturbed samples. On the basis of the above triaxial test results, the initiation of static liquefaction is presented for long-runout loess landslides caused by rise in groundwater table, that is, the loess slope soil gradually transfer from unsaturated to saturated state under the infiltration of irrigation. A contractive failure occurs in the local region at very small strain by increasing the pore-water pressure at constant deviator stresses under drained conditons. It is the contractive failrue resulting from rise of pore pressure that leads to high excess pore pressure in the neighbour soil which reduces shear resistance of soil. The neighbour soils also fail due to the rapid increase in pore-water pressure. Thus a connected failure surface is developed quickly and a flowslide occurs. Based on the saturated-unsaturated seepage theory, transient seepage is computed using the finite element method on loess slope under groundwater table rise. Pore-water pressure distribution for every time step after irrigation are obtained. The phreatic surface in the slope increases with the groundwater table. Pore-water pressure distribution within 8m above the phreatic surface changes very quickly，but the water content and pore water pressure in the region ranging from 8m above the phreatic surface up to ground surface is almost not affected and the matric suction usually is kept at 100～120 kPa. Based on the results of laboratory tests and seepage flow analysis, the development process of loess landslide is modeled considering groundwater table rise. The shearing plastic zone first occurs at the slope toe where the soil is soaked for long term during rise in groundwater table. As irrigation continues, the shearing plastic zone gradually extends to the interior soils, with the results that the tensile plastic zone occurs at the slope crown. As time goes on, both the shearing plastic zone and tensile plastic zone continue to extend. Then a connected plastic zone is formed and fowslide occurs. In comparision to laboratory test results, the results of numerical simulation quite well verify the presented mechanism of static liquefaction of long-runout loess landslides caused by rise in groundwater table.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

Measurement of thermal rise-time of a laser diode based on spectrally resolved waveforms

Thermal resistance and thermal rise-time are two basic parameters that affect most of the performances of a laser diode greatly. By measuring waveforms received after a spectroscope at wavelengths varied step-by-step, the spectrally resolved waveforms can be converted to calculate the thermal rise-time. Basic formulas for the spectrum variation of a laser diode and the measurement set-up by using a Boxcar are described in the paper. As an example, the thermal rise-time of a p-side up packaged short-pulse laser diode was measured by the method to be 390 mu s. The method will be useful in characterizing diode lasers and LID modules in high-power applications. (c) 2005 Elsevier B.V. All rights reserved.

Enhancement of photoluminescence intensity of GaInNAs/GaAs quantum wells by two-step rapid thermal annealing

We investigate the effect of rapid thermal annealing on InGaNAs/GaAs quantum wells. At optimized annealing temperatures and times, the greatest enhancement of the photoluminescence intensity is obtained by a special two-step annealing process. To identify the mechanism affecting the material quality during the rapid thermal annealing, differential temperature analysis is applied, and temperature- and power-dependent photoluminescence is carried out on the samples annealed under different conditions. Our experiment reveals that some composition redistribution or other related ordering process may occur in the quantum-well layer during annealing. Annealing at a lower temperature for a long time primarily can remove defects and dislocations while annealing at a higher temperature for a short time primarily homogenizes the composition in the quantum wells.

The structure and current-voltage characteristics of multi-sheet InGaN quantum dots grown by a new multi-step method

Multi-sheet InGaN/GaN quantum dots (QDs) were grown successfully by surface passivation processing and low-temperature growth in metalorganic chemical vapor deposition. This method based on the principle of increasing the energy barrier of adatom hopping by surface passivation and low-temperature growth, is quite different from present methods. The InGaN quantum dots in the first layer of about 40-nm-wide and 15-nm-high grown by this method were revealed by atomic force microscopy. The InGaN QDs in upper layer grew bigger. To our knowledge, the current-voltage characteristics of multi-sheet InGaN/GaN QDs were measured for the fist time. Two kinds of resonance-tunneling-current features were observed which were attributed to the low-dimensional localization effect. Some current peaks only appeared in positive voltage for sample due to the non-uniformity of the QDs in the structure. (C) 2002 Elsevier Science B.V. All rights reserved.

An automatic step adjustment method for average power analysis technique used in fiber amplifiers

An automatic step adjustment (ASA) method for average power analysis (APA) technique used in fiber amplifiers is proposed in this paper for the first time. In comparison with the traditional APA technique, the proposed method has suggested two unique merits such as a higher order accuracy and an ASA mechanism, so that it can significantly shorten the computing time and improve the solution accuracy. A test example demonstrates that, by comparing to the APA technique, the proposed method increases the computing speed by more than a hundredfold under the same errors. By computing the model equations of erbium-doped fiber amplifiers, the numerical results show that our method can improve the solution accuracy by over two orders of magnitude at the same amplifying section number. The proposed method has the capacity to rapidly and effectively compute the model equations of fiber Raman amplifiers and semiconductor lasers. (c) 2006 Optical Society of America

Ultrafast Kerr rotations and zero-field dephasing time of electron spins in InAs/GaAs quantum disks

Time-resolved Kerr rotation (TRKR) measurements based on pump-probe arrangement were carried out at 5 K on the monolayer fluctuation induced InAs/GaAs quantum disks grown on GaAs substrate without external magnetic field. The lineshape of TRKR signals shows an unusual dependence on the excitation wavelength, especially antisymmetric step-shaped structures appearing when the excitation wavelength was resonantly scanned over the heavy- and light-hole subbands. Moreover, these step structures possess an almost identical decay time of similar to 40 Ps which is believed to be the characteristic spin dephasing time of electrons in the extremely narrow InAs/GaAs quantum disks.

Application of step-scan FTIR to the research of quantum cascade lasers

The principle of step-scan Fourier transform infrared （FTIR） spectroscopy is introduced. Double modulation step-scan FTIR technique is used to obtain the quantum cascade laser＇s stacked emission spectra in the time domain. Optical property and thermal accumulation of devices due to large drive current are analyzed.

One-step aqueous solvothermal synthesis of In2O3 nanocrystals

Highly crystalline and nearly monodisperse In2O3 nanocrystals with both cube and flower shapes were successfully synthesized in one step through a facile aqueous solvothermal method for the first time, free of any surfactant or template. X-ray diffraction (XRD), transmission electron microscopy (TEM), selective area electron diffraction (SAED), and high-resolution transmission electron microscopy (HRTEM) were used to characterize the samples. In our work, the use of diethylene glycol (DEG) is a crucial factor for the formation of the In2O3 phase.