Generalized detection of a turbulent front generated by an oscillating grid

This report presents experimental results on the propagation of a turbulent front induced by an oscillating grid starting from rest. The purpose of this preliminary investigation is to implement and validate detection methods of the turbulent/non-turbulent interface, which are based on flow measurements (velocity and vorticity) and scalar intensity, for oscillating grid turbulence. This is done using particle image velocimetry (PIV) and fluorescent dye visualization, separately. The results of both techniques describe the spreading of the turbulent front, confirming the known dependency of the front location, H, on time, t. It is demonstrated, that the level-based detection of a turbulent front can be applied to an unsteady flow, such as grid turbulence advancing into a fluid at rest.     

Numerical simulation of expansion/compression effect on shock induced turbulent flow

The interaction of homogeneous and isotropic turbulence with a shock wave is observed by solving the Reynolds-averaged Navier–Stokes equations with the k–? turbulence model. All turbulent fluctuations are measured at the period of expansion in the turbulent field and during compression by the reflected shock on turbulent field, and it is observed that the longitudinal turbulent velocity fluctuation is enhanced more at the period of expansion due to incident shock wave movement far from the turbulent field. The amplification of the turbulent kinetic energy (TKE) level in the shock/turbulence interaction depends on the shock wave strength and the longitudinal velocity difference across the shock wave. On decreasing the longitudinal velocity difference across the shock, the turbulent kinetic energy (TKE) level is less amplified. The TKE level is amplified by the factor of 1.5–1.8 in the shock/turbulence interaction where the dissipation rate of TKE decreases in all cases of shock/turbulence interaction. After the shock/turbulence interaction, the turbulent dissipative-length scale is amplified slightly and the amplification of the length scales decreases when increasing the shock strength. To cite this article: M.A. Jinnah, K. Takayama, C. R. Mecanique 333 (2005).     

振荡潜流带沉积层?水界面污染物输运的研究

潜流带中污染物质交换与输运特性是影响水资源环境的重要问题之一. 本文对底部为高渗透沉积层的三维槽道振荡流高Schmidt数传质问题进行了大涡模拟研究. 采用动力学亚格子模型来封闭滤波后的三维不可压缩Navier-Stokes方程以及污染物输运方程, 同时采用修正的Darcy-Brinkman-Forcheimer模型来描述沉淀有锌离子污染溶质的可渗透沉积层. 通过对沉积层内外流场和浓度场的统计特性以及瞬态结构的分析, 探究了上覆水体中振荡流驱动作用对污染物输运的动力学影响以及扩散率随振荡周期和振荡角的变化规律. 研究结果表明, 浓度通量中的湍流浓度分量在垂向物质交换中起主导作用, 流向、展向速度, 湍流强度和污染物浓度的波动跟随振荡驱动力呈现准周期变化, 同时发现沉积层?水交界面处的湍流浓度通量与法向湍流强度两者之间的变化具有明确的相关性. 并且在较大振荡角和低频振荡的情况下, 沉积层?水交界面处的有效扩散率增大, 这主要是来自于沉积层?水交界面处流体的猝发行为有效促进了湍流混合和物质交换, 进而增强了污染物的垂向输运.      

稳定分层二层流室内实验早期格栅湍流结构

采用室内实验混合箱和粒子图像测速技术,本文研究了稳定分层无平均剪切二层流(上层淡水、下层盐水)振动湍流结构.对实验录像进行粒子图像测速技术处理,获得垂向二维流场(垂直于格栅平面)瞬时速度和涡量,并用于计算:①时均速度和时均涡量;②均方根速度;③均匀程度和各向同性程度;④平均流强度;⑤时均泰勒的欧拉积分长度尺度;⑥时均湍...     

Turbulent thermal diffusion in a multi-fan turbulence generator with imposed mean temperature gradient

