Thermo-mechanical modeling of turbulent heat transfer in gas–solid flows including particle collisions

A thermo-mechanical turbulence model is developed and used for predicting heat transfer in a gas–solid flow through a vertical pipe with constant wall heat flux. The new four-way interaction model makes use of the thermal k_θ–τ_θ equations, in addition to the hydrodynamic k–τ transport, and accounts for the particle–particle and particle–wall collisions through a Eulerian/Lagrangian formulation. The simulation results indicate that the level of thermal turbulence intensity and the heat transfer are strongly affected by the particle collisions. Inter-particle collisions attenuate the thermal turbulence intensity near the wall but somewhat amplify the temperature fluctuations in the pipe core region. The hydrodynamic-to-thermal times-scale ratio and the turbulent Prandtl number in the region near the wall increase due to the inter-particle collisions. The results also show that the use of a constant or the single-phase gas turbulent Prandtl number produces error in the thermal eddy diffusivity and thermal turbulent intensity fields. Simulation results also indicate that the inter-particle contact heat conduction during collision has no significant effect in the range of Reynolds number and particle diameter studied.     

可压缩各向同性衰减湍流直接数值模拟研究

采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性 和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等. 研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大 尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量(压缩)耗散率所占比率也显著增长. 标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高(p ≥ 5)时,纵向速度结构函数的扩展自相似性则不再成立. 对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的 间歇性要强于横向湍流脉动.     

金属丝壁面加热湍流标度律的实验研究

通过在固壁表面的平板湍流边界层沿流向平行放置若干通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的流向涡结构。对壁湍流小尺度结构标度律统计特性的研究表明,金属丝加热后形成的规则流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,增强了湍流小尺度结构的层次结构相似性,减小了壁湍流中小尺度结构的间歇性和奇异性,抑制了壁湍流中奇异的湍涡结构。     

The experimental results of micro-structure in grid-produced turbulence

The variation of main turbulent quantities in an isotropic turbulent flow, such as the decay of turbulent energy and the variation of Taylor microscale of turbulence with time are obtained, by employing a hot-wire anemometer and a nearly isotropic turbulent flow which is produced by a gridscreen located at the entrance of the test section in a low-level turbulence and low-speed wind tunnel in Peking University. The experimental results of the decay of turbulent energy and the variation of Taylor microscale of turbulence with time at the whole period from initial to final stage, normalized in an non-dimensional form, are consistent quite well with the computational results by the theory of the statistical vorticity structure^[1]. The experimental results presented in this paper also agree with Townsend's results obtained in earlier years^[2] as well as with Bennett's in the seventy's^[3].     

Interaction of a normal shock wave with a compressible turbulent flow

A speckle photographic method, which is sensitive to changes of gradients in fluid density, is applied for analyzing a compressible turbulent air flow with density fluctuations. Spatial correlation coefficients, turbulent length scales, and energy spectra are determined under the assumption of homogeneous isotropic turbulence. The experiments are performed in a shock tube where the flow is passed through a turbulence grid. Measurements are taken before and after the turbulent regime interacts with the normal shock wave reflected from the tube's end wall. Amplification of the turbulence intensity by the shock interaction process is verified quantitatively and is shown to be restricted to the lower wave numbers in the spectrum.A version of this paper was presented at the 11th Symposium on Turbulence, University of Missouri-Rolla, Oct. 17–19, 1988.To Professor Dr.-Ing. Klaus Gersten on the occasion of his 60th birthday     

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

Symmetries, Invariance and Scaling-Laws in Inhomogeneous Turbulent Shear Flows

An approach to derive turbulent scaling laws based on symmetry analysis is presented. It unifies a large set of scaling laws for the mean velocity of stationary parallel turbulent shear flows. The approach is derived from the Reynolds averaged Navier–Stokes equations, the fluctuation equations, and the velocity product equations, which are the dyad product of the velocity fluctuations with the equations for the velocity fluctuations. For the plane case the results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the centre of a Couette flow and in the centre of a rotating channel flow, and a new exponential mean velocity profile that is found in the mid-wake region of high Reynolds number flat-plate boundary layers. The algebraic scaling law is confirmed in both the centre and the near wall regions in both experimental and DNS data of turbulent channel flows. For a non-rotating and a moderately rotating pipe about its axis an algebraic law was found for the axial and the azimuthal velocity near the pipe-axis with both laws having equal scaling exponents. In case of a rapidly rotating pipe, a new logarithmic scaling law for the axial velocity is developed. The key elements of the entire analysis are two scaling symmetries and Galilean invariance. Combining the scaling symmetries leads to the variety of different scaling laws. Galilean invariance is crucial for all of them. It has been demonstrated that two-equation models such as the k– model are not consistent with most of the new turbulent scaling laws.     

