Aspect ratio effects on three-dimensional incompressible flow in a two-sided non-facing lid-driven parallelepiped cavity

A numerical study of the three-dimensional fluid flow has been carried out to determine the effects of the transverse aspect ratio, Ay, on the flow structure in two-sided non-facing lid-driven cavities. The flow is complex, unstable and can undergo bifurcation. The numerical method is based on the finite volume method and multigrid acceleration. Computations have been investigated for several Reynolds numbers and various aspect ratio values. At a fixed Reynolds number, Re=500, the three-dimensional flow characteristics are analyzed considering four transverse aspect ratios, Ay=1,0.75,0.5 and 0.25. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. An analysis of the flow evolution shows that, at Ay=0.75, the flow bifurcates to a periodic regime at (Re=600) with a frequency f=0.093 less than the predicted value in the cubical cavity. A correlation is established when Ay=0.5 and gives the critical Reynolds number value. At Ay=0.25, the periodic regime occurs at high Re value beyond 3500, after which the flow becomes chaotic. It is shown that, when increasing Ay over the unit, the flow in the cavity exhibits a complex behavior. The kinetic energy transmission from the driven walls to the cavity center is reduced at low Ay values.     

Lattice Boltzmann simulation of natural convection in open end cavity with inclined hot wall

Natural convection in an open end cavity with a hot inclined wall is simulated based on the lattice Boltzmann method (LBM). The physics of flow and energy transfer in open end cavities are addressed when the hot wall is inclined. The combination of the two topics (open cavity and inclined walls) is the main novelty of the present study. The effects of the angle of the hot inclined wall on the flow field and heat transfer are thoroughly investigated. The Prandtl number is fixed to 0.71 (air). The Rayleigh number and the angle of the hot inclined wall are varied in the range of 10⁴ to 10⁶ and 60? to 85?, respectively. The results are presented for two different aspect ratios, i.e., A = 1 and 2. The results obtained with the LBM are also compared with those of the finite volume method (FVM). The predicted results of the LBM conform to those of the FVM. The results show that by increasing the angle of the hot inclined wall and the aspect ratio of the cavity, the average Nusselt number decreases. The trend of the local Nusselt number on the inclined wall is also discussed.     

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

Lattice Boltzmann simulation of MHD natural convection in a cavity with porous media and sinusoidal temperature distribution

The lattice Boltzmann method (LBM) is used to simulate the effect of magnetic field on the natural convection in a porous cavity. The sidewalls of the cavity are heated sinusoidally with a phase derivation, whereas the top and bottom walls are thermally insulated. Numerical simulation is performed, and the effects of the pertinent parameters, e.g., the Hartmann number, the porosity, the Darcy number, and the phase deviation, on the fluid flow and heat transfer are investigated. The results show that the heat transfer is affected by the temperature distribution on the sidewalls clearly. When the Hartmann number is 0, the maximum average Nusselt number is obtained at the phase deviation 90°. Moreover, the heat transfer enhances when the Darcy number and porosity increase, while decreases when the Hartman number increases.     

不同壁面绕流特性的格子Boltzmann模拟研究

为了探讨不同壁面的绕流特性,针对粘性流场中,不同壁面诱导的涡脱落现象以及升阻力系数等流场特性进行了格子Boltzmann数值研究。利用基于分子动理论的格子Boltzmann方法(LBM)求解Navier-Stokes方程,实现对流体运动的描述,针对不同的壁面条件,分别采用不同的格子Boltzmann流-固壁面处理方法。采用Half-way反弹边界条件来处理平直壁面,而曲壁面则采用LBM与有限差分法相结合的形式进行处理,计入了壁面与标准网格不重合对结果造成的影响。开发相应的计算程序,计算结果与已发表文献结果吻合良好,验证了数值模型的正确性。同时,探讨了进出口边界与钝体中心的距离对结果的影响。对比分析了不同壁面的绕流模型中升阻力系数、斯托罗哈数和涡量云图等,并进一步研究了雷诺数条件的影响。结果表明,不同壁面的绕流特性具有明显差异,且同时受雷诺数的显著影响;一般地,平直壁面对于来流作出的响应更迅速。     

用格子Boltzmann方法模拟Belousov-Zhabotinsky反应中的靶型波

构造了用于Belosov-Zhabotinsky反应的格子Boltzmann模型。通过对使用多组分的分布函数满足的格子Boltzmann方程,进行多重尺度Knudsen数展开,得到了模型的平衡态分布函数的各向同性解。作为算例,给出随机初始条件下反应区域内的靶型波的模拟结果,再现了Belousov-Zhabotinsky反应的经典结果。     

Lattice Boltzmann methods for shallow water flow applications

