Numerical solution of Rosseland model for transient thermal radiation in non-grey optically thick media using enriched basis functions

Abstract Heat radiation in optically thick non-grey media can be well approximated with the Rosseland model which is a class of nonlinear diffusion equations with convective boundary conditions. The optical spectrum is divided into a set of finite bands with constant absorption coefficients but with variable Planckian diffusion coefficients. This simplification reduces the computational costs significantly compared to solving a full radiative heat transfer model. Therefore, the model is very popular for industrial and engineering applications. However, the opaque nature of the media often results in thermal boundary layers that requires highly refined meshes, to be recovered numerically. Such meshes can significantly hinder the performance of numerical methods. In this work we explore for the first time using enriched basis functions for the model in order to avoid using refined meshes. In particular, we discuss the finite element method when using basis functions enriched with a combination of exponential and hyperbolic functions. We show that the enrichment can resolve thermal boundary layers on coarse meshes and with few elements. Comparisons to the standard finite element method for thermal radiation in non-grey optically thick media with multi-frequency bands show the efficiency of the approach. Although we mainly study the enriched basis functions in glass cooling applications but the substantial saving in the computational requirements makes the approach highly relevant to a large number of engineering applications that involve solving the Rosseland model.

[1]  Axel Klar,et al.  A comparison of approximate models for radiation in gas turbines , 2004 .

[2]  M. S. Mohamed,et al.  A partition of unity finite element method for nonlinear transient diffusion problems in heterogeneous materials , 2019, Comput. Appl. Math..

[3]  D. Mihalas,et al.  Foundations of Radiation Hydrodynamics , 1985 .

[4]  R. Gama,et al.  An upper bound for the steady-state temperature for a class of heat conduction problems wherein the thermal conductivity is temperature dependent , 2013 .

[5]  Ralf Hiptmair,et al.  Dispersion analysis of plane wave discontinuous Galerkin methods , 2014 .

[6]  M. S. Mohamed,et al.  Mixed enrichment for the finite element method in heterogeneous media , 2015 .

[7]  Peter Monk,et al.  The Ultra Weak Variational Formulation Using Bessel Basis Functions , 2012 .

[8]  S. G. Martyushev,et al.  Characteristics of Rosseland and P-1 approximations in modeling nonstationary conditions of convection-radiation heat transfer in an enclosure with a local energy source , 2012 .

[9]  Owe Axelsson.  The method of diagonal compensation of reduced matrix entries and multilevel iteration , 1991 .

[10]  Jon Trevelyan,et al.  Wave boundary elements: a theoretical overview presenting applications in scattering of short waves , 2004 .

[11]  Sj Steven Hulshoff,et al.  The partition‐of‐unity method for linear diffusion and convection problems: accuracy, stabilization and multiscale interpretation , 2003 .

[12]  Mohammed Seaïd,et al.  A partition of unity finite element method for three-dimensional transient diffusion problems with sharp gradients , 2019, J. Comput. Phys..

[13]  Mohammed Seaïd,et al.  A three-dimensional enriched finite element method for nonlinear transient heat transfer in functionally graded materials , 2020 .

[14]  C. Farhat,et al.  A discontinuous Galerkin method with Lagrange multipliers for spatially-dependent advection-diffusion problems , 2017 .

[15]  M. Nawaz.  Role of hybrid nanoparticles in thermal performance of Sutterby fluid, the ethylene glycol , 2020 .

[16]  Mohammed Seaïd,et al.  Time-independent hybrid enrichment for finite element solution of transient conduction-radiation in diffusive grey media , 2013, J. Comput. Phys..

[17]  Axel Klar,et al.  Numerical Solvers for Radiation and Conduction in High Temperature Gas Flows , 2005 .

[18]  E. E. Anderson,et al.  Heat Transfer in Semitransparent Solids , 1975 .

[19]  A. Collin,et al.  Glass sagging simulation with improved calculation of radiative heat transfer by the optimized reciprocity Monte Carlo method , 2014 .

[20]  S. Saleem,et al.  Combined effects of partial slip and variable diffusion coefficient on mass and heat transfer subjected to chemical reaction , 2020, Physica Scripta.

[21]  W. Malalasekera,et al.  An Analysis of Thermal Radiation in Porous Media Under Local Thermal Non-equilibrium , 2020, Transport in Porous Media.

[22]  S. Duerigen,et al.  The simplified P3 approach on a trigonal geometry of the nodal reactor code DYN3D , 2012 .

[23]  J. Sanders.  Nonlinear, Transient Conduction Heat Transfer Using A Discontinuous Galerkin Hierarchical Finite Element Method , 2004 .

[24]  W. H. McCrea,et al.  Theoretical astrophysics : atomic theory and the analysis of stellar atmospheres and envelopes , 1936 .

[25]  Omar Laghrouche,et al.  Short wave modelling using special finite elements , 2000 .

[26]  Charbel Farhat,et al.  A discontinuous Galerkin method with Lagrange multipliers for the solution of Helmholtz problems in the mid-frequency regime , 2003 .

[27]  Carlos Armando Duarte,et al.  Transient analysis of sharp thermal gradients using coarse finite element meshes , 2011 .

[28]  Carlos Armando Duarte,et al.  Generalized finite element analysis of three-dimensional heat transfer problems exhibiting sharp thermal gradients , 2009 .

[29]  Rainer Koch,et al.  DISCRETE ORDINATES QUADRATURE SCHEMES FOR MULTIDIMENSIONAL RADIATIVE TRANSFER , 1995 .

[30]  Axel Klar,et al.  Simplified P N approximations to the equations of radiative heat transfer and applications , 2002 .

[31]  Mohammed Seaïd,et al.  An enriched finite element model with q-refinement for radiative boundary layers in glass cooling , 2014, J. Comput. Phys..

[32]  R. Al-Khoury,et al.  Time‐dependent shape functions for modeling highly transient geothermal systems , 2009 .

[33]  Carlos Armando Duarte,et al.  A two-scale generalized finite element method for parallel simulations of spot welds in large structures , 2018, Computer Methods in Applied Mechanics and Engineering.

[34]  Axel Klar,et al.  Adaptive solutions of SPN-approximations to radiative heat transfer in glass☆ , 2005 .

[35]  O. Laghrouche,et al.  A partition of unity FEM for time‐dependent diffusion problems using multiple enrichment functions , 2013 .

[36]  P. Weidler,et al.  Optical parameters for characterization of thermal radiation in ceramic sponges – Experimental results and correlation , 2014 .

[37]  Charbel Farhat,et al.  A higher‐order discontinuous enrichment method for the solution of high péclet advection–diffusion problems on unstructured meshes , 2010 .

[38]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[39]  P. S. Datti,et al.  MHD visco-elastic fluid flow over a non-isothermal stretching sheet , 2004 .

[40]  Mohammed Seaid.  Multigrid Newton-Krylov method for radiation in diffusive semitransparent media , 2007 .

[41]  Charbel Farhat,et al.  A discontinuous enrichment method for variable‐coefficient advection–diffusion at high Péclet number , 2011 .

[42]  M. Modest.  Radiative heat transfer , 1993 .

[43]  Fuqiang Wang,et al.  Thermal performance analyses of porous media solar receiver with different irradiative transfer models , 2014 .

[44]  M. S. Mohamed,et al.  Pollution studies for high order isogeometric analysis and finite element for acoustic problems , 2019, Computer Methods in Applied Mechanics and Engineering.

[45]  Abderrahman Bouhamidi,et al.  Iterative solvers for generalized finite element solution of boundary‐value problems , 2018, Numer. Linear Algebra Appl..