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

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

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.     

Numerical predictions of the bifurcation of confined swirling flows

The bifurcation of confined swirling flows was numerically investigated by employing both the k-? and algebraic stress turbulence models. Depending upon the branch solution examined, dual flow patterns were predicted at certain swirl levels. In the lower-branch solution which is obtained by gradually increasing the swirl level from a low-swirl flow, the flow changes with increasing swirl number from the low-swirl flow pattern to a high-swirl flow pattern. In the upper-branch solution which is acquired by gradually decreasing the swirl level from a high-swirl flow, on the other hand, the flow can maintain itself in the high-swirl flow pattern at the swirl levels where it exhibits the low-swirl flow pattern in the lower branch. The bifurcation of confined swirling flows was predicted with either the k-? model or the algebraic stress model being employed. Both the k-? and algebraic stress models result in comparable and sufficiently good predictions for confined swirling flows if high-order numerical schemes are used. The reported poor performance of the k-? model was clarified to be mainly attributable to the occurrence of the bifurcation and the use of low-order numerical schemes.     

On the effect of the nozzle design on the performances of gas–liquid cylindrical cyclone separators

This paper constitutes an experimental study of the separation performances of a gas–liquid cylindrical cyclone (GLCC) separator that interests the oil industry. The global hydrodynamics behavior in the GLCC is characterized by flow visualization under various inflow operating conditions. The effect of the inlet nozzle design on the performances of the separator is studied by using three different nozzles, and it proves to be a key parameter. With an insufficient nozzle restriction, low swirl intensity is imparted to the flow. Due to inadequate centrifugal effects, liquid is prematurely carried over by the gas as flooding occurs in the separator upper part. High amounts of gas are also carried under by the liquid stream. On the other hand, with a too severe nozzle convergence, the important drag applied by the gas leads to liquid “short circuiting” the cyclone toward the gas outlet. In addition to the nozzle design, the separator performances are influenced by phenomena such as liquid bridging or the occurrence of the slug flow regime at the cyclone inlet. This paper leads to a better understanding of the links between the hydrodynamics in the GLCC and its operational limits, which is necessary to enable reliable scaling up tools.     

Numerical simulation of a dense solid particle flow inside a cyclone separator using the hybrid Euler–Lagrange approach

This paper presents a numerical simulation of the flow inside a cyclone separator at high particle loads. The gas and gas–particle flows were analyzed using a commercial computational fluid dynamics code. The turbulence effects inside the separator were modeled using the Reynolds stress model. The two phase gas–solid particles flow was modeled using a hybrid Euler–Lagrange approach, which accounts for the four-way coupling between phases. The simulations were performed for three inlet velocities of the gaseous phase and several cyclone mass particle loadings. Moreover, the influences of several submodel parameters on the calculated results were investigated. The obtained results were compared against experimental data collected at the in-house experimental rig. The cyclone pressure drop evaluated numerically underpredicts the measured values. The possible reason of this discrepancies was disused.     

Fluid Dynamical Analysis of Atmospheric Reacting and Isothermal Swirling Flows

A series of unconfined swirling premixed natural gas/air flames was investigated. Reynolds-numbers spanned from 10,000 to 42,300. Respective isothermal flows were studied additionally to gain insight into changes of fluid dynamical features caused by combustion. Statistical moments, Reynolds-stresses, temporal time scales, spatial length scales, and power spectral densities were deduced from one- and two-point laser Doppler velocimetry (LDV) data. Properties of the turbulent flows and dependencies on Reynolds-number, swirl number, and chemical reactions are discussed. Most distinct differences between combusting and isothermal flows were precessing vortex cores (PVC) occurring only for the latter cases. The study is aimed to serve as a database of a generic flame geometry featuring important characteristics of industrial applications for validation of numerical simulations. Therefore, nozzle exit profiles as important inlet conditions to numerical simulations are thoroughly documented.     

