Flow and axial dispersion in a sinusoidal-walled tube: Effects of inertial and unsteady flows

Abstract In this work, we consider a sinusoidal-walled tube (a three-dimensional tube with sinusoidally-varying diameter) as a simplified conceptualization of flow in porous media. Direct numerical simulation using computational fluid dynamics (CFD) methods was used to compute velocity fields by solving the Navier–Stokes equations, and also to numerically solve the volume averaging closure problem, for a range of Reynolds numbers (Re) spanning the low-Re to inertial flow regimes, including one simulation at Re = 449 for which unsteady flow was observed. The longitudinal dispersion observed for the flow was computed using a random walk particle tracking method, and this was compared to the longitudinal dispersion predicted from a volume-averaged macroscopic mass balance using the method of volume averaging; the results of the two methods were consistent. Our results are compared to experimental measurements of dispersion in porous media and to previous theoretical results for both the low-Re, Stokes flow regime and for values of Re representing the steady inertial regime. In the steady inertial regime, a power-law increase in the effective longitudinal dispersion ( D L ) with Re was found, and this is consistent with previous results. This rapid rate of increase is caused by trapping of solute in expansions due to flow separation (eddies). One unsteady (but non-turbulent) flow case ( Re = 449 ) was also examined. For this case, the rate of increase of D L with Re was smaller than that observed at lower Re. Velocity fluctuations in this regime lead to increased rates of solute mass transfer between the core flow and separated flow regions, thus diminishing the amount of tailing caused by solute trapping in eddies and thereby reducing longitudinal dispersion. The observed tailing was further explored through analysis of concentration skewness (third moment) and its assymptotic convergence to conventional advection–dispersion behavior (skewness = 0). The method of volume averaging was applied to develop a skewness model, and demonstrated that the skewness decreases as a function of inverse square root of time. Our particle tracking simulation results were shown to conform to this theoretical result in most of the cases considered.

[1]  C. Wieselsberger.  New Data on the Laws of Fluid Resistance , 1922 .

[2]  Howard Brenner,et al.  Dispersion and reaction in two‐dimensional model porous media , 1993 .

[3]  Evangelos Keramaris,et al.  Turbulent flow over and within a porous bed , 2003 .

[4]  Message P Forum,et al.  MPI: A Message-Passing Interface Standard , 1994 .

[5]  L. Beghin,et al.  On the Solutions of Linear Odd-Order Heat-Type Equations with Random Initial Conditions , 2007 .

[6]  D. Bolster,et al.  The impact of inertial effects on solute dispersion in a channel with periodically varying aperture , 2012 .

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

[8]  Francisco J. Valdés-Parada,et al.  Volume averaging: Local and nonlocal closures using a Green’s function approach , 2013 .

[9]  Adam Lambert.  Advection-diffusion in inertial flow through periodically converging-diverging tubes , 2010 .

[10]  Stephen Whitaker,et al.  Dispersion in pulsed systems—II: Theoretical developments for passive dispersion in porous media , 1983 .

[11]  V. C. Patel,et al.  Large eddy simulation of turbulent open-channel flow with free surface simulated by level set method , 2005 .

[12]  John H. Cushman,et al.  A primer on upscaling tools for porous media , 2002 .

[13]  G. Taylor.  Dispersion of soluble matter in solvent flowing slowly through a tube , 1953, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[14]  G. Bretschko,et al.  Interstitial flow through preferential flow paths in the hyporheic zone of the Oberer Seebach, Austria , 2002, Aquatic Sciences.

[15]  Martin J. Blunt,et al.  Pore‐scale modeling and continuous time random walk analysis of dispersion in porous media , 2006 .

[16]  Chiang C. Mei,et al.  Some Applications of the Homogenization Theory , 1996 .

[17]  R. Maier,et al.  Hydrodynamic dispersion in confined packed beds , 2003 .

[18]  D. Tritton.  Experiments on the flow past a circular cylinder at low Reynolds numbers , 1959, Journal of Fluid Mechanics.

[19]  C. W. Gardiner,et al.  Handbook of stochastic methods - for physics, chemistry and the natural sciences, Second Edition , 1986, Springer series in synergetics.

[20]  B. Fornberg.  A numerical study of steady viscous flow past a circular cylinder , 1980, Journal of Fluid Mechanics.

