方腔内耦合传热传质非对称流动数值模拟

对带有质热源的方腔内流体传热传质进行数值研究。针对不同Ra、Nc、Sr和Df,探究对称方腔内流体传热传质的分岔特性。结果表明：存在临界Ra_c使流体流动形态发生转变,当Ra<Ra_c时,流体流线、温度场和浓度场对称分布;当Ra>Ra_c时,流体发生偏斜。增大浮升力,流体更易发生分岔现象。增强Soret和Dufour效应可增强传热对称性并增大流体发生分岔的临界Rayleigh数。     

具有质热源的方腔内对流传热传质及其非线性特性

以火灾为背景,对一个带有质热源的方腔进行了研究。针对不同Ra,浮升力比Nc,S_r和D_f,探讨了腔内流体的传热传质规律及其非线性特性。结果表明,存在临界R{a}_{\mathrm{c}}使流体流动由导热驱动流动转变为对流驱动流动,随着Ra增大质热源表面传热传质强度增加,流体由稳态流动转为振荡流动。当Nc>-1时,对流传热系数和对流传质系数增大,流体由稳态转化为倍周期振荡,最后转化为混沌。增大S_r和D_f可增强传热传质能力。     

Natural convection heat transfer enhancement of different nanofluids by adding dimple fins on a vertical channel wall

In this numerical study, natural flow and heat transfer of nanofluids with Al_2O_3, TiO_2, Cu and CNT nanoparticles in a vertical channel with dimpled fins at Rayleigh number(Ra) of Ra = 3.25 × 10~7 to Ra = 1 × 10~8 are investigated by using the finite volume method. The obtained results revealed that, using CNT in volume fractions of 2% and 4%leads to significant heat transfer and at φ = 6%, using TiO_2 nanoparticles has a great effect on Nu number enhancement. Also, using solid nanoparticles in base fluid causes more uniform heat transfer distribution,especially in areas close to heated surface and by adding more volume fraction in base fluid, temperature level reduces. In general, according to temperature contours, reduction of wall temperature depends on the increase of Ra and volume fraction and the type of solid nanoparticles.     

ROUTE TO CHAOS IN CHEMICALLY ENHANCED THERMAL CONVECTION IN POROUS MEDIA

Exothermic chemical reactions can influence natural convection effects in a porous medium. Such phenomena may occur in tubular reactors, oxidation of solid materials in large containers, chemical vapor deposition systems, liquid explosives, and others. Experimental evidence indicates that the influence of natural convection in many chemically reacting systems cannot be neglected.

In the present work, transient effects of a two-dimensional convection generated and sustained by an exothermic chemical reaction and a constant boundary wall temperature are studied. The Darcy-Boussinesq equations are used to describe fluid flow through porous media with a zero-order chemical reaction. A recently developed method is used to compute the singular points (such as limit and symmetry breaking bifurcation points) precisely as a function of Rayleigh (Ra) and Frank-Kamenetskii (Fk) numbers. Fold curves are drawn as a function of Ra, and Fk. Flow behavior is governed by two natural parameters, namely, Ra and Fk. Multiple stationary solutions arise over an intermediate range of Ra. The solution structure, i.e., the interconnections between the various branches appear quite complicated. Determination of linear stability on the these solution branches reveal that all two-dimensional, stationary solutions develop some form of instability at higher values of Ra. Transient simulations reveal the emergence of time periodic solutions at higher Ra. The nature and frequency of these periodic solutions depend on the route followed in the parameter space Fk and Ra. The various routes to chaos are identified in this parameter space.     

耦合表面张力的封闭腔体内管外自然对流传热特性

