BRIEF RESEARCH REPORT article
Localization in Flow of Non-Newtonian Fluids Through Disordered Porous Media
- 1Institute of Terrestrial Ecosystems, Department of Environmental Systems Science, ETH Zürich, Zurich, Switzerland
- 2Departamento de Física, Universidade Federal do Ceará, Campus do Pici, Fortaleza, Brazil
- 3Institute of Environmental Engineering, Department of Civil, Environmental and Geomatic Engineering, ETH Zürich, Zürich, Switzerland
- 4Instituto Federal do Ceará, Campus Tianguá, Ceará, Brasil
- 5EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Dübendorf, Switzerland
- 6Instituto de Física, Universidade Federal da Bahia, Salvador, Brazil
- 7Swiss 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.
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  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 , shear-thickening fluids being used as active dampers  and hybrid fluids as components of enhanced body armors . 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 .
However, quantitative experiments with non-Newtonian materials which go beyond simple bulk measurements [39, 40] are scarce  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  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 ϱ,
Here the apparent viscosity
Our analysis is based on experimental results from the setup presented recently in Eberhard et al. . 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) . 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  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 , 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 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 . In the horizontal plane, we first created an unstructured quadrilateral mesh with an average cell size of
FIGURE 2. Comparison of the experimentally measured  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 qin = 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 qin 0.05 µL/min and qin 5 µL/min are labeled with yellow and blue stars, respectively.
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
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
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.
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.
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.
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
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
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).
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
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
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
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
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
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, 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
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, email@example.com