Pore-scale simulation of fluid flow and solute dispersion in three-dimensional porous media

Impact of hydraulic tortuosity on microporous and nanoporous media flow

Using two-dimensional porous structures made up of homogeneously arranged solid obstacles, we examine the effects of rarefaction on the hydraulic tortuosity in the slip and early transition flow regimes via extended lattice Boltzmann method. We observed that modification in either the obstacle's arrangement or the porosity led to a power-law relation between the porosity-tortuosity. Along with this, we also found that in the slip-flow regime, the exponent of this relation contains the effect of finite Knudsen number (Kn). In addition, we observed that on properly scaling Kn with porosity and hydraulic tortuosity, a generalized correlation can be obtained for apparent permeability.

Modeling flow of Carreau fluids in porous media

Carreau fluids occur routinely in porous medium systems for a range of applications, and the dependence of the viscosity for such fluids on the rate of strain tensor poses challenges to modeling at an averaged macroscale. Traditional approaches for macroscale modeling such flows have relied upon experimental observations of flows for generalized Newtonian fluids (GNFs) and a phenomenological approach referred to herein as the shift factor. A recently developed approach based upon averaging conservation and thermodynamic equations from the microscale for Cross model GNFs is extended to the case of Carreau fluids and shown to predict the flow through both isotropic and anisotropic media accurately without the need for GNF-flow experiments. The model is formulated in terms of rheological properties, a standard Newtonian resistance tensor, and a length-scale tensor, which does require estimation. An approach based upon measures of the morphology and topology of the pore space is developed to approximate this length-scale tensor. Thus, this work provides the missing components needed to predict Carreau GNF macroscale flow with only rheological information for the fluid and analysis of the pore morphology and topology independent of any fluid flow experiments. Accuracy of predictions based upon this approach is quantified, and extension to other GNFs is straightforward.

Analysis of Particle-Resolved CFD Results for Dispersion in Packed Beds

Dispersion is the spreading of a solute while it is moved by a flowing medium. The study of dispersion in catalytic chemical reactors is fundamental to their design, since dispersion influences the reactant and product transport within the bed. In this paper, longitudinal and transverse dispersion of an inert tracer in slender packed beds of spheres and spherocylinders is studied using Computational Fluid Dynamics simulations. The focus is on the analysis of dispersion from full field data. The purpose is to develop a methodology that can later also be used to characterize dispersion from full field experimental data such as MRI measurements. Results obtained by means of particle-resolved CFD simulations are discussed. Spatial distributions and residence times are analyzed and the results are interpreted by comparison to results obtained from 1D and 2D convection-diffusion equations.

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

Microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. A two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

Generalized Newtonian fluid flow in porous media

Single-fluid-phase porous medium systems are typically modeled at an averaged length scale termed the macroscale, and Darcy's law is typically relied upon as an approximation of the momentum equation under Stokes flow conditions. Standard approaches for modeling macroscale single-fluid-phase flow of generalized Newtonian fluids (GNFs) extend the standard Newtonian model based upon Darcy's law using an effective viscosity and assuming that the intrinsic permeability is invariant with respect to fluid properties. This approach results in a need to perform an experiment for a non-Newtonian fluid, the introduction of effective parameters that are not tied to known microscale physics, and uncertainty regarding the dependencies of the fitted empirical parameter on system properties. We use the thermodynamically constrained averaging theory (TCAT) to examine the formulation and closure of a macroscale model for GNF flow that is consistent with microscale conservation principles and the second law of thermodynamics. A direct connection between microscale and macroscale quantities is used to formulate an expression for interphase momentum transfer for GNF flow in porous medium systems. Darcy's law is shown to approximate momentum transfer from the fluid phase to the solid phase. Momentum transfer is found to depend on the viscosity at the solid surface, which is only invariant for Newtonian flow. TCAT is used to derive a macroscale equation for the hydraulic resistance based on accessible fluid and solid properties. This hydraulic resistance may be used in the same way that hydraulic conductivity is typically used to model flow at the macroscale, and it includes parameters that can be calculated a priori, without the need to carry out microscale simulations, or experiments, for any GNF. The TCAT approach is validated for four model isotropic and anisotropic media and five Cross-model fluids. The traditional shift factor and effective viscosity are related to the newly derived TCAT model, shedding new light on this common empirical approach. The results from this work form a basis for the modeling of GNF flow in porous medium systems under Stokes flow, which is predictive given the rheological properties of the GNF and the resistance observed for Newtonian flow.

