A lattice Boltzmann method for axisymmetric thermocapillary flows

In this work, we develop a two-phase lattice Boltzmann method (LBM) to simulate axisymmetric thermocapil- lary flows. This method simulates the immiscible axisymmetric two-phase flow by an improved color-gradient model, in which the single-phase collision, perturbation and recoloring operators are all presented with the axisymmetric effect taken into account in a simple and computational consistent manner. An additional lattice Boltzmann equation is introduced to describe the evolution of the axisymmetric temperature field, which is coupled to the hydrodynamic equations through an equation of state. This method is first validated by simulations of Rayleigh-B ́enard convection in a vertical cylinder and thermocapillary migration of a de- formable droplet at various Marangoni numbers. It is then used to simulate the thermocapillary migration of two spherical droplets in a constant applied temperature gradient along their line of centers, and the influence of the Marangoni number (Ca), initial distance between droplets (S0), and the radius ratio of the leading to trailing droplets (Λ) on the migration process is systematically studied. As Ma increases, the thermal wake behind the leading droplet strengthens, resulting in the transition of the droplet migration from coalescence to non-coalescence; and also, the final distance between droplets increases with Ma for the non-coalescence cases. The variation of S0 does not change the final state of the droplets although it has a direct impact on the migration process. In contrast, Λ can significantly influence the migration process of both droplets and their final state: at low Ma, decreasing Λ favors the coalescence of both droplets; at high Ma, the two droplets do not coalesce eventually but migrate with the same velocity for the small values of Λ, and decreasing Λ leads to a shorter equilibrium time and a faster migration velocity

A lattice Boltzmann method for axisymmetric multicomponent flows with high viscosity ratio

A color-gradient lattice Boltzmann method (LBM) is proposed to simulate ax- isymmetric multicomponent flows. This method uses a collision operator that is a combination of three separate parts, namely single-component collision op- erator, perturbation operator, and recoloring operator. A source term is added into the single-component collision operator such that in each single-component region the axisymmetric continuity and momentum equations can be exactly re- covered. The interfacial tension effect is realized by the perturbation operator, in which an interfacial force of axisymmetric form is derived using the concept of continuum surface force. A recoloring operator proposed by Latva-Kokko and Rothman is extended to the axisymmetric case for phase segregation and maintenance of the interface. To enhance the method’s numerical stability for handling binary fluids with high viscosity ratio, a multiple-relaxation-time mod- el is used for the collision operator. Several numerical examples, including static droplet test, oscillation of a viscous droplet, and breakup of a liquid thread, are presented to test the capability and accuracy of the proposed color-gradient LB- M. It is found that the present method is able to accurately capture the phase interface and produce low spurious velocities. Also, the LBM results are all in good agreement with the analytical solutions and/or available experimental data for a very broad range of viscosity ratios

Lattice Boltzmann modeling of contact angle and its hysteresis in two-phase flow with large viscosity difference

Contact angle hysteresis is an important physical phenomenon omnipresent in nature and various industrial processes, but its effects are not considered in many existing multiphase flow simulations due to modeling complexity. In this work, a multiphase lattice Boltzmann method (LBM) is developed to simulate the contact-line dynamics with consideration of the contact angle hysteresis for a broad range of kinematic viscosity ratios. In this method, the immiscible two-phase flow is described by a color-fluid model, in which the multiple-relaxation-time collision operator is adopted to increase numerical stability and suppress unphysical spurious currents at the contact line. The contact angle hysteresis is introduced using the strategy proposed by Ding and Spelt [Ding and Spelt, J. Fluid Mech. 599, 341 (2008)], and the geometrical wetting boundary condition is enforced to obtain the desired contact angle. This method is first validated by simulations of static contact angle and dynamic capillary intrusion process on ideal (smooth) surfaces. It is then used to simulate the dynamic behavior of a droplet on a nonideal (inhomogeneous) surface subject to a simple shear flow. When the droplet remains pinned on the surface due to hysteresis, the steady interface shapes of the droplet quantitatively agree well with the previous numerical results. Four typical motion modes of contact points, as observed in a recent study, are qualitatively reproduced with varying advancing and receding contact angles. The viscosity ratio is found to have a notable impact on the droplet deformation, breakup, and hysteresis behavior. Finally, this method is extended to simulate the droplet breakup in a microfluidic T junction, with one half of the wall surface ideal and the other half nonideal. Due to the contact angle hysteresis, the droplet asymmetrically breaks up into two daughter droplets with the smaller one in the nonideal branch channel, and the behavior of daughter droplets is significantly different in both branch channels. Also, it is found that the contact angle hysteresis is strengthened with decreasing the viscosity ratio, leading to an earlier droplet breakup and a decrease in the maximum length that the droplet can reach before the breakup. These simulation results manifest that the present multiphase LBM can be a useful substitute to Ba et al. [Phys. Rev. E 88, 043306 (2013)] for modeling the contact angle hysteresis, and it can be easily implemented with higher computational efficiency

Optimized detector tomography for photon-number resolving detectors with hundreds of pixels

Photon-number resolving detectors with hundreds of pixels are now readily available, while the characterization of these detectors using detector tomography is computationally intensive. Here, we present a modified detector tomography model that reduces the number of variables that need optimization. To evaluate the effectiveness and accuracy of our model, we reconstruct the photon number distribution of optical coherent and thermal states using the expectation-maximization-entropy algorithm. Our results indicate that the fidelity of the reconstructed states remains above 99%, and the second and third-order correlations agree well with the theoretical values for a mean number of photons up to 100. We also investigate the computational resources required for detector tomography and find out that our approach reduces the solving time by around a half compared to the standard detector tomography approach, and the required memory resources are the main obstacle for detector tomography of a large number of pixels. Our results suggest that detector tomography is viable on a supercomputer with 1~TB RAM for detectors with up to 340 pixels