We studied experimentally the effect of turbulent thermal diffusion in a multi-fan turbulence generator which produces a nearly homogeneous and isotropic flow with a small mean velocity. Using particle image velocimetry and image processing techniques, we showed that in a turbulent flow with an imposed mean vertical temperature gradient (stably stratified flow) particles accumulate in the regions with the mean temperature minimum. These experiments detected the effect of turbulent thermal diffusion in a multi-fan turbulence generator for relatively high Reynolds numbers. The experimental results are in compliance with the results of the previous experimental studies of turbulent thermal diffusion in oscillating grid turbulence (Buchholz et al. 2004; Eidelman et al. 2004). We demonstrated that the turbulent thermal diffusion is an universal phenomenon. It occurs independently of the method of turbulence generation, and the qualitative behavior of particle spatial distribution in these very different turbulent flows is similar. Competition between turbulent fluxes caused by turbulent thermal diffusion and turbulent diffusion determines the formation of particle inhomogeneities.     

细颗粒泥沙净冲刷和输移的大涡模拟研究

在传统水沙输移数值模拟研究中一般采用雷诺时均模拟技术(Reynolds-averaged simulation,RANS).与RANS相比,大涡模拟技术(large eddy simulation,LES)能够更加精确反映细部流动结构,计算机的发展使得采用LES探讨水流和泥沙运动规律成为可能.本文尝试给出净冲刷条件下悬沙计算的边界条件,采用动态亚格子模式对循环槽道和长槽道中的水流运动和泥沙输移进行了三维大涡模拟研究.利用直接数值模拟(direct numerical simulation,DNS)结果对LES模型进行了率定,计算结果符合良好,在此基础上初步探讨了泥沙浓度、湍动强度和湍动通量等的分布特征.结果表明,净冲刷条件下输沙平衡时泥沙浓度符合Rouse公式分布,单向流动中泥沙浓度沿着流向逐渐增大.泥沙浓度湍动强度和湍动通量都在近底部达到最大值,沿着垂向迅速减小.湍动黏性系数和扩散系数基本上在水深中间处达到最大.湍动Schmidt数沿着水深方向不是常数,在近底部和自由水面附近较大,在水深中间处较小.     

Modeling droplet dispersion and interphase turbulent kinetic energy transfer using a new dual-timescale Langevin model

Dispersion of spray droplets and the modulation of turbulence in the ambient gas by the dispersing droplets are two coupled phenomena that are closely linked to the evolution of global spray characteristics, such as the spreading rate of the spray and the spray cone angle. Direct numerical simulations (DNS) of turbulent gas flows laden with sub-Kolmogorov size particles, in the absence of gravity, report that dispersion statistics and turbulent kinetic energy (TKE) evolve on different timescales. Furthermore, each timescale behaves differently with Stokes number, a non-dimensional flow parameter (defined in this context as the ratio of the particle response time to the Kolmogorov timescale of turbulence) that characterizes how quickly a particle responds to turbulent fluctuations in the carrier or gas phase. A new dual-timescale Langevin model (DLM) composed of two coupled Langevin equations for the fluctuating velocities, one for each phase, is proposed. This model possesses a unique feature that the implied TKE and velocity autocorrelation in each phase evolve on different timescales. Consequently, this model has the capability of simultaneously predicting the disparate Stokes number trends in the evolution of dispersion statistics, such as velocity autocorrelations, and TKE in each phase. Predictions of dispersion statistics and TKE from the new model show good agreement with published DNS of non-evaporating and evaporating droplet-laden turbulent flow.     

Experiments on the spreading of shear-free turbulence under the influence of confinement and rotation

From Lie-group (symmetry) analysis of the multi-point correlation equation Oberlack and Günther (Fluid Dyn Res 33:453–476, 2003) found three different solutions for the behavior of shear-free turbulence: (i) a diffusion like solution, in which turbulence diffuses freely into the adjacent calm fluid, (ii) a deceleration wave like solution when there is an upper bound for the integral length scale and (iii) a finite domain solution for the case when rotation is applied to the system. This paper deals with the experimental validation of the theory. We use an oscillating grid to generate turbulence in a water tank and Particle Image Velocimetry (PIV) to determine the two-dimensional velocity and out-of-plane vorticity components. The whole setup is placed on a rotating table. After the forcing is initiated, a turbulent layer develops which is separated from the initially irrotational fluid by a sharp interface, the so-called turbulent/non-turbulent interface (TNTI). The turbulent region grows in time through entrainment of surrounding fluid. We measure the propagation of the TNTI and find quantitative agreement with the predicted spreading laws for case one and two. For case three (system rotation), we observe that there is a sharp transition between a 3D turbulent flow close to the source of energy and a more 2D-like wavy flow further away. We measure that the separation depth becomes constant and in this sense, we confirm the theoretical finite domain solution.     

