Packed Glass Beads Research Articles

Influence of stagnant zones on solute transport in heterogeneous porous media at the pore scale

Solute transport in porous media is sensitive to heterogeneity at all scales. However, the pore-scale solute transport behavior may considerably affect the behavior at larger scales. Here, a multi-relaxation-time lattice Boltzmann method with Flekkøy's mass transfer scheme is employed for simulating the fluid flow and solute transport in three-dimensional porous media obtained from high-resolution micro focus x-ray computed tomography, namely, randomly packed glass beads and four consolidated sandstones with an increasing level of heterogeneity, i.e., Fontainebleau, Berea, Takoh, and Shirahama. The flow field heterogeneity is carefully resolved for each porous media in terms of streamlines, Eulerian velocity fields, and the ratio of stagnant zones, which is consistent with the sequence of coordination numbers. Dispersion results show that Fick's law is satisfied well for glass beads, whereas early arrivals and late-time tailings are observed for heterogeneous rocks from the residence time distribution. Then, the dispersion coefficient is calculated using the time moment method, indicating that more heterogeneous porous media exhibit larger dispersion coefficients. The scalar dissipation rate (SDR) is resolved to characterize the mixing state. Two distinctive time regimes are recognized for heterogeneous rocks, separating at around 10 convective time scales. At a later time, a universal power-law scaling of SDR with time is observed, with the power-law γ being approximately 1.5 for glass beads (indicative of Fickian dispersion) and 2–3 for heterogeneous rocks. Finally, the significance of the mass transfer rate between the mobile and stagnant zones on the mechanical dispersion is evaluated in terms of the Damhöhler (Da) number.

Effect of selected sampling media, flow rate, and time on the sampling efficiency of a liquid impinger packed with glass beads for the collection of airborne viruses.

The liquid impingers can be used for sampling of viral aerosols, such as COVID-19 virus, influenza, and measles. However, the lowest cutoff diameter of commercially available liquid impingers was about 0.3 μm, and the physical collection efficiency for nano-bioaerosol is only about 10–20%. Here, we enhanced the impinger’s collection efficiency and recovery of viable viral aerosols by using packed glass beads and selected sampling media (1% peptone and lysogeny broth, LB). Single-stranded RNA (ssRNA) MS2 bacteriophage with uranine (as a physical tracer) was used as model viral aerosols. The effects of different sampling flow rates (4, 6, and 12.5 L per minute) and different sampling time (10, 20, and 30 min) on the collection efficiency and recovery of MS2 aerosols were also tested. Collection efficiency and recovery of viable viral aerosols were analyzed as a function of sampling media, flow rate, and sampling time and packed glass beads by using a general linear model. Although the packed glass beads considerably enhanced the collection efficiency of the liquid impinger for MS2 aerosols, the recovery of viable MS2 becomes lower due to the higher pressure drop across the impinger. Using peptone or LB as sampling media, reducing sampling flow rate, and decreasing sampling time was proven to improve the recovery of viable MS2. Conclusively, this study provides some practical methods to improve the collection efficiency of liquid impinger for viral aerosols and preserve their viability.

Packing simulations of complex-shaped rigid particles using FDEM: An application to catalyst pellets

A new component of the combined finite-discrete element method (FDEM) is employed to estimate the effects of geometrical features, friction and energy dissipation parameters on the bulk properties of rigid pellet packs. This work constitutes the first systematic validation of the Solidity FDEM code for rigid particles. The experimental and numerical axial and radial packing density profiles and orientation distributions have been compared, confirming that the numerical simulations of packing of cylindrical catalyst supports, glass beads and trilobe pellets deposited in a cylindrical container match the corresponding emergent bulk properties obtained from X-Ray CT scans. The presented results are a first confirmation of the applicability of FDEM based methods to the simulation of this class of multi-body problems. The assessment of the accuracy of the simulated topology of the pellet pack that is established in this work is encouraging for further investigations of the multi-physical engineering systems of interest for catalyst pellets involving hydro-thermo-mechanical, fracturing and fragmentation interactions using coupled FDEM formulations.

Experimental and numerical upscaling of foam flow in highly permeable porous media

