板壳式换热器波纹倾角对换热和阻力性能影响

研究了波纹倾角（α=45°，60°，75°）对板壳式换热器单流道内单相流动与换热过程的影响，分析了波纹流道内的速度分布、湍动能分布、压力分布和温度分布；基于范宁摩擦阻力因子f和Nu与Re的关系，提出了板壳式换热器不同波纹倾角下换热特性和流动阻力特性的预测关联式；通过计算不同波纹倾角下的综合性能评价指标（performance evaluation criteria,PEC）和面积质量因子（j/f）综合评价了流动与换热性能。结果表明：波纹倾角是影响圆形板间流体流动形态的因素之一，随着波纹倾角增大，会出现“十字交叉流”向“曲折流”转变；所提出的关联式能够很好地预测板壳式换热器内阻力与换热性能，阻力特性偏差在±14%范围内，换热特性偏差在±7%范围内；α为45°时，j/f较大，流道内阻力相对小；α为60°时，PEC较大，流道内换热性能相对强。     

壁面温度在多参数下对平板边界层的影响研究

为了得到壁面温度在不同来流速度、不同湍流强度条件下对边界层转捩与减阻的影响规律,本文采用Transitionk-kl-ω模型对低来流速度下无压力梯度的光滑平板进行了数值模拟。结果表明,随着来流速度的升高,壁温升高所起到的减阻效果更好,即高来流速度对壁面温度更为敏感。当来流处于中高湍流强度下时,壁温升高能起到推迟转捩的作用,且随着湍流强度的升高,转捩推迟的效果越好,但减阻效果正好相反;当来流处于低湍流强度下时,壁温升高会使得转捩提前发生。壁温升高抑制了边界层内流体的脉动程度,使得层流的稳态不易被破坏,流动更加稳定;同时,壁温升高使得边界层内流体的速度梯度减小,从而降低了壁面摩擦系数,故壁温升高能起到推迟边界层转捩与减阻的作用。     

不同润湿性纳米通道内库埃特流动的模拟

利用非平衡分子动力学模拟方法, 模拟了两无限大平行平板组成的纳米通道内的库埃特流动, 并给出了壁面润湿性和速度对流场密度、速度分布及壁面滑移的影响规律.数值模拟中, 统计系综采用微正则系综, 势能函数选用LJ/126模型, 壁面设为刚性原子壁面, 温度校正使用速度定标法, 牛顿运动方程的求解则采用文莱特算法.结果表明, 纳米通道内流体密度呈对称的衰减振荡分布, 且随壁面润湿性的降低, 振荡幅度减小, 振荡周期保持不变;滑移量随壁面润湿性的提高而降低, 甚至在亲水壁面时出现负滑移现象;随壁面速度的增加滑移速度逐渐增大, 且在流体呈现非线性流动阶段其增幅显著加大.另外, 还发现当壁面设置为超疏水性时, 壁面滑移呈现出随润湿性降低而减小的反常现象, 并基于杨氏方程对其进行了解释.     

微型矩形凹槽的长方体通道的减阻和传热特性

在长方体通道底面沿展向方向间隔设置了微型矩形凹槽，凹槽的深度与边界层尺度相当。采用大涡模拟方法对长方体通道内的流动及传热特性进行研究。数值计算结果表明：在长方体通道内设置的微矩形凹槽可以诱导“突出效应”及二次涡，二次涡的作用类似于微型空气滚动轴承，因而可减小流阻，并使传热性能略有提高。研究表明：微凹槽导致了速度滑移，从而有效降低了通道底部附近流体速度梯度；造成低速条纹变宽，使高低速流体的混合受到抑制。微凹槽内产生的二次涡增加了黏性底层的厚度，且二次涡与微凹槽上方流体之间的滚动摩擦代替了壁面与流体之间的滑动摩擦。与没有布置微型矩形凹槽的长方体通道相比，布置微凹槽的长方体通道可在不影响传热效果的前提下达到6%以上的减阻率。     

表面粗糙度对微细管内气体流动特性的影响

采用了表面粗糙度粘性系数模型对微细管内的气体流动进行数值模拟，以研究微管内壁表面粗糙度对微管内气体流动的影响。运用本文改进的表面粗糙度粘性系数模型，数值模拟与实验数据十分吻合。计算结果表明，进出口压力一定时，表面粗糙度对流场的压力、密度及温度分布的影响不大，但是对速度场的影响十分显著，表面粗糙度使气体流动速度减小，并使壁面附近的速度梯度减小，从而使通过管道的气体质量流量减小，在微管内的气体流动中，表面粗糙度的影响是不能被忽略的。     

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