Effect of Particle Clusters on Carrier Flow Turbulence: A Direct Numerical Simulation Study

Experiments indicate that particle clusters that form in fluidized–bed risers can enhance gas-phase velocity fluctuations. Direct numerical simulations (DNS) of turbulent flow past uniform and clustered configurations of fixed particle assemblies at the same solid volume fraction are performed to gain insight into particle clustering effects on gas-phase turbulence, and to guide model development. The DNS approach is based on a discrete-time, direct-forcing immersed boundary method (IBM) that imposes no-slip and no-penetration boundary conditions on each particle’s surface. Results are reported for mean flow Reynolds number Re _p ?=?50 and the ratio of the particle diameter d _p to Kolmogorov scale is 5.5. The DNS confirm experimental observations that the clustered configurations enhance the level of fluid-phase turbulent kinetic energy (TKE) more than the uniform configurations, and this increase is found to arise from a lower dissipation rate in the clustered particle configuration. The simulations also reveal that the particle-fluid interaction results in significantly anisotropic fluid-phase turbulence, the source of which is traced to the anisotropic nature of the interphase TKE transfer and dissipation tensors. This study indicates that when particles are larger than the Kolmogorov scale (d _p ?>?η), modeling the fluid-phase TKE alone may not be adequate to capture the underlying physics in multiphase turbulence because the Reynolds stress is anisotropic. It also shows that multiphase turbulence models should consider the effect of particle clustering in the dissipation model.     

Large eddy simulation of free surface shallow‐water flow

Shallow‐water flow with free surface frequently occurs in ambient water bodies, in which the horizontal scale of motion is generally two orders of magnitude greater than the water depth. To accurately predict this flow phenomenon in more detail, a three‐dimensional numerical model incorporating the method of large eddy simulation (LES) has been developed and assessed. The governing equations are split into three parts in the finite difference solution: advection, dispersion and propagation. The advection part is solved by the QUICKEST scheme. The dispersion part is solved by the central difference method and the propagation part is solved implicitly using the Gauss–Seidel iteration method. The model has been applied to free surface channel flow for which ample experimental data are available for verification. The inflow boundary condition for turbulence is generated by a spectral line processor. The computed results compare favourably with the experimental data and those results obtained by using a periodic boundary condition. The performance of the model is also assessed for the case in which anisotropic grids and filters with horizontal grid size of the order of the water depth are used for computational efficiency. The coarse horizontal grid was found to cause a significant reduction in the large‐scale turbulent motion generated by the bottom turbulence, and the turbulent motion is predominately described by the sub‐grid scale (SGS) terms. The use of the Smagorinsky model for SGS turbulence in this situation is found inappropriate. A parabolic mixing length model, which accounts for the filtered turbulence, is then proposed. The new model can reproduce more accurately the flow quantities. Copyright © 2000 John Wiley & Sons, Ltd.     

Turbulent kinetic energy balance measurements in the wake of a low-pressure turbine blade

The turbulent kinetic energy budget in the wake generated by a high lift, low-pressure two-dimensional blade cascade of the T106 profile was investigated experimentally using hot-wire anemometry. The purpose of this study is to examine the transport mechanism of the turbulent kinetic energy and provide validation data for turbulence modeling. Point measurements were conducted on a high spatial resolution, two-dimensional grid that allowed precise derivative calculations. Positioning of the probe was achieved using a high accuracy traversing mechanism. The turbulent kinetic energy (TKE) convection, production, viscous diffusion and turbulent diffusion were all obtained directly from experimental measurements. Dissipation and pressure diffusion were calculated indirectly using techniques presented and validated by previous investigators. Results for all terms of the turbulent kinetic energy budget are presented and discussed in detail in the present work.     

Stereoscopic multi-planar PIV measurements of in-cylinder tumbling flow

