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波成像；选用场分离技术、位移数字成像技术重新处理了重力和航磁资料，通过联合反演来建立深部二维构造剖面，对剖面所揭示的基底构造特征和地壳结构特征、主要断裂构造特征以及莫霍面的起伏特征进行了精细分析和细致研究．从地球物理平面场特征出发建立了岩石圈构造三维框架。研究结果表明，这样的研究思路准确、方法得当．结论可靠；沿该剖面，把可解译的断裂分为Ⅲ级：Ⅰ级为超岩石圈断裂；Ⅱ级为壳断裂；Ⅲ级为盖层断裂。依据地球物理特征，准确揭示了库区的断裂构造特征，达到了预期效果。     

Synthetic SH seismograms in a laterally varying medium by the discrete wavenumber method

Summary. We present a new method to calculate the SH wavefield produced by a seismic source in a half-space with an irregular buried interface. The diffracting interface is represented by a distribution of body forces. The Green's functions needed to solve the boundary conditions are evaluated using the discrete wavenumber method. Our approach relies on the introduction of a periodicity in the source-medium configuration and on the discretization of the interface at regular spacing. The technique developed is applicable to boundaries of arbitrary shapes and is valid at all frequencies. Some examples of calculation in simple configurations are presented showing the capabilities of the method.     

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.     

A Critical review of theoretical models of negatively polarized light scattered by atmosphereless solar system bodies

About a dozen physical mechanisms and models aspire to explain the negative polarization of light scattered by atmosphereless celestial bodies. This is too large a number for the reliable interpretation of observational data. Through a comparative analysis of the models, our main goal is to answer the question: Does any one model have an advantage over the others? Our analysis is based on new laboratory polarimetric and photometric data as well as on theoretical results. We show that the widely used models due to Hopfield and Wolff cannot realistically explain the phase-angle dependence of the degree of polarization observed at small phase angles. The so-called interference or coherent backscattering mechanism is the most promising model. Models based on that mechanism use well-defined physical parameters to explain both negative polarization and the opposition effect. They are supported by laboratory experiments, particularly those showing enhancement of negative polarization with decreasing particle size down to the wavelength of light. According to the interference mechanism, pronounced negative branches of polarization, like those of C-class asteroids, may indicate a high degree of optical inhomogeneity of light-scattering surfaces at small scales. The mechanism also seems appropriate for treating the negative polarization and opposition effects of cometary dust comae, planetary rings, and the zodiacal light.     

CygX—1硬态高能光子的时延

硬X射线和软射线光子的时延时研究主能辐射的一个重要方面，傅里叶交叉谱常常被用于教育处硬软光子之间的时延，但是交叉说示能在高于30Hz的傅里叶频率上从CygX-1测量到有统计意义的快速光变，由李惕培在时域上发展的交叉相关函方法能有效测量不同时间尺度上的时延，利用此交叉相关分析方法，讨论了不同观测时期CygX-1硬态高能光子时延 的性质，分析结果表明，CygX-1硬态在短的尺度（＜0．1秒）上存存在有意义的时间延迟，并将短时间尺度上的观测结果和各种CygX-1硬态模型进行了比较讨论。     

Laue diffraction lenses for astrophysics: From theory to experiments

Based on the laws of X-ray diffraction in crystals, Laue lenses offer a promising way to achieve the sensitivity and angular resolution leap required for the next generation of hard X-ray and gamma-ray telescopes.The present paper describes the instrumental responses of Laue diffraction lenses designed for nuclear astrophysics. Different possible geometries are discussed, as well as the corresponding spectral and imaging capabilities. These theoretical predictions are then compared with Monte-Carlo simulations and experimental results (ground and stratospheric observations from the CLAIRE project). PACS 95.55.Ka, 61.50.Ah, 61.10.−i, 41.50.+h     

Techniques to measure grain-size distributions of loamy sediments: a comparative study of ten instruments for wet analysis

Ten instruments for measuring the grain‐size distribution of loamy sediments in water were tested in a comparative study. Instruments tested were: the Malvern Mastersizer S, the Coulter LS 200, the Fritsch Analysette 22 (version C), the Horiba Partica LA‐950, the Sedigraph 5100, the Atterberg cylinder, the Coulter Multisizer 3, the CIS‐100, the EyeTech and the image analysis software Histolab. The first four instruments use laser diffraction. The Sedigraph and Atterberg cylinder are based on sedimentation. The Coulter Multisizer 3 uses impedance measurements and the last three instruments use optical techniques. Grain‐size analyses were carried out on four sediments with median grain diameters of approximately 35, 30, 12 and 9 μm, respectively, and with no particles >90 μm. The laser diffraction instruments produced the best results for the various criteria considered in this study, followed by the Coulter Multisizer 3 (electrical sensing zone technique), the sedimentation techniques and the optical methods. However, the study also showed that there is no such thing as an ‘ideal’ method because techniques may score differently according to the criterion tested. An appropriate evaluation thus requires that various criteria be considered. Issues to be taken into account are: type of sediment, quantity of sediment available, speed of measurement, complexity of the measurement protocol, processing of the data, reproducibility of the results and the specific aims of the study undertaken.     

Dynamic recrystallization near the brittle-plastic transition in naturally and experimentally deformed quartz aggregates

The electron backscattering diffraction technique (EBSD) was used to analyze bulging recrystallization microstructures from naturally and experimentally deformed quartz aggregates, both of which are characterized by porphyroclasts with finely serrated grain boundaries and grain boundary bulges set in a matrix of very fine recrystallized grains. For the Tonale mylonites we investigated, a temperature range of 300–380 °C, 0.25 GPa confining pressure, a flow stress range of ~ 0.1–0.2 GPa, and a strain rate of ~ 10^− 13 s^− 1 were estimated. Experimental samples of Black Hills quartzite were analyzed, which had been deformed in axial compression at 700 °C, 1.2–1.5 GPa confining pressure, a flow stress of ~ 0.3–0.4 GPa, a strain rate of ~ 10^− 6 s^− 1, and to 44% to 73% axial shortening. Using orientation imaging we investigated the dynamic recrystallization microstructures and discuss which processes may contribute to their development. Our results suggest that several deformation processes are important for the dismantling of the porphyroclasts and the formation of recrystallized grains. Grain boundary bulges are not only formed by local grain boundary migration, but they also display a lattice misorientation indicative of subgrain rotation. Dynamic recrystallization affects especially the rims of host porphyroclasts with a hard orientation, i.e. with an orientation unsuitable for easy basal slip. In addition, Dauphiné twins within porphyroclasts are preferred sites for recrystallization. We interpret large misorientation angles in the experimental samples, which increase with increasing strain, as formed by the activity of fluid-assisted grain boundary sliding.     

实验条件对X射线衍射物相定量分析结果的影响

通过实验验证了多物相试样X射线衍射定量分析的基体效应,着重研究了参比强度K值、衍射强度、混合物各物相粒度、X射线管功率等实验条件对物相定量分析结果的影响。结果表明,物相定性分析是定量分析的前提,采用积分强度、使各物相颗粒尽可能细而且均匀、采用较高的X射线管功率以及适合的扫描速率均有利于提高物相X射线衍射定量分析的准确度。