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

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

横向紊动射流的数值与实验研究进展

横向紊动射流作为流体运动的一种重要类型,广泛存在于如: 燃气轮机气膜冷却、锅炉燃烧室等的燃烧控制, V/STOL(垂直或短距离起落)飞机、废气排放的控制等工程实际应用中.由于射流的存在,增加了流场的复杂性,流场中同时存在射流剪切层涡、马蹄形涡系、反向旋涡对和尾迹涡等4种涡系结构,这对流体力学理论研究具有重要意义.长期以来,研究人员从理论分析、实验测量和数值模拟方面对横向紊动射流进行了大量的研究工作,目前已经认识了流场中的许多流动特性和流动机理.从数值模拟和实验研究两个方面,比较并分析了国内外横向紊动射流研究的现状和研究结果,评述了不同湍流模型以及不同的实验测量方法对横向紊动射流的预测能力,讨论了存在的问题并对该领域的研究方向进行了展望.      

展向喷嘴湍流横向射流大涡模拟及其统计分析

采用大涡模拟方法数值模拟了展向椭圆喷嘴的湍流横向射流,对其大尺度结构的时空演化和湍流脉动速度场的时间序列分析、频谱分析、PDF分析以及时、空截面上的统计平均特性进行分析.结果表明,在射流出口附近的下游核心区中速度脉动剧烈,显现出明显的湍流特征.除了三维涡环脱落、扭曲、变形、摆动所对应频率之外,还存在很宽的湍流基频,它与在喷嘴出口附近产生的三维涡环的时空演化过程密切相关.由于展向椭圆喷嘴的湍流横向射流中的三维涡环快速脱落和强相互作用导致射流尾迹中的强湍流脉动,展向椭圆喷嘴湍流横向射流的PDF空间演化特征结构复杂.在射流核心区的湍流偏应力变化平缓,其统计平均值分布接近左右对称.展向椭圆喷嘴的湍流横向射流脉动速度场具有极为复杂的统计行为,与流向椭圆喷嘴相比具有更好的掺混能力.     

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

RNG湍流模型在冲击射流数值计算中的应用

 将3种版本的RNG,k-ε湍流模型应用于轴对称冲击射流场的数 值模拟. 通过与实验结果的比较，考察了模型系数的改变和非线性附 加项对数值结果的影响，评价了三种模型对冲击射流的数值预测能力.     

表面活性剂减阻溶液湍流流动研究进展

与聚合物添加剂相比,表面活性剂具有寿命长,不受机械力和高温影响发生降解的特点,目前被认为是最具有适用价值的减阻添加剂．针对表面活性剂湍流减阻机理,尽管已经开展了大量的实验研究、理论分析和数值模拟,但仍处于探索阶段,尚未定论．文中对目前国内外有关表面活性剂减阻溶液湍流流动的研究如平均速度、湍流强度、雷诺应力和相关系数等湍流统计量以及取得的成果进行了归纳和总结．目前普遍认为阻力减小是由流动垂直方向的湍流强度受抑和脉动速度分量解耦导致雷诺应力极大降低引起的．分析了目前研究存在的问题,对今后继续开展研究提出了自己的观点．      

添加剂湍流减阻流动与换热研究综述

添加剂湍流减阻是指在液体的管道湍流中添加少量的高分子聚合物或某种表面活性剂从而使湍流阻力大大降低的现象.从其被发现至今,经过近半个世纪的研究(实验研究、理论分析、数值模拟和实际系统的应用研究),尽管对这一现象及其实际应用价值已有了较为深入的认识,但仍有许多方面尚有欠缺,例如对湍流减阻的机理仍然在探索中.本文归纳评述了高分子聚合物或表面活性剂添加剂湍流减阻流动与换热现象的研究现状,从湍流减阻剂的特性、减阻剂的湍流减阻机理、湍流减阻发生时的换热机理、减阻流动速度场分布和换热控制等几个方面综述了添加剂湍流减阻流动与换热特性,并综述了湍流减阻剂在实际工业系统中的应用情况,在对添加剂湍流减阻机理、有湍流减阻发生时的对流换热机理等的理解方面进行了新的总结.      

基于CFD流态分析技术的座便器虹吸管道优化设计

在Fluent软件中应用RNG $k$-$\varepsilon$湍流模型及流体体积函数(VOF)对座便 器内部三维湍流流动进行了雷诺平均N-S方程的数值模拟，得到 了座便器内流场三维流动形态，研究了座便器虹吸管内流动规律，分析了虹吸管形状对流速 分布、压力分布及虹吸性能的影响，以三维湍流场的分析结果为依据，实现了座便器虹吸管 道的优化设计. 通过PIV测试数据验证了三维湍流数值模拟结果的准确性.     

半封闭狭缝湍流冲击射流的数值模拟