The non-reacting flow field within the combustion chamber of a motored direct-injection spark-ignition engine with tumble intake port is measured. The three-dimensionality of the flow necessitates the measurement of all three velocity components via stereoscopic particle-image velocimetry in multiple planes. Phase-locked stereoscopic PIV is applied at 15 crank angles during the intake and compression strokes, showing the temporal evolution of the flow field. The flow fields are obtained within a set of 14 axial planes, covering nearly the complete cylinder volume. The stereoscopic PIV setup applied to engine in-cylinder flow and the arising problems and solutions are discussed in detail. The three-dimensional flow field is reconstructed and analyzed using vortex criteria. The tumble vortex is the dominant flow structure, and this vortex varies significantly regarding shape, strength, and position throughout the two strokes. The tumble vortex center moves clockwise through the combustion chamber. At first, the tumble has a c-shape which turns into an almost straight tube at the end of the compression. Small-scale structures are analyzed by the distribution of the turbulent kinetic energy. It is evident that the symmetry plane only represents the 3D flow field after 100 CAD. For earlier crank angles, both kinetic energy (KE) and turbulent kinetic energy (TKE) in the combustion chamber are well below the KE and TKE in the symmetry plane. This should be taken into account when the injection and breakup of the three-dimensional fuel jet are studied. The mean kinetic energy is conserved until late compression by the tumble motion. This conservation ensures through the excited air motion an enhancement of the initial air-fuel mixture which is of interest for direct-injection gasoline engines.     

Effect of bubbles on turbulent kinetic energy transport in downward flow measured by time-resolved PTV

A time-resolved particle tracking velocimetry (PTV) system and a shape projection imaging system were used to investigate the turbulence modifications by bubbles in a downward bubbly flow. Two bubble sizes and three mean void fractions were tested at a Reynolds number of about 20,000. The strong modifications in the mean velocity, turbulent kinetic energy (TKE) budget, and velocity spectra are observed in the central region of the pipe that has a high local void fraction. In particular, kinetic energy production decreased, whereas the TKE dissipation rate increased. This suggests that the transfer of energy due to bubbles has a very large effect on the TKE budget. Moreover, velocity spectra reveal that the presence of bubbles modifies the length scales of turbulent eddies, which contain, transfer, and dissipate energy.     

Turbulence manipulation in air–water flows on a stepped chute: An experimental study

In a stepped channel operating with large flow rates, the flow skims over the pseudo-bottom formed by the step edges as a coherent stream. Intense three-dimensional recirculation is maintained by shear stress transmission from the mainstream to the step cavities, while significant free-surface aeration takes place. The interactions between free-surface aeration and cavity recirculation are investigated herein with seven step cavity configurations. The experiments were conducted in a large stepped channel operating at large Reynolds numbers. For some experiments, triangular vanes, or longitudinal ribs, were placed across the step cavities to manipulate the flow turbulence to enhance the interactions between the mainstream flow and the cavity recirculation region. The results showed a strong influence of the vanes on the air–water flow properties in both free-stream and cavity flows. The findings demonstrate some passive turbulence manipulation in highly turbulent air–water flows.     

定床弯道内水沙两相运动的数值模拟

在适体同位网格中采用非正交曲线坐标系下的三维k-ε-kp固液两相双流体湍流模型研究弯道内水流和悬浮泥沙运动,主要计算了试验室S型水槽内清水流动的三维流场、120°弯道内水沙两相流动中底沙与底流的运动轨迹以及S型水槽内水沙两相流动的两相流场和泥沙浓度场. 对于S型水槽内清水流动,数值结果与试验结果吻合良好. 120°弯道内水沙两相流动中固液两相的运动轨迹在弯道直线段基本重合,在弯道内泥沙轨迹逐步偏离水体轨迹,其偏离程度随泥沙粒径增大而增大. 从S型水槽内水沙两相流动计算结果中发现泥沙纵向流速在壁面附近比水流纵向速度大,在远离壁面区域比水流纵向速度小;弯道内泥沙横向流速比水流横向流速小;垂向流速在直线段和泥沙沉速相当,在弯道内受螺旋水流影响而变化;两相流速差别随泥沙粒径增大而变大;泥沙浓度呈现下浓上稀的分布,在弯道内横向断面上呈现凸岸大凹岸小的分布,泥沙浓度随泥沙粒径增大而减小.      

