基于LBM的部分热活跃边界下多孔腔体内自然对流换热

采用Boltzmann方法模拟部分热活跃边界下的多孔腔体内自然对流,探讨不同热边界布置方案、孔隙度、Da数及Ra数对其流动传热的影响.数值计算表明：Da=10^-4时,腔体内中央出现一个循环流模式,只在Ra数很大时孔隙度才对传热有影响; Da=10^-2时,腔体内出现两个循环流,在Ra数很小时孔隙度对传热产生强烈的的影响.热活跃边界位置影响腔体内流体对流传热的强度,加热边界布置在底部、而冷却边界布置在顶部（Bottom-Top布置方式）,对多孔腔体内对流传热最有利,优于全热边界布置方式的传热效果.     

Analysis of transport of collimated radiation in a participating media using the lattice Boltzmann method

Application of the lattice Boltzmann method (LBM) recently proposed by Asinari et al. [Asinari P, Mishra SC, Borchiellini R. A lattice Boltzmann formulation to the analysis of radiative heat transfer problems in a participating medium. Numer Heat Transfer B 2010; 57:126–146] is extended to the analysis of transport of collimated radiation in a planar participating medium. To deal with azimuthally symmetric radiation in planar medium, a new lattice structure for the LBM is used. The transport of the collimated component in the medium is analysed by two different, viz., flux splitting and direct approaches. For different angles of incidence of the collimated radiation, the LBM formulation is tested for the effects of the extinction coefficient, the anisotropy factor, and the boundary emissivities on heat flux and emissive power distributions. Results are compared with the benchmark results obtained using the finite volume method. Both the approaches in LBM provide accurate results.     

Lattice Boltzmann method for nanofluid flow in a porous cavity with heat sources and magnetic field

In this article, Lattice Boltzmann method (LBM) has been applied to investigate the influences of magnetic field and heat sources on water based nanofluid natural convection inside a porous cavity with three square heat sources. Koo–Kleinstreuer–Li (KKL) model is applied to study Brownian motion impact on nanofluid flow. Effects of Rayleigh number (Ra), Darcy number (Da), nanofluid volume fraction (ϕ), and Hartmann number (Ha) on heat transfer characteristics are analyzed. From the obtained results we observe a decrease in the temperature gradient with increasing Ha; while quite the opposite effect is true with increasing Da and Ra. In the absence of magnetic field, for higher values of Darcy and Rayleigh numbers, thermal plumes are generated and the temperature gradient is enhanced. Moreover, small eddies are generated near the vertical centerline. However, in the presence of magnetic field, the number of thermal plumes decreases.     

Radiative and conductive heat transfer in a nongrey semitransparent medium. Application to fire protection curtains

This paper deals with heat transfer in nongrey media which scatter, absorb and emit radiation. Considering a two dimensional geometry, radiative and conductive phenomena through the medium have been taken into account. The radiative part of the problem was solved using the discrete ordinate method with classical S_n quadratures. The absorption and scattering coefficients involved in the radiative transfer equation (RTE) were obtained from the Mie theory. Conduction inside the medium was linked to the RTE through the energy conservation. Validation of the model has been achieved with several simulation of water spray curtains used as fire protection walls.     

Time-dependent radiative transfer using Chapman-Enskog maximum entropy method

The time-dependent problems of radiative transfer involve a coupling between radiation and material energy fields and are nonlinear because of proposed temperature dependence of the medium characteristics in semi-infinite medium with Rayleigh anisotropic scattering. By means of the limited flux, Chapman-Enskog and maximum entropy technique the time-dependent radiative transfer equation has been solved explicitly. The maximum entropy method is used to solve the resulting differential equation for radiative energy density. The calculations are carried out for temperature (normalized dimensionless) Θ(x,τ), radiative energy density and net flux with Rayleigh and anisotropic scattering for different space at different times.     

A Lattice Boltzmann Kinetic Model for Microflow and Heat Transfer