将Yakhot和Orszag提出的RNGk-ε模型推广应用于半封闭狭缝冲击射流场的数值模拟，以评价该模型对这种复杂湍流的预测能力。将计算得到的流场平均速度分布、湍流强度分布和流函数分布与标准k—ε模型的预测结果以及相应的实验数据进行了比较，结果表明：RNGk—ε模型的预测结果总体上要好于标准k—ε模型，但与实验值相比，所有预测结果都还存在不同程度的误差，尤其是近壁区和滞止点较远下游处的湍流强度分布。说明RNG模型虽然已在某些湍流的预测中取得了一定的成功，但要定量准确地预测冲击射流场，还必须针对其流动特征对模型加以改进。     

不同参数下环形通道内湍流旋流流动的模拟

应用一种合理考虑湍流一旋流相互作用及湍流脉动各向异性的新的代数ＲｅｙｎｏｌｄＳ应力模型,对环形通道内的湍流旋流流动进行了数值模拟．研究了旋流数、进口轴向速度和内外半径比等参数对环形通道内湍流旋流流动的影响,以及由此产生的流场变化对强化环形通道内传热的作用．     

Computations for a jet impinging obliquely on a flat surface

A SIMPLE-C algorithm and Jones-Launder k-ε two-equation turbulence model are used to simulate a two-dimensional jet impinging obliquely on a flat surface. Both the continuity and momentum equations for the unsteady state are cast into suitable finite difference equations. The pressure, velocity, turbulent kinetic energy and turbulent energy dissipation rate distributions are solved and show good agreement with various experimental data. The calculations show that the flow field structure of the jet impinging obliquely on a flat surface is strongly affected by the oblique impingement angle. The maximum pressure zone of the obliquely impinging jet flow field moves towards the left as the oblique impingement angle is decreased.     

On the analysis of an impinging jet on ground effects

Laser Doppler measurements and flow visualization are presented for a turbulent circular jet emerging into a low-velocity cross-stream and, then, impinging on a flat surface perpendicular to the jet-nozzle axis. The experiments were performed for a Reynolds number based on the jet-exit conditions of 6 × 10⁴, a jet-to-crossflow velocity ratio of 30 and for an impingement height of 5 jet diameters and include mean and turbulent velocity characteristics along the two normal directions contained in planes parallel to the nozzle axis. The results, which have relevance to flows found beneath VSTOL aircraft in ground effect, show the presence of a complex 3-D scarf vortex formed around the impinging jet. In zones where measurement data are not available, the flow details are numerically-visualized using a solution of the finite difference form of the fully threedimensional Reynolds-averaged Navier-Stokes equations, incorporating the turbulence viscosity concept. The turbulent structure of the flow is affected by flow distortion at the impinging zone, which results in an unconventional behaviour of the dimensionless structure parameters that determine the empirical constants in engineering models of turbulence. The relative magnitude of the terms involved in the transport equations for the turbulent stresses is quantified from the experimental data in order to assess the importance of these effects and show the extent to which the turbulent structure of the impingement zone is affected by extra rates of strain.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

Jet impingement onto a cylindrical cavity: consideration of annular nozzle cone angles,and cavity diameter

In the present study, gas jet emerging from an annular nozzle and impinging onto a cylindrical cavity is considered. The geometric configuration of the nozzle is varied in the simulations. Air is used as impinging gas while stainless steel is considered as workpiece material. Reynolds turbulence model is accommodated to account for the turbulence. A numerical scheme employing a control volume approach is used to simulate the flow field. Heat transfer characteristic and shear stress distribution around the cavity are computed. It is found that outer cone angle of the annular nozzle influences the heat transfer rates from the cavity wall. The flow structure around the cavity changes significantly with increasing cavity diameter. Moreover, increasing cavity depth results in stagnation zone moving into the cavity.     

An Investigation of SGS Stress Anisotropy Modeling in Complex Turbulent Flow Fields

An anisotropy-resolving subgrid-scale (SGS) model for large eddy simulation was investigated. Primary attention was given to the predictive performance of the SGS model in the case of complex turbulence with flow impingement and/or flow separation. The SGS model was constructed by combining an isotropic linear eddy-viscosity model with an extra anisotropic term. Since the extra anisotropic term was modeled to prevent undesirable energy transfer between the grid-scale and SGS parts, the model is expected not to seriously affect computational stability. To validate the model performance for complex turbulent flow fields, the SGS model was applied to numerical simulations of a plane impinging jet and 3-D diffuser flow as well as fundamental plane channel flows. The SGS model provided reasonable predictions for these test cases. Furthermore, the predicted SGS stress components were decomposed into linear and anisotropic parts and their roles were investigated in detail. The usefulness of the present anisotropy-resolving SGS model in practical engineering applications was thus described.     

An experimental investigation on flow structures of confined and unconfined impinging air jets

The flow characteristics of both confined and unconfined air jets, impinging normally onto a flat plate have been experimentally investigated. The mean and turbulence velocities, and surface pressures were measured for Reynolds numbers ranging from 30,000 to 50,000 and the nozzle-to-plate spacings in range of 0.2–6. Smoke-wire technique is used to visualize the flow behavior. The effects of Reynolds number, nozzle-to-plate spacing and flow confinement on the flow structure are reported. In the case of confined jet, subatmospheric regions occur on both impingement and confinement surfaces at nozzle-to-plate spacings up to 2 for all Reynolds numbers in consideration and they lie up to nearly the same radial location at both surfaces. However, there is no evidence of the subatmospheric region in unconfined jet. It is concluded that there exists a linkage among the subatmospheric region, turbulence intensity and the peaks in heat transfer coefficients for low spacings in impinging jets.     