Laminar/turbulent oscillating flow in circular pipes

A two-dimensional oscillating flow analysis was conducted simulating the gas flow inside Stirling engine heat exchangers. Both laminar and turbulent oscillating pipe flow were investigated numerically for Re_max = 1,920 (Va = 80), 10,800 (Va = 272), 19,300 (Va = 272), and 60,800 (Va = 126). The results are here compared with experimental results of previous investigators. Predictions of the flow regime on present oscillating flow conditions are also checked by comparing velocity amplitudes and phase difference with those from laminar theory and quasi-steady profile. A high Reynolds number k-ε turbulence model was used for turbulent oscillating pipe flow. Finally, the performance of the k-ε model was evaluated to explore the applicability of quasi-steady turbulent models to unsteady oscillating flow analysis.     

Resuspension and Sedimentation of Particles from a Sediment Bed by Turbulent Jets

An experimental study was carried out to investigate the resuspension of heavy particles of a sediment bed by axi-symmetric turbulent jets issuing from below. The case of point jets and loosely-held monodispersed particles was considered, in particular, the cases of a single jet and an array of jets located in the intersection points of a square grid. The aim was to map the flow field and the nature of sediment distribution as a function of governing parameters and to parameterize salient observables such as the maximum height of rise of sediments. The results were extrapolated to study sediment resuspension in karstic lakes, especially lake Banyoles located in Catalonia, northeastern Spain.     

Turbulence measurements in the bubbly flow region of hydraulic jumps

A hydraulic jump is characterized by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air–water flow measurements were performed in a large-size facility using two types of phase-detection intrusive probes: i.e. single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent integral time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of Chanson [H. Chanson, Bubbly flow structure in hydraulic jump, Eur. J. Mech. B/Fluids 26 (3) (2007) 367–384], providing the turbulent scales of the eddy structures advecting the air bubbles in the developing shear layer. The length scale L_xz is an integral air–water turbulence length scale which characterized the transverse size of the large vortical structures advecting the air bubbles. The experimental data showed that the dimensionless integral turbulent length scale L_xz/d₁ was closely related to the inflow depth: i.e. L_xz/d₁ = 0.2–0.8, with L_xz increasing towards the free-surface.     

Experimental investigation of streamwise velocity fluctuation based on the Reynolds-number dependency in turbulent viscoelastic-fluid flow

This paper describes an experimental verification of energy supply mechanisms for the streamwise component of the turbulent kinetic energy (TKE) at different Reynolds numbers in viscoelastic-fluid flow. We investigated the characteristics of the streamwise turbulent velocity fluctuation by analyzing the production and turbulent diffusion terms in the TKE transport equation. In addition, we reported on the Reynolds-number dependency in a high Reynolds-number regime where direct numerical simulation cannot demonstrate changes in fluid properties. Based on the experimental verification, we proposed a conceptual model of the energy-exchange term between the TKE and the elastic energy, with focusing on the dependency of the fluid properties on the shear stress. This model is indirectly reflected in the streamwise TKE, the instantaneous velocity field, and the wave number relevant to energy-containing eddies. The main gain term of the TKE switches between the energy-exchange term and the production term dependently on the Reynolds number: as the Reynolds number exceeds the value which provides the maximum drag reduction rate, the production term becomes dominant and the magnitude of streamwise TKE becomes high compared to the water flow case.     

Total pressure fluctuations and two-phase flow turbulence in hydraulic jumps

The large-scale turbulence and high air content in a hydraulic jump restrict the application of many traditional flow measurement techniques. This paper presents a physical modelling of hydraulic jump, where the total pressure and air–water flow properties were measured simultaneously with intrusive probes, namely a miniature pressure transducer and a dual-tip phase-detection probe, in the jump roller. The total pressure data were compared to theoretical values calculated based upon void fraction, water depth and flow velocity measured by the phase-detection probe. The successful comparison showed valid pressure measurement results in the turbulent shear region with constant flow direction. The roller region was characterised by hydrostatic pressure distributions, taking into account the void fraction distributions. The total pressure fluctuations were related to both velocity fluctuations in the air–water flow and free-surface dynamics above the roller, though the time scales of these motions differed substantially.