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.     

Relict mountain slope deposits of northern Portugal: facies,sedimentogenesis and environmental implications

This paper deals with the formation processes and the palaeoenvironmental significance of relict slope deposits located on the uppermost part of the north Portugal mountains. For this purpose, seven key sites representative of the different lithofacies have been selected and analysed in detail. The data show that three main dynamic processes are responsible for the emplacement of regional fossil slope deposits: runoff, debris flows and dry grain flows. The ubiquity of these processes and the lack of frost‐related features or landforms do not support the existence of severe Pleistocene climates in this part of the lberian Peninsula as postulated by previous work. Pedological data gathered at one of the study sites show that a subalpine environment was probably present at 700–800 m altitude between 29 and 14 kyr. Using data from the Pyrenees Mountains, a 6.5 to 12°C depression in mean annual temperature has been tentatively postulated for this Pleniglacial period. Copyright © 2003 John Wiley & Sons, Ltd.     

LRS Bianchi Type V bulk viscous fluid string dust cosmological model in general relativity

An LRS Bianchi Type V bulk viscous fluid dust distribution string cosmological model in General Relativity is investigated. It has been shown that if coefficient of bulk viscosity (ζ) is proportional to the expansion (θ) in the model then string cosmological model for Bianchi Type V space-time is possible. In absence of bulk viscosity(ζ) i.e. when ζ → 0 then there is no string cosmological model for Bianchi Type V space-time. The physical and geometrical aspects are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Plane Symmetric Domain Wall in Lyra Geometry

In this paper general solutions are found for domain walls in Lyra geometry in the plane symmetric spacetime metric given by Taub. Expressions for the energy density and pressure of domain walls are derived in both cases of uniform and time varying displacement field β. It is also shown that the results obtained by Rahaman et al [IJMPD, 10, 735 (2001)] are particular case of our solutions. Finally, the geodesic equations and acceleration of the test particle are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

Abundance and isotopic composition of noble gases in metal and graphite of the Bohumilitz IAB iron meteorite

Abstract— Abundances and isotopic compositions of noble gases in metal and graphite of the Bohumilitz IAB iron meteorite were measured. The abundance ratios of spallogenic components in metal reveal a ³He deficiency which is due to the diffusive loss of parent isotopes, that is, tritium (Tilles, 1963; Schultz, 1967). The diffusive loss likely has been induced by thermal heating by the Sun during cosmic‐ray exposure (～160 Ma; Lavielle et al, 1999). Thermal process such as impact‐induced partial loss may have affected the isotopic composition of spallogenic Ne. The ¹²⁹Xe/¹³¹Xe ratio of cosmogenic components in the metal indicates an enhanced production of epi‐thermal neutrons. The abundance ratios of spallogenic components in the graphite reveal that it contained small amounts of metal and silicates. The isotopic composition of heavy noble gases in graphite itself was obtained from graphite treated with HF/HCl. The isotopic composition of the etched graphite shows that it contains two types of primordial Xe (i.e., Q‐Xe and El Taco Xe). The isotopic heterogeneity preserved in the Bohumilitz graphite indicates that the Bohumilitz graphite did not experience any high‐temperature event and, consequently, must have been emplaced into the metal at subsolidus temperatures. This situation is incompatible with an igneous model as well as the impact melting models for the IAB‐IIICD iron meteorites as proposed by Choi et al. (1995) and Wasson et al (1980).     

Radar observations of Itokawa in 2004 and improved shape estimation

Abstract— We present June 2004 radar images of asteroid 25143 Itokawa (1998 SF36) that improve upon the longitude‐latitude coverage of images obtained in 2001 by Ostro et al. (2004) and use the 2001–2004 data to refine that paper's constraints on Itokawa's shape. The 2004 images, the first of the asteroid's southern side, look distinctly different from the 2001 images, revealing leading edges that are much more curved and rugged than the nearly convex leading edges seen at northern latitudes in 2001. Itokawa is shaped like a slightly asymmetrical, bent, lumpy ellipsoid with dimensions along the principal axes within 10% of 594 times 320 times 288 m. To illustrate the uncertainty space associated with shape reconstruction from images with suboptimal orientational coverage, we present two alternative three‐dimensional models of the object.     

The onset of shock wave in an electrically conducting and radiating gas

The present work applies the method of characteristics to study the behaviour of planar and cylindrical wave-heads propagating through a perfectly electrically conducting and thermally radiating inviscid gas under the optically thin limit in the presence of a transverse magnetic field. The true nonlinear progress of the flow variable gradients at the wavefront is predicted and the critical distance at which the characteristics pile up at the wavefront to form a shock wave is obtained. It is investigated as to how the effects of radiative flux, the magnetic field strength and the specific heat ratio influence the process of steepening or flattening of the characteristic wavefront.