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.     

Flows in clumpy planetary nebulae

The dynamics of clumps observed in planetary nebulae are considered. The possibility that SiO maser spots in evolved stars and the planetary nebula clumps are formed by the Parker instability behind shocks in pulsating stars' atmospheres is raised. Molecular observations of the clumps are suggested. The effects of the ablation of clumps on the global flow structure of a more tenuous plasma in which they are embedded are reviewed.     

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.     

Hydrogeological framework and water balance studies in parts of Krishni–Yamuna interstream area, Western Uttar Pradesh, India

The Krishni–Yamuna interstream area is a micro-watershed in the Central Ganga Plain and a highly fertile track of Western Uttar Pradesh. The Sugarcane and wheat are the major crops of the area. Aquifers of Quaternary age form the major source of Irrigation and municipal water supplies. A detailed hydrogeological investigation was carried out in the study area with an objective to assess aquifer framework, groundwater quality and its resource potential. The hydrogeological cross section reveals occurrence of alternate layers of clay and sand. Aquifer broadly behaves as a single bodied aquifer down to the depth of 100 m bgl (metre below ground level) as the clay layers laterally pinch out. The depth to water in the area varies between 5 and 16.5 m bgl. The general groundwater flow direction is from NE to SW with few local variations. An attempt has been made to evaluate groundwater resources of the area. The water budget method focuses on the various components contributing to groundwater flow and groundwater storage changes. Changes in ground water storage can be attributed to rainfall recharge, irrigation return flow and ground water inflow to the basin minus baseflow (ground water discharge to streams or springs), evapotranspiration from ground water, pumping and ground water outflow from the basin. The recharge is obtained in the study area using Water table fluctuation and Tritium methods. The results of water balance study show that the total recharge in to the interstream region is of the order of 185.25 million m³ and discharge from the study area is of the order of 203.24 million m³, leaving a deficit balance of −17.99 million m³. Therefore, the present status of groundwater development in the present study area has acquired the declining trend. Thus, the hydrogeological analysis and water balance studies shows that the groundwater development has attained a critical state in the region.     

Risk assessment strategies for the reactivation of earth flows in the Northern Apennines (Italy)

The aim of this paper is to analyze the reactivation mechanism of ancient earth flows, with a view to gleaning information that can subsequently be utilized to formulate a risk-reduction strategy. All considerations made herein are the result of direct experience and observation of actual events which have occurred over the past few decades in the Northern Apennines. Particular attention has been paid to the analysis of the evolution of landslides during actual reactivation, acknowledging a typical, recurring succession of events that precede the failure of the slope. The hazard assessment of these large landslide bodies, which are of slope scale, constitutes a thorny problem, especially in view of the inapplicability of traditional deterministic models such as limit equilibrium stability analysis. Nevertheless, a site-specific assessment of probability of reactivation of these large and ancient earth flows is fundamental to effective land-use planning.     

Deep-seated gravitational slope deformations and their influence on consequent mass movements (case studies from the highest part of the Czech Carpathians)

Gravitational spreading of mountain ridges displays primary disequilibrium of flysch mountain areas of the Czech Carpathians. The progression of various types of mass movements is a product of long-term ridge disintegration and is predisposed by the geological structure of the area and the upper Tertiary-Quaternary morphogenesis of the mountain area. Deep-seated slope deformations are spatially interconnected by the occurrence of some other types of slope deformations (e.g. debris flows, debris slides, slumps, rock avalanches, etc.), which pose a considerable risk for the existence of human society. An important causative factor in these dynamically developing hazardous processes is, among other factors, the way in which land has been used in the last three centuries. Therefore, the occurrence of various types of slope deformations is studied in terms of their relation to deep-seated gravitational deformations and in terms of other limiting factors (structural geological, morphological and climatic factors, manmade impacts, etc.). The paper presents several case studies of slope deformations (Velká Čantoryje Mt, Lysá hora Mt, Ropice Mt and Smrk Mt) in the area of the Outer Carpathians within the territory of the Czech Republic and also adverts to some consequences in terms of the socioeconomic structure of the landscape.     

Rain-induced debris and mudflow triggering factors assessment in the Santiago cordilleran foothills,Central Chile

Debris and mudflows are some of the main geological hazards in the mountain foothills of the Chilean capital city of Santiago. There, the risk of flows triggered in the basins of ravines that drain the range into the city increases with time due to the city growth. A multivariate statistical study based on the logistic regression method is presented. The model provides equations that allow the computation of combined meteorological triggering factors associated with a probability of rain-induced flow occurrence. Daily rainfall, accumulated rainfall and the snowfall level, traditionally considered as the relevant factors, are analysed for a 25-year period. The results show a strong relevance of the rainfall on the day of the flow event over the other factors. However, the relatively low probabilities returned for some real flow events suggest that the model does not capture all the significant variables and the problem is more complex than as it has been traditionally assumed, and further investigations are needed to develop predictive models of flow triggering.     

The spatial implications of the functional proximity deriving from air passenger flows between European metropolitan urban regions

Until recently the traditional spatial configuration of the European geography was based upon the core-periphery model. The ‘pentagon’, broadly defined as lying between London, Paris, Milan, Munich and Hamburg, was seen as the core area characterised by having the highest concentration of economic development in the European Union (EU), with the remainder of the European territory viewed as peripheral, albeit to varying degrees. In a number of cases such peripheral areas equated with clear regional disparities. The elaboration of the European Spatial Development Perspective (ESDP) (CEC, European spatial development perspective, towards balanced and sustainable development in the territory of the European Union, 1999) challenged this core-periphery model. European spatial planning policies, aimed at encouraging social and economic, and with ever increasing importance, territorial cohesion, seek amongst other aspects to encourage the development of a balanced and polycentric urban system. This paper adopts a network analysis approach to the analysis of air passenger flows between some 28 principal European metropolitan urban regions. The evaluation of these flows contributes to an enhanced comprehension of the spatial dynamics of the European metropolitan territory which goes beyond that deriving from the more standard analyses of the individual components of the urban system. Several indicators are used, deriving from gravitational modelling techniques, to analyse the complexity of the air passenger flows. A multidimensional scaling (MDS) technique is introduced in order to interpret and visualise the resulting spatial configuration and positioning of the different metropolitan centres within the conceptual European ‘space of air passenger flows’, thereby contrasting with the more traditional map-based geographical image of Europe, based upon Cartesian coordinates.