An improved axisymmetric interfacial lattice Boltzmann flux solver for large-density-ratio multiphase flows

In this paper, an improved axisymmetric interfacial lattice Boltzmann flux solver abandoning the previous predictor-corrector scheme is proposed. Unlike the previous model starting from the two-dimensional standard lattice Boltzmann method (LBM), the present method is developed using the axisymmetric LBM, which directly incorporates the axisymmetric effects into the distribution functions. As a result, the proposed solver does not need the corrector step involving complicated space derivatives. It makes this method simpler and more computationally efficient. In the present solver, the resultant governing equation is globally resolved by the finite volume method, while the fluxes are reconstructed by local application of the axisymmetric LBM. Therefore, the inconsistency between the local reconstruction and the global governing equation can be eliminated because the global equation can be strictly derived from the axisymmetric LBM, which holds stronger physical basis than the previous method. Numerical experiments about the interface capturing and the multiphase flows are conducted to test the proposed model. Results show that the present method is superior to the fractional step solver in terms of the accuracy, stability, and computational efficiency. In addition, this solver has the capacity of simulating large-density-ratio and complex interfacial change.