格子Boltzmann方法分布函数重构算子的应用

给出了格子Boltzmann方法(LBM)密度分布函数和温度分布函数的两个重构算子,解决了LBM与其它方法耦合的关键难题.构造了LBM与有限容积法(FVM)以及LBM与分子动力模拟方法(MD)的两个耦合模型.通过方腔自然对流以及平板通道内的泊肃叶流动对耦合模型进行验证.结果表明,基于重构算子的耦合模型可以正确地应用于L...     

过渡区微尺度流动的有效黏性多松弛系数格子Boltzmann模拟

为提高采用二维九速离散速度模型的格子Boltzmann方法 (LBM)模拟微尺度流动中非线性现象的精度和效率,引入Dongari等提出的有效平均分子自由程对黏性进行修正(Dongari N,Zhang Y H,Reese J M2011 J.Fluids Eng.133 071101);并针对以往研究微尺度流动时采用边界处理格式含有离散误差的问题,采用多松弛系数格子Boltzmann方法结合二阶滑移边界条件,对微尺度Couette流动和周期性Poiseuille流动进行模拟,并将速度分布以及质量流量等模拟结果与直接模拟蒙特卡罗方法模拟数据、线性Boltzmann方程的数值解以及现有的LBM模型模拟结果进行对比.结果表明,相对于现有的LBM模型,引入新的修正函数所建立的有效黏性多松弛系数LBM模型有效提高了LBM模拟过渡区的微尺度流动中的非线性现象的能力.     

超燃冲压发动机等截面隔离段中流动的LBM模拟

采用耦合的双分布函数格子Boltzmann方法模拟超燃冲压发动机中隔离段的流动.计算结果反映了隔离段内流场的基本结构,表明LBM在隔离段的模拟中是-种具有潜力的数值计算方法.     

不同精度格式的格子Boltzmann热模型的传热分析

通过在格子Boltzmann (LBM)热模型中添加参数项,使得在对应的宏观传热方程中,消除了一阶非线性误差项,具备二阶精度.通过Rayleigh-Benard对流数值试算,初步探索该二阶精度格式及其对应的一阶精度格式三个热模型的传热特征和适应性,并做出相应对比分析.针对一二阶精度模型在Ra数极高或热传导系数极大时,Nu数的计算与经验值相比出现较大偏差,分析LBM对应宏观热传导方程的截断误差后,在平衡分布函数中引进一个调节因子.通过调节对应宏观传热方程的截断误差项系数,校正Nu数的计算偏差,提高模拟精度,拓展模拟范围,增强了LBM作为一个数值方法在传热中的适应性.     

POD方法在扁管管翅式换热器研究中的应用

针对复杂流动与传热现象的扁管管翅式换热器板芯结构构建特征正交分解（POD）低阶模型,数值分析不同参数组合的传热特性并与有限容积法（FVM）计算结果进行比对.结果表明：POD计算结果与FVM计算结果间的相对偏差平均值在5.25%以内,但计算速度比FVM提高了164.78~2 021倍.所得结论对于扁管管翅式换热器结构参数设计和运行优化有理论参考价值.     

电热水器能耗数值计算及保温性能优化分析

电热水器中的能量传递过程主要包括电加热丝对内胆中水的加热过程,内胆中水的自然对流过程以及通过保温层的散热过程.本文采用多尺度数值模拟方法:用LBM确定保温材料导热系数作为基于连续介质的Fluent的输入数据,对发生在电热水器中的流动和耦合传热过程进行了数值模拟,并采用UDF来控制电加热丝的动态间歇加热工作过程,建立了模拟电热水器24 h固有能耗系数测试过程的数值模拟方法.采用该方法模拟得到的能耗系数与实验值符合良好,证明了该模型的准确性。同时模拟获得了电热水器中水在加热和散热阶段的流动特性、温度分布,以及各个区域的散热占比和热流密度。提出了优化电热水器保温隔热性能的判据,为电热水器保温结构的优化提供了依据.     

矩形管道载氚热磁流体输运特性的数值模拟研究

聚变堆液态金属包层矩形管道中的氚输运过程与磁流体动力学(MHD)流动传热过程耦合在一起，形成了复杂的载氚热磁流体输运特性。基于开发的MHD流动与传热数值模拟程序对矩形管道中液态金属MHD流动传热特性及其氚输运的影响进行了数值模拟。该程序首先求解了动量守恒方程，并与理论解进行了对比验证，然后与能量守恒方程耦合求解，得到了温度影响下矩形管道中的液态金属流场分布，在此基础上对强磁场高核热梯度影响下的氚浓度分布进行了数值模拟，得到了氚浓度在管道中的分布特性。结果显示，液态金属在矩形管道中的流动传热对氚输运过程产生了显著影响。     

格子Boltzmann方法应用于气动声学研究

