Pore-Scale Modeling of Non-Newtonian Flow in Porous Media

The thesis investigates the flow of non-Newtonian fluids in porous media using pore-scale network modeling. Non-Newtonian fluids show very complex time and strain dependent behavior and may have initial yield stress. Their common feature is that they do not obey the simple Newtonian relation of proportionality between stress and rate of deformation. They are generally classified into three main categories: time-independent, time-dependent and viscoelastic. Two three-dimensional networks representing a sand pack and Berea sandstone were used. An iterative numerical technique is used to solve the pressure field and obtain the flow rate and apparent viscosity. The time-independent category is investigated using two fluid models: Ellis and Herschel-Bulkley. The analysis confirmed the reliability of the non-Newtonian network model used in this study. Good results are obtained, especially for the Ellis model, when comparing the network model results to experimental data sets found in the literature. The yield-stress phenomenon is also investigated and several numerical algorithms were developed and implemented to predict threshold yield pressure of the network. An extensive literature survey and investigation were carried out to understand the phenomenon of viscoelasticity with special attention to the flow in porous media. The extensional flow and viscosity and converging-diverging geometry were thoroughly examined as the basis of the peculiar viscoelastic behavior in porous media. The modified Bautista-Manero model was identified as a promising candidate for modeling the flow of viscoelastic materials which also show thixotropic attributes. An algorithm that employs this model was implemented in the non-Newtonian code and the initial results were analyzed. The time-dependent category was examined and several problems in modeling and simulating the flow of these fluids were identified.

[1]  Kenneth Stuart Sorbie,et al.  Polymer-improved oil recovery , 1991 .

[2]  C. Tiu,et al.  Flow of non-Newtonian fluids through tubes with abrupt expansions and contractions (square wave tubes) , 1997 .

[3]  W. L. Wilkinson,et al.  Non-Newtonian Flow and Heat Transfer , 1972 .

[4]  C. Petrie Dynamics of polymeric liquids. Volume 1. Fluid mechanics (2nd ed) : R.B. Bird, R.C. Armstrong and O. Hassager, Wiley-Interscience, New York, NY, 1987, 649 + xxi pages, ISBN 0-471-80245-X (V.1), price US$ 69.95 , 1988 .

[5]  Yasuki Nakayama,et al.  Flow in pipes , 1998 .

[6]  Pavel Kuzhir,et al.  Flow of magnetorheological fluid through porous media , 2003 .

[7]  J. Mewis,et al.  Thixotropy : Build-up and breakdown curves during flow , 2005 .

[8]  M. Blunt,et al.  Predicting the Impact of Non-Newtonian Rheology on Relative Permeability Using Pore-Scale Modeling , 2004 .

[9]  J. C. Godfrey Steady shear measurement of thixotropic fluid properties , 1973 .

[10]  O. Manero,et al.  Dynamics of worm-like micelles: the Cox–Merz rule , 2002 .

[11]  H. Barnes Thixotropy—a review , 1997 .

[12]  C. Pérez-García,et al.  Linear instability in viscoelastic fluid convection , 1990 .

[13]  A. Beris,et al.  Applications of Domain Decomposition Spectral Collocation Methods in Viscoelastic Flows Through Model Porous Media , 1992 .

[14]  D. V. Boger A highly elastic constant-viscosity fluid , 1977 .

[15]  G. Vradis,et al.  Macroscopic conductivities for flow of bingham plastics in porous media , 1994 .

[16]  W. B. Gogarty,et al.  Viscoelastic Effects in Polymer Flow Through Porous Media , 1972 .

[17]  A. Garrouch A Viscoelastic Model for Polymer Flow in Reservoir Rocks , 1999 .

[18]  O. Manero,et al.  On the shear banding flow of elongated micellar solutions , 2000 .

[19]  Yu-Shu Wu,et al.  Flow and Displacement of Bingham Non-Newtonian Fluids in Porous Media , 1992 .

[20]  M. Balhoff,et al.  A macroscopic model for shear-thinning flow in packed beds based on network modeling , 2006 .

[21]  Karsten E. Thompson,et al.  Modeling the steady flow of yield‐stress fluids in packed beds , 2004 .

[22]  Stig Bakke,et al.  Extending Predictive Capabilities to Network Models , 1998 .

[23]  M. Letelier,et al.  Pulsating flow of viscoelastic fluids in straight tubes of arbitrary cross-section—Part I: longitudinal field , 2002 .

[24]  F. Civan,et al.  Experimental and theoretical investigation of the Zaoyuan field heavy oil flow through porous media , 2006 .

