孔板空调风口送风射流的数值模拟

介绍Ｎ点风口模型用于数值模拟室内空气流动时描述孔板类送风口的入流边界条件．然后采用该风口模型对不同的孔板风口出流条件算例进行数值计算,并就轴心速度衰减、射流扩展角以及断面流速分布等射流特性与实验数据进行了对比．比较结果表明,Ｎ点风口模型用于描述数值模拟室内空气流动的孔板类风口入流边界条件,可以获得工程上足够满意的结果．     

燃气射流冲击传热特性的数值模拟

针对射流传热问题,利用基于RNGk-ε湍流模型的数值方法模拟了射流垂直冲击平板的流动过程,并与实验数据比较,验证了模型的可行性。在此基础上,以火箭喷管入口参数为入口条件,建立了超音速燃气射流垂直冲击平板和冲击浸没平板的计算模型,分析了不同冲击条件下努塞尔数分布规律和温度分布规律, 论述了超音速射流传热的特性及影响传热特性的因素。得到了冲击距离为(14~18)D的努塞尔数取值范围,并表明冲击距离和射流温度是影响传热效率的关键因素;冲击距离增加,传热效率降低,冲击平板表面的射流温度越高,传热效率越高。     

方型散流器空调室内空气流动的数值模拟

用N点风口动量模型和一个新零方程湍流模型对某办公室方型散流器空调的室内温度场和速度场进行了模拟,并和实验数据进行对比.结果表明,计算所得速度和温度分布与实测值吻合得很好,所用的风口模型和湍流模型能快速地将方型散流器空调通风的温度和速度场合理地模拟出来,可用于指导和优化同类空调通风气流组织设计.     

倾斜管道内圆柱状颗粒运动的数值模拟

本文采用计算流体力学–离散单元法对倾斜管道内圆柱颗粒在气体中的运动进行数值模拟,分别分析了颗粒密度、颗粒长径比、吸入速度、管道倾角、管道长径比等因素对颗粒输运效率的影响。结果表明,颗粒密度越大,颗粒的输运效率越低;输运效率随着颗粒长径比的增大而减小;输运效率随着进口速度的增大而增大。管道的倾角对输运效率存在影响,但二者不是单调关系;颗粒输运效率与管道长径比成反比。

     

室内空气流动数值模拟的N点风口动量模型

为适应工程应用中快速、准确模拟室内空气流动的需要,提出N点风口动量模型,以简化描述利用计算流体动力学CFD方法模拟室内空气流动时百叶、多孔板类送风口的入流边界条件。百叶和多孔板风口的等温自由射流算例以及HESCO孔板类散流器在室内送风的算例和实验数据对比表明,N点风口动量模型可以较好地解决数值模拟室内空气流动的风口入流边界条件描述问题。     

带自由表面水流与离散颗粒耦合模型研究

基于考察泥沙运动的细观行为特征,采用离散单元法(DEM)模拟泥沙颗粒运动,结合带自由表面的水动力学计算模型,建立了CFD-DEM耦合数值模型。计算程序开发基于Fortran语言来实现。耦合模型中实现了硬球模型和软球模型两种颗粒碰撞模型,应用范围较广。作为自由表面水流与泥沙颗粒流数值模型的初步研究,在模型建立的基础上,对模型做了基本的验证。分别通过单颗粒静水沉降和混合颗粒群分选两个计算工况,验证了模型的正确性及模拟精度。该耦合模型可进一步丰富带自由表面水流条件下泥沙运动的研究手段。     

基于离散元方法的颗粒系统休止角数值模拟研究

