Comment on “Frequency-dependent dispersion in porous media”

In a recent paper, Valdes-Parada and Alvarez-Ramirez [Phys. Rev. E 84 031201 (2011)] used the technique of volume averaging to derive a frequency-dependent dispersion tensor, D*(w), the goal of which is to describe solute transport in porous media undergoing periodic processes. We describe two issues with this dispersion tensor that become particularly apparent in the time domain. Firstly, we demonstrate that the definition of D*(w), and equivalently of D*(t), is erroneous, and derive its correct counterpart, D*c(t). With this modification, the approach becomes strictly equivalent to the one devised by Moyne [Adv. Water Res. 20 63 (1997)]. Secondly, we show that the term "frequency-dependent dispersion" is misleading in this case, because D*(t) and D*c(t) do not depend on the process operating frequency. The study carried out by Valdes-Parada and Alvarez-Ramirez represents a frequency-analysis of the relaxation of the dispersion coefficients towards their steady-state, independently of any periodic operation or excitation

A two-pressure model for slightly compressible single phase flow in bi-structured porous media

Problems involving flow in porous media are ubiquitous in many natural and engineered systems. Mathematical modeling of such systems often relies on homogenization of pore-scale equations and macroscale continuum descriptions. For single phase flow, Stokes equations at the pore-scale are generally approximated by Darcy’s law at a larger scale. In this work, we develop an alternative model to Darcy’s law that can be used to describe slightly compressible single phase flow within bi-structured porous media. We use the method of volume averaging to upscale mass and momentum balance equations with the fluid phase split into two fictitious domains. The resulting macroscale model combines two coupled equations for average pressures with regional Darcy’s laws for velocities. In these equations, effective parameters are expressed via integrals of mapping variables that solve boundary value problems over a representative unit cell. Finally, we illustrate the behaviour of these equations in a two-dimensional model porous medium and validate our approach by comparing solutions of the homogenized equations with computations of the exact microscale problem

Correspondence between one- and two-equation models for solute\ud transport in two-region heterogeneous porous media

In this work, we study the transient behavior of upscaled models for solute transport in two-region porous media. We focus on the following three models: (1) a time non-local, two-equation model (2eq-nlt). This model does not rely on time constraints and, therefore, is particularly useful in the short-time regime, when the time scale of interest (t) is smaller than the characteristic time (T1) for the relaxation of the effective macroscale parameters (i.e., when t ≤ T1); (2) a time local, two-equation model (2eq). This model can be adopted when (t) is significantly larger than (T1) (i.e., when t » T1); and (3) a one-equation, time-asymptotic formulation (1eq∞). This model can be adopted when (t) is significantly larger than the time scale (T2) associated with exchange processes between the two regions (i.e., when t » T2). In order to obtain some physical insight into this transient behavior, we combine a theoretical approach based on the analysis of spatial moments with numerical and analytical results in simple cases. The main result of this paper is to show that there is weak long-time convergence of the solution of (2eq) toward the solution of (1eq∞) in terms of standardized moments but, interestingly, not in terms of centered moments. Physically, our interpretation of this result is that the spreading of the solute is dominating higher order non-zero perturbations in the asymptotic regime

Validity of the Cauchy-Born rule applied to discrete cellular-scale models of biological tissues

The development of new models of biological tissues that consider cells in a discrete manner is becoming increasingly popular as an alternative to PDE-based continuum methods, although formal relationships between the discrete and continuum frameworks remain to be established. For crystal mechanics, the discrete-to-continuum bridge is often made by assuming that local atom displacements can be mapped homogeneously from the mesoscale deformation gradient, an assumption known as the Cauchy-Born rule (CBR). Although the CBR does not hold exactly for non-crystalline materials, it may still be used as a first order approximation for analytic calculations of effective stresses or strain energies. In this work, our goal is to investigate numerically the applicability of the CBR to 2-D cellular-scale models by assessing the mechanical behaviour of model biological tissues, including crystalline (honeycomb) and non-crystalline reference states. The numerical procedure consists in precribing an affine deformation on the boundary cells and computing the position of internal cells. The position of internal cells is then compared with the prediction of the CBR and an average deviation is calculated in the strain domain. For centre-based models, we show that the CBR holds exactly when the deformation gradient is relatively small and the reference stress-free configuration is defined by a honeycomb lattice. We show further that the CBR may be used approximately when the reference state is perturbed from the honeycomb configuration. By contrast, for vertex-based models, a similar analysis reveals that the CBR does not provide a good representation of the tissue mechanics, even when the reference configuration is defined by a honeycomb lattice. The paper concludes with a discussion of the implications of these results for concurrent discrete/continuous modelling, adaptation of atom-to-continuum (AtC) techniques to biological tissues and model classification

Solute transport within porous biofilms: diffusion or dispersion?

