Researcher profile

K. H. Luo

K. H. Luo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2016arXiv

Thermal cascaded lattice Boltzmann method

In this paper, a thermal cascaded lattice Boltzmann method (TCLBM) is developed in combination with the double-distribution-function (DDF) approach. A density distribution function relaxed by the cascaded scheme is employed to solve the flow field, and a total energy distribution function relaxed by the BGK scheme is used to solve temperature field, where two distribution functions are coupled naturally. The forcing terms are incorporated by means of central moments, which is consistent with the previous force scheme [Premnath \emph{et al.}, Phys. Rev. E \textbf{80}, 036702 (2009)] but the derivation is more intelligible and the evolution process is simpler. In the method, the viscous heat dissipation and compression work are taken into account, the Prandtl number and specific-heat ratio are adjustable, the external force is considered directly without the Boussinesq assumption, and the low-Mach number compressible flows can also be simulated. The forcing scheme is tested by simulating a steady Taylor-Green flow, and the new TCLBM is then tested by numerical experiments of the thermal Couette flow, low-Mach number shock tube problem, and Rayleigh-Bénard convection. The simulation results agree well with the analytical solutions and/or results given by previous researchers.

preprint2015arXiv

Lattice Boltzmann methods for multiphase flow and phase-change heat transfer

Over the past few decades, tremendous progress has been made in the development of particle-based discrete simulation methods versus the conventional continuum-based methods. In particular, the lattice Boltzmann (LB) method has evolved from a theoretical novelty to a ubiquitous, versatile and powerful computational methodology for both fundamental research and engineering applications. It is a kinetic-based mesoscopic approach that bridges the microscales and macroscales, which offers distinctive advantages in simulation fidelity and computational efficiency. Applications of the LB method have been found in a wide range of disciplines including physics, chemistry, materials, biomedicine and various branches of engineering. The present work provides a comprehensive review of the LB method for thermofluids and energy applications, focusing on multiphase flows, thermal flows and thermal multiphase flows with phase change. The review first covers the theoretical framework of the LB method, revealing the existing inconsistencies and defects as well as common features of multiphase and thermal LB models. Recent developments in improving the thermodynamic and hydrodynamic consistency, reducing the spurious currents, enhancing the numerical stability, etc., are highlighted. These efforts have put the LB method on a firmer theoretical foundation with enhanced LB models that can achieve larger liquid-gas density ratio, higher Reynolds number and flexible surface tension. Examples of applications are provided in fuel cells and batteries, droplet collision, boiling heat transfer and evaporation, and energy storage. Finally, further developments and future prospect of the LB method are outlined for thermofluids and energy applications.

preprint2015arXiv

Lattice Boltzmann modeling of boiling heat transfer: The boiling curve and the effects of wettability

A hybrid thermal lattice Boltzmann (LB) model is presented to simulate thermal multiphase flows with phase change based on an improved pseudopotential LB approach [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 87, 053301 (2013)]. The present model does not suffer from the spurious term caused by the forcing-term effect, which was encountered in some previous thermal LB models for liquid-vapor phase change. Using the model, the liquid-vapor boiling process is simulated. The boiling curve together with the three boiling stages (nucleate boiling, transition boiling, and film boiling) is numerically reproduced in the LB community for the first time. The numerical results show that the basic features and the fundamental characteristics of boiling heat transfer are well captured, such as the severe fluctuation of transient heat flux in the transition boiling and the feature that the maximum heat transfer coefficient lies at a lower wall superheat than that of the maximum heat flux. Furthermore, the effects of the heating surface wettability on boiling heat transfer are investigated. It is found that an increase in contact angle promotes the onset of boiling but reduces the critical heat flux, and makes the boiling process enter into the film boiling regime at a lower wall superheat, which is consistent with the findings from experimental studies.

preprint2014arXiv

Contact angles in the pseudopotential lattice Boltzmann modeling of wetting

