Abstract
We combine results of high-resolution microfluidic experiments with extensive numerical simulations to show how the flow patterns inside a “swiss-cheese” type of pore geometry can be systematically controlled through the intrinsic rheological properties of the fluid. Precisely, our analysis reveals that the velocity field in the interstitial pore space tends to display enhanced channeling under certain flow conditions. This observed flow “localization”, quantified by the spatial distribution of kinetic energy, can then be explained in terms of the strong interplay between the disordered geometry of the pore space and the nonlinear rheology of the fluid. Our results disclose the possibility that the constitutive properties of the fluid can enhance the performance of chemical reactors and chromatographic devices through control of the channeling patterns inside disordered porous media.
Introduction
Flow through porous media is of great interest in chemical engineering, physics, and biology [1–3]. Previous studies have shown that the disordered characteristics of the pore structure naturally leads to heterogeneous flow patterns [4–8] and preferential channeling [9–11]. Understanding how to control and manipulate these flow patterns can help to optimize catalysts [12, 13] or chromatographic devices [14, 15], and allows to steer chemical reactions inside the porous medium itself [16–18].
In order to understand the physics of important problems like, for example, blood flow through the kidney [19] or oil flow through porous rocks [20, 21], one must also consider the nonlinear constitutive behavior of the fluids involved in these processes. Technological applications which make use of non-Newtonian fluids are ubiquitous nowadays [22–25]. It is, for instance, the case of shear-thinning solvents that are present in dropless paints [26], shear-thickening fluids being used as active dampers [27] and hybrid fluids as components of enhanced body armors [28]. While Newtonian flows in irregular media have been extensively investigated theoretically and confirmed by many experiments, the study of non-Newtonian fluids lack a generalized framework due to their diverse constitutive nature. Non-Newtonian flows through porous media have mainly been studied theoretically [29, 30] and through numerical simulations [31, 32], where the main focus of interest was to find non-Darcian models for the flow of generalized Newtonian fluids [30, 33–37]. In the particular case of power-law fluids, it has been shown that, in spite of the nonlinear nature of the fluid’s rheology and the geometrical complexity of the pore volume, the general behavior of the system can still be quantified in terms of a universal permeability extending over a broad range of Reynolds conditions and power-law exponents [38].
However, quantitative experiments with non-Newtonian materials which go beyond simple bulk measurements [39, 40] are scarce [41] because the design of the experimental pore geometry and the operating conditions need to be adjusted in order to match the nonlinear constitutive regime of the fluid’s rheology. Here we combine the results of microfluidic experiments [42] with fluid dynamics simulations to demonstrate how the nonlinear rheological properties of a fluid can be effectively exploited in order to control the macroscopic transport properties of a flow through the external operational flow conditions. These results have important consequences for the design of chemical reactors and chromatographic systems as well as for the enhancement of oil recovery and transport in porous media in general.
Under steady-state conditions the motion of an incompressible fluid through the interstitial space of a porous medium is described by mass and momentum conservation, respectively,together with appropriate boundary conditions. The variables ϱ, and p are the fluid’s density, velocity and pressure, and is the deviatoric stress tensor which depends on the fluid’s rheology. For many fluids, this constitutive relation is well described by a simple linear rheology , where is the shear strain rate tensor [43] and the proportionality constant μ defines the kinematic viscosity. Examples of these so-called Newtonian fluids are water, light oil and most diluted gases. However, many fluids present in industrial products, biology and environmental flows obey much more complex nonlinear constitutive laws [23, 25, 44]. These fluids are called non-Newtonian fluids. The constitutive behavior of most non-Newtonian fluids can be described by a generalization of the Newtonian relation, namely,Here the apparent viscosity is a nonlinear function of the second principal invariant of the shear strain rate tensor alone [43]. Examples of fluids—which are often called generalized Newtonian fluids—are colloidal suspensions, protein or polymeric mixtures, heavy petroleum, blood or debris flows, only to mention a few [45–48].
Our analysis is based on experimental results from the setup presented recently in Eberhard et al. [42]. Their main purpose was to map the local viscosity of a non-Newtonian flow in a porous microfluidic channel by means of a high-resolution technique of image velocimetry, namely, Ghost Particle Velocimetry (GVP) [49]. The geometry of the microfluidic chip is shown in Figure 1. As non-Newtonian fluid we used a 0.5 wt% xanthan gum solution, which is a polysaccharide mainly found in food industry [50] and enhanced oil recovery [22, 24]. It has a shear-thinning rheology which changes its apparent viscosity over several orders of magnitude. While polymeric solutions often show viscoelastic behavior [51], the concentration of xanthan gum in our experimental solution was so low that no measurable elastic behavior could be observed during the experiment. The rheology of xanthan gum closely follows a Carreau model,approaching the viscosity of the solvent (water) in the limit of very high shear [52]. Conversely, for low shear, Eq. 4 reduces to , corresponding to a constant viscosity . At an intermediate range of shear rates, the fluid follows a power-law relation with in our specific case. The remaining Carreau parameter was determined using a nonlinear least-square fit to the experimentally measured values [42]. As shown in Figure 1, the experimental pore structure consisted of a microfluidic device containing a quasi-2D porous medium of size 30 mm × 15 mm and depth of 100 µm. It contains pillars of radius 100 µm that are randomly allocated and can overlap, forming a “swiss-cheese” pore geometry with void fraction approximately equal to 0.8.
FIGURE 1
Figure 2 compares the velocimetry measurements obtained from ref. 42 in a section of the mid plane of the microfluidic chip with those obtained from numerical simulations calculated with exactly the same pore geometry, fluid properties and flow conditions. More precisely, the flow rates for the presented cases are qin = 0.05 µL/min (Figure 2A) and qin = 5 µL/min (Figure 2B) for the xanthan case, and qin = 5 µL/min for the measurement with water (Figure 2C). The color scale has been normalized to the 95% quantile of the velocity distribution to facilitate the comparison of the flow fields at different flow rates. Although differences between the experimentally measured and simulated velocity fields can be visually detected, they are mostly local and could be explained by the natural difficulties of exactly reproducing in the mathematical model the detailed features of the flow, the fluid rheology, and the flow operational conditions. To perform numerical simulations, the computational mesh was generated by capturing the two-dimensional technical drawing of the device geometry into Ansys’ meshing module [53]. In the horizontal plane, we first created an unstructured quadrilateral mesh with an average cell size of . The three-dimensional structure was obtained by extruding ten vertical layers, resulting in a computational mesh composed of approximately 30 million unstructured hexahedral cells (corresponding to roughly 20 µm3/cell), which was then imported into Ansys FluentTM [53]. Non-slip boundary conditions were applied on all solid walls of the microfluidic chip, which is a reasonable assumption on surfaces without hyper-hydrophilic coatings and at scales much larger than the polymer coil size [54]. The fluid was injected via a constant velocity inlet corresponding to the inflow rate reported in the experiment. The density of the non-Newtonian fluid used in the computational simulations matches exactly the experimental value of the xanthan solution, namely, ρxan = 1.0 g/cm3. As for the rheology of the fluid, we used Eq. 4 to interpolate the local viscosity values obtained from the independent rheometer experiments [42]. Finally, the steady-state flow solution in terms of the velocity and pressure fields was calculated using a second-order integration scheme and convergence was achieved if the residuals reached a threshold of .
FIGURE 2
In order to check the variability of our results with respect to the disorder level of the pore space, numerical simulations have also been performed with three additional realizations of the swiss-cheese geometry, but keeping the same physico-chemical properties of the fluid, operational parameters of the flow, and boundary conditions. Considering the independence of the rheometry and velocimetry measurements, the excellent agreement between results from the numerical model and experiments (Figures 2A–C) clearly demonstrates the global consistency of our methodological approach. In order to highlight the differences between the Newtonian and non-Newtonian flows and the tendency for stronger localization in the non-Newtonian flow, we show in Figure 3A the contour plot of the ratio between the local velocity magnitudes measured with the xanthan solution and water normalized by their respective mean velocities. In both cases, the applied flow rate was set to qin = 5 µL/min.
FIGURE 3
Analysis
The channeling effect present in the flow fields shown in Figure 2 can be statistically quantified in terms of their spatial distributions of kinetic energy . This is performed here in terms of a measure utilized in previous studies on localization of vibrational modes in harmonic chains [55], namely, the participation number defined as,where the total volume of the fluid in a domain is given by . Thus the participation ratio varies between , corresponding to a limiting state of equal partition of kinetic energy and the value for a sufficiently large system , indicating strong localization, namely, the presence of intense channeling effects in the flow field [56]. For reference, the Reynolds number is defined here as,where v is the mean velocity at the entrance of the pore zone and is the diameter of the solid obstacles. Figure 3B shows how the participation ratio obtained with our computational model varies as a function of Reynolds number. The thick black curve is obtained by averaging the results from the pore geometry used in the microfluidic experiment (dark gray) with four additional realizations of the “swiss-cheese” pore geometry (light gray markers) having the same porosity. The spread of the markers therefore gives a good indication for the statistical variability of the participation ratio for different pore geometries. The values of the participation ratios obtained from the two experimentally measured velocity fields are marked with yellow and blue stars, respectively, and are in good agreement with the numerical calculations.
At low Reynolds numbers, the participation ratio is practically constant at , similar to a Stokes flow with fixed viscosity . By increasing the flow rate, the local shear within the interstitial pore space also generally increases to eventually reach a point where its range of variability overlaps with the range in which the rheology of the xanthan gum solution follows a power-law behavior. At this point, the local viscosity spans over a wide range of values, leading to a drop in the participation ratio by almost 20% at a Reynolds numbers around .
Interestingly, while the absolute value of the participation ratio varies slightly from realization to realization, the location of the participation ratio minimum is determined by the fluid’s rheology and does not seem to be influenced by the details of the pore geometry. More precisely, we find for the two experimentally measured velocity fields a participation ratio of for qin = 0.05 µL/min and for qin = 5 µL/min. For the simulations, the mean participation ratios averaged over four realizations yield for qin = 0.05 µL/min and for qin = 5 µL/min. The participation ratios of the two simulations which share the same pore geometry as the experimental setup are (qin = 0.05 µL/min) and (qin = 5 µL/min). After reaching the minimum, where the flow is most heterogeneously distributed in the porous medium, the participation ratio increases again as the flow rate pushes the fluid’s rheology beyond the power-law regime, ultimately approaching the flow pattern of water. In this limiting case, the Newtonian behavior is recovered, but inertial effects on the flow should prevent the participation ratio to reach the same value obtained for very low Reynolds numbers, namely, .
A question that naturally arises is how rheology influences the flow’s heterogeneity in space. At low and high Reynolds number, the Carreau fluid has an almost constant viscosity equal to the low and high shear limits and , respectively. In the intermediate regime, however, the local viscosity covers a broad spectrum of values [42]. In this case, both experimental and simulation results reveal that the interplay between the disordered geometry of the pore space and the fluid rheology leads to a larger flow heterogeneity and therefore to a stronger localization pattern.
Conclusion
It is well accepted that flow and transport processes in porous media are fundamentally controlled by the complex interplay between the fluid and the pore space structure. Here we showed that these processes can be tailored by tuning the rheology of the fluid. Precisely, the heterogeneity of the pore scale structure of the medium causes a high variability of shear rates that, when matched with the nonlinear viscosity window of the non-Newtonian fluid, can substantially enhance macroscopic properties of the system like the participation ratio. These effects may be exploited to improve filters and catalysts or to enhance chemical reactions by spreading the transported chemicals more uniformly throughout the pore space. In particular, localization should have a deleterious influence on the effectiveness of catalysts subjected to flow, for example, in a packed bed chemical reaction. Precisely, the preferential channeling at the minimum of the participation number should be avoided to maximize the activity of the surface area available for reaction in the system.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
HS and JA designed the research. UE carried out the experiments with input from ES, MH, and JJ-M All authors discussed the results. HS and JA wrote the paper with input from all other authors.
Acknowledgments
We acknowledge Prof. R. Stocker and E. Burmeister from the Department of Civil, Environmental and Geomatic Engineering at ETH Zurich for kindly providing access to laboratory equipment and material. ES acknowledges SNSF PRIMA Grant 179834. JA acknowledges financial support from the Brazilian agencies CNPq, CAPES and FUNCAP. HS and JA acknowledge financial support from Petrobras (“Física do Petróleo em Meios Porosos”, Project Number: F0185. RA acknowledges financial support from the Brazilian agency CNPq and the INCT-SC Project.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
DullienFAL. Porous media: fluid transport and pore structure. San Diego, NY: Academic Press (2012).
2.
SahimiM. Flow and transport in porous media and fractured rock: from classical methods to modern approaches. Weinheim, Germany: John Wiley & Sons (2011).
3.
KhaledARAVafaiK. The role of porous media in modeling flow and heat transfer in biological tissues. Int J Heat Mass Tran (2003) 46:4989–5003. 10.1016/s0017-9310(03)00301-6
4.
WarrenJEPriceHS. Flow in heterogeneous porous media. Soc Petrol Eng J (1961) 1:153–69.
5.
BearJ. Dynamics of fluids in porous materials. New York: Courier Corporation (2013).
6.
DavidC. Geometry of flow paths for fluid transport in rocks. J Geophys Res (1993) 98:12267–78.
7.
AndradeJSStreetDAShinoharaTShibusaYAraiY. Percolation disorder in viscous and nonviscous flow through porous media. Phys Rev E (1995) 51:5725–31.
8.
SeyboldHJCarmonaHALeonardo FilhoFAAraújoADNepomuceno FilhoFAndradeJS. Flow through three-dimensional self-affine fractures. Phys Rev Fluids (2020) 5:104101. 10.1103/physrevfluids.5.104101
9.
Bruderer-WengCCowiePBernabéYMainI. Relating flow channelling to tracer dispersion in heterogeneous networks. Adv Water Resour (2004) 27:843–55. 10.1016/j.advwatres.2004.05.001
10.
AndradeJSAlmeidaMPMendes FilhoJHavlinSSukiBStanleyHE. Fluid flow through porous media: the role of stagnant zones. Phys Rev Lett (1997) 79:3901.
11.
TsangCFNeretnieksI. Flow channeling in heterogeneous fractured rocks. Rev Geophys (1998) 36:275–98.
12.
WangGJohannessenEKleijnCRde LeeuwSWCoppensMO. Optimizing transport in nanostructured catalysts: a computational study. Chem Eng Sci (2007) 62:5110–6. 10.1016/j.ces.2007.01.046
13.
DavisME. Ordered porous materials for emerging applications. Nature (2002) 417:813–21. 10.1038/nature00785
14.
BillenJDesmetG. Understanding and design of existing and future chromatographic support formats. J Chromatogr A (2007) 1168:73–2. 10.1016/j.chroma.2007.07.069
15.
TennikovMBGazdinaNVTennikovaTBSvecF. Effect of porous structure of macroporous polymer supports on resolution in high-performance membrane chromatography of proteins. J Chromatogr A (1998) 798:55–64. 10.1016/s0021-9673(97)00873-x
16.
RubinJ. Transport of reacting solutes in porous media: relation between mathematical nature of problem formulation and chemical nature of reactions. Water Resour Res (1983) 19:1231–52. 10.1029/wr019i005p01231
17.
KeilFJ. Diffusion and reaction in porous networks. Catal Today (1999) 53:245–58. 10.1016/s0920-5861(99)00119-4
18.
VafaiK. Handbook of porous media. Boca Raton: CRC Press (2015).
19.
BrimbleKSMcFarlaneAWinegardNCrowtherMChurchillDN. Effect of chronic kidney disease on red blood cell rheology. Clin Hemorheol Microcirc (2006) 34:411–20.
20.
McCainWDJr. Properties of petroleum fluids. Tusla: PennWell Corporation (2017).
21.
XieCLvWWangM. Shear-thinning or shear-thickening fluid for better eor?–a direct pore-scale study. J Petrol Sci Eng (2018) 161:683–91. 10.1016/j.petrol.2017.11.049
22.
SorbieKS. Polymer-improved oil recovery. Edinburgh: Springer Science & Business Media (2013).
23.
LaiSKWangYYWirtzDHanesJ. Micro- and macrorheology of mucus. Adv Drug Deliv Rev (2009) 61:86–100. 10.1016/j.addr.2008.09.012
24.
SandvikEIMaerkerJM. Application of xanthan gum for enhanced oil recovery. Am Chem Soc (1977) 45:242–64. 10.1021/bk-1977-0045.ch019
25.
LópezOVCastilloLANinagoMDCiolinoAEVillarMA. Modified starches used as additives in enhanced oil recovery (EOR) In: GoyanesSND’AccorsoNB, editors Industrial applications of renewable biomass products: past, present and future. Cham: Springer International Publishing (2017).
26.
ReuversAJ. Control of rheology of water-borne paints using associative thickeners. Prog Org Coating (1999) 35:171–81.
27.
WangFX. Flow analysis and modeling of field-controllable, electro-and magneto-rheological fluid dampers. J Appl Mech (2007) 74:13–22. 10.1115/1.2166649
28.
MajumdarAButolaSBSrivastavaA. Optimal designing of soft body armour materials using shear thickening fluid. Mater Des (2013) 46:191–8. 10.1016/j.matdes.2012.10.018
29.
SahimiM. Nonlinear transport processes in disordered media. AIChE J (1993) 39:369–86. 10.1002/aic.690390302
30.
ShahCBYortsosYC. Aspects of flow of power-law fluids in porous media. AIChE J (1995) 41:1099–112. 10.1002/aic.690410506
31.
DeSKuipersJAMPetersEAJFPaddingJT. Viscoelastic flow past mono- and bidisperse random arrays of cylinders: flow resistance, topology and normal stress distribution. Soft Matter (2017) 13:9138–46. 10.1039/c7sm01818e
32.
MokhtariZZippeliusA. Dynamics of active filaments in porous media. Phys Rev Lett (2019) 123:028001. 10.1103/physrevlett.123.028001
33.
CannellaWJHuhCSerightRS. Prediction of xanthan rheology in porous media. In: SPE annual technical conference and exhibition; 1988 Oct 2–5; Houston TX. Society of Petroleum Engineers (1988). 10.2118/18089-MS
34.
TsakiroglouCD. A methodology for the derivation of non-darcian models for the flow of generalized Newtonian fluids in porous media. J Non Newton Fluid (2002) 105:79–110. 10.1016/s0377-0257(02)00066-6
35.
SochiTBluntMJ. Pore-scale network modeling of ellis and herschel–bulkley fluids. J Petrol Sci Eng (2008) 60:105–24. 10.1016/j.petrol.2007.05.009
36.
BergSvan WunnikJ. Shear rate determination from pore-scale flow fields. Transp Porous Med (2017) 117:229–46. 10.1007/s11242-017-0830-3
37.
EberhardUSeyboldHJFloriancicMBertschPJiménez-MartínezJAndradeJSJret alDetermination of the effective viscosity of non-Newtonian fluids flowing through porous media. Front Phys (2019) 7:71. 10.3389/fphy.2019.00071
38.
MoraisAFSeyboldHHerrmannHJAndradeJS. Non-Newtonian fluid flow through three-dimensional disordered porous media. Phys Rev Lett (2009) 103:194502. 10.1103/PhysRevLett.103.194502
39.
PerrinCLTardyPMSorbieKSCrawshawJC. Experimental and modeling study of Newtonian and non-Newtonian fluid flow in pore network micromodels. J Colloid Interface Sci (2006) 295:542–50. 10.1016/j.jcis.2005.09.012
40.
de CastroARRadillaG. Non-darcian flow of shear-thinning fluids through packed beads: experiments and predictions using forchheimer’s law and ergun’s equation. Adv Water Resour (2017) 100:35–47. 10.1016/j.advwatres.2016.12.009
41.
HopkinsCCHawardSJShenAQ. Tristability in viscoelastic flow past side-by-side microcylinders. arXiv:2010 14749v1 (2020).
42.
EberhardUSeyboldHJSecchiEJiménez-MartínezJRhüsPOfnerAet alMapping the local viscosity of heterogeneous non- Newtonian flows. Sci Rep (2020). 10:11733. 10.1038/s41598-020-68545-7
43.
BatchelorGK. An introduction to fluid dynamics. Cambridge: Cambridge University Press (2000).
44.
AstaritaGMarrucciG. Principles of non-Newtonian fluid mechanics. London, New York: McGraw-Hill Companies (1974).
45.
RoyerJRBlairDLHudsonSD. Rheological signature of frictional interactions in shear thickening suspensions. Phys Rev Lett (2016) 116:188301. 10.1103/PhysRevLett.116.188301
46.
ChhabraRP. Bubbles, drops, and particles in non-Newtonian fluids. Boca Raton: CRC Press (2006).
47.
BaskurtOKMeiselmanHJ. Blood rheology and hemodynamics. New York (2003).
48.
CoussotPMeunierM. Recognition, classification and mechanical description of debris flows. Earth Sci Rev (1996) 40:209–27. 10.1016/0012-8252(95)00065-8
49.
BuzzaccaroSSecchiEPiazzaR. Ghost particle velocimetry: accurate 3d flow visualization using standard lab equipment. Phys Rev Lett (2013) 111, 048101. 10.1103/PhysRevLett.111.048101
50.
KatzbauerB. Properties and applications of xanthan gum. Polym Degrad Stabil (1998) 59:81–4. 10.1016/s0141-3910(97)00180-8
51.
WalkamaDMWaisbordNGuastoJS. Disorder suppresses chaos in viscoelastic flows. Phys Rev Lett (2020) 124:164501. 10.1103/PhysRevLett.124.164501
52.
BewersdorffHWSinghRP. Rheological and drag reduction characteristics of xanthan gum solutions. Rheol Acta (1988) 27:617–27. 10.1007/bf01337457
53.
AnsysA. Workbench user manual. Canonsburg: ANSYS (2019).
54.
HaaseASWoodJASprakelLMLammertinkRG. Inelastic non-Newtonian flow over heterogeneously slippery surfaces. Phys Rev E (2017) 95:023105. 10.1103/PhysRevE.95.023105
55.
RussSSapovalB. Anomalous viscous damping of vibrations of fractal percolation clusters. Phys Rev Lett (1994) 73:1570. 10.1103/PhysRevLett.73.1570
56.
AndradeJSCostaUMSAlmeidaMPMakseHAStanleyHE. Inertial effects on fluid flow through disordered porous media. Phys Rev Lett (1999) 82:5249–52. 10.1103/physrevlett.82.5249
Summary
Keywords
localization, microfluidics, particle velocimetry, non-newtonian fluids, porous media
Citation
Seybold HJ, Eberhard U, Secchi E, Cisne Jr. RLC, Jiménez-Martínez J, Andrade RFS, Araújo AD, Holzner M and Andrade Jr. JS (2021) Localization in Flow of Non-Newtonian Fluids Through Disordered Porous Media. Front. Phys. 9:635051. doi: 10.3389/fphy.2021.635051
Received
29 November 2020
Accepted
04 January 2021
Published
16 February 2021
Volume
9 - 2021
Edited by
Ferenc Kun, University of Debrecen, Hungary
Reviewed by
Taotao Fu, Tianjin University, China
Simon Haward, Okinawa Institute of Science and Technology Graduate University, Japan
Sandro Longo, University of Parma, Italy
Updates
Copyright
© 2021 Seybold, Eberhard, Secchi, Cisne, Jiménez-Martínez, Andrade, Araújo, Holzner and Andrade.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: H. J. Seybold, hseybold@ethz.ch
This article was submitted to Interdisciplinary Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.