Dynamic Pore-Scale Model of Drainage in Granular Porous Media: The Pore-Unit Assembly Method.

Dynamics of drainage is analyzed for packings of spheres, using numerical experiments. For this purpose, a dynamic pore-scale model was developed to simulate water flow during drainage. The pore space inside a packing of spheres was extracted using regular triangulation, resulting in an assembly of grain-based tetrahedra. Then, pore units were constructed by identifying and merging tetrahedra that belong to the same pore, resulting in an assembly of pore units. Each pore unit was approximated by a volume-equivalent regular shape (e.g., cube and octahedron), for which a local capillary pressure-saturation relationship was obtained. To simulate unsaturated flow, a pore-scale version of IMPES (implicit pressure solver and explicit saturation update) was employed in order to calculate pressure and saturation distributions as a function of time for the assembly of pore units. To test the dynamic model, it was used on a packing of spheres to reproduce the corresponding measured quasi-static capillary pressure-saturation curve for a sand packing. Calculations were done for a packing of spheres with the same grain size distribution and porosity as the sand. We obtained good agreement, which confirmed the ability of the dynamic code to accurately describe drainage under low flow rates. Simulations of dynamic drainage revealed that drainage occurred in the form of finger-like infiltration of air into the pore space, caused by heterogeneities in the pore structure. During the finger-like infiltration, the pressure difference between air and water was found to be significantly higher than the capillary pressure. Furthermore, we tested the effects of the averaging, boundary conditions, domain size, and viscosity on the dynamic flow behavior. Finally, the dynamic coefficient was determined and compared to experimental data.

Modelling the Evolution of Pore Structure during the Disintegration of Pharmaceutical Tablets.

Pharmaceutical tablet disintegration is a critical process for dissolving and enabling the absorption of the drug substance into the blood stream. The tablet disintegration process consists of multiple connected and interdependent mechanisms: liquid penetration, swelling, dissolution, and break-up. One key dependence is the dynamic change of the pore space in a tablet caused by the swelling of particles while the tablet takes up liquid. This study analysed the changes in the pore structure during disintegration by coupling the discrete element method (DEM) with a single-particle swelling model and experimental liquid penetration data from terahertz-pulsed imaging (TPI). The coupled model is demonstrated and validated for pure microcrystalline cellulose (MCC) tablets across three porosities (10, 15, and 22%) and MCC with three different concentrations of croscarmellose sodium (CCS) (2, 5, and 8% w/w). The model was validated using experimental tablet swelling from TPI. The model captured the difference in the swelling behaviour of tablets with different porosities and formulations well. Both the experimental and modelling results showed that the swelling was lowest (i.e., time to reach the maximum normalised swelling capacity) for tablets with the highest CCS concentration, cCCS = 8%. The simulations revealed that this was caused by the closure of the pores in both the wetted volume and dry volume of the tablet. The closure of the pores hinders the liquid from accessing other particles and slows down the overall swelling process. This study provides new insights into the changes in the pore space during disintegration, which is crucial to better understand the impact of porosity and formulations on the performance of tablets.

The impact of water saturation on the infiltration behaviour of elemental mercury DNAPL in heterogeneous porous media.

Industrial use has led to the presence of liquid elemental mercury (Hg0) worldwide in the subsurface as Dense NonAqueous Phase Liquid (DNAPL), resulting in long lasting sources of contamination, which cause harmful effects on human health and detrimental consequences on ecosystems. However, to date, insight into the infiltration behaviour of elemental mercury DNAPL is lacking. In this study, a two-stage flow container experiment of elemental mercury DNAPL infiltration into a variably water saturated stratified sand is described. During the first stage of the experiment, 16.3 ml of liquid Hg0 infiltrated and distributed into an initially partially water saturated system. Afterwards, during the second stage of the experiment, consisting of the simulation of a "rain event" to assess whether the elemental mercury already infiltrated could be mobilized due to local increases in water saturation, a significant additional infiltration of 4.7 ml of liquid mercury occurred from the remaining DNAPL source. The experiment showed that, under conditions similar to those found in the field, Hg0 DNAPL infiltration is likely to occur via fingers and is strongly controlled by porous medium structure and water saturation. Heterogeneities within the porous medium likely determined preferential pathways for liquid Hg0 infiltration and distribution, as also suggested by dual gamma ray measurements. Overall, this study highlights that the infiltration behaviour of mercury DNAPL is strongly impacted by water saturation. In the field, this may result in a preferential infiltration of Hg0 DNAPL in wetter areas or in its mobilization due to wetting during a rain event, as indicated by this study, or a groundwater table rise. This should be considered when assessing the likely distribution pathways of historic mercury DNAPL contamination as well as the remediation efforts.

Infiltration behaviour of elemental mercury DNAPL in fully and partially water saturated porous media.

