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.     

Thick and thin models of the evolution of carbon dioxide on Mars

We present results from a new simulation code that accounts for the evolution of the reservoirs of carbon dioxide on Mars, from its early years to the present. We establish a baseline model parameter set that produces results compatible with the present (i.e., Patm?6.5 mbar with permanent CO₂ ice cap) for a wide range of initial inventories. We find that the initial inventory of CO₂ broadly determines the evolutionary course of the reservoirs of CO₂. The reservoirs include the atmosphere, ice cap, adsorbed CO₂ in the regolith, and carbonate rocks. We track the evolution of the free inventory: the atmosphere, ice cap and regolith. Simulations begin at 4.53 Gyr before present with a rapid loss of free inventory to space in the early Noachian. Models that assume a relatively small initial inventory (?5 bar) have pronounced minima in the free inventory of CO₂ toward the end of the Noachian. Under baseline parameters, initial inventories below ∼4.5 bar result in a catastrophic loss of the free inventory to space. The current free inventory would be then determined by the balance between outgassing, sputtering losses and chemical weathering following the end of the late bombardment. We call these “thin” models. They generically predict small current free inventories in line with expectations of a small present CO₂ ice cap. For “thick” models, with initial inventories ?5 bar, a surplus of 300-700 mbar of free CO₂ remains during the late-Noachian. The histories of free inventory in time for thick models tend to converge within the last 3.5 Gyr toward a present with an ice cap plus atmospheric inventory of about 100 mbar. For thick models, the convergence is largely due to the effects of chemical weathering, which draws down higher free inventories more rapidly than the low. Thus, thick models have ?450 mbar carbonate reservoirs, while thin models have ?200 mbar. Though both thick and thin scenarios can reproduce the current atmospheric pressure, the thick models imply a relatively large current CO₂ ice cap and thin models, little or none. While the sublimation of a massive cap at a high obliquity would create a climate swing of greenhouse warming for thick models, under the thin model, mean temperatures and pressures would be essentially unaffected by increases in obliquity.     

Exact analysis of the combined data of SNO and Super-Kamiokande

Comparison of solar-neutrino signals in SNO [Phys. Rev. Lett. 87 (2001) 071301] and Super-Kamiokande (SK) [Phys. Rev. Lett. 86 (2001) 5651] detectors results in discovery of ν_e→ν_μ,τ oscillations at level 3.1–3.3σ [Phys. Rev. Lett. 87 (2001) 071301]. This comparison involves the assumption of neutrino spectrum and a choice for the thresholds of detection in both experiments. In this paper we obtain an exact formula for the comparison of the signals which is valid for arbitrary spectra and thresholds. We find that the no-oscillation hypothesis (astrophysical solutions) is excluded at 3.3σ. If the energy-dependent component of the survival probability for electron neutrinos is small as compared with the average value, i.e. in the case of small distortion of the observed spectrum, the oscillation hypothesis can also be tested to similar accuracy. The oscillation to sterile neutrino only, is excluded at 3.3σ level, and oscillation to active neutrinos is confirmed at 2.8σ.     

Central regions of the early-type galaxies in the NGC 3169 group

We have investigated the central regions of the galaxies in the NGC 3169/NGC 3166/NGC 3156 group with the multipupil fiber spectrograph of the 6-m telescope; the first (central) galaxy in the group is a spiral (Sa) one and the other two galaxies are lenticular ones. The group is known to have an extended HI cloud with a size of more than 100 kpc that is associated in its position, orientation, and rotation with the central galaxy NGC 3169. The mean age of the stellar populations in the centers of all three galaxies has been found to be approximately the same, ～1 Gyr. Since the galaxies are early-type ones and since NGC 3166 and NGC 3156 show no global star formation, we are dealing here with a synchronous star formation burst in the centers of all three galaxies.     

Size-Flux Relation in Solar Active Regions

We present a study of the relationship between integral area and corresponding total magnetic flux for solar active regions. It is shown that some of these relationships are satisfied to simple power laws. Fractal examination showed that some of these power laws can not be justified inside the simple models of stationary magnetic flux tube aggregation. All magnetic fluxes and corresponding areas were calculated using the data measured with the Solar Magnetic Field Telescope of the Huairou Solar Observing Station in Beijing.     

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.     

A mechanism of formation of globular clusters in the corona of a giant elliptical galaxy

Leningrad State University. Translated from Astrofizika, Vol. 32, No. 2, pp. 267–275, March–April, 1990.     

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