浮力对混合对流流动及换热特性的影响

用热线和冷线相结合的技术测量垂直圆管内逆混合对流流体的平均速度、 温度以及它们的脉动. 较详细地研究了浮力对逆混合对流的流动特性和传热特性的影响. 评 估了实验中采用的冷线测量温度补偿速度探头温度敏感的影响. 逆混合对流的传热结果用无 量纲参数Ω (Ω= Grd / Red2 )来表示，其中，基于管道直 径的雷诺数Red变化范围为900~18000, 浮力参数Ω变化范围为 0.004899~0.5047. 研究结果表明，浮力对逆混合对流的换热有强化作用. 随着葛拉晓夫数Grd的增加，温度脉动，流向雷诺正应力和流向温度通量增 大，并且在靠近壁面的流体区域尤其明显. 热线与冷线相结合的技术适合于研究非绝热的流 动测量，可以用于研究浮力对流动和换热特性的影响.     

Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method

An efficient immersed boundary-lattice Boltzmann method (IB-LBM) is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows. Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor. A boundary thickening-based direct forcing IB method is adopted to solve the particle–fluid interactions with high accuracy for non-slip boundary conditions. A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computational parameters. An asynchronous calculation strategy is adopted to further improve the computing efficiency. The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods. The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed. Subsequently, 504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh–Taylor instabilities was observed. At last, simulations of a dense flow involving 11001 particles, the largest number of particles to date, were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid. The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time. These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.     

壁面约束对柱状粒子在牛顿流体中沉降影响的研究

格子Boltzmann方法(LBM)直接数值模拟了在有壁面约束的流场中柱状粒子的沉降。结果说明粒子在壁面的影响下，会离开管壁向管中心移动，移动速度先是增大，然后随着与管壁距离的增加，逐渐减小，粒子中心最后并不一定停留在管中心。随着粒子与壁面距离的增加，粒子的沉降速度将增大，可见壁面对粒子沉降有阻碍作用，平行于粒子轴线的壁面对粒子沉降的影响比垂直于粒子轴线壁面的影响大。     

A numerical study of micropolar flow inside a lid-driven triangular enclosure

A study is carried out to analyze the mixed convection flow and heat transfer inside a lid-driven triangular conduit under the effects of micro-gyration boundary conditions. The micropolar constitutive equation characterizes the fluid inside the cavity. The lower boundary is at a uniform temperature and sliding in its plane with constant velocity u₀, while the inclined walls are cold. Dual cases are considered here, namely the intense concentration (d) and the weak concentration of microelements (\(m = 0.5\)). The governing nonlinear equations are simulated employing the Galerkin finite element method, where the pressure term is handled via the Penalty approach. Using the numerical data, graphical results are produced to illustrate the effects of physical parameters. Specifically, this refers to the effects of the Grashof number (Gr), Prandtl number (Pr), Reynolds number (Re) and vortex viscosity parameter (K) on the streamlines, mid-section velocity profiles, temperature contours, and local and average Nusselt numbers on the cold and heated boundaries of the conduit. Particular emphasis is given on the identification of the set of parameters for which simultaneous symmetry in streamlines and isotherms prevails. The grid independence test is also performed by comparing the average Nusselt numbers (on the hot and cold boundaries of the conduit) for various mesh sizes, and the optimal solution is found. Moreover, the results are also benchmarked with the previously published data.     

Non-Darcy effects on convective boundary layer flow past a semi-infinite vertical plate in saturated porous media

The composite effects of viscosity, porosity, buoyancy parameter, thermal conductivity ratio and non-Darcy effects of Brinkman friction and Forscheimmer quadratic drag on the mixed convection boundary layer flow past a semi-infinite plate in a fully-saturated porous regime are theoretically and numerically investigated using Keller’s implicit finite-difference technique and a double-shooting Runge-Kutta method. The Brinkman Forcheimer-extended Darcy model is implemented in the hydrodynamic boundary layer equation. The effects of the various non-dimensional thermofluid parameters, viz Grashof number, Darcy number, and Forchheimer number, and also porosity, thermal conductivity and viscosity parameters on the velocity and temperature fields are discussed. Computations for both numerical schemes are made where possible and found to be in excellent agreement.     

Effect of the coflow stream on a plane wall jet

