Numerical Approach of the Permeability of a Macroporous Bioceramic with Interconnected Spherical Pores

Manufacturing a hybrid bone substitute requires a dynamic culture of the cells preliminarily seeded in a scaffold through a flow of physiological fluid. The velocity, pressure, and the distribution of fluid flow in this kind of macroporous medium are the important keys. Because of the difficulties in determining these parameters by experiment, a numerical approach has been chosen. One of the primary step of this study consists in the determination of permeability K. In this article, two types of structure of macroporous bioceramics are concerned. One is the interconnected pore spheres arranged either simple cubic, body-centered cubic or face-centered cubic systems. The other is the interconnected pore spheres randomly arranged. Based on Darcy??s law, the permeability K was calculated for many cases (type, porosity) by simulating the fluid flow through a small representative volume. These results are compared with some previous models such as Ergun, Carman?CKozeny, Rumpf?CGupte, and Du Plessis. The limits of Darcy??s law and the above-mentioned models have been determined using numerical simulation. The result showed that the porous media with spherical interconnected pores of BCC systems can be used to replace a complex random system in a range of porosity from 0.71 to 0.76 (i.e., porosity of our scaffolds). This assumption is validated for a pressure gradient lower the 1,000?Pa m^?C1 and a simple polynomial relation linking permeability and porosity (0.71?C0.76) has been established.     

Grain Shape Effects on Permeability,Formation Factor,and Capillary Pressure from Pore-Scale Modeling

We invoke pore-scale models to evaluate grain shape effects on petrophysical properties of three-dimensional (3D) images from micro-CT scans and consolidated grain packs. Four sets of grain-packs are constructed on the basis of a new sedimentary algorithm with the following shapes: exact angular grain shapes identified from micro-CT scans, ellipsoids fitted to angular grains, and spheres with volume and surface-to-volume ratio equal to original angular grains on a grain-by-grain basis. Subsequently, a geometry-based cementation algorithm implements pore space alteration due to diagenesis. Eight micro-CT scans and 144 grain-pack images with $500 \times 500 \times 500$ voxels (the resolution units of 3D images) are analyzed in this study. Absolute permeability, formation factor, and capillary pressure are calculated for each 3D image using numerical methods and compared to available core measurements. Angular grain packs give rise to the best agreement with experimental measurements. Cement volume and its spatial distribution in the pore space significantly affect all calculated petrophysical properties. Available empirical permeability correlations for non-spherical grains underestimate permeability between 30 and 70 % for the analyzed samples. Kozeny–Carman’s predictions agree with modeled permeability for spherical grain packs but overestimate permeability for micro-CT images and non-spherical grain packs when volume-based radii are used to calculate the average grain size in a pack. We identify surface-to-volume ratio and grain shape as fundamental physical parameters that control fluid distribution and flow in porous media for equivalent porosity samples.     

Remarks on the Rayleigh-Bénard Convection on Spherical Shells

The main objective of this article is to study the effect of spherical geometry on dynamic transitions and pattern formation for the Rayleigh-Bénard convection. The study is mainly motivated by the importance of spherical geometry and convection in geophysical flows. It is shown in particular that the system always undergoes a continuous (Type-I) transition to a 2l _c -dimensional sphere ${S^{2l_c}}$ , where l _c is the critical wave number corresponding to the critical Rayleigh number. Furthermore, it has shown in Ma and Wang (Physica D 239:3–4, 167–189, 2010) that it is critical to add nonisotropic turbulent friction terms in the momentum equation to capture the large-scale atmospheric and oceanic circulation patterns. We show in particular that the system with turbulent friction terms added undergoes the same type of dynamic transition, and obtain an explicit formula linking the critical wave number (pattern selection), the aspect ratio, and the ratio between the horizontal and vertical turbulent friction coefficients.     

Connectivity of Pore Space as a Control on Two-Phase Flow Properties of Tight-Gas Sandstones