Parametric Study of Alternating Flow Patterns in Non-Reacting Multiple-Swirl Flows

Multiple nozzle combustors, under certain conditions, may result in flowfields that differ between nozzles in an alternating pattern. Previous work has provided some clues on the parameters which govern the appearance of this behavior, but there is a lack of systematic studies. A series of non-reacting simulations of adjacent swirling flows is used to investigate the effect of nozzle exit flare angle and swirl number on the presence of the alternating flow pattern. Two-nozzle simulations are shown to accurately predict if an asymmetric flow characteristic appears and are therefore used in the parametric investigation. Alternating flow patterns are predicted at nozzle exit flare angles of 105 degrees (for a swirl number of 0.79) and 120 degrees (for a swirl number of 0.69 and 0.79). Under conditions close to the stability boundary between symmetric and asymmetric flows, the nozzle exit flare and increased swirl number push the shear layers against the dome wall so that the flows between each nozzle are largely opposite in direction. An increase in nozzle exit flare above 120^° results in separated flows exiting from the inlet and a return to a symmetric flow state. This is consistent with a proposed physical mechanism based on hydrodynamic stability in turbulent opposed jets.     

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

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

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.     

强旋湍流气粒两相流动的PDPA研究

采用相多普勒颗粒分析仪（ＰＤＰＡ）对切向进气，轴向缩口出口的旋风筒内强旋单相和气粒两相流动进行了实验研究，给出了强旋流场中，两相湍流的运动及相互作用规律     

Large Eddy Simulation of Turbulent Confined Highly Swirling Annular Flows

The present paper discusses the Large Eddy Simulation of a confined non-reacting annular swirling jet. The configuration corresponds to a well investigated flow studied experimentally by Sheen (1993). The flow field is characterised by a high swirl number resulting in relatively complex features. The challenging behaviour of the flow is governed by the interaction of several recirculation zones. The central recirculation zone formed by the swirling jet is strongly affected by the cylindrical centre body which acts as a bluff body. The flow features coherent structures such as Precessing Vortex Cores (PVCs), which create regions with high velocity fluctuations. The simulations presented comprise a detailed investigation of the parameters controlling the inert flow and a thorough comparison with the experimental data.     

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.     

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.     

Experimental investigation of helical structures in swirling flows

In this investigation the flow in a generic swirl tube with a tangential double-inlet swirl generator and variable exit orifices was experimentally investigated. Using magnetic resonance velocimetry (MRV) three-dimensional, three-component velocity fields were measured for two different Reynolds numbers: 10,000 and 15,000, and for three different exit orifices. The swirl generator had a fixed geometry producing an initial swirl number of 1.6 for all cases. One major observation is the occurrence of a three-layered flow structure. An annular main flow was surrounded by a recirculation zone, as reported in previous literature. However, this recirculation zone – also of an annular shape – exhibited a third layer inside: a thin, high speed jet in the center of the tube with the same flow direction as the main flow. Therefore, the conventional classification of swirling flows into ring and recirculation zone, has to be extended by a core zone. This three-layered flow structure develops independently of the exit configuration. Helical structures were observed in the near-wall region for all cases investigated. Applying an eccentric exit orifice results in the occurrence of strong stationary helical structures not only in the near-wall region but also in the center of the tube. The results, deviating significantly from previous results in the literature, underline the need for more detailed research on the topic of cyclone type flows.     

Large-eddy structures of turbulent swirling flows and methane-air swirling diffusion combustion

