Search results

The purpose of this study is to validate whether the local thermal equilibrium for unsteady state is an appropriate assumption for the porous media with closed pores. It also compares the transient temperatures between the pore scale and volume averaged approaches to prove that the volume averaged method is an appropriate technique for the heat transfer in closed-cell porous media. The interfacial heat transfer coefficient for the closed-cell porous media is also discussed in details.

The governing equations for the pore scale and continuum domains are given. They are solved numerically for the pore scale and volume-averaged domains. The results are compared and discussion was done. The performed discussions and explanations are supported with figure and graphics.

A local thermal non-equilibrium exits for the closed-cell porous media in which voids are filled with water during the unsteady heat transfer process. Local thermal non-equilibrium condition exists in the cells under high temperature gradient and it disappears when the heat transfer process becomes steady-state. Although a local thermal equilibrium exists in the porous media in which the voids are filled with air, a finite value for heat transfer coefficient is found. The thermal diffusivity of air and solid phase are close to each other and hence a local thermal equilibrium exists.

The study is done only for the closed-cell porous media and for Rayleigh number till 10⁵. Two common working fluids as water and air are considered.

There are many applications of porous media with closed pores particularly in the industry, such as the closed-cell metal foam or the closed cells in porous materials such as foods and plastic-based insulation material. The obtained results are important for transient heat transfer in closed-cell porous materials.

The obtained results are important from the transient application of heat transfer in the closed-cell material existing in nature and industry.

The authors’ literature survey shows that it is the first time the closed-cell porous media is discussed from local thermal non-equilibrium point of view and it is proved that the local thermal non-equilibrium can exist in the closed-cell porous media. Hence, two equations as solid and fluid equations should be used for unsteady heat transfer in a closed-cell porous medium.

We present a numerical methodology for the study of convective pore‐fluid, thermal and mass flow in fluid‐saturated porous rock basins. In particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore‐fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore‐fluid flow and mass transport in fluid‐saturated hydrothermal basins; (2) Convective pore‐fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence on the distribution patterns of convective pore‐fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore‐fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

Porous titanium is used in many bioimplant and electrode applications because of its interconnected pore structure and good corrosion resistance. The purpose of this paper is to study the anodic polarization behavior of porous titanium in different electrolytes and clarify the influences of the porosity and macro‐pore size on the corrosion resistance.

The porous titanium with 10‐70% porosities and average macro‐pore sizes in the range of 100‐500 μm was prepared by the powder metallurgy method using polymethyl methacrylate (PMMA) as a space holder. Electrochemical corrosion tests were performed on porous titanium as well as solid titanium (with the same irregular and isolated micro‐pore structures as that on the interconnected spheroidal macro‐pore walls of porous titanium) in the 0.1 M H₂SO₄, 1 M NaOH and 0.9% NaCl (37 °C) solutions.

It was found that porous titanium exhibited an active‐passive transition behavior in the 1 M NaOH and 0.1 M H₂SO₄ solutions. In contrast, a self‐passivation transition behavior was observed in the 0.9% NaCl solution (37 °C).

The paper demonstrates that both the porosity and macro‐pore size of porous titanium play an important role in determining the corrosion rate, rather than the corrosion potential.

The purpose of this paper is to provide a summarization and review of the present author's main investigations on failure modes of reticular metal foams under different loadings in engineering applications.

With the octahedral structure model proposed by the present authors themselves, the fundamentally mechanical relations have been systematically studied for reticular metal foams with open cells in their previous works. On this basis, such model theory is continually used to investigate the failure mode of this kind of porous materials under compression, bending, torsion and shearing, which are common loading forms in engineering applications.

The pore-strut of metal foams under different compressive loadings will fail in the tensile breaking mode when it is brittle. While it is ductile, it will tend to the shearing failure mode when the shearing strength is half or nearly half of the tensile strength for the corresponding dense material and to the tensile breaking mode when the shearing strength is higher than half of the tensile strength to a certain value. The failure modes of such porous materials under bending, torsional and shearing loads are also similarly related to their material species.

This paper presents a distinctive method to conveniently analyze and estimate the failure mode of metal foams under different loadings in engineering applications.

The knowledge of structural parameters of nonwovens media is poorly understood. The pores size distribution (PSD) function is one of those parameters. The difficulty is not only the understanding of the distribution of pores but also the identification of pores geometry distribution (PGD) and their behaviour concerning the dynamic fluid transportation. The purpose of this paper is to present an efficient and reliable method based on image analysis which on one hand, performs the estimation of the PSD function and takes into account the geometric aspect of pores, and on the other hand, analyses liquid wicking in very thin filter media.