Foam in porous media has been studied as a tool for various applications. Recently, the technology has become relevant for contaminated-aquifer remediation, where porous media are highly permeable. Therefore, the behavior of foam flow in high permeability porous media still raises numerous questions. In particular, upscaling of the foam flow from pore to Darcy scale is still under debate. Since the behavior of bulk foam has been studied principally in the food and cosmetics industries, and foam flow in porous media has mainly been investigated in the oil industry, the link between bulk-foam behavior and foam flow in porous media is still missing. The upscaling of foam flow from the pore scale to the laboratory scale could give valuable insight for understanding foam flow in aquifers.We studied the behavior of pre-generated foam with different foam qualities through the rheological characterization of bulk foam using a rheometer and also when flowing in a porous medium composed of 1 mm glass beads. Foam was formed by co-injecting surfactant solution and nitrogen gas through a porous column filled by fine sand. The homogenization method is used to study macroscopic foam flow properties in porous media by solving the non-linear boundary value problem. The rheology of bulk foam is then used as an input in the upscaling procedure for foam flow in different periodic model 2D and 3D unit cells.From our experiments, we found that the bulk foam is a yield-stress fluid and that the yield-stress values increase with foam quality. Moreover, the rheology of bulk foam corresponds well to the yield stress (Herschel-Bulkley-Papanastasiou) model. We found that foam behaves as a continuous yield-stress fluid in highly permeable porous media. It was also shown that the apparent foam viscosity in porous media increases with the foam quality at the same total flow rate. The results obtained from the rheometer successfully match the outcomes of apparent foam viscosity obtained by flow in porous media by a shifting parameter for the same foam quality. The apparent foam viscosity found in 1 mm glass-bead packing was much higher than bulk foam viscosity.Experimental results were compared to numerical results on simple unit cells. Although we observed considerable differences between the experimental and numerical results of upscaling, the general trend was identical. The differences can be explained by the complexity of the foam flow in porous media, especially foam compressibility. We found that foam flow at low capillary numbers is influenced by the trapping effect and at high pressure gradients by the compressibility. Compressibility was estimated for foam flow in 1 mm glass-bead packing. When foam compressibility is insignificant, the upscaling model can predict foam-flow behavior well at the Darcy scale.

Direct characterization of solute transport in unsaturated porous media using fast X-ray synchrotron microtomography

Solute transport in unsaturated porous materials is a complex process, which exhibits some distinct features differentiating it from transport under saturated conditions. These features emerge mostly due to the different transport time scales at different regions of the flow network, which can be classified into flowing and stagnant regions, predominantly controlled by advection and diffusion, respectively. Under unsaturated conditions, the solute breakthrough curves show early arrivals and very long tails, and this type of transport is usually referred to as non-Fickian. This study directly characterizes transport through an unsaturated porous medium in three spatial dimensions at the resolution of 3.25 μm and the time resolution of 6 s. Using advanced high-speed, high-spatial resolution, synchrotron-based X-ray computed microtomography (sCT) we obtained detailed information on solute transport through a glass bead packing at different saturations. A large experimental dataset (>50 TB) was produced, while imaging the evolution of the solute concentration with time at any given point within the field of view. We show that the fluids' topology has a critical signature on the non-Fickian transport, which yet needs to be included in the Darcy-scale solute transport models. The three-dimensional (3D) results show that the fully mixing assumption at the pore scale is not valid, and even after injection of several pore volumes the concentration field at the pore scale is not uniform. Additionally, results demonstrate that dispersivity is changing with saturation, being twofold larger at the saturation of 0.52 compared to that at the fully saturated domain.

Core-Shell Nanohydrogels with Programmable Swelling for Conformance Control in Porous Media.

Conformance control during waterflooding in an oil reservoir is utilized to redistribute water and increase the sweep efficiency and hence oil production. Using preformed gel particles can effectively redirect the flow by blocking the high-permeability zones and forcing water into low-permeability zones where the oil is trapped. However, the size of such gel particles can limit their applications deeper within the reservoir and can result in shear-induced degradation near the well bore. Here, we fabricate core-shell nanohydrogels with delayed swelling behavior; their volume increases by a factor of 200 after about 30 days in brine under reservoir conditions. We study their effect on the flow behavior in a three-dimensional porous medium micromodel consisting of randomly packed glass beads. Using confocal microscopy, we directly visualize the spatial variations of flow in the micromodel before and after nanohydrogel injection and swelling. The swollen nanohydrogels block some pores reducing the permeability of the micromodel and diverting the water into low-permeability regions. A core flood experiment further confirms that the nanohydrogels can significantly reduce the permeability of a reservoir sample and divert the fluid flow. Our results demonstrate that these core-shell nanohydrogels might be useful for flow control in porous media and can be used as a conformance control agent.