In this paper, we aim to investigate the implementation of contact angles in the pseudopotential lattice Boltzmann modeling of wetting at a large density ratio. The pseudopotential lattice Boltzmann model [X. Shan and H. Chen, Phys. Rev. E 49, 2941 (1994)] is a popular mesoscopic model for simulating multiphase flows and interfacial dynamics. In this model, the contact angle is usually realized by a fluid-solid interaction. Two widely used fluid-solid interactions: the density-based interaction and the pseudopotential-based interaction, as well as a modified pseudopotential-based interaction formulated in the present paper, are numerically investigated and compared in terms of the achievable contact angles, the maximum and the minimum densities, and the spurious currents. It is found that the pseudopotential-based interaction works well for simulating small static (liquid) contact angles, however, is unable to reproduce static contact angles close to 180 degrees. Meanwhile, it is found that the proposed modified pseudopotential-based interaction performs better in light of the maximum and the minimum densities and is overall more suitable for simulating large contact angles as compared with the other two types of fluid-solid interactions. Furthermore, the spurious currents are found to be enlarged when the fluid-solid interaction force is introduced. Increasing the kinematic viscosity ratio between the vapor and liquid phases is shown to be capable of reducing the spurious currents caused by the fluid-solid interactions.

preprint2014arXiv

Effect of the forcing term in the pseudopotential lattice Boltzmann modeling of thermal flows

The pseudopotential lattice Boltzmann (LB) model is a popular model in the LB community for simulating multiphase flows. Recently, several thermal LB models, which are based on the pseudopotential LB model and constructed within the framework of the double-distribution-function LB method, were proposed to simulate thermal multiphase flows [G. Házi and A. Márkus, Phys. Rev. E 77, 026305 (2008); L. Biferale et al., Phys. Rev. Lett. 108, 104502 (2012); S. Gong and P. Cheng, Int. J. Heat Mass Transfer 55, 4923 (2012)]. The objective of the present paper is to show that the effect of the forcing term on the temperature equation must be eliminated in the pseudopotential LB modeling of thermal flows. First, the effect of the forcing term on the temperature equation is shown via the Chapman-Enskog analysis. For comparison, alternative treatments that are free from the forcing-term effect are provided. Subsequently, numerical investigations are performed for two benchmark tests. The numerical results clearly show that the existence of the forcing-term effect will lead to significant numerical errors in the pseudopotential LB modeling of thermal flows.

preprint2013arXiv

Achieving tunable surface tension in the pseudopotential lattice Boltzmann modeling of multiphase flows

In this paper, we aim to address an important issue about the pseudopotential lattice Boltzmann (LB) model, which has attracted much attention as a mesoscopic model for simulating interfacial dynamics of complex fluids, but suffers from the problem that the surface tension cannot be tuned independently of the density ratio. In the literature, a multi-range potential was devised to adjust the surface tension [Sbragaglia et al., Phys. Rev. E 75, 026702 (2007)]. However, it was recently found that the density ratio of the system will be changed when the multi-range potential is employed to adjust the surface tension. A new approach is therefore proposed in the present work. The basic strategy is to add a source term to the LB equation so as to tune the surface tension of the pseudopotential LB model. The proposed approach can guarantee that the adjustment of the surface tension does not affect the mechanical stability condition of the pseudopotential LB model, and thus provides a separate control of the surface tension and the density ratio. Meanwhile, it still retains the mesoscopic feature and the computational simplicity of the pseudopotential LB model. Numerical simulations are carried out for stationary droplets, capillary waves, and droplet splashing on a thin liquid film. The numerical results demonstrate that the proposed approach is capable of achieving a tunable surface tension over a very wide range and can keep the density ratio unchanged when adjusting the surface tension.

preprint2013arXiv

Lattice Boltzmann modeling of multiphase flows at large density ratio with an improved pseudopotential model

Owing to its conceptual simplicity and computational efficiency, the pseudopotential multiphase lattice Boltzmann (LB) model has attracted significant attention since its emergence. In this work, we aim to extend the pseudopotential LB model to simulate multiphase flows at large density ratio and relatively high Reynolds number. First, based on our recent work [Li et al., Phys. Rev. E. 86, 016709 (2012)], an improved forcing scheme is proposed for the multiple-relaxation-time pseudopotential LB model in order to achieve thermodynamic consistency and large density ratio in the model. Next, through investigating the effects of the parameter a in the Carnahan-Starling equation of state, we find that the interface thickness is approximately proportional to 1/sqrt(a). Using a smaller a will lead to a wider interface thickness, which can reduce the spurious currents and enhance the numerical stability of the pseudopotential model at large density ratio. Furthermore, it is found that a lower liquid viscosity can be gained in the pseudopotential model by increasing the kinematic viscosity ratio between the vapor and liquid phases. The improved pseudopotential LB model is numerically validated via the simulations of stationary droplet and droplet oscillation. Using the improved model as well as the above treatments, numerical simulations of droplet splashing on a thin liquid film are conducted at a density ratio in excess of 500 with Reynolds numbers ranging from 40 to 1000. The dynamics of droplet splashing is correctly reproduced and the predicted spread radius is found to obey the power law reported in the literature.

