Revisiting the Lattice Boltzmann Method Through a Nonequilibrium Thermodynamics Perspective

2021 ◽  
Vol 143 (5) ◽  
Author(s):  
Anirudh Jonnalagadda ◽  
Atul Sharma ◽  
Amit Agrawal

Abstract In this paper, we incorporate a nonequilibrium thermodynamics perspective that is consistent with the Onsager reciprocity principle into the lattice Boltzmann framework to propose a novel regularized lattice Boltzmann formulation for modeling the Navier–Stokes–Fourier equations. The new method is applied to one-dimensional (1D) isothermal situations wherein the advantages of incorporating such a nonequilibrium perspective can be explicitly appreciated. In such situations, the nonequilibrium contribution of the lattice populations obtained by the new method completely vanishes, and the lattice update is entirely reduced to evaluating the equilibrium distribution function. Such a counterintuitive 1D mesoscopic description is not obtained in any other existing lattice Boltzmann scheme. We therefore numerically test the proposed formulation on two complex problems, namely, shockwave and nonlinear wave propagation, and compare results with analytical results along with six existing lattice Boltzmann schemes; it is found that the new method indeed yields results that are more stable and accurate. These results highlight the potency of the nonequilibrium thermodynamics-based approach for obtaining accurate and stable lattice Boltzmann computations, and provide new insights into established lattice Boltzmann simulation methods.

2016 ◽  
Vol 9 (1) ◽  
pp. 1-22 ◽  
Author(s):  
Z. Chen ◽  
C. Shu ◽  
Y. Wang ◽  
L. M. Yang ◽  
D. Tan

AbstractIn this paper, a simplified lattice Boltzmann method (SLBM) without evolution of the distribution function is developed for simulating incompressible viscous flows. This method is developed from the application of fractional step technique to the macroscopic Navier-Stokes (N-S) equations recovered from lattice Boltzmann equation by using Chapman-Enskog expansion analysis. In SLBM, the equilibrium distribution function is calculated from the macroscopic variables, while the non-equilibrium distribution function is simply evaluated from the difference of two equilibrium distribution functions. Therefore, SLBM tracks the evolution of the macroscopic variables rather than the distribution function. As a result, lower virtual memories are required and physical boundary conditions could be directly implemented. Through numerical test at high Reynolds number, the method shows very nice performance in numerical stability. An accuracy test for the 2D Taylor-Green flow shows that SLBM has the second-order of accuracy in space. More benchmark tests, including the Couette flow, the Poiseuille flow as well as the 2D lid-driven cavity flow, are conducted to further validate the present method; and the simulation results are in good agreement with available data in literatures.


1998 ◽  
Vol 09 (08) ◽  
pp. 1383-1391 ◽  
Author(s):  
Yu Chen ◽  
Shulong Teng ◽  
Takauki Shukuwa ◽  
Hirotada Ohashi

A model with a volumetric stress tensor added to the Navier–Stokes Equation is used to study two-phase fluid flows. The implementation of such an interface model into the lattice-Boltzmann equation is derived from the continuous Boltzmann BGK equation with an external force term, by using the discrete coordinate method. Numerical simulations are carried out for phase separation and "dam breaking" phenomena.


Author(s):  
Joris C. G. Verschaeve

By means of the continuity equation of the incompressible Navier–Stokes equations, additional physical arguments for the derivation of a formulation of the no-slip boundary condition for the lattice Boltzmann method for straight walls at rest are obtained. This leads to a boundary condition that is second-order accurate with respect to the grid spacing and conserves mass. In addition, the boundary condition is stable for relaxation frequencies close to two.


2014 ◽  
Vol 554 ◽  
pp. 665-669
Author(s):  
Leila Jahanshaloo ◽  
Nor Azwadi Che Sidik

The Lattice Boltzmann Method (LBM) is a potent numerical technique based on kinetic theory, which has been effectively employed in various complicated physical, chemical and fluid mechanics problems. In this paper multi-relaxation lattice Boltzmann model (MRT) coupled with a Large Eddy Simulation (LES) and the equation are applied for driven cavity flow at different Reynolds number (1000-10000) and the results are compared with the previous published papers which solve the Navier stokes equation directly. The comparisons between the simulated results show that the lattice Boltzmann method has the capacity to solve the complex flows with reasonable accuracy and reliability. Keywords: Two-dimensional flows, Lattice Boltzmann method, Turbulent flow, MRT, LES.


2017 ◽  
Vol 2017 ◽  
pp. 1-11 ◽  
Author(s):  
L. El Moutaouakil ◽  
Z. Zrikem ◽  
A. Abdelbaki

Laminar natural convection in a water filled square enclosure containing at its center a horizontal hexagonal cylinder is studied by the lattice Boltzmann method. The hexagonal cylinder is heated while the walls of the cavity are maintained at the same cold temperature. Two orientations are treated, corresponding to two opposite sides of the hexagonal cross-section which are horizontal (case I) or vertical (case II). For each case, the results are presented in terms of streamlines, isotherms, local and average convective heat transfers as a function of the dimensionless size of the hexagonal cylinder cross-section (0.1≤B≤0.4), and the Rayleigh number (103≤Ra≤106).