The proposed methods, in this paper, are applied on thin filter media made of polyester. The samples have not sudden any treatment. The authors set up an optical test bed in order to observe the dynamic properties of the samples. Dynamic raw data about the liquid wicking are extracted directly from video sequences using the appropriate test bed. The structural parameters are extracted from the non‐wetted samples.

Obtained results allow a better understanding of the liquid wicking in very thin filter media. In addition to the PSD function, the PGD function adds informations about the shape of pores. The dynamic data of the liquid wicking explains that pores have different behaviour when liquid reached them. It can be deduced from this study that the fluid transport in the pore network is defined by three main parameters: geometric parameter (size, shape), capillary action and pores connection in the network.

The led back‐lighting system is not sufficient to observe precisely the liquid wicking. An additional front‐lighting will be added in further studies.

The extraction of dynamic properties from video sequences, by performing image analysis is an original method to characterise the porosity in thin media filter.

This research study aims to develop regular cylindrical pore network models (RCPNMs) to calculate topology and geometry properties of explosively created fractures along with their resulting hydraulic permeability. The focus of the investigation is to define a method that generates a valid geometric and topologic representation from a computational modelling point of view for explosion-generated fractures in rocks. In particular, extraction of geometries from experimentally validated Eulerian smoothed particle hydrodynamics (ESPH) approach, to avoid restrictions for image-based computational methods.

Three-dimensional stabilized ESPH solution is required to model explosively created fracture networks, and the accuracy of developed ESPH is qualitatively and quantitatively examined against experimental observations for both peak detonation pressures and crack density estimations. SPH simulation domain is segmented to void and solid spaces using a graphical user interface, and the void space of blasted rocks is represented by a regular lattice of spherical pores connected by cylindrical throats. Results produced by the RCPNMs are compared to three pore network extraction algorithms. Thereby, once the accuracy of RCPNMs is confirmed, the absolute permeability of fracture networks is calculated.

The results obtained with RCPNMs method were compared with three pore network extraction algorithms and computational fluid dynamics method, achieving a more computational efficiency regarding to CPU cost and a better geometry and topology relationship identification, in all the cases studied. Furthermore, a reliable topology data that does not have image-based pore network limitations, and the effect of topological disorder on the computed absolute permeability is minor. However, further research is necessary to improve the interpretation of real pore systems for explosively created fracture networks.

Although only laboratory cylindrical rock specimens were tested in the computational examples, the developed approaches are applicable for field scale and complex pore network grids with arbitrary shapes.

It is often desirable to develop an integrated computational method for hydraulic conductivity of explosively created fracture networks which segmentation of fracture networks is not restricted to X-ray images, particularly when topologic and geometric modellings are the crucial parts. This research study provides insight to the reliable computational methods and pore network extraction algorithm selection processes, as well as defining a practical framework for generating reliable topological and geometrical data in a Eulerian SPH setting.

Natural gas leak from underground pipelines could lead to serious damage and global warming, whose spreading in soil should be systematically investigated. This paper aims to propose a three-dimensional numerical model to analyze the methane–air transportation in soil. The results could help understand the diffusion process of natural gas in soil, which is essential for locating leak source and reducing damage after leak accident.

A numerical model using finite element method is proposed to simulate the methane spreading process in porous media after leaking from an underground pipe. Physical models, including fluids transportation in porous media, water evaporation and heat transfer, are taken into account. The numerical results are compared with experimental data to validate the reliability of the simulation model. The effects of methane leaking direction, non-uniform soil porosity, leaking pressure and convective mass transfer coefficient on ground surface are analyzed.

The methane mole fraction distribution in soil is significantly affected by the leaking direction. Horizontally and vertically non-uniform soil porosity has a stronger effect. Increasing leaking pressure causes increasing methane mole flux and flow rate on the ground surface.

Most existing gas diffusion models in porous media are for one- or two-dimensional simulation, which is not enough for predicting three-dimensional diffusion process after natural gas leak in soil. The heat transfer between gas and soil was also neglected by most researchers, which is very important for predicting the gas-spreading process affected by the soil moisture variation because of water evaporation. In this paper, a three-dimensional numerical model is proposed to further analyze the methane–air transportation in soil using finite element method, with the presence of water evaporation and heat transfer in soil.