preprint2012arXiv

Forcing scheme in pseudopotential lattice Boltzmann model for multiphase flows

The pseudo-potential lattice Boltzmann (LB) model is a widely used multiphase model in the LB community. In this model, an interaction force, which is usually implemented via a forcing scheme, is employed to mimic the molecular interactions that cause phase segregation. The forcing scheme is therefore expected to play an important role in the pseudo-potential LB model. In this paper, we aim to address some key issues about forcing schemes in the pseudo-potential LB model. Firstly, theoretical and numerical analyses will be made for Shan-Chen's forcing scheme and the exact-difference-method (EDM) forcing scheme. The nature of these two schemes and their recovered macroscopic equations will be shown. Secondly, through a theoretical analysis, we will reveal the physics behind the phenomenon that different forcing schemes exhibit different performances in the pseudo-potential LB model. Moreover, based on the analysis, we will present an improved forcing scheme and numerically demonstrate that the improved scheme can be treated as an alternative approach for achieving thermodynamic consistency in the pseudo-potential LB model.

preprint2012arXiv

Lattice Boltzmann method for relativistic hydrodynamics: Issues on conservation law of particle number and discontinuities

In this paper, we aim to address several important issues about the recently developed lattice Boltzmann (LB) model for relativistic hydrodynamics [M. Mendoza et al., Phys. Rev. Lett. 105, 014502 (2010); Phys. Rev. D 82, 105008 (2010)]. First, we study the conservation law of particle number in the relativistic LB model. Through the Chapman-Enskog analysis, it is shown that in the relativistic LB model the conservation equation of particle number is a convection-diffusion equation rather than a continuity equation, which makes the evolution of particle number dependent on the relaxation time. Furthermore, we investigate the origin of the discontinuities appeared in the relativistic problems with high viscosities, which were reported in a recent study [D. Hupp et al., Phys. Rev. D 84, 125015 (2011)]. A multiple-relaxation-time (MRT) relativistic LB model is presented to examine the influences of different relaxation times on the discontinuities. Numerical experiments show the discontinuities can be eliminated by setting the relaxation time $τ_e$ (related to the bulk viscosity) to be sufficiently smaller than the relaxation time $τ_v$ (related to the shear viscosity). Meanwhile, it is found that the relaxation time $τ_\varepsilon$, which has no effect on the conservation equations at the Navier-Stokes level, will affect the numerical accuracy of the relativistic LB model. Moreover, the accuracy of the relativistic LB model for simulating moderately relativistic problems is also investigated.

preprint2011arXiv

Additional interfacial force in lattice Boltzmann models for incompressible multiphase flows

The existing lattice Boltzmann models for incompressible multiphase flows are mostly constructed with two distribution functions, one is the order parameter distribution function, which is used to track the interface between different phases, and the other is the pressure distribution function for solving the velocity field. In this brief report, it is shown that in these models the recovered momentum equation is inconsistent with the target one: an additional interfacial force is included in the recovered momentum equation. The effects of the additional force are investigated by numerical simulations of droplet splashing on a thin liquid film and falling droplet under gravity. In the former test, it is found that the formation and evolution of secondary droplets are greatly affected, while in the latter the additional force is found to increase the falling velocity and limit the stretch of the droplet.

preprint2011arXiv

Coupling lattice Boltzmann model for simulation of thermal flows on standard lattices

In this paper, a coupling lattice Boltzmann (LB) model for simulating thermal flows on the standard 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 basic advantages 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.

preprint2011arXiv

Enhanced lattice Bhatnagar-Gross-Krook method for fluid dynamics simulation

A generalization of the lattice Bhatnagar-Gross-Krook (LBGK) model for the simulation of hydrodynamics is presented, which takes into account the difference and the frame-independence of the relaxation of non-hydrodynamic modes. The present model retains the computationally efficient standard LBGK form with the generalized equilibrium explicitly derived. The two-dimensional realization on the standard lattice is discussed in detail. Performance of the model is assessed through a shear layer simulation and enhanced stability and accuracy with respect to the standard LBGK are reported. The results demonstrate that the present model is a useful upgrade of the standard LBGK without compromising its computational efficiency and accuracy.