孤立波斜坡堤堤顶及后坡越浪流数值分析

孤立波往往被用来模拟海啸等浅水大波.运用基于RANS方程、VOF方法以及Goring造波方法建立了二维数值波浪水槽,进行了孤立波在简单斜坡堤上越浪过程的数值模拟.通过数值模拟,研究了不同相对堤顶宽及不同相对波高条件下对孤立波越堤流基本特征,包括堤顶前缘水层厚度及速度、堤顶末端水层厚度及速度、堤顶流厚度及速度沿程分布、后坡流厚度及速度沿程分布.依据数值实验结果,建立了堤顶流最大厚度分布及最大速度分布的经验公式.计算结果表明,堤顶流厚度沿程减小,而堤顶流速度则沿程增大;后坡流厚度先增大,随后再沿程缓慢减小;后坡流速度先减小,随后再沿程缓慢增大.     

二维水平通道内流动沸腾换热的格子Boltzmann模拟

利用格子Boltzmann方法模拟二维水平通道内水的流动沸腾过程,获得不同壁面过热度下流型特点和不同因素对换热过程的影响规律。结果表明,随着壁面过热度升高,流道内流型依次经历从泡状流、弹状流到反环流的转变,平均热流密度和平均换热系数先增大后减小。入口流速降低会使流道内出现受限气泡流,核态沸腾受到抑制。提高入口流速能够有效促进气泡脱离,壁面平均换热系数随入口流速增大而增大,但增长速率有所减小。减小通道宽度有利于汽化现象发生,核态沸腾得到强化,壁面平均换热系数有所提高。     

弓形折流板开孔位置的数值模拟研究

在相同的开孔面积和开孔数目前提下,本文数值研究在开孔位置不同时圆缺高度为0.2D(D为折流板直径)的单弓型折流板换热器的壳程换热与流动阻力问题。数值结果表明:在折流板压强高的“高压区”,均匀地开孔对减小换热器壳程压降效果更好,其综合换热性能的性能评价因子E值也更好;在壳程入口速度为0.1m/s时,通过改变开孔位置,壳程压降可优化3%左右。     

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

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

Turbulent mixed convection flow over a backward-facing step––the effect of the step heights

Effect of the backward-facing step heights on turbulent mixed convection flow along a vertical flat plate is examined experimentally. The step geometry consists of an adiabatic backward-facing step, an upstream wall and a downstream wall. Both the upstream and downstream walls are heated to a uniform and constant temperature. Laser–Doppler velocimeter and cold wire anemometer were used, respectively, to measure simultaneously the time-mean velocity and temperature distributions and their turbulent fluctuations. The experiment was carried out for step heights of 0, 11, and 22 mm, at a free stream air velocity, u_∞, of 0.41 m/s, and a temperature difference, ΔT, of 30 °C between the heated walls and the free stream air. The present results reveal that the turbulence intensity of the streamwise and transverse velocity fluctuations and the intensity of temperature fluctuations downstream of the step increase as the step height increases. Also, it was found that both the reattachment length and the heat transfer rate from the downstream heated wall increase with increasing step height.     

Direct numerical simulation of supersonic pipe flow at moderate Reynolds number

We study compressible turbulent flow in a circular pipe at computationally high Reynolds number. Classical related issues are addressed and discussed in light of the DNS data, including validity of compressibility transformations, velocity/temperature relations, passive scalar statistics, and size of turbulent eddies. Regarding velocity statistics, we find that Huang’s transformation yields excellent universality of the scaled Reynolds stresses distributions, whereas the transformation proposed by Trettel and Larsson (2016) yields better representation of the effects of strong variation of density and viscosity occurring in the buffer layer on the mean velocity distribution. A clear logarithmic layer is recovered in terms of transformed velocity and wall distance coordinates at the higher Reynolds number under scrutiny (Re_τ ≈ 1000), whereas the core part of the flow is found to be characterized by a universal parabolic velocity profile. Based on formal similarity between the streamwise velocity and the passive scalar transport equations, we further propose an extension of the above compressibility transformations to also achieve universality of passive scalar statistics. Analysis of the velocity/temperature relationship provides evidence for quadratic dependence which is very well approximated by the thermal analogy proposed by Zhang et al. (2014). The azimuthal velocity and scalar spectra show an organization very similar to canonical incompressible flow, with a bump-shaped distribution across the flow scales, whose peak increases with the wall distance. We find that the size growth effect is well accounted for through an effective length scale accounting for the local friction velocity and for the local mean shear.     

