1
|
A Non-Isothermal Chemical Lattice Boltzmann Model Incorporating Thermal Reaction Kinetics and Enthalpy Changes. COMPUTATION 2017. [DOI: 10.3390/computation5030037] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
2
|
Zhang L, Yang S, Zeng Z, Chen J, Yin L, Chew JW. Forcing scheme analysis for the axisymmetric lattice Boltzmann method under incompressible limit. Phys Rev E 2017; 95:043311. [PMID: 28505753 DOI: 10.1103/physreve.95.043311] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2016] [Indexed: 11/07/2022]
Abstract
Because the standard lattice Boltzmann (LB) method is proposed for Cartesian Navier-Stokes (NS) equations, additional source terms are necessary in the axisymmetric LB method for representing the axisymmetric effects. Therefore, the accuracy and applicability of the axisymmetric LB models depend on the forcing schemes adopted for discretization of the source terms. In this study, three forcing schemes, namely, the trapezium rule based scheme, the direct forcing scheme, and the semi-implicit centered scheme, are analyzed theoretically by investigating their derived macroscopic equations in the diffusive scale. Particularly, the finite difference interpretation of the standard LB method is extended to the LB equations with source terms, and then the accuracy of different forcing schemes is evaluated for the axisymmetric LB method. Theoretical analysis indicates that the discrete lattice effects arising from the direct forcing scheme are part of the truncation error terms and thus would not affect the overall accuracy of the standard LB method with general force term (i.e., only the source terms in the momentum equation are considered), but lead to incorrect macroscopic equations for the axisymmetric LB models. On the other hand, the trapezium rule based scheme and the semi-implicit centered scheme both have the advantage of avoiding the discrete lattice effects and recovering the correct macroscopic equations. Numerical tests applied for validating the theoretical analysis show that both the numerical stability and the accuracy of the axisymmetric LB simulations are affected by the direct forcing scheme, which indicate that forcing schemes free of the discrete lattice effects are necessary for the axisymmetric LB method.
Collapse
Affiliation(s)
- Liangqi Zhang
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
| | - Shiliang Yang
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
| | - Zhong Zeng
- Department of Engineering Mechanics, Colledge of Aerospace Engineering, Chongqing University, Chongqing 400044, People's Republic of China.,State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, People's Republic of China
| | - Jie Chen
- State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, People's Republic of China
| | - Linmao Yin
- Department of Engineering Mechanics, Colledge of Aerospace Engineering, Chongqing University, Chongqing 400044, People's Republic of China
| | - Jia Wei Chew
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore.,Singapore Membrane Technology Center, Nanyang Environment and Water Research Institute, Nanyang Technological University, Singapore 637141, Singapore
| |
Collapse
|
3
|
Zhang L, Yang S, Zeng Z, Yin L, Zhao Y, Chew JW. Consistent lattice Boltzmann methods for incompressible axisymmetric flows. Phys Rev E 2016; 94:023302. [PMID: 27627407 DOI: 10.1103/physreve.94.023302] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2016] [Indexed: 06/06/2023]
Abstract
In this work, consistent lattice Boltzmann (LB) methods for incompressible axisymmetric flows are developed based on two efficient axisymmetric LB models available in the literature. In accord with their respective original models, the proposed axisymmetric models evolve within the framework of the standard LB method and the source terms contain no gradient calculations. Moreover, the incompressibility conditions are realized with the Hermite expansion, thus the compressibility errors arising in the existing models are expected to be reduced by the proposed incompressible models. In addition, an extra relaxation parameter is added to the Bhatnagar-Gross-Krook collision operator to suppress the effect of the ghost variable and thus the numerical stability of the present models is significantly improved. Theoretical analyses, based on the Chapman-Enskog expansion and the equivalent moment system, are performed to derive the macroscopic equations from the LB models and the resulting truncation terms (i.e., the compressibility errors) are investigated. In addition, numerical validations are carried out based on four well-acknowledged benchmark tests and the accuracy and applicability of the proposed incompressible axisymmetric LB models are verified.
Collapse
Affiliation(s)
- Liangqi Zhang
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
- Department of Engineering Mechanics, College of Aerospace Engineering, Chongqing University, Chongqing 400044, People's Republic of China
| | - Shiliang Yang
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
| | - Zhong Zeng
- Department of Engineering Mechanics, College of Aerospace Engineering, Chongqing University, Chongqing 400044, People's Republic of China
- State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, People's Republic of China
| | - Linmao Yin
- Department of Engineering Mechanics, College of Aerospace Engineering, Chongqing University, Chongqing 400044, People's Republic of China
| | - Ya Zhao
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
| | - Jia Wei Chew
- School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore 637459, Singapore
- Singapore Membrane Technology Center, Nanyang Environment and Water Research Institute, Nanyang Technological University, Singapore 637141, Singapore
| |
Collapse
|
4
|
Karani H, Huber C. Lattice Boltzmann formulation for conjugate heat transfer in heterogeneous media. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:023304. [PMID: 25768633 DOI: 10.1103/physreve.91.023304] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/17/2014] [Indexed: 06/04/2023]
Abstract
In this paper, we propose an approach for studying conjugate heat transfer using the lattice Boltzmann method (LBM). The approach is based on reformulating the lattice Boltzmann equation for solving the conservative form of the energy equation. This leads to the appearance of a source term, which introduces the jump conditions at the interface between two phases or components with different thermal properties. The proposed source term formulation conserves conductive and advective heat flux simultaneously, which makes it suitable for modeling conjugate heat transfer in general multiphase or multicomponent systems. The simple implementation of the source term approach avoids any correction of distribution functions neighboring the interface and provides an algorithm that is independent from the topology of the interface. Moreover, our approach is independent of the choice of lattice discretization and can be easily applied to different advection-diffusion LBM solvers. The model is tested against several benchmark problems including steady-state convection-diffusion within two fluid layers with parallel and normal interfaces with respect to the flow direction, unsteady conduction in a three-layer stratified domain, and steady conduction in a two-layer annulus. The LBM results are in excellent agreement with analytical solution. Error analysis shows that our model is first-order accurate in space, but an extension to a second-order scheme is straightforward. We apply our LBM model to heat transfer in a two-component heterogeneous medium with a random microstructure. This example highlights that the method we propose is independent of the topology of interfaces between the different phases and, as such, is ideally suited for complex natural heterogeneous media. We further validate the present LBM formulation with a study of natural convection in a porous enclosure. The results confirm the reliability of the model in simulating complex coupled fluid and thermal dynamics in complex geometries.
Collapse
Affiliation(s)
- Hamid Karani
- School of Earth and Atmospheric Sciences, Georgia Institute of Technology, 311 Ferst Drive, Atlanta, Georgia 30332, USA
| | - Christian Huber
- School of Earth and Atmospheric Sciences, Georgia Institute of Technology, 311 Ferst Drive, Atlanta, Georgia 30332, USA
| |
Collapse
|
5
|
Li Q, Luo KH, He YL, Gao YJ, Tao WQ. Coupling lattice Boltzmann model for simulation of thermal flows on standard lattices. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:016710. [PMID: 22400704 DOI: 10.1103/physreve.85.016710] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/18/2011] [Revised: 11/14/2011] [Indexed: 05/31/2023]
Abstract
In this paper, a coupling lattice Boltzmann (LB) model for simulating thermal flows on the standard two-dimensional nine-velocity (D2Q9) lattice is developed in the framework of the double-distribution-function (DDF) approach in which the viscous heat dissipation and compression work are considered. In the model, a density distribution function is used to simulate the flow field, while a total energy distribution function is employed to simulate the temperature field. The discrete equilibrium density and total energy distribution functions are obtained from the Hermite expansions of the corresponding continuous equilibrium distribution functions. The pressure given by the equation of state of perfect gases is recovered in the macroscopic momentum and energy equations. The coupling between the momentum and energy transports makes the model applicable for general thermal flows such as non-Boussinesq flows, while the existing DDF LB models on standard lattices are usually limited to Boussinesq flows in which the temperature variation is small. Meanwhile, the simple structure and general features of the DDF LB approach are retained. The model is tested by numerical simulations of thermal Couette flow, attenuation-driven acoustic streaming, and natural convection in a square cavity with small and large temperature differences. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.
Collapse
Affiliation(s)
- Q Li
- Key Laboratory of Thermal Fluid Science and Engineering of MOE, National Key Laboratory of Multiphase Flow in Power Engineering, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi 710049, China
| | | | | | | | | |
Collapse
|
6
|
Zhou JG. Axisymmetric lattice Boltzmann method revised. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:036704. [PMID: 22060526 DOI: 10.1103/physreve.84.036704] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/19/2011] [Revised: 07/19/2011] [Indexed: 05/31/2023]
Abstract
A reformulated lattice Boltzmann model is described for incompressible axisymmetric flows with or without swirling. It is a further development and improvement on the author's axisymmetric lattice Boltzmann method. The main features of the revised scheme are (a) all the macroscopic variables such as velocities are determined in the same formulas as those in the conventional lattice Boltzmann approach to the Navier-Stokes equations and (b) the added sink or source and force terms are simple with no calculation for a derivative. Such features distinguish the present method from the other existing simplified schemes, leading to a simple and efficient model. The scheme is naturally suitable for generic incompressible axisymmetric rotational flows involving more physical phenomena. The numerical solutions to Womersley and cylindrical cavity flows are presented to demonstrate its accuracy and capability.
Collapse
Affiliation(s)
- Jian Guo Zhou
- School of Engineering, University of Liverpool, Liverpool L69 3GQ, United Kingdom
| |
Collapse
|
7
|
Alloui Z, Guiet J, Vasseur P, Reggio M. Natural convection of nanofluids in a shallow rectangular enclosure heated from the side. CAN J CHEM ENG 2011. [DOI: 10.1002/cjce.20523] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
8
|
Li Q, He YL, Tang GH, Tao WQ. Improved axisymmetric lattice Boltzmann scheme. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:056707. [PMID: 20866354 DOI: 10.1103/physreve.81.056707] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/30/2009] [Indexed: 05/29/2023]
Abstract
This paper proposes an improved lattice Boltzmann scheme for incompressible axisymmetric flows. The scheme has the following features. First, it is still within the framework of the standard lattice Boltzmann method using the single-particle density distribution function and consistent with the philosophy of the lattice Boltzmann method. Second, the source term of the scheme is simple and contains no velocity gradient terms. Owing to this feature, the scheme is easy to implement. In addition, the singularity problem at the axis can be appropriately handled without affecting an important advantage of the lattice Boltzmann method: the easy treatment of boundary conditions. The scheme is tested by simulating Hagen-Poiseuille flow, three-dimensional Womersley flow, Wheeler benchmark problem in crystal growth, and lid-driven rotational flow in cylindrical cavities. It is found that the numerical results agree well with the analytical solutions and/or the results reported in previous studies.
Collapse
Affiliation(s)
- Q Li
- National Key Laboratory of Multiphase Flow in Power Engineering, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi 710049, China
| | | | | | | |
Collapse
|
9
|
Li Q, He YL, Tang GH, Tao WQ. Lattice Boltzmann model for axisymmetric thermal flows. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:037702. [PMID: 19905256 DOI: 10.1103/physreve.80.037702] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/17/2009] [Indexed: 05/28/2023]
Abstract
A thermal lattice Boltzmann (LB) model is presented for axisymmetric thermal flows in the incompressible limit. The model is based on the double-distribution-function LB method, which has attracted much attention since its emergence for its excellent numerical stability over the multispeed LB method. Compared with the existing axisymmetric thermal LB models, the present model is simpler and retains the inherent features of the standard LB method. Numerical simulations are carried out for the thermally developing laminar flows in circular ducts and the natural convection in an annulus between two coaxial vertical cylinders. The Nusselt number obtained from the simulations agrees well with the analytical solutions and/or the results reported in previous studies.
Collapse
Affiliation(s)
- Q Li
- National Key Laboratory of Multiphase Flow in Power Engineering, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi 710049, China
| | | | | | | |
Collapse
|
10
|
Chen S, Tölke J, Krafczyk M. Simple lattice Boltzmann subgrid-scale model for convectional flows with high Rayleigh numbers within an enclosed circular annular cavity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:026702. [PMID: 19792276 DOI: 10.1103/physreve.80.026702] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/20/2009] [Indexed: 05/28/2023]
Abstract
Natural convection within an enclosed circular annular cavity formed by two concentric vertical cylinders is of fundamental interest and practical importance. Generally, the assumption of axisymmetric thermal flow is adopted for simulating such natural convections and the validity of the assumption of axisymmetric thermal flow is still held even for some turbulent convection. Usually the Rayleigh numbers (Ra) of realistic flows are very high. However, the work to design suitable and efficient lattice Boltzmann (LB) models on such flows is quite rare. To bridge the gap, in this paper a simple LB subgrid-scale (SGS) model, which is based on our recent work [S. Chen, J. Tölke, and M. Krafczyk, Phys. Rev. E 79, 016704 (2009); S. Chen, J. Tölke, S. Geller, and M. Krafczyk, Phys. Rev. E 78, 046703 (2008)], is proposed for simulating convectional flow with high Ra within an enclosed circular annular cavity. The key parameter for the SGS model can be quite easily and efficiently evaluated by the present model. The numerical experiments demonstrate that the present model works well for a large range of Ra and Prandtl number (Pr). Though in the present study a popularly used static Smagorinsky turbulence model is adopted to demonstrate how to develop a LB SGS model for simulating axisymmetric thermal flows with high Ra, other state-of-the-art turbulence models can be incorporated into the present model in the same way. In addition, the present model can be extended straightforwardly to simulate other axisymmetric convectional flows with high Ra, for example, turbulent convection with internal volumetric heat generation in a vertical cylinder, which is an important simplified representation of a nuclear reactor.
Collapse
Affiliation(s)
- Sheng Chen
- State Key Laboratory of Coal Combustion, Huazhong University of Science and Technology, Wuhan 430074, China
| | | | | |
Collapse
|
11
|
Chen S, Tölke J, Krafczyk M. Numerical simulation of fluid flow and heat transfer inside a rotating disk-cylinder configuration by a lattice Boltzmann model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:016702. [PMID: 19658833 DOI: 10.1103/physreve.80.016702] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/16/2009] [Indexed: 05/28/2023]
Abstract
A simple lattice Boltzmann model for numerical simulation of fluid flow and heat transfer inside a rotating disk-cylinder configuration, which is of fundamental interest and practical importance in science as well as in engineering, is proposed in this paper. Unlike existing lattice Boltzmann models for such flows, which were based on "primitive-variable" Navier-Stokes equations, the target macroscopic equations of the present model for the flow field are vorticity-stream function equations, inspired by our recent work designed for nonrotating flows [S. Chen, J. Tölke, and M. Krafczyk, Phys. Rev. E 79, 016704 (2009); S. Chen, J. Tölke, S. Geller, and M. Krafczyk, Phys. Rev. E 78, 046703 (2008)]. The flow field and the temperature field both are solved by the D2Q5 model. Compared with the previous models, the present model is more efficient, more stable, and much simpler. It was found that, even though with a relatively low grid resolution, the present model can still work well when the Grashof number is very high. The advantages of the present model are validated by numerical experiments.
Collapse
Affiliation(s)
- Sheng Chen
- State Key Laboratory of Coal Combustion, Huazhong University of Science and Technology, Wuhan 430074, China
| | | | | |
Collapse
|