网格湍流微结构的实验研究

本文试验研究了网格湍流从前期到后期整个连续衰变过程即湍能和Taylor微尺度随时间的变化规律,以及高阶速度相关系数的变化规律,试验结果是在一个低湍流度、低速风洞内,用TSI热线风速仪测得的,而拟均匀各向同性湍流是用在风洞试验段入口处加网格产生的,本文的试验结果与文献[1]提出的涡旋结构理论的计算结果做了比较,发现理论计算的和曲线与本文实测值非常吻合,本文的实测结果与Townsendt早年的试验以及Beanett七十年代末的试验也做了比较。结果表明,这些试验结果彼此也很一致,因而,所有这些试验结果与理论计算值都相互获得了验证。     

EXPERIMENTAL STUDY OF MEASUREMENT FOR DISSIPATION RATE SCALING EXPONENT IN HEATED WALL TURBULENCE

Experimental investigations have been devoted to the study of scaling law of coarse-grained dissipation rate structure function for velocity and temperature fluctuation of non-isotropic and inhomogeneous turbulent flows at moderate Reynolds number. Much attention has been paid to the case of turbulent boundary layer, which is typically the non-istropic and inhomogeneous trubulence because of the dynamically important existence of organized coherent structure burst process in the near wall region . Longitudinal velocity and temperature have been measured at different vertical positions in turbulent boundary layer over a heated and unheated flat plate in a wind tunnel using hot wire anemometer. The influence of non-isotropy and inhomogeneity and heating the wall on the scaling law of the dissipation rate structure function is studied because of the existence of organized coherent structure burst process in the near wall region . The scaling law of coarse-grained dissipation rate structure function is foun     

DNS of heat transfer in transitional, accelerated boundary layer flow over a flat plate affected by free-stream fluctuations

Direct numerical simulations (DNS) of flow over and heat transfer from a flat plate affected by free-stream fluctuations were performed. A contoured upper wall was employed to generate a favourable streamwise pressure gradient along a large portion of the flat plate. The free-stream fluctuations originated from a separate LES of isotropic turbulence in a box. In the laminar portions of the accelerating boundary layer flow the formation of streaks was observed to induce an increase in heat transfer by the exchange of hot fluid near the surface of the plate and cold fluid from the free-stream. In the regions where the streamwise pressure gradient was only mildly favourable, intermittent turbulent spots were detected which relaminarised downstream as the streamwise pressure gradient became stronger. The relaminarisation of the turbulent spots was reflected by a slight decrease in the friction coefficient, which converged to its laminar value in the region where the streamwise pressure gradient was strongest.     

Predicting Turbulent Spectra in Drag-reduced Flows

In the present work we describe how turbulent skin-friction drag reduction obtained through near-wall turbulence manipulation modifies the spectral content of turbulent fluctuations and Reynolds shear stress with focus on the largest scales. Direct Numerical Simulations (DNS) of turbulent channels up to Re _τ = 1000 are performed in which drag reduction is achieved either via artificially removing wall-normal turbulent fluctuations in the vicinity of the wall or via streamwise-travelling waves of spanwise wall velocity. This near-wall turbulence manipulation is shown to modify turbulent spectra in a broad range of scales throughout the whole channel. Above the buffer layer, the observed changes can be predicted, exploiting the vertical shift of the logarithmic portion of the mean streamwise velocity profile, which is a classic performance measure for wall roughness or drag-reducing riblets. A simple model is developed for predicting the large-scale contribution to turbulent fluctuation and Reynolds shear stress spectra in drag-reduced turbulent channels in which a flow control acts at the wall. Any drag-reducing control that successfully interacts with large scales should deviate from the predictions of the present model, making it a useful benchmark for assessing the capability of a control to affect large scales directly.     

