Direct numerical simulations of viscous suspensions with variably shaped crystals

Abstract Our ability to numerically model and understand the complex flow behavior of solid-bearing suspensions has increased significantly over the last couple of years, partly due to direct numerical simulations that compute flow around individual immersed particles. Direct numerical simulations hence resolve suspension dynamics in unprecedented detail. While most previous studies focus on spherical particles, we develop a direct numerical approach to also capture rectangular particles. Our approach uses distributed Lagrange multipliers to enforce rigid-body motion in the solid domain in combination with an immersed boundary method to correctly enforce the no-slip constraint on the solid-fluid interfaces. An important component of our model is an efficient particle collision scheme that prevents overlap between particles of different shapes and allows for the transfer of both translational and angular momentum during particle collision. We verify and validate our numerical method through several benchmark cases. Applied to particulate suspensions, we test the hypothesize that particle rotations alter the aggregated dynamics of the suspension even if the relative rotational energy of the particles remains small as compared to the translational energy. At low solid fraction, we reproduce experimental observations of strongly nonlinear coupling between rectangular particles that is reminiscent of particle aggregation in the inertial regime but occurs at zero Reynolds number as a result of the long-range interaction between non-spherical particles. At intermediate solid fraction, we show that particle rotations can destabilize force chains. The dynamic consequences include the delayed onset of jamming and strong nonlinear coupling to the flow field in the fluid domain, which channelizes more strongly for rectangular as compared to spherical particles. Our model was motivated specifically by the need to better understand hazardous, crystal-bearing lava flows, but our insights generalize to viscous suspension flow in other scientific or engineering contexts.

[1]  P. Mucha,et al.  Fluid dynamics: That sinking feeling , 2001, Nature.

[2]  Renaud Toussaint,et al.  The Mechanical Coupling of Fluid-Filled Granular Material Under Shear , 2011 .

[3]  Timothy A. Davis,et al.  Direct methods for sparse linear systems , 2006, Fundamentals of algorithms.

[4]  W. Carlson,et al.  Plagioclase-chain networks in slowly cooled basaltic magma , 1999 .

[5]  Marcelo Horacio Garcia,et al.  Experiments on Saltation of Sand in Water , 1998 .

[6]  I. Orlanski.  A Simple Boundary Condition for Unbounded Hyperbolic Flows , 1976 .

[7]  A. Burgisser,et al.  On the kinematics and dynamics of crystal‐rich systems , 2017 .

[8]  D. Joseph,et al.  Nonlinear mechanics of fluidization of beds of spherical particles , 1987, Journal of Fluid Mechanics.

[9]  G. B. Jeffery.  The motion of ellipsoidal particles immersed in a viscous fluid , 1922 .

[10]  Allen T. Chwang,et al.  Hydromechanics of low-Reynolds-number flow. Part 3. Motion of a spheroidal particle in quadratic flows , 1975, Journal of Fluid Mechanics.

[11]  J. Fröhlich,et al.  Collision modelling for the interface-resolved simulation of spherical particles in viscous fluids , 2012, Journal of Fluid Mechanics.

[12]  C. Petrie,et al.  The rheology of fibre suspensions , 1999 .

[13]  S. Tulaczyk,et al.  Discrete element modeling of subglacial sediment deformation , 2013 .

[14]  Julian Simeonov,et al.  Modeling mechanical contact and lubrication in Direct Numerical Simulations of colliding particles , 2012 .

[15]  W. Breugem,et al.  Turbulent channel flow of dense suspensions of neutrally buoyant spheres , 2014, Journal of Fluid Mechanics.

[16]  A. Rempel,et al.  Gas accumulation in particle-rich suspensions and implications for bubble populations in crystal-rich magma , 2010 .

[17]  W. R. Briley,et al.  Solution of the multidimensional compressible Navier-Stokes equations by a generalized implicit method , 1977 .

[18]  Hideya Nakamura,et al.  Numerical modeling of particle fluidization behavior in a rotating fluidized bed , 2007 .

[19]  E. Balaras.  Modeling complex boundaries using an external force field on fixed Cartesian grids in large-eddy simulations , 2004 .

[20]  G. Hulme,et al.  The Interpretation of Lava Flow Morphology , 1974 .

[21]  Gregory H. Wannier,et al.  A contribution to the hydrodynamics of lubrication , 1950 .

[22]  R. Rangel,et al.  Numerical investigation of particle–particle and particle–wall collisions in a viscous fluid , 2008, Journal of Fluid Mechanics.

[23]  E. W. Llewellin,et al.  The rheology of suspensions of solid particles , 2010, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[24]  C. Chaffey,et al.  Rotation of an isolated triaxial ellipsoid suspended in slow viscous flow , 1978 .

[25]  David Baraff,et al.  Analytical methods for dynamic simulation of non-penetrating rigid bodies , 1989, SIGGRAPH.

[26]  Paul M. Chaikin,et al.  Long-range correlations in sedimentation , 1997 .

[27]  T. Y. Wu,et al.  Hydromechanics of low-Reynolds-number flow. Part 2. Singularity method for Stokes flows , 1975, Journal of Fluid Mechanics.

[28]  A. Philpotts,et al.  PHYSICAL PROPERTIES OF PARTLY MELTED THOLEIITIC BASALT , 1996 .

[29]  S. G. Mason,et al.  The kinetics of flowing dispersions: I. Concentrated suspensions of rigid particles , 1966 .