Mercury is a contaminant of global concern due to its harmful effects on human health and for the detrimental consequences of its release in the environment. Sources of liquid elemental mercury are usually anthropogenic, such as chlor-alkali plants. To date insight into the infiltration behaviour of liquid elemental mercury in the subsurface is lacking, although this is critical for assessing both characterization and remediation approaches for mercury DNAPL contaminated sites. Therefore, in this study the infiltration behaviour of elemental mercury in fully and partially water saturated systems was investigated using column experiments. The properties affecting the constitutive relations governing the infiltration behaviour of liquid Hg0, and PCE for comparison, were determined using Pc(S) experiments with different granular porous media (glass beads and sands) for different two- and three-phase configurations. Results showed that, in water saturated porous media, elemental mercury, as PCE, acted as a non-wetting fluid. The required entry head for elemental mercury was higher (from about 5 to 7 times). However, due to the almost tenfold higher density of mercury, the required NAPL entry heads of 6.19cm and 12.51cm for mercury to infiltrate were 37.5% to 20.7% lower than for PCE for the same porous media. Although Leverett scaling was able to reproduce the natural tendency of Hg0 to be more prone than PCE to infiltrate in water saturated porous media, it considerably underestimated Hg0 infiltration capacity in comparison with the experimental results. In the partially water saturated system, in contrast with PCE, elemental mercury also acted as a nonwetting fluid, therefore having to overcome an entry head to infiltrate. The required Hg0 entry heads (10.45 and 15.74cm) were considerably higher (68.9% and 25.8%) than for the water saturated porous systems. Furthermore, in the partially water saturated systems, experiments showed that elemental mercury displaced both air and water, depending on the initial water distribution within the pores. This indicates that the conventional wettability hierarchy, in which the NAPL has an intermediate wetting state between the air and the water phases, is not valid for liquid elemental mercury. Therefore, for future modelling of elemental mercury DNAPL infiltration behaviour in variably water saturated porous media, a different formulation of the governing constitutive relations will be required.

The Effects of Swelling and Porosity Change on Capillarity: DEM Coupled with a Pore-Unit Assembly Method.

In this study, a grain-scale modelling technique has been developed to generate the capillary pressure-saturation curves for swelling granular materials. This model employs only basic granular properties such as particles size distribution, porosity, and the amount of absorbed water for swelling materials. Using this model, both drainage and imbibition curves are directly obtained by pore-scale simulations of fluid invasion. This allows us to produce capillary pressure-saturation curves for a large number of different packings of granular materials with varying porosity and/or amount of absorbed water. The algorithm is based on combining the Discrete Element Method for generating different particle packings with a pore-unit assembly approach. The pore space is extracted using a regular triangulation, with the centres of four neighbouring particles forming a tetrahedron. The pore space within each tetrahedron is referred to as a pore unit. Thus, the pore space of a particle packing is represented by an assembly of pore units for which we construct drainage and imbibition capillary pressure-saturation curves. A case study on Hostun sand is conducted to test the model against experimental data from literature and to investigate the required minimum number of particles to have a Representative Elementary Volume. Then, the capillary pressure-saturation curves are constructed for Absorbent Gelling Material particles, for different combinations of porosity values and amounts of absorbed water. Each combination yields a different configuration of pore units, and thus distinctly different capillary pressure-saturation curves. All these curves are shown to collapse into one curve for drainage and one curve for imbibition when we normalize capillary pressure and saturation values. We have developed a formula for the Van Genuchten parameter [Formula: see text] (which is related to the inverse of the entry pressure) as a function of porosity and the amount of absorbed water.

The transport behaviour of elemental mercury DNAPL in saturated porous media: analysis of field observations and two-phase flow modelling.

Mercury is a contaminant of global concern. The use of elemental mercury in various (former) industrial processes, such as chlorine production at chlor-alkali plants, is known to have resulted in soil and groundwater contaminations worldwide. However, the subsurface transport behaviour of elemental mercury as an immiscible dense non-aqueous phase liquid (DNAPL) in porous media has received minimal attention to date. Even though, such insight would aid in the remediation effort of mercury contaminated sites. Therefore, in this study a detailed field characterization of elemental mercury DNAPL distribution with depth was performed together with two-phase flow modelling, using STOMP. This is to evaluate the dynamics of mercury DNAPL migration and the controls on its distribution in saturated porous media. Using a CPT-probe mounted with a digital camera, in-situ mercury DNAPL depth distribution was obtained at a former chlor-alkali-plant, down to 9 m below ground surface. Images revealing the presence of silvery mercury DNAPL droplets were used to quantify its distribution, characteristics and saturation, using an image analysis method. These field-observations with depth were compared with results from a one-dimensional two-phase flow model simulation for the same transect. Considering the limitations of this approach, simulations reasonably reflected the variability and range of the mercury DNAPL distribution. To further explore the impact of mercury's physical properties in comparison with more common DNAPLs, the migration of mercury and PCE DNAPL in several typical hydrological scenarios was simulated. Comparison of the simulations suggest that mercury's higher density is the overall controlling factor in controlling its penetration in saturated porous media, despite its higher resistance to flow due to its higher viscosity. Based on these results the hazard of spilled mercury DNAPL to cause deep contamination of groundwater systems seems larger than for any other DNAPL.