In this paper, we propose a lattice Boltzmann BGK model for simulation of micro flows with heat transfer based on kinetic theory and the thermal lattice Boltzmann method (He et al., J. Comp. Phys. 146:282, 1998). The relaxation times are redefined in terms of the Knudsen number and a diffuse scattering boundary condition (DSBC) is adopted to consider the velocity slip and temperature jump at wall boundaries. To check validity and potential of the present model in modelling the micro flows, two two-dimensional micro flows including thermal Couette flow and thermal developing channel flow are simulated and numerical results obtained compare well with previous studies of the direct simulation Monte Carlo (DSMC), molecular dynamics (MD) approaches and the Maxwell theoretical analysis     

Non-Fourier heat conduction with radiation in an absorbing, emitting, and isotropically scattering medium

This article numerically analyses the combined conductive and radiative heat transfer in an absorbing, emitting, and isotropically scattering medium. The non-Fourier heat conduction equation, which includes the time lag between heat flux and the temperature gradient, is used to model the conductive heat transfer in the medium. It predicts that a temperature disturbance will propagate as a wave at finite speed. The radiative heat transfer is solved using the P₃ approximation method. In addition, the MacCormack's explicit predictor-corrector scheme is used to solve the non-Fourier problem. The effects of radiation including single scattering albedo, conduction-to-radiation parameter, and optical thickness of the medium on the transient and steady state temperature distributions are investigated in detail. Analysis results indicate that the internal radiation in the medium significantly influences the wave nature. The thermal wave nature in the combined non-Fourier heat conduction with radiation is more obvious for large values of conduction-to-radiation parameter, small values of optical thickness and higher scattering medium. The results from non-Fourier-effect equation are also compared to those obtained from the Fourier equation. Non-Fourier effect becomes insignificant as either time increases or the effect of radiation increases.     

A finite element-spherical harmonics radiation transport model for photon migration in turbid media

In this paper, we solve the steady-state form of the Boltzmann transport equation in homogeneous and heterogeneous tissue-like media with a finite element-spherical harmonics (FE-P_N) radiation transport method. We compare FE-transport and diffusion solutions in terms of the ratio of absorption to reduced scattering coefficient, (μ_a/μ_s′) and the anisotropy factor g. Two different scattering phase function formulas are employed to model anisotropic scattering in the slab media with high g-value. Influence of void-like heterogeneities, and of their boundaries with the surrounding medium on the transport of photons are also examined.     

基于总能形式的耦合的双分布函数热晶格玻尔兹曼数值方法

双分布函数热晶格玻尔兹曼数值方法在微尺度热流动系统中得到广泛的应用. 本文基于晶格玻尔兹曼平衡分布函数低阶Hermite展开式, 创新性地提出了包含黏性热耗散和压缩功的耦合的双分布函数热晶格玻尔兹曼数值方法, 将能量场内温度的变化以动量源的形式引入晶格波尔兹曼动量演化方程, 实现了能量场与动量场之间的耦合. 研究了考虑黏性热耗散和压缩功的和不考虑的两种热自然对流模型, 重点分析了不同瑞利数和普朗特数下流场内的流动情况以及温度、速度和平均努赛尔数的变化趋势. 本文实验结果与文献结果一致, 验证了本文数值方法的可行性和准确性. 研究结果表明: 随着瑞利数和普朗特数的增大, 方腔内对流传热作用逐渐增强, 边界处形成明显的边界层; 考虑黏性热耗散和压缩功的模型对流作用相对增强, 黏性热耗散和压缩功对自然对流的影响在微尺度流动过程中不能忽略.     

利用格子Boltzmann方法模拟矩形腔内纳米流体Raleigh-Benard对流

利用格 子 Boltzmann方法模拟矩形腔内纳米流体Rayleigh-Benard对流, 得到温度场和流线分布, 比较分析不同Ra数、体积分数、粒径下纳米流体对流换热的变化情况. 结果表明: 在相同的Ra 数和体积分数下, 纳米流体的对流换热随着粒径的增大而减弱; 在相同的Ra数和粒径下, 纳米流体的对流换热随着体积分数增大而增强. 关键词： 纳米流体 Raleigh-Benard 多相流 格子Boltzmann方法     