[30]  Markus Uhlmann,et al.  Interface-resolved direct numerical simulation of vertical particulate channel flow in the turbulent regime , 2008, 1108.6233.

[31]  Jochen Fröhlich,et al.  Direct Numerical Simulations of spherical bubbles in vertical turbulent channel flow , 2015 .

[32]  H. A. Becker.  The effects of shape and reynolds number on drag in the motion of a freely oriented body in an infinite fluid , 1959 .

[33]  Chaikin,et al.  Diffusion, dispersion, and settling of hard spheres. , 1992, Physical review letters.

[34]  Roger H. Rangel,et al.  Collision of multi-particle and general shape objects in a viscous fluid , 2008, J. Comput. Phys..

[35]  G. G. Stokes.  On the Effect of the Internal Friction of Fluids on the Motion of Pendulums , 2009 .

[36]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[37]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[38]  G. Bergantz,et al.  On the hydrodynamics of crystal clustering , 2018, Philosophical Transactions of the Royal Society A.

[39]  K. R. Schwindinger.  Particle dynamics and aggregation of crystals in a magma chamber with application to Kilauea Iki olivines , 1999 .

[40]  Z. Qin,et al.  Direct numerical simulations of gas–solid–liquid interactions in dilute fluids , 2017 .

[41]  Jianren Fan,et al.  A modified immersed boundary method for simulations of fluid–particle interactions , 2007 .

[42]  Gene H. Golub,et al.  Matrix computations , 1983 .

[43]  E. Lajeunesse,et al.  Bed load transport in turbulent flow at the grain scale: Experiments and modeling , 2010 .

[44]  Wim-Paul Breugem,et al.  A second-order accurate immersed boundary method for fully resolved simulations of particle-laden flows , 2012, J. Comput. Phys..

[45]  R. Glowinski,et al.  A distributed Lagrange multiplier/fictitious domain method for particulate flows , 1999 .

[46]  Allen T. Chwang,et al.  Hydromechanics of low-Reynolds-number flow. Part 4. Translation of spheroids , 1976, Journal of Fluid Mechanics.

[47]  J. Westerweel,et al.  Collision model for fully resolved simulations of flows laden with finite-size particles. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  Bruce J. Ackerson,et al.  Analogies between Colloidal Sedimentation and Turbulent Convection at High Prandtl Numbers , 1998 .

[49]  George W. Bergantz,et al.  Open-system dynamics and mixing in magma mushes , 2015 .

[50]  S. Elghobashi,et al.  Modulation of isotropic turbulence by particles of Taylor length-scale size , 2010, Journal of Fluid Mechanics.

[51]  A. Burgisser,et al.  Lubrication effects on magmatic mush dynamics , 2019, Journal of Volcanology and Geothermal Research.

[52]  Alexandre J. Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations , 1969 .

[53]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[54]  H. Ishibashi.  Non-Newtonian behavior of plagioclase-bearing basaltic magma: Subliquidus viscosity measurement of the 1707 basalt of Fuji volcano, Japan , 2009 .

[55]  Pradipta Kumar Panigrahi,et al.  Experimental Investigation of Flow Past a Square Cylinder at an Angle of Incidence , 2008 .

[56]  J. Happel,et al.  Low Reynolds number hydrodynamics: with special applications to particulate media , 1973 .

[57]  R. Glowinski,et al.  A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow , 2001 .

[58]  L. G. Leal,et al.  Rotation of small non-axisymmetric particles in a simple shear flow , 1979, Journal of Fluid Mechanics.

[59]  Mathieu Martin,et al.  A numerical method for fully resolved simulation (FRS) of rigid particle-flow interactions in complex flows , 2009, J. Comput. Phys..

[60]  J. Sethian,et al.  Crystals stirred up: 1. Direct numerical simulations of crystal settling in nondilute magmatic suspensions , 2012 .

[61]  Howard Brenner,et al.  Rheology of a dilute suspension of axisymmetric Brownian particles , 1974 .

[62]  James P. Kauahikaua,et al.  Cooling and crystallization of lava in open channels, and the transition of Pāhoehoe Lava to 'A'ā , 1999 .

[63]  Allen T. Chwang,et al.  Hydromechanics of low-Reynolds-number flow. Part 1. Rotation of axisymmetric prolate bodies , 1974, Journal of Fluid Mechanics.

[64]  R. Glowinski,et al.  A new formulation of the distributed Lagrange multiplier/fictitious domain method for particulate flows , 2000 .

[65]  P. Umbanhowar,et al.  An effective gravitational temperature for sedimentation , 2001, Nature.

[66]  M. Vergassola,et al.  Particles and fields in fluid turbulence , 2001, cond-mat/0105199.

[67]  S. Pannala,et al.  Open-source MFIX-DEM software for gas-solids flows: Part I – verification studies , 2012 .

[68]  Edward Biegert,et al.  A collision model for grain-resolving simulations of flows over dense, mobile, polydisperse granular sediment beds , 2016, J. Comput. Phys..

[69]  SCREENED AND UNSCREENED PHASES IN SEDIMENTING SUSPENSIONS , 1998, cond-mat/9801164.

[70]  Pui-Kuen Yeung,et al.  LAGRANGIAN INVESTIGATIONS OF TURBULENCE , 2003 .

[71]  Jochen Fröhlich,et al.  On the relevance of collision modeling for interface-resolving simulations of sediment transport in open channel flow , 2014 .