Efficient projection kernels for discontinuous Galerkin simulations of disperse multiphase flows on arbitrary curved elements

Abstract In this work, we develop projection kernels for Euler-Lagrange simulations of particle-laden flows on arbitrary curved elements. These kernels are employed under a high-order discontinuous Galerkin framework for projecting the action of the particles to the Eulerian mesh. Instead of commonly used isotropic kernels, such as a Gaussian-type kernel, we construct an anisotropic polynomial-based smoothing function that preserves the compactness of the discontinuous Galerkin method on high-aspect-ratio elements and maintains an acceptable computational cost. At the same time, it mitigates inaccuracies associated with larger particles and numerical instabilities arising from the Dirac delta function and low-order kernels. Specifically, the geometric mapping of a physical element to the reference element is exploited to construct a kernel that is elliptical in 2D and ellipsoidal in 3D. We also develop a strategy to conserve interphase transfer near boundaries, particularly on curved elements. This strategy employs a high-order polynomial approximation to appropriately rescale the interphase source terms in an efficient manner. The compatibility of the proposed methodology with different types of meshes is investigated. We then apply it to a number of particle-laden flow configurations, including supersonic dusty flow over a flat plate, moving shocks interacting with clouds of particles, and hypersonic dusty flows over blunt bodies.

[1]  S. Balachandar,et al.  Dynamics of rapidly depressurized multiphase shock tubes , 2019, Journal of Fluid Mechanics.

[2]  Wai-Sun Don,et al.  A high-order WENO-Z finite difference based particle-source-in-cell method for computation of particle-laden flows with shocks , 2009, J. Comput. Phys..

[3]  Rodolfo Bermejo,et al.  A Generalized Particle Search-Locate Algorithm for Arbitrary Grids , 1997 .

[4]  Per-Olof Persson,et al.  High-order, linearly stable, partitioned solvers for general multiphysics problems based on implicit–explicit Runge–Kutta schemes , 2018, Computer Methods in Applied Mechanics and Engineering.

[5]  Clayton T. Crowe,et al.  Review—Numerical Models for Dilute Gas-Particle Flows , 1982 .

[6]  C. Munz,et al.  A Particle-in-Cell solver based on a high-order hybridizable discontinuous Galerkin spectral element method on unstructured curved meshes , 2019, Computer Methods in Applied Mechanics and Engineering.

[7]  王柏懿 Compressible laminar boundary-layer flows of a dusty gas over a semi-infinite flat-plate , 1988 .

[8]  S. Elghobashi,et al.  On the two-way interaction between homogeneous turbulence and dispersed solid particles. I: Turbulence modification , 1993 .

[9]  A. Osiptsov,et al.  Experimental and numerical study of heat transfer on a blunt body in dusty hypersonic flow , 1999 .

[10]  C. B. Henderson,et al.  Drag Coefficients of Spheres in Continuum and Rarefied Flows , 1976 .

[11]  M. Ihme,et al.  Sensitivity study of high-speed dusty flows over blunt bodies simulated using a discontinuous Galerkin method , 2019, AIAA Scitech 2019 Forum.

[12]  Matthias Ihme,et al.  A two-way coupled Euler-Lagrange method for simulating multiphase flows with discontinuous Galerkin schemes on arbitrary curved elements , 2020, J. Comput. Phys..

[13]  Krzysztof J. Fidkowski,et al.  A local sampling approach to anisotropic metric-based mesh optimization , 2016 .

[14]  S. Rebay,et al.  GMRES Discontinuous Galerkin Solution of the Compressible Navier-Stokes Equations , 2000 .

[15]  D. A. Dunavant High degree efficient symmetrical Gaussian quadrature rules for the triangle , 1985 .

[16]  Jan S. Hesthaven,et al.  High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids , 2006, J. Comput. Phys..

[17]  Xavier Pennec,et al.  A Riemannian Framework for Tensor Computing , 2005, International Journal of Computer Vision.

[18]  K. Squires,et al.  Particle response and turbulence modification in isotropic turbulence , 1990 .

[19]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[20]  E. V. Vorozhtsov,et al.  Interaction of a shock wave with a particle cloud of finite size , 2006 .

[21]  Farzad Mashayek,et al.  Spectral-based Simulations of Particle-laden Turbulent Flows ⋆ , 2009 .

[22]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[23]  Olivier Desjardins,et al.  An Euler-Lagrange strategy for simulating particle-laden flows , 2013, J. Comput. Phys..

[24]  A. Osiptsov,et al.  Heat Exchange on the Front Surface of a Blunt Body in a High‐Speed Flow Containing Low‐Inertia Particles , 2001 .

[25]  J. Nicholls,et al.  Shock wave ignition of magnesium powders , 1978 .

[26]  A. L. Stasenko,et al.  Velocity recovery factors of a particle repelled from a solid surface , 2007 .

[27]  A. Stock,et al.  Three-Dimensional Numerical Simulation of a 30-GHz Gyrotron Resonator With an Explicit High-Order Discontinuous-Galerkin-Based Parallel Particle-In-Cell Method , 2012, IEEE Transactions on Plasma Science.

[28]  David Zwick,et al.  A scalable Euler–Lagrange approach for multiphase flow simulation on spectral elements , 2020, Int. J. High Perform. Comput. Appl..

[29]  M. W. Chase NIST-JANAF thermochemical tables , 1998 .

[30]  Scott M. Murman,et al.  An Adjoint Elasticity Solver for High-Order Mesh Deformation , 2021 .

[31]  V. Fomin,et al.  Shock wave interaction with a cloud of particles , 1997 .

[32]  Peter A. Gnoffo,et al.  Shock capturing for discontinuous Galerkin methods with application to predicting heat transfer in hypersonic flows , 2019, J. Comput. Phys..

[33]  Sanjay Ranka,et al.  Dynamic Load Balancing for Compressible Multiphase Turbulence , 2018, ICS.

[34]  S. Balachandar,et al.  Pairwise-interaction extended point-particle model for particle-laden flows , 2017, J. Comput. Phys..

[35]  Rodney O. Fox,et al.  A volume-filtered description of compressible particle-laden flows , 2020 .