# Localization in Flow of Non-Newtonian Fluids Through Disordered Porous Media

^{1}Institute of Terrestrial Ecosystems, Department of Environmental Systems Science, ETH Zürich, Zurich, Switzerland^{2}Departamento de Física, Universidade Federal do Ceará, Campus do Pici, Fortaleza, Brazil^{3}Institute of Environmental Engineering, Department of Civil, Environmental and Geomatic Engineering, ETH Zürich, Zürich, Switzerland^{4}Instituto Federal do Ceará, Campus Tianguá, Ceará, Brasil^{5}EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Dübendorf, Switzerland^{6}Instituto de Física, Universidade Federal da Bahia, Salvador, Brazil^{7}Swiss Federal Research Institute WSL, Birmensdorf, Switzerland

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 *ϱ*, *p* are the fluid’s density, velocity and pressure, and

Here the apparent viscosity

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)

**FIGURE 1**. Sketch of the experimental setup showing the whole pore geometry together with the simulated flow field in the mid plane for the non-Newtonian case at qin 5 µL/min. The porous region is L 30 mm long, B 15 mm wide and has a height of h 100 µm. The radius of the circular pillars is also 100 mµ. The gray shaded region marks the part of the porous device where the experimental flow velocity measurement has been performed.

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 *q*_{in} = 0.05 µL/min (Figure 2A) and *q*_{in} = 5 µL/min (Figure 2B) for the xanthan case, and *q*_{in} = 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 ^{3}/cell), which was then imported into *Ansys Fluent*^{TM} [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/cm^{3}. 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**. Comparison of the experimentally measured [42] and simulated velocity fields in the mid plane of the porous device for the two non-Newtonian cases at **(A)** qin 0.05 µL/min, **(B)** qin 5 µL/min, and for the Newtonian case **(C)**. Note the different velocity ranges between **(A–C)**. For comparison, the color scales has been normalized to the 95% quantile of the velocity distribution.

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 *q*_{in} = 5 µL/min.

**FIGURE 3**. **(A)** Contour plot of the ratio between the local velocity magnitudes measured with the non-Newtonian (xanthan) and Newtonian (water) fluids, normalized by the corresponding mean velocities in the observation mid planes, both obtained at qin 5 µL/min. **(B)** Participation ratio as a function of Reynolds number. Gray markers label the participation values obtained from the simulations of four different realizations of the pore structure. The realization which corresponds to the experimental device is marked in darker gray. The black solid line was obtained by averaging the four realizations of the “swiss-cheese” pore geometry. The participation ratio calculated from the experimentally measured velocity fields of the xanthan gum solution at q_{in} 0.05 µL/min and q_{in} 5 µL/min are labeled with yellow and blue stars, respectively.

## 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

where the total volume of the fluid in a domain is given by

where *v* is the mean velocity at the entrance of the pore zone and

At low Reynolds numbers, the participation ratio is practically constant at

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 *q*_{in} = 0.05 µL/min and *q*_{in} = 5 µL/min. For the simulations, the mean participation ratios averaged over four realizations yield *q*_{in} = 0.05 µL/min and *q*_{in} = 5 µL/min. The participation ratios of the two simulations which share the same pore geometry as the experimental setup are *q*_{in} = 0.05 µL/min) and *q*_{in} = 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

## 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.

## 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.

## 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.

## 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.

## References

1. Dullien FAL. *Porous media: fluid transport and pore structure*. San Diego, NY: Academic Press (2012).

2. Sahimi M. *Flow and transport in porous media and fractured rock: from classical methods to modern approaches*. Weinheim, Germany: John Wiley & Sons (2011).

3. Khaled ARA, Vafai K. The role of porous media in modeling flow and heat transfer in biological tissues. *Int J Heat Mass Tran* (2003) 46:4989–5003. doi:10.1016/s0017-9310(03)00301-6

7. Andrade JS, Street DA, Shinohara T, Shibusa Y, Arai Y. Percolation disorder in viscous and nonviscous flow through porous media. *Phys Rev E* (1995) 51:5725–31.

8. Seybold HJ, Carmona HA, Leonardo Filho FA, Araújo AD, Nepomuceno Filho F, Andrade JS. Flow through three-dimensional self-affine fractures. *Phys Rev Fluids* (2020) 5:104101. doi:10.1103/physrevfluids.5.104101