We predict capillary-pressure (drainage) curves in tight-gas sandstones which have little matrix or microporosity using a quantitative grain-scale model. The model accounts for the geometric results of some depositional and diagenetic processes important for porosity and permeability reduction in tight-gas sandstones, such as deformation of ductile grains during burial and quartz cementation. The model represents the original sediment as a dense, disordered packing of spheres. We simulated the evolution of this model sediment into a low-porosity sandstone by applying different amounts of ductile grains and quartz precipitation. A substantial fraction of original pore throats in the sediment is closed by the simulated diagenetic alteration. Because the percolation threshold corresponds to closure of half of the pore throats, the pore space in this type of tight-gas sandstone is poorly connected and is often close to being completely disconnected. The drainage curve for different model rocks was computed using invasion percolation in a network taken directly from the grain-scale geometry and topology of the model. Some general trends follow classical expectations and were confirmed by experimental measurements: increasing the amount of cement shifts the drainage curve to larger pressures. This is related to reduction of the connectivity of pore space resulting from closure of throats. Existence of ductile grains in the ductile grain model also reduces the connectivity of pore space but it treats the throats distribution differently causing the drainage curves to be shifted to larger irreducible water saturation when cement is added to the model. The range of porosities in which these connectivity effects are important corresponds to the range of porosities common for tight gas sandstones. Consequently these rocks can exhibit small effective permeability to gas even at large gas saturations. This problem occurs at larger porosities in rocks with significant content of ductile grains because ductile deformation blocks a significant fraction of pore throats even before cementation begins. Predicted drainage curves agree with measurements on two samples with little microporosity, one dominated by rigid grains, the other containing a significant fraction of ductile grains. We conclude that connectivity of the matrix pore space is an important factor for an understanding of flow properties of tight-gas sandstones.     

Poisson–Nernst–Planck Systems for Ion Flow with Density Functional Theory for Hard-Sphere Potential: I–V Relations and Critical Potentials. Part II: Numerics

We consider a one-dimensional steady-state Poisson–Nernst–Planck type model for ionic flow through membrane channels. Improving the classical Poisson–Nernst–Planck models where ion species are treated as point charges, this model includes ionic interaction due to finite sizes of ion species modeled by hard sphere potential from the Density Functional Theory. The resulting problem is a singularly perturbed boundary value problem of an integro-differential system. We examine the problem and investigate the ion size effect on the current–voltage (I–V) relations numerically, focusing on the case where two oppositely charged ion species are involved and only the hard sphere components of the excess chemical potentials are included. Two numerical tasks are conducted. The first one is a numerical approach of solving the boundary value problem and obtaining I–V curves. This is accomplished through a numerical implementation of the analytical strategy introduced by Ji and Liu in [Poisson–Nernst–Planck systems for ion flow with density functional theory for hard-sphere potential: I–V relations and critical potentials. Part I: Analysis, J. Dyn. Differ. Equ. (to appear)]. The second task is to numerically detect two critical potential values V _c and V ^c .The existence of these two critical values is first realized for a relatively simple setting and analytical approximations of V _c and V ^c are obtained in the above mentioned reference. We propose an algorithm for numerical detection of V _c and V ^c without using any analytical formulas but based on the defining properties and numerical I–V curves directly. For the setting in the above mentioned reference, our numerical values for V _c and V ^c agree well with the analytical predictions. For a setting including a nonzero permanent charge in which case no analytic formula for the I–V relation is available now, our algorithms can still be applied to find V _c and V ^c numerically.     

Quantifying the dynamics of flow within a permeable bed using time-resolved endoscopic particle imaging velocimetry (EPIV)

This paper presents results of an experimental study investigating the mean and temporal evolution of flow within the pore space of a packed bed overlain by a free-surface flow. Data were collected by an endoscopic PIV (EPIV) technique. EPIV allows the instantaneous velocity field within the pore space to be quantified at a high spatio-temporal resolution, thus permitting investigation of the structure of turbulent subsurface flow produced by a high Reynolds number freestream flow (Re _s in the range 9.8?×?10³?C9.7?×?10⁴). Evolution of coherent flow structures within the pore space is shown to be driven by jet flow, with the interaction of this jet with the pore flow generating distinct coherent flow structures. The effects of freestream water depth, Reynolds and Froude numbers are investigated.     

Two-Phase Flow Properties Prediction from Small-Scale Data Using Pore-Network Modeling