[25]  A. B. Metzner,et al.  Thermodynamic and heat transport considerations for viscoelastic fluids , 1965 .

[26]  Timothy Nigel Phillips,et al.  Viscoelastic flow in an undulating tube using spectral methods , 2004 .

[27]  N. Phan-Thien,et al.  Flow of an oldroyd-type fluid through a sinusoidally corrugated tube , 1987 .

[28]  Matthew D. Jackson,et al.  Detailed physics, predictive capabilities and macroscopic consequences for pore-network models of multiphase flow. , 2002 .

[29]  G. B. The Dynamical Theory of Gases , 1916, Nature.

[30]  R. Larson The Structure and Rheology of Complex Fluids , 1998 .

[31]  M. Denn ISSUES IN VISCOELASTIC FLUID MECHANICS , 1990 .

[32]  W. R. Schowalter,et al.  Modeling the flow of viscoelastic fluids through porous media , 1981 .

[33]  Franz Durst,et al.  The nature of flows through porous media , 1987 .

[34]  M. Nicolet The Earth: its Origin, History, and Physical Constitution , 1929, Nature.

[35]  Raanan Fattal,et al.  Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation , 2005 .

[36]  H. Pascal Nonsteady flow through porous media in the presence of a threshold gradient , 1981 .

[37]  A. Collyer Time dependent fluids , 1974 .

[38]  Gary A. Pope,et al.  ANALYSIS OF FACTORS INFLUENCING MOBILITY AND ADSORPTION IN THE FLOW OF POLYMER SOLUTION THROUGH POROUS MEDIA. , 1974 .

[39]  H. Kharabaf,et al.  INVASION PERCOLATION WITH MEMORY , 1997 .

[40]  A. Hansen,et al.  A new algorithm to extract the backbone in a random resistor network , 1987 .

[41]  J. G. Savins NON-NEWTONIAN FLOW THROUGH POROUS MEDIA , 1969 .

[42]  Krishnamurthy Jayaraman,et al.  The effect of polymer extensibility on crossflow of polymer solutions through cylinder arrays , 1992 .

[43]  Stanley Middleman,et al.  Flow of Viscoelastic Fluids through Porous Media , 1967 .

[44]  Matthew T. Balhoff,et al.  Modeling the flow of non-Newtonian fluids in packed beds at the pore scale , 2005 .

[45]  Martin J. Blunt,et al.  Predictive pore‐scale modeling of two‐phase flow in mixed wet media , 2004 .

[46]  Marco Dressler,et al.  Computational Rheology , 2002 .

[47]  J. Billingham,et al.  Laminar, unidirectional flow of a thixotropic fluid in a circular pipe , 1993 .

[48]  F. Dullien,et al.  Single phase flow through porous media and pore structure , 1975 .

[49]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[50]  O. Manero,et al.  Understanding thixotropic and antithixotropic behavior of viscoelastic micellar solutions and liquid crystalline dispersions. I. The model , 1999 .

[51]  R. Keunings MICRO-MACRO METHODS FOR THE MULTISCALE SIMULATION OF VISCOELASTIC FLOW USING MOLECULAR MODELS OF KINETIC THEORY , 2004 .

[52]  M. Escudier,et al.  Pipe flow of a thixotropic liquid , 1996 .

[53]  P. M. J. Tardy,et al.  Models for flow of non-Newtonian and complex fluids through porous media , 2002 .

[54]  B. Khomami,et al.  Application of higher order finite element methods to viscoelastic flow in porous media , 1992 .

[55]  William R. Rossen,et al.  Percolation theory of creation and mobilization of foams in porous media , 1990 .

[56]  R. V. Craster,et al.  Geophysical Aspects of Non-Newtonian Fluid Mechanics , 2001 .

[57]  Antony N. Beris,et al.  Calculations of steady-state viscoelastic flow in an undulating tube , 1989 .

[58]  T. Papanastasiou,et al.  Viscous Fluid Flow , 1999 .

[59]  Howard A. Barnes,et al.  The yield stress—a review or ‘παντα ρει’—everything flows? , 1999 .

[60]  T. Al-Fariss Flow of polymer solutions through porous media , 1990 .

[61]  Howard A. Barnes,et al.  A review of the slip (wall depletion) of polymer solutions, emulsions and particle suspensions in viscometers: its cause, character, and cure , 1996 .

[62]  V. Chaplain,et al.  Model for the flow of a yield fluid through a porous medium , 1992 .

[63]  A. Beris,et al.  Viscoelastic flow in a periodically constricted tube: The combined effect of inertia, shear thinning, and elasticity , 1991 .