Event-based contact angle measurements inside porous media using time-resolved micro-computed tomography

HypothesisCapillary-dominated multiphase flow in porous materials is strongly affected by the pore walls’ wettability. Recent micro-computed tomography (mCT) studies found unexpectedly wide contact angle distributions measured on static fluid distributions inside the pores. We hypothesize that analysis on time-resolved mCT data of fluid invasion events may be more directly relevant to the fluid dynamics. ExperimentWe approximated receding contact angles locally in time and space on time-resolved mCT datasets of drainage in a glass bead pack and a limestone. Whenever a meniscus suddenly entered one or more pores, geometric and thermodynamically consistent contact angles in the surrounding pores were measured in the time step just prior to the displacement event. We introduced a new force-based contact angle, defined to recover the measured capillary pressure in the invaded pore throat prior to interface movement. FindingsUnlike the classical method, the new geometric and force-based contact angles followed plausible, narrower distributions and were mutually consistent. We were unable to obtain credible results with the thermodynamically consistent method, likely because of sensitivity to common imaging artifacts and neglecting dissipation. Time-resolved mCT analysis can yield a more appropriate wettability characterization for pore scale models, despite the need to further reduce image analysis uncertainties.

An automated testing device for continuous measurement of the hysteretic water retention curve of granular media

The water retention curve (WRC), representing an important key for the modelling of hydro-mechanical behaviour of unsaturated soils, is still not fully understood, because it originates from microscopic hydraulic and capillary phenomena. Furthermore, its experimental measurement, especially for cyclic drainage and imbibition paths, is challenging and time-consuming. In this contribution, a recently developed low-cost easy-to-use miniature testing device for the investigation of the WRC of unsaturated granular soils, such as coarse-grained sand and a packing of glass beads, is presented. With the help of the new device, that can be controlled by a Raspberry Pi single-board computer, the hysteretic WRC can be investigated in a conventional macroscopic approach by plotting the macroscopic specimen degree of saturation versus measured matric suction. The test set-up allows an automatic measurement of the WRC which is measured continuously following a programmed test procedure. In addition to the technical realisation of the new device, this contribution focuses on macroscopic results of water retention tests. Moreover, the testing device has been designed in a miniaturised size, in order to obtain microscopic insights into the phase distribution during cyclic drainage and imbibition paths with the help of computed tomography in future applications.

Applications of critical path analysis to uniform grain packings with narrow conductance distributions: I. Single-phase permeability

Accurate estimation of single-phase permeability (k) has broad application in numerous areas, particularly modeling flow and transport in porous materials. Various techniques have been proposed in the literature to estimate k from other medium's properties, such as porosity, grain and/or pore size distribution, and pore connectivity. Among them critical path analysis (CPA) from statistical physics − first developed to model fluid flow in media with broad conductance distributions − has been successfully applied to heterogeneous soils and rocks. However, its application to uniform sphere and/or glass bead packs that represent homogeneous porous media with narrow conductance distributions needs to be investigated. In this study, we invoke concepts from CPA and estimate k from average grain diameter and formation factor in uniform sand and glass bead packings. By comparing theory with eight datasets including 105 packs from the literature, we demonstrate that CPA estimates permeability in homogeneous media accurately. We also compare our CPA-based model estimations with those from the Kozeny–Carman, Revil and Cathles, and RGPZ models. Results indicate that the CPA approach estimates k more precisely than other three models studied here.

Slow dynamic nonlinearity in unconsolidated glass bead packs.

Slow dynamic nonlinearity describes a poorly understood, creeplike phenomena that occurs in brittle composite materials such as rocks and cement. It is characterized by a drop in stiffness induced by a mechanical conditioning, followed by a log(time) recovery. A consensus theoretical understanding of the behavior has not been developed. Here we introduce an alternative experimental venue with which to inform theory. Unconsolidated glass bead packs are studied rather than rocks or cement because the structure and internal contacts of bead packs are less complex and better understood. Slow dynamics has been observed in such systems previously. However, the measurements to date tend to be irregular. Particular care is used here in the experimental design to overcome the difficulties inherent in bead pack studies. This includes the design of the bead pack support, the use of low-frequency conditioning, and the use of ultrasonic waves as a probe with coda wave interferometry to assess changes. Slow dynamics is observed in our system after three different methods for low-frequency conditioning, one of which has not been reported in the literature previously.

