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.     

Simulation of natural convection melting in a cavity with fin using lattice Boltzmann method

In this study, a numerical investigation of melting phenomenon with natural convection in a cavity with fin has been performed using enthalpy‐based lattice Boltzmann method. The lattice D2Q9 model was applied to determine the density and velocity fields, and the D2Q5 model for the temperature field. The effect of vertical position and length of the fin on the melting rate was studied. The simulations were carried out for Stefan number of 10, Rayleigh number of 10 ⁵ and relative thermal conductivity (k_fin∕k_fluid) ranging from 5 to 30. The obtained results show that the rate of melting increases when the relative thermal conductivity and the length of the fin become greater. We also found that the variation of vertical position of the fin from bottom to middle has an insignificant effect on melting while it causes the increase of full melting time when the fin is mounted on the top of the cavity. Copyright © 2011 John Wiley & Sons, Ltd.     

Lattice Boltzmann simulation of two-sided lid-driven flow in a staggered cavity

In this article, the lattice Boltzmann method is employed in order to explore incompressible fluid flow inside a two-sided lid-driven staggered cavity. Results of the lattice Boltzmann simulation for antiparallel motion of lids are compared with the data from existing literature. For parallel motion of lids, the characteristics of flow pattern for a variety of Re numbers (50–3200) are presented. An asymmetric steady-state flow pattern for parallel motion of lids is obtained.     

Free convection of nanofluid filled enclosure using lattice Boltzmann method （LBM）

The lattice Boltzmann method (LBM) is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanoparticles: TiO₂, Ag, Cu, and Al₂O₃. The Brinkman and Maxwell-Garnetts models are used to simulate the viscosity and the effective thermal conductivity of nanofluids, respectively. Results from the performed numerical analysis show good agreement with those obtained from other numerical methods. A variety of the Rayleigh number, the nanoparticle volume fraction, and the aspect ratio are examined. According to the results, choosing copper as the nanoparticle leads to obtaining the highest enhancement for this problem. The results also indicate that the maximum value of enhancement occurs at λ = 2.5 when Ra = 10⁶ while at λ = 1.5 for other Rayleigh numbers.     

LES-based filter-matrix lattice Boltzmann model for simulating turbulent natural convection in a square cavity

In this paper, a novel thermal filter-matrix lattice Boltzmann model based on large eddy simulation (LES) is proposed for simulating turbulent natural convection. In this study, the Vreman subgrid-scale eddy-viscosity model is introduced into the present framework of LES to accurately predict the flow in near-wall region. Two dimensional numerical simulations of natural convection in a square cavity were performed at high Rayleigh number varying from 10⁷ to 10¹⁰ with a fixed Prandtl number of Pr = 0.71. The influences of the higher-order terms upon the present results at high Rayleigh numbers are examined, taking Ra = 10⁷ and 10⁸ as the example, revealing that the proper minimization of the higher-order terms can improve numerical accuracy of present model for high Rayleigh convective flow. For the turbulent convective flow, the time-averaged quantities in the median lines are presented and compared against those available results from previous studies. The general structure of turbulent boundary layers is well predicted. All numerical results exhibit good agreement with the benchmark solutions available in the previous literatures.     

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.     

矩形倾斜腔体中的两种对流斑图和分区

为了研究矩形倾斜腔体中普朗特数Pr=0.72的流体对流斑图和分区,本文基于流体力学方程组进行了数值模拟。在相对瑞利数r=6.0的情况下,观察了倾角θ=10°和θ=60°时对流斑图随着时间的发展,发现系统存在单圈型对流和多圈型对流两种斑图。流线随着倾角的变化说明:随着倾角增加,对流圈数逐渐减少,对流波长逐渐增加,对流波数减小;然后,随着对流圈数显著地减少,系统由多圈型对流过渡到单圈型对流。根据模拟计算结果,给出了多圈型对流过渡到单圈型对流的临界倾角θc随着相对瑞利数r变化的关系曲线。对流在θ-r平面上分为两个区域:θ<θc时系统是单圈型对流,θ>θc时系统是多圈型对流。对流最大振幅A和努塞尔数Nu随着倾角θ的变化曲线被临界倾角θc分成两段,它们有不同的变化规律。因此,临界倾角也可以由对流最大振幅A或努塞尔数Nu的变化曲线来确定。     

