含不同形状内热源的圆管内纳米流体自然对流及换热数值研究

采用格子玻尔兹曼方法对有三种恒温热源（圆形、三角形、方形）参与的圆管内纳米流体（铜-水）自然对流进行数值研究。主要研究瑞利（Ra）数,纳米颗粒体积分数以及热源几何形状等控制参数对纳米流体的流动与传热的影响。结果发现纳米颗粒体积分数的增加有利于强化传热,且在Ra数较小时,平均努塞尔（Nu）数增加的幅度要优于Ra数较大的情况。在所研究的控制参数范围内,方形热源的平均Nu数最大。根据数值结果给出不同热源表面的平均Nu数、纳米颗粒体积分数、Ra数三者之间的函数关系式,该函数关系可为此类工程的设计提供理论指导。     

Fully Developed Convective Heat Transfer of Power Law Fluids in a Circular Tube

We present a theoretical analysis for fully developed convective beat transfer in a circular tube for power law fluids by assuming that the thermal diffusivity is a function of temperature gradient. The analytical eolution is obtained and the heat transfer behaviour is investigated under a constant heat flux boundary condition. It is shown that the Nusselt number strongly depends on the value of power law index n. The Nusselt number sharply decreases in the range of 0 〈 n 〈 0.1. However, for n 〉 0.5, the Nusselt number decreases monotonically with the increasing n, and for n 〉 20, the values of Nusselt number approach a constant.     

Numerical simulations of electro-thermo-convection and heat transfer in 2D cavity

In this article the electro-thermo-convective phenomena in a dielectric liquid enclosed in a 2D cavity and subjected to the simultaneous action of an electric field and a thermal gradient is studied. We solved directly the full set of coupled equations of Electro-Hydro-Dynamic (EHD) and energy equation using a finite volume method. In order to characterize the influence of the electric field on heat transfer the liquid is first heated (from a lateral wall) till the thermal steady state is obtained and then the electric potential and injection of electric charge is applied. Two cases of injection are considered: from the lower electrode and from a lateral wall (left or right). The flow pattern and Nusselt number strongly depend on the non-dimensional characteristic parameters: electrical parameter, Rayleigh number, Prandtl number and mobility parameter M. The convective motion passing from a purely thermal convection to a purely electrical convection and the number of electro-thermo-convective rolls patterns are investigated.As a consequence of the analysis of the combined effect of electric and thermal fields on the flow structure and on Nusselt number, we have also evaluated the heat transfer enhancement due to electroconvection. It is shown that the injection of electric charge increases the heat transfer and Nusselt number is independent of Rayleigh number for high enough values of T.     

微方肋冷却系统的换热特性

采用去离子水为冷却介质,对自行设计的不同结构微方肋散热器内的换热特性进行实验研究,结果表明:在进口温度为20 ℃、进口流量为57.225 L/h、底面平均温度为73.4 ℃时,散热器散热量可达2.83106 W/cm2,可以满足当前高热流密度散热需求;当散热面温度一定时,散热量随着散热器进口流量的增加而增加,但增速随散热器底面温度的增加变缓;努谢尔特数随雷诺数的增加而成幂次方增加,常规针肋结构和微针肋结构换热关系式不满足微方肋散热器特性。为了更好地表达微方肋散热器内的换热特性,拟合了微方肋散热器内对流换热关系式。     

格子Boltzmann方法分析加热尺寸和瑞利数对可变形开口腔内自然对流的影响

采用格子Boltzmann方法对可变形腔体内自然对流问题进行数值研究，给出平均努赛尔数的经验关系式.腔体左壁加热长度分为左壁面的整个区域（H）和左壁面的中间区域（0.5H）两种情况，右壁向外界环境开放，上下边界绝热且可以上下移动，以此调节右出口尺寸.主要研究瑞利数（10⁴ ≤ Ra ≤ 10⁶），右出口尺寸（1.0H ≤ L ≤ 2.0H），左壁加热尺寸（L_h=0.5H或L_h=H）对腔体内等温线、流线、局部努塞尔数和平均努赛尔数的影响.结果表明：腔体内换热随着瑞利数的增大越来越强烈，表现为椭圆形准静止区域更加靠近上绝热壁，且热分层厚度逐渐变小，平均努赛尔数增加.而右出口尺寸的增加，对于两种加热尺寸下腔内的换热效果有不同程度影响，其中与加热尺寸为左壁面的全部区域L_h=H相比，加热尺寸为左壁面的中间情况L_h=0.5H时，右侧开口尺寸的增加对换热效果的影响不显著.此外，左壁加热尺寸为0.5H时显示出比加热尺寸为H时更高的平均传热效率.最后，针对不同的加热尺寸，提出加热面平均努赛尔数与Ra数及右壁面开口尺寸L^*之间函数关系的经验预测，拟合效果满足工程实践与设计需要.     