Response of the streaks, the active and passive eddies in an unsteady channel flow

Spanwise space–time correlations of the wall shear stress and the longitudinal velocity fluctuations in the low buffer layer of an unsteady channel flow are reported. The imposed amplitude is 20% of the centerline velocity and the imposed frequency covers a large range going from the quasi-steady limit to the bursting frequency of the corresponding steady flow. The unsteady spanwise correlation coefficient is investigated both through its own modulation characteristics (amplitude and phase shifts) and those of the resulting streak spacing. A good correspondence is found between the modulation of the streak spacing and that of the ejection period. The data is further analyzed by temporal filtering of the wall shear stress and streamwise velocity fluctuations. It is shown that the large outer-layer structures play a “passive” role in the unsteady response of the near wall turbulence. The inner wall eddies, in return, are amply responsible for the unsteady reaction of both the turbulent wall shear stress and the streamwise velocity intensities in the buffer layer.     

壁湍流边界层奇异标度律的实验研究

采用热线风速仪对平板湍流边界层的流向速度进行测量,用速度结构函数研究不同尺度结构标度律的变化规律,结果显示小尺度区的概率密度曲线尾部明显偏离高斯型,说明高幅值间歇性事件占的份额较多;惯性子区的曲线向高斯型靠近,间歇性事件所占份额减少;大尺度结构的曲线趋于高斯型,间歇性事件所占份额最小。在耗散区、惯性子区和较大的尺度结构区存在大小不同的绝对标度指数,越靠近壁面这些区域的标度指数均越偏离p/3而逐渐变小。绝对标度指数与边界层位置有关,在缓冲层各阶标度指数与线性标度律偏差很大,显示较强的奇异性,当过渡到对数层及外区,标度指数逐渐增大,接近均匀各向同性湍流的状态。缓冲层、对数层及外区具有各异的绝对标度指数增长率,与各层的不同湍流结构特征和运动形式有关。     

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.     

Rotation effects on a fully-developed turbulent pipe flow

A fully-developed turbulent pipe flow is allowed to pass through a rotating pipe section, whose axis of rotation coincides with the pipe axis. At the exit end of the rotating section, the flow passes into a stationary pipe. As a result of the relaxation of surface rotation, the turbulent flow near the pipe wall is affected by extra turbulence production created by the large circumferential shear strain set up by the rapid decrease of the rotational velocity to zero at the wall. However, the flow in the most part of the pipe is absent of this extra turbulence production because the circumferential strain is zero as a result of the solid-body rotation imparted to the flow by the rotating pipe section. The combined effect of these two phenomena on the flow is investigated in detail using hot-wire anemometry techniques. Both mean and turbulence fields are measured, together with the wall shear and the turbulent burst behavior at the wall. A number of experiments at different rotational speeds are carried out. Therefore, the effects of rotation on the behavior of wall shear, turbulent burst at the wall, turbulence production and the near-wall flow can be documented and analysed in detail.     

Experimental characterization of non-premixed turbulent jet propane flames

This paper reports an experimental study conducted on turbulent jet propane flames aiming at further understanding of turbulent structure in non-premixed slow-chemistry combustion systems. Measurements of mean and fluctuating velocity and temperature fields, mean concentration of major chemical species, correlation between velocity and temperature fluctuations, and dissipation of temperature fluctuations are reported in a turbulent round jet non-premixed propane flame, Re=20 400 and 37 600, issuing vertically in still air. The experimental conditions were designed to provide a complete definition of the upstream boundary conditions in the measurement domain for the purpose of validating computational models. The measured data depicts useful flow field information for describing turbulent non-premixed slow-chemistry flames. Velocity–temperature correlation measurements show turbulent heat fluxes tended to be restricted to the mixing layer where large temperature gradients occurred. Observations of non-gradient diffusion of heat at x/D=10 were verified. Temperature fluctuation dissipation, χ, showed the highest values in the shear layer, where the variance of temperature fluctuations was maximum and combustion occurred. The isotropy between the temperature dissipation in the radial and tangential directions was confirmed. By contrast, the observed anisotropy between axial and radial directions of dissipation suggests the influence of large structures in the entrainment shear layer on the production of temperature fluctuations in the flame region. The value of the normalized scalar dissipation at the stoichiometric mixture fraction surface, χ_st, was calculated, and ranges between 2 and 4 s⁻¹. The measured data were used to estimate the budgets in the balance equations for turbulent kinetic energy, Reynolds shear stresses, turbulent heat flux and temperature variance, quantifying the mechanisms involved in the generation of turbulence as well as in the transport of the temperature.     