[22]  A. Straatman,et al.  Closure of non-equilibrium volume-averaged energy equations in high-conductivity porous media , 2011 .

[23]  Peter K. Kitanidis,et al.  Pore‐scale dilution of conservative solutes: An example , 1998 .

[24]  Brian D. Wood,et al.  Inertial effects in dispersion in porous media , 2007 .

[25]  S. Vanka,et al.  Convective heat transfer in periodic wavy passages , 1995 .

[26]  John L. Wilson,et al.  The influence of ambient groundwater discharge on exchange zones induced by current-bedform interactions , 2006 .

[27]  M. Cardenas.  Direct simulation of pore level Fickian dispersion scale for transport through dense cubic packed spheres with vortices , 2009 .

[28]  W. Shyy,et al.  Regular Article: An Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries , 1999 .

[29]  M. Darwish,et al.  TVD schemes for unstructured grids , 2003 .

[30]  Martin J. Blunt,et al.  Pore‐scale modeling of longitudinal dispersion , 2004 .

[31]  H. Fischer.  Mixing in Inland and Coastal Waters , 1979 .

[32]  R. Aris.  On the dispersion of a solute in a fluid flowing through a tube , 1956, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[33]  J. Delgado.  A critical review of dispersion in packed beds , 2006 .

[34]  S. Vanka,et al.  Numerical Study of Developing Flow and Heat Transfer in a Wavy Passage , 1999 .

[35]  M. Darwish,et al.  NORMALIZED VARIABLE AND SPACE FORMULATION METHODOLOGY FOR HIGH-RESOLUTION SCHEMES , 1994 .

[36]  S. Balachandar,et al.  Effect of three‐dimensionality on the lift and drag of nominally two‐dimensional cylinders , 1995 .

[37]  Jarek Nieplocha,et al.  Advances, Applications and Performance of the Global Arrays Shared Memory Programming Toolkit , 2006, Int. J. High Perform. Comput. Appl..

[38]  John B. McLaughlin,et al.  Direct Numerical Simulation of a Fully Developed Turbulent Flow over a Wavy Wall , 1998 .

[39]  S. Dennis,et al.  Numerical solutions for steady flow past a circular cylinder at Reynolds numbers up to 100 , 1970, Journal of Fluid Mechanics.

[40]  Richard A. Ketcham,et al.  Effects of inertia and directionality on flow and transport in a rough asymmetric fracture , 2009 .

[41]  Tanguy Le Borgne,et al.  Solute dispersion in channels with periodically varying apertures , 2009 .

[42]  M. Ralph.  Steady Flow Structures and Pressure Drops in Wavy-Walled Tubes , 1987 .

[43]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[44]  Steven B. Yabusaki,et al.  Building conceptual models of field‐scale uranium reactive transport in a dynamic vadose zone‐aquifer‐river system , 2008 .

[45]  V. Nassehi,et al.  Modeling of contaminants mobility in underground domains with multiple free/porous interfaces , 2003 .

[46]  Michel Quintard,et al.  Volume averaging for determining the effective dispersion tensor: Closure using periodic unit cells and comparison with ensemble averaging , 2003 .

[47]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[48]  J. Ultman,et al.  Axial dispersion through tube constrictions , 1980 .

[49]  J J Fredberg,et al.  Axial dispersion in respiratory bronchioles and alveolar ducts. , 1988, Journal of applied physiology.

[50]  M. Bayani Cardenas,et al.  Three‐dimensional vortices in single pores and their effects on transport , 2008 .

[51]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[52]  Helio Pedro Amaral Souto,et al.  Dispersion in two-dimensional periodic porous media. Part II. Dispersion tensor , 1997 .

[53]  D. S. Henn,et al.  Large-eddy simulation of flow over wavy surfaces , 1999, Journal of Fluid Mechanics.

[54]  C. Williamson.  Vortex Dynamics in the Cylinder Wake , 1996 .

[55]  Robert J. Harrison,et al.  Global arrays: A nonuniform memory access programming model for high-performance computers , 1996, The Journal of Supercomputing.

[56]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[57]  Robert W. Zimmerman,et al.  Creeping flow through a pipe of varying radius , 2001 .