Temperature and velocity field characteristics of turbulent natural convection in a vertical parallel-plate channel with asymmetric heating

Turbulent natural convection in an asymmetrically heated vertical parallel-plate channel has been studied experimentally and numerically using LDA and CFD. Simultaneous velocity and temperature measurements across the channel at different elevations have been carried out. Three different Ra(b/h) values of 1.91 × 10⁷, 2.74 × 10⁷ and 3.19 × 10⁷ are considered with the channel aspect ratio of b/h = 1/20. Experimental and numerical data are presented in the form of streamwise direction heated wall surface temperature, mean velocity, mean temperature, Reynolds shear stress and turbulent kinetic energy profiles along the channel for one case. These profiles exhibit the flow field development along the channel emphatically. The numerical technique used predicts temperature field fairly well, considerably over-estimating velocity field in the core region.     

Numerical simulation of compressible gas flow and heat transfer in a microchannel surrounded by solid media

In this study, a numerical model is developed to investigate the coupled compressible gas flow and heat transfer in a microchannel surrounded by solid media. To accommodate the varying flow cross-section, the compressible gas flow model is established in a non-orthogonal curvilinear coordinate system. An iterative numerical procedure is employed to solve the coupled heat transfer and gas flow equations. The computer code for the compressible gas flow is first validated against two test problems, and then extended by including the heat conduction in the solid media. The effect of the inlet Mach number on the Nusselt number is examined. It is found that the pressure difference from the pyrolysis front to the heated surface is induced essentially by the gas addition from the channel wall, instead from the pyrolysis front. The necessity of accounting for the gas compressibility is clearly demonstrated when severe heating is applied. The pressure distribution obtained along the channel axial direction is useful for further structural analysis of composite materials.     

Direct numerical simulation of a turbulent axisymmetric jet with buoyancy induced acceleration

Direct numerical simulations of an axisymmetric jet with off-source volumetric heat addition are presented in this paper. The system solved here involves a three-way coupling between velocity, concentration and temperature. The computations are performed on a spherical coordinate system, and application of a traction free boundary condition at the lateral edges allows physical entrainment into the computational domain. The Reynolds and Richardson numbers based on local scales employed in the simulations are 1000 and 12 respectively. A strong effect of heat addition on the jet is apparent. Heating causes acceleration of the jet, and an increased dilution due to an increase in entrainment. Further, the streamwise velocity profile is distorted, and the cross-stream velocity is inward for all radial locations for the heated jet. Interestingly, the maximum temperature is realized off-axis and a short distance upstream of the exit of the heat injection zone (HIZ). The temperature width is intermediate between the scalar and velocity widths in the HIZ. Normalized rms of the concentration and temperature increases in the HIZ, whereas that of streamwise, cross-stream and tangential velocities increases rapidly after decreasing. Both mass flux and entrainment are larger for the heated jet as compared to their unheated counterparts. The buoyancy flux increases monotonically in the HIZ, and subsequently remains constant.     

An experimental study into the effects of streamwise and spanwise acceleration in a turbulent boundary layer

We compare two turbulent boundary layers produced in a low-speed water channel experiment. Both are subjected to an identical streamwise pressure gradient generated via a lateral contraction of the channel, and an additional spanwise pressure gradient is imposed on one of the layers by curving the contraction walls. Despite a relatively high streamwise acceleration, hot-film probe measurements of the mean-velocity distributions show that the Reynolds number increases whilst the coefficient of friction decreases downstream. Visualization of the viscous layers using hydrogen bubbles reveal an increase in the non-dimensional streak spacing in response to the acceleration. Changes in statistical moments of the streamwise velocity near the wall suggest an increased dominance of high-velocity fluctuations. The near-wall streaks and velocity statistics have little sensitivity to the boundary layer three-dimensionality induced by the spanwise pressure gradient, with the boundary-layer crossflow velocity reaching 11 % that of the local freestream velocity.     

Non-invasive measurement of void fraction and liquid temperature in microchannel flow boiling

