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 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.     

Fluid production and isograd reactions at contacts of carbonate-rich and carbonate-poor layers during progressive metamorphism

ABSTRACT With increasing temperature during prograde metamorphism reactions will occur first at the lithological contacts of mixed pelite and calcsilicate terranes. At these interfaces, a fluid of lower chemical potential of H₂O and CO₂ than that required to produce a fluid in either layer can be produced whether reaction is caused by fluid infiltration or is initially fluid absent. If the interface region does not allow fluid transport then as temperature increases, a fluid pressure greater than lithostatic can develop. At some degree of over-pressure relative to rock pressure, the fluid hydraulically fractures the rock and a gradient in fluid composition away from the contact can be produced. These phenomena occur at the compositional interfaces whenever univariant reactions in the differing layers cross on a temperature vs. mole fraction of CO₂ diagram with slopes of opposite sign. The first occurrence of these reaction products at lithological contacts delineates an isograd that defines temperature as well as the mole fraction of CO₂ at constant pressure in systems open to fluid transport. These isograds can be contrasted with fluid-producing isograds in closed systems. As an illustration of possible effects, the reactions quartz + clinozoisite + muscovite = anorthite + K-feldspar + H₂O and phlogopite + quartz + calcite = tremolite + K-feldspar + H₂O + CO₂ at 4 kbar are analysed and equations for fluid production and transport are developed.     

The mass ratio of Nova Muscae 1991

We report intermediate resolution H spectroscopy of the black hole candidate Nova Muscae 1991 during quiescence. We classify the companion star as a K3-K4V which contributes 85±6 percent to the total flux from the binary. The photospheric absorption lines are broadened by 106±13 kms⁻¹ with respect to template field stars, leading to a system mass ratio of q =M₁/M₂ = 7.8_−2.0^+3.4. Doppler imaging of the H line shows strong emission coming from the secondary star (EW=3.1±0.6Å) which we associate with chromospheric activity. However, the hot-spot is not detected and this may indicate a lower mass transfer rate than in other X-ray transients of comparable orbital periods. The surface brightness distribution of the accretion disk in H follows a relation I∝R^−1.1, less steep than typically observed in cataclysmic variables. Updated system parameters are also presented.     

Relaxation of the Vela pulsar’s angular velocity following its jumps

The dynamics of the vortex lattice in the inner crust of a neutron star is considered. A general equation of motion is obtained and solved under the assumption that there are regions of pinned and of free vortices. By comparing these solutions with observational data for the Vela pulsar, the relative moments of inertia of regions of relaxation with the corresponding characteristic times are calculated for two model stars with different equations of state. It is shown that the theory can be reconciled with observations of the relaxation of pulsar angular velocity only for model stars with extremely stiff equations of state. Translated from Astrofizika, Vol. 40, No. 1, pp. 67–76, January–March, 1997.     

Near-Infrared Imaging Spectrophotometry Of Comet C/1995 O1 (Hale-Bopp) Near Perihelion

We present 1- to 5-μm broadband and CVF images of comet Hale-Bopp taken 1997 February 10.5 UT, 50 days before perihelion. All the images exhibit a nonspherical coma with a bright “ridge” in the direction of the dust tail approximately 10″ from the coma. Synthetic aperture spectrophotometry implies that the optically important grains are of a radius ≤0.4 μm; smallest radius for any comet seen to date. The variation of the integrated surface brightness with radial distance from the coma (ρ) in all the images closely follows the “steady state” ρ⁻¹ model for comet dust ablation (Gehrz and Ney, 1992). The near-infrared colors taken along the dust tail are not constant implying the dust grain properties vary with coma distance. This revised version was published online in July 2006 with corrections to the Cover Date.