采用格子-玻尔兹曼方法（LBM）构建了考虑表面张力影响下封闭腔体内加热圆柱外自然对流数学模型。在分别对自然对流与表面张力模型模拟验证的基础上,探究了液固表面张力与重力场下浮升力同时作用下加热圆柱外流体自然对流的传热特性。结果表明,在给定Ra数下（10³≤Ra≤10⁶）,随着表面张力的增大（Oh数减小）,封闭腔体内管外自然对流扰动会加剧,流型会变复杂,壁面换热效率有明显提高。在Ra数为10⁵,长度比（加热圆柱半径和方腔长度之比）为1/10情况下,加入表面张力σ=0.076302N/mm（Oh=0.122）,腔体左侧壁面平均努塞尔数和加热圆柱壁面与未加表面张力相比分别提高了93.5%和60.35%。当表面张力大小和自然对流的浮升力相当时,此时的流场波动更加明显剧烈,在一定范围内增大浮升力反而减弱换热。     

COMPREHENSIVE, THEORETICALLY BASED, CORRELATING EQUATIONS FOR FREE CONVECTION FROM ISOTHERMAL SPHERES

Correlating equations are developed for the local and mean Nusselt number for free convection from an isothermal sphere as a function of the Rayleigh and Prandtl numbers. These expressions are based primarily on theoretical solutions for limiting cases, and hence are presumed to be more reliable than purely empirical correlations. The predictions of the proposed expressions are, however, validated by comparisons with prior experimental data. The expressions for the mean Nusselt number are shown to be applicable for all Ra and Pr. The expressions for the local Nusselt number are limited in applicability to the laminar boundary layer regime. The same equations are applicable to mass transfer and to combined heat and mass transfer in terms of the Sherwood, Schmidt and appropriately modified Rayleigh numbers.     

方形微通道内超临界CO2流动换热特性研究

采用SST k-ω湍流模型对加热条件下超临界CO₂在方形微通道内的流动换热特性进行了数值模拟。通过对比三种壁面平均传热系数、浮升力参数和二次流强度的沿程变化研究了管型、热通量、质量流量和倾斜角度对微通道内流动换热性能的影响。结果表明：水平方形微通道的整体换热性能优于相同水力直径的半圆形微通道。流体域典型截面的温度分布、速度分布和湍动能分布等信息可以很好地解释水平方向流动时上、下壁面传热差异的现象。减小热通量、增大质量流量或减小流体流动方向与重力方向之间的夹角,可提高方形微通道的整体换热水平。该模拟结果对以超临界CO₂为工质的微通道换热器的设计和优化具有一定的理论指导意义。     

A numerical study of the Dufour and Soret effects on unsteady natural convection flow past an isothermal vertical cylinder

The Dufour and Soret effects on the unsteady laminar free convective flow with mass transfer flow past a semi-infinite isothermal vertical cylinder were studied numerically. The governing partial differential equations were converted into a non-dimensional form and solved numerically by applying a Crank-Nicolson type of implicit finite-difference method with a tri-diagonal matrix manipulation and an iterative procedure. For the hydrogen-air mixture, which is a non-chemical reacting fluid, the profiles of the unsteady dimensionless velocity, temperature and concentration are shown graphically for the different values of thermal and mass Grashof numbers, thermal diffusion parameters (Soret numbers) and diffusion-thermo parameters (Dufour numbers). Finally, the simulated values of the average skin-friction coefficient, the average Nusselt number and the average Sherwood number are presented. The numerical results reveal that for an increasing Soret number or decreasing Dufour number, the time to reach the temporal maximum and the steady-state decreases for the flow variables. As the Soret number increases or the Dufour number decreases, both the skin friction and the Sherwood number increase, whereas the Nusselt number decreases.     

高热流低流速条件下超临界CO2在小圆管内的对流传热特性

