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.     

Nature of the Martian uplands: Effect on Martian meteorite age distribution and secondary cratering

Abstract— Martian meteorites (MMs) have been launched from an estimated 5–9 sites on Mars within the last 20 Myr. Some 80–89% of these launch sites sampled igneous rock formations from only the last 29% of Martian time. We hypothesize that this imbalance arises not merely from poor statistics, but because the launch processes are dominated by two main phenomena: first, much of the older Martian surface is inefficient in launching rocks during impacts, and second, the volumetrically enormous reservoir of original cumulate crust enhances launch probability for 4.5 Gyr old rocks. There are four lines of evidence for the first point, not all of equal strength. First, impact theory implies that MM launch is favored by surface exposures of near‐surface coherent rock (≤10² m deep), whereas Noachian surfaces generally should have ≥10² m of loose or weakly cemented regolith with high ice content, reducing efficiency of rock launch. Second, similarly, both Mars Exploration Rovers found sedimentary strata, 1–2 orders of magnitude weaker than Martian igneous rocks, favoring low launch efficiency among some fluvial‐derived Hesperian and Noachian rocks. Even if launched, such rocks may be unrecognized as meteorites on Earth. Third, statistics of MM formation age versus cosmic‐ray exposure (CRE) age weakly suggest that older surfaces may need larger, deeper craters to launch rocks. Fourth, in direct confirmation, one of us (N. G. B.) has found that older surfaces need larger craters to produce secondary impact crater fields (cf. Barlow and Block 2004). In a survey of 200 craters, the smallest Noachian, Hesperian, and Amazonian craters with prominent fields of secondaries have diameters of ?45 km, ?19 km, and ?10 km, respectively. Because 40% of Mars is Noachian, and 74% is either Noachian or Hesperian, the subsurface geologic characteristics of the older areas probably affect statistics of recognized MMs and production rates of secondary crater populations, and the MM and secondary crater statistics may give us clues to those properties.