统一二阶矩模型用于模拟旋流湍流两相流动

用统一二阶矩模型（ＵＳＭ）模拟了旋流数为０４７和１５的气粒两相流动，并和实验结果以及ｋ ε ｋｐ模型的模拟结果进行了对比．研究结果表明，提高旋流数减小了轴向速度反流区，增大了切向速度似固核区．ＵＳＭ和ｋ ε ｋｐ模型预报旋流数为０４７时的两相速度场差别不大，并都和实验结果接近，但前者预报的旋流数为１５的两相速度场比后者有改进，在两种情况下，前者都能揭示出后者无法预报的两相湍流各向异性规律．     

环形通道内湍流旋流流动的数值模拟

本文应用一种考虑湍流－旋流相互作用及湍流脉动各向异性的新的代数Ｒｅｙｎｏｌｄｓ应力模型,对环形通道内的湍流旋流流动进行了数值模拟,研究了改主为旋流流数,进口轴向速度及半径比等参数对环形通道内湍流流动的影响,以及对强化环形通道内传热的作用。     

Particle statistics in a gas–solid coaxial strongly swirling flow: A direct numerical simulation

A direct numerical simulation of a strongly coaxial swirling particle-laden flow is conducted with reference to a previous experiment. The carrier phase is simulated as a coaxial swirling flow through a short nozzle injecting into a large container. The particle phase is carried by the primary jet, and simulated in the Lagrangian approach. The drag force, slip-shear force and slip-rotation force experienced by particles are calculated. A partial validation of the results is followed. The results are analyzed in Eulerian approach focusing on the statistical behavior of particle motion. The relative importance of the drag, slip-shear and slip-rotation forces under different Stokes numbers is indicated quantitatively. The particle velocity profiles, fluctuations, Reynolds stress, and turbulence intensity are demonstrated and analyzed respectively. An important “choke” behavior for large particles within the mainstream is found and interpreted. Additionally, the patterns of particle distribution and the helical structures of particle motion under different Stokes numbers are demonstrated qualitatively and analyzed quantitatively.     

Characterization of particle-laden, confined swirling flows by phase-doppler anemometry and numerical calculation

The particle dispersion characteristics in a confined swirling flow with a swirl number of approx. 0.5 were studied in detail by performing measurements using phase-Doppler anemometry (PDA) and numerical predictions. A mixture of gas and particles was injected without swirl into the test section, while the swirling airstream was provided through a co-flowing annular inlet. Two cases with different primary jet exit velocities were considered. For these flow conditions, a closed central recirculation bubble was established just downstream of the inlet.

The PDA measurements allowed the correlation between particle size and velocity to be obtained and also the spatial change in the particle size distribution throughout the flow field. For these results, the behaviour of different size classes in the entire particle size spectrum, ranging from about 15 to 80 μm, could be studied, and the response of the particles to the mean flow and the gas turbulence could be characterized. Due to the response characteristics of particles with different diameters to the mean flow and the flow turbulence, a considerable separation of the particles was observed which resulted in a streamwise increase in the particle mean number diameter in the core region of the central recirculation bubble. For the lower particle inlet velocity (i.e. low primary jet exit velocity), this effect is more pronounced, since here the particles have more time to respond to the flow reversal and the swirl velocity component. This also gave a higher mass of recirculating particle material.

The numerical predictions of the gas flow were performed by solving the time-averaged Navier-Stokes equations in connection with the well known kε turbulence model. Although this turbulence model is based on the assumption of isotropic turbulence, the agreement of the calculated mean velocity profiles compared to the measured gas velocities is very good. The gas-phase turbulent kinetic energy, however, is considerably underpredicted in the initial mixing region. The particle dispersion characteristics were calculated by using the Lagrangian approach, where the influence of the particulate phase on the gas flow could be neglected, since only very low mass loadings were considered. The calculated results for the particle mean velocity and the mass flux are also in good agreement with the experiments. Furthermore, the change in the particle mean diameter throughout the flow field was predicted approximately, which shows that the applied simple stochastic dispersion model also gives good results for such very complex flows. The variation of the gas and particle velocity in the primary inlet had a considerable impact on the particle dispersion behaviour in the swirling flow and the particle residence time in the central recirculation bubble, which could be determined from the numerical calculations. For the lower particle inlet velocity, the maximum particle size-dependence residence time within the recirculation region was considerably shifted towards larger particles.     

