三维精细积分法地震正演及边界条件

 将任意差分精细积分法用于三维波动方程地震正演，关键在于如何消除数值计算中有限波场区域边界引起的边界反射。文中采用Berenger给出的电磁波完全匹配层吸收边界条件，推导出三维波动方程任意差分精细积分法地震正演的完全匹配层吸收边界条件计算公式，并给出了完全匹配层吸收边界条件算例。计算结果表明，此方法压制边界反射效果明显。三维波动方程地震正演模拟实例表明，完全匹配层吸收边界条件的任意差分精细积分法为复杂区地震波传播规律研究提供了一种实用的正演模拟工具     

Effect of pH on the Crystallization of Homogeneously Precipitated Nanosized PZT Powder

Nanosized lead zirconate titanate (PZT) powder with Zr:Ti ratio in the morphotropic phase boundary region was synthesized by homogeneous precipitation of metal ions. The powder precipitated at 90°C and at pH 6.7 resulted single-phase perovskite lead zirconate titanate powder when calcined at 550°C and above for 4 hours in air. The solution pH and the precipitation temperature strongly affect the composition of the calcined powder. The results obtained by structural characterization of homogeneously precipitated powder were compared with that obtained from the conventional precipitation method using ammonia in terms of crystallization, homogeneity, and microstructure. The homogeneously precipitated powder showed smaller particle size, minimum agglomeration and uniform shape on calcination and annealing. Powdered samples that precipitated by homogeneous precipitation crystallized directly to perovskite PZT, without any intermediate pyrochlore phase formation. In contrast, the NH₃ precipitated powder converted to perovskite PZT via metastable pyrochlore and it showed phase segregation upon annealing at higher temperatures. The reaction kinetics has been studied by X-ray diffraction, differential thermal analysis, and differential scanning calorimetry.     

Flow Control Through the Use of Topography

In this work, optimal shaft shapes for flow in the annular space between a rotating shaft with axially-periodic radius and a fixed coaxial outer circular cylinder, are investigated. Axisymmetric steady flows in this geometry are determined by solving the full Navier-Stokes equations in the actual domain. A measure of the flow field, a weighted convex combination of the volume averaged square of the L₂-norm of the velocity and vorticity vectors, is employed. It has been demonstrated that boundary shape can be used to influence the characteristics of the flow field, such as its velocity component distribution, kinetic energy, or even vorticity. This ability to influence flow fields through boundary shape may be employed to improve microfluidic mixing or, possibly, to minimize shear in biological applications.     

Symmetric BEM formulations for 2-D steady-state and transit potential problems

This paper presents a new concept for symmetric boundary element method (SBEM) applicable to 2-D steady-state and transit potential problems. Two kinds of SBEM formulations are derived. Symmetry is obtained simply through matrix manipulation, and no hypersingularity appears. Therefore, SBEM is much easier than the traditional symmetric Galerkin BEM. Compared with the traditional asymmetric BEM, the present SBEM can reduce the computational cost for time domain problems only. However, when applied to BEM/FEM coupling procedure, SBEM can reduce the computational cost for both steady-state and time domain problems. Three numerical examples are included to illustrate the effectiveness and accuracy of the present formulations.     

Transition to turbulence in rotating-disk boundary layers—convective and absolute instabilities

This work is an experimental study of mechanisms for transition to turbulence in the boundary layer on a rotating disk. In one case, the focus was on a triad resonance between pairs of traveling cross-flow modes and a stationary cross-flow mode. The other was on the temporal growth of traveling modes through a linear absolute instability mechanism first discovered by Lingwood (1995, J Fluid Mech 314:373–405). Both research directions made use of methods for introducing controlled initial disturbances. One used a distributed array of ink dots placed on the disk surface to enhance a narrow band of azimuthal and radial wave numbers of both stationary and traveling modes. The size of the dots was small so that the disturbances they produce were linear. Another approach introduced temporal disturbances by a short-duration air pulse from a hypodermic tube located above the disk and outside the boundary layer. Hot-wire sensors primarily sensitive to the azimuthal velocity component, were positioned at different spatial (r,θ) locations on the disk to document the growth of disturbances. Spatial correlation measurements were used with two simultaneous sensors to obtain wavenumber vectors. Cross-bicoherence was used to identify three-frequency phase locking. Ensemble averages conditioned on the air pulses revealed wave packets that evolved in time and space. The space–time evolution of the leading and trailing edges of the wave packets were followed past the critical radius for the absolute instability, r _{c A} . With documented linear amplitudes, the spreading of the disturbance wave packets did not continue to grow in time as r _{c A} was approached. Rather, the spreading of the trailing edge of the wave packet decelerated and asymptotically approached a constant. This result supports the linear DNS simulations of Davies and Carpenter (2003, J Fluid Mech 486:287–329) who concluded that the absolute instability mechanism does not result in a global mode, and that linear-disturbance wave packets are dominated by the convective instability. In contrast, wave-number matching between traveling cross-flow modes confirmed a triad resonance that lead to the growth of a low azimuthal number (n = 4) stationary mode. At transition, this mode had the largest amplitude. Signs of this mechanism can be found in past flow visualization of transition to turbulence in rotating disk flows.     

