A three-dimensional color-fluid lattice Boltzmann model for immiscible two-phase flows is developed in the framework of a three-dimensional 27-velocity (D3Q27) lattice. The collision operator comprises the D3Q27 versions of three suboperators: a multiple-relaxation-time (MRT) collision operator, a generalized Liu–Valocchi–Kang perturbation operator, and a Latva-Kokko–Rothman recoloring operator. A D3Q27 version of an enhanced equilibrium distribution function is also incorporated into this model to improve the Galilean invariance. Three types of numerical tests, namely, a static droplet, an oscillating droplet, and the Rayleigh–Taylor instability, show a good agreement with analytical solutions and numerical simulations. Following these numerical tests, this model is applied to liquid-jet-breakup simulations. The simulation conditions are matched to the conditions of the previous experiments. In this case, numerical stability is maintained throughout the simulation, although the kinematic viscosity for the continuous phase is set as low as $1.8\times {10}^{-4}$, in which case the corresponding Reynolds number is $3.4\times {10}^3$; the developed lattice Boltzmann model based on the D3Q27 lattice enables us to perform the simulation with parameters directly matched to the experiments. The jet's liquid column transitions from an asymmetrical to an axisymmetrical shape, and entrainment occurs from the side of the jet. The measured time history of the jet's leading-edge position shows a good agreement with the experiments. Finally, the reproducibility of the regime map for liquid-liquid systems is assessed. The present lattice Boltzmann simulations well reproduce the characteristics of predicted regimes, including varicose breakup, sinuous breakup, and atomization.

When boiling occurs in a liquid flow field, the phenomenon is known as forced-convection boiling. We numerically investigate such a boiling system on a cylinder in a flow at saturated conditions. To deal with the complicated liquid–vapor phase-change phenomenon, we develop a numerical scheme based on the pseudopotential lattice Boltzmann method (LBM). The collision stage is performed in the space of central moments (CMs) to enhance numerical stability for high Reynolds numbers. The adopted forcing scheme, consistent with the CM-based LBM, leads to a concise yet robust algorithm. Furthermore, additional terms required to ensure thermodynamic consistency are derived in a CM framework. The effectiveness of the present scheme is successfully tested against a series of boiling processes, including nucleation, growth, and departure of a vapor bubble for Reynolds numbers varying between 30 and 30 000. Our CM-based LBM can reproduce all the boiling regimes, i.e., nucleate boiling, transition boiling, and film boiling, without any artificial input such as initial vapor phase. We find that the typical boiling curve, also known as the Nukiyama curve, appears even though the focused system is not the pool boiling but the forced-convection system. Also, our simulations support experimental observations of intermittent direct solid–liquid contact even in the film-boiling regime. Finally, we provide quantitative comparison with the semi-empirical correlations for the forced-convection film boiling on a cylinder on the Nu-Ja diagram.

We develop a lattice Boltzmann (LB) model for immiscible two-phase flow simulations with central moments (CMs). This successfully combines a three-dimensional nonorthogonal CM-based LB scheme [De Rosis, Phys. Rev. E 95, 013310 (2017)] with our previous color-gradient LB model [Saito, Abe, and Koyama, Phys. Rev. E 96, 013317 (2017)]. Hydrodynamic melt-jet breakup simulations show that the proposed model is significantly more stable, even for flow with extremely high Reynolds numbers, up to $O\left({10}^6\right)$. This enables us to investigate the phenomena expected under actual reactor conditions.

To better understand the fundamental interactions between melt jet and coolant during a core-disruptive accident at a sodium-cooled fast reactor, the jet breakup and droplet formation in immiscible liquid-liquid systems were studied experimentally. Experiments using two different pairs of test fluids were carried out at isothermal conditions. The observed jet breakup behavior was classified into characteristic regimes based on the classical Ohnesorge classification in liquid-gas systems. The variation in breakup length obtained in the present liquid-liquid system was similar to that in a liquid-gas system. The droplet size distribution in each breakup regime was analyzed using image processing and droplet formation via pinch-off, satellite formation, and entrainment was observed. The measured droplet size was compared with those available from melt jet experiments. Based on the observation and analysis results, the breakup regimes were organized on a dimensionless operating diagram, with the derived correlations representing the criteria for regime boundaries of a liquid-liquid system. Finally, the experimental data were extrapolated to the expected conditions of a sodium-cooled fast reactor. From this, it was implied that most of the hydrodynamic conditions during an accident would be close to the atomization regime, in which entrainment is the dominant process for droplet formation