A characterization of the coupled evolution of grain fabric and pore space using complex networks: Pore connectivity and optimized flows in the presence of shear bands

A framework for the multiscale characterization of the coupled evolution of the solid grain fabric and its associated pore space in dense granular media is developed. In this framework, a pseudo-dual graph transformation of the grain contact network produces a graph of pores which can be readily interpreted as a pore space network. Survivability, a new metric succinctly summarizing the connectivity of the solid grain and pore space networks, measures material robustness. The size distribution and the connectivity of pores can be characterized quantitatively through various network properties. Assortativity characterizes the pore space with respect to the parity of the number of particles enclosing the pore. Multiscale clusters of odd parity versus even parity contact cycles alternate spatially along the shear band: these represent, respectively, local jamming and unjamming regions that continually switch positions in time throughout the failure regime. Optimal paths, established using network shortest paths in favor of large pores, provide clues on preferential paths for interstitial matter transport. In systems with higher rolling resistance at contacts, less tortuous shortest paths thread through larger pores in shear bands. Notably the structural patterns uncovered in the pore space suggest that more robust models of interstitial pore flow through deforming granular systems require a proper consideration of the evolution of in situ shear band and fracture patterns – not just globally, but also inside these localized failure zones.     

Study of the effect of capillary pressure on the permeability of porous media embedded with a fractal-like tree network

In this paper, the analytical expressions for permeability of (both saturated and unsaturated) porous media embedded with a fractal-like tree network are presented based on fractal theory and technique when the capillary pressure is taken into account. Both the dimensionless effective permeability and the relative permeability of the composites, which are defined as porous media embedded with a fractal-like tree network in this work, are derived and found to be a function of saturation, the capillary pressure and microstructural parameters of the networks. The relative permeabilities predicted by the present fractal model are compared with the available experimental data and a fair agreement between them is found.     

Predicting the Permeability of Very Low Porosity Sandstones

The permeability predictions of two geometric pore-scale models, one being predominantly granular and the other consolidated with tube-like pores, are compared with experimental results for Fontainebleau sandstones and the results interpreted. Percolation thresholds are determined from experimental data and applied in the modelling exercise by means of cut-off asymptotes on porosity. It is found that, although both granular and foamlike models yield plausible results, the granular model appears to be superior, at least for the sets of data considered. The Klinkenberg correction is analytically derived and incorporated into the models to relate gas and liquid permeabilities and an analytical expression for the Klinkenberg factor is proposed for each model. The permeability predictions are promising and yield an effective manner to correlate sandstone percolation data.     

Characterization of bubble mobility in channel flow with fibrous porous media walls

In composites processing, resin is introduced into a fibrous domain to cover all the empty spaces between the fibers. It is important to extract air bubbles from the domain before the resin solidifies. Failure to do so will entrap these voids in the final part, which is detrimental to its performance. Hence, there is a need to understand bubble motion in a fibrous porous domain in which the bubbles move with the resin in channels surrounded by fibrous walls. A rising bubble model is presented that consists of a single spherical void in a cylindrical axisymmetric two-phase domain of resin and air surrounded by porous media boundaries. The motion of a bubble in a channel flow with porous boundaries is modeled by replacing the walls with a slip velocity. Focus is on how the porous media permeability influences the bubble motion. A parameter called bubble mobility is defined as the ratio of bubble rise velocity to the resin free surface velocity. Results suggest that fabric permeability and fluid properties can be optimized to increase bubble mobility and ultimately lead to reduction in void content during composites processing.     

Series solution for steady flow of a third grade fluid through porous space

In this paper steady flow of a third grade fluid through porous space is considered. Modified Darcy’s law for third grade fluid in a porous space has been introduced. The governing non-linear equation is first modelled and then solved using homotopy analysis method (HAM). The convergence of the obtained series solution is discussed. The effects of the emerging parameters on the velocity field are seen. It is noted that meaningful solution exists only in the case of suction.     

Rotating flow of a third grade fluid in a porous space with Hall current

This study investigates the rotating magnetohydrodynamic (MHD) flow of a third-grade fluid in a porous space. Modified Darcy's law has been utilized for the flow modeling. The Hall effects are taken into consideration. The basic equations governing the flow are reduced to a highly nonlinear ordinary differential equation. This equation has been solved analytically by employing the homotopy analysis method (HAM). The effects of the various interesting parameters on the velocity distribution have been discussed.     

Numerical simulation of non-Darcian flow through a porous medium