The objective of this work is to evaluate the prediction accuracy of network modeling to calculate transport properties of porous media based on the interpretation of mercury invasion capillary pressure curves only. A pore-scale modeling approach is used to model the multi-phase flow and calculate gas/oil relative permeability curves. The characteristics of the 3-D pore-network are defined with the requirement that the network model satisfactorily reproduces the capillary pressure curve (P_c curve), the porosity and the permeability. A sensitivity study on the effect of the input parameters on the prediction of capillary pressure and gas/oil relative permeability curves is presented. The simulations show that different input parameters can lead to similarly good reproductions of the experimental P_c, although the predicted relative permeabilities K_r are somewhat widespread. This means that the information derived from a mercury invasion P_c curve is not sufficient to characterize transport properties of a porous medium. The simulations indicate that more quantitative information on the wall roughness and the node/bond aspect ratio would be necessary to better constrain the problem. There is also evidence that in narrow pore size distributions pore body volume and pore throat radius are correlated while in broad pore size distributions they would be uncorrelated.     

Flow in a spherical cavity due to a stokeslet

The earth is considered as a rigid spherical cavity of radius a filled with a highly viscous incompressible fluid of viscosity η. The non-axisymmetric problem of flow due to a stokeslet of strength F/8πη located at (0, 0, c), c < a with its axis along OX, O being the centre of the sphere, is discussed for small Reynolds numbers. The expressions for velocity and pressure are obtained in terms of A(r, θ, φ) and B(r, θ, φ), biharmonic and harmonic functions respectively, using a sphere theorem for non-axisymmetric flow inside a sphere. The drag on the sphere exerted by the fluid is found to be independent of the location of the stokeslet.     

Slow flow through a periodic array of spheres

Slow flow through a periodic array of spheres is studied theoretically, and the drag force by the fluid on a sphere forming the periodic array is calculated using a modification of the method developed by Hashimoto (1959). Results for the complete range of volume fraction $\text{c}$ of spheres are given for simple cubic, body-centered cubic, and face-centered cubic arrays and these agree well with the corresponding values reported by previous investigators. Also, series expansions for the drag force to 0($\text{c}{}^{10}$) are derived for each of these cubic arrays. The method is also applied to determine the drag force to 0($\text{c}{}^3$) on infinitely long cylinders in square and hexagonal arrays.     

Relative permeability and capillary pressure functions of porous media as related to the displacement growth pattern

Visualization experiments of the unsteady immiscible displacement of a fluid by another are performed on glass-etched pore networks of well-controlled morphology by varying the fluid system and flow conditions. The measured transient responses of the fluid saturation and pressure drop across the porous medium are introduced into numerical solvers of the macroscopic two-phase flow equations to estimate the non-wetting phase, k_rnw, and wetting phase, k_rw, relative permeability curves and capillary pressure, P_c, curve. The correlation of k_rnw, k_rw, and P_c with the displacement growth pattern is investigated. Except for the capillary number, wettability, and viscosity ratio, the immiscible displacement growth pattern in a porous medium may be governed by the shear-thinning rheology of the injected or displaced fluid, and the porous sample length as compared to the thickness of the frontal region. The imbibition k_rnw increases as the flow pattern changes from compact displacement to viscous fingering or from viscous to capillary fingering. The imbibition k_rw increases as the flow pattern changes from compact displacement or capillary fingering to viscous fingering. As the shear-thinning behaviour of the NWP strengthens and/or the contact angle decreases, then the flow pattern is gradually dominated by irregular interfacial configurations, and the imbibition k_rnw increases. The imbibition P_c is a decreasing function of the capillary number or increasing function of the injected phase viscosity in agreement with the linear thermodynamic theory.     

On the modeling of miscible flow in multi-component porous media

