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.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

Wolf-Rayet stars

Summary Recent literature on Population I Wolf-Rayet star research extending from the Milky Way to blue compact dwarf galaxies is reviewed, broken down into inventory, basic parameters and galactic distribution, atmospheres, binaries, intrinsic variability, mass loss, enrichment and evolution. Also the incidence of Wolf-Rayet stars with variable non-thermal radio emission, excess X-ray fluxes, and episodic/periodic IR excesses is reviewed. These phenomena appear to be associated with wind-wind interaction in wide long-period WR+OB binaries and with wind-compact object interaction in WR+c binaries, with orbit sizes of the order of magnitude of the WR radio photosphere sizes or larger.     

Evidence of climatic warming in the southern Urals region derived from borehole temperatures and meteorological data

Thirty borehole temperature–depth profiles in the central and southern Urals, Russia were scrutinized for evidence of ground surface temperature histories. We explored two inversion schemes: a simple ramp inversion in which solutions are parameterized in terms of an onset time and magnitude of change and a more sophisticated functional space inverse algorithm in which the functional form of the solution is left unspecified. To enhance and potentially identify latitudinal differences in the ground surface temperature signal, we subdivided the data into three groups based on geographic proximity and simultaneously inverted the borehole temperature–depth logs. The simultaneous inversions highlighted 13 temperature–depth logs that could not both fit a common ground surface temperature history and a priori models within reasonable bounds. Our results confirm that this is an effective way to reduce site-specific noise from an ensemble of boreholes. Each inversion scheme gives comparable results indicating locally variable warming on the order of 1°C starting between 1800 and 1900 AD. Similarly surface air temperature records from 12 nearby meteorological stations exhibit locally variable warming also on the order of 1°C of warming during the 20th century. To explore the degree to which borehole temperatures and surface air temperature (SAT) time series are responding to the same signal, we average the SAT data into the same three groups and used these averages as a forcing function at the Earth's surface to generate synthetic transient temperature profiles. Root mean square (RMS) misfits between these synthetic temperature profiles and averaged temperature–depth profiles are low, suggesting that first-order curvature in borehole temperatures and variations in SAT records are correlated.     

Detecting the first objects in the mid-infrared with the Next Generation Space Telescope

We calculate the expected mid-infrared (MIR) molecular hydrogen line emission from the first objects in the Universe. As a result of their low masses, the stellar feedback from massive stars is able to blow away their gas content and collect it into a cooling shell where H₂ rapidly forms and IR roto-vibrational (as for example the rest-frame 2.12 μm) lines carry away a large fraction (up to 10 per cent) of the explosion energy. The fluxes from these sources are in the range 10⁻²¹–10⁻¹⁷ erg s⁻¹ cm⁻² . The highest number counts are expected in the 20-μm band, where about 10⁵ sources deg⁻² are predicted at the limiting flux of 3×10⁻¹⁸ erg s⁻¹ cm⁻². Among the planned observational facilities, we find that the best detection perspectives are offered by the Next Generation Space Telescope ( NGST ), which should be able to reveal about 200 first objects in one hour observation time at its limiting flux in the above band. Therefore, mid-IR instruments appear to represent perfect tools to trace star formation and stellar feedback in the high ( z ≳5) redshift Universe.     

About Darcy's law in non-Galilean frame

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.