颗粒系统堆积的休止角是颗粒学中一个基础性的科学问题,深入认识其影响因素和变化规律对于土木、化工等工程优化具有重要意义.目前基于实验测量的研究受限于可用的颗粒种类和测量手段,难以完备地揭示各种物理参数对休止角的影响规律.离散元(Discrete element method,DEM)是一种可直接对复杂的大型颗粒系统的运动进行仿真的数值方法.本文借助DEM方法对颗粒系统的休止角问题开展了高精度数值研究,揭示了影响颗粒堆积休止角的重要颗粒物性参数.在DEM建模方面,为了表征非球形度对仿真的影响,本研究采用了滚动摩擦系数,并利用既有的实验数据进行了模型验证.同时,本文利用堆积体在垂直平面内投影的二值化图像来合理计算休止角.仿真结果表明,滑动和滚动摩擦系数均与休止角的大小呈正相关,滑动摩擦系数的增加可将休止角提升一倍,而滚动摩擦系数对休止角的影响存在一个上限阈值,达到该阈值后将不能继续增加休止角的大小;颗粒的杨氏模量、回弹系数对休止角的影响不显著.此外,本文利用DEM方法研究了堆积体的传热特性,发现在加热过程中,堆积体内部具有最高热流率的区域从堆积体的底部迁移到上部.     

流化床生物质燃烧过程中传热特性的CFD–DEM数值模拟

本文基于CFD–DEM方法,将多种传热模型与稠密气固流动与化学反应耦合,对鼓泡流化床内生物质燃烧过程中的传热特性进行了数值模拟研究。结果表明,对于燃料颗粒,反应热占据主导地位,辐射和对流传热同样起到重要作用,导热的影响相对较小。随着燃料挥发分含量的增加,对流传热和颗粒–颗粒导热作用增强,颗粒–壁面导热减弱。床内热点温度由燃料挥发分燃烧引起,从床层表面产生并上升到自由空域。生物质颗粒的温度则随局部颗粒浓度的减小而逐渐增加。

     

碎石料直剪实验的组合颗粒单元数值模拟

通过构造三维组合颗粒单元来描述颗粒间的互锁效应,对非规则颗粒材料的力学行为进行了离散元数值模拟,并通过碎石料的直剪实验进行了验证.该组合颗粒的质量与碎石块具有相同的概率分布特性,其几何形态则由不同数目、镶嵌尺寸、组合方位和粒径的球形颗粒进行随机构造.组合颗粒单元在局部与整体坐标之间的转动、力矩和方位关系通过四元素方法进行确定;颗粒之间的作用力采用具有Mohr-Coulomb摩擦定侓的Hertz-Mindlin 非线性接触模型,并考虑了非线性法向粘滞力的影响.在不同的法向应力下,对碎石料在直剪实验中的剪切应力和剪胀现象进行了离散元模拟,计算结果与实测结果相吻合;此外,在不同的法向应力和接触摩擦系数下,对碎石料的有效摩擦系数进行了计算和讨论.本文工作验证了组合颗粒单元在非规则颗粒材料的离散元模拟中的可行性.     

基于离散单元法和人工神经网络的近壁颗粒动力学特征研究

颗粒与壁面的相互作用往往对颗粒流动具有显著影响. 为研究颗粒与壁面作用机理, 对滚筒内颗粒流动过程进行离散单元法(DEM)数值模拟. 基于模拟结果统计分析靠近壁面处颗粒的运动特征, 结果表明, 小摩擦系数时颗粒平动和旋转速度均近似满足正态分布, 但由于壁面影响, 摩擦系数增大时颗粒沿滚筒轴向的旋转速度偏离正态分布, 颗粒动力学理论推导壁面边界条件时应考虑速度正态分布的修正及速度脉动的各向异性. 采用人工神经网络(ANN)构建了颗粒无因次旋转温度、滑移速度和平动温度之间的函数模型, 进而可以在常规双流模型壁面边界条件中考虑颗粒旋转的影响. 基于DEM模拟及结果分析可以为壁面边界条件的理论构造和半经验修正提供基础数据和封闭模型.      

Streamline upwind finite element method for conjugate heat transfer problems