An analytical model of miscible flow in multi-component porous media is presented to demonstrate the influence of pore capacitance in extending diffusive tailing. Solute attenuation is represented naturally by accommodating diffusive and convective flux components in macropores amd micropores as elicited by the local solute concentration and velocity fields. A set of twin, coupled differential equations result from the Laplace transform and are solved simultaneously using a differential operator for one-dimensional flow geometry. The solutions in real space are achieved using numeric inversion. In addition, to represent more faithfully the dominant physical processes, this approach enables efficient and stable semi-analytical solution procedure of the coupled system that is significantly more complex than current capacitance type models. Parametric studies are completed to illustrate the ability of the model to represent sharp breakthrough and lengthy tailing, as well as investigating the form of the nested heterogeneity as a result of solute exchange between macropores and micropores. Data from a laboratory column experiment is examined using the present model and satisfactory agreement results.Roman Letters a rate coefficient of internal flow - b velocity ratio (v ₁/v ₂) - h dispersion ratio (D ₂/D ₁) - c ₁ macropore concentration - c ₂ micropore concentration - ¯c ₁ macropore concentration in Laplace space - ¯c ₂ micropore concentration in Laplace space - c ₁ ⁰ macropore concentration at source location - c ₂ ⁰ micropore concentration at source location - D ₁ macropore dispersion coefficient - D ₂ micropore dispersion coefficient - f fraction of pore space occupied by fluid in primary channel - L length of laboratory sample column - K mass exchange rate - t time from initial stage - v ₁ primary flow channel velocity - v ₂ micropore interstitial velocity - x distance from source - y dimensionless distance Greek Letters equivalent Péclet number - dimensionless time, or injected pore volume     

Connectivity in Vuggy Carbonates, New Experimental Methods and Applications

The length and spatial distribution of the touching-vugs channels affect the degree of permeability variations and is the main contributor to heterogeneity in vuggy carbonates. Hence, this article focuses on vug connectivity characterization and its impact on fluid flow. A whole core sample was scanned by X-ray computed tomography (CT). Image segmentation was used to obtain a binarized three-dimensional (3D) model of the vuggy pore space. Analysis of the binarized 3D model is used to calculate the correlation function and correlation length for the vuggy pore space. Connectivity analysis of the binarized 3D model shows that 79% of the vugs connected network spanning along the sample. The remaining 21% vug porosity exists in a large number of isolated vugs. The correlation length for the connected vug network is found to be larger than for vugs in general. NMR T ₂ measurements at increasing capillary pressure is tested on the vuggy material and used to investigate the amount of connected- and isolated vugs. The results verify the large fraction of connected vugs. Application of NMR T ₂ measurements in combination with capillary pressure experiments can also reveal matrix properties that play an important role in recovery processes. The transition between non-Fickian and Fickian regimes for tracer/solute transport is studied by laboratory experiments performed at various sample lengths, from cm to m scale. For the largest sample measured in our experiments show that effluent concentration curve conform to the CDE solution, suggesting that the Fickian regime has been established.     

Stability of viscoelastic flow around periodic arrays of cylinders

Low Reynolds number flow of Newtonian and viscoelastic Boger fluids past periodic square arrays of cylinders with a porosity of 0.45 and 0.86 has been studied. Pressure drop measurements along the flow direction as a function of flow rate as well as flow visualization has been performed to investigate the effect of fluid elasticity on stability of this class of flows. It has been shown that below a critical Weissenberg number (We_c), the flow in both porosity cells is a two-dimensional steady flow, however, pressure fluctuations appear above We_c which is 2.95±0.25 for the 0.45 porosity cell and 0.95±0.08 for the higher porosity cell. Specifically, in the low porosity cell as the Weissenberg number is increased above We_c a transition between a steady two-dimensional to a transient three-dimensional flow occurs. However, in the high porosity cell a transition between a steady two-dimensional to a steady three-dimensional flow consisting of periodic cellular structures along the length of the cylinder in the space between the first and the second cylinder occurs while past the second cylinder another transition to a transient three-dimensional flow occurs giving rise to time- dependent cellular structures of various wavelengths along the length of the cylinder. Overall, the experiments indicate that viscoelastic flow past periodic arrays of cylinders of various porosities is susceptible to purely elastic instabilities. Moreover, the instability observed in lower porosity cells where a vortex is present between the cylinders in the base flow is amplifieds spatially, that is energy from the mean flow is continuously transferred to the disturbance flow along the flow direction. This instability gives rise to a rapid increase in flow resistance. In higher porosity cells where a vortex between the cylinders is not present in the base flow, the energy associated with the disturbance flow is not greatly changed along the flow direction past the second cylinder. In addition, it has been shown that in both flow cells the instability is a sensitive function of the relaxation time of the fluid. Hence, the instability in this class of flows is a strong function of the base flow kinematics (i.e., curvature of streamlines near solid surfaces), We and the relaxation time of the fluid.     