This study aims to numerically study the buoyant convective flow of two different nanofluids in a porous annular domain. A uniformly heated inner cylinder, cooled outer cylindrical boundary and adiabatic horizontal surfaces are considered because of many industrial applications of this geometry. The analysis also addresses the comparative study of different porous media models governing fluid flow and heat transport.

The finite difference method has been used in the current simulation work to obtain the numerical solution of coupled partial differential equations. In particular, the alternating direction implicit method is used for solving transient equations, and the successive line over relaxation iterative method is used to solve time-independent equation by choosing an optimum value for relaxation parameter. Simpson’s rule is adopted to estimate average Nusselt number involving numerical integration. Various grid sensitivity checks have been performed to assess the sufficiency of grid size to obtain accurate results. In this analysis, a general porous media model has been considered, and a comparative study between three different models has been investigated.

Numerical simulations are performed for different combinations of the control parameters and interesting results are obtained. It has been found that the an increase in Darcy and Rayleigh numbers enhances the thermal transport rate and strengthens the nanofluid movement in porous annulus. Also, higher flow circulation rate and thermal transport has been detected for Darcy model as compared to non-Darcy models. Thermal mixing could be enhanced by considering a non-Darcy model.

The present results could be effectively used in many practical applications under the limiting conditions of two-dimensionality and axi-symmetry conditions. The only drawback of the current study is it does not include the three-dimensional effects.

The results could be used as a first-hand information for the design of any thermal systems. This will help the design engineer to have fewer trial-and-run cases for the new design.

A pioneering numerical investigation on the buoyant convective flow of two different nanofluids in an annular porous domain has been carried out by using a general Darcy–Brinkman–Forchheimer model to govern fluid flow in porous matrix. The results obtained from current investigation are novel and original, with numerous practical applications of nanofluid saturated porous annular enclosure in the modern industry.

This paper aims to numerically investigate the natural convection heat transfer of a hybrid nanofluid into a porous cavity exposed to a variable magnetic field.

The non-linear elliptical governing equations have been solved numerically using control volume based finite element method. The effects of different governing parameters including Rayleigh number (Ra = 10³ − 10⁶), Hartman number (Ha = 0 − 50), volume fraction of nanoparticles (φ = 0 − 0.02), curvature of horizontal isolated wall (a = 0.85 − 1.15), porosity coefficient (ε = 0.1 − 0.9) and Darcy number (Da = 10⁻⁵ − 10⁻¹) have been studied.

The results indicate that at low Darcy numbers close to 0, the average Nusselt number Nu_a enhances as porosity coefficient increases. For a = 1 and a = 1.15 in comparison with a = 0.85, the stretching of the isothermal lines is maintained from the left side to the right side and vice versa, which indicates increased natural convection heat transfer for this configuration of the top and bottom walls. In addition, at higher Rayleigh numbers, by increasing the Hartmann number, a significant decrease is observed in the Nusselt number, which can be attributed to the decreased power of the flow.

The authors believe that all the results, both numerical and asymptotic, are original and have not been published elsewhere.

This paper aims to investigate the behaviour of geosynthetic reinforced deep cement mixed (DCM) column-supported embankments constructed over soft soils.

Coupled consolidation analyses based on the finite element method are carried out assuming that the soil and DCM columns are fully saturated porous mediums. In the first part of the paper, a case study of an embankment constructed over a very soft soil deposit in Finland is presented. Two- and three-dimensional finite element models for the case study are developed including isolated and attached DCM columns beneath the embankment to capture the arching mechanism between DCM columns. The model simulations were carried out considering the actual staged construction procedure adopted in the field. Finite element predictions show good agreement with field data and confirm that the load transfer is mainly between attached columns beneath the embankment. Next, the significance of geosynthetic reinforcement on the load transfer mechanism is investigated. Finally, the influence of permeability of columns and soft soil on the performance of geosynthetic reinforcement column-supported embankments is studied.

Results demonstrate that the excess pore pressure dissipation rate is fast in DCM column-improved ground compared to the same case without any columns, although the same permeability is assigned to both DCM columns and surrounding soft soil. When DCM column permeability exceeds soil permeability, excess pore pressure dissipation rate shows a remarkable increase compared to that observed when the DCM column permeability is less than or equal to the permeability of surrounding soft soil. [ ]

This paper investigates the contribution of permeability and geosynthetic layer on the vertical load transfer mechanism of the embankment and modelling issues related to application of the embankment load and the properties of the cement-improved columns.