应用格子Boltzmann方法(LBM)对不同类型的气动声学问题进行数值研究.通过模拟一维平面声波和二维点源声波的传播,得到沿传播方向的声压脉动,其振荡幅值和衰减趋势与理论值相吻合.其次进行声波衍射和干涉现象的数值模拟.最后,模拟处于流场中运动声源辐射声场的多普勒效应.模拟结果说明LBM方法能较好地模拟低马赫数下的声学问题,包括声压脉动的传播,声波的波动特性以及流动与声波间的相互作用.     

采用双分布函数格子Boltzmann方法的直流氩等离子体流动与传热分析

本文采用双分布函数格子Boltzmann方法(DDF-LB)模拟了直流氩电弧等离子体中的流动与传热问题。针对非局域热平衡状态下热等离子体中的原子、离子和电子三种组分,推导给出了不同组分的格子Boltzmann方程,通过耦合迭代求解得到各组分的温度和速度分布。文中采用随温度变化的松弛时间用以提高数值计算的稳定性。通过将本文的计算结果与文献结果进行比较,验证了DDF-LB方法求解等离子体中能量输运问题的稳定性和数值精度。     

L型混沌流道流动与传热特性的数值分析

采用周期性模型对L型混沌流道内流体流动与传热性能进行数值模拟,并与普通平直流道相对比,分析得到了两种流道内流体流场与温度场分布的细观信息,同时还讨论了不同Re数对混沌流道内流体传热与流动的影响。设计制作了测试流体流动性能的冷模试验装置,采用激光多普勒测速仪对L型混沌流道中流道内特殊点的速度进行测量,并将试验测量结果与数值模拟计算结果进行对比。结果标明:试验测量结果与数值模拟计算结果总体变化趋势基本一致,主流方向上速度相对误差最大为13.8%,在工程允许的误差范围之内,由此验证了数值模拟计算的可靠性。     

Lattice Boltzmann simulation of viscous-fluid flow and conjugate heat transfer in a rectangular cavity with a heated moving wall

In the present work, conjugate heat transfer in a rectangular cavity with a heated moving lid is investigated using the lattice Boltzmann method (LBM). The simulations are performed for incompressible flow, with Reynolds numbers ranging from 100 to 500, thermal diffusivity ratios ranging from 1 to 100, and Prandtl numbers ranging from 0.7 to 7. A uniform heat flux through the top of the lid is assumed. Results show that LBM is suitable for the study of heat transfer in conjugate problems. Effects of the Reynolds number, the Prandtl number and the thermal diffusivity ratio on hydrodynamic and thermal characteristics are investigated and discussed. The streamlines and temperature distribution in flow field, dimensionless temperature and Nusselt number along the hot wall are illustrated. The results indicate that increase of thermal diffusivity yields the removal of a higher quantity of energy from lid and its temperature decreases when increasing the Reynolds and the Prandtl numbers.     

格子玻尔兹曼方法对不同张角聚焦声束的建模

高强度聚焦超声(HIFU)是一种新型的无创治疗肿瘤新技术,其中换能器声场数值计算能够为HIFU治疗提供重要的依据。传统非线性KZK和SBE模型广泛应用于换能器声场数值计算,但依然存在某些不足。我们采用一种介观尺度的新型流体力学方法,即格子Boltzmann方法(LBM),基于2维9离散速度(D2Q9)格子构建了轴对称多弛豫参数LBM模型,并通过调节弛豫参数分析其对模型的影响;利用该模型对两个具有不同张角的球面聚焦换能器的声场进行数值模拟,并与KZK和SBE模型的计算结果进行比较。结果表明LBM模型能够很好地描述超声波的激发和传播机制,从流体力学的角度描述聚焦声场的分布,具有清晰的物理意义,且计算过程不受换能器张角的限制,在换能器声场的理论分析和模拟计算及其在HIFU治疗中的应用有着积极的意义。      

耦合不可压流场输运方程的格子Boltzmann方法研究

基于格子Boltzmann方法,提出了求解耦合不可压缩流场输运方程的一种改进数值方法. 该方法使用格子Boltzmann方法求解流场方程,并根据流场格子模型的密度分布函数构建了输运方程的二阶离散格式. 通过二维平板通道流场输运系统验证了该方法的有效性.数值结果表明,该方法可以有效地减少计算过程中出现的非物理耗散, 并克服了传统模型所需巨大存储量的缺点.     

An improved bounce-back scheme for complex boundary conditions in lattice Boltzmann method

In recent years, the lattice Boltzmann method (LBM) has been widely adopted to simulate various fluid systems, and the boundary treatment has been an active topic during the LBM development. In this paper, we present a novel approach to improve the bounce-back boundary treatment for moving surfaces with arbitrary configurations. We follow the framework originally proposed by Ladd [A.J.C. Ladd, Numerical simulations of particulate suspensions via a discretized Boltzman equation. Part 1. Theoretical foundation, Journal of Fluid Mechanics 271 (1994) 285–309]; however, the adjustment in the density distribution during the bouncing-back process at the boundary is calculated using the midpoint velocity inter-/extrapolated from the boundary and fluid velocities, instead of the real boundary velocity in the Ladd method. This modification ensures that the bouncing-back process and the density distribution adjustment both take place at a same location: the midpoint of a boundary lattice link, and thus removes the discrepancy of bouncing-back at the midpoint but density distribution adjustment at the boundary point in the original Ladd method. When compared with other existing boundary models, this method involves a simpler algorithm and exhibits a comparable or even better accuracy in describing flow field and flow-structure interaction, as demonstrated by several test simulations. Therefore, this boundary method could be considered as a competitive alternative for boundary treatment in LBM simulations, especially for particulate and porous flows with large fluid–solid interfacial areas.     