Past thermometry research for two-phase microfluidic systems made much progress regarding wall temperature distributions, yet the direct measurement of fluid temperature has received little attention. This paper uses a non-invasive two-dye/two-color fluorescent technique to capture fluid temperature along with local liquid fraction in a two-phase microflow generated by injecting air into a heated microchannel. The fluorescent emission of Rhodamine 110 and Rhodamine B, measured with photodiodes, is used to obtain local liquid temperature (±3°C) and void fraction (±2% full-scale) over a temperature range from 45 to 100°C. Arrays of these sensors can significantly expand the set of measurable flow parameters to include bubble/slug frequency, size, velocity, and growth rates in addition to mapping the local liquid temperature and void fraction.     

Measurement in a Zero-Pressure Gradient Turbulent Boundary Layer with Forced Thermal Convection

A subsonic zero-pressure gradient turbulent boundary layer developing on a uniformly heated surface at a Reynolds number in the range of 3, 560?≤?Re _θ ?≤ 5,360 was investigated. Particle-image velocimetry measurements were performed at various positions in the streamwise direction for several wind-tunnel speeds and for different wall excess temperatures to show the thermal convection effects to expand the boundary-layer thickness δ _0.99 and to enlarge the turbulence intensities in the log-law and wake region. The mean velocity profiles are found to be self-preserving. The inclination of large-scale ramp-like vortex packets increases to higher characteristic angles, i.e., the mean angles are enlarged by approximately 5–10°. Hairpin-like vortex structures originating from the near-wall region seem to undergo higher climbing rates in the wall-normal direction causing the above mentioned significant changes in the boundary-layer thickness δ _0.99 and the strongly increased distributions of turbulence intensities in the wake region of the boundary layer. Changes in the distributions of the skewness and flatness of the probability density function (PDF) of the streamwise fluctuations corroborate these findings. The two-point correlation distribution of the streamwise velocity fluctuations R _uu is increased for wall distances y/δ _0.99?=?0.1 to y/δ _0.99?=?0.75 indicating the existence of coherent structures in higher regions of the boundary layer.     

Turbulent boundary layer on a mildly curved convex surface

Extensive single point turbulence measurements made in the boundary layer on a mildly curved heated convex wall show that the turbulence heat fluxes and Stanton number are more sensitive to a change in wall curvature than the Reynolds stresses and skinfriction coefficient, and that downstream, as the flow adjusts to new curved conditions, the St/c _f ratio of Reynolds analogy is appreciably lower than in plane wall flow for the same conditions. Details of the turbulence structure in unheated flow have been documented in an earlier paper; temperature field measurements now described comprise mean temperature distributions, the streamwise variation of wall heat flux, profiles of the temperature variance, transverse and streamwise heat fluxes, and triple correlations. Turbulent diffusion of heat flux is drastically reduced even by mild curvature; changes in the heat fluxes are of the same order as changes in the shear stress, that is, an order of magnitude greater than the ratio of boundary layer thickness to wall radius of curvature. The data include plane flow measurements taken in a developed boundary layer upstream of a change in wall curvature.     

Using digital holographic microscopy for simultaneous measurements of 3D near wall velocity and wall shear stress in a turbulent boundary layer

A digital holographic microscope is used to simultaneously measure the instantaneous 3D flow structure in the inner part of a turbulent boundary layer over a smooth wall, and the spatial distribution of wall shear stresses. The measurements are performed in a fully developed turbulent channel flow within square duct, at a moderately high Reynolds number. The sample volume size is 90 × 145 × 90 wall units, and the spatial resolution of the measurements is 3–8 wall units in streamwise and spanwise directions and one wall unit in the wall-normal direction. The paper describes the data acquisition and analysis procedures, including the particle tracking method and associated method for matching of particle pairs. The uncertainty in velocity is estimated to be better than 1 mm/s, less than 0.05% of the free stream velocity, by comparing the statistics of the normalized velocity divergence to divergence obtained by randomly adding an error of 1 mm/s to the data. Spatial distributions of wall shear stresses are approximated with the least square fit of velocity measurements in the viscous sublayer. Mean flow profiles and statistics of velocity fluctuations agree very well with expectations. Joint probability density distributions of instantaneous spanwise and streamwise wall shear stresses demonstrate the significance of near-wall coherent structures. The near wall 3D flow structures are classified into three groups, the first containing a pair of counter-rotating, quasi streamwise vortices and high streak-like shear stresses; the second group is characterized by multiple streamwise vortices and little variations in wall stress; and the third group has no buffer layer structures.