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.     

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

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

Radio Sources and Scintillation

A review is given of the interplay between studies of compact radio sources and the scattering and scintillations that occur as the signals travel through the irregular refractive index of the interstellar and interplanetary plasmas. This revised version was published online in July 2006 with corrections to the Cover Date.     

CCD Polarimetric Imaging Of Comet Hale-Bopp (C/1995 O1)

Polarimetric observations of the light scattered by dust have been carried out at Pic-du-Midi Observatory with the 2 m telescope in June and September–October 1996, and at Haute-Provence Observatory with the 0.80 m telescope in April 1997. They cover a total number of 11 nights and a large (6.9°–47.7°) phase angle range. The spatial resolution allows to underline structures in the coma, as well in the brightness images as in the polarization maps, with a correlation between the regions of bright structures and the regions of higher polarization. A clear difference appears between the sunward and antisunward side, with higher polarization on the antisunward side. The phase angle coverage allows us to obtain a polarimetric phase curve for the whole coma and to compare it with other cometary phase curves. The degree of polarization is higher for Hale-Bopp than for the comets previously observed. This revised version was published online in July 2006 with corrections to the Cover Date.     

The first World Atlas of the artificial night sky brightness

We present the first World Atlas of the zenith artificial night sky brightness at sea level. Based on radiance-calibrated high-resolution DMSP satellite data and on accurate modelling of light propagation in the atmosphere, it provides a nearly global picture of how mankind is proceeding to envelop itself in a luminous fog. Comparing the Atlas with the United States Department of Energy (DOE) population density data base, we determined the fraction of population who are living under a sky of given brightness. About two-thirds of the World population and 99 per cent of the population in the United States (excluding Alaska and Hawaii) and European Union live in areas where the night sky is above the threshold set for polluted status. Assuming average eye functionality, about one-fifth of the World population, more than two-thirds of the United States population and more than one half of the European Union population have already lost naked eye visibility of the Milky Way. Finally, about one-tenth of the World population, more than 40 per cent of the United States population and one sixth of the European Union population no longer view the heavens with the eye adapted to night vision, because of the sky brightness.     

Scattered light models of the dust shell around HD 179821

The dust shell around the evolved star HD 179821 has been detected in scattered light in near-IR imaging polarimetry observations. Here, we subtract the contribution of the unpolarized stellar light to obtain an intrinsic linear polarization of between 30 and 40 per cent in the shell that seems to increase with radial offset from the star. The J - and K -band data are modelled using a scattering code to determine the shell parameters and dust properties. We find that the observations are well described by a spherically-symmetric distribution of dust with an r ⁻² density law, indicating that when mass-loss was occurring, the mass-loss rate was constant. The models predict that the detached nature of a spherically-symmetric, optically-thin dust shell, with a distinct inner boundary, will only be apparent in polarized flux. This is in accordance with the observations of this and other optically-thin circumstellar shells, such as IRAS 17436+5003. By fitting the shell brightness we derive an optical depth to the star that is consistent with V -band observations and that, assuming a distance of 6 kpc, gives an inner-shell radius of     , a dust number density of     at r _in and a dust mass of     . We have explored axisymmetric shell models but conclude that any deviations from spherical symmetry in the shell must be slight, with an equator-to-pole density contrast of less than 2:1. We have not been able to fit simultaneously the high linear polarizations and the small     colour excess of the shell and we attribute this to the unusual scattering properties of the dust. We suggest that the dust grains around HD 179821 either are highly elongated or consist of aggregates of smaller particles.     

关于非球形粒子光散射的T-矩阵数值计算方法的研究

首先介绍了无规取向轴对称非球形粒子光散射的T－矩阵数值方法。然后，进一步研究了T－矩阵计算的收敛问题，提出了一种新的关于无规取向无吸收非线形散射物体的物理收敛方法，我们研究了椭球粒子的收敛问题，结果表明我们的数值方法和T－矩阵计算程序是有效的， 收敛精度与粒子的尺度和形状有很强的依赖性，在一定的条件下，我们的物理收敛速度优于NASA Mishchenko的数学收敛速序。     

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.