Heat transfer and fluid flow of benard-cell convection in rectangular container with free surface sensed by infrared thermography

The natural convection flow phenomena that occur inside an enclosed space are very interesting examples of complex fluid systems that may yield to analytical, empirical and numerical solutions, and many reports have looked into this basic problem. In the present study, heat transfer and fluid flow for natural convection in a horizontal rectangular container with a free surface are investigated using infrared thermography. The temperature field was measured and visualized at a gas-liquid (air — silicon oil) interface using infrared thermography. The heat transfer phenomena were also investigated by statistically analyzing the temperature data. The applicability of the infrared thermography to quantitative heat transfer measurement at the gas-liquid interface was evaluated. It is revealed that infrared thermography is effective not only in visualization of a gas-liquid interface but also in heat transfer measurement. A new heat transfer correlation is proposed for the gas-liquid interface of this flow system. The coefficient of heat transfer can be summarized by a specific heat transfer correlation formula regardless of several conditions, including container aspect ratio, fluid viscosity and fluid layer depth.     

Numerical modeling of the effect of number of electrodes on natural convection in an EHD fluid

Numerical modeling of the electrohydrodynamic effect on natural convection in enclosures is investigated. The interactions between electric field, flow field, and temperature field are analyzed by using the computational fluid dynamics technique. Flow pattern and temperature distribution are substantially affected by the voltage supplied at the wire electrodes, especially at low Rayleigh number. It can be concluded that the fluid velocity and heat transfer coefficient in the presence of an electric field are significantly increased with when a large number of electrodes is used. Surprisingly, a minimum value of augmented heat transfer occurs with an intermediate number of electrodes. The optimized condition between the aspect ratio of an enclosure and number of electrodes which leads to maximum heat transfer enhancement is expressed in terms of the parameters concerned.     

Heat transfer analysis for stretching flow over a curved surface with magnetic field

The analysis of a viscous fluid flow and heat transfer is carried out under the influence of a constant applied magnetic field over a curved stretching sheet. Heat transfer analysis is carried out for two heating processes, namely, prescribed surface temperature (PST) and prescribed heat flux (PHF). The equations governing the flow are modeled in a curvilinear coordinate system (r, s, z). The nonlinear partial differential equations are then transformed to nonlinear ordinary differential equations by using similarity transformations. The obtained system of equations is solved numerically by a shooting method using Runge-Kutta algorithm. The interest lies in determining the influence of dimensionless radius of curvature on the velocity, temperature, skin friction, and rate of heat transfer at the wall prescribed by the Nusselt number. The effects of Hartmann number are also presented for the fluid properties of interest.     

Experimental Investigations on Latent Heat Storage Unit using Paraffin Wax as Phase Change Material

Thermal performance of a latent heat storage unit is evaluated experimentally. The latent heat thermal energy storage system analyzed in this work is a shell-and-tube type of heat exchanger using paraffin wax (melting point between 58°C and 60°C) as the phase change material. The temperature distribution in the phase change material is measured with time. The influence of mass flow rate and inlet temperature of the heat transfer fluid on heat fraction is examined for both the melting and solidification processes. The mass flow rate of heat transfer fluid (water) is varied in the range of 0.0167 kg/s to 0.0833 kg/s (1 kg/min to 5 kg/min), and the fluid inlet temperature is varied between 75°C and 85°C. The experimental results indicate that the total melting time of the phase change material increases as the mass flow rate and inlet temperature of heat transfer fluid decrease. The fluid inlet temperature influences the heat fraction considerably as compared to the mass flow rate of heat transfer fluid during the melting process of the phase change material.     

Lattice Boltzmann Simulation for Complex Flow in a Solar Wall

In this letter,we present a lattice Boltzmann simulation for complex flow in a solar wall system which includes porous media flow and heat transfer,specifically for solar energy utilization through an unglazed transpired solar air collector(UTC).Besides the lattice Boltzmann equation(LBE) for time evolution of particle distribution function for fluid field,we introduce an analogy,LBE for time evolution of distribution function for temperature.Both temperature fields of fluid(air) and solid(porous media) are modeled.We study the effects of fan velocity,solar radiation intensity,porosity,etc.on the thermal performance of the UTC.In general,our simulation results are in good agreement with what in literature.With the current system setting,both fan velocity and solar radiation intensity have significant effect on the thermal performance of the UTC.However,it is shown that the porosity has negligible effect on the heat collector indicating the current system setting might not be realistic.Further examinations of thermal performance in different UTC systems are ongoing.The results are expected to present in near future.