This work reports on fluid flow in a fluid-saturated porous medium, accounting for the boundary and inertial effects in the momentum equation. The flow is simulated by Brinkman-Forchheimer-extended Darcy formulation (DFB), using MAC (Marker And Cell) and Chorin pressure iteration method. The method is validated by comparison with analytic results. The effect of Reynolds number, Darcy number, porosity and viscosity ratio on velocity is investigated. As a result, it is found that Darcy number has a decisive influence on pressure as well as velocity, and the effect of viscosity ratio on velocity is very strong given the Darcy number. Additional key findings include unreasonable choice of effective viscosity can involve loss of important physical information.     

Permeability of microscale fibrous porous media using the lattice Boltzmann method

The permeabilities of microscale fibrous porous media were calculated using the multiple-relaxation-time (MRT) lattice Boltzmann method (LBM). Two models of the microscale fibrous porous media were constructed based on overlapping fibers (simple cubic, body-centered cubic). Arranging the fibers in skew positions yielded two additional models comprising non-overlapping fibers (skewed simple cubic, skewed body-centered cubic). As the fiber diameter increased, the fibers acted as granular inclusions. The effects of the overlapping fibers on the media permeability were investigated. The overlapping fibers yielded permeability values that were a factor of 2.5 larger than those obtained from non-overlapping fibers, but the effects of the fiber arrangement were negligible. Two correlations were obtained for the overlapping and non-overlapping fiber models, respectively. The effects of the rarefaction and slip flow are also discussed. As the Knudsen number increased, the dimensionless permeability increased; however, the increase differed depending on the fiber arrangement. In the slip flow regime, the fiber arrangement inside the porous media became an important factor.     

Poromechanical Modelling and Inverse Approach of Drying Tests on Weakly Permeable Porous Rocks

The present paper deals with the determination of permeability in partially saturated conditions for weakly permeable porous rocks such as argillites or deep clayey formations. The level of permeability can be obtained via the measurements of transient weight loss of a sample submitted to a decrease in relative humidity imposed by saline solution in a hermetic chamber. An identification method based on simplified uncoupled 1D-linear and 1D-non-linear modelling was presented in a previous paper (Giraud et al. Trans Porous Media 69(2):259–280, 2006). The present paper takes into account generalized mass transfer phenomena such as Darcean advective transport of liquid and gas mixtures and Fickean diffusive transport of the vapour specie inside a gas mixture. Poromechanical coupling as well as 3D effects due to the geometry and finite dimensions of the tested samples are also covered by this approach. The coupled THM finite element computer code Code_Aster is then used to model the forward problem. The parameter identification procedure is based upon the solution of an inverse problem. The Levenberg–Marquardt algorithm was used for the problem of minimization. Comparisons between previous simplified 1D modelling and 2D-axisymmetrical coupled modelling show that the former method efficiently provides the correct order of magnitude of the level of permeability or the equivalent storage coefficients. Due to the boundary condition, the real 2D-axisymmetrical geometry of the sample must not be neglected if we are to obtain accurate results.     

Flow of oil–water emulsions through a constricted capillary

The flow of oil-in-water emulsions through quartz micro-capillary tubes was analyzed experimentally. The capillaries were used as models of connecting pore-throats between adjacent pore body pairs in high-permeability media. Pressure drop between the inlet and outlet ends of the capillary was recorded as a function of time, for several values of the volumetric flow rate. Several distinct emulsions were prepared using synthetic oils in deionized water, stabilized by a surfactant (Triton X-100). Two oils of different viscosity values were used to prepare the emulsions, while two distinct drop size distributions were obtained by varying the mixing procedure. The average oil drop size varied from smaller to larger than the neck radius. The results are presented in terms of the extra-pressure drop due to the presence of the dispersed phase, i.e. the difference between the measured pressure drop and the one necessary to drive the continuous phase alone at the same flow rate. For emulsions with drops smaller than the capillary throat diameter, the extra-pressure drop does not vary with capillary number and it is a function of the viscosity ratio, dispersed phase concentration and drop size distribution. For emulsions with drops larger than the constriction, the large oil drops may partially block the capillary, leading to a high extra pressure difference at low capillary numbers. Changes in the local fluid mobility by means of pore-throat blockage may help to explain the additional oil recovery observed in laboratory experiments and the sparse data on field trials.     

Permeability identification of a weakly permeable partially saturated porous rock