[64]  R. Larson Constitutive equations for polymer melts and solutions , 1988 .

[65]  R. Bird Dynamics of Polymeric Liquids , 1977 .

[66]  R. Keunings FINITE ELEMENT METHODS FOR INTEGRAL VISCOELASTIC FLUIDS , 2003 .

[67]  A. S. Lodge Elastic Liquids: An Introductory Vector Treatment of Finite-Strain Polymer Rheology , 1964 .

[68]  T. Al-Fariss A new correlation for non-newtonian flow through porous media , 1989 .

[69]  Y. Yortsos,et al.  The flow and displacement in porous media of fluids with yield stress , 2005 .

[70]  J. Pearson,et al.  Viscous fingering of a thixotropic fluid in a porous medium or a narrow fracture , 2006 .

[71]  M. A. Hulsen A numerical method for solving steady 2D and axisymmetrical viscoelastic flow problems with an application to inertia effects in contraction flows , 1990 .

[72]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[73]  P. Saffman,et al.  A theory of dispersion in a porous medium , 1959, Journal of Fluid Mechanics.

[74]  S. Middleman,et al.  Power-Law Flow through a Packed Tube , 1965 .

[75]  Hans Muhlhaus,et al.  A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials , 2003 .

[76]  K. Jayaraman,et al.  Elastic instability in crossflow of polymer solutions through periodic arrays of cylinders , 1993 .

[77]  Marcel Crochet,et al.  Numerical Methods in Non-Newtonian Fluid Mechanics , 1983 .

[78]  J.J. Shea,et al.  Rheology of Polymeric Systems , 1998, IEEE Electrical Insulation Magazine.

[79]  J. Schieber,et al.  Viscoelastic flow through fibrous media using the CONNFFESSIT approach , 1998 .

[80]  M. Blunt Flow in porous media — pore-network models and multiphase flow , 2001 .

[81]  D. D. Joye,et al.  Characteristics of Thixotropic Behavior , 1971 .

[82]  N. Kohler,et al.  Two-Phase Flow Through Porous Media: Effect of an Adsorbed Polymer Layer , 1988 .

[83]  J. Mewis,et al.  Thixotropy - a general review , 1979 .

[84]  N. Mori,et al.  Measurement of pressure loss and observation of the flow field in viscoelastic flow through an undulating channel , 2000 .

[85]  P. Royer,et al.  Derivation of macroscopic filtration law for transient linear viscoelastic fluid flow in porous media , 2000 .

[86]  D. V. Boger Viscoelastic Flows Through Contractions , 1987 .

[87]  George G. Chase,et al.  Incompressible Cake Filtration of a Yield Stress Fluid , 2003 .

[88]  A. Beris,et al.  Spectral collocation/domain decomposition method for viscoelastic flow simulations in model porous geometries , 1996 .

[89]  J. Bear Dynamics of Fluids in Porous Media , 1975 .

[90]  D. Cheng,et al.  Phenomenological characterization of the rheological behaviour of inelastic reversible thixotropic and antithixotropic fluids , 1965 .

[91]  R. Dharmadhikari,et al.  Flow of non-Newtonian fluids through porous media , 1985 .

[92]  J. Sherwood,et al.  Oscillation superimposed on steady shearing: Measurements and predictions for wormlike micellar solutions , 2006 .

[93]  Eugene H. Wissler,et al.  Viscoelastic Effects in the Flow of Non-Newtonian Fluids through a Porous Medium , 1971 .

[94]  Martin J Blunt,et al.  Predictive network modeling of single-phase non-Newtonian flow in porous media. , 2003, Journal of colloid and interface science.

[95]  L. G. Leal,et al.  The influence of rheological properties on the slow flow past spheres , 1987 .

[96]  P. Valvatne Predictive Pore-Scale Modelling of Multiphase Flow , 2009 .

[97]  D. Wang,et al.  Prediction of IPR Curve of Oil Wells in Visco-Elastic Polymer Solution Flooding Reservoirs , 2001 .

[98]  R. P. Chhabra,et al.  Flow of non-Newtonian fluids in fixed and fluidised beds , 2001 .

[99]  Stig Bakke,et al.  Reconstruction of Berea sandstone and pore-scale modelling of wettability effects , 2003 .

[100]  Spectral element methods for transient viscoelastic flow problems , 2004 .

[101]  Rossen,et al.  Minimal path for transport in networks. , 1993, Physical review. B, Condensed matter.

[102]  D. Piotrowski,et al.  [Introduction to rheology]. , 1982, Acta haematologica Polonica.