We propose in this work to study an isothermal and a non-isothermal laminar plane wall jet emerging in a coflow steam. The numerical solution of the governing equations was performed by a finite difference method. In this work, we are interested in the study of the influence of Grashof numbers on the wall jet emerging in a medium at rest. Further, we will examine the effect of the coflow stream on the behavior of the dynamic and thermal properties of the wall jet subjected to a constant temperature. A comparison with a simple wall jet is carried out. The results show that for a buoyant wall jet, two parameters can influence the flow: the inertial and buoyancy forces. The velocity effect indicates that the potential core length increases with the velocity ratio. We are also showed that when using a momentum length scale, the normalized longitudinal maximum velocity can reach an asymptotic curve at different velocity ratios.     

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 particle sedimentation in 3D rectangular channel

The 3D lattice Boltzmann method is used to simulate particle sedimentation in a rectangular channel. The results of single particle sedimentation indicate that the last position of the particle is along the center line of the channel regardless of the initial position, the particle diameter, and the particle Reynolds number. The wall effect on the terminal velocity is in good agreement with experimental results quantitatively. The drafting, kissing, and tumbling (DKT) process is reproduced and analyzed by simulating two-particle cluster sedimentation. The effects of the diameter ratio, initial position, and wall on the DKT process are investigated. When the two particles have equal diameter sediment in the rectangular channel, a periodical DKT process and the spiraling trajectory are found. The last equilibrium configuration is obtained from the simulation results. The interesting regular sedimentation phenomena are found when 49 particles fall down under gravity.     

Lower bound for transient growth of inclined buoyancy layer

The relationship between stabilities of the buoyancy boundary layers along an inclined plate and a vertical plate immersed in a stratified medium is studied theoretically and numerically. The eigenvalue problem of energy stability is solved with the method of descending exponentials. The disturbance energy is found to be able to grow to 11.62 times as large as the initial disturbance energy for P r = 0.72 when the Grashof number is between the critical Grashof numbers of the energy stability and the linear stability.We prove that, with a weighted energy method, the basic flow of the vertical buoyancy boundary layer is stable to finite-amplitude streamwise-independent disturbances.     

Direct-forcing immersed boundary lattice Boltzmann simulation of particle/fluid interactions for spherical and non-spherical particles

The lattice Boltzmann method (LBM) is a useful technique for simulating multiphase flows and modeling complex physics. Specifically, we use LBM combined with a direct-forcing (DF) immersed boundary (IB) method to simulate fluid–particle interactions in two-phase particulate flows. Two grids are used in the simulation: a fixed uniform Eulerian grid for the fluid phase and a Lagrangian grid that is attached to and moves with the immersed particles. Forces are calculated at each Lagrangian point. To exchange numerical information between the two grids, discrete delta functions are used. The resulting DF IB-LBM approach is then successfully applied to a variety of reference flows, namely the sedimentation of one and two circular particles in a vertical channel, the sedimentation of one or two spheres in an enclosure, and a neutrally buoyant prolate spheroid in a Couette flow. This last application proves that the developed approach can be used also for non-spherical particles. The three forcing schemes and the different factors affecting the simulation (added mass effect, corrected radius) are also discussed.     

Stratified magnetohydrodynamic flow of tangent hyperbolic nanofluid induced by inclined sheet

This paper studies stratified magnetohydrodynamic(MHD) flow of tangent hyperbolic nanofluid past an inclined exponentially stretching surface. The flow is subjected to velocity, thermal, and solutal boundary conditions. The partial differential systems are reduced to ordinary differential systems using appropriate transformations.The reduced systems are solved for convergent series solutions. The velocity, temperature,and concentration fields are discussed for different physical parameters. The results indicate that the temperature and the thermal boundary layer thickness increase noticeably for large values of Brownian motion and thermophoresis effects. It is also observed that the buoyancy parameter strengthens the velocity field, showing a decreasing behavior of temperature and nanoparticle volume fraction profiles.     

A direct‐forcing pressure‐based lattice Boltzmann method for solving fluid–particle interaction problems

A direct‐forcing immersed boundary‐lattice Boltzmann method (IB–LBM) is developed to simulate fluid–particle interaction problems. This method uses the pressure‐based LBM to solve the incompressible flow field and the immersed boundary method to handle the fluid–particle interactions. The pressure‐based LBM uses the pressure distribution functions instead of the density distribution functions as the independent dynamic variables. The main idea is to explicitly eliminate the compressible effect due to the density fluctuation. In the IB method, a direct‐forcing method is introduced to capture the particle motion. It directly computes an IB force density at each lattice grid from the differences between the pressure distribution functions obtained by the LBM and the equilibrium pressure distribution functions computed from the particle velocity. By applying this direct‐forcing method, the IB–LBM becomes a purely LBM version. Also, by applying the Gauss theorem, the formulas for computing the force and the torque acting on the particle from the flows are derived from the volume integrals over the particle volume instead of from the surface integrals over the particle surface. The order of accuracy of the IB–LBM is demonstrated on the errors of velocity field, wall stress, and gradients of velocity and pressure. As a demonstration of the efficiency and capabilities of the new method, sedimentation of a large number of spherical particles in an enclosure is simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

