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.     

MATLAB在运用系统建模处理南海气象数据过程中应用

介绍了MATLAB语言特点和系统建模方法的基本理论.根据南海气象数据的实际建模处理过程,给出了建模的详细步骤及其MATLAB实现过程以及MATLABTM的主要实现程序.试验讨论和结果表明利用MATLAB语言可以方便地对南海气象数据用系统建模方法进行建模和处理,MATLAB在运用系统建模法处理南海气象数据方面具有明显的优越性.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

LU Vel (GJ 375): A M3.5Ve Flare and Double-Lined Spectroscopic Binary

High resolution echelle spectroscopic observations taken with the FEROS spectrograph at the 2.2 m telescope ESO confirm the binary nature of the flare M3.5V star LU Vel (GJ 375, RE J0958-462) previously reported by Christian and Mathioudakis (2002). Emission of similar intensity from both components is detected in the Balmer, Na i D₁&D₂, He i D₃, Ca ii H&K, and Ca ii IRT lines. We have determined precise radial velocities by cross correlation with radial velocity standard stars, which have allowed us to obtain for the first time the orbital solution of the system. The binary consists of two near-equal M3.5V components with an orbital period shorter than 2 days. We have analyzed the behaviour of the chromospheric activity indicators (variability and possible flares). In addition, we have determined its rotational velocity and kinematics.     

OSIRIS: AO-assisted integral-field spectroscopy at the Keck Observatory

OSIRIS (OH-Suppressing Infra-Red Integral-field Spectrograph) is a new facility instrument for the Keck Observatory. After seeing first light in February 2005, OSIRIS is currently undergoing commissioning. OSIRIS provides the capability of performing three-dimensional spectroscopy in the near-infrared z, J, H, and K bands at the resolution limit of the Keck II telescope, which is equipped with adaptive optics and a laser guide star. The science case for OSIRIS is summarized, and the instrument and associated data reduction software are described.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

An atlas of the optical spectrum of supernova 1987A in the large magellanic cloud

Photographic spectra of SN1987A in the LMC have been obtained from 1987 February 25 to 1988 June 30. Microdensitometer tracings of these have been reduced to intensity and corrections for instrumental response have been applied to the spectra. This paper presents these data in an atlas format, discusses the reduction procedures in detail, and presents radial velocity measurements of selected lines in the spectra     

利用自回归模型检测地球极移重力效应

同时采用4个台站的国际超导重力仪长期连续重力观测资料和国际地球自转服务中心提供的同步地球自转参数，研究了极移引起的地球重力场变化特征。利用自回归模型估计了各序列的功率谱密度和积谱密度，结果表明极移导致的重力效应的主要能量集中在Chandler摆动和周年项附近，叠积后实际重力观测与极移重力信号理论值之间的差异分别为0．4％和3．9％，说明超导重力仪可有效监测极移导致的重力变化。     

多条太阳光球FeI线在磁场中形成及形成深度的综合分析(Ⅰ)

针对多通道滤光器太阳磁场望远镜的磁场观测定标及掌握(所采用的)谱线特征之需要,取VAL宁静太阳大气模型计算了7条Fel光球线的Stokes轮廓、形成深度、贡献函数分布,从而较为系统地对多条反常及正常Zeeman线的特征及性质作出了分析与总结,解释了在磁光效应影响下Stokes Q,U参量的形成深度曲线在近线心区域处出现陡峭峰值的原因。     

Modelling Gas And Dust Release From Comet Hale–Bopp

Numerical simulations of the evolving activity of comet Hale-Bopp are presented, assuming a porous, spherical nucleus, 20 km in radius, made of dust and gas-laden amorphous ice. The main effects included are: crystallization of amorphous ice and release of occluded gas, condensation, sublimation and flow of gases through the pores, changing pore sizes, and flow of dust grains. The model parameters, such as initial pore size and porosity, emissivity, dust grain size, are varied in order to match the observed activity. In all cases, a sharp rise in the activity of the nucleus occurs at a large heliocentric distance pre-perihelion, marked by a few orders of magnitude increase in the CO and the CO₂ fluxes and in the rate of dust emission. This is due to the onset of crystallization, advancing down to a few meters below the surface, accompanied by release of the trapped gases. A period of sustained, but variable, activity ensues. The emission of water molecules is found to surpass that of CO at a heliocentric distance of 3 AU. Thereafter the activity is largely determined by the behaviour of the dust. If a dust mantle is allowed to build up, the water production rate does not increase dramatically towards perihelion; if most of the dust is ejected, the surface activity increases rapidly, producing a very bright comet. This revised version was published online in July 2006 with corrections to the Cover Date.