Slow dynamics in a single glass bead.

Slow dynamic nonlinearity is ubiquitous amongst brittle materials, such as rocks and concrete, with cracked microstructures. A defining feature of the behavior is the logarithmic-in-time recovery of stiffness after a mechanical conditioning. Materials observed to exhibit slow dynamics are sufficiently different in microstructure, chemical composition, and scale (ranging from the laboratory to the seismological) to suggest some kind of universality. A consensus of theoretical understanding of the universality in general and the log(time) recovery in particular is lacking. Seminal studies were focused on sandstones and other natural rocks, but in recent years other experimental venues have been introduced with which to inform theory. One such system is unconsolidated glass bead packs. However, bead packs still contain many contact points. The force distribution amongst the contacts is unknown. Here, we present slow dynamics measurements on a yet simpler system-a single glass bead confined between two large glass plates. The system is designed with a view towards rapid control of the contact zone environment. Ultrasonic waves are used as a probe of the system, and changes are assessed with coda wave interferometry. Three different methods of low-frequency conditioning are applied; all lead to slow dynamic recoveries. Results imply that force chains do not play an essential role in granular media slow dynamics, as they are absent in our system.

An Experimental Study of Emulsion Flow in Alkaline/Solvent Coinjection with Steam for Heavy-Oil/Bitumen Recovery

SummaryWater is the dominant component in steam-injection processes, such as steam-assisted gravity drainage (SAGD). The central hypothesis in this research is that in-situ oil transport can be enhanced by generating oil-in-water emulsion, where the water-continuous phase acts as an effective oil carrier. As part of the research project, this paper presents an experimental study of how oil-in-water emulsion can improve oil transport in porous media at elevated temperatures.Diethylamine (DEA) was selected as the organic alkali that generates oil-in-water emulsions with Athabasca bitumen at a 1,000-ppm NaCl brine and a 0.5-wt% alkali concentration. This aqueous composition had been confirmed to be an optimum in terms of oil content in the water-external emulsion phase at a wide range of temperatures. Then, flow experiments with a glass-bead pack were conducted to measure the effective viscosities of emulsion samples at shear rates from 5 to 29 seconds−1 at 35 bar and temperatures from 373 to 443 K.Results show that the oil-in-water emulsions were more than 15 times less viscous than the original bitumen at temperatures from 373 to 443 K. At the shear rate of 5 seconds−1, for example, the emulsion viscosity was 12 cp at 373 K, at which the bitumen viscosity was 206 cp. The efficiency of in-situ bitumen transport was evaluated by calculating the bitumen molar flow rate under gravity drainage with the new experimental data. Results show that oil-in-water emulsion can enhance the in-situ molar flow of bitumen by a factor of 273 at 403 K and 345 at 373 K, in comparison with the two-phase flow of oil and water in conventional SAGD. At 443 K, only a fraction of bitumen is emulsified in water, but the bitumen transport by both oil-in-water emulsion and an excess oil phase in DEA-SAGD can enhance the molar flow of bitumen by a factor of 19 in comparison to SAGD. This is mainly because the mobility of the bitumen-containing phase is enhanced by the reduced viscosity and increased effective permeability. A marked difference between alkaline solvents and conventional hydrocarbon solvents is that only a small amount of an alkaline solvent enables enhancing the in-situ transport of bitumen.

Plasma-assisted abatement of SF6 in a dielectric barrier discharge reactor: investigation of the effect of packing materials