Lattice Boltzmann simulation of 3-dimensional natural convection heat transfer of CuO/water nanofluids

The present study investigated fluid flow and natural convection heat transfer in an enclosure embedded with isothermal cylinder. The purpose was to simulate the three-dimensional natural convection by thermal lattice Boltzmann method based on the D3Q19 model. The effects of suspended nanoparticles on the fluid flow and heat transfer analysis have been investigated for different parameters such as particle volume fraction, particle diameters, and geometry aspect ratio. It is seen that flow behaviors and the average rate of heat transfer in terms of the Nusselt number (Nu) are effectively changed with different controlling parameters such as particle volume fraction (5 % ≤ φ ≤ 10 %), particle diameter (d _p = 10 nm to 30 nm) and aspect ratio (0.5 ≤ AR ≤ 2) with fixed Rayleigh number, Ra = 10⁵. The present results give a good approximation for choosing an effective parameter to design a thermal system.     

Double MRT thermal lattice Boltzmann method for simulating convective flows

A two-dimensional double Multiple Relaxation Time-Thermal Lattice Boltzmann Equation (2-MRT-TLBE) method is developed for predicting convective flows in a square differentially heated cavity filled with air (Pr=0.71). In this Letter, we propose a numerical scheme to solve the flow and the temperature fields using the MRT-D2Q9 model and the MRT-D2Q5 model, respectively. Thus, the main objective of this study is to show the effectiveness of such model to predict thermodynamics for heat transfer. This model is validated by the numerical simulations of the 2-D convective square cavity flow. Excellent agreements are obtained between numerical predictions. These results demonstrate the accuracy and the effectiveness of the proposed methodology.     

Mesoscopic analysis of heatline and massline during double-diffusive MHD natural convection in an inclined cavity

Double-diffusive natural convection in an inclined cavity with the presence of magnetic field is studied numerically via heatline and massline approach. The governing equations are discretized using the Lattice Boltzmann Method (LBM). In this investigation, the controlling parameters involved are Rayleigh number (10³ ≤ Ra ≤ 10⁵), buoyancy ratio (−5 ≤ N ≤ 5), cavity inclination angle (0° ≤ Ø ≤ 180°), Lewis number (2 ≤ Le ≤ 10), Prandtl number (Pr = 5.0) and Hartmann number (0 ≤ Ha ≤ 50). The numerical results obtained by the effect of parameters mentioned above are reported as contours of streamlines, isotherms, isoconcentrations, heatlines, and masslines. The obtained results are compared with the existing literature to validate the coding. The heat and mass transfer rate decrease with increasing the magnetic field and increase with an increase in Rayleigh number. The impact of Ø is maximum for higher Ra and negligible for lower Ra (10³). Increasing Le influences the mass transfer rate to increase and heat transfer to reduce for both opposing and aiding flow.     

Lattice Boltzmann simulation for temperature-sensitive magnetic fluids in a porous square cavity

A lattice Boltzmann method is developed to simulate temperature-sensitive magnetic fluids in a porous cavity. In the simulation, the magnetic force, efficient gravity, viscous loss term and geometric loss term in porous medium are imported to the momentum equation. To test the reliability of the method, a validation with water in porous cavity is carried out. Good agreements with the previous results verify that the present lattice Boltzmann method is promising for simulation of magnetic fluids in porous medium. In this study, we investigate the change of magnetization with external magnetic field, and we present numerical results for the streamlines, isotherms, and magnetization at vertical or horizontal mid-profiles for different values of Ram. In addition, Nusselt numbers changing with magnetic Rayleigh numbers are also investigated.     

Surface radiation influence on the regimes of conjugate natural convection in an enclosure with local energy source