HEAT TRANSFER CHARACTERISTICS OF AN IMPINGING BUTANE/AIR FLAME JET OF LOW REYNOLDS NUMBER

An experimental study was performed to determine the heat transfer characteristics of a premixed butane/air round flame jet, of low Reynolds number, impinging upwards normally on a flat rectangular plate. The effects of the exit Reynolds number and equivalence ratio of the flame jet, and the distance between the nozzle and the impingement plate, on the thermal performance of the jet were examined. The range of Reynolds numbers was selected to cover the laminar to the transitional flow conditions. The investigations were conducted with equivalence ratios corresponding to the fuel-rich, stoichiometric, and fuel-lean conditions. The nozzle-plate distance was varied from 1d to 8d. Within the range of Reynolds numbers investigated, the highest Nusselt numbers were obtained at the equivalence ratio of φ = 0.85 when the nozzle-plate distance was maintained at 5d. At the stoichiometric condition, the highest Nusselt number was obtained at the nozzle- plate distance of 6d. Nondimensional correlations were obtained from the experimental results and presented to predict the maximum Nusselt number and average Nusselt number for laminar flame jets as a function of the nozzle-plate distance, Reynolds number, and equivalence ratio.     

Axisymmetric Stagnation-Point Flow of Nanofluid over a Stretching Surface

This paper investigates the laminar boundary layer flow of nanofluid induced by a radially stretching sheet. Nanofluid model exhibiting Brownian motion and thermophoresis is used. Series solutions for a reduced system of nonlinear ordinary differential equations are obtained by homotopy analysis method (HAM). Comparative study between the HAM solutions and previously published numerical results shows an excellent agreement. Velocity, temperature and mass fraction are displayed for various values of parameters. The local skin friction coefficient, the local Nusselt number and the local Sherwood number are computed. It is observed that the presence of nanoparticles enhances the thermal conductivity of base fluid. It is found that the convective heat transfer coefficient (Nusselt number) is decreased with an increase in concentration of nanoparticles whereas Sherwood number increases when concentration of nanoparticles in the base fluid is increased.     

A numerical study of laminar natural convective heat transfer inside a closed cavity with different aspect ratio

Two-dimensional steady-state laminar natural convection was studied numerically for differentially heated air-filled closed cavity with adiabatic top and bottom walls. The temperature of the left heated wall and cooled right wall was assumed to be constant. The governing equations were iteratively solved using the control volume approach. In this paper, the effects of the Rayleigh number and the aspect ratio were examined. Flow and thermal fields were exhibited by means of streamlines and isotherms, respectively.Variations of the maximum stream function and the average heat transfer coefficient were also shown. The average Nusselt number and was correlated to the Rayleigh number based on curve fitting for each aspect ratio. The investigation covered the range 10⁴ ≤ RA ≤ 10⁷ and is done at Prandtl number equal to 0.693. The result shows the average Nusselt number is the increasing function of Rayleigh number. As the aspect ratio increases, Nusselt number decreases along the hot wall of the cavity. As Rayleigh number increases, Nusselt number increases. Result indicates that at constant aspect ratio, with increase in Rayleigh number the heat transfer rate increases.     

Influence of the Shape of the Orifice on the Local Heat Transfer Distribution of a Flat Plate Impinged by Under-Expanded Sonic Jets

Experiments are carried out for the circular, square, and equilateral triangular orifices of the same contraction ratio in order to explore the heat transfer characteristics for nozzle pressure ratios of 2.36, 3.04, 3.72, 4.4, and 5.08 at z/d = 2, 4, 6, and 8. The presence of vena contracta and shock cells and its locations are estimated from the impact pressure distribution along the axial direction. The pressure drop during the initial expansion is minimal for the equilateral triangular orifice. However, this pressure drop is comparable for the square and circular orifices. The heat transfer characteristics of the circular orifice are nearly axisymmetric, but that for non-circular jets are asymmetric, and the three dimensional effects are clearly observed. This is inferred through the spatial gradients of the Nusselt number and adiabatic wall temperature distributions. The linear correlations for the average Nusselt number are proposed for the circular, square, and equilateral triangular orifices.     

