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.     

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.     

Can the big-bang singularity be avoided in the scale-covariant theory?

The Raychaudhuri equation is written in terms of atomic units in the scale-covariant theory of Canutoet al. (1977) and it is pointed out that even dust-filled nonrotating cosmological models may not have an observable singularity of infinite density in the beginning. Simple explicit solutions of Einstein's equations are also constructed in terms of atomic units and its is shown that the Einstein-de Sitter solution in gravitational units may appear to be a static or an oscillating universe in atomic units.     

Range of validity of seismic ray and beam methods in general inhomogeneous media – I. General theory

Summary. The limitations of asymptotic wave theory and its geometrical manifestations are newly formalized and scrutinized. Necessary and sufficient conditions for the existence of acoustic and seismic rays and beams in general inhomogeneous media are expressed in terms of new physical parameters: the threshold frequency ω₀ associated with the P/S decoupling condition, the cut-off frequency ω_c associated with the radiation-zone condition, the total curvature of the wavefront and the Fresnel-zone radius.
The analysis is facilitated with the introduction of a new ancillary functional – the hypereikonal which is capable of representing ordinary as well as evanescent waves. The hypereikonal is the natural extension of the eikonal theory.
With the aid of the above new parameters, simple conditions are obtained for the decoupled far field, the decoupled near field, two point dynamic ray tracing, paraxial wavefields and Gaussian beams.     

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.