Non-linear current–voltage characteristics related to native defects in SrTiO3 and ZnO bicrystals

SrTiO₃ and ZnO bicrystals with various types of boundaries were fabricated in order to examine their current–voltage characteristics across single grain boundaries. Their grain boundary structures were also investigated by high-resolution transmission electron microscopy. In Nb-doped SrTiO₃, electron transport behaviors depend on the type of boundaries. Random type boundaries exhibit highly non-linear current–voltage characteristics, while low angle boundaries show a slight non-linearity. On the contrary, undoped ZnO does not exhibit non-linear current–voltage characteristics in any type of boundaries including random ones. It is suggested that the differences observed in current–voltage properties between the two systems are mainly due to the difference in the accumulation behavior of acceptor-like native defects at grain boundaries. A clear non-linearity is obtained by means of Co-doping even for the highly coherent Σ1 boundary in a ZnO bicrystal. This is considered to result from the production of acceptor-like native defects by Co-doping.     

Exact integration in the boundary element method for two-dimensional elastostatic problems

This paper derives the exact integrations for the integrals in the boundary element analysis of two-dimensional elastostatics. For facilitation, the derivation is based on the simple forms of the fundamental functions by taking constant, discontinuous linear and discontinuous quadratic elements as examples. The efficiency and accuracy of the derived exact integrations are verified against five benchmark problems; the results indicate that the derived exact integrations significantly reduces the CPU time for forming the matrices of the boundary element analysis and solving the internal displacements.     

有限元数值试井解释技术在复杂断块油气藏评价中的应用

有限元数值模拟试井技术,可通过描述不同网格的瞬变压力对测试范围内各个小单元进行精细描述.综合运用解析试井、数值试井联带分析方法,将测试动态成果与地质静态成果有机结合、相互补充,可实现对复杂断块油气藏的精细描述,达到一口井搞清一个小型复杂断块油气藏的目的.     

An inverse solution methodology for estimating the diffusion coefficient of gases in Athabasca bitumen from pressure-decay data

An inverse solution methodology is developed for the estimation of diffusion coefficient of gases in highly viscous, oil-sands bitumens from isothermal, pressure-decay measurements. The approach involves modeling the rate of change in pressure using the diffusion equation for the liquid phase coupled with a mass balance equation for the gas phase. The inverse solution framework is utilized to arrive at two graphical techniques for estimating the diffusion coefficient. Both techniques involve the determination of the slope of a straight line resulting from plotting the experimental data in accordance with the developed model. An advantage of the proposed techniques is that the diffusion coefficient is estimated directly, i.e. without making it an adjustable parameter. The novelty of the proposed method is in its simplicity as well as its ability to isolate portions of the pressure-decay data affected by experimental fluctuations. The effect of the initial pressure on the predicted diffusion coefficient and pressure-decay profile was also investigated. The diffusion coefficients of CO₂, CH₄, C₂H₆ and N₂ in Athabasca bitumen at 50–90 °C and about 8 MPa were estimated and compared with literature values.     

地下水运动数值模拟过程中边界条件问题探讨

本文对地下水运动数值模拟过程中边界条件的涵义 和处理方法进行了分析和讨论。阐述了边界条件所包含的双重意义。指出随着人类活动影响 强度的日益增大，边界条件的处理要面临一些新的更为复杂的问题。在模型预报之前必须首 先对边界条件做出预报。边界条件的预报既要考虑自然因素的作用，同时也要考虑人类活动 (人工开采和人工补给)的影响及由于邻区水流条件变化而产生的耦合效应。之后，给出了两 个应用实例。