An octree-based immersogeometric approach for modeling inertial migration of particles in channels

Abstract In this paper, we develop a scalable, adaptively refined, octree-based finite element approach with immersogeometric analysis to track inertial migration of particles in microchannels. Fluid physics is modeled using a residual-based variational multiscale method, and the particle movement is modeled as rigid body motion. A parallel, hierarchically refined octree mesh is employed as the background mesh, on which a variationally consistent immersogeometric formulation is adopted for tracking the particle motion in the fluid. Adaptations of immersogeometric analysis on an octree background mesh are developed to enable efficient searching of background element for a given surface quadrature point, as well as a distribution of surface quadrature points over processors to reduce memory overhead and better parallelize the surface assembly. An octree-based adaptive mesh refinement algorithm adapted to in-out test in the immersogeometric approach is also developed. The validation of our octree-based immersogeometric approach is carried out using a benchmark case of a sphere settling in quiescent fluid, with good agreement presented. In addition, good strong (and weak) scalability on supercomputing resources for this benchmark case up to 16,384 processes is demonstrated. The proposed method is further deployed for exploring particle migration in straight and converging-diverging channels. This example illustrates the potential of the octree-based immersogeometric approach for efficiently tracking particle motion in complex channel flows – a problem with a diverse array of applications.

[1]  Ming-Chen Hsu,et al.  An immersogeometric formulation for free-surface flows with application to marine engineering problems , 2020 .

[2]  Cedric Thieulot,et al.  Adaptive octree-based finite element analysis of two- and three-dimensional indentation problems , 2008 .

[3]  Marco S. Pigazzini,et al.  Optimizing fluid–structure interaction systems with immersogeometric analysis and surrogate modeling: Application to a hydraulic arresting gear , 2017 .

[4]  Lexing Ying,et al.  A New Parallel Kernel-Independent Fast Multipole Method , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[5]  Ming-Chen Hsu,et al.  Immersogeometric analysis of compressible flows with application to aerodynamic simulation of rotorcraft , 2019, Mathematical Models and Methods in Applied Sciences.

[6]  J. D. Teresco,et al.  Parallel structures and dynamic load balancing for adaptive finite element computation , 1998 .

[7]  T. Hughes,et al.  Isogeometric variational multiscale modeling of wall-bounded turbulent flows with weakly enforced boundary conditions on unstretched meshes , 2010 .

[8]  O. Ghattas,et al.  Parallel Octree-Based Finite Element Method for Large-Scale Earthquake Ground Motion Simulation , 2005 .

[9]  Victor M. Calo,et al.  Weak Dirichlet Boundary Conditions for Wall-Bounded Turbulent Flows , 2007 .

[10]  Kumar Saurabh,et al.  Solving PDEs in space-time: 4D tree-based adaptivity, mesh-free and matrix-free approaches , 2019, SC.

[11]  Adarsh Krishnamurthy,et al.  Direct immersogeometric fluid flow analysis using B-rep CAD models , 2016, Comput. Aided Geom. Des..

[12]  Xingshi Wang,et al.  Modified Immersed Finite Element Method For Fully-Coupled Fluid-Structure Interations. , 2013, Computer methods in applied mechanics and engineering.

[13]  Baskar Ganapathysubramanian,et al.  Immersogeometric analysis of moving objects in incompressible flows , 2019, Computers & Fluids.

[14]  Barry F. Smith,et al.  PETSc Users Manual , 2019 .

[15]  Wing Kam Liu,et al.  Mathematical foundations of the immersed finite element method , 2006 .

[16]  Hari Sundar,et al.  Machine and Application Aware Partitioning for Adaptive Mesh Refinement Applications , 2017, HPDC.

[17]  Hari Sundar,et al.  Parallel geometric-algebraic multigrid on unstructured forests of octrees , 2012, 2012 International Conference for High Performance Computing, Networking, Storage and Analysis.

[18]  Thomas J. R. Hughes,et al.  Weak imposition of Dirichlet boundary conditions in fluid mechanics , 2007 .

[19]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[20]  David Eppstein,et al.  Parallel Construction of Quadtrees and Quality Triangulations , 1993, Int. J. Comput. Geom. Appl..

[21]  Adarsh Krishnamurthy,et al.  Rapid B-rep model preprocessing for immersogeometric analysis using analytic surfaces , 2017, Comput. Aided Geom. Des..

[22]  Roland Becker,et al.  Multigrid techniques for finite elements on locally refined meshes , 2000 .

[23]  R. Glowinski,et al.  Distributed Lagrange multipliers based on fictitious domain method for second order elliptic problems , 2007 .