This paper presents a combined finite element method for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow. The streamline upwind finite element method is used for the analysis of thermal viscous flow in the fluid region, whereas the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the proposed method is to consistently couple heat transfer along the fluid-solid interface. Three test cases, i.e. conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the efficiency of the present method. The English text was polished byYunming Chen.     

Experiments on the local heat transfer characteristics of a circulating fluidized bed

The role of particle diameter in the heat transfer of a gas–solid suspension to the walls of a circulating fluidized bed was studied for particles of uniform size. This work reports and analyzes new experimental results for the local bed to wall heat transfer coefficient, not including the radiation component, in a long active heat transfer surface length laboratory bed, which extend previous findings and clear up some divergences. The research included determining the effects of extension and location of the heat transfer surface, circulating solids mass flux and average suspension density. An experimental set-up was built, with a 72.5 mm internal diameter riser, 6.0 m high, composed of six double pipe heat exchangers, 0.93 m high, located one above the other. Five narrow sized diameter quartz sand particles − 179, 230, 385, 460 and 545 μm − were tested. Temperature was kept approximately constant at 423 K and the superficial gas velocity at 10.5 m/s. The major influence of suspension density on the wall heat transfer was confirmed, and contrary to other authors, a significant effect of particle size was found, which becomes more relevant for smaller particles and increasing suspension density. It was observed that the extension of the heat transfer surface area did not influence the heat transfer coefficient for lengths greater than 0.93 m.The heat transfer surface location did not show any effect, except for the exchanger at the botton of the riser. A simple correlation was proposed to calculate the heat transfer coefficient as a function of particle diameter and suspension density.     

Experimental study of heat transfer characteristics due to confined impinging two-dimensional jets

Experimental results are presented for characteristics of impingement heat transfer caused by three slot jets. Experimental values were obtained for the dimensionless distance H = 0.5−3, dimensionless pitch P = 6−16, and Reynolds number Re = 500−8000. For laminar impinging flow, they were compared with numerical results. For turbulent impinging flow, two peaks of the local Nusselt number were obtained behind the second nozzle. The position of the second peak approached the nozzle as the space between nozzle and impinged surface decreased. The average Nusselt number between the central and second nozzles was determined from the ratio P/H and the Reynolds number based on the pitch of the nozzles.     

A tube-by-tube reduction method for simultaneous heat and mass transfer characteristics for plain fin-and-tube heat exchangers in dehumidifying conditions

This study proposed a new method, namely a tube-by-tube reduction method to analyze the performance of fin-and-tube heat exchangers having plain fin configuration under dehumidifying conditions. The mass transfer coefficients which seldom reported in the open literature, are also presented. For fully wet conditions, it is found that the reduced results for both sensible heat transfer performance and the mass transfer performance by the present method are insensitive to change of inlet humidity. Unlike those tested in fully dry condition, the sensible heat transfer performance under dehumidification is comparatively independent of fin pitch. The ratio of the heat transfer characteristic to mass transfer characteristic (h_c,o/h_d,o C_p,a) is in the range of 0.6~1.0, and the ratio is insensitive to change of fin spacing at low Reynolds number. However, a slight drop of the ratio of (h_c,o/h_d,o C_p,a) is seen with the decrease of fin spacing when the Reynolds number is sufficient high. This is associated with the more pronounced influence due to condensate removal by the vapor shear. Correlations are proposed to describe the heat and mass performance for the present plate fin configurations. These correlations can describe 89% of the Chilton Colburn j-factor of the heat transfer (j_h) within 15% and can correlate 81% of the Chilton Colburn j-factor of the mass transfer (j_m) within 20%.     

A note on flow and heat transfer of a viscoelastic fluid over a stretching sheet

This paper presents a study of the flow and heat transfer of an incompressible homogeneous second grade fluid past a stretching sheet. The governing partial differential equations are converted into ordinary differential equations by a similarity transformation. The effects of viscous dissipation and work due to deformation are considered in the energy equation and the variations of dimensionless surface temperature and dimensionless surface temperature gradient with various parameters are graphed and tabulated. Two cases are studied, namely, (i) the sheet with constant surface temperature (CST case) and (ii) the sheet with prescribed surface temperature (PST case).     