The present paper deals with the determination of permeability in partially saturated conditions for weakly permeable porous continua such as argillites or deep clayey formations. The permeability can be deduced from measurements of transient weight loss of a sample submitted to a laboratory drying test: a decrease of relative humidity is imposed by saline solution in an hermetic chamber. Assumptions of constant gas pressure equal to atmospheric pressure and of negligible Fickean diffusive transport of vapour are adopted. The only transport phenomenon taken into account inside the sample is the Darcean advective transport of the water liquid. The forward problem is solved by following two modelling approaches: a linear one and a nonlinear one. The parameter identification procedure is based upon the solution of corresponding inverse problems. In the two cases, the Levenberg–Marquardt algorithm has been used for the minimization problem. In the linear approach, the solution of the forward problem is explicit. In the non linear approach, finite volume method for the spatial discretization combined with a Newton–Raphson algorithm has been used to solve the non linear forward problem. The identification method enables variations of permeability and capillary capacity to be estimated. Comparisons between linear and non linear approaches show that the first one is useful to give mean values and order of magnitude of permeability and capacity. A more complete information is deduced from the non linear approach as variations of equivalent capacity and permeability during a test are significant in most cases. The analysis of the obtained results shows that the basic modelling assumption of constant gas pressure inside the sample would not be relevant for lower range of relative humidities and liquid permeability than those investigated.     

Gas cluster growth by solute diffusion in porous media. Experiments and automaton simulation on pore network

The study presented in this paper deals with the liquid–gas phase change by pressure decline of supersaturated CO₂ solutions in 2D porous media. The growth of the gas phase is studied experimentally and numerically as a function of supersaturation, wettability and gravity. Experiments are performed on a transparent etched network (micromodel) and simulations with a specific numerical automaton.In the experiments, the nucleation process, i.e. the occurrence of the gas bubbles, as well as the growth of these bubbles are visualised and analysed by means of a micro video camera and an image processing apparatus. The observations confirm the heterogeneous nature of nucleation and the disordered growth pattern of the gas phase. The analysis of the growth rate of a single gas cluster shows that this phenomenon is different from the compact growth of an isolated single bubble in the bulk. As previously predicted, the bubble growth by mass transfer and volume expansion in porous media is characterised by a pattern of the invasion percolation type under normal laboratory conditions.Numerical simulations of the growth pattern and the growth rate of a single gas cluster are performed with a numerical automaton. Based on a pore network modelling technique and a set of hypotheses derived from the observations, this automaton is first validated by comparing the numerical results with the experiments. Then, the automaton is used to conduct a sensitivity study. In particular, the influences of the Jakob number, pressure decline rate, Bond number, wettability and characteristics of the microstructure are investigated.     

A Finite Element-based Heuristic Estimation of Local Preform Permeability for Resin Transfer Molding

Uniformity of fabrics significantly affects the resin flow behavior in the resin transfer molding (RTM) process. Due to fabric defects or improper fiber preform preparation/loading, non-uniformity in fabric structure frequently occurs in RTM processing and creates local permeability variations. Such variations often lead to unbalanced resin flow patterns and thus result in defects of finished composite parts. In RTM process modeling, an accurate estimation of the whole field permeability profile of the fiber preform is critical for predicting resin flow pattern correctly. In this article, a finite element-based heuristic computing method is introduced for estimating the in situ whole-field isotropic permeability profile of the preform using a steady flow of gas. Compared with conventional approaches, this method is effective in measuring local permeability variations and applicable to molds with complex 2-D geometries, as well as diverse injection strategies. Several case studies were presented with experimental designs and numerical computations to illustrate the effectiveness and efficiency of the method.     

Combustion in a horizontal channel partially filled with a porous media

Experiments were carried out to investigate the combustion propagation phenomenon in a horizontal channel partially filled with ceramic-oxide spherical beads. A 1.22 m long, 43 mm nominally thick layer of spherical beads is located at the ignition end of a 2.44 m long, 76 mm square channel. Tests were performed with 6.4 and 12.7 mm diameter beads. A flame is ignited at the bead end wall by an automotive spark ignition system. Flame propagation and pressure measurements are obtained via ionization probes and piezoelectric pressure transducers mounted on the top and bottom surfaces of the channel. High-speed schlieren video was used to visualize the structure of the explosion front. Experiments were performed with a 31% nitrogen diluted stoichiometric methane–oxygen mixture at room temperature and at an initial pressure in the range of 15–50 kPa. For initial pressures of 15 and 20 kPa the flame accelerates to a velocity close to the speed of sound in the combustion products. For initial pressure of 30 kPa and higher DDT occurs in the gap above the bead layer. An explosion front propagating at a velocity just under the CJ detonation velocity is detected in the bead layer even though the bead layer pore size is much smaller than the detonation cell size. It is demonstrated that flame propagation within the bead layer is the driving force behind the very rapid flame acceleration observed, however the DDT event occurring in the gap above the bead layer is not affected by the bead layer porosity. Schlieren video indicates that the structure of the explosion front varies across the channel height and with propagation distance down the channel.     