[24]  Ming-Chen Hsu,et al.  A residual-based variational multiscale method with weak imposition of boundary conditions for buoyancy-driven flows , 2019, Computer Methods in Applied Mechanics and Engineering.

[25]  Hari Sundar,et al.  Bottom-Up Construction and 2: 1 Balance Refinement of Linear Octrees in Parallel , 2008, SIAM J. Sci. Comput..

[26]  David R. O'Hallaron,et al.  Scalable Parallel Octree Meshing for TeraScale Applications , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[27]  Yuri Bazilevs,et al.  Aerodynamic and FSI Analysis of Wind Turbines with the ALE-VMS and ST-VMS Methods , 2014 .

[28]  C. Bona-Casas,et al.  A NURBS-based immersed methodology for fluid–structure interaction , 2015 .

[29]  Dinggang Shen,et al.  Robust Computation of Mutual Information Using Spatially Adaptive Meshes , 2007, MICCAI.

[30]  Charles S. Peskin,et al.  Flow patterns around heart valves: a digital computer method for solving the equations of motion , 1973 .

[31]  Lucy T. Zhang,et al.  On computational issues of immersed finite element methods , 2009, J. Comput. Phys..

[32]  A. Oberai,et al.  The residual-based variational multiscale formulation for the large eddy simulation of compressible flows , 2012 .

[33]  D. Baraff.  An Introduction to Physically Based Modeling: Rigid Body Simulation I—Unconstrained Rigid Body Dynamics , 1997 .

[34]  Carsten Burstedde,et al.  p4est: Scalable Algorithms for Parallel Adaptive Mesh Refinement on Forests of Octrees , 2011, SIAM J. Sci. Comput..

[35]  Xiaowei Deng,et al.  Numerical Simulations of Two Back-To-Back Horizontal Axis Tidal Stream Turbines in Free-Surface Flows , 2020 .

[36]  H. Stone,et al.  Particle segregation and dynamics in confined flows. , 2009, Physical review letters.

[37]  A. K. Patra,et al.  Data structures and load balancing for parallel adaptive hp finite-element methods☆ , 2003 .

[38]  Dino Di Carlo,et al.  Nonlinear Microfluidics. , 2018, Analytical chemistry.

[39]  Lucy T. Zhang,et al.  Immersed finite element method , 2004 .

[40]  Ming-Chen Hsu,et al.  The tetrahedral finite cell method for fluids: Immersogeometric analysis of turbulent flow around complex geometries , 2016 .

[41]  Michael S. Warren,et al.  A parallel hashed oct-tree N-body algorithm , 1993, Supercomputing '93. Proceedings.

[42]  J. Nitsche.  Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind , 1971 .

[43]  Tayfun E. Tezduyar,et al.  Multiscale space-time methods for thermo-fluid analysis of a ground vehicle and its tires , 2015 .

[44]  Yuri Bazilevs,et al.  An immersogeometric variational framework for fluid-structure interaction: application to bioprosthetic heart valves. , 2015, Computer methods in applied mechanics and engineering.

[45]  Hari Sundar,et al.  A scalable framework for adaptive computational general relativity on heterogeneous clusters , 2019, ICS.

[46]  Y. Bazilevs,et al.  Isogeometric analysis of multi-phase flows with surface tension and with application to dynamics of rising bubbles , 2019, Computers & Fluids.

[47]  Hari Sundar,et al.  Massively Parallel Simulations of Binary Black Hole Intermediate-Mass-Ratio Inspirals , 2018, SIAM J. Sci. Comput..

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

[49]  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 .

[50]  H. Amini,et al.  Label-free cell separation and sorting in microfluidic systems , 2010, Analytical and bioanalytical chemistry.

[51]  T. Hughes,et al.  Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows , 2007 .

[52]  Christos Davatzikos,et al.  Low-constant parallel algorithms for finite element simulations using linear octrees , 2007, Proceedings of the 2007 ACM/IEEE Conference on Supercomputing (SC '07).

[53]  Jos Derksen,et al.  Particle imaging velocimetry experiments and lattice-Boltzmann simulations on a single sphere settling under gravity , 2002 .

[54]  David Sondak,et al.  A new class of finite element variational multiscale turbulence models for incompressible magnetohydrodynamics , 2014, J. Comput. Phys..

[55]  Jinhui Yan,et al.  Residual-based variational multi-scale modeling for particle-laden gravity currents over flat and triangular wavy terrains , 2019, Computers & Fluids.

[56]  Christos Davatzikos,et al.  Biomechanically-Constrained 4D Estimation of Myocardial Motion , 2009, MICCAI.

[57]  Grégory Legrain,et al.  On the use of the extended finite element method with quadtree/octree meshes , 2011 .