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.     

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.     

Wave runup on a concentric twin perforated circular cylinder

The regular wave interaction with a twin concentric porous circular cylinder system consisting of an inner impermeable cylinder and an outer perforated cylinder was studied through physical model and numerical model studies. The experiments were carried out on the twin concentric cylinder model in a wave flume to study the wave runup and rundown at the leading and trailing edges of the perforated cylinder. It was found that the maximum wave runup on the perforated cylinder is almost same as the incident wave height. The experimental results were used to develop the predictive formulae for the wave runup and rundown on the perforated cylinder, which can be easily used for design applications. The wave runup profiles around the perforated cylinder for different values of ka and porosities were studied numerically using Green's Identity Method. The results of the numerical study are presented and compared with the experimental measurements.     

Threshold velocity for wind erosion: the effects of porous fences

Porous fence is a kind of artificial windbreak that has many practical applications. The threshold wind velocities at different distances downwind from porous fences were measured and the corresponding characteristics of particle movement observed to assess their shelter effect. It is found that the fence’s porosity is the key factor that determines the resulting shelter effect. The area near a fence can be typically classified into five regions, each with a different mode of particle movement. Dense fences, and especially solid fences, favor the accumulation of sand upwind of the fences. Fences with porosities of 0.3–0.4 produce the maximum threshold wind velocity; those with porosities of 0.3–0.6 (depending on the fence height) provide the maximum effective shelter distance. It is confirmed that the fence porosities of 0.3–0.4 that have been proposed for practical application in previous research are the most effective for abating wind erosion.     

Static deformation of two welded monoclinic elastic half-spaces due to a long inclined strike-slip fault

Static deformation of two monoclinic elastic half-spaces in welded contact due to a long inclined strike-slip fault situated in one of the half-spaces is studied analytically and numerically. Closed-form algebraic expressions for the displacement at any point of the medium are obtained. The variation of the displacement at the interface with the horizontal distance from the fault is studied. The effect of anisotropy on the displacement field is examined. It is found that while the anisotropy of the source half-space has a significant effect on the displacement at the interface, the anisotropy of the other half-space has only a marginal effect.     

Surface vibration due to a sequence of high speed moving harmonic rectangular loads

The transmission of vibrations over the surface of the ground, due to high-speed moving, vertical harmonic rectangular loads, is investigated theoretically. The problem is three-dimensional and the interior of the ground is modelled as an elastic half-space or a multilayered ground. The transformed solutions are obtained using the Fourier transform on the space variable. A new damping model in the spatial wavenumber domain, presented in Lefeuve-Mesgouez et al. [J. Sound. Vibr. 231 (2000) 1289] is used. Numerical results for the displacements on the surface are presented for loads moving with speeds up to and beyond the Rayleigh wave speed of the half-space.