No Arabic abstract
In this paper, a diffuse-interface lattice Boltzmann method (DI-LBM) is developed for fluid-particle interaction problems. In this method, the sharp interface between the fluid and solid is replaced by a thin but nonzero thickness transition region named diffuse interface, where the physical variables varies continuously. In order to describe the diffuse interface, we introduce a smooth function, which is similar to the order parameter in phase-field model or the volume fraction of solid phase in the partially saturated lattice Boltzmann method (PS-LBM). In addition, to depict the fluid-particle interaction more accurately, a modified force term is also proposed and included in the evolution equation of the DI-LBM. Some classical problems are used to test the DI-LBM, and the results are in good agreement with some available theoretical and numerical works. Finally, it is also found that the DI-LBM is more efficient and accurate than the PS-LBM with the superposition model.
We present a novel method for fluid structure interaction (FSI) simulations where an original 2nd-order curved space lattice Boltzmann fluid solver (LBM) is coupled to a finite element method (FEM) for thin shells. The LBM can work independently on a standard lattice in curved coordinates without the need for interpolation, re-meshing or an immersed boundary. The LBM distribution functions are transformed dynamically under coordinate change. In addition, force and momentum can be calculated on the nodes exactly in any geometry. Furthermore, the FEM shell is a complete numerical tool with implementations such as growth, self-contact and strong external forces. We show resolution convergent error for standard tests under metric deformation. Mass and volume conservation, momentum transfer, boundary-slip and pressure maintenance are verified through specific examples. Additionally, a brief deformation stability analysis is carried out. Next, we study the interaction of a square fluid flow channel to a deformable shell. Finally, we simulate a flag at moderate Reynolds number, air flow channel. The scheme is limited to small deformations of O(10%) relative to domain size, by improving its stability the method can be naturally extended to multiple applications without further implementations.
We propose a novel approach to the numerical simulation of thin film flows, based on the lattice Boltzmann method. We outline the basic features of the method, show in which limits the expected thin film equations are recovered and perform validation tests. The numerical scheme is applied to the viscous Rayleigh-Taylor instability of a thin film and to the spreading of a sessile drop towards its equilibrium contact angle configuration. We show that the Cox-Voinov law is satisfied, and that the effect of a tunable slip length on the substrate is correctly captured. We address, then, the problem of a droplet sliding on an inclined plane, finding that the Capillary number scales linearly with the Bond number, in agreement with experimental results. At last, we demonstrate the ability of the method to handle heterogenous and complex systems by showcasing the controlled dewetting of a thin film on a chemically structured substrate.
An Onsager-like relation is proposed as a new criterion for constructing and analysing the lattice Boltzmann (LB) method. For LB models obeying the relation, we analyse their linearized stability, establish their diffusive limit, and find new constraints for those with free parameters. The new relation seems of fundamental importance for the LB method.
In this paper, we develop a simplified hybrid weighted essentially non-oscillatory (WENO) method combined with the modified ghost fluid method (MGFM) [28] to simulate the compressible two-medium flow problems. The MGFM can turn the two-medium flow problems into two single-medium cases by defining the ghost fluids status in terms of the predicted the interface status, which makes the material interface invisible. For the single medium flow case, we adapt between the linear upwind scheme and the WENO scheme automatically by identifying the regions of the extreme points for the reconstruction polynomial as same as the hybrid WENO scheme [50]. Instead of calculating their exact locations, we only need to know the regions of the extreme points based on the zero point existence theorem, which is simpler for implementation and saves computation time. Meanwhile, it still keeps the robustness and has high efficiency. Extensive numerical results for both one and two dimensional two-medium flow problems are performed to demonstrate the good performances of the proposed method.
In this paper, we first present a unified framework for the modelling of generalized lattice Boltzmann method (GLBM). We then conduct a comparison of the four popular analysis methods (Chapman-Enskog analysis, Maxwell iteration, direct Taylor expansion and recurrence equations approaches) that have been used to obtain the macroscopic Navier-Stokes equations and nonlinear convection-diffusion equations from the GLBM, and show that from mathematical point of view, these four analysis methods are equivalent to each other. Finally, we give some elements that are needed in the implementation of the GLBM, and also find that some available LB models can be obtained from this GLBM.