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 comparison of theoretical Sv emission line strengths with extreme ultraviolet observations of a sunspot

Electron impact excitation rates for transitions in the S v ion, calculated with theR-matrix code, are used to derive the electron temperature sensitive emission line ratiosR ₁ =I(854.8 Å)/I(786.9 Å),R ₂ =I(852.2 Å)/I(786.9 Å),R ₃ =I(849.2 Å)/I(786.9 Å), andR ₄ =I(1199.1 Å)/I(786.9 Å), which are found to be significantly different from previous estimates. A comparison of the present results with observational data for a sunspot obtained with the Harvard S-055 spectrometer on boardSkylab reveals generally good agreement between theory and experiment, except in the case ofR ₁, which is probably due to blending in the 854.8 Å feature. The possible effects of Lyman continuum absorption on the observed line ratios is briefly discussed.     

The dust coma of periodic comet Churyumov-Gerasimenko (1982 VIII)

The dust coma of Comet P/Churyumov-Gerasimenko was monitored in the infrared (1–20 μm) from September 1982 to March 1983. Maximum dust production rate of ～2 × 10⁵ g/sec occured in December, 1 month postperihelion. The ratio of dust/gas production was higher than that in other short-period comets. No silicate feature was visible in the 8- to 13-μm spectrum on 23 October. The mean geometric albedo of the grains was ～0.04 at 1.25 μm and ～0.05 at 2.2 μm.     

Merging of low-mass systems and the origin of the fundamental plane

We present the preliminary results of a study of how small stellar systems merge to form larger ones. As we display the families of galaxies in the μ_e - R_e plane (effective surface brightness versus effective radius) we realize that different morphological types occupy different loci, evidencing the different physical mechanisms operating in each family. As proposed by Capaccioli et al. (1992) this diagram is the logical equivalent of the HR diagram for stars. Here we take some initial steps in understanding of how we can establish the evolutionary tracks, solely due to dynamical processes, in the μ_e - R_e plane, ultimately making a dwarf elliptical to turn into a normal elliptical galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

Isotope fractionation of neon during stepheating extraction?: a comment on ‘Re‐interpretation of the existence of a primitive plume under Australia based on neon isotope fractionation during step heating’ by Gautheron and Moreira (2003)

Using a Rayleigh distillation fractionation model, we calculate that the maximum isotope fractionation potentially achievable is less than 5% during the early stages of gas release from a sample. Our calculation corrects the erroneous conclusions of Gautheron and Moreira (2003), who re‐interpreted the plume‐like neon isotopic compositions found in metasomatic apatite from a south‐eastern Australian xenolith (Matsumoto et al., 1997) to be the result of Rayleigh‐type isotope fractionation of originally MORB‐type neon during stepheating gas extraction. We stress that the modelling of neon isotopic fractionation by Gautheron and Moreira (2003) is incorrect, and that the finding of a plume‐like neon isotopic composition in the apatite by Matsumoto et al. (1997) remains a quite valid and robust conclusion.     

Chandrasekhar's perturbation method-oriented theories up to third order for uniformly and differentially-rotating polytropic stars: Error-removing techniques

In this paper we analyse error-removing techniques used without sufficient theoretical support in a previous paper, where Chandrasekhar's higher-order perturbation theories were developed for either uniformly or differentially-rotating polytropic stars.