Distribution of liquid flow in a pore network during evaporation

The variation of the distribution of the liquid flow in porous media during evaporation is still a puzzle. We resolve it with the pore network modeling approach. The distribution of the evaporation-induced liquid flow in a pore network composed of about 2.5 million pores is determined. The probability density function of the magnitude of the normalized liquid flow rate is obtained. For the low normalized liquid flow rate, the probability density function is power-lawlike. The power-law exponent depends on both the liquid saturation and the location of the moving meniscus in the main liquid cluster. The evaporation-induced liquid flow in the pores in the pore network can be correlated. Whether the liquid flow distributions in various zones in the pore network are similar or not relies significantly on the location of the moving meniscus in the main liquid cluster. The functions depicting the relation between the power-law exponent and the local liquid saturation for the zones adjacent to and away from the open side of the pore network are different. These findings from the pore scale studies provide insights into developing the accurate continuum model for evaporation in porous media.

Evaporation-Induced Diffusion Acceleration in Liquid-Filled Porous Materials

Liquid-filled porous materials exist widely in nature and engineering fields, with the diffusion of substances in them playing an important role in system functions. Although surface evaporation is often inevitable in practical scenarios, the evaporation effects on diffusion behavior in liquid-filled porous materials have not been well explored yet. In this work, we performed noninvasive diffusion imaging experiments to observe the diffusion process of erioglaucine disodium salt dye in a liquid-filled nitrocellulose membrane under a wide range of relative humidities (RHs). We found that evaporation can significantly accelerate the diffusion rate and alter concentration distribution compared with the case without evaporation. We explained the accelerated diffusion phenomenon by the mechanism that evaporation would induce a weak flow in liquid-filled porous materials, which leads to convective diffusion, i.e., evaporation-induced flow and diffusion (EIFD). Based on the EIFD mechanism, we proposed a convective diffusion model to quantitatively predict the diffusion process in liquid-filled porous materials under evaporation and experimentally validated the model. Introducing the dimensionless Peclet (P _e) number to measure the relative contribution of the evaporation effect to pure molecular diffusion, we demonstrated that even at a high RH of 95%, where the evaporation effect is usually assumed negligible in common sense, the evaporation-induced diffusion still overwhelms the molecular diffusion. The flow velocity induced by evaporation in liquid-filled porous materials was found to be 0.4-5 μm/s, comparable to flow in many biological and biomedical systems. The present analysis may help to explain the driving mechanism of tissue perfusion and provide quantitative analysis or inspire new control methods of flow and material exchange in numerous cutting-edge technologies, such as paper-based diagnostics, hydrogel-based flexible electronics, evaporation-induced electricity generation, and seawater purification.

Reconstruction and seepage simulation of a coal pore-fracture network based on CT technology

The distribution of multiscale pores and fractures in coal and rock is an important basis for reflecting the capacity of fluid flow in coal seam seepage passages. Accurate extraction and qualitative and quantitative analysis of pore-fracture structures are helpful in revealing the flow characteristics of fluid in seepage channels. The relationship between pore and fracture connectivity can provide a scientific reference for optimizing coal seam water injection parameters. Therefore, to analyse the change in permeability caused by the variability in the coal pore-fracture network structure, a CT scanning technique was used to scan coal samples from the Leijia District, Fuxin. A total of 720 sets of original images were collected, a median filter was used to filter out the noise in the obtained images, and to form the basis of a model, the reconstruction and analysis of the three-dimensional pore-fracture morphology of coal samples were carried out. A pore-fracture network model of the coal body was extracted at different scales. Using the maximum sphere algorithm combined with the coordination number, the effect of different quantitative relationships between pore size and pore throat channel permeability was studied. Avizo software was used to simulate the flow path of fluid in the seepage channels. The change trend of the fluid velocity between different seepage channels was discussed. The results of the pore-fracture network models at different scales show that the pore-fracture structure is nonuniform and vertically connected, and the pores are connected at connecting points. The pore size distribution ranges from 104 μm to 9425 μm. The pore throat channel length distribution ranges from 4206 μm to 48073 μm. The size of the coordination number determines the connectivity and thus the porosity of the coal seam. The more connected pore channels there are, the larger the pore diameters and the stronger the percolation ability. During flow in the seepage channels of the coal, the velocity range is divided into a low-speed region, medium-speed region and high-speed region. The fluid seepage in the coal seam is driven by the following factors: pore connectivity > pore and pore throat dimensions > pore and pore throat structure distribution. Ultimately, the pore radius and pore connectivity directly affect the permeability of the coal seam.