Effective elasticity tensor of a periodic composite

The effective elasticity tensor of a composite is defined to be the four-tensor C which relates the average stress to the average strain. We determine it for an array of rigid spheres centered on the points of a periodic lattice in a homogeneous isotropic elastic medium. We first express C in terms of the traction exerted on a single sphere by the medium, and then derive an integral equation for this traction. We solve this equation numerically for simple, body-centered and face-centered cubic lattices with inclusion concentrations up to 90% of the close-packing concentration. For lattices with cubic symmetry the effective elasticity tensor involves just three parameters, which we compute from the solution for the traction. We obtain approximate asymptotic formulas for low concentrations which agree well with the numerical results. We also derive asymptotic results for C at high inclusion concentrations for arbitrary lattice geometries. We find them to be in good agreement with the numerical results for cubic lattices. For low and moderate concentrations the approximate results of Nemat-Nasseret al., also agree well with the numerical results for cubic lattices.     

Extensional rheology of concentrated emulsions as probed by capillary breakup elongational rheometry (CaBER)

Elongational flow behavior of w/o emulsions has been investigated using a capillary breakup elongational rheometer (CaBER) equipped with an advanced image processing system allowing for precise assessment of the full filament shape. The transient neck diameter D(t), time evolution of the neck curvature κ(t), the region of deformation l _def and the filament lifetime t _c are extracted in order to characterize non-uniform filament thinning. Effects of disperse volume fraction ?, droplet size d _sv , and continuous phase viscosity η _c on the flow properties have been investigated. At a critical volume fraction ? _c , strong shear thinning, and an apparent shear yield stress τ _y,s occur and shear flow curves are well described by a Herschel–Bulkley model. In CaBER filaments exhibit sharp necking and t _c as well as κ _max ?=?κ (t?=?t _c ) increase, whereas l _def decreases drastically with increasing ?. For ? <?? _c , D(t) data can be described by a power-law model based on a cylindrical filament approximation using the exponent n and consistency index k from shear experiments. For ??≥?? _c , D(t) data are fitted using a one-dimensional Herschel–Bulkley approach, but k and τ _y,s progressively deviate from shear results as ? increases. We attribute this to the failure of the cylindrical filament assumption. Filament lifetime is proportional to η _c at all ?. Above ? _c, κ _max as well as t _c /η _c scale linearly with τ _y,s . The Laplace pressure at the critical stretch ratio ε _c which is needed to induce capillary thinning can be identified as the elongational yield stress τ _y,e , if the experimental parameters are chosen such that the axial curvature of the filament profile can be neglected. This is a unique and robust method to determine this quantity for soft matter with τ _y ?< 1,000 Pa. For the emulsion series investigated here a ratio τ _y,e /τ _y,s = 2.8 ± 0.4 is found independent of ?. This result is captured by a generalized Herschel–Bulkley model including the third invariant of the strain-rate tensor proposed here for the first time, which implies that τ _y,e and τ _y,s are independent material parameters.     

Permeability of Pore Networks Under Compaction

The characteristic pore length fixes the scale of permeability of a porous medium. For pore networks undergoing strong random compaction, this length becomes singular at transition porosities, revealing a change in the microstructure of the porespace controlling the transport. Nodal balances and lattice Boltzmann simulations of flow in pore networks under compaction show that the scaling between permeability and porosity changes near the transition porosities. Simulation results are compared with experimental permeability data from transparent two-dimensional micromodels of networks decorated with the same pore size distribution. Permeability?Cporosity data of media undergoing smooth compaction is well described by a single power law. Under strong compaction, however, the scaling between permeability and porosity is possible by traits only, the scaling exponent changes notably at given transition porosities. These behaviors are common to a wealth of permeability?Cporosity data thus far unexplained.     

On the Domain of Validity of Brinkman’s Equation

