圆管中湍旋流的精细模拟与总流机械能损失

湍旋流广泛存在于自然界与工程实际中,其流动特性及机械能损失的研究具有很高的理论意义与工程价值.通过开发数值旋流生成装置,采用雷诺应力模式封闭雷诺时均方程及有针对性的网格布置,实现了对圆管中湍旋流的三维精细模拟.计算成果表明:1)湍旋流的旋流数沿程呈分段指数衰减,其衰减系数有随雷诺数的增加而减小的趋势;2)湍旋流的切向时均流速沿径向的变化在旋流数较大时可近似用组合涡的速度分布来描述;3)与一般恒定均匀圆管湍流的总流机械能损失系数相比,带有旋流的圆管湍流的总流机械能损失系数要普遍大于前者,且旋流数越大,总流机械能损失系数越大.

Stabilization of Pattern in Complex Ginzburg-Landau Equation with Spatial Perturbation Scheme

In this paper, a new scheme of spatial perturbation is proposed to stabilize the pattern in the oscillatory media, which could be described with the complex Ginzburg-Landau equation. The numerical simulation results confirm that the spiral wave, antispiral wave and spatiotemporal chaos in the complex Ginzburg-Landau equation could be suppressed, and the scheme is also discussed with the conservative field theory. Furthermore, the spatiotemporal noise is also introduced into the whole media, it just confirms that it is robust to the spatiotemporal noise.

Effects of the Blade Shape on the Trailing Vortices in Liquid Flow Generated by Disc Turbines

Particle image velocimetry technique was used to analyze the trailing vortices and elucidate their rela-tionship with turbulence properties in a stirred tank of 0.48 m diameter,agitated by four different disc turbines,in-cluding Rushton turbine,concaved blade disk turbine,half elliptical blade disk turbine,and parabolic blade disk turbine.Phase-averaged and phase-resolved flow fields near the impeller blades were measured and the structure of trailing vortices was studied in detail.The location,size and strength of vortices were determined by the simplified λ2-criterion and the results showed that the blade shape had great effect on the trailing vortex characteristics.The larger curvature resulted in longer residence time of the vortex at the impeller tip,bigger distance between the upper and lower vortices and longer vortex life,also leads to smaller and stronger vortices.In addition,the turbulent ki-netic energy and turbulent energy dissipation in the discharge flow were determined and discussed.High turbulent kinetic energy and turbulent energy dissipation regions were located between the upper and lower vortices and moved along with them.Although restricted to single phase flow,the presented results are essential for reliable de-sign and scale-up of stirred tank with disc turbines.

Large eddy simulation of the rotation effect on the ocean turbulence kinetic energy budget in the surface mixed layer

A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.

Turbulence Characteristics of Swirling Reacting Flow in a Combustor with Staged Air Injection

This paper presents an experimental investigation of the turbulent reacting flow in a swirl combustor with staged air injection. The air injected into the combustor is composed of the primary swirling jet and the secon-dary non-swirling jet. A three dimension-laser particle dynamic analyzer (PDA) was employed to measure the in-stantaneous gas velocity. The probability density functions (PDF) for the instantaneous gas axial and tangential ve-locities at each measuring location, as well as the radial profiles of the root mean square of fluctuating gas axial and tangential velocities and the second-order moment for the fluctuating gas axial and tangential velocities are ob-tained. The measured results delineate the turbulence properties of the swirling reacting flow under the conditions of staged combustion.

Effects of Coordinate Rotation on Turbulent Flux Measurements during Wintertime Haze Pollution in Beijing, China

Eddy-covariance observations from the Beijing 325-m meteorological tower are used to evaluate the effects of coordinate rotation on the turbulent exchange of momentum and scalars during wintertime haze pollution(January-February 2013). Two techniques are used in the present evaluation; namely, the natural wind coordinate(NWC) and the planar fit coordinate(PFC), with the latter being applied by means of two methods for linear regression(i.e., overall and sector-wise). The different techniques show a general agreement in both turbulent fluxes and transport efficiencies, especially evident at the lower, 140-m level above the ground(compared to the higher, 280-m level), perhaps implying that the selection of a technique for coordinate rotation(NWC or PFC) is less of a concern for a sufficiently low level, despite the complexities of urban terrain. Additionally, sector-wise regression is a recommended approach for practical application of the PFC in a complex urban environment subjected to particulate pollution, because this method is found to produce a better correlation between the mean vertical velocity at the 140- and 280-m heights.

Reconstruction of Turbulent Swirling Flow in a Dump Combustor

Experimental data of the continuous evolution of fluid flow characteristics in a dump combustor is very useful and essential for better and optimum designs of gas turbine combustors and ramjet engines. Unfortunately, experimental techniques such as 2D and/or 3D LDV （Laser Doppler Velocimetry） measurements provide only limited discrete information at given points; especially, for the cases of complex flows such as dump combustor swirling flows. For this type of flows, usual numerical interpolating schemes appear to be unsuitable. Recently, neural networks have emerged as viable means of expanding a finite data set of experimental measurements to enhance better understanding of a particular complex phenomenon. This study showed that generalized feed forward network is suitable for the prediction of turbulent swirling flow characteristics in a model dump combustor. These techniques are proposed for optimum designs of dump combustors and ramjet engines.