The greenhouse gas hexafluoride (SF6) and its degradation is a research hotspot in the environmental field. In a double-layer dielectric barrier discharge reactor, the degradation process of SF6 was studied when packed in glass beads or γ-Al2O3 pellets as well as with no packing system. The results show that the packing materials changed the physical characteristics of the discharge, increasing the equivalent capacitance and discharge power and leading to the generation of more active particles and high-energy electrons. Additionally, the two packing materials effectively improved the degradation removal efficiency (DRE) and energy yield (EY) of SF6. In the γ-Al2O3 system, 2% SF6 gas obtained a maximum DRE of 85.97% with a 110 W input power and the EY reached 9.17 g kWh−1, far exceeding the case without packing. Moreover, the glass bead packing had no obvious effect on the degradation products, while the packing of γ-Al2O3 particles promoted SO2 production and reduced the formation of SOF2, SO2F2 and SOF4. Meanwhile, in the γ-Al2O3 system, the yield of SO2 increased with increasing input power and some undesirable products such as SO2F2 were reduced. These results show that the addition of suitable packing materials such as γ-Al2O3 pellets can effectively promote the harmless degradation of SF6.

Analysis of the Kozeny–Carman model based on pore networks

Abstract The Kozeny–Carman (KC) model is commonly used to predict permeability (k) but sometimes presents obvious predictive deviations. Comparatively, k prediction based on pore networks could more effectively describe the dependence of k on porous structure and accordingly shows high accuracy and reliability. This triggers us to examine the rationality of the KC assumptions about the pore structure and analyze the model applicability according to the pore networks. Two glass bead packs, a sintered glass bead plate and a Berea core are measured using computed tomography imaging and their pore parameters are accordingly calculated. Their pore parameters are obviously distinct, generally reflecting the progressively stronger roles of particle size, compaction-alike sintering and weak cementation in reducing porosity (ϕ), k, pore and pore throat radii (rp and rt). When correlating the pore parameters of the KC model to those of the pore networks, it is found that the specific surface area (av) has no specific relation to rp and rt and that pore structures could be distinct despite the equivalent av. Thus, av is insufficient to distinguish the differences in pore geometry and reflects their influence on k. An analysis of the modified KC equations of our cores and the ϕ−k relationships of some relative homogeneous reservoir rocks (e.g. silty sandstone and Fontainebleau sands) indicates that the variety of the ϕ−k relationships induced by such factors as particle size and pore radius could not be fully predicted by the KC model in some cases, especially when the ϕ is relatively low.

Application of Non-toxic Yield Stress Fluids Porosimetry Method and Pore-Network Modelling to Characterize the Pore Size Distribution of Packs of Spherical Beads

With X-ray computed tomography still being flawed as a result of limitations in terms of spatial resolution and cost, toxic mercury intrusion porosimetry (MIP) is nowadays the prevailing technique to determine PSDs of most porous media. Recently, yield stress fluids porosimetry method (YSM) has been identified as a promising clean alternative to MIP. This technique is based on the particular percolation patterns followed by yield stress fluids, which only flow through certain pores when injected at a given pressure gradient. In previous works, YSM was used to characterize natural and synthetic porous media, and the results were compared with MIP showing reasonable agreement. However, considerable uncertainty still remains regarding the characterized pore dimension with each method arising from the highly complex geometry of the interstices in real porous media. Therefore, a critical stage for the validation of YSM consists in achieving successful characterization of model porous media with well-known pore morphology and topology. With this objective in mind, a set of four packs of glass beads each with a given monodisperse bead size were characterized in the present work using different porosimetry methods: experimental YSM, numerically simulation of MIP and pore-network extraction from a 3D image. The results provided by these techniques were compared, allowing the identification of the pore dimensions being characterized in each case. The results of this research elucidate the causes of the discrepancies between the considered porosimetry methods and demonstrate the usefulness of the PSD provided by YSM when predicting flow in porous media.

A two-region transport model for interpreting T1-T2 measurements in complex systems.

A 1D two region coupled pore model with discrete pore coupling is developed to elucidate the eigenmode interactions in regions with different surface relaxivity. Numerical solution of the model and simulation of the correlation experiment for varying surface relaxivity, pore connectivity and pore size ratio indicate the role of negative eigenmodes and overlap of T1 and T2 eigenmodes in generating a time domain signal increase with inversion recovery time, t1. The eigenmodes and eigenfunctions are considered in detail providing connection between the mathematical model and the diffusion dynamics and spin physics of the system. Physical systems, i.e. a microporous glass bead pack, a cyclopentane/water hydrate former, and beeswax, showing experimentally measured T1-T2 time domain signal rise are considered within the limitations of the model.

Acoustic emission generated by glass beads in compression and shearing