An increasing number of articles are adopting Brinkman’s equation in place of Darcy’s law for describing flow in porous media. That poses the question of the respective domains of validity of both laws, as well as the question of the value of the effective viscosity μ ^e which is present in Brinkman’s equation. These two topics are addressed in this article, mainly by a priori estimates and by recalling existing analyses. Three main classes of porous media can be distinguished: “classical” porous media with a connected solid structure where the pore surface S _p is a function of the characteristic pore size l _p (such as for cylindrical pores), swarms of low concentration fixed particles where the pore surface is a function of the characteristic particle size l _s , and fiber-made porous media at low solid concentration where the pore surface is a function of the fiber diameter. If Brinkman’s 3D flow equation is valid to describe the flow of a Newtonian fluid through a swarm of fixed particles or fibrous media at low concentration under very precise conditions (Lévy 1983), then we show that it cannot apply to the flow of such a fluid through classical porous media.     

A Numerical Model of Tracer Transport in a Non-isothermal Two-Phase Flow System for {\text{ CO}}_2 Geological Storage Characterization

For the purpose of characterizing geologically stored $\text{ CO}_{2}Air sparging is an in situ soil/groundwater remediation technology, which involves the injection of pressurized air through air sparging well below the zone of contamination. To investigate the rate-dependent flow properties during multistep air sparging, a rule-based dynamic two-phase flow model was developed and applied to a 3D pore network which is employed to characterize the void structure of porous media. The simulated dynamic two-phase flow at the pore scale or microscale was translated into functional relationships at the continuum-scale of capillary pressure?Csaturation (P ^c?CS) and relative permeability??saturation (K _r?CS) relationships. A significant contribution from the air injection pressure step and duration time of each air injection pressure on both of the above relationships was observed during the multistep air sparging tests. It is observed from the simulation that at a given matric potential, larger amount of water is retained during transient flow than that during steady flow. Shorter the duration of each air injection pressure step, there is higher fraction of retained water. The relative air/water permeability values are also greatly affected by the pressure step. With large air injection pressure step, the air/water relative permeability is much higher than that with a smaller air injection pressure step at the same water saturation level. However, the impact of pressure step on relative permeability is not consistent for flows with different capillary numbers (N _ca). When compared with relative air permeability, relative water permeability has a higher scatter. It was further observed that the dynamic effects on the relative permeability curve are more apparent for networks with larger pore sizes than that with smaller pore sizes. In addition, the effect of pore size on relative water permeability is higher than that on relative air permeability.     

Estimation of Shale Intrinsic Permeability with Process-Based Pore Network Modeling Approach

A multi-scale pore network model is developed for shale with the process-based method (PBM). The pore network comprises three types of sub-networks: the \(\upmu \)m-scale sub-network, the nm-scale pore sub-network in organic matter (OM) particles and the nm-scale pore sub-network in clay aggregates. Process-based simulations mimic shale-forming geological processes and generate a \(\upmu \)m-scale sub-network which connects interparticle pores, OM particles and clay aggregates. The nm-scale pore sub-networks in OM and clay are extracted from monodisperse sphere packing. Nm-scale throats in OM and clay are simplified to be cylindrical and cuboid-shaped, respectively. The nm-scale pore sub-networks are inserted into selected OM particles and clay aggregates in the \(\upmu \)m-scale sub-network to form an integrated multi-scale pore network. No-slip permeability is evaluated on multi-scale pore networks. Permeability calculations verify that shales permeability keeps decreasing when nm-scale pores and throats replace \(\upmu \)m-scale pores. Soft shales may have higher porosity but similar range of permeability with hard shales. Small compaction leads to higher permeability when nm-scale pores dominate a pore network. Nm-scale pore networks with higher interconnectivity contribute to higher permeability. Under constant shale porosity, the shale matrix with cuboid-shaped nm-scale throats has lower no-slip permeability than that with cylindrical throats. Different from previous reconstruction processes, the new reconstruction process first considers the porous OM and clay distribution with PBM. The influence of geological processes on the multi-scale pore networks is also first analyzed for shale. Moreover, this study considers the effect of OM porosities and different pore morphologies in OM and clay on shale permeability.