为获取高热流、低流速条件下超临界CO₂的传热规律，开展了超临界CO₂在内径2 mm水平小圆管内对流传热试验研究，并重点探讨了变物性、浮升力和热加速等效应对传热过程的影响。试验参数范围：系统压力7.6~8.4 MPa，质量流速400~500 kg/(m²?s)，热通量0~200 kW/m²，流体温度20~60℃，Reynolds数1.2×10⁴~4.3×10⁴。分别采用Gr/Re ²和Kv作为浮升力效应和热加速效应的判别因子。结果显示，在高热流低流速工况下，浮升力效应显著（Gr/Re ² > 10^-3），同一个截面处的上壁面传热系数始终小于下壁面传热系数。浮升力效应是高热流低流速工况下传热恶化的主要诱发因素。试验中热加速因子较小（Kv < 8.5×10^-7），其效应可以忽略。将试验数据与典型的传热经验关联式作对比，结果表明Liao-Zhao关联式的计算结果与试验结果最吻合。     

A NEW CLASS OF SIMILARITY SOLUTIONS OF AN UNSTEADY ELECTRICALLY CONDUCTING FREE-FORCED CONVECTIVE FLOW IN A VERTICAL POROUS SURFACE WITH DUFOUR AND SORET EFFECTS

The 2-D unsteady magnetohydrodynamic free-forced convective boundary layer flow of a viscous incompressible fluid is studied numerically taking into account heat and mass transfer. The fluid is subjected to uniform heat and mass fluxes embedded in a porous medium by the presence of coupled Dufour and Soret effects. A new class of similarity equations has been obtained by introducing a time-dependent length scale and a corresponding similarity variable. The resulting equations are then integrated numerically using the Nachtsheim-Swigert shooting iteration technique along with the sixth-order Runge-Kutta integration scheme. By developing locally similar solutions of the fluid flow, the behavior of the velocity, temperature, and concentration fields as well as the rate of heat transfer, wall temperature gradient, rate of mass transfer, and skin friction coefficient have been investigated. The effects of Grashof number (Gr), modified Grashof number (Gm), combined effects of the porous and magnetic parameter (S), suction/injection parameter Fw, Brinkman number (Br), Soret number (Sr), and Dufour number (D_f) have been observed on the flow field and discussed.     

Thermal performance and entropy generation for nanofluid jet injection on a ribbed microchannel with oscillating heat flux: Investigation of the first and second laws of thermodynamics

In current numerical study, forced flow and heat transfer of water/NDG (Nitrogen-doped graphene) nanofluid in nanoparticles mass fractions (φ) of 0, 2% and 4% at Reynolds numbers (Re) of 10, 50, 100 and 150 are simulated in steady states. Studied geometry is a two-dimensional microchannel under the influence of nanofluid jet injection. Temperature of inlet fluid equals with T_c=293 K and hot source of microchannel is under the influence of oscillating heat flux. Also, in this research, the effect of the variations of attack angle of triangular rib (15°, 30°, 45° and 60°) on laminar nanofluid flow behavior inside the studied rectangular geometry with the ratio of L/H=28 and nanofluid jet injection is investigated. Obtained results indicate that the increase of Reynolds number, nanoparticles mass fraction and attack angle of rib leads to the increase of pressure drop. By increasing fluid viscosity, momentum depreciation of fluid in collusion with microchannel surfaces enhances. Also, the increase of attack angle of rib at higher Reynolds numbers has a great effect on this coefficient. At low Reynolds numbers, due to slow motion of fluid, variations of attack angle of rib, especially in angles of 30°, 45° and 60° are almost similar. By increasing fluid velocity, the effect of the variations of attack angle on pressure drop becomes significant and pressure drop figures act differently. In general, by using heat transfer enhancement methods in studied geometry, heat transfer increases almost 25%.     

Numerical study on non-uniform heat transfer deterioration of supercritical RP-3 aviation kerosene in a horizontal tube