Pore-Scale Mixing and the Evolution of Hydrodynamic Dispersion in Porous Media

We study the interplay of pore-scale mixing and network-scale advection through heterogeneous porous media, and its role for the evolution and asymptotic behavior of hydrodynamic dispersion. In a Lagrangian framework, we identify three fundamental mechanisms of pore-scale mixing that determine large scale particle motion, namely, the smoothing of intrapore velocity contrasts, the increase of the tortuosity of particle paths, and the setting of a maximum time for particle transitions. Based on these mechanisms, we derive a theory that predicts anomalous and normal hydrodynamic dispersion in terms of the characteristic pore length, Eulerian velocity distribution, and Péclet number.

The Complexity of Porous Media Flow Characterized in a Microfluidic Model Based on Confocal Laser Scanning Microscopy and Micro-PIV

Abstract

In this study, the complexity of a steady-state flow through porous media is revealed using confocal laser scanning microscopy (CLSM). Micro-particle image velocimetry (micro-PIV) is applied to construct movies of colloidal particles. The calculated velocity vector fields from images are further utilized to obtain laminar flow streamlines. Fluid flow through a single straight channel is used to confirm that quantitative CLSM measurements can be conducted. Next, the coupling between the flow in a channel and the movement within an intersecting dead-end region is studied. Quantitative CLSM measurements confirm the numerically determined coupling parameter from earlier work of the authors. The fluid flow complexity is demonstrated using a porous medium consisting of a regular grid of pores in contact with a flowing fluid channel. The porous media structure was further used as the simulation domain for numerical modeling. Both the simulation, based on solving Stokes equations, and the experimental data show presence of non-trivial streamline trajectories across the pore structures. In view of the results, we argue that the hydrodynamic mixing is a combination of non-trivial streamline routing and Brownian motion by pore-scale diffusion. The results provide insight into challenges in upscaling hydrodynamic dispersion from pore scale to representative elementary volume (REV) scale. Furthermore, the successful quantitative validation of CLSM-based data from a microfluidic model fed by an electrical syringe pump provided a valuable benchmark for qualitative validation of computer simulation results.

Graphic Abstract

Microscale modeling of nondilute flow and transport in porous medium systems

Nondilute transport occurs routinely in porous medium systems. Experimental observations have revealed effects that seemingly depend upon density, viscosity, velocity, and chemical activity. Macroscale models based upon averaged behavior over many pores have been relied upon to describe such systems to date, which require parametrization of important physical phenomena in material coefficients. To advance fundamental understanding of these complex systems, we examine nondilute transport from a fundamental microscale, or pore-scale, continuum modeling perspective. We approximate the solution of a model based upon the variable-density Navier-Stokes equations and a nondilute species transport equation. Known dependencies of the densities, viscosities, chemical activity, and diffusion for a salt solution on chemical composition are included in the model. Microscale model solutions are averaged to the macroscale and compared with extant experimental observations. Investigation of the effects of various physical phenomena on the microscale velocity distribution and the observed macroscale dispersion are considered using dimensional analysis and constrained simulations. Simulation results are used to explain observed experimental results in light of underlying mechanisms. Conditions under which the various physicochemical effects investigated are important are revealed.

Plume deformation, mixing, and reaction kinetics: An analysis of interacting helical flows in three-dimensional porous media

Heterogeneity and macroscopic anisotropy of porous media play an important role for dilution and reaction enhancement of conservative and reactive plumes. In this study, we perform numerical simulations to investigate steady-state flow and transport in three-dimensional heterogeneous porous media. We consider two macroscopic anisotropic inclusions resulting in helical flows with twisting streamlines in a three-dimensional flow-through domain. The inclusions are obtained by alternating two layers of angled slices of coarse and fine porous media with different hydraulic conductivity. We investigate flow and transport scenarios considering different geometry and relative position of the two anisotropic inclusions yielding helical flow fields with different extent of interaction. We use metrics of stretching and folding to characterize the flow field and entropy-based metrics for the analysis of the conservative and reactive transport problems. The outcomes show that the two helices result in different patterns of twisting streamlines, which cause distinct deformation of the plumes. However, mixing and reaction enhancement could not be directly related to the extent of the flow field deformation: Configurations with strong deformation can result in only moderate mixing enhancement, whereas configurations with limited deformation of the flow field can lead to significant mixing of the solute plume. Finally, we explore the impact of different degradation rates on reactive transport and the role of reaction kinetics on the entropy balance for a reactant undergoing transport and mixing-controlled degradation in the twisting flow fields. The results show that strong mixing enhancement due to helical flow increases the importance of the reaction kinetics that becomes the rate-limiting process for solute reactive transport.