泄洪洞内具有阻塞与扩散段的水平旋流消能方式水力特性的试验与数值模拟(2)

本文采用Realizable k-ε双方程紊流模型和VOF法,并结合模型试验,对水平旋流复合式内消能泄洪洞的流速分布、旋流角、湍动能、耗散率、以及消能率进行了分析和计算。结果表明:起旋室段存在径向速度,切向入流条件较差;轴向流速沿径向呈对称分布,最大值靠近洞壁;切向流速沿径向呈反对称分布,最大值靠近气水交界面,可视为组合涡分布;旋流角在起旋室段与阻塞扩散段变化较大,中间旋流洞段变化较平缓;紊动能、耗散率在起旋室出口附近及阻塞扩散段较大;计算总消能率约为83%,其中竖井段、旋流洞段和阻塞扩散段各占4.58%、28.45%和49.51%。旋流阻塞的设置不仅减小了旋流洞内发生空化空蚀的可能性,而且增大了洞内消能率,提高了适用水头。

Scaling properties of Navier-Stokes turbulence

The property of the velocity field and the cascade process of the fluid flow are key problems in turbulence research. This study presents the scaling property of the turbulent velocity field and a mathematical description of the cascade process, using the following methods: (1) a discussion of the general self-similarity and scaling invariance of fluid flow from the viewpoint of the physical mechanism of turbulent flow; (2) the development of the relationship between the scaling indices and the key parameters of the She and Leveque (SL) model in the inertial range; (3) an investigation of the basis of the fractal model and the multi-fractal model of turbulence; (4) a demonstration of the physical meaning of the flowing field scaling that is related to the real flowing vortex. The results illustrate that the SL model could be regarded as an approximate mathematical solution of Navier-Stokes (N-S) equations, and that the phenomena of normal scaling and anomalous scaling is the result of the mutual interactions among the physical factors of nonlinearity, dissipation, and dispersion. Finally, a simple turbulent movement conceptional description model is developed to show the local properties and the instantaneous properties of turbulence.

Large eddy simulation of unconfined turbulent swirling flow

Large eddy simulations(LES) were performed to study the non-reacting flow fields of a Cambridge swirl burner. The dynamic Smagorinsky eddy viscosity model is used as the sub-grid scale turbulence model. Comparisons of experimental data show that the LES results are capable of predicting mean and root-mean-square velocity profiles. The LES results show that the annular swirling flow has a minor impact on the formation of the bluff-body recirculation zone. The vortex structures near the shear layers, visualized by the iso-surface of Q-criterion, display ring structures in non-swirling flow and helical structures in swirling flow near the burner exit. Spectral analysis was employed to predict the occurrence of flow oscillations induced by vortex shedding and precessing vortex core(PVC). In order to extract accurately the unsteady large-scale structures in swirling flow, a three-dimensional proper orthogonal decomposition(POD) method was developed to reconstruct turbulent fluctuating velocity fields. POD analysis reveals that flow fields contain co-existing helical and toroidal shaped coherent structures. The helical structure associated with the PVC is the most energetic dynamic flow structure. The latter toroidal structure associated with vortex shedding has lower energy content which indicates that it is a secondary structure.

Experimental Study of the Flow in Helical Circular Pipes:Torsion Effect on the Flow Velocity and Turbulence

An objective of the present paper is to experimentally clarify the torsion effect on the flow in helical circular pipes. We have made six helical circular pipes having different pitches and common non-dimensional curvature δ of about 0.1. The torsion parameter β0, which is defined by β0 = τ/（2δ）1/2 with non-dimensional torsion r, are taken to be 0.02, 0.45, 0.69, 1.01, 1.38 and 1.89 covering from small to very large pitch. The velocity distributions and the turbulence of the flow are measured using an X-type hot-wire anemometer in the range of the Reynolds number from 200 to 20000. The results obtained are summarized as follows： The mean secondary flow pattern in a cross section of the pipe changes from an ordinary twin-vortex type as is seen in a curved pipe without torsion （toroidal pipe） to a single vortex type after one of the twin-vortex gradually disappears as β0 increases. The circulation direction of the single vortex is the same as the direction of torsion of the pipe. The mean velocity distribution of the axial flow is similar to that of the toroidal pipe at small β0, but changes its shape as β0 increases, and attains the shape similar to that in a straight circular pipe when ,β0 = 1.89. It is also found that the critical Reynolds number, at which the flow shows a marginal behavior to turbulence, decreases as ,β0 increases for small ,β0, and then increases after taking a minimum at ,β0 ≈ 1.4 as ,β0 increases. The minimum of the critical Reynolds number experimentally obtained is about 400 at ,β0 ≈ 1.4.

