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.     

Axial convergence in a well-mixed estuary

Transverse secondary circulations involving surface convergence, observed in a well-mixed estuary in North Wales, are made visible by the collection of surface material along an axial line which extends continuously for many kilometres through the estuary. The circulation and axial convergence, however, are seen only during the flood phase of the tide and no similar behaviour has been observed during the ebb phase.Convergent circulations in the estuary are associated with small but steady transverse density gradients in the cross-section, produced by non-uniform advection of the longitudinal gradient through the channel. A diagnostic model, using measured mean distributions of cross-sectional density, indicates surface transverse velocities (～0.1 ms^?1) similar to those observed in the estuary. The model further predicts appreciable transverse divergent currents at a fractional depth of 0.75: a prediction which has been tested in the estuary using a vertical array of accurately resolving current direction indicators.     

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.     

On spatial characteristics of five-minute oscillations in the sunspot umbra

Using a differential method we have carried out observations of oscillations in six sunspots. Spectral lines Fe i 5434 Å and Fe i 5576 Å were used. Horizontal waves are not observed in the sunspot umbra photosphere. Results obtained indicate that, at least, the sunspot umbra oscillates as a single whole.     

A WKBJ spectral method for computation of SV synthetic seismograms in a cylindrically symmetric medium

Summary. A method of synthetic seismogram computation for teleseismic SV -waves is developed in order to treat quantitatively SV -waves in problems of body wave source inversion and source—receiver structure studies. The method employs WKBJ theory for a generalized ray in a vertically inhomogeneous half-space and the propagator matrix technique for waves in near-surface homogeneous layers. Wavenumber integration is done along the real axis of the wavenumber plane and anelasticity is included by using complex velocity in all regions of the earth model. The near-surface source structure is taken into account in the computation for the case of the shallow source by allowing a point source to be located in the homogeneous layers. Source and receiver area structures are also allowed to differ. A general moment tensor point source is considered.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite

Two distinct age estimates for eclogite-facies metamorphism in the Sanbagawa belt have been proposed: (i) c.  120–110 Ma based on a zircon SHRIMP age for the Western Iratsu unit and (ii) c.  88–89 Ma based on a garnet–omphacite Lu–Hf isochron age from the Seba and Kotsu eclogite units. Despite the contrasting estimates of formation ages, petrological studies suggest the formation conditions of the Western Iratsu unit are indistinguishable from those of the other two units—all ∼20 kbar and 600–650 °C. Studies of the associated geological structures suggest the Seba and Western Iratsu units are parts of a larger semi-continuous eclogite unit. A combination of geochronological and petrological studies for the Western Iratsu eclogite offers a resolution to this discrepancy in age estimates. New Lu–Hf dating for the Western Iratsu eclogite yields an age of 115.9 ± 0.5 Ma that is compatible with the zircon SHRIMP age. However, petrological studies show that there was significant garnet growth in the Western Iratsu eclogite before eclogite facies metamorphism, and the early core growth is associated with a strong concentration of Lu. Pre-eclogite facies garnet (Grt1) includes epidote–amphibolite facies parageneses equilibrated at 550–650 °C and ∼10 kbar, and this is overgrown by prograde eclogite facies garnet (Grt2). The Lu–Hf age of c.  116 Ma is strongly skewed to the isotopic composition of Grt1 and is interpreted to reflect the age of the pre-eclogite phase. The considerable time gap ( c.  27 Myr) between the two Lu–Hf ages suggests they may be related to separate tectonic events or distinct phases in the evolution of the Sanbagawa subduction zone.