Solute transport in porous media studied by lattice Boltzmann simulations at pore scale and x-ray tomography experiments

With the aid of nondestructive microfocus x-ray computed tomography (CT), we performed three-dimensional (3D) tracer dispersion experiments on randomly unconsolidated packed beds. Plumes of nonreactive sodium iodide solution were point injected into a sodium chloride solvent as a tracer for the evaluation of the dispersion process. The asymptotic dispersion coefficient was obtainable within the experimental scale and was summarized over Péclet numbers from 11.7 to ∼860. Then, the lattice Boltzmann method and moment propagation method were used to elucidate the mechanisms embedded in the dispersion phenomenon. The methods were rigorously verified against the classical Taylor dispersion problem and extended to simulate fluid flow and tracer dispersion in high-resolution 3D digital porous structures from CT. The method of moments, Lagrangian velocity correction function, and dilution index were thoroughly analyzed to evaluate the dispersion behaviors. Numerical simulations revealed ballistic and superdiffusive regimes at the transient times, whereas asymptotic dispersion behaviors appear at longer characteristic times. Besides, the observed transient times unanimously persist over convective length scales of around 12 particles transversely and 16 particles longitudinally. The estimated dispersion coefficients from simulation are in consistence with the experimental result. Furthermore, the simulation also enabled the identification of regimes, including diffusive, power law, and mechanical dispersion. Thus, the proposed experimental and computational schemes are of practical means to study dispersion behaviors by direct pore scale imaging and modeling.

Ubiquity of anomalous transport in porous media: Numerical evidence, continuous time random walk modelling, and hydrodynamic interpretation

Anomalous transport in porous media is commonly believed to be induced by the highly complex pore space geometry. However, this phenomenon is also observed in porous media with rather simple pore structure. In order to answer how ubiquitous can anomalous transport be in porous media, we in this work systematically investigate the solute transport process in a simple porous medium model with minimal structural randomness. The porosities we consider range widely from 0.30 up to 0.85, and we find by lattice Boltzmann simulations that the solute transport process can be anomalous in all cases at high Péclet numbers. We use the continuous time random walk theory to quantitatively explain the observed scaling relations of the process. A plausible hydrodynamic origin of anomalous transport in simple porous media is proposed as a complement to its widely accepted geometric origin in complex porous media. Our results, together with previous findings, provide evidence that anomalous transport is indeed ubiquitous in porous media. Consequently, attentions should be paid when modelling solute transport by the classical advection-diffusion equation, which could lead to systematic error.

Pore-scale statistics of flow and transport through porous media

Flow in porous media is known to be largely affected by pore morphology. In this work, we investigate the effects of pore geometry on the transport and spatial correlations of flow through porous media in two distinct pore structures arising from three-dimensional assemblies of overlapping and nonoverlapping spheres. Using high-resolution direct numerical simulations (DNS), we perform Eulerian and Lagrangian analysis of the flow and transport characteristics in porous media. We show that the Eulerian velocity distributions change from nearly exponential to Gaussian distributions as porosity increases. A stretched exponential distribution can be used to represent this behavior for a wide range of porosities. Evolution of Lagrangian velocities is studied for the uniform injection rule. Evaluation of tortuosity and trajectory length distributions of each porous medium shows that the model of overlapping spheres results in higher tortuosity and more skewed trajectory length distributions compared to the model of nonoverlapping spheres. Wider velocity distribution and higher tortuosity for overlapping spheres model give rise to non-Fickian transport while transport in nonoverlapping spheres model is found to be Fickian. Particularly, for overlapping spheres model our analysis of first-passage time distribution shows that the transport is very similar to those observed for sandstone. Finally, using three-dimensional (3D) velocity field obtained by DNS at the pore-scale, we quantitatively show that despite the randomness of pore-space, the spatially fluctuating velocity field and the 3D pore-space distribution are strongly correlated for a range of porous media from relatively homogeneous monodisperse sphere packs to Castlegate sandstone.

Discrete unified gas kinetic scheme for all Knudsen number flows. III. Binary gas mixtures of Maxwell molecules