We apply the lattice Boltzmann (LB) method for solving the shallow water equations with source terms such as the bed slope and bed friction. Our aim is to use a simple and accurate representation of the source terms in order to simulate practical shallow water flows without relying on upwind discretization or Riemann problem solvers. We validate the algorithm on problems where analytical solutions are available. The numerical results are in good agreement with analytical solutions. Furthermore, we test the method on a practical problem by simulating mean flow in the Strait of Gibraltar. The main focus is to examine the performance of the LB method for complex geometries with irregular bathymetry. The results demonstrate its ability to capture the main flow features. Copyright © 2007 John Wiley & Sons, Ltd.     

脉冲流在Y型和倒Y型微反应器中强化传质的格子Boltzmann模拟

连续流微反应器的迅速发展为化学合成技术提供了一条可精准控制的路径。微反应器中流体的流动、混合和传质是反应的物理基础,因此强化传递就能够合理改善混和效果,增加相接触面积,减小微通道尺寸,以缩短分子扩散距离。由于微反应器中对流传递数量级非常低,不容易控制。通过改进微混合器构型和引入脉冲流动等主被动强化措施,可以有效改善对流和传质,影响化合反应。如何量化分析这些影响,是微反应器强化传质和优化控制反应过程的基础。本文基于有效的虚拟串行竞争反应格子Boltzmann模型,通过对Y型和倒Y型的微反应器的流场结构、混合传质和化合反应进行数值研究,定量分析了脉冲流动在不同流场构型下传质和化学反应的影响,所得结论可为连续流微反应器设计以及微反应精准控制提供有意义的参考。     

Lattice Boltzmann simulation of flow past a circular cylinder near a moving wall

The two‐dimensional flows past a circular cylinder near a moving wall are simulated by lattice Boltzmann method. The wall moves at the inlet velocity and the Reynolds number ranges from 300 to 500. The influence of the moving wall on the flow patterns is demonstrated and the corresponding mechanism is illustrated by using instability theory. The correlations among flow features based on gap ratio are interpreted. Force coefficients, velocity profile and vortex structure are analyzed to determine the critical gap ratio. Copyright © 2011 John Wiley & Sons, Ltd.     

The simulation of turbulent particle‐laden channel flow by the Lattice Boltzmann method

We perform direct numerical simulation of three‐dimensional turbulent flows in a rectangular channel, with a lattice Boltzmann method, efficiently implemented on heavily parallel general purpose graphical processor units. After validating the method for a single fluid, for standard boundary layer problems, we study changes in mean and turbulent properties of particle‐laden flows, as a function of particle size and concentration. The problem of physical interest for this application is the effect of water droplets on the turbulent properties of a high‐speed air flow, near a solid surface. To do so, we use a Lagrangian tracking approach for a large number of rigid spherical point particles, whose motion is forced by drag forces caused by the fluid flow; particle effects on the latter are in turn represented by distributed volume forces in the lattice Boltzmann method. Results suggest that, while mean flow properties are only slightly affected, unless a very large concentration of particles is used, the turbulent vortices present near the boundary are significantly damped and broken down by the turbulent motion of the heavy particles, and both turbulent Reynolds stresses and the production of turbulent kinetic energy are decreased because of the particle effects. We also find that the streamwise component of turbulent velocity fluctuations is increased, while the spanwise and wall‐normal components are decreased, as compared with the single fluid channel case. Additionally, the streamwise velocity of the carrier (air) phase is slightly reduced in the logarithmic boundary layer near the solid walls. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of laminar jet-forced flow using lattice Boltzmann method

In the paper, a numerical study on symmetrical and asymmetrical laminar jet-forced flows is carried out by using a lattice Boltzmann method （LBM） with a special boundary treatment. The simulation results are in very good agreement with the available numerical prediction. It is shown that the LBM is a competitive method for the laminar jet-forced flow in terms of computational efficiency and stability.     

Investigation of incompressible flow within ½ circular cavity using lattice Boltzmann method

A wall‐driven incompressible viscous flow in a ½ circular cavity is simulated, based on the lattice Boltzmann method (LBM). The treatment of curved boundary with second‐order accuracy is used. The force evaluation is based on the momentum‐exchange method. The streamlines and vorticity contours and the velocity component along the central line of a semi‐circular cavity are obtained for different Reynolds numbers. The numerical results show that the LBM can capture the formation of primary, secondary and tertiary vortices exactly as the Reynolds number increases and has a great agreement with those of current literatures. Copyright © 2008 John Wiley & Sons, Ltd.     

微通道人工泡状流的Lattice Boltzmann数值模拟

采用高频电控热激发汽泡的方式构造微通道人工泡状流,可以有效抑制微通道沸腾流动的不稳定性和强化传热。本文基于Lattice Boltzmann大密度比多相流复合模型,数值研究了通道内人工泡状流的流动和传热,通过比较分析不同发泡频率的泡状流,量化分析了汽泡运动和增长对微通道流动与传热的相互影响。一方面着重分析了汽泡运动对微通道运动边界层以及汽泡相变增长对热边界层的影响,另一方面也研究了边界层对汽泡动力行为的影响,所得结论对研究抑制微通道沸腾流动不稳定性和强化传热有参考意义。     

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

