Great advances have been made with the lattice Boltzmann (LB) method for complicated fluid phenomena and fundamental thermal processes over the past three decades. This paper presents a systematic overview of the LB method from 1990 to 2018, based on bibliometric analysis and the Science Citation Index Expanded (SCI-E) database. The results show that China took the leading position in this field, followed by the USA and UK. The Chinese Academy of Sciences had the most publications, while the Los Alamos National Laboratory was first as far as highest average citation per paper and h-index are concerned. Physical Review E was the most productive journal and "Mechanics" was the most frequently used subject category. Keyword analysis indicated that recent research has focused on the natural convection and heat transfer of nanofluid or multiphase flow in complex porous media. Hydrothermal treatment of nanofluid with shape factor on the conditions, such as variable magnetic fields, thermal radiation and slipping boundary, were the research hotspots. Further research perspectives mainly explore the multiscale models for coupling multiple transport phenomena, morphology optimization of porous parameters, new nanoparticles with shape factor, multicomponent LB method considering Knudsen diffusion effect, LB-based hybrid methods, radiation performance or boiling-heat transfer of nanofluid, and the active control of droplets, may continue to attract more attention. Moreover, some new applications, such as phase change of metal foam, erosion induced by nanaofluid, anode circulating, 3D modeling in thermal systems with vibration, and magnetohydrodynamics microfluid devices, could be of interest going forward.

The bubble bursting process existing in the particle flow is a complex gas-liquid-solid three-phase coupling dynamic problem. The bubble bursting mechanism, including dynamic characteristics and wall effects, is not clear. To address the above matters, we present a modeling method for the piecewise linear interface calculation-volume of fluid (PLIC-VOF) based bubble burst. The bubble bursting process near or on the wall is analyzed to reveal the dynamic characteristics of bubble bursting and obtain the effect of a bubble bursting on the surrounding flow field. Then a particle image velocimetry (PIV) based self-developed experimental observation platform is established, and the effectiveness of the proposed method is verified. Research results indicate that, in the high-speed turbulent environment, a large pressure difference existed in the bubble tail, which induces the bubble burst to occur; the distance between the wall and the bubble decreases; the higher the flow velocity is, the less time is acquired for bubble bursting, but when the flow velocity exceeds the critical velocity 50 m/s, more time is needed; the coalescence-burst process of double bubbles increases the bubble bursting time, which causes the acceleration of particle motion to reduce.

Microreactors have the advantages of high heat and mass transfer efficiency, strict control of reaction parameters, easy amplification, and good safety performance, and have been widely used in various fields such as chip manufacturing, fine chemicals, and biomanufacturing. However, narrow microchannels in microreactors often become filled with catalyst particles, leading to blockages. To address this challenge, this study proposes a multiphase flow heat transfer model based on the lattice Boltzmann method (LBM) to investigate the dynamic changes during the bubble collapse process and temperature distribution regularities. Based on the developed three-phase flow dynamics model, this study delves into the shock dynamic evolution process of bubble collapse and analyzes the temperature distribution regularities. Then, the flow patterns under different particle density conditions are explored. The study found that under the action of shock wave, the stable structure of the liquid film of the bubble is destroyed, and the bubble deforms and collapses. At the moment of bubble collapse, energy is rapidly transferred from the potential energy of the bubble to the kinetic energy of the flow field. Subsequently, the kinetic energy is converted into pressure waves. This results in the rapid generation of extremely high pressure in the flow field, creating high-velocity jets and intense turbulent vortices, which can enhance the mass transfer effects of the multiphase flows. At the moment of bubble collapse, a certain high temperature phenomenon will be formed at the collapse, and the high temperature phenomenon in this region is relatively chaotic and random. The pressure waves generated during bubble collapse have a significant impact on the motion trajectories of particles, while the influence on high-density particles is relatively small. The results offer a theoretical basis for understanding mass transfer mechanisms and particle flow patterns in three-phase flow. Moreover, these findings have significant practical implications for advancing technologies in industrial applications, including chip manufacturing and chemical process transport.

Basic oxygen furnace (BOF) steelmaking is widely used in the metallurgy field. The slagging reaction is a necessary process that oxidizes C, Mn, Si, P, S, and other impurities and therefore directly affects the quality of the resultant steel. Relevant research has suggested that intensifying the stirring effect can accelerate the slagging reaction and that the dynamic characteristics of the top blow are the key factor in exploring the related complex physical and chemical phenomena. To address the issue, the standard k-$\omega$ turbulence model and level-set method were adopted in the present work and a fluid dynamics model was developed for a BOF. Accordingly, the slag–metal–gas emulsion interaction and stirring effect were investigated, and the interference mechanism of a multi-nozzle supersonic coherent jet was revealed. Finally, a self-adjustment method based on fuzzy control is proposed for the oxygen lance. The results indicate that the transfer efficiency of jet kinetic energy at the gas–liquid interface is the critical factor for the slagging reaction and that multi-nozzle oxygen lances with a certain twisted angle have important advantages with respect to stirring effects and splashing inhibition. The fuzzy control method predicts that the optimal nozzle twist angle is within the range of 7.2° to 7.8°. The results presented herein can provide theoretical support and beneficial reference information for BOF steelmaking.