PTV investigation of phase interaction in dispersed liquid–liquid two-phase turbulent swirling flow

An investigation of dispersed liquid–liquid two-phase turbulent swirling flow in a horizontal pipe is conducted using a particle tracking velocimetry (PTV) technique and a shadow image technique (SIT). Silicone oil with a low specific gravity is used as immiscible droplets. A swirling motion is given to the main flow by an impeller installed in the pipe. Fluorescent tracer particles are applied to flow visualization. Red/green/blue components extracted from color images taken with a digital color CCD camera are used to simultaneously estimate the liquid and droplet velocity vectors. Under a relatively low swirl motion, a large number of droplets with low specific gravity tend to accumulate in the central region of the pipe. With increasing droplet volume fraction, the liquid turbulence intensity in the axial direction increases while that in the wall-normal direction decreases in the central region of the pipe. In addition, the turbulence modification in the present flow is strongly dependent on the droplet Reynolds number; however, the interaction of droplet-induced turbulences is significant due to vortex shedding, particularly at high droplet Reynolds numbers and higher droplet volume fraction.     

双柱单锥型液-液旋流管内流场的激光诊断

应用激光测速仪，对一种双柱单锥型液 液旋流管内的流动结构，进行了全场范围内的多工况流动诊断研究．揭示出其切向速度由内旋流区和外旋流区构成，其中内旋流区中的速度分布符合准强制涡关系，外旋流区中的速度分布符合准自由涡关系；轴向速度由上行流动区和下行流动区构成，两者之间在直管段以零速点作分界，在锥体段则以零速区作过渡并伴随有一定的回流出现，且该过渡区或回流区的大小随锥体截面的收缩而减小，直到进入直管段后消失；各湍流量的分布以管芯处最大向外逐渐减小，沿轴向是直管段中的湍流度大于锥体段中的湍流度，而且湍流度在旋流管内具有各向异性的特性．     

Buoyancy effects on turbulent swirling flows

A three-parameter model of turbulence applicable to free boundary layers has been developed and applied for the prediction of axisymmetric turbulent swirling flows in uniform and stagnant surroundings under the action of buoyancy forces. The turbulent momentum and heat fluxes appearing in the time-averaged equations for the mean motion have been determined from algebraic expressions, derived by neglecting the convection and diffusion terms in the differential transport equations for these quantities, which relate the turbulent fluxes to the kinetic energy of turbulence, k, the dissipation length scale of turbulence, L, and the temperature covariance, T ′². Differential transport equations have been used to determine these latter quantities. The governing equations have been solved using fully implicit finite difference schemes. The turbulence model is capable of reproducing the gross features of pure jet flows, buoyant flows and swirling flows for weak and moderate swirl. The behaviour of a turbulent buoyant swirling jet has been found to depend solely on exit swirl and Froude numbers. The predicted results indicate that the incorporation of buoyancy can cause significant changes in the behaviour of a swirling jet, particularly when the buoyancy strength is high. The jet exhibits similarity behaviour in the initial region for weak swirl and weak buoyancy strengths only, and the asymptotic case of a swirling jet under the action of buoyancy forces is a pure plume in the far field. The predicted results have been found to be in satisfactory agreement with the available experimental data and in good qualitative agreement with other predicted results.     

Numerical study of hysteresis in annular swirling jets with a stepped‐conical nozzle

