A coupled high-order implicit-explicit flux reconstruction lattice Boltzmann method for nearly incompressible thermal flows. (15th May 2022)

• Record Type: Journal Article
• Authors: Ma, Chao; Wu, Jie; Yang, Liming; Dong, Hao
• Is Part Of: International journal of heat and mass transfer. Volume 187(2022)
• Publication Date: 2022-05-15
• Subjects: Lattice Boltzmann method -- Double distribution function model -- Flux reconstruction scheme -- Implicit-explicit Runge-Kutta scheme -- High order accuracy -- Incompressible thermal flow; Heat -- Transmission -- Periodicals; Mass transfer -- Periodicals; Chaleur -- Transmission -- Périodiques; Transfert de masse -- Périodiques; Electronic journals; 621.4022
• DOI: 10.1016/j.ijheatmasstransfer.2022.122575
• Languages: English
• ISSNs: 0017-9310

Abstract

Highlights: A high-order implicit-explicit flux reconstruction thermal lattice Boltzmann method (FRTLBM) is proposed. FRTLBM can easily achieve high-order accuracy (>2). FRTLBM has a good adaptability of complicated geometries. FRTLBM is a reliable and efficient tool for incompressible thermal flows. Abstract: In this paper, based on the double distribution function (DDF) model and the Boussinesq approximation, a high-order implicit-explicit flux reconstruction thermal lattice Boltzmann method (FRTLBM) is proposed for simulating nearly incompressible thermal flows. Two discrete velocity Boltzmann equations (DVBEs), one governing the flow field and the other governing the temperature field, are solved by using a high-order flux reconstruction scheme for spatial discretization and a second-order implicit-explicit Runge-Kutta scheme (IMEX RK) for time discretization in the generalized curvilinear coordinates with uniform and non-uniform grids. Numerical validations of the proposed method are implemented by simulating (a) porous plate problem, (b) natural convection in a square cavity, (c) Rayleigh-Bénard convection, (d) natural convection in a concentric annulus and (e) mixed heat transfer from a heated circular cylinder. Numerical results demonstrate that the present method can achieve the high-order accuracy and it is stable and accurate even at the relatively high Rayleigh number (Ra = 10 7 and 10 8 ). The present method is efficient due to the small amount of mesh, the simple algorithm and the time step unlimited by the relaxation time. Since the rescontruction is conducted in a single cell, the present method is compact for parallel computing. Moreover, the present method has a good adaptability of complicated geometries.