New studying of Lattice Boltzmann method for two-phase driven in porous media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｗｏ＿ｐｈａｓｅｄｒｉｖｅｎｉｎｐｏｒｏｕｓｍｅｄｉａｉｓａｎｉｍｐｏｒｔａｎｔｓｕｂｊｅｃｔｏｆｗｉｄｅｉｎｔｅｒｄｉｓｃｉｐｌｉｎａｒｙｃｏｎｃｅｒｎ ,ｓｕｃｈａｓｐｅｔｒｏｌｅｕｍｉｎｄｕｓｔｒｙ ,ｓｕｂｔｅｒｒａｎｅａｎ ,ｃｈｅｍｉｃａｌｐｒｏｊｅｃｔｅｔｃ..Ｉｔ’ｓｖｅｒｙｄｉｆｆｉｃｕｌｔｉｆｗｅａｄｏｐｔｎｕｍｅｒｉｃａｌｉｎｔｅｇｒａｌｍｅｔｈｏｄｆｏｒｔｈｅｃａｓｅｏｆｍｕｌｔｉｐｈａｓｅａｎｄｍｕｌｔｉｃｏｍｐｏｓｉｔｉｏｎ ,ｔｈｅｃｒｕｘｏｆ…     

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.     

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.     

Numerical investigation of natural convection inside an inclined parallel‐walled channel

Steady two‐dimensional natural convection in an inclined parallel‐walled channel was investigated numerically. The full elliptic forms of conservation equations were solved together and the velocity vectors, temperature contours and local and average Nusselt number distribution were obtained. The comparisons of local and average Nusselt number with published experimental and numerical results indicate very good agreement. Results are presented for a single aspect ratio, L/b=24, over the range of Rayleigh number of 3–1000 and angle of inclination 0–90°. The results indicate that the overall channel average Nusselt number is reduced as the inclination angle is increased. Significant reductions in the overall Nusselt number are exhibited at high angle of channel inclination. Copyright © 2005 John Wiley & Sons, Ltd.     

Lattice Boltzmann study of the interaction between a single solid particle and a thin liquid film

The isothermal single-component multi-phase lattice Boltzmann method(LBM) combined with the particle motion model is used to simulate the detailed process of liquid film rupture induced by a single spherical particle.The entire process of the liquid film rupture can be divided into two stages.In Stage 1,the particle contacts with the liquid film and moves into it due to the interfacial force and finally penetrates the liquid film.Then in Stage 2,the upper and lower liquid surfaces of the thin fi...     

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.     

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

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

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

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

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

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

Transient natural convection in a cavity with heat input and a constant temperature wall on opposite sides

The transient natural convection of a fluid with Prandtl number of order 200 in a two-dimensional square cavity has been numerically studied. One of the vertical walls of the cavity is kept at a constant (ambient) temperature and a constant heat flux is applied on the opposite wall. The other walls are adiabatic. Initially, a boundary layer is formed near the heated wall; subsequently, a large vortical structure is generated, together with an upper intrusion layer. As time progresses, the average temperature in the cavity increases, and a descending boundary layer is formed near the constant temperature wall. During the transition to the steady-state regime, a thermal stratification pattern is formed. The results are compared with the scale analysis presented by Patterson and Imberger (1980).     

Natural convection in an inclined quadrantal cavity heated and cooled on adjacent walls

Natural convective flow and heat transfer in an inclined quadrantal cavity is studied experimentally and numerically. The particle tracing method is used to visualize the fluid motion in the enclosure. Numerical solutions are obtained via a commercial CFD package, Fluent. The working fluid is distilled water. The effects of the inclination angle, ? and the Rayleigh number, Ra on fluid flow and heat transfer are investigated for the range of angle of inclination between 0° ? ? ? 360°, and Ra from 10⁵ to 10⁷. It is disclosed that heat transfer changes dramatically according to the inclination angle which affects convection currents inside, i.e. flow physics inside. A fairly good agreement is observed between the experimental and numerical results.     

Coupling of finite volume method and thermal lattice Boltzmann method and its application to natural convection

On the basis of the existing density distribution function reconstruction operator, the temperature distribution operator was derived to calculate heat transfer by coupling the lattice Boltzmann method (LBM) with the finite volume method. The present coupling model was validated by two‐dimensional natural convection flows with and without an isolated internal vertical plate. The results from the coupling model agree well with those from the pure‐finite volume method, pure‐LBM and references, and all the physical quantities cross the coupled interface smoothly. On the basis of residual history curves, it is likely that the convergence property and the numerical stability of the present model are better than those of the pure‐LBM at fine grid numbers and high Rayleigh numbers. Copyright © 2011 John Wiley & Sons, Ltd.