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.     

Wavelet Analysis of Space Solar Telescope Images

The scientific satellite SST (Space Solar Telescope) is an important research project strongly supported by the Chinese Academy of Sciences. Every day, SST acquires 50 GB of data (after processing) but only 10GB can be transmitted to the ground because of limited time of satellite passage and limited channel volume. Therefore, the data must be compressed before transmission. Wavelets analysis is a new technique developed over the last 10 years, with great potential of application. We start with a brief introduction to the essential principles of wavelet analysis, and then describe the main idea of embedded zerotree wavelet coding, used for compressing the SST images. The results show that this coding is adequate for the job.     

三峡重庆库区深部地球物理特征与断裂构造

为了深入研究三峡重庆库区岩石圈动力学特征及其对断裂构造活动的控制和影响．并为之提供基础资料，为三峡重庆库区地震、地质灾害的监测与防治提供基础依据，在已有地球物理资料的基础上，从综合地球物理研究角度出发，通过实测地震测线资料的再解释．采用新的处理技术方法，对本区的东西向主剖面和南北向支测线的地震测深资料进行二维射线追踪处理、Pg波成像；选用场分离技术、位移数字成像技术重新处理了重力和航磁资料，通过联合反演来建立深部二维构造剖面，对剖面所揭示的基底构造特征和地壳结构特征、主要断裂构造特征以及莫霍面的起伏特征进行了精细分析和细致研究．从地球物理平面场特征出发建立了岩石圈构造三维框架。研究结果表明，这样的研究思路准确、方法得当．结论可靠；沿该剖面，把可解译的断裂分为Ⅲ级：Ⅰ级为超岩石圈断裂；Ⅱ级为壳断裂；Ⅲ级为盖层断裂。依据地球物理特征，准确揭示了库区的断裂构造特征，达到了预期效果。     

Real-time difference imaging analysis of MOA Galactic bulge observations during 2000

We describe observations carried out by the MOA group of the Galactic bulge during 2000 that were designed to detect efficiently gravitational microlensing of faint stars in which the magnification is high and/or of short duration. These events are particularly useful for studies of extrasolar planets and faint stars. Approximately 17 deg² were monitored at a sampling rate of up to six times per night. The images were analysed in real time using a difference imaging technique. 20 microlensing candidates were detected, of which eight were alerted to the microlensing community whilst in progress. Approximately half of the candidates had high magnifications (≳10), at least one had very high magnification (≳50), and one exhibited a clear parallax effect. The details of these events are reported here, together with details of the on-line difference imaging technique. Some nova-like events were also observed and these are described, together with one asteroid.     

NORTHERN AND SOUTHERN HEMISPHERIC POLAR LOWS——A COMPARATIVE STUDY

Northern Hemispheric“polar lows”,i.e.maritime,small,but fairly intense cyclonic systems have been studied formore than two decades.Recently researchers have turned their attention towards the occurrence of polar lows and othersmall-scale cyclonic disturbances on the Southern Hemisphere.In this paper a short account of our present knowledgeof Northern Hemispheric polar lows will be given,followed by some preliminary results of an investigation of the natureand occurrence of polar lows on the Southern Hemisphere.     

Particle-mesh simulations of the Lyα forest

We have analysed short-exposure high-resolution images obtained on the William Herschel Telescope (WHT), Big Telescope Azimuthal (BTA) and CanadaFranceHawaii Telescope (CFHT) in order to investigate mechanical vibrations that are capable of reducing the high angular resolution allowed by speckle-interferometric methods. After filtering the photon noise we have computed power spectra of the image centres of gravity. In these spectra we have found vibration features that differ from one telescope to another and that vary in particular with the zenith angle. In the case of the WHT we discuss these results in terms of possible causes and in terms of degradation of the transfer function. We present some means to improve the data.     

Sequences in the Hardness Ratio-Peak Energy Plane of Gamma-Ray Bursts

The narrowness of the distribution of the peak energy of the vFv spectrum of gamma-ray bursts (GRBs) and the unification of GRB populations are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB populations, the long GRBs, the short GRBs, and the X-ray flashes (XRFs), in the HR - Ep plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that the long GRBs and the XRFs observed by HETE-2 seem to follow the same sequence in the HR-Ep plane, with the XRFs at the low end of this sequence. We fit the sequence by a universal Band function, and find that this sequence is mainly defined by the low energy index α, and is insensitive to the high energy index,β. With fixed β=-5, a best fit is given by α=-1.00 with X^2min/dof=2.2. The long and short GRBs observed by BATSE follow significantly different sequences in the HR-Ep plane, with most of the short GRBs having a larger hardness ratio than the long GRBs at a given Ep. For the long GRBs a best-fit yields α=-0.30 and β=-2.05. For the short GRBs, a best fit gives α=-0.60 with x^2min=1.1 (withβfixed at-2.0 because it is numerically unstable). The α value for the short GRBs is significantly greater than that for the long GRBs. These results indicate that the global spectral behaviors of the long GRB sample and the XRF sample are similar, while that of the short GRBs is different. The short GRBs seem to be a unique subclass of GRBs, and they are not the higher energy extension of the long GRBs.