No Arabic abstract
In this brief report, a thermal lattice-Boltzmann (LB) model is presented for axisymmetric thermal flows in the incompressible limit. The model is based on the double-distribution-function LB method, which has attracted much attention since its emergence for its excellent numerical stability. Compared with the existing axisymmetric thermal LB models, the present model is simpler and retains the inherent features of the standard LB method. Numerical simulations are carried out for the thermally developing laminar flows in circular ducts and the natural convection in an annulus between two coaxial vertical cylinders. The Nusselt number obtained from the simulations agrees well with the analytical solutions and/or the results reported in previous studies.
In this paper, a coupling lattice Boltzmann (LB) model for simulating thermal flows on the standard D2Q9 lattice is developed in the framework of the double-distribution-function (DDF) approach in which the viscous heat dissipation and compression work are considered. In the model, a density distribution function is used to simulate the flow field, while a total energy distribution function is employed to simulate the temperature field. The discrete equilibrium density and total energy distribution functions are obtained from the Hermite expansions of the corresponding continuous equilibrium distribution functions. The pressure given by the equation of state of perfect gases is recovered in the macroscopic momentum and energy equations. The coupling between the momentum and energy transports makes the model applicable for general thermal flows such as non-Boussinesq flows, while the existing DDF LB models on standard lattices are usually limited to Boussinesq flows in which the temperature variation is small. Meanwhile, the simple structure and basic advantages of the DDF LB approach are retained. The model is tested by numerical simulations of thermal Couette flow, attenuation-driven acoustic streaming, and natural convection in a square cavity with small and large temperature differences. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.
In this article, a coupled Two-relaxation-time Lattice Boltzmann-Volume penalization (TRT-LBM-VP) method is presented to simulate flows past obstacles. Two relaxation times are used in the collision operator, of which one is related to the fluid viscosity and the other one is related to the numerical stability and accuracy. The volume penalization method is introduced into the TRT-LBM by an external forcing term. In the procedure of the TRT-LBM-VP, the processes of interpolating velocities on the boundaries points and distributing the force density to the Eulerian points are unneeded. Performing the TRT-LBM-VP on a certain point, only the variables of this point are needed. As a consequence, the TRT-LBM-VP can be conducted parallelly. From the comparison between the result of the cylindrical Couette flow solved by the TRT-LBM-VP and that solved by the Single-relaxation-time LBM-VP (SRT-LBM-VP), the accuracy of the TRT-LBM-VP is higher than that of the SRT-LBM-VP. Flows past a single circular cylinder, a pair of cylinders in tandem and side-by-side arrangements, two counter-rotating cylinders and a NACA-0012 airfoil are chosen as numerical experiments to verify the present method further. Good agreements between the present results and those in the previous literatures are achieved.
Non-Newtonian fluid flows, especially in three dimensions (3D), arise in numerous settings of interest to physics. Prior studies using the lattice Boltzmann method (LBM) of such flows have so far been limited to mainly to two dimensions and used less robust collision models. In this paper, we develop a new 3D cascaded LBM based on central moments and multiple relaxation times on a three-dimensional, nineteen velocity (D3Q19) lattice for simulation of generalized Newtonian (power law) fluid flows. The relaxation times of the second order moments are varied locally based on the local shear rate and parameterized by the consistency coefficient and the power law index of the nonlinear constitutive relation of the power law fluid. Numerical validation study of the 3D cascaded LBM for various benchmark problems, including the complex 3D non-Newtonian flow in a cubic cavity at different Reynolds numbers and power law index magnitudes encompassing shear thinning and shear thickening fluids, are presented. Furthermore, numerical stability comparisons of the proposed advanced LBM scheme against the LBM based on other collision models, such as the SRT model and MRT model based on raw moments, are made. Numerical results demonstrate the accuracy, second order grid convergence and significant improvements in stability of the 3D cascaded LBM for simulation of 3D non-Newtonian flows of power law fluids.
Current multi-component, multiphase pseudo-potential lattice Boltzmann models have thermodynamic inconsistencies that prevent them to correctly predict the thermodynamic phase behavior of partially miscible multi-component mixtures, such as hydrocarbon mixtures. This paper identifies these inconsistencies and attempts to design a thermodynamically consistent multi-component, multiphase pseudo-potential lattice Boltzmann model that allows mass transfer across the phase interfaces and is capable to predict the phase behavior of typically partially miscible hydrocarbon mixtures. The designed model defines the total interaction force for the entire phase and split the force into individual components. Through a properly derived force split factor associated with the volatility of each component, the model can achieve precise thermodynamic consistency in multi-component hydrocarbon mixtures, which is described by the iso-fugacity rule.
This paper proposes an improved lattice Boltzmann scheme for incompressible axisymmetric flows. The scheme has the following features. First, it is still within the framework of the standard lattice Boltzmann method using the single-particle density distribution function and consistent with the philosophy of the lattice Boltzmann method. Second, the source term of the scheme is simple and contains no velocity gradient terms. Owing to this feature, the scheme is easy to implement. In addition, the singularity problem at the axis can be appropriately handled without affecting an important advantage of the lattice Boltzmann method: the easy treatment of boundary conditions. The scheme is tested by simulating Hagen-Poiseuille flow, three-dimensional Womersley flow, Wheeler benchmark problem in crystal growth, and lid-driven rotational flow in cylindrical cavities. It is found that the numerical results agree well with the analytical solutions and/or the results reported in previous studies.