Large-eddy simulation in a mixing tee junction: High-order turbulent statistics analysis

This study analyses the mixing and thermal fluctuations induced in a mixing tee junction with circular cross-sections when cold water flowing in a pipe is joined by hot water from a branch pipe. This configuration is representative of industrial piping systems in which temperature fluctuations in the fluid may cause thermal fatigue damage on the walls. Implicit large-eddy simulations (LES) are performed for equal inflow rates corresponding to a bulk Reynolds number Re = 39,080. Two different thermal boundary conditions are studied for the pipe walls; an insulating adiabatic boundary and a conducting steel wall boundary. The predicted flow structures show a satisfactory agreement with the literature. The velocity and thermal fields (including high-order statistics) are not affected by the heat transfer with the steel walls. However, predicted thermal fluctuations at the boundary are not the same between the flow and the solid, showing that solid thermal fluctuations cannot be predicted by the knowledge of the fluid thermal fluctuations alone. The analysis of high-order turbulent statistics provides a better understanding of the turbulence features. In particular, the budgets of the turbulent kinetic energy and temperature variance allows a comparative analysis of dissipation, production and transport terms. It is found that the turbulent transport term is an important term that acts to balance the production. We therefore use a priori tests to evaluate three different models for the triple correlation.     

Chaotic advection in the near-wake of a low-rise building

Properties of the mean-streamlines of a turbulent wake flow behind a low-rise building have been investigated in numerical simulations. For some near wind directions, the solutions to the equations for the mean-streamlines in the recirculating region of flow are found to be chaotic. This is shown to be analogous to the occurrence of chaotic advection in laminar flows and consistent with the current understanding of chaotic dynamics in Hamiltonian systems. The transitions between regular and chaotic solutions for the mean-streamlines are found to occur abruptly over incremental changes, Δθ<2°, in the mean wind direction, θ. These transitions are shown to be accompanied by abrupt changes in the characteristics of plumes formed from tracer-particles undergoing turbulent dispersion. Small fluctuations in the mean wind direction at the transition between regular and chaotic solutions for the mean-streamlines are predicted to result in strongly intermittent turbulent dispersion.     

利用子波分析对平壁湍流猝发现象的研究

利用槽道湍流直接数值模拟的数据库,采用子波分析的方法。对平壁湍流猝发现象的多尺度特性进行了研究,在不同惊讶上对猝发平均周期进行了统计,并利用局部标度指数研究了猝发过程的奇异性。     

Some problems and issues in geophysical turbulence

Velocity fluctuations in the large scales of the atmosphere's meso-scale have turbulent characteristics of random fluctuations and a scale-size distribution near k^−5/3 (Gage, 1979. J. Atmos. Sci. 36, 1950–1954; Lilly, 1983. J. Atmos. Sci. 40, 749–761; Lilly et al., 1998. Theoret. Comput. Fluid Dyn. 11, 139–153). Explanations of this motion field have ranged from inverse cascading quasi-geostrophic (i.e. quasi-two-dimensional) turbulence, to gravity waves (VanZandt, 1982. Geophys. Res. Lett. 9, 575–578). We describe efforts to relate observational spectra to various theories ranging from quasi-geostrophic turbulence to gravity waves. We note that at the larger scales quasi-geostrophic theory may suffice, but at smaller scales, a quasi-geostrophic explanation becomes untenable because the importance of rotation becomes progressively weaker as scales of the flow becomes smaller (the Rossby number approaches unity). We then discuss numerical simulations designed to discriminate between alternative explanations of the flow. Several simulations are reviewed, starting with those of Herring and Métais (1989. J. Fluid Mech. 202, 97–115), and finally those described by Lilly et al. (1998. Theoret. Comput. Fluid Dyn. 11, 139–153).