This study investigates the experimentally observed hysteresis in the mean flow field of an annular swirling jet with a stepped‐conical nozzle. The flow is simulated using the Reynolds‐averaged Navier–Stokes (RANS) approach for incompressible flow with a k–ε and a Reynolds stress transport (RSTM) turbulence model. Four different flow structures are observed depending on the swirl number: ‘closed jet flow’, ‘open jet flow low swirl’, ‘open jet flow high swirl’ and ‘coanda jet flow’. These flow patterns change with varying swirl number and hysteresis at low and intermediate swirl numbers is revealed when increasing and subsequently decreasing the swirl. The influence of the inlet velocity profile on the transitional swirl numbers is investigated. When comparing computational fluid dynamics with experiments, the results show that both turbulence models predict the four different flow structures and the associated hysteresis and multiple solutions at low and intermediate swirl numbers. Therefore, a good agreement exists between experiments and numerics. Copyright © 2006 John Wiley & Sons, Ltd.     

Experimental study of pre-swirl flow effect on the heat transfer process in the entry region of a convergent pipe

Experiments have been performed to study the heat transfer process of swirling flow issued into a heated convergent pipe with a convergent angle of 5° with respect to the pipe axis. A flat vane swirler situated at the entrance of the pipe is used to generate the swirling flow. During the experiments, the Reynolds number ranges from 7970 to 47,820, and the swirl number from 0 to 1.2. It is found that the convergence of the pipe can accelerate the flow which has an effect to suppress the turbulence generated in the flow and reduce the heat transfer. However, in the region of weak swirl (S = 0-0.65), the Nusselt numbers increase with increasing swirl numbers until S = 0.65, where turbulence intensity is expected to be large enough and not suppressible. In the region of strong swirl (S > 0.65), where recirculation flow is expected to be generated in the core of the swirling flow, the heat transfer characteristic can be altered significantly. At very high swirl (S ? 1.0), the accelerated flow in the circumferential direction is expected to be dominant, which leads to suppress the turbulence and reduce the heat transfer. The Nusselt number is found proportional to the swirl number. Correlations of the Nusselt numbers in terms of the swirl number, the Reynolds number and the dimensionless distance are attempted and are very successful in both the weak and the strong swirl regions.     

An experimental study of fluid flow and heat transfer in decaying swirl through a heated annulus

The mean and turbulent structures of turbulent swirling flow in a heated annulus have been measured. Both forced and free vortex swirling flows were generated, and the outer wall of the test section was heated uniformly. The maximum swirl number was 1.39, Reynolds numbers were up to 200000, and heat input was 10.5 kW. Mean and turbulent velocity components, air and wall temperatures, and wall static pressures were all measured. Hot-film techniques were developed to measure turbulence. From these parameters, the flow and temperature fields, pressure distribution, and heat transfer coefficients were determined. The mechanisms of heat transfer were identified.     

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

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

Prediction of Pressure Oscillations in a Premixed Swirl Combustor Flow and Comparison to Measurements

The most common and reliable technique used for flame stabilization of industrial combustors with high thermal loads is the application of strongly swirling flows. In addition to stabilization, swirl flames offer the possibility to influence emission characteristics by simply changing the swirl intensity or the type of swirl generation. Despite of these major advantages, swirling flows tend to evolve flow instabilities, that considerably constitute a significant source of noise. In general, noise generation is substantially enhanced, when such a swirling flow is employed for flames. Thus, the minimization of the resulting noise emissions under conservation of the benefit of high ignition stability is one major design challenge for the development of modern swirl stabilized combustion devices. The present investigation makes an attempt to determine mechanisms and processes to influence the noise generation of flames with underlying swirling flows. Therefore, a new burner has been designed, that offers the possibility to vary geometrical parameters as well as the type of swirl generation, typically applied in industrial devices. Experimental data has been acquired for the isothermal flow as well as swirl flames by means of 3-D-LDV-diagnostics comprising the components of long-time averaged mean and rms-velocities as well as spectrally resolved velocity fluctuations for all components. The noise emission data was acquired with microphone probes resulting in sound pressure levels outside the zone of the perceptible fluid flow. Along to the experiments, numerical simulations using RANS and LES have been carried out for isothermal cases with different burner outlet geometries. The results of the measurements show a distinct rise of the sound pressure level, obtained by changing both the test setup from the isothermal into the flame configuration as well as the geometrical parameters. This is also resembled by the LES simulation results. Furthermore, a physical model has been developed from experiments and verified by the LES simulation, that explains the formation of coherent flow structures and allows to separate their contribution to the overall noise emission from ordinary turbulent noise sources.     