Finite element analysis of convective heat transfer in porous media

The finite element method is used to analyse convective heat transfer in a porous medium. Convection past a vertical surface embedded in the medium and convection in a confined porous medium enclosure are analysed using the above method. The results are compared with those available in the literature and the agreement is found to be good. The method is applicable for two-dimensional analysis in a porous body of any arbitrary shape. The restriction of the boundary layer assumption is relaxed.     

Inertia dominated flow and heat transfer in liquid drop spreading on a hot substrate

The present work deals with computational modeling of the fluid flow and heat transfer taking place in the process of impact of a cold liquid drop (T_d = 20-25 °C) onto a dry heated substrate characterized by different thermophysical properties. The computational model, based on the volume-of-fluid method for the free-surface capturing, is validated by simulating the configurations accounting for the conjugate heat transfer. The simulations were performed in a range of impact Reynolds numbers (Re = 2000-4500), Weber numbers (We = 27-110) and substrate temperatures (T_s = 100-120 °C). The considered temperature range of the drop-surface, i.e. liquid-solid system does not account for the phase change, that is boiling and evaporation. The model performances are assessed by contrasting the results to the reference database originating from the experimental and complementary numerical investigations by Pasandideh-Fard et al. [Pasandideh-Fard, M., Aziz, S., Chandra, S., Mostaghimi, J., 2001. Cooling effectiveness of a water drop impinging on a hot surface. International Journal of Heat and Fluid Flow, 22, 201-210] and Healy et al. [Healy, W., Hartley, J., Abdel-Khalik, S., 2001. On the validity of the adiabatic spreading assumption in droplet impact cooling. International Journal of Heat and Mass Transfer, 44, 3869-3881]. In addition, the thermal field obtained is analyzed along with the corresponding asymptotic analytical solution proposed by Roisman [Roisman, I.V., 2010. Fast forced liquid film spreading on a substrate: flow, heat transfer and phase transition. Journal of Fluid Mechanics, 656, 189-204]. Contrary to some previous numerical studies, the present computational model accounts for the air flow surrounding the liquid drop. This model feature enables a small air bubble to be resolved in the region of the impact point. The reported results agree reasonably well with experimental and theoretical findings with respect to the drop spreading pattern and associated heat flux and temperature distribution.     

Numerical study of steady/unsteady flow and heat transfer in porous media using a characteristics-based matrix-free implicit FV method on unstructured grids

In this study, a high-resolution characteristic-based finite-volume (FV) method on unstructured grids [Int. J. Numer. Method Eng. 50 (2001) 11; Int. J. Heat Fluid Flow 21 (2000) 432] is extended by a matrix-free implicit dual-time stepping scheme for the numerical simulation of steady and unsteady flow and heat transfer with porous media. The method has been used to study the characteristics of a complex problem: flow and heat transfer in a channel with multiple discrete porous blocks, which was originally proposed by Huang and Vafai [J. Thermophys. Heat Transfer 8 (3) (1994) 563]. In addition, flow and heat transfer in a channel partially or fully filled with porous layers and containing solid protruding blocks with constant heat flux on its lower surface are also investigated in details. Hydrodynamic and heat transfer results are reported for both steady and transient flow cases. In particular, the effects of Darcy and Reynolds numbers on heat transfer augmentation and pressure loss are studied. An in-depth discussion of the formation and variation of recirculation is presented and the existence of optimum porous insert is demonstrated. At high Reynolds numbers the flow in the porous channel exhibits a cyclic characteristics although unlike the non-porous channel flow, the cyclic vortex development is only restricted to a small area behind the last solid block, while temperature changes more slowly and does not exhibit cyclic variations over a long period of time. It is shown that for all the cases studied altering some parametric values can have significant and interesting effects on both flow pattern as well as heat transfer characteristics.