侧加热腔体内重力波演化过程的数值模拟

利用二维数值模拟的方法研究了侧加热腔体内的自然对流．基于数值模拟结果,描述了水平热入侵流(intrusion)的整个演化过程,并对该过程的物理机制进行了讨论．结果表明：当热入侵流抵达腔体冷壁后,由于冷壁无法卷入所有的热入侵流,热入侵流在冷上角堆积并产生一个反向流动,在冷壁边界层附近形成一个顺时针涡,该涡在浮力效应驱动下可返回热壁,并在腔体的冷热壁之间形成了腔体尺度的流体振荡,即内重力波．     

Exact solution of the kinetic equation in the problem of nonisothermal flow in the neighborhood of a weakly curved surface

An analytic method for solving the half-space boundary value problem for the inhomogeneous Boltzmann equation with the collision operator in the form of an elliptico-statistical model (the ES-model of the Boltzmann equation) is proposed for the problem of nonisothermal rarefied gas flow in the neighborhood of a curved surface. An exact analytic expression is derived for the thermal slip of a monatomic gas along the surface of a rigid spherical aerosol particle. A numerical value of the gas-kinetic coefficient which takes into account the effect of the curvature of the surface on the thermal slip coefficient is obtained. A comparison with published data is carried out. Moscow, Arkhangelsk. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 165–173, March–April, 1998.     

Particle flow simulations with homogenised lattice Boltzmann methods

An alternative approach to simulating arbitrarily shaped particles submersed in viscous fluid in two dimensions is proposed, obtained by adapting the velocity parameter of the equilibrium distribution function of a standard lattice Boltzmann method (LBM). Comparisons of exemplifying simulations to results in the literature validate the approach as well as the convergence analysis. Pressure fluctuations occurring in Ladd’s approach are greatly reduced. In comparison with the immersed boundary method, this approach does not require cost intensive interpolations. The parallel efficiency of LBM is retained. An intrinsic momentum transfer is observed during particle–particle collisions. To demonstrate the capabilities of the approach, sedimentation of particles of several shapes is simulated despite omitting an explicit particle collision model.     

Transient natural convection along a vertical clyinder with variable surface temperature and mass diffusion

A numerical study is carried out for thermal and concentration driven transient natural convection adjacent to a vertical cylinder. The temperature and concentration level at the cylinder surface are assumed to vary as power-law type functions, with exponents n and m respectively in the streamwise co-ordinate. The governing boundary layer equations are converted into a non-dimensional form. A Crank-Nicolson type of implicit finite-difference method is used to solve the governing non-linear set of equations. Numerical results are obtained and presented with various thermal and mass Grashof numbers and power law variations. Transient effects of velocity, temperature and concentration are analyzed. Local and average skin-friction, Nusselt number and Sherwood number are shown graphically.     

Kinetic theory of interacting particles in dilute suspensions of mineral waste

A kinetic theory of interacting spherical particles in dilute suspension is developed which results in a Boltzmann transport equation. This equation is solved in the relaxation time approximation to calculate the settling velocity of fine particulates in the steady state. The theory is applied to the suspended Jamaican bauxite waste and kaolinite particles. The experimental settling velocity compares well with the calculated values at low concentrations. This treatment can form the basis for a more rigorous theory applicable to denser systems and non-spherical particulates.     

幂律流体圆柱绕流的格子波尔兹曼模拟简

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明：插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.     

Study on the interaction of sedimenting cylindrical particles in still fluid

The sedimentations of two cylindrical particles in three different initial relative positions are numerically simulated using the lattice Boltzmann method. The movement characteristics and particle interactions during their sedimentation are presented and discussed in detail. The results show that, (i) if the two particles are released parallel but separated horizontally, they push away each other, rotate inwards and separate horizontally as they fall; (ii) if the two particles are released parallel but separated vertically, the sedimentation behavior can be classified into three stages: trailing, tumbling and separating; (iii) if the two particles are released perpendicular but separated vertically, the sedimentation behavior can be characterized as: trailing and rotating, touching and sliding. In order to validate our simulation, experiments were also conducted and the results agree well with the numerical ones. The project supported by the National Natural Science Foundation of China for Distinguished Scholars (19925210)