On the stability of the flow and heat transfer over a moving thin needle with prescribed surface heat flux

The steady flow and heat transfer over a moving thin needle with prescribed surface heat flux is studied. The similarity equations are obtained by using similarity transformation technique. The problem is solved numerically using the boundary value problem solver (bvp4c) in Matlab software. The plots of the skin friction coefficient and the local Nusselt number as well as the velocity and temperature profiles are presented and their behaviors are discussed for different values of the needle size and the velocity ratio parameter. Results show that the decreasing of the needle size enhance the skin friction coefficient and the local Nusselt number on the needle surface. It is found that dual solutions exist (upper and lower branches) for a certain range of the velocity ratio parameter. A stability analysis of the solutions are performed and it shows that the upper branch solution is stable, while the lower branch solution is unstable.     

不连续冷源布局多孔介质内热流耦合LBM数值模拟

为研究不连续冷源边界对内置发热体多孔介质方腔内传热及流动的影响，采用格子Boltzmann方法对REV尺度下多孔介质方腔内的自然对流进行计算，并研究瑞利数（Ra）、达西数（Da）、孔隙度对多孔介质方腔内传热流动的影响.发现Da对方腔内的流体流型影响很大，Da为10^-4时，多孔介质方腔内只有一个涡流，而Da为10^-2时，方腔内有两个涡流.增大Ra、Da、孔隙度可以提高冷源壁面的平均努赛尔数（Nu），增强散热效果，孔隙率对平均Nu影响程度和Da的大小有关.当冷源布置在壁面上方，壁面的平均Nu随Ra的增加剧烈变化，方腔处于高Ra条件下时，将冷源布置在边界的上方可以提高散热效果.6种布置方案中Case 6的散热效果最好.     

裂解吸热反应对乙烷超临界传热的影响

本文对圆管内超临界状态下乙烷的湍流传热过程进行了数值模拟,分析了乙烷高温裂解对壁面热流密度、平均温度以及对流换热Nusselt数的影响。计算结果表明:考虑乙烷高温裂解吸热反应,管内平均温度会降低,出口处温度降低可达135 K;壁面热流密度则会显著增加,热出口处可增加近2倍;热裂解反应改善了超临界对流换热效果,Nusselt数可提高约20%。在本文计算条件下,经典的对流换热关系式可准确适用于不考虑热裂解反应的情况,而对于考虑裂解吸热反应的超临界乙烷传热情况则误差较大。     

Experimental Analysis of Gas Micro-Convection Through Commercial Microtubes

In this article, an experimental campaign devoted to analyzing the forced micro-convection features of heated gas flows through commercial stainless-steel microtubes having inner diameters of 172 μm and 750 μm is described. The experimental results obtained by heating the microtubes with an imposed uniform heat flux (H-boundary condition) at the external wall, in terms of Nusselt numbers, are compared to the predictions of the classical correlations validated for conventional pipes and to the correlations proposed for gas flows through microtubes under laminar and transitional conditions (100 < Re < 4,000). The cross-sections of the tested microtubes enabled the analysis of the effects of wall axial heat conduction on the Nusselt number. It was observed that the Nusselt number is strongly dependent on the Reynolds number in the laminar regime, and this fact is explained in the article with the effects of wall axial heat conduction and the difficulties in the experimental determination of the right exit bulk temperature of the gas flow, which cannot be ignored in the thermal analysis. The agreement between the Gnielinski correlation and the experimental Nusselt number is poor, especially for low Reynolds numbers, if one uses the average gas bulk temperature, obtained as the arithmetic mean between the inlet and outlet gas bulk temperature, in the definition of the experimental Nusselt number. On the contrary, the agreement with the Gnielinski correlation improves if the local wall-gas temperature difference near the exit of the microtube is used instead. The experimental results presented in the article demonstrate that the criteria for the design of accurate micro-convection tests can be quite different from those for the analysis of forced convection through conventional pipes.     

Fully HOC Scheme for Mixed Convection Flow in a Lid-Driven Cavity Filled with a Nanofluid