Revisiting poromechanics in the context of microporous materials

Poromechanics offers a consistent theoretical framework for describing the mechanical response of porous solids, fully or partially saturated with a fluid phase. When dealing with fully saturated microporous materials, which exhibit pores of the nanometre size, aside from the fluid pressure acting on the pore walls additional effects due to adsorption and confinement of the fluid molecules in the smallest pores must be accounted for. From the mechanical point of view, these phenomena result into volumetric deformations of the porous solid: the so-called “swelling” phenomenon. The present work investigates how the poromechanical theory should be refined in order to describe adsorption and confinement induced swelling in microporous solids. Firstly, we report molecular simulation results that show that the pressure and density of the fluid in the smallest pores are responsible for the volumetric deformation of the material. Secondly, poromechanics is revisited in the context of a microporous material with a continuous pore size distribution. Accounting for the thermodynamic equilibrium of the fluid phase in the overall pore space, the new formulation introduces an apparent porosity and an interaction free energy. We use a prototype constitutive relation relating these two quantities to the Gibbs adsorption isotherm, and then calculate the induced deformation of the solid matrix. Agreement with experimental data found in the literature is observed. As an illustrating example, we show the predicted strains in the case of adsorption of methane on activated carbon.     

考虑摩擦特性的柔性关节空间机器人自适应CMAC神经网络鲁棒控制

讨论了关节摩擦力矩影响下,具有柔性铰关节的漂浮基空间机器人系统的动力学控制问题.设计了基于高斯基函数的小脑神经网络(CMAC)鲁棒控制器和摩擦力矩补偿器.用奇异摄动理论对系统的动力学模型进行快慢变子系统分解,针对快变子系统,设计力矩微分反馈控制器来抑制机械臂关节柔性引起的振动;对于慢变子系统,设计了基于自适应CMAC神...     

On the Non-Isentropic Sound Waves Propagation in a Cylindrical Tube Filled with Saturated Porous Media

A rigid frame, cylindrical capillary theory of sound propagation in porous media that includes the nonlinear effects of the Forchheimer type is laid out by using variational solutions. It is shown that the five main parameters governing the propagation of sound waves in a fluid contained in rigid cylindrical tubes filled with a saturated porous media are: the shear wave number, , the reduced frequency parameter, , the porosity, ε, Darcy number, , and Forchheimer number, . The manner in which the flow influences the attenuation and the phase velocities of the forward and backward propagating non-isentropic acoustic waves is deduced. It is found that the inclusion of the solid matrix increases wave’s attenuations and phase velocities for both forward and backward sound waves, while increasing the porosity and the reduced frequency number decreased attenuation and increased phase velocities. The effect of the steady flow is found to decrease the attenuation and phase velocities for forward sound waves, and enhance them for the backward sound waves. This work is done during a sabbatical leave year granted form the University of Jordan to Dr. Hamzeh Duwairi for the academic year 2007/2008 at the German Jordanian University.     

Numerical and experimental investigation of convective drying in unsaturated porous media with bound water

The heat and mass transfer in an unsaturated wet cylindrical porous bed packed with quartz particles was investigated theoretically for relatively low convective drying rates. Local thermodynamic equilibrium was assumed in the mathematical model describing the multi-phase flow in the unsaturated porous media using the energy and mass conservation equations to describe the heat and mass transfer during the drying. The drying model included convection and capillary transport of the free water, diffusion of bound water, and convection and diffusion of the gas. The numerical results indicated that the drying process could be divided into three periods, the temperature rise period, the constant drying rate period and the decreasing drying rate period. The numerical results agreed well with the experimental data verifying that the mathematical model can evaluate the drying performance of porous media for low drying rates. The effects of drying conditions such as the ambient temperature, the relative humidity, and the velocity of the drying air, on the drying process were evaluated by numerical solution.