9. Bruderer-Weng C, Cowie P, Bernabé Y, Main I. Relating flow channelling to tracer dispersion in heterogeneous networks. *Adv Water Resour* (2004) 27:843–55. doi:10.1016/j.advwatres.2004.05.001

10. Andrade JS, Almeida MP, Mendes Filho J, Havlin S, Suki B, Stanley HE. Fluid flow through porous media: the role of stagnant zones. *Phys Rev Lett* (1997) 79:3901.

11. Tsang CF, Neretnieks I. Flow channeling in heterogeneous fractured rocks. *Rev Geophys* (1998) 36:275–98.

12. Wang G, Johannessen E, Kleijn CR, de Leeuw SW, Coppens MO. Optimizing transport in nanostructured catalysts: a computational study. *Chem Eng Sci* (2007) 62:5110–6. doi:10.1016/j.ces.2007.01.046

13. Davis ME. Ordered porous materials for emerging applications. *Nature* (2002) 417:813–21. doi:10.1038/nature00785

14. Billen J, Desmet G. Understanding and design of existing and future chromatographic support formats. *J Chromatogr A* (2007) 1168:73–2. doi:10.1016/j.chroma.2007.07.069

15. Tennikov MB, Gazdina NV, Tennikova TB, Svec F. Effect of porous structure of macroporous polymer supports on resolution in high-performance membrane chromatography of proteins. *J Chromatogr A* (1998) 798:55–64. doi:10.1016/s0021-9673(97)00873-x

16. Rubin J. 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. doi:10.1029/wr019i005p01231

17. Keil FJ. Diffusion and reaction in porous networks. *Catal Today* (1999) 53:245–58. doi:10.1016/s0920-5861(99)00119-4

19. Brimble KS, McFarlane A, Winegard N, Crowther M, Churchill DN. Effect of chronic kidney disease on red blood cell rheology. *Clin Hemorheol Microcirc* (2006) 34:411–20.

21. Xie C, Lv W, Wang M. Shear-thinning or shear-thickening fluid for better eor?–a direct pore-scale study. *J Petrol Sci Eng* (2018) 161:683–91. doi:10.1016/j.petrol.2017.11.049

23. Lai SK, Wang YY, Wirtz D, Hanes J. Micro- and macrorheology of mucus. *Adv Drug Deliv Rev* (2009) 61:86–100. doi:10.1016/j.addr.2008.09.012

24. Sandvik EI, Maerker JM. Application of xanthan gum for enhanced oil recovery. *Am Chem Soc* (1977) 45:242–64. doi:10.1021/bk-1977-0045.ch019

25. López OV, Castillo LA, Ninago MD, Ciolino AE, Villar MA. Modified starches used as additives in enhanced oil recovery (EOR) In: SN Goyanes, and NB D’Accorso, editors *Industrial applications of renewable biomass products: past, present and future*. Cham: Springer International Publishing (2017).

26. Reuvers AJ. Control of rheology of water-borne paints using associative thickeners. *Prog Org Coating* (1999) 35:171–81.

27. Wang FX. Flow analysis and modeling of field-controllable, electro-and magneto-rheological fluid dampers. *J Appl Mech* (2007) 74:13–22. doi:10.1115/1.2166649

28. Majumdar A, Butola SB, Srivastava A. Optimal designing of soft body armour materials using shear thickening fluid. *Mater Des* (2013) 46:191–8. doi:10.1016/j.matdes.2012.10.018

29. Sahimi M. Nonlinear transport processes in disordered media. *AIChE J* (1993) 39:369–86. doi:10.1002/aic.690390302

30. Shah CB, Yortsos YC. Aspects of flow of power-law fluids in porous media. *AIChE J* (1995) 41:1099–112. doi:10.1002/aic.690410506

31. De S, Kuipers JAM, Peters EAJF, Padding JT. Viscoelastic flow past mono- and bidisperse random arrays of cylinders: flow resistance, topology and normal stress distribution. *Soft Matter* (2017) 13:9138–46. doi:10.1039/c7sm01818e

32. Mokhtari Z, Zippelius A. Dynamics of active filaments in porous media. *Phys Rev Lett* (2019) 123:028001. doi:10.1103/physrevlett.123.028001

33. Cannella WJ, Huh C, Seright RS. 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). doi:10.2118/18089-MS

34. Tsakiroglou CD. 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. doi:10.1016/s0377-0257(02)00066-6