Author(s):  
R. Kamali ◽  
A. H. Tabatabaee Frad

It is known that the Lattice Boltzmann Method is not very effective when it is being used for the high speed compressible viscous flows; especially complex fluid flows around bodies. Different reasons have been reported for this unsuccessfulness; Lacking in required isotropy in the employed lattices and the restriction of having low Mach number in Taylor expansion of the Maxwell Boltzmann distribution as the equilibrium distribution function, might be mentioned as the most important ones. In present study, a new numerical method based on Li et al. scheme is introduced which enables the Lattice BoltzmannMethod to stably simulate the complex flows around a 2D circular cylinder. Furthermore, more stable implementation of boundary conditions in Lattice Boltzmann method is discussed.


2015 ◽  
pp. 157-168
Author(s):  
Natasa Lukic ◽  
Predrag Tekic ◽  
Jelena Radjenovic ◽  
Ivana Sijacki

The present study is concerned with two-sided lid-driven incompressible flow in rectangular, deep cavities applying lattice Boltzmann method. After validating the code for the square cavity, solutions for cavities with an aspect ratio 1.5 and 4 were obtained for the Reynolds numbers of 100, 400, 1000 and 3200. The influence of the Reynolds number and aspect ratio on the flow pattern and on the characteristics of vortices inside the cavity was studied. Symmetric flow pattern was obtained for all investigated cases. The middle of the cavity is mostly influenced by the increase in the aspect ratio. Critical aspect ratio, at which the birth of a primary vortex in the middle of the cavity takes place, was determined to be between 2.7 and 2.725.


2021 ◽  
Vol 321 ◽  
pp. 01014
Author(s):  
Makoto Sugimoto ◽  
Tatsuya Miyazaki ◽  
Zelin Li ◽  
Masayuki Kaneda ◽  
Kazuhiko Suga

Stator coils of automobiles in operation generate heat and are cooled by a coolant poured from above. Since the behavior characteristic of the coolant poured on the coils is not clarified yet due to its complexity, the three-dimensional two-phase flow simulation is conducted. In this study, as a steppingstone to the simulation of the liquid falling on the actual coils, the coils are modelled with horizontal rectangular pillar arrays whose governing parameters can be easily changed. The two-phase flows are simulated using the lattice Boltzmann method and the phase-field model, and the effects of the governing parameters, such as the physical properties of the cooling liquid, the wettability, and the gap between the pillars, on the wetting area are investigated. The results show that the oil tends to spread across the pillars because of its high viscosity. Moreover, the liquid spreads quickly when the contact angle is small. In the case that the pillars are stacked, the wetting area of the inner pillars is larger than that of the exposed pillars.


Author(s):  
Felipe A. Valenzuela ◽  
Amador M. Guzmán ◽  
Andrés J. Díaz

During the last years the aerodynamics characteristics of airfoils have been studied solving numerically the Navier-Stokes (NS) equations. These calculations require a significant computational cost due to both the second order and the nonlinear characteristics of the NS partial differential equations. Therefore, efforts have been devoted to reduce this cost and increase the accuracy of the numerical methods. The Lattice-Boltzmann Method (LBM) has become a great alternative to simulate this problem and a variety of fluid flows. In this method, the convective operator is linear and the pressure is calculated directly by the equation of state without implementing iterative methods. This work represents a preliminary investigation of a laminar flow over airfoils under low Reynolds number conditions (Re = 500). Solutions are obtained using a Multi-Block mesh refinement method. In order to validate the computational code, calculations are performed on a SD7003 airfoil at an angle of attack of 4° and 30°, which corresponds to the available numerical and experimental results. The results of this study agree well with previous experimental and numerical studies demonstrating the capabilities of the LBM to simulate accurately laminar flows over airfoils as well as capturing and predicting the laminar separation bubbles.


Author(s):  
Minglei Shan ◽  
Yu Yang ◽  
Hao Peng ◽  
Qingbang Han ◽  
Changping Zhu

Understanding the dynamic characteristic of the cavitation bubble near a solid wall is a fundamental issue for the bubble collapse application and prevention. In the present work, an improved three-dimensional multi-relaxation-time pseudopotential lattice Boltzmann model is adopted to investigate the cavitation bubble collapse near the solid wall. With respect to thermodynamic consistency, Laplace law verification, the three-dimensional pseudopotential multi-relaxation-time lattice Boltzmann model is investigated. By the theoretical analysis, it is proved that the model can be regarded as a solver of the Rayleigh–Plesset equation, and confirmed by comparing the results of the lattice Boltzmann simulation and the Rayleigh–Plesset equation calculation for the case of cavitation bubble collapse in the infinite medium field. The bubble collapse near the solid wall is modeled using the improved pseudopotential multi-relaxation-time lattice Boltzmann model. We find the lattice Boltzmann simulation and the experimental results have the same dynamic process by comparing the bubble profiles evolution. Form the pressure field and the velocity field evolution it is found that the tapered higher pressure region formed near the top of the bubble is a crucial driving force inducing the bubble collapse. This exploratory research demonstrates that the lattice Boltzmann method is an alternative tool for the study of the interaction between collapsing cavitation bubble and matter.


Sign in / Sign up

Export Citation Format

Share Document