An internal penalty discontinuous Galerkin method for simulating conjugate heat transfer in a closed cavity

Summary Using the discontinuous Galerkin (DG) method for conjugate heat transfer problems can provide improved accuracy close to the fluid-solid interface, localizing the data exchange process, which may further enhance the convergence and stability of the entire computation. This paper presents a framework for the simulation of conjugate heat transfer problems using DG methods on unstructured grids. Based on an existing DG solver for the incompressible Navier-Stokes equation, the fluid advection-diffusion equation, Boussinesq term, and solid heat equation are introduced using an explicit DG formulation. A Dirichlet-Neumann partitioning strategy has been implemented to achieve the data exchange process via the numerical flux of interface quadrature points in the fluid-solid interface. Formal h and p convergence studies employing the method of manufactured solutions demonstrate that the expected order of accuracy is achieved. The algorithm is then further validated against 3 existing benchmark cases, including a thermally driven cavity, conjugate thermally driven cavity, and a thermally driven cavity with conducting solid, at Rayleigh numbers from 1000 to 100 000. The computational effort is documented in detail demonstrating clearly that, for all cases, the highest-order accurate algorithm has several magnitudes lower error than first- or second-order schemes for a given computational effort.

[1]  Christopher J. Roy,et al.  Review of code and solution verification procedures for computational simulation , 2005 .

[2]  B. Rivière,et al.  Esaim: Mathematical Modelling and Numerical Analysis a Splitting Method Using Discontinuous Galerkin for the Transient Incompressible Navier-stokes Equations , 2022 .

[3]  Béatrice Rivière,et al.  Discontinuous finite element methods for incompressible flows on subdomains with non-matching interfaces , 2006 .

[4]  Tanmay Basak,et al.  Heatline analysis on thermal management with conjugate natural convection in a square cavity , 2013 .

[5]  Theodore J. Heindel,et al.  CONJUGATE NATURAL CONVECTION FROM AN ARRAY OF PROTRUDING HEAT SOURCES , 1996 .

[6]  C. Gu,et al.  Application of a discontinuous Galerkin method on the compressible flow in the transonic axial compressor , 2016 .

[7]  S. Rebay,et al.  An implicit high-order discontinuous Galerkin method for steady and unsteady incompressible flows , 2007 .

[8]  Taehun Lee,et al.  A spectral-element discontinuous Galerkin lattice Boltzmann method for nearly incompressible flows , 2011, J. Comput. Phys..

[10]  Béatrice Rivière,et al.  A Discontinuous Subgrid Eddy Viscosity Method for the Time-Dependent Navier-Stokes Equations , 2005, SIAM J. Numer. Anal..

[11]  M. Wheeler.  An Elliptic Collocation-Finite Element Method with Interior Penalties , 1978 .

[12]  Carlos A. Felippa,et al.  Staggered transient analysis procedures for coupled mechanical systems: Formulation , 1980 .

[13]  Xiaodong Ren,et al.  An investigation of the conjugate heat transfer in an intercooled compressor vane based on a discontinuous Galerkin method , 2017 .

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

[15]  John M. House,et al.  EFFECT OF A CENTERED CONDUCTING BODY ON NATURAL CONVECTION HEAT TRANSFER IN AN ENCLOSURE , 1990 .

[16]  G. Chavent,et al.  A finite-element method for the 1-D water flooding problem with gravity , 1982 .

[17]  Antonio Huerta,et al.  Discontinuous Galerkin methods for the Navier–Stokes equations using solenoidal approximations , 2009 .

[18]  David E. Keyes,et al.  Smooth and robust solutions for Dirichlet boundary control of fluid-solid conjugate heat transfer problems , 2015, J. Comput. Phys..

[19]  David L. Darmofal,et al.  An adaptive simplex cut-cell method for high-order discontinuous Galerkin discretizations of elliptic interface problems and conjugate heat transfer problems , 2014, J. Comput. Phys..

[20]  Pierre Lallemand,et al.  Lattice Boltzmann simulations of thermal convective flows in two dimensions , 2013, Comput. Math. Appl..

[21]  Meng-Sing Liou,et al.  Application of Chimera/unstructured hybrid grids for conjugate heat transfer , 1997 .

[22]  J. Douglas,et al.  Interior Penalty Procedures for Elliptic and Parabolic Galerkin Methods , 1976 .

[23]  Thomas J. R. Hughes,et al.  A comparison of discontinuous and continuous Galerkin methods bases on error estimates, conservation, robustness and efficiency , 2000 .

[24]  Hester Bijl,et al.  Analysis and application of high order implicit Runge-Kutta schemes for unsteady conjugate heat transfer: A strongly-coupled approach , 2014, J. Comput. Phys..

[25]  Saumil Patel,et al.  A spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for conjugate heat transfer applications , 2016 .

[26]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[27]  A. Veeraragavan,et al.  Use of the method of manufactured solutions for the verification of conjugate heat transfer solvers , 2016, Journal of Computational Physics.

[28]  B. Rivière,et al.  Improved energy estimates for interior penalty, constrained and discontinuous Galerkin methods for elliptic problems. Part I , 1999 .

[29]  Michael B. Giles,et al.  Stability analysis of numerical interface conditions in fluid-structure thermal analysis , 1997 .