35. Sochi T, Blunt MJ. Pore-scale network modeling of ellis and herschel–bulkley fluids. *J Petrol Sci Eng* (2008) 60:105–24. doi:10.1016/j.petrol.2007.05.009

36. Berg S, van Wunnik J. Shear rate determination from pore-scale flow fields. *Transp Porous Med* (2017) 117:229–46. doi:10.1007/s11242-017-0830-3

37. Eberhard U, Seybold HJ, Floriancic M, Bertsch P, Jiménez-Martínez J, Andrade JS, et al. Determination of the effective viscosity of non-Newtonian fluids flowing through porous media. *Front Phys* (2019) 7:71. doi:10.3389/fphy.2019.00071

38. Morais AF, Seybold H, Herrmann HJ, Andrade JS. Non-Newtonian fluid flow through three-dimensional disordered porous media. *Phys Rev Lett* (2009) 103:194502. doi:10.1103/PhysRevLett.103.194502

39. Perrin CL, Tardy PM, Sorbie KS, Crawshaw JC. Experimental and modeling study of Newtonian and non-Newtonian fluid flow in pore network micromodels. *J Colloid Interface Sci* (2006) 295:542–50. doi:10.1016/j.jcis.2005.09.012

40. de Castro AR, Radilla G. 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. doi:10.1016/j.advwatres.2016.12.009

41. Hopkins CC, Haward SJ, Shen AQ. Tristability in viscoelastic flow past side-by-side microcylinders. arXiv:2010 14749v1 (2020).

42. Eberhard U, Seybold HJ, Secchi E, Jiménez-Martínez J, Rhüs P, Ofner A, et al. Mapping the local viscosity of heterogeneous non- Newtonian flows. *Sci Rep* (2020). 10:11733. doi:10.1038/s41598-020-68545-7

44. Astarita G, Marrucci G. *Principles of non-Newtonian fluid mechanics*. London, New York: McGraw-Hill Companies (1974).

45. Royer JR, Blair DL, Hudson SD. Rheological signature of frictional interactions in shear thickening suspensions. *Phys Rev Lett* (2016) 116:188301. doi:10.1103/PhysRevLett.116.188301

46. Chhabra RP. *Bubbles, drops, and particles in non-Newtonian fluids*. Boca Raton: CRC Press (2006).

48. Coussot P, Meunier M. Recognition, classification and mechanical description of debris flows. *Earth Sci Rev* (1996) 40:209–27. doi:10.1016/0012-8252(95)00065-8

49. Buzzaccaro S, Secchi E, Piazza R. Ghost particle velocimetry: accurate 3d flow visualization using standard lab equipment. *Phys Rev Lett* (2013) 111, 048101. doi:10.1103/PhysRevLett.111.048101

50. Katzbauer B. Properties and applications of xanthan gum. *Polym Degrad Stabil* (1998) 59:81–4. doi:10.1016/s0141-3910(97)00180-8

51. Walkama DM, Waisbord N, Guasto JS. Disorder suppresses chaos in viscoelastic flows. *Phys Rev Lett* (2020) 124:164501. doi:10.1103/PhysRevLett.124.164501

52. Bewersdorff HW, Singh RP. Rheological and drag reduction characteristics of xanthan gum solutions. *Rheol Acta* (1988) 27:617–27. doi:10.1007/bf01337457

54. Haase AS, Wood JA, Sprakel LM, Lammertink RG. Inelastic non-Newtonian flow over heterogeneously slippery surfaces. *Phys Rev E* (2017) 95:023105. doi:10.1103/PhysRevE.95.023105

55. Russ S, Sapoval B. Anomalous viscous damping of vibrations of fractal percolation clusters. *Phys Rev Lett* (1994) 73:1570. doi:10.1103/PhysRevLett.73.1570

Keywords: localization, microfluidics, particle velocimetry, non-newtonian fluids, porous media

Citation: Seybold HJ, Eberhard U, Secchi E, Cisne RLC, Jiménez-Martínez J, Andrade RFS, Araújo AD, Holzner M and Andrade J (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.

Edited by:

Ferenc Kun, University of Debrecen, HungaryReviewed by:

Taotao Fu, Tianjin University, ChinaSimon Haward, Okinawa Institute of Science and Technology Graduate University, Japan

Sandro Longo, University of Parma, Italy

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