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.     

The exact theory of inductive conductivity sensors for oceanographic application

The results of a general theoretical investigation of three commonly used types of inductive conductivity sensors, i.e., the single transformer, the double transformer, and the double transformer with an additional loop, are presented. The resulting formulas describe the dependence of the sensor output signal not only on the conductivity of the seawater but also on the parameters of the electrical circuit, among them the permeability of the transformer core(s), which-unlike the other parameters-shifts considerably during oceanographic in situ measurements. A mathematical discussion of these formulas shows that for certain circuit configurations, the sensor output is independent of changes in permeability. Most of these configurations form the basis of existing oceanographical conductivity sensors, among them the "classical" sensors developed by H. Hinkelmann [3], [4], and by N. L. Brown [14], while some others make evident further possibilities for eliminating the unwanted effects of shifting permeability. In the era of microelectronics, the latter might lead to a reassessment, especially of the single transformer-type sensor.     

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.     

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.     

Zircon U–Pb age, trace element and Hf isotope composition of Kongling terrane in the Yangtze Craton: refining the timing of Palaeoproterozoic high-grade metamorphism

U–Pb age, trace element and Hf isotope compositions of zircon were analysed for a metasedimentary rock and two amphibolites from the Kongling terrane in the northern part of the Yangtze Craton. The zircon shows distinct morphological and chemical characteristics. Most zircon in an amphibolite shows oscillatory zoning, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high formation temperature, high trace element contents, clear negative Eu anomaly, as well as HREE-enriched patterns, suggesting that it is igneous. The zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2857 ± 8 Ma, representing the age of the magmatic protolith. The zircon in the other two samples is metamorphic. It has low Th/U ratios, low trace element concentrations, variable HREE contents (33.8 ≥ Lu_N≥2213; 14.7 ≤ Lu_N/Sm_N ≤ 354) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.000030–0.001168). The data indicate that the zircon formed in the presence of garnet and under upper amphibolite facies conditions. The metamorphic zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2010 ± 13 Ma. These results combined with previously obtained Palaeoproterozoic metamorphic ages suggest a c. 2.0 Ga Palaeoproterozoic collisional event in the Yangtze Craton, which may result from the assembly of the supercontinent Columbia. The zircon in two samples yields weighted mean two-stage Hf model ( T _DM2) ages of 3217 ± 110 and 2943 ± 50 Ma, respectively, indicating that their protoliths were mainly derived from Archean crust.