Acoustic emission (AE) monitoring offers the potential to sense particle-scale interactions that lead to macro-scale responses of granular materials. This paper presents results from a programme of drained triaxial tests performed on densely packed glass beads to establish quantitative interpretation of AE during isotropic compression, shearing and associated stick–slip events. Relationships have been quantified between: AE and boundary work (i.e. AE generated per Joule) for a unit volume of glass beads under isotropic compression and shear; AE and shear displacement rate; and the amplitude of deviator stress cycles and AE activity during stick–slip events. In shear, AE generation increased with shear strain and reached peak values that were maintained from volume minimum (i.e. the transition from contractive to dilative behaviour) to peak dilatancy, whereupon AE generation gradually reduced and then remained around a constant mean value with further increments of shear strain. In each stick–slip event, AE activity increased during shear strength mobilisation, particle climbing and dilation, and then reduced with the subsequent deviator stress drop during particle sliding and contraction. The amplitude of these cycles in AE activity were governed by the amplitude of deviator stress cycles during stick–slip events, which were also proportional to the imposed stress level and inversely proportional to particle size.

Pore Scale Visualization of Drainage in 3D Porous Media by Confocal Microscopy

We visualize the dynamics of immiscible displacement of a high viscosity wetting phase by a low viscosity non-wetting phase in a three-dimensional (3D) glass bead packing using confocal microscopy. Both phases were doped with two different fluorescent dyes, which enabled visualization of both phases simultaneously and quantification of the phase volumes without the need of image subtraction operations. The transient results show details of the displacement process and how pores are invaded by the non-wetting displacing phase. The static images at the end of the displacement process reveal how the trapped ganglia volume and morphology change with capillary number. The wetting phase is trapped as pendular rings spanning one to multiple pore necks. Details of the pore scale flow of oil wet media revealed with the experimental methods presented here can lead to better fundamental understanding of the physical processes and optimized enhanced oil recovery methods, CO2 sequestration and aquifer remediation.

Experimental study of the dynamic mechanism on gas bubbles migration, fragment, coalescence and trapping in a porous media

Immiscible gas injection recovery is one of the most important methods used for EOR. However, it is very difficult to directly observe the flow of gas bubbles in the reservoir stratum. In our experimental work, we investigated visually the morphological configurations of gas bubbles migration, fragment, coalescence and trapping in different liquid-simulated porous media packs of glass beads. Also, the volume fraction of the trapped gas bubbles in the drainage process were measured. A fast video recording technique was used for the visual observations of the morphological characteristics of the rising gas bubbles in the saturated porous media. The obtained experimental images recorded reveal three types of the rising migration bubbles namely: free migration in the pore throat, attached wall migration on the beads and deformation migration. Additionally, coalescence and fragment have an interdependency for a rising gas bubble. The adsorption and capillary trappings are both two kinds of important trapping mechanisms, which determine the volume fraction of the trapped bubbles. The volume fraction of the trapped bubbles increases with an increment in the gas injection time and rates with the trapped bubbles fluctuating upward and downward due to the more coalescence and fragment after a long gas injection time/period. Lastly, the volume fraction of the trapped bubbles increases with a decrease in the inclination angles and porosities, which influences the areal sweep efficiency for enhanced liquid recovery in a porous media.

Ultrasonic tracking of a sinking ball in a vibrated dense granular suspension

Observing and understanding the motion of an intruder through opaque dense suspensions such as quicksand remains a practical and conceptual challenge. Here we use an ultrasonic probe to monitor the sinking dynamics of a steel ball in a dense glass bead packing (3D) saturated by water. We show that the frictional model developed for dry granular media can be used to describe the ball motion induced by horizontal vibration. From this rheology, we infer the static friction coefficient and effective viscosity that decrease when increasing the vibration intensity. Our main finding is that the vibration-induced reduction of the yield stress and increase of the sinking depth are presumably due to micro-slips induced at the grain contacts but without visible plastic deformation due to macroscopic rearrangements, in contrast to dry granular packings. To explain these results, we propose a mechanism of acoustic lubrication that reduces the inter-particle friction and leads to a decrease of the yield stress. This scenario is different from the mechanism of liquefaction usually invoked in loosely packed quicksands where the vibration-induced compaction increases the pore pressure and decreases the confining pressure on the solid skeleton, thus reducing the granular resistance to external load.