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.     

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.     

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.     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.