The convective heat transfer of supercritical-pressure RP-3(Rocket Propellant 3) aviation kerosene in a horizontal circular tube has been numerically studied, focusing mainly on the non-uniform heat transfer deterioration along the circumferential direction. The governing equations of mass, momentum and energy have been solved using the pressure-based segregated solver based on the finite volume method. The re-normalization group(RNG) k-ε turbulence model with an enhanced wall treatment was selected. Considering the heat conduction in the solid wall, the mechanism of heat transfer deterioration and the buoyancy effect on deteriorated heat transfer were discussed. The evolution of secondary flow was analyzed. Effects of the outer-wall heat flux,mass flux, pressure and tube thermal conductivity on heat transfer were investigated. Moreover, the buoyancy criterion and the heat transfer correlation were obtained. Results indicate that the poor flow performance of near-wall fluid causes the pseudo-film boiling, further leads to the heat transfer deterioration. The strong buoyancy has an effect of enhancing the heat transfer at the bottom of tube, and weakening the heat transfer at the top of tube, which results in the non-uniform inner-wall temperature and heat flux distributions. Decreasing the ratio of outer-wall heat flux and mass flux, increasing the pressure could weaken the heat transfer difference along the circumferential direction, while the effect of thermal conductivity of tube on the circumferential parameters distributions is more complicated. When the buoyancy criterion of(Gr_q/Gr_(th))_(max)≤ 0.8 is satisfied, the effect of buoyancy could be ignored. The new correlations work well for non-uniform heat transfer predictions.     

封闭腔内水自然对流换热数值模拟

为了揭示封闭腔内非Boussinesq流体在浮力驱动下所特有的流动换热现象和形成机理，采用CFD软件Fluent对封闭腔内水的自然对流进行数值模拟，得到矩形封闭腔高宽比、Rayleigh数、倾斜角度、壁面温度差对流动和传热的影响规律。研究结果表明：由于水的密度在3.98℃达到最大，两竖壁面温度跨越这一点时会引起流动图像反转；具有流动反转的双涡结构降低了对流换热平均Nusselt数；相同Rayleigh数下,高宽比为1对应对流换热平均Nusselt数最大值；倾斜角度对平均Nusselt数影响与Rayleigh数和温度边界条件有关。     

微通道内超临界氮气三维热流场实验与数值模拟

利用实验与数值模拟相结合的方法研究了超临界氮气（SCN₂）三维热流场特性，在用实验数据验证数值方法可靠性基础上，分析了压力（3.6~7 MPa）和质量流速［800~1200 kg/（m²?s）］对SCN₂对流传热特性的影响，揭示了微通道圆管不同圆周方向上SCN₂热流场规律。结果表明：在低压力和高质量流速下，同一轴向位置处，径向内壁温最大值出现在圆管90°处；质量流速越大，内壁温最大值和对流传热系数最小值由圆管180°向90°处发生了转移；当浮升力系数Gr*/Re²>1时，浮升力有利于强化圆管底部区域流体传热能力；基于获得的数据，提出了一个新的适合预测微通道圆管内SCN₂对流传热特性的无量纲换热关联式，预测误差小于20%。研究结果为微通道换热器优化设计提供了参考。     

绝热圆腔内4根冷热微管阵列振动强化传热实验

目前针对微细管管外同时有自然对流和周期横向振动存在情况下的混合对流问题鲜有研究。为了探究这种现象, 进而开发一种新型的换热器, 本文提出了一种由4根冷热微细管组成的微管阵列, 将其平行放置于有中间流体的圆形腔体内, 并对管束施加微小横向振动的混合对流传热机理进行研究。详细描述了实验系统及实验方法, 以及后期对实验数据的处理方法, 得到了传热Nu数随管间自然对流强度Ra数, 强迫振动Re_v数等量纲为1参数的变化曲线;揭示了该结构形式下传热特性随冷热管温差、振动频率、微管振幅等传热规律。结果表明, 横向振动在小温差情况下传热强化效果更为明显。将本文实验数据与Lemlich单管振动实验中的数据进行对比, 表明了振动对4根微管阵列的强化传热效果比单根管更明显。     

冷热圆管在封闭方腔内不同垂直位置的自然对流数值研究