Recently a discrete unified gas kinetic scheme (DUGKS) in a finite-volume formulation based on the Boltzmann model equation has been developed for gas flows in all flow regimes. The original DUGKS is designed for flows of single-species gases. In this work, we extend the DUGKS to flows of binary gas mixtures of Maxwell molecules based on the Andries-Aoki-Perthame kinetic model [P. Andries et al., J. Stat. Phys. 106, 993 (2002)JSTPBS0022-471510.1023/A:1014033703134. A particular feature of the method is that the flux at each cell interface is evaluated based on the characteristic solution of the kinetic equation itself; thus the numerical dissipation is low in comparison with that using direct reconstruction. Furthermore, the implicit treatment of the collision term enables the time step to be free from the restriction of the relaxation time. Unlike the DUGKS for single-species flows, a nonlinear system must be solved to determine the interaction parameters appearing in the equilibrium distribution function, which can be obtained analytically for Maxwell molecules. Several tests are performed to validate the scheme, including the shock structure problem under different Mach numbers and molar concentrations, the channel flow driven by a small gradient of pressure, temperature, or concentration, the plane Couette flow, and the shear driven cavity flow under different mass ratios and molar concentrations. The results are compared with those from other reliable numerical methods. The results show that the proposed scheme is an effective and reliable method for binary gas mixtures in all flow regimes.

A robust upscaling of the effective particle deposition rate in porous media

In the upscaling from pore to continuum (Darcy) scale, reaction and deposition phenomena at the solid-liquid interface of a porous medium have to be represented by macroscopic reaction source terms. The effective rates can be computed, in the case of periodic media, from three-dimensional microscopic simulations of the periodic cell. Several computational and semi-analytical models have been studied in the field of colloid filtration to describe this problem. They typically rely on effective deposition rates defined by complex fitting procedures, neglecting the advection-diffusion interplay, the pore-scale flow complexity, and assuming slow reactions (or large Péclet numbers). Therefore, when these rates are inserted into general macroscopic transport equations, they can lead to several conceptual inconsistencies and significant errors. To study more accurately the dependence of deposition on the flow parameters, in this work we advocate a clear distinction between the surface processes (that altogether defines the so-called attachment efficiency), and the pore-scale processes. With this approach, valid when colloidal particles are small enough, we study Brownian and gravity-driven deposition on a face-centred cubic (FCC) arrangement of spherical grains, and define a robust upscaling based on a linear effective reaction rate. The case of partial deposition, defined by an attachment probability, is studied and the limit of perfect sink is retrieved as a particular case. We introduce a novel upscaling approach and a particularly convenient computational setup that allows the direct computation of the asymptotic stationary value of effective rates. This allows to drastically reduce the computational domain down to the scale of the single repeating periodic unit. The savings are ever more noticeable in the case of higher Péclet numbers, when larger physical times are needed to reach the asymptotic regime and thus, equivalently, much larger computational domain and simulation time would be needed in a traditional setup. We show how this new definition of deposition rate is more robust and extendable to the whole range of Péclet numbers; it also is consistent with the classical heat and mass transfer literature.

Longitudinal dispersion coefficients for numerical modeling of groundwater solute transport in heterogeneous formations

Most recent research on hydrodynamic dispersion in porous media has focused on whole-domain dispersion while other research is largely on laboratory-scale dispersion. This work focuses on the contribution of a single block in a numerical model to dispersion. Variability of fluid velocity and concentration within a block is not resolved and the combined spreading effect is approximated using resolved quantities and macroscopic parameters. This applies whether the formation is modeled as homogeneous or discretized into homogeneous blocks but the emphasis here being on the latter. The process of dispersion is typically described through the Fickian model, i.e., the dispersive flux is proportional to the gradient of the resolved concentration, commonly with the Scheidegger parameterization, which is a particular way to compute the dispersion coefficients utilizing dispersivity coefficients. Although such parameterization is by far the most commonly used in solute transport applications, its validity has been questioned. Here, our goal is to investigate the effects of heterogeneity and mass transfer limitations on block-scale longitudinal dispersion and to evaluate under which conditions the Scheidegger parameterization is valid. We compute the relaxation time or memory of the system; changes in time with periods larger than the relaxation time are gradually leading to a condition of local equilibrium under which dispersion is Fickian. The method we use requires the solution of a steady-state advection-dispersion equation, and thus is computationally efficient, and applicable to any heterogeneous hydraulic conductivity K field without requiring statistical or structural assumptions. The method was validated by comparing with other approaches such as the moment analysis and the first order perturbation method. We investigate the impact of heterogeneity, both in degree and structure, on the longitudinal dispersion coefficient and then discuss the role of local dispersion and mass transfer limitations, i.e., the exchange of mass between the permeable matrix and the low permeability inclusions. We illustrate the physical meaning of the method and we show how the block longitudinal dispersivity approaches, under certain conditions, the Scheidegger limit at large Péclet numbers. Lastly, we discuss the potential and limitations of the method to accurately describe dispersion in solute transport applications in heterogeneous aquifers.

Recirculation zones induce non-Fickian transport in three-dimensional periodic porous media

In this work, the influence of pore space geometry on solute transport in porous media is investigated performing computational fluid dynamics pore-scale simulations of fluid flow and solute transport. The three-dimensional periodic domains are obtained from three different pore structure configurations, namely, face-centered-cubic (fcc), body-centered-cubic (bcc), and sphere-in-cube (sic) arrangements of spherical grains. Although transport simulations are performed with media having the same grain size and the same porosity (in fcc and bcc configurations), the resulting breakthrough curves present noteworthy differences, such as enhanced tailing. The cause of such differences is ascribed to the presence of recirculation zones, even at low Reynolds numbers. Various methods to readily identify recirculation zones and quantify their magnitude using pore-scale data are proposed. The information gained from this analysis is then used to define macroscale models able to provide an appropriate description of the observed anomalous transport. A mass transfer model is applied to estimate relevant macroscale parameters (hydrodynamic dispersion above all) and their spatial variation in the medium; a functional relation describing the spatial variation of such macroscale parameters is then proposed.

Experimental investigation of transverse mixing in porous media under helical flow conditions

Plume dilution and transverse mixing can be considerably enhanced by helical flow occurring in three-dimensional heterogeneous anisotropic porous media. In this study, we perform tracer experiments in a fully three-dimensional flow-through chamber to investigate the effects of helical flow on plume spiraling and deformation, as well as on its dilution. Porous media were packed in angled stripes of materials with different grain sizes to create blocks with macroscopically anisotropic hydraulic conductivity, which caused helical flows. Steady-state transport experiments were carried out by continuously injecting dye tracers at different inlet ports. High-resolution measurements of concentration and flow rates were performed at 49 outlet ports. These measurements allowed quantifying the spreading and dilution of the solute plumes at the outlet cross section. Direct evidence of plume spiraling and visual proof of helical flow was obtained by freezing and slicing the porous media at different cross sections and observing the dye-tracer distribution. We simulated flow and transport to interpret our experimental observations and investigate the effects of helical flow on mixing-controlled reactive transport. The simulation results were evaluated using metrics of reactive mixing such as the critical dilution index and the length of continuously injected steady-state plumes. The results show considerable reaction enhancement, quantified by the remarkable decrease of reactive plume lengths (up to four times) in helical flows compared to analogous scenarios in uniform flows.

Power-exponential velocity distributions in disordered porous media

Velocity distribution functions link the micro- and macro-level theories of fluid flow through porous media. Here we study them for the fluid absolute velocity and its longitudinal and lateral components relative to the macroscopic flow direction in a model of a random porous medium. We claim that all distributions follow the power-exponential law controlled by an exponent γ and a shift parameter u_{0} and examine how these parameters depend on the porosity. We find that γ has a universal value 1/2 at the percolation threshold and grows with the porosity, but never exceeds 2.

Experimental Evidence of Helical Flow in Porous Media

Helical flow leads to deformation of solute plumes and enhances transverse mixing in porous media. We present experiments in which macroscopic helical flow is created by arranging different materials to obtain an anisotropic macroscopic permeability tensor with spatially variable orientation. The resulting helical flow entails twisting streamlines which cause a significant increase in lateral mass exchange and thus a large enhancement of plume dilution (up to 235%) compared to transport in homogenous media. The setup may be used to effectively mix solutes in parallel streams similarly to static mixers, but in porous media.

Nanosafety forum for young scientists: a meeting report

In October of 2014, a meeting jointly organized by the EU Nanosafety Cluster and the COST Action TD 1204, was held on the beautiful island of Ortygia in Syracuse (Sicily). The meeting was specifically conceived to give the opportunity to young researchers in the field of nanotoxicology to present and discuss the results of their research. Presentations were divided into eight sessions over 2 days, reflecting the eight working groups of the Nanosafety Cluster. This report gives a description of the meeting activities and a summary of the data presented there.