Impinging sweeping jet and convective heat transfer on curved surfaces

The application of an impinging sweeping jet, which oscillates periodically with a large angle, to convective heat transfer has received attention owing to its capability to provide a more spatially uniform and enhanced heat removal rate when compared to a steady jet. Herein, we study how the surface curvature affects the heat transfer performance of a sweeping jet and couple it with the representative flow characteristics. Heat transfer measurement and quantitative flow visualization are conducted experimentally for concave and convex surfaces as well as a flat surface. Whereas concave surfaces have a better heat transfer rate than a flat surface, the enhancement of the heat transfer is relatively small for a convex surface. For both concave and convex surfaces, the Nusselt number does not increase monotonically with the curvature magnitude but has a peak for a moderate curvature. The variation in heat transfer performance with the surface curvature is correlated with the phase-averaged velocity profile of the wall jet deflected after an impingement and the turbulence kinetic energy inside the jet. For both concave and convex surfaces, the wall jet becomes thinner than a flat surface in general, which contributes to improved heat transfer. However, whereas the turbulence kinetic energy is significantly larger for a concave surface of a moderate curvature than that of a flat surface, the turbulence kinetic energy for a convex surface is reduced from that of a flat surface, resulting in degradation of the heat transfer performance.     

The low Reynolds number turbulent flow and mixing in a confined impinging jet reactor

Turbulent mixing takes an important role in chemical engineering, especially when the chemical reaction is fast compared to the mixing time. In this context a detailed knowledge of the flow field, the distribution of turbulent kinetic energy (TKE) and its dissipation rate is important, as these quantities are used for many mixing models. For this reason we conduct a direct numerical simulation (DNS) of a confined impinging jet reactor (CIJR) at Re = 500 and Sc = 1. The data is compared with particle image velocimetry (PIV) measurements and the basic flow features match between simulation and experiment. The DNS data is analysed and it is shown that the flow is dominated by a stable vortex in the main mixing duct. High intensities of turbulent kinetic energy and dissipation are found in the impingement zone which decrease rapidly towards the exit of the CIJR. In the whole CIJR the turbulence is not in equilibrium. The strong mixing in the impingement zone leads to a rapid development of a monomodal PDF. Due to the special properties of the flow field, a bimodal PDF is generated in cross-sections downstream the impingement zone, that slowly relaxes under relaminarising conditions. The time required for meso-mixing is dominating the overall mixing performance.     

Drag reduction and heat transfer in surfactant solutions with excess counterion

We present the results of a study of turbulent drag reduction in a small circulating loop using surfactant solutions with excess counterion. In addition, these solutions were used in measurements of heat transfer, both in pipe flow and in an impinging jet. Both frictional drag and heat transfer were reduced in the pipe flow experiments. Measurements of heat transfer in the impinging jet revealed a dependence on the molar concentration ratio of the counterion. When the counterion was added at a molar concentration 30 times higher than that of the surfactant, the resulting surfactant solution did not reduce the rate of heat transfer in the impinging jet. By using this surfactant system in an impinging jet, we show both a reduction in pipe friction and normal heat transfer potential in a circulating heat exchange system. In order to investigate this difference in heat transfer between pipe flows and impinging jet flows, a comparison was made of the wall shear stress between these two flow regimes. The estimated wall shear stress was of the same order in both flows, and thus was not considered to be the primary cause of the difference in heat transfer. It is instead suggested that the micellar structure of the surfactant is influenced by a compressive deformation of the impinging flow in a manner that is different from the shear deformation observed in pipe flow.     

Stereo-PIV measurements of vapor-induced flow modifications in confined jet impingement boiling

A single subcooled jet of water which undergoes boiling upon impingement on a discrete heat source is studied experimentally using time-resolved stereo particle image velocimetry (PIV). The impinging jet issues from a 3.75 mm diameter sharp-edged orifice in a confining orifice plate positioned 4 orifice diameters from the target surface. The behavior at jet Reynolds numbers of 5,000 and 15,000 is compared for a constant jet inlet subcooling of 10 °C. Fluorescent illumination allows for simultaneous imaging of both the flow tracers and the vapor bubbles in the flow. Flow structure, time-averaged velocities, and turbulence statistics are reported for the liquid regions within the confinement gap for a range of heat inputs at both Reynolds numbers, and the effect of the vapor generation on the flow is discussed. Vapor generation from boiling is found to modify the liquid velocities and turbulence fluctuations in the confinement gap. Flow in the confinement gap is dominated by vapor flow, and the vapor bubbles disrupt both the vertical impinging jet and horizontal wall jet flow. Moreover, vapor bubbles are a significant source of turbulence kinetic energy and dissipation, with the bubbly regions above the heated surface experiencing the most intense turbulence modification. Spectral analysis indicates that a Strouhal number of 0.023 is characteristic of the interaction between bubbles and turbulent liquid jets.