A note on the consistency of hybrid Eulerian/Lagrangian approach to multiphase flows

The aim of the present paper is to introduce and to discuss inconsistency errors that may arise when Eulerian and Lagrangian models are coupled for the simulations of turbulent poly-dispersed two-phase flows. In these hybrid models, two turbulence models are implicitly used at the same time and it is important to check that they are consistent, in spite of their apparent different formulations. This issue is best revealed in the case of very small particles, or tracer-limit particles, where it is assessed that coupling inconsistent turbulence models (Eulerian and Lagrangian) can result in non-physical results, notably for second-order fluid velocity moments. Computations for fluid particles in a turbulent channel flow using several coupling strategies are presented to illustrate this question.     

Considerations for the design of swirl chambers for the cyclone cooling of turbine blades and for other applications with high swirl intensity

The focus of this study lies on turbulent incompressible swirling flows with high swirl intensity. A systematic parameter study is conducted to examine the sensitivity of the mean velocity field in a swirl chamber to changes in the Reynolds number, swirl intensity and channel outlet geometry. The investigated parameter range reflects the typical kinematic flow conditions found in heat transfer applications, such as the cooling of the turbine blade known as cyclone cooling. These applications require a swirl intensity, which is typically much higher than necessary for vortex breakdown. The resulting flows are known as flow regime II and III. In comparison to flow regime I, which denotes a swirling flow without vortex breakdown, these flow regimes are characterized by a subcritical behavior. In this context, subcritical means that the flow is affected by the downstream channel section. Based on mean velocity field measurements in various swirl chamber configurations, it is shown that flow regime III is particularly sensitive to these effects. The channel outlet geometry becomes a determining parameter and, therefore, small changes at the outlet can produce entirely different flow patterns in the swirl chamber. In contrast, flow regime II, as well as flow regime I and axial channel flows, are much less sensitive to changes at the channel outlet. The knowledge about the sensitivity of the flow in different flow regimes is highly relevant for the design of a cyclone cooling system. Cooling systems employing flow regime III can result in a weakly robust flow system that may change completely over the operating range. As a remedy, the swirl intensity needs to be decreased so that flow regime III cannot be reached, which, however, reduces the maximum achievable heat transfer in the cooling system. Alternatively, the flow has to transition back from flow regime III to flow regime II or I before the flow leaves the swirl chamber. Two practical methods are presented. These findings can be directly applied in the design processes of future cyclone cooling systems, and other applications of swirling flow.     

Cubic eddy‐viscosity turbulence models for strongly swirling confined flows with variable density

An investigation on the predictive performance of cubic eddy‐viscosity turbulence models for strongly swirling confined flows with variable density is presented. Comparisons of the prediction with the experiments show some improvements of cubic models over the linear k–ε model. The linear k–ε model does not contain any mechanism to represent the interaction of swirl and density variation and as a consequence it performs poorly. With appropriate modelling, two‐equation cubic turbulence models can capture the subcritical nature of the flow, represent the azimuthal velocity profiles of combined forced‐free vortex motion, and predict the combined effects of swirl and density variation fairly well. However, the calibration of model coefficients is still a topic of investigation. Further amendments are also needed for the equations of k and ε to take into account the effects of swirl and density gradients correctly. Copyright © 2004 John Wiley & Sons, Ltd.     

Axial plus tangential entry swirling jet

This paper presents an experimental investigation on swirling jets with well-defined initial conditions. The axial, radial, and azimuthal velocity components, with their respective fluctuations were measured using high spatial–resolution particle image velocimetry. These detailed measurements allow the initial conditions of the swirling jets to be established and the jets to be characterized using various swirl number definitions. The significance of each term in the swirl number calculations are quantified, and the effect of the common assumptions and simplifications are examined. The characteristics of the jets in relation to the initial conditions are then investigated and compared with the previous studies using similar characterization parameters. Jets with Reynolds number of approximately 5700 and swirl conditions ranging from a non-swirling reference case to high swirl are studied. General properties of swirling jets such as higher spreading rate, higher centerline velocity decay, and higher turbulence level are observed. When the degree of swirl is sufficiently high, vortex breakdown occurs. A swirl number of 0.94 is recorded for a high swirl case prior to vortex breakdown, much higher than the critical swirl number reported in the literature. This behavior is attributed to the effect of the initial conditions on the swirl number calculation.     

Turbulent pipe flow downstream a 90° pipe bend with and without superimposed swirl

In the present work, the turbulent flow downstream a 90° pipe bend is investigated by means of stereoscopic particle image velocimetry. In particular, the three dimensional flow field at the exit of the curved pipe is documented for non-swirling and swirling flow conditions, with the latter being generated through a unique axially rotating pipe flow facility. The non-swirling flow was examined through snapshot proper orthogonal decomposition (POD) with the aim to reveal the unsteady behaviour of the Dean vortices under turbulent flow conditions, the so-called “swirl-switching” phenomenon. In respect to the swirling turbulent pipe flow, covering a wide range of swirl strengths, POD has been employed to study the effect of varying strength of swirl on the Dean vortices as well as the interplay of swirling motion and Dean cells. Furthermore, the visualised large scale structures in turbulent swirling flows through the bend are found to incline and tear up with increasing swirl intensity. The present time-resolved, three component, experimental velocity field data will provide a unique and useful database for future studies; in particular for the CFD community.     

旋流燃烧室内湍流燃烧速度场的实验研究

建立了采用分级进风方式的同轴射流旋流燃烧室实验装置，选用耐高温的氧化铝细粉作为示踪粒子，实现了用三维激光粒子动态分析仪(PDA)测量湍流旋流燃烧的热态瞬时速度场．在分级进风比率和旋流致不同的3组实验工况条件下，得到了气体时均轴向与切向速度、轴向与切向脉动速度均方根值和轴向—切向脉动速度二阶关联量的分布．     

An experimental study and CFD analysis towards heat transfer and fluid flow characteristics of decaying swirl pipe flow generated by axial vanes

In this presentation, influences of axial vane swirler on heat transfer augmentation and fluid flow are investigated both experimentally and numerically. The swirl generator is installed at the inlet of the annular duct to generate decaying swirling pipe flow. Three different blade angels of 30°, 45° and 60° were examined. Meanwhile, flow rate was adjusted at Reynolds numbers ranging from 10000 to 30000. Study has been done under uniform heat flux condition and air was used as working fluid. Experimental results confirm that the use of vane swirler leads to a higher heat transfer compared with those obtained from plain tubes. Depending on blade angle, overall Nusselt augmentation is found from 50% to 110% while friction factor increases by the range of 90–500%. Thermal Performance evaluation has been done for test section and test section together with swirler. In both cases, thermal performance increases as vane angle is raised and decreases by growth of Re number. When increasing the blade angle, higher decay rate has been observed for local Nusselt number. In CFD analysis, time-averaged governing equations were solved numerically and RSM model was applied as the turbulence model. Here, the simulation results of axial and tangential velocities, turbulent kinetic energy, wall stresses and swirl intensity are provided. They illustrate the effect of swirling pattern on mean flow and turbulence structure, as well as on improving heat transfer enhancement in the annular duct.