Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

A study on quality of aquatic environment in Tumen River Area

After the survey of pollution sources, a study on surface water quality assessment and forecast is given by means of grey system method, fuzzy mathematical method and multiple-index method. Based on it, aquatic environment quality features, treatment measures and environmental strategies of the area are proposed. The quality of aquatic environment of 5 rivers in the Tumen River area is studied. The results show that the pollution of surface water is serious; water quality of most rivers is between grade IV and V except the Hunchun River, being higher than grade IV standard; pollution levels of most rivers have been basically controlled except the of Burhatong River, which is deteriorating gradually. Pollutants of the rivers are comparatively regular, mainly are SS, COD, BOD, AR-OH, NH₃-N. The main pollution trades are chemical fibre industry, pulp and paper making industry and mining industry. If the growth rate of gross industrial product is higher than 25 percent under the encouraging-model of regional exploitation, the pollutants’ load will overtake the bearing capacity of aquatic environment. Thus some protection program against pollution must be worked out in order to achieve the harmonious development of economy society and environment. A project from “Studies on Earlier Stage of Regional Development” of the Chinese Academy of Sciences     

基于RDBMS的积件技术的应用研究

积件技术是教学软件发展的方向，探讨了用Power Builder前台结合SQL Server关系数据库后台进行多媒体教学积件管理的应用模式，并结合具体教学示例，进行了分析说明。     

Search for groups of extragalactic radio sources

We explore the possibility of searching for groups of radio sources from the FIRST catalog on angular scales 1′–5′. We developed an efficient method of searching for such groups that takes into account the need for combining the components of extended sources represented in the catalog by separate objects. We found 31 groups of radio sources with angular sizes <5′ that contain no fewer than five sources with flux densities ≥3 mJy. This number is at least triple the expected number of such groups for a random Poisson distribution of radio sources in the sky. The prospects for using groups of radio sources to detect and study distant systems of galaxies are discussed.     

天文(重力)水准的新手段及其计算方法

介绍当今确定天文大地垂偏差的新仪器－CCD自动天体测量仪和确定两点间高程异常差的新方法，该方法用单极坐标代替繁琐的双极坐标进行计算，讨论了天文重力水准的误差及垂线偏差非线性影响等问题，估计在不久的将来，用这一手段施测山区似大地水准面的精度可望达到厘米级。     

RXTE observations of the transient sources IGR J17091-3624 and IGR J18539+0727

We present the results of our analysis of the RXTE observations for two transient sources, IGR J17091-3624 and IGR J18539+0727, in April 2003. The derived energy spectra of the sources and the power-density spectra of their light curves make it possible to classify them as low/hard-state X-ray binaries. The parameters of the power spectrum for IGR J18539+0727 lead us to tentatively conclude that the compact object in this binary is a black hole.     

Hard X-ray spectrum of Mkn 421 during the active phase

Spectral measurement of Mkn 421 were made in the hard X-ray energy band of 20–200 keV using a high sensitivity, large area scintillation counter telescope on November 21, 2000 and these coincided with the onset of an active X-ray phase as seen in the ASM counting rates on board RXTE. The observed spectrum can not be fitted to a single power law similar to the PDS data of BeppoSAX. The data can be fitted both by a two component power-law function or a combination of an exponential function with a power law component at the high energies above 80 keV. We identify these components with those arising from the synchrotron self compton and the high energy power-law tail arising from the upgrading of the thermal photons due to multiple Compton scattering a la Cyg X-1. A comparison with the earlier data clearly suggests a spectral variability in the hard X-ray spectrum of the source. We propose a continuously flaring geometry for the source as the underlying mechanism for energy release.     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

地壳演化与成矿作用——以川滇地洼系“四层楼”铜矿床序列为例

提出川滇地洼系“四层楼”铜矿床序列的形成与陆壳演化的成生联系,是与本区陆壳由前地槽—地槽—地台—地洼演化各阶段与之相匹配的成矿作用的产物.与此同时,并总结了本区“四层楼”铜矿床序列的成矿作用具有明显的继承性、新生性、旋回性及层控性四大特点和多因复成矿床的成矿模式.