Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Overview of In - Situ Biodegradation and Enhancement

Microbial degradation technologies have been developed to restore ground water quality in aquifers polluted by organic contaminants effectively in recent years. However, in course of the degradation, the formation of biofilms in ground water remediation technology can be detrimental to the effectiveness of a ground water remediation project. Several alternatives are available to a remedial design engineer, such as Permeable Reactive Barriers (PRBs) and in -situ bioremediation, Hydrogen Releasing Compounds (HRCs) barrier, Oxygen Releasing Compounds (ORCs) barrier etc. which are efficient and cost- effective technologies. Excessive biomass formation renders a barrier ineffective in degrading the contaminants, Efforts are made to develop kinetics models which accurately determine bio - fouling and bio - filn formation and to control excessive biomass formation.     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

Discovery of 13 New Variable Stars in the Field of the Open Cluster NGC 2168 （M35）

A wide-field time-series CCD photometric survey of variable stars in the field of the open cluster NGC 2168 was carried out using the BATC Schmidt telescope. In total 13 new variable stars are discovered with three W UMa systems, one EA type and two EB type eclipsing binaries (one of them could be a W UMa system), and seven pulsating stars including three candidates of δScuti stars.     

Multi-wavelength Diagnostics of Bombardment by Non-thermal Particles in Solar Flares

Chromospheric lines, including Ha, Lyα, Lyβand CaⅡK, CaⅡ8542, are systemically and quantitatively investigated with respect to the non-thermal excitation and ionization due to particle beam bombardment for a series of solar semi-empirical atmospheric models. As a result we propose to use the contrast in the integrated intensity of hydrogen lines to estimate the total energy flux of the bombarding beam during the solar flare impulsive phase. Partial frequency redistribution is considered in the Lyαline calculation and a smaller intensity enhancement in the Ha line-centers is found than in the previous results of Fang et al.     

Influence of the Choice of Core-Envelope Transition Point on the Binary Merger of Two Main-sequence Components

We have studied the influence of different choices of core-envelope transition point on the final merger of contact binaries with two main-sequence components. A binary of 1.00 + 0.90M⊙ with an initial orbital period of 0.35d is examined. The mass fraction of the primary mixed with the matter of the secondary, qmix, determined by the chosen core-envelope transition point, ranges from 0.04 to 1.00 in our analysis. If as qmix< 0.8, none of the     

The Difference between the α-disks of Seyfert 1 Galaxiesand Quasars

In a previous paper, it was suggested that contamination of the nuclear luminosity by the host galaxy plays an important role in determining the parameters of the standard a disk of AGNs. Using the nuclear absolute B band magnitude instead of the total absolute B band magnitude, we have recalculated the central black hole masses, accretion rates and disk inclinations for 20 Seyfert 1 galaxies and 17 Palomar-Green (PG) quasars. It is found that a small value of a is needed for the Seyfert 1 galaxies than for the PG quasars. This difference in a possibly leads to the different properties of Seyfert 1 galaxies and quasars. Furthermore, we find most of the objects in this sample are not accreting at super-Eddington rates if we adopt the nuclear optical luminosity in our calculation.     

Cygnus A

Cygnus A was the first hyper-active galaxy discovered, and it remains by far the closest of the ultra-luminous radio galaxies. As such, Cygnus A has played a fundamental role in the study of virtually all aspects of extreme activity in galaxies. We present a review of jet theory for powering the double-lobed radio emitting structures in powerful radio galaxies, followed by a review of observations of Cygnus A in the radio, optical, and X-ray relevant to testing various aspects of jet theory. Issues addressed include: jet structure from pc- to kpc-scales, jet stability, confinement, composition, and velocity, the double shock structure for the jet terminus and the origin of multiple radio hotspots, the nature of the filamentary structure in the radio lobes, and the hydrodynamic evolution of the radio lobes within a dense cluster atmosphere, including an analysis of pressure balance between the various gaseous components. Also discussed are relativistic particle acceleration and loss mechanisms in Cygnus A, as well as magnetic field strengths and geometries both within the radio source, and in the intracluster medium. We subsequently review the classification, cluster membership, and the emission components of the Cygnus A galaxy. The origin of the activity is discussed. Concentrating on the nuclear regions of the galaxy, we review the evidence for an obscured QSO, also given the constraints on the orientation of the radio source axis with respect to the sky plane. We present an overview of models of central engines in AGN and observations of Cygnus A which may be relevant to testing such models. We conclude with a brief section concerning the question of whether Cygnus A is representative of powerful high redshift radio galaxies. Received October 10, 1995