Simple lattice Boltzmann approach for turbulent buoyant flow simulation

In the present work, a simple large eddy simulation （LES）-based lattice Boltz- mann model （LBM） is developed for thermal turbulence research. This model is validated by some benchmark tests. The numerical results demonstrate the good performance of the present model for turbulent buoyant flow simulation.     

Lattice Boltzmann simulations for flow and heat/mass transfer problems in a three‐dimensional porous structure

The lattice Boltzmann method (LBM) for a binary miscible fluid mixture is applied to problems of transport phenomena in a three‐dimensional porous structure. Boundary conditions for the particle distribution function of a diffusing component are described in detail. Flow characteristics and concentration profiles of diffusing species at a pore scale in the structure are obtained at various Reynolds numbers. At high Reynolds numbers, the concentration profiles are highly affected by the flow convection and become completely different from those at low Reynolds numbers. The Sherwood numbers are calculated and compared in good agreement with available experimental data. The results indicate that the present method is useful for the investigation of transport phenomena in porous structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical study of Prandtl effect on MHD flow at a lid‐driven porous cavity

In this paper, the lattice Boltzmann method is used to study the Prandtl number effect on flow structure and heat transfer rates in a magnetohydrodynamic flow mixed convection in a lid‐driven cavity filled with a porous medium. The right and left walls are at constant but different temperatures (θ_h and θ_c), while the other walls are adiabatic. Gallium and salt water (0.02 < Pr < 13.4) are used as samples of the electroconducting fluids in the cavity. Typical sets of streamlines and isotherms are presented to analyze the flow patterns set up by the competition among the forced flow created by the lid‐driven wall, the buoyancy force of the fluid and the magnetic force of the applied magnetic field. Mathematical formulations in the porous media were constructed based on the Brinkman–Forchheimer model, while the multidistribution‐function model was used for the magnetic field effect. Numerical results were obtained and the effects of the Prandtl number and the other effective parameters such as Richardson, Hartman, and Darcy numbers were investigated. It was found that the fluid fluctuations within the cavity were reduced by increasing the Hartman number. A similar pattern was observed for the Darcy number reduction. Heat transfer was essentially dominated by the conduction for the low Prandtl number and forced convection dominated as the Prandtl number increased. Also, the average Nusselt number was raised by increasing the Prandtl number. It was discovered that a remarkable heat transfer enhancement of up to 28% could be reached by increasing the Prandtl number (from 0.02 to 13.4) at constant Richardson and Darcy numbers. Copyright © 2011 John Wiley & Sons, Ltd.     

Intercomparison of boundary schemes in Lattice Boltzmann method for flow simulation in porous media

The Lattice Boltzmann method has been widely adopted to simulate flow in porous media. The choice of appropriate boundary schemes is essential to achieve simulation accuracy; however, the criteria for the most suitable boundary treatment in the simulation of flow in porous media flow remain unresolved. Here, three types of the most commonly used boundary conditions are tested: interpolation bounce back (IBB), partial saturated method (PSM), and immersed boundary method (IBM). The dimensionless drag of face-centered cubic (FCC) sphere array and the dimensionless permeability of a random closely packed (RCP) sphere array are calculated and compared at different viscosities and resolutions. In the FCC sphere array case where spheres are not contacted, the IBB and PSM exhibit the same accuracy and both are of the second-order convergence rate. The IBM is less accurate and is of the first-order convergence rate. In the RCP sphere array case where the spheres are contacted, the IBB shows finer results and a second-order convergence rate. PSM underestimates the dimensionless permeability and increases resolution only slightly improved the results of PSM. The IBM overestimates the dimensionless permeability. These results indicate that among the three methods, the IBB is the most accurate. The PSM has the same accuracy as the IBB when sediments are not contacted; however, it loses its accuracy in the simulation of flow in closely packed porous media. This work could serve as a benchmark for further research in choosing the most appropriate method in the simulation of flow in porous media.     

A highly accurate GPU Lattice Boltzmann method with directional interpolation for the probability distribution functions

In this article, a highly accurate and graphics processing unit (GPU)-accelerated Lattice Boltzmann Method (LBM) is presented. The methodology is derived from a combination of conventional and recent LBM algorithms, mainly focusing on reducing the computational time, memory allocation, and complexity of existing algorithms. The general implementation focuses on accelerating the overall methodology using GPGPU technology based on Compute Shaders from OpenGL and avoids the storage of the distribution function components to reduce the memory allocation size. Furthermore, an efficient spatial interpolation of the probability distribution function components is described, based on a directional interpolation, without unnecessary control points for the reconstruction of virtual nodes data. The present methodology, tested for spatial accuracy via two- and three-dimensional Lid-Driven Cavity benchmark cases, shows excellent agreement with the results reported in the literature. Additionally, time efficiency is analyzed by comparing different configurations for the construction of virtual streaming points.