Numerical Simulation of Melting Problems Using the Lattice Boltzmann Method with the Interfacial Tracking Method

International Journal of Numerical Methods for Heat &Amp; Fluid Flow

2016

Self Cite

Purpose: The purposes of this paper are testing an efficiency algorithm based on LBM and using it to analyze two-dimensional natural convection with low Prandtl number. Design/methodology/approach: Steady state or oscillatory results are obtained using double multiple-relaxation-time thermal lattice Boltzmann method. The velocity and temperature fields are solved using D2Q9 and D2Q5 models, respectively. Findings: With different Rayleigh number, the tested natural convection can either achieve to steady state or oscillatory. With fixed Rayleigh number, lower Prandtl number leads to a weaker convection effect, longer oscillation period and higher oscillation amplitude for the cases reaching oscillatory solutions. At fixed Prandtl number, higher Rayleigh number leads to a more notable convection effect and longer oscillation period. Originality/value: Double multiple-relaxation-time thermal lattice Boltzmann method is applied to simulate the low Prandtl number (0.001 -0.01) fluid natural convection. Rayleigh number and Prandtl number effects are also investigated when the natural convection results oscillate.

Section: Introductionmentioning

confidence: 99%

“…It has more possibility to reach oscillatory results. Li et al (2015) discussed low Prandtl number melting problem with double LBGK model. Kosec and Sarler (2013) reported the solution of a low Prandtl number natural convection benchmark problem.…”

Section: Introductionmentioning

confidence: 99%

Double MRT thermal lattice Boltzmann method for simulating natural convection of low Prandtl number fluids

International Journal of Numerical Methods for Heat &Amp; Fluid Flow

2016

Self Cite

“…Different models exist in LBM to solve heat transfer problems. Li et al [13,14] solve melting problems using FVM-LBM hybrid method and double distribution model in LBM with interfacial tracking method. All of the above applications of the interfacial tracking method are based on 2-D problems.…”

Section: Introductionmentioning

confidence: 99%

Lattice Boltzmann Method Simulation of 3-D Melting Using Double MRT Model With Interfacial Tracking Method

Volume 2: Heat Transfer in Multiphase Systems; Gas Turbine Heat Transfer; Manufacturing and Materials Processing; Heat Transfer

2016

Self Cite

Three-dimensional melting problems are investigated numerically with Lattice Boltzmann method (LBM). Regarding algorithm's accuracy and stability, Multiple-Relaxation-Time (MRT) models are employed to simplify the collision term in LBM. Temperature and velocity fields are solved with double distribution functions, respectively. 3-D melting problems are solved with double MRT models for the first time in this article.The key point for the numerical simulation of a melting problem is the methods to obtain the location of the melting front and this article uses interfacial tracking method. The interfacial tracking method combines advantages of both deforming and fixed grid approaches. The location of the melting front was obtained by calculating the energy balance at the solid-liquid interface. Various 3-D conduction controlled melting problems are solved firstly to verify the numerical method. Liquid fraction tendency and temperature distribution obtained from numerical methods agree with the analytical results well. The proposed double MRT model with interfacial tracking method is valid to solve 3-D melting problems. Different 3-D convection controlled melting problems are then solved with the proposed numerical method. Various locations of the heat surface have different melting front moving velocities, due to the natural convection effects. Rayleigh number's effects to the 3-D melting process is discussed.

“…Fusegi et al [9] used SIMPLE algorithm with QUICK scheme to simulate natural convection in a three-dimensional square cavity. Li et al [10] discussed the nonlinear results of natural convection melting problems in a square cavity. Ishida et al [11] and Benouaguef et al [12] studied nonlinear characteristics of natural convection in inflatable square cavity using different numerical methods.…”

Section: Introductionmentioning

confidence: 99%

Effects of slotted structures on the nonlinear characteristics of natural convection in a cylinder with an internal concentric slotted annulus

Shen

Numerical Heat Transfer, Part A: Applications

et al. 2016

Self Cite

Natural convection in a cylinder with an internally slotted annulus was solved by SIMPLE algorithm, and the effects of different slotted structures on nonlinear characteristics of natural convection were investigated. The results show that the equivalent thermal conductivity Keq increases with Rayleigh number, and reaches the maximum in the vertical orientation. Nonlinear results were obtained by simulating the fluid flow at different conditions. With increasing Rayleigh number, heat transfer is intensified and the state of heat transfer changes from the steady to unsteady. We investigated different slotted structures effects on natural convection, and analyze the corresponding nonlinear characteristics. Keywords