Many microorganisms live within surface-associated consortia, termed biofilms, that can form intricate porous structures interspersed with a network of fluid channels. In such systems, transport phenomena, including flow and advection, regulate various aspects of cell behaviour by controllling nutrient supply, evacuation of waste products and permeation of antimicrobial agents. This study presents multiscale analysis of solute transport in these porous biofilms. We start our analysis with a channel-scale description of mass transport and use the method of volume averaging to derive a set of homogenized equations at the biofilmscale. We show that solute transport may be described via two coupled partial differential equations for the averaged concentrations, or telegrapher’s equations. These models are particularly relevant for chemical species, such as some antimicrobial agents, that penetrate cell clusters very slowly. In most cases, especially for nutrients, solute penetration is faster, and transport can be described via an advection-dispersion equation. In this simpler case, the effective diffusion is characterised by a second-order tensor whose components depend on: (1) the topology of the channels’ network; (2) the solute’s diffusion coefficients in the fluid and the cell clusters; (3) hydrodynamic dispersion effects; and (4) an additional dispersion term intrinsic to the two-phase configuration. Although solute transport in biofilms is commonly thought to be diffusion-dominated, this analysis shows that dispersion effects may significantly contribute to transport

Assessing the accuracy of energy turbulent diffusion dispersion correlation in a porous two-fluid model dedicated to PWR core simulations

International audienceCATHARE is a 2-fluid thermal-hydraulic code, capable of simulating thermal and mechanical phenomena occurring in the primary and secondary circuits of Pressurized Water Reactor under a wide variety of accidental situations. One of the medium-term objectives of system code CATHARE-3 is modeling a PWR core at assembly scale to simulate various accidental situations such as the loss of coolant accident (LOCA) and steam line break accident. This requires the monophasic and two-phase models that adapted to the assembly scale. However, there exists 3D models for the whole core and sub-channel scale models, which have a certain degree of validation. For more macroscopic three-dimensional models, we only have global validations without local measurements, which is necessary for the validations of each closure law's separate effects. The objective of my PhD project is improving the sub-channel scale models and developing the assembly scale models in CATHARE-3 system code with the sub-channel scale simulations and experiments results

Hydrodynamic dispersion within porous biofilms

Many microorganisms live within surface-associated consortia, termed biofilms, that can form intricate porous structures interspersed with a network of fluid channels. In such systems, transport phenomena, including flow and advection, regulate various aspects of cell behavior by controlling nutrient supply, evacuation of waste products, and permeation of antimicrobial agents. This study presents multiscale analysis of solute transport in these porous biofilms. We start our analysis with a channel-scale description of mass transport and use the method of volume averaging to derive a set of homogenized equations at the biofilm-scale in the case where the width of the channels is significantly smaller than the thickness of the biofilm. We show that solute transport may be described via two coupled partial differential equations or telegrapher's equations for the averaged concentrations. These models are particularly relevant for chemicals, such as some antimicrobial agents, that penetrate cell clusters very slowly. In most cases, especially for nutrients, solute penetration is faster, and transport can be described via an advection-dispersion equation. In this simpler case, the effective diffusion is characterized by a second-order tensor whose components depend on (1) the topology of the channels' network; (2) the solute's diffusion coefficients in the fluid and the cell clusters; (3) hydrodynamic dispersion effects; and (4) an additional dispersion term intrinsic to the two-phase configuration. Although solute transport in biofilms is commonly thought to be diffusion dominated, this analysis shows that hydrodynamic dispersion effects may significantly contribute to transport

Polymer Flow Through Porous Media: Numerical Prediction of the Contribution of Slip to the Apparent Viscosity.

The flow of polymer solutions in porous media is often described using Darcy’s law with an apparent viscosity capturing the observed thinning or thickening effects. While the macroscale form is well accepted, the fundamentals of the pore-scale mechanisms, their link with the apparent viscosity, and their relative influence are still a matter of debate. Besides the complex effects associated with the rheology of the bulk fluid, the flow is also deeply influenced by the mechanisms occurring close to the solid/liquid interface, where polymer molecules can arrange and interact in a complex manner. In this paper, we focus on a repulsive mechanism, where polymer molecules are pushed away from the interface, yielding a so-called depletion layer in the vicinity of the wall. This depletion layer acts as a lubricating film that may be represented by an effective slip boundary condition. Here, our goal is to provide a simple mean to evaluate the contribution of this slip effect to the apparent viscosity. To do so, we solve the pore-scale flow numerically in idealized porous media with a slip length evaluated analytically in a tube. Besides its simplicity, the advantage of our approach is also that it captures relatively well the apparent viscosity obtained from core-flood experiments, using only a limited number of inputs. Therefore, it may be useful in many applications to rapidly estimate the influence of the depletion layer effect over the macroscale flow and its relative contribution compared to other phenomena, such as non-Newtonian effects