Fully-implicit orthogonal reconstructed Discontinuous Galerkin method for fluid dynamics with phase change

A new reconstructed Discontinuous Galerkin (rDG) method, based on orthogonal basis/test functions, is developed for fluid flows on unstructured meshes. Orthogonality of basis functions is essential for enabling robust and efficient fully-implicit Newton-Krylov based time integration. The method is designed for generic partial differential equations, including transient, hyperbolic, parabolic or elliptic operators, which are attributed to many multiphysics problems. We demonstrate the method's capabilities for solving compressible fluid-solid systems (in the low Mach number limit), with phase change (melting/solidification), as motivated by applications in Additive Manufacturing (AM). We focus on the method's accuracy (in both space and time), as well as robustness and solvability of the system of linear equations involved in the linearization steps of Newton-based methods. The performance of the developed method is investigated for highly-stiff problems with melting/solidification, emphasizing the advantages from tight coupling of mass, momentum and energy conservation equations, as well as orthogonality of basis functions, which leads to better conditioning of the underlying (approximate) Jacobian matrices, and rapid convergence of the Krylov-based linear solver.

[1]  G. D. Davis.  Natural convection of air in a square cavity: A bench mark numerical solution , 1983 .

[2]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[3]  Hong Luo,et al.  A Parallel Reconstructed Discontinuous Galerkin Method for the Compressible Flows on Aritrary Grids , 2010 .

[4]  J. J. Moré,et al.  Estimation of sparse jacobian matrices and graph coloring problems , 1983 .

[5]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[6]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[7]  Hong Luo,et al.  An implicit Hermite WENO reconstruction-based discontinuous Galerkin method on tetrahedral grids , 2014 .

[8]  Robert Youngblood,et al.  Development, Selection, Implementation and Testing of Architectural Features and Solution Techniques for Next Generation of System Simulation Codes to Support the Safety Case if the LWR Life Extension , 2010 .

[9]  Rainald Löhner,et al.  A discontinuous Galerkin method based on a Taylor basis for the compressible flows on arbitrary grids , 2008, J. Comput. Phys..

[10]  Ben Q. Li.  Discontinuous Finite Elements in Fluid Dynamics and Heat Transfer , 2005 .

[11]  Rainald Löhner,et al.  A p-multigrid discontinuous Galerkin method for the Euler equations on unstructured grids , 2006 .

[12]  Michael Dumbser,et al.  Efficient implementation of high order unstructured WENO schemes for cavitating flows , 2013 .

[13]  R. Aris.  Vectors, Tensors and the Basic Equations of Fluid Mechanics , 1962 .

[14]  H. J.,et al.  Hydrodynamics , 1924, Nature.

[15]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[16]  L. M. Albright.  Vectors , 2003, Current protocols in molecular biology.

[17]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[18]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[19]  Hong Luo,et al.  A Reconstructed Discontinuous Galerkin Method Based on a Hierarchical Hermite WENO Reconstruction for Compressible Flows on Tetrahedral Grids , 2012 .

[20]  Michael Dumbser,et al.  Very high order PNPM schemes on unstructured meshes for the resistive relativistic MHD equations , 2009, J. Comput. Phys..

[21]  S. Khairallah,et al.  Mesoscopic Simulation Model of Selective Laser Melting of Stainless Steel Powder , 2014 .

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

[23]  Y. Mukaigawa,et al.  Large Deviations Estimates for Some Non-local Equations I. Fast Decaying Kernels and Explicit Bounds , 2022 .

[24]  Vicente Hernández,et al.  SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems , 2005, TOMS.

[25]  V. Voller,et al.  A fixed grid numerical modelling methodology for convection-diffusion mushy region phase-change problems , 1987 .

[26]  Robert R. Nourgaliev.  Modeling and analysis of heat and mass transfer processes during in-vessel melt progression stage of light water reactor (LWR) severe accidents , 1998 .

[27]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[28]  Bram Van Leer,et al.  Discontinuous Galerkin for Diffusion , 2005 .

[29]  Hong Luo,et al.  A third‐order implicit discontinuous Galerkin method based on a Hermite WENO reconstruction for time‐accurate solution of the compressible Navier–Stokes equations , 2015 .

[30]  Chia-Jung Hsu.  Numerical Heat Transfer and Fluid Flow , 1981 .

[31]  Theo G. Theofanous,et al.  Direct Numerical Simulation of Interfacial Flows: Implicit Sharp-Interface Method (I-SIM) , 2008 .

[32]  D. Kershaw.  Differencing of the diffusion equation in Lagrangian hydrodynamic codes , 1981 .

[33]  Hong Luo,et al.  A set of parallel, implicit methods for a reconstructed discontinuous Galerkin method for compressible flows on 3D hybrid grids , 2014 .

[34]  Hong Luo,et al.  A Hermite WENO reconstruction-based discontinuous Galerkin method for the Euler equations on tetrahedral grids , 2012, J. Comput. Phys..

[35]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[36]  A. Mccarthy.  Development , 1996, Current Opinion in Neurobiology.

[37]  A. A. Amsden,et al.  A numerical fluid dynamics calculation method for all flow speeds , 1971 .

[38]  Robert Nourgaliev,et al.  Recovery Discontinuous Galerkin Jacobian-Free Newton-Krylov Method for Multiphysics Problems , 2010 .

[39]  Robert A. van de Geijn,et al.  Elemental: A New Framework for Distributed Memory Dense Matrix Computations , 2013, TOMS.

[40]  A. A. Amsden,et al.  Numerical calculation of almost incompressible flow , 1968 .

[41]  Vincent Mousseau,et al.  A reconstructed discontinuous Galerkin method for the compressible Navier-Stokes equations on arbitrary grids , 2010, J. Comput. Phys..

[42]  Hong Luo,et al.  Fully-Implicit Orthogonal Reconstructed Discontinuous Petrov-Galerkin Method for Multiphysics Problems , 2015 .

[43]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[44]  Jonathan A. Dantzig,et al.  MODELLING LIQUID-SOLID PHASE CHANGES WITH MELT CONVECTION , 1989 .

[45]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[46]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one-dimensional case , 2004 .

[47]  Michael Dumbser,et al.  Arbitrary high order PNPM schemes on unstructured meshes for the compressible Navier–Stokes equations , 2010 .

[48]  Richard C. Martineau,et al.  On physics-based preconditioning of the Navier-Stokes equations , 2009, J. Comput. Phys..

[49]  Hong Luo,et al.  A Reconstructed Discontinuous Galerkin Method for the Compressible Euler Equations on Arbitrary Grids , 2009 .

[50]  Hester Bijl,et al.  Implicit Time Integration Schemes for the Unsteady Compressible Navier–Stokes Equations: Laminar Flow , 2002 .

[51]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[52]  R. I. Issa,et al.  AN IMPROVED PISO ALGORITHM FOR THE COMPUTATION OF BUOYANCY-DRIVEN FLOWS , 2001 .

[53]  Richard D. Hornung,et al.  The RAJA Portability Layer: Overview and Status , 2014 .

[54]  Philip M. Gresho,et al.  On the theory of semi‐implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. Part 1: Theory , 1990 .

[55]  Bram van Leer,et al.  Analysis and Implementation of Recovery-Based Discontinuous Galerkin for Diffusion , 2009 .

[56]  Hong Luo,et al.  A Reconstructed Discontinuous Galerkin Method for the Euler Equations on Arbitrary Grids , 2012 .

[57]  Hester Bijl,et al.  Fourth-Order Runge–Kutta Schemes for Fluid Mechanics Applications , 2005, J. Sci. Comput..