利用有限体积法对冷热圆管在封闭方腔内不同垂直位置的自然对流现象进行了数值研究。讨论了瑞利数Ra和冷热圆管间距δ对方腔内自然对流流动与换热的影响, 其中瑞利数的变化范围为10³～10⁶, 圆管间距变化范围为0.3～0.6。为了揭示冷热圆管间的相互作用和圆管与方腔间的相互作用对自然对流换热与流动的影响规律, 比较分析了热圆管在上、冷圆管在下和热圆管在下、冷圆管在上两种情形下冷热圆管、方腔的自然对流换热能力的差异。研究表明:瑞利数的改变, 对方腔内温度场分布和涡流结构有显著影响;热圆管在下、冷圆管在上这种情形更有利于自然对流换热的进行;增加圆管间距δ, 热圆管和方腔的换热能力增强, 但冷圆管的换热能力却有所减弱。研究结果为核电站安全壳非能动余热排出系统的性能研究提供了理论依据。     

Power-law shear-thinning flow around a heated square bluff body under aiding buoyancy at low Reynolds numbers

Power-law shear-thinning fluid flow over a heated square bluff body is numerically investigated under aiding buoyancy mixed convection at low Reynolds numbers. Semi-explicit finite volume code is developed to solve the governing equations along with the appropriate boundary conditions. Both aiding buoyancy and shear-thinning natures are found to augment the heat transfer rate from the surface of the long square bar. In aiding buoyancy, the total drag coefficient is found to be more for the square cylinder than that of the circular cylinder, whereas the average cylinder Nusselt number for the square cylinder is found to be lower than the circular one on equal side/diameter basis. Maximum augmentation in heat transfer is found to be approximately 20% with respect to forced convection. Finally, a heat transfer correlation is established by using the Colburn heat transfer factor.     

圆管束中导流器对其自然对流换热的影响

对无限大空间中5根同规格圆管组成的圆管束在管间放置导流器的自然对流换热进行了数值模拟研究。考虑了Ra在10³～10⁴范围内,导流器偏转角为0°～60°,圆管间距为2～4倍圆管直径的自然对流换热,分析了5根圆管的局部Nusselt数（Nu_loc）和平均Nusselt数（Nu_ave）。研究结果表明,导流器对管束结构的自然对流换热影响体现在两方面：一是对导流器下方圆管而言相当于障碍物削弱其换热,二是对导流器上方圆管而言隔绝了下方羽状流的影响,从而增强系统的整体换热。在圆管间距较大时,与无导流器的圆管束相比,C1～C4各圆管换热下降趋势明显放缓,并且从C4、C5圆管换热趋势上升,系统整体换热增强。圆管间距S=2D,Ra=10³,导流器主要起障碍物的作用,削弱C1～C4各圆管换热,导致系统整体换热下降。当导流器偏转角度为45°时系统换热达到最大值,较无导流器时换热最多有着21%的提升。     

Heat and mass transfer enhancement in a double diffusive mixed convection lid cavity under pulsating flow

This work assess the effectiveness of the pulsating flow regime for the enhancement of heat and mass transfer in double diffusive problems by studying a two-dimensional, heated, lid cavity. The research characterizes the influence of pulsating parameters such as temporal frequency, wave number and amplitude in the process. Results show that the pulsating regime enhances heat/mass transfer within a square cavity up to a 14%/38% respectively with respect to the non-pulsating case, due to the promotion of additional shear stress fields. As Richardson number, Brownian diffusion or the solute/thermal buoyancy ratio increase or the cavity becomes narrower, heat/mass transfer increases.     

The Soret and Dufour Effects in Non-thermal Equilibrium Packed Beds with Forced Convection and Endothermic Reactions

To study the influence of the Soret and Dufour effects on the reactive characteristics of a porous packed bed with endothermic reactions and forced convection, a two-dimensional mathematical model cons...