A fully higher-order compact (HOC) finite difference scheme on the 9-point two-dimensional (2D) stencil is formulated for solving the steady-state laminar mixed convection flow in a lid-driven inclined square enclosure filled with water-$Al_2O_3$ nanofluid. Two cases are considered depending on the direction of temperature gradient imposed (Case I, top and bottom; Case II, left and right). The developed equations are given in terms of the stream function-vorticity formulation and are non-dimensionalized and then solved numerically by a fourth-order accurate compact finite difference method. Unlike other compact solution procedure in literature for this physical configuration, the present method is fully compact and fully higher-order accurate. The fluid flow, heat transfer and heat transport characteristics were illustrated by streamlines, isotherms and averaged Nusselt number. Comparisons with previously published work are performed and found to be in excellent agreement. A parametric study is conducted and a set of graphical results is presented and discussed to elucidate that significant heat transfer enhancement can be obtained due to the presence of nanoparticles and that this is accentuated by inclination of the enclosure at moderate and large Richardson numbers.     

Correlation for Nusselt number in pure magnetic convection ferrofluid flow in a square cavity by a numerical investigation

Magnetic convection heat transfer in a two-dimensional square cavity induced by magnetic field gradient is investigated numerically using a semi-implicit finite volume method. The side walls of the cavity are heated with different temperatures, the top and bottom walls are isolated, and a permanent magnet is located near the bottom wall. Thermal buoyancy-induced flow is neglected due to the nongravity condition on the plane of the cavity. Conditions for the different values of non-dimensional variables in a variety of ferrofluid properties and magnetic field parameters are studied. Based on this numerical analysis, a general correlation for the overall Nusselt number on the side walls is introduced for a wide range of effective parameters. Results showed that maximum error produced by use of this correlation is about 6 percent.     

Dust in Plasma I. Particle Size and Ion‐Neutral Collision Effects

A fully kinetic self‐consistent model of an absorbing particle immersed in stationary isotropic weakly collisional plasma has been developed. The combined effects of particle size and ion‐neutral charge exchange collisions have been investigated for intermediate regimes, where no analytic theories are available. It is shown that collisional effects related to the ion orbital destruction (presence of extrema in ion flux collected on the particle surface and in particle potential and charge) are important for small particles, while they are totally absent for large particles. The potential distribution around the particle is quite well represented by a Yukawa form, but with an effective screening length that shows different dependences from the gas pressure for small and large particle size. Analytical fitting formulas of particle charge and potential and screening length depending on the particle radius parameter and on the Knudsen number have been obtained (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Calculation of the drag and heat transfer from a sphere in the gas flow in a cylindrical channel

A numerical experiment on the simulation of heat transfer from a sphere to a gas flow in a cylindrical channel in the Stokes and transient flow regimes has been described. Radial and axial profiles of the gas temperature and the dependences of drag coefficient C_d of the body and Nusselt number Nu on Reynolds number Re have been calculated and analyzed. The problem of the influence of the early drag crisis for a sphere on its heat transfer to the gas flow has been considered. The estimation of this phenomenon has shown that the early drag crisis of the sphere in a strongly turbulent flow causes a reduction in heat transfer from the sphere to the gas by three to six times (in approximately the same proportion as for its drag coefficient).     

Air ionization in a near-critical electric field

A set of experimental dependences of the air ionization effective rate on the electric field strength is presented. The concept of the critical breakdown field is discussed. It is indicated that experimental data are quantitatively inconsistent with analytical results based on this concept. This inconsistency is eliminated if the ionization balance takes into account not only dissociative adhesion of electrons to oxygen molecules but also their detachment from the molecules that gained a charge during the charge exchange process. Based on the results obtained, a new physical interpretation of the critical field is suggested. A formula for the effective rate of air ionization in near-critical fields is derived.     

油基钻屑螺纹推进式换热器流动传热特性

对油基钻屑在螺纹推进式换热器内的流动换热过程进行了数值模拟,研究了螺杆转速、油基钻屑雷诺数Re和螺纹截面形状对流动换热的影响。结果表明:随着螺杆转速增加,传热系数、油基钻屑出口温度均增大;同时发现,当雷诺数Re<250时,壳侧Nusselt数随雷诺数Re增大而迅速增大,此后雷诺数对Nusselt数影响较小;Nusselt数随曲率比d_i/D增大而增大。为方便工程设计,利用数值结果给出了油基钻屑的流动换热关系式。     

高电荷态Ar^q＋与Ne碰撞中入射离子电荷交换截面的研究

报道Ar^q Ne(q=8，9，11，12)碰撞体系中多电子转移过程，得到了多组实验测量电荷交换截面数据，讨论入射离子电荷交换截面、反冲离子产生截面与入射离子电荷态、能量以及散射离子电荷态的关系，并且将实验结果与Ar^q Ar碰撞体系进行对比研究。在修正分子库仑过垒模型的基础上，对实验现象做了合理的解释。