Evolution of the vortex structures and turbulent spots at the late-stage of transitional boundary layers

The nonlinear evolution process of new vortex structures at the late-stage of the transition, including the 3-D spatial structure of barrel-shaped vortex and "dark spots" structure observed by experiment research, has been confirmed by our computational results. The formation mechanisms of these structures have been explored. It is revealed that the new vortex structures, the ring-like vortex chain and induced disturbance velocities play a dominant role in the generation of turbulent spots.

Drive Control of Spiral Wave and Turbulence by a Target Wave in CGLE

Suppression of spiral wave and turbulence in the complex Cinzburg-Landau equation （CCLE） plays a prominent role in nonlinear science and complex dynamical system. In this paper, the nonlinear behavior of the proposed drive-response system, which consists of two coupled OGLEs, is investigated and controlled by a state error feedback controller with the lattice Boltzmann method. First, spiral wave appropriate parameters of the response system under the no-flux and turbulence are, respectively, generated by selecting boundary and perpendicular gradient initial conditions. Then, based on the random initial condition, the target wave yielded by introducing spatially localized inhomogeneity into the drive system is applied on the above response system. The numerical simulation results show that the spiral wave and turbulence existing in the response system could be successfully eliminated by the target wave in the drive system during a short evolution time. Furthermore, it turns out that the transient time for the drive course is related to the control intensity imposed on the whole media.

Comparison of Reynolds Averaged Navier-Stokes Based Simulation and Large-eddy Simulation for One Isothermal Swirling Flow

The flow structure of one isothermal swirling case in the Sydney swirl flame database was studied using two numerical methods. Results from the Reynolds-averaged Navier-Stokes (RANS) approach and large eddy simulation (LES) were compared with experimental measurements. The simulations were applied in two different Cartesian grids which were investigated by a grid independence study for RANS and a post-estimator for LES. The RNG k-ε turbulence model was used in RANS and dynamic Smagorinsky-Lilly model was used as the sub-grid scale model in LES. A validation study and cross comparison of ensemble average and root mean square (RMS) results showed LES outperforms RANS statistic results. Flow field results indicated that both approaches could capture dominant flow structures, like vortex breakdown (VB), and precessing vortex core (PVC). Streamlines indicate that the formation mechanisms of VB deducted from the two methods were different. The vorticity field was also studied using a velocity gradient based method. This research gained in-depth understanding of isothermal swirling flow.

Numerical Investigation of Heat Transfer Coefficient in a Low Speed 1.5-Stage Turbine

The paper numerically investigated the heat transfer coefficients over the rotating blades in a 1.5-stage turbine. The hexahedral structured grids and k-ε turbulence model were applied in the simulation. A film hole with diameter of 0.004 m, angled 36°and 28° tangentially to the suction side and pressure side in streamwise respectively, was set in the middle span of the rotor blade. Simulations are done at three different rotating numbers of 0.0239, 0.0265 and 0.0280 with the blowing ratio varying from 0.5 to 2.0. The effects of mainstream Reynolds number and density ratio are also compared. Results show that increasing blowing ratio can increase the heat transfer coefficient ratio on the pressure side, but the rule is parabola on the suction side. Besides, increasing rotating number and Reynolds number is positive while increasing density ratio is negative to the heat transfer on both the pressure side and the suction side.

Numerical Simulation and Experiment Analyses for the Gas-liquid Two-phase Vortex Pump

Based on the gas-liquid two-phase mixture transportation test, the k-c-A; turbulence model was applied to simulate the two-phase turbulent flow in a vortex pump. By comparing the simulation and experiment results, inner flow features were revealed. The bubbles in the channel distribute mainly at the pressure side of the blades, and the aggregation degree of the bubbles is enhanced with an increase in inlet gas volume fraction. Experimental results indicate that the influence of the gas phase on vortex pump performance is small when the gas volume fraction is less than 10%. When the gas volume fraction contiuuously increases to 15%, the characteristic curves abruptly drop due to the gas blocking phenomenon.

Effect of rib height on turbulent flow and heat transfer of kerosene in rectangular duct

The coupling effects of rib heights and fluid properties on turbulent convective heat transfer of kerosene flow through the rectangular duct on the ribbed bottom wall are studied numerically in this paper.The numerical simulation is based on the ten components surrogate model of kerosene and the Reynolds average method combined with the re-normalized group(RNG)k-εturbulence model.The turbulent vortex structures and heat transfer characteristics of kerosene flowing over rectangular ribs of different heights are obtained.The results show that three dimensional vortices are generated by the ribs,and the vortices alter local flow significantly,leading to both enhanced and reduced convective heat transfer at different locations near the ribs.In addition,it is found that with the increase of rib height,the average Nusselt number and the wall friction factor on the ribbed wall also increase.For the present study,the maximum heat transfer enhancement rate of kerosene flow is 72.16%,and the ratio of rib-to-duct height is 0.75.