A numerical analysis of the unsteady regimes of the natural convection and thermal radiation in a square enclosure with heat-conducting walls has been carried out in the presence of a heat source of finite sizes located in the base zone under the conditions of convective exchange with the ambient medium. The mathematical model formulated in the dimensionless variables “stream function — vorticity — temperature” has been implemented numerically by a finite difference method. The influence of the reduced emissivity of internal surfaces of bounding walls on the local characteristics (the streamlines and temperature fields) and on the integral complex (the mean Nusselt number on typical boundaries) has been analyzed in detail for 0 ≤ ? < 1, the location of the heat source 0.1 ≤ l/L ≤ 0.4, and its length 0.2 ≤ l _hs /L ≤ 0.6 for Ra = 10⁶, Pr = 0.7. The approximation relations have been derived for the mean convective and radiation Nusselt numbers depending on the reduced emissivity of the internal surfaces of bounding walls and the energy source location relative to vertical walls.     

Transient,combined conduction and radiation in an absorbing,emitting, and isotropically scattering solid sphere

Transient, combined conduction and radiation is solved in an absorbing, emitting, and isotropically scattering solid sphere with a black boundary initially at a uniform temperature and for times t > 0 subjected to a constant temperature at the spherical surface. The collocation method is used to solve the radiation part of this problem and the implicit finite difference scheme is used to solve the conduction part. The effects of the conduction-to-radiation parameter, the single scattering albedo, the optical thickness of the medium on the temperature distribution and the heat flux in the medium are examined.     

Étude numérique de l'influence du transfert radiatif sur la convection naturelle laminaire,bidimensionnelle, permanente,dans un espace annulaire d'axe horizontal,délimité par deux cylindres circulaires isothermes

Combined radiation and natural convection in a participating medium between concentric or vertically eccentric horizontal cylinders is investigated numerically. The annular medium is considered as a gray, emitting, absorbing, and isotropically scattering gas. The equations of steady, laminar, two-dimensional, thermal, natural convection are written by using a two-cylindrical coordinate system, the stream function, and the vorticity. The finite volumes method is used to discretize the coupled equations of momentum, energy, and radiative transfer. To solve the global nonlinear algebraic equations the successive-over-relaxation iterating scheme is applied. Numerical solutions are obtained for a Rayleigh number in the range 10³–10⁵ and radiation-conduction parameter ranging from 0 to ∞. The influences of radiation-conduction parameter, Rayleigh number and other parameters on flow and temperature distributions and heat transfer are discussed.     

An inverse analysis of a transient 2-D conduction-radiation problem using the lattice Boltzmann method and the finite volume method coupled with the genetic algorithm

This article deals with the simultaneous estimation of parameters in a 2-D transient conduction-radiation heat transfer problem. The homogeneous medium is assumed to be absorbing, emitting and scattering. The boundaries of the enclosure are diffuse gray. Three parameters, viz. the scattering albedo, the conduction-radiation parameter and the boundary emissivity, are simultaneously estimated by the inverse method involving the lattice Boltzmann method (LBM) and the finite volume method (FVM) in conjunction with the genetic algorithm (GA). In the direct method, the FVM is used for computing the radiative information while the LBM is used to solve the energy equation. The temperature field obtained in the direct method is used in the inverse method for simultaneous estimation of unknown parameters using the LBM-FVM and the GA. The LBM-FVM-GA combination has been found to accurately predict the unknown parameters.     

Radiative transfer in one-dimensional hollow sphere with anisotropic scattering and variable medium properties

The effects of variable medium properties on radiation transfer in participating and anisotropically scattering one-dimensional spherical medium were investigated by Tsai et al. (JQSRT 42(3) (1989) 187). The discrete ordinates method solutions they provided for hollow spherical medium cases are incorrect. The correct DOM S₈ and the integral transfer equation solutions are provided.     

Natural Convection of Temperature-Sensitive Magnetic Fluids in Porous Media

In this article, natural convection of a temperature-sensitive magnetic fluid in a porous media is studied numerically by using lattice Boltzmann method. Results show that the heat transfer decreases when the ball numbers increase. When the magnetic field is increased, the heat transfer is enhanced; however, the average wall Nusselt number increases at small ball numbers but decreases at large ball numbers due to the induced flow being more likely confined near the bottom walls with a high number of obstacles.