海底沉积层中水合物成藏过程模拟

水合物成藏是一种动态演变过程，从生成和分解速度相对大小角度提出了天然气水合物的动态成藏理论；分析指出天然气水合物储层状态主要存在三种类型：成长型（包括渗透型与扩散型）、成熟型和消退型，必须结合经济型指标对这三种类型水合物藏进行开采可行性研究。 通过在自制的透明鼓泡反应釜内模拟气体迁移与水合物形成实验，解释了气体迁移对体系内传热与传质、溶解性质、过压等的影响因，进一步指出需要研究新的水合物相图。 从传热角度建立了水合物生长模型，计算并比较了南海北部陆坡琼东南盆地里甲烷水合物在渗漏系统和扩散系统下的生长速度，指出前者比后者快约20-40倍。 建立了水-气（游离气）-水合物-盐反应的多相流模型，演绎了在四个不同时刻随着水合物的形成，溶解度、渗透率、毛细压力、饱和度及盐度等的联动演变关系，并分析和比较了在沉积层水合物稳定区域内不同游离气饱和度和甲烷渗漏通量情况下水合物的形成过程和三种成藏类别。 反演了神狐海域水合物的形成，指出该区很可能是属于发育阶段的成长型渗漏系统水合物储藏，按照甲烷渗漏通量为0.5kg/m2·a计算，该地区已演化了约4800a，还需要约2700ka孔隙中水合物可到最大体积饱和度，达约75%，并在此后游离气将溢出沉积层进入海水中；计算得最终的甲烷水合物储量约为1.47×109m3，即该区最终将储存着约2200亿立方米甲烷气，是现在预测的约13倍！

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

Structure of Vortex Dislocations in a Wake-type Flow

Structure and dynamical processes of vortex dislocations in a kind of wake-type flow are described clearly by vortex lines, which are directly constructed from data of three-dimensional direct numerical simulations of the flow evolution.

Investigation of Mn-Implanted n-Type Ge

Mn+ irons were implanted to n-type Ge(1 1 1) single crystal at room temperature with an energy of 100 keV and a dose of 3 x 10(16) cm(-2). Subsequently annealing was performed at 400degreesC for 1 h under flowing nitrogen gas. X-ray diffraction measurements show that as-implanted sample is amorphous and the structure of crystal is restored after annealing. Polycrystalline germanium is formed in annealed sample. There are no new phases found except germanium. The samples surface morphologies indicate that annealed sample has island-like feature while there is no such kind of characteristic in as-implanted sample. The elemental composition of annealed sample was analyzed by Auger electron spectroscopy. It shows that manganese ions are deeply implanted into germanium substrate and the highest manganese atomic concentration is 8% at the depth of 120 nm. The magnetic properties of samples were investigated by an alternating gradient magnetometer. The annealed sample shows ferromagnetic behavior at room temperature.

Low-Dimensional Dynamical Systems for a Wake-Type Shear Flow with Vortex Dislocations

Low-dimensional systems are constructed to investigate dynamics of vortex dislocations in a wake-type shear flow. High-resolution direct numerical simulations are employed to obtain flow snapshots from which the most energetic modes are extracted using proper orthogonal decomposition (POD). The first 10 modes are classified into two groups. One represents the general characteristics of two-dimensional wake-type shear flow, and the other is related to the three-dimensional properties or non-uniform characteristics along the span. Vortex dislocations are generated by these two kinds of coherent structures. The results from the first 20 three-dimensional POD modes show that the low- dimensional systems have captured the basic properties of the wake-type shear flow with vortex dislocation, such as two incommensurable frequencies and their beat frequency.

Magnetic Properties of Mn-Implanted n-Type Ge

Mn+ ions were implanted into n-type Ge(111) single crystal at room temperature at an energy of 100 keV with a dose of 3 x 1016 cm-2. Subsequent annealing was performed on the samples at 400 °C and 600 °C in a flowing nitrogen atmosphere. The magnetic properties of the samples have been investigated by alternating gradient magnetometer at room temperature. The compositional properties of the annealed samples were studied by Auger electron spectroscopy and the structural properties were analyzed by X-ray diffraction measurements. Magnetization measurements reveal room-temperature ferromagnetism for the annealed samples. The magnetic analysis supported by compositional and structural properties indicates that forming the diluted magnetic semiconductor (DMS) MnxGe1-x after annealing may account for the ferromagnetic behavior in the annealed samples.

Quantifying Cell Binding Kinetics Mediated By Surface-Bound Blood Type B Antigen To Immobilized Antibodies

Cell adhesion is crucial to many biological processes, such as inflammatory responses, tumor metastasis and thrombosis formation. Recently a commercial surface plasmon resonance (SPR)-based BIAcore biosensor has been extended to determine cell binding mediated by surface-bound biomolecular interactions. How such cell binding is quantitatively governed by kinetic rates and regulating factors, however, has been poorly understood. Here we developed a novel assay to determine the binding kinetics of surface-bound biomolecular interactions using a commercial BIAcore 3000 biosensor. Human red blood cells (RBCs) presenting blood group B antigen and CM5 chip bearing immobilized anti-B monoclonal antibody (mAb) were used to obtain the time courses of response unit, or sensorgrams, when flowing RBCs over the chip surface. A cellular kinetic model was proposed to correlate the sensorgrams with kinetic rates. Impacts of regulating factors, such as cell concentration, flow duration and rate, antibody-presenting level, as well as pH value and osmotic pressure of suspending medium were tested systematically, which imparted the confidence that the approach can be applied to kinetic measurements of cell adhesion mediated by surface-bound biomolecular interactions. These results provided a new insight into quantifying cell binding using a commercial SPR-based BIAcore biosensor.

Lattice type of fracture model for concrete

Concrete is usually described as a three-phase material, where matrix, aggregate and interface zones are distinguished. The beam lattice model has been applied widely by many investigators to simulate fracture processes in concrete. Due to the extremely large computational effort, however, the beam lattice model faces practical difficulties. In our investigation, a new lattice called generalized beam (GB) lattice is developed to reduce computational effort. Numerical experiments conducted on a panel subjected to uniaxial tension show that the GB lattice model can reproduce the load-displacement curves and crack patterns in agreement to what are observed in tests. Moreover, the effects of the particle overlay on the fracture process are discussed in detail. (C) 2007 Elsevier Ltd. All rights reserved.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

A new type of boundary integral-equation for plane problems of elasticity including rotational forces

Based on the idea proposed by Hu [Scientia Sinica Series A XXX, 385-390 (1987)], a new type of boundary integral equation for plane problems of elasticity including rotational forces is derived and its boundary element formulation is presented. Numerical results for a rotating hollow disk are given to demonstrate the accuracy of the new type of boundary integral equation.

Vortex dislocation in wake-type flow caused by local spanwise nonuniformity

DNS of spatiotemporal evolution of a wake-type flow is performed. In the incoming flow, a local spanwise nonuniformity in momentum defect is initially imposed. Results show that the spanwise nonuniformity leads to a series of symmetric twist vortex dislocation in downstream of the flow. Vortex line variations and substantial transition of vorticity from spanwise to the streamwise and vertical directions clearly feature the generation of a vortex dislocation and the real vortex linking in the dislocation. Dynamical process and the mechanism responsible for the vortex dislocation are described.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.