Turbulent swirling flows and methane-air swirling diffusion combustion are studied by large-eddy simulation (LES) using a Smagorinsky-Lilly subgrid scale turbulence model and a second-order moment (SOM) SGS combustion model, and also by RANS modeling using the Reynolds Stress equation model with the IPCM+wall and IPCM pressure-strain models and SOM combustion model. The LES statistical results for swirling flows give good agreement with the experimental results, indicating that the adopted subgrid-scale turbulence model is suitable for swirling flows. The LES instantaneous results show the complex vortex shedding pattern in swirling flows. The initially formed large vortex structures soon break up in swirling flows. The LES statistical results of combustion modeling are near the experimental results and are as good as the RANS-SOM modeling results. The LES results show that the size and range of large vortex structures in swirling combustion are different from those of isothermal swirling flows, and the chemical reaction is intensified by the large-eddy vortex structures. The project supported by the Special Funds for Major State Basic Research (G-1999-0222-07). The English text was polished by Keren Wang.     

LDV measurements in complex swirling flows,their physics and a database for CFD evaluations

In an earlier paper in this journal (Mocikat et al. in Exp Fluids 34:442–448, 2003) LDV measurements in a simple geometry but for a complex flow have been provided as a database for CFD evaluation purposes. With special inflow devices swirl could now be added to the flow. By changing the exit position of the test section in order to get a non-symmetric flow field, a steady swirling flow without instability induced precessing motions could be established. This flow can be interpreted as a superposition of a swirling motion to an otherwise swirl-free flow by introducing “swirl influence factors” for various aspects of the flow field. With a modified inflow device a periodically unsteady flow with swirl emerged. The turbulence features of this flow are distinctively different from the steady flow case with swirl. For all flows under consideration the three time-averaged components of the velocity vector and all components of the Reynolds stress tensor are measured in selected cross sections and provided as a data base for CFD calculations.     

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

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

Turbulent flow in converging nozzles,part one: boundary layer solution

The boundary layer integral method is used to investigate the development of the turbulent swirling flow at the entrance region of a conical nozzle. The governing equations in the spherical coordinate system are simplified with the boundary layer assumptions and integrated through the boundary layer. The resulting sets of differential equations are then solved by the fourth-order Adams predictor-corrector method. The free vortex and uniform velocity profiles are applied for the tangential and axial velocities at the inlet region, respectively. Due to the lack of experimental data for swirling flows in converging nozzles, the developed model is validated against the numerical simulations. The results of numerical simulations demonstrate the capability of the analytical model in predicting boundary layer parameters such as the boundary layer growth, the shear rate, the boundary layer thickness, and the swirl intensity decay rate for different cone angles. The proposed method introduces a simple and robust procedure to investigate the boundary layer parameters inside the converging geometries.     

静电旋风分离器气相流场的数值模拟

对一内部安装电晕极的切向进口旋风分离器,以三维贴体坐标为基础,应用Bradshaw的修正k-ε湍流模型,用非交错的SIMPLE算法对静电旋风分离器的气相流场进行求解,计算结果与实验结果进行了对比。分析了电晕极的不同安装位置对旋风分离器流场的影响。从流场的角度来看,电晕极安装在筒体与排气管之间并靠近排气管的位置有利于提高静电旋风分离器的分离效率。     

Investigation of the “TECFLAM” Non-premixed Flame Using Large Eddy Simulation and Proper Orthogonal Decomposition

The confined turbulent swirling non-premixed TECFLAM S09c flame has been investigated using Large Eddy Simulation and a pre-calculated mixture-fraction-based flamelet model in a preliminary step to simulate this flame. The simplified formulation, despite the fact that it cannot capture localized extinctions, is found to reproduce the experimentally-observed vortex breakdown and the results agree well with available experimental data for velocity and mixture fraction. The unsteady flow features before the burner exit and inside the combustion chamber are analyzed with spectral analysis, correlations, and Proper Orthogonal Decomposition. The results show the presence of longitudinal vortices whose axes rotate with the swirl and which cause separation inside the inlet pipe. With combustion, some of the structures are damped and the spectral peaks shift to higher frequencies. Finally, simulations with the full 3D Conditional Moment Closure equation, which allows spatial and temporal fluctuations of the conditionally-filtered reacting scalars, reproduce successfully the flame lift-off observed experimentally.