Turbulent Flow Produced by Piston Motion in a Spark-ignition Engine

Turbulence produced by the piston motion in spark-ignition engines is studied by 2D axisymmetric numerical simulations in the cylindrical geometry as in the theoretical and experimental work by Breuer et al. (Flow Turbul Combust 74:145, 2005). The simulations are based on the Navier–Stokes gas-dynamic equations including viscosity, thermal conduction and non-slip at the walls. Piston motion is taken into account as a boundary condition. The turbulent flow is investigated for a wide range of the engine speed, 1,000–4,000 rpm, assuming both zero and non-zero initial turbulence. The turbulent rms-velocity and the integral length scale are investigated in axial and radial directions. The rms-turbulent velocity is typically an order-of-magnitude smaller than the piston speed. In the case of zero initial turbulence, the flow at the top-dead-center may be described as a combination of two large-scale vortex rings of a size determined by the engine geometry. When initial turbulence is strong, then the integral turbulent length demonstrates self-similar properties in a large range of crank angles. The results obtained agree with the experimental observations of Breuer et al. (Flow Turbul Combust 74:145, 2005).     

Experimental study of turbulent thermal diffusion in oscillating grids turbulence

We have experimentally detected a new effect, turbulent thermal diffusion, as predicted theoretically by Elperin et al. (Phys Rev Lett (1996) 76:224–228) and associated with the turbulent transport of inertial particles. The essence of this effect is an appearance of a non-diffusive mean flux of particles in the direction of the mean heat flux. This results in formation of large-scale inhomogeneities in the spatial distribution of inertial particles that are accumulated in regions of minimum mean temperature in the surrounding fluid. The experiments were performed in oscillating grids turbulence with an imposed mean temperature gradient. We used Particle Image Velocimetry to determine the turbulent velocity field, and an Image Processing Technique based on the analysis of the intensity of the Mie scattering to determine the spatial distribution of tracer particles. Analysis of the intensity of laser light Mie scattering by tracer particles showed that the tracer particles accumulate in the vicinity of the minimum of the mean temperature. The latter finding confirms the existence of the effect of turbulent thermal diffusion.     

Full 3D correlation tensor computed from double field stereoscopic PIV in a high Reynolds number turbulent boundary layer

The turbulence structure near a wall is a very active subject of research and a key to the understanding and modeling of this flow. Many researchers have worked on this subject since the fifties Hama et al. (J Appl Phys 28:388–394, 1957). One way to study this organization consists of computing the spatial two-point correlations. Stanislas et al. (C R Acad Sci Paris 327(2b):55–61, 1999) and Kahler (Exp Fluids 36:114–130, 2004) showed that double spatial correlations can be computed from stereoscopic particle image velocimetry (SPIV) fields and can lead to a better understanding of the turbulent flow organization. The limitation is that the correlation is only computed in the PIV plane. The idea of the present paper is to propose a new method based on a specific stereoscopic PIV experiment that allows the computation of the full 3D spatial correlation tensor. The results obtained are validated by comparison with 2D computation from SPIV. They are in very good agreement with the results of Ganapthisubramani et al. (J Fluid Mech 524:57–80, 2005a).     

Development and characterization of a variable turbulence generation system

Experimental turbulent combustion studies require systems that can simulate the turbulence intensities [u′/U ₀ ~ 20–30% (Koutmos and McGuirk in Exp Fluids 7(5):344–354, 1989)] and operating conditions of real systems. Furthermore, it is important to have systems where turbulence intensity can be varied independently of mean flow velocity, as quantities such as turbulent flame speed and turbulent flame brush thickness exhibit complex and not yet fully understood dependencies upon both U ₀ and u′. Finally, high pressure operation in a highly pre-heated environment requires systems that can be sealed, withstand high gas temperatures, and have remotely variable turbulence intensity that does not require system shut down and disassembly. This paper describes the development and characterization of a variable turbulence generation system for turbulent combustion studies. The system is capable of a wide range of turbulence intensities (10–30%) and turbulent Reynolds numbers (140–2,200) over a range of flow velocities. An important aspect of this system is the ability to vary the turbulence intensity remotely, without changing the mean flow velocity. This system is similar to the turbulence generators described by Videto and Santavicca (Combust Sci Technol 76(1):159–164, 1991) and Coppola and Gomez (Exp Therm Fluid Sci 33(7):1037–1048, 2009), where variable blockage ratio slots are located upstream of a contoured nozzle. Vortical structures from the slots impinge on the walls of the contoured nozzle to produce fine-scale turbulence. The flow field was characterized for two nozzle diameters using three-component Laser Doppler velocimetry (LDV) and hotwire anemometry for mean flow velocities from 4 to 50 m/s. This paper describes the key design features of the system, as well as the variation of mean and RMS velocity, integral length scales, and spectra with nozzle diameter, flow velocity, and turbulence generator blockage ratio.     

A simple model for the effect of peak-locking on the accuracy of boundary layer turbulence statistics in digital PIV

A simple model was constructed to study the effect of peak-locking on the accuracy of particle image velocimetry (PIV) turbulence statistics. A crucial parameter is the ratio between the root-mean-square (rms) velocity and the discretization velocity, which reflects the number of peaks distributed over the velocity probability density functions. When the ratio of the discretization velocity, which is set by the PIV setup parameters, to the rms, given by the flow, is larger than two, the maximum errors introduced in the mean and rms values become significant (larger than 1%). The errors introduced also depend on the amplitude, or severity, of the peak-locking, and whether the mean displacement corresponds to an integer or a fractional number of pixels. The peak-locking affects the statistical moments of different order in such a way that the errors are phase shifted. The proposed model can be used to predict errors in the turbulence statistics in a laboratory PIV experiment. According to our model predictions, the most significant influence of peak-locking in a boundary layer type of flow is an overall underestimation of the wall-normal rms. Our predictions are in good agreement with our experimental results from turbulent boundary layers and the recent experimental results from a turbulent channel flow by Christensen (Exp Fluids 36:484–497, 2004) for a case of moderate peak-locking.

---

Hysteresis phenomenon in turbulent convection

Coherent large-scale circulations of turbulent thermal convection in air have been studied experimentally in a rectangular box heated from below and cooled from above using Particle Image Velocimetry. The hysteresis phenomenon in turbulent convection was found by varying the temperature difference between the bottom and the top walls of the chamber (the Rayleigh number was changed within the range of 10⁷–10⁸). The hysteresis loop comprises the one-cell and two-cells flow patterns while the aspect ratio is kept constant (A=2–2.23). We found that the change of the sign of the degree of the anisotropy of turbulence was accompanied by the change of the flow pattern. The developed theory of coherent structures in turbulent convection (Phys Rev E 66:1–15, 2002, Boundary-Layer Meteorol, 2005) is in agreement with the experimental observations. The observed coherent structures are superimposed on a small-scale turbulent convection. The redistribution of the turbulent heat flux plays a crucial role in the formation of coherent large-scale circulations in turbulent convection.     

A stochastic Langevin model of turbulent particle dispersion in the presence of thermophoresis

Direct numerical simulation (DNS) and experimental data have shown that inertial particles exhibit concentration peaks in isothermal turbulent boundary layers, whereas tracer-like particles remain well mixed in the domain. It is therefore expected that the interactions between turbulence and thermophoresis will be strong in particle-laden flows where walls and carrier fluid are at significantly different temperatures. To capture turbulent particle dispersion with active thermophoresis, a coupled CFD-Lagrangian continuous random walk (CRW) model is developed. The model uses 3D mean flow velocities obtained from the Fluent 6.3 CFD code, to which are added turbulent fluid velocities derived from the normalized Langevin equation which accounts for turbulence inhomogeneities. The mean thermophoretic force is included as a body force on the particle following the Talbot formulation. Validation of the model is performed against recent integral thermophoretic deposition data in long pipes as well as the TUBA TT28 test with its detailed local deposition measurements. In all cases, the agreement with the data is very good. In separate parametric studies in a hypothetical cooled channel flow, it is found that turbulence strongly enhances thermophoretic deposition of particles with dimensionless relaxation times τ⁺ of order 1 or more. On the other hand, the thermophoretic deposition of very small inertia particles (τ⁺ < 0.2) in the asymptotic region far from the injection point tends to that which characterizes stagnant flow conditions, in agreement with the DNS results of Thakurta et al.     

Turbulent film condensation on a horizontal elliptical tube

An investigation has been made of turbulent film condensation on a horizontal elliptical tube. The present study is based on Colburn analogy [1] and potential flow theory to determine the high tangential velocity of vapor flow at the boundary layer and to define the local interfacial shear owing to high velocity vapor flow across the tube surface. The condensate film flow and local/or mean heat transfer characteristics from a horizontal elliptical tube with variable ellipticities, e, under the influence of Froude number, sub-cooling parameter and system pressure have been performed. The present result for dimensionless mean heat transfer coefficient reduces to the same result obtained by Sarma et al.s [2] e=0 (circular tube). Compared with laminar model by Yang and Hsu [3], the present turbulent model shows in better agreement with Michaels experimental data [4] (for e=0). The dependence of mean Nusselt coefficient on the effect of n (power of Reynolds) [1] is also discussed.     

Roughness effects on turbulent plane wall jets in an open channel

This paper reports the effects of surface roughness on the mean flow characteristics for a turbulent plane wall jet created in an open channel. The velocity measurements were obtained using a laser Doppler anemometer over smooth and transitionally rough surfaces. The power law proposed by George et al. (2000) was used to determine the friction velocity. Both conventional scaling and the momentum–viscosity scaling proposed by Narasimha et al. (1973) were used to analyze the streamwise evolution of the flow. The results show that surface roughness increases the skin friction coefficient and the inner layer thickness, but the jet half-width is nearly independent of surface roughness.     

分层流体中栅格湍流的特性

本文考察了分层流体中栅格湍流衰减和演化过程的细节,通过对空间场的信息二维图象处理,获得了湍流动能、耗散率、功率谱及多种湍流尺度等主要湍流特征量。结果显示层结加速了湍流垂直动能的衰减,而水平动能的衰减几乎很少受影响。流场各特征量与前人结果也符合得较好。     

Flow temporal reconstruction from non time-resolved data part II: practical implementation,methodology validation,and applications

This paper proposes a method to sort experimental snapshots of a periodic flow using information from the first three POD coefficients. Even in presence of turbulence, phase-average flow fields are reconstructed with this novel technique. The main objective is to identify and track traveling coherent structures in these pseudo periodic flows. This provides a tool for shedding light on flow dynamics and allows for dynamical contents comparison, instead of using mean statistics or traditional point-based correlation techniques. To evaluate the performance of the technique, apart from a laminar test on the relative strength of the POD modes, four additional tests have been performed. In the first of these tests, time-resolved PIV measurements of a turbulent flow with an externally forced main frequency allows to compare real phase-locked average data with reconstructed phase obtained using the technique proposed in the paper. The reconstruction technique is then applied to a set of non-forced, non time-resolved Stereo PIV measurements in an atmospheric burner, under combustion conditions. Besides checking that the reconstruction on different planes matches, there is no indication of the magnitude of the error for the proposed technique. In order to obtain some data regarding this aspect, two additional tests are performed on simulated non-externally forced laminar flows with the addition of a digital filter resembling turbulence (Klein et al. in J Comput Phys 186:652–665, 2003). With this information, the limitation of the technique applicability to periodic flows including turbulence or secondary frequency features is further discussed on the basis of the relative strength of the Proper Orthogonal Decomposition (POD) modes. The discussion offered indicates coherence between the reconstructed results and those obtained in the simulations. In addition, it allows defining a threshold parameter that indicates when the proposed technique is suitable or not. For those researchers interested on the background and possible generalizations of the technique, part I of this work (Legrand et al. in Exp Fluid (submitted in 2010) 2011) offers the mathematic fundamentals of the general space–time reconstruction technique using POD coefficients. Noteworthy, the involved computational time is relatively small: all the reconstructions have been performed in the order of minutes.     

PIV study of the influence of large-scale streamwise vortices on a turbulent boundary layer

Vortex generator jets were used to generate large-scale longitudinal vortices embedded in a flat-plate turbulent boundary layer. The investigation was performed in a water tunnel, measuring instantaneous flow fields in planes parallel and normal to the flat plate, using particle image velocimetry. The objective of the research was to observe the response of near-wall turbulence to the imposed perturbing flow. It was shown that a small-amplitude forcing vortical flow had significant influence on the mean and fluctuating velocity profiles. Moreover, particle image velocimetry permitted speculation upon the behaviour of the wall velocity streaks under the action of the perturbing forcing vortical flow.     

Experiments on the Flow Field and Acoustic Properties of a Mach number 0·75 Turbulent Air Jet at a Low Reynolds Number

In this paper we present the experimental results of a detailed investigation of the flow and acoustic properties of a turbulent jet with Mach number 0·75 and Reynolds number 3·5 10³. We describe the methods and experimental procedures followed during the measurements, and subsequently present the flow field and acoustic field. The experiment presented here is designed to provide accurate and reliable data for validation of Direct Numerical Simulations of the same flow. Mean Mach number surveys provide detailed information on the centreline mean Mach number distribution, radial development of the mean Mach number and the evolution of the jet mixing layer thickness both downstream and in the early stages of jet development. Exit conditions are documented by measuring the mean Mach number profile immediately above the nozzle exit. The fluctuating flow field is characterised by means of a hot-wire, which produced radial profiles of axial turbulence at several stations along the jet axis and the development of flow fluctuations through the jet mixing layer. The axial growth rate of the jet instabilities are determined as function of Strouhal number, and the axial development of several spectral components is documented. The directivity of the overall sound pressure level and several spectral components were investigated. The spectral content of the acoustic far field is shown to be compatible with findings of hot-wire experiments in the mixing layer of the jet. In addition, the measured acoustic spectra agree with Tam’s large-scale similarity and fine-scale similarity spectra (Tam et al., AIAA Pap 96, 1996).     

On improvement of PIV image interrogation near stationary interfaces

In this paper the problem posed by interfaces when present in PIV measurements is addressed. Different image pre-processing, processing and post-processing methodologies with the intention to minimize the interface effects are discussed and assessed using Monte Carlo simulations. Image treatment prior to the correlation process is shown to be incapable of fully removing the effects of the intensity pedestal across the object edge. The inherent assumption of periodicity in the signal causes the FFT-based correlation technique to perform the worst when the correlation window contains a signal truncation. Instead, an extended version of the masking technique introduced by Ronneberger et al. (Proceedings of the 9th international symposium on applications of laser techniques to fluid mechanics, Lisbon, 1998) is able to minimize the interface-correlation, resolving only the particle displacement peak. Once the displacement vector is obtained, the geometric center of the interrogation area is not the correct placement. Instead, the centre of mass position allows an unbiased representation of the wall flow (Usera et al. in Proceedings of the 12th international symposium on applications of laser techniques to fluid mechanics, Lisbon, 2004). The aforementioned concepts have been implemented in an adaptive interrogation methodology (Theunissen et al. in Meas Sci Technol 18:275–287, 2007) where additionally non-isotropic resolution and re-orientation of the correlation windows is applied near the interface, maximizing the wall-normal spatial resolution. The increase in resolution and robustness are demonstrated by application to a set of experimental images of a flat-plate, subsonic, turbulent boundary layer and a hypersonic flow over a double compression ramp.     

Experimental and numerical investigation of turbulent cross-flow in a staggered tube bundle

This paper presents the results of measurements and numerical predictions of turbulent cross-flow in a staggered tube bundle. The bundle consists of transverse and longitudinal pitch-to-diameter ratios of 3.8 and 2.1, respectively. The experiments were conducted using a particle image velocimetry technique, in a flow of water in a channel at a Reynolds number of 9300 based on the inlet velocity and the tube diameter. A commercial CFD code, ANSYS CFX V10.0, is used to predict the turbulent flow in the bundle. The steady and isothermal Reynolds–Averaged Navier–Stokes (RANS) equations were used to predict the turbulent flow using each of the following four turbulence models: a k-epsilon, a standard k-omega, a k-omega-based shear stress transport, and an epsilon-based second moment closure. The epsilon-based models used a scalable wall function and the omega-based models used a wall treatment that switches automatically between low-Reynolds and standard wall function formulations.

The experimental results revealed extremely high levels of turbulence production by the normal stresses, as well as regions of negative turbulence production. The convective transport by mean flow and turbulent diffusion were observed to be significantly higher than in classical turbulent boundary layers. As a result, turbulence production is generally not in equilibrium with its dissipation rate. In spite of these characteristics, it was observed that the Reynolds normal stresses approximated from the k-based two-equation models were in a closer agreement with experiments than values obtained from the second moment closure. The results show that none of the turbulence models was able to consistently reproduce the mean and turbulent quantities reasonably well. The omega-based models predicted the mean velocities better in the developing region while the epsilon-based models gave better results in the region where the flow is becoming spatially periodic.     

Effect of thermophoresis and other parameters on the particle deposition on a tilted surface

In this study a modified version of v2-f turbulence model (φ-α), is applied to simulate a non-isothermal air-flow. The φ-α model and a two-phase Eulerian approach complement each other to predict the rate of particle deposition on a tilted surface. The φ-α model can accurately calculate the normal fluctuations, which mainly represent the non-isotropic nature of turbulence regime near the wall. The Eulerian model was modified considering the most important mechanism in the particle deposition rate when compared to the experimental data. The model performance is examined by comparing the rate of particle deposition on a vertical surface with the experimental data in a turbulent channel flow available in the literature. The effects of lift force, turbophoretic force, thermophoreric force, electrostatic force, gravitational force and Brownian/turbulent diffusion were examined on the particle deposition rate. The results show that, using the φ-α model predicts the rate of deposition with reasonable accuracy. The results of modified particle model are in good agreement with the experimental data. This study highlights the paramount effect of thermophoretic force on the particle deposition rate and clearly shows that when the temperature difference exceeds a certain limit, the electrostatic force has insignificant effect on the particle deposition rate. Furthermore, it is indicated that even at small temperature differences, the effect of tilt angle on the particle deposition rate for intermediate-size particles is negligible.     

Numerical modeling of a two-dimensional vertical turbulent two-phase jet

The results of numerically modeling two-dimensional two-phase flow of the “gas-solid particles” type in a vertical turbulent jet are presented for three cases of its configuration, namely, descending, ascending, and without account of gravity. Both flow phases are modeled on the basis of the Navier-Stokes equations averaged within the framework of the Reynolds approximation and closed by an extended k-? turbulence model. The averaged two-phase flow parameters (particle and gas velocities, particle concentration, turbulent kinetic energy, and its dissipation) are described using the model of mutually-penetrating continua. The model developed allows for both the direct effect of turbulence on the motion of disperse-phase particles and the inverse effect of the particles on turbulence leading to either an increase or a decrease in the turbulent kinetic energy of the gas. The model takes account for gravity, viscous drag, and the Saffman lift. The system of equations is solved using a difference method. The calculated results are in good agreement with the corresponding experimental data which confirms the effect of solid particles on the mean and turbulent characteristics of gas jets.     

Approach to local axisymmetry in a turbulent cylinder wake

The objective of this experimental study is to characterise the small-scale turbulence in the intermediate wake of a circular cylinder using measured mean-squared velocity gradients. Seven of the twelve terms which feature in ε, the mean dissipation rate of the turbulent kinetic energy, were measured throughout the intermediate wake at a Reynolds number of Re _d  ≈ 3000 based on the cylinder diameter (d). Earlier measurements of the nine major terms of ε by Browne et al. (J Fluid Mech 179: 307–326 1987) at a downstream distance (x) of x = 420d and Re _d  ≈ 1170 are also used. Whilst departures from local isotropy are significant at all locations in the wake, local axisymmetry of the small-scale turbulence with respect to the mean flow direction is first satisfied approximately at x = 40d. The approach towards local axisymmetry is discussed in some detail in the context of the relative values of the mean-squared velocity gradients. The data also indicate that axisymmetry is approximately satisfied by the large scales at x/d ≥ 40, suggesting that the characteristics of the small scales reflect to a major extent those of the large scales. Nevertheless, the far-wake data of Browne et al. (1987) show a discernible departure from axisymmetry for both small and large scales.     

Scalar turbulence model investigation with variable turbulent Prandtl number in heated jets and diffusion flames

Traditional turbulence models using constant turbulent Prandtl number fail to predict the experimentally observed anisotropies in the thermal eddy diffusivity and thermal turbulent intensity fields. Accurate predictions depend strongly on the turbulence model employed. Consequently, the objective of this paper is to assess the performance of turbulence model with variable turbulent Prandtl number in predicting of thermal and scalar fields quantities. The model is applied to axisymmetric turbulent round jet with variable density and in turbulent hydrogen diffusion flames using the flamelet concept. The k − ɛ turbulence model is used in conjunction with thermal field; the model involves solving supplemental scalar equations for the temperature variance and its dissipation rate. The model predictions are compared with available experimental data for the purpose of validating model. In reacting cases, velocity and scalar (including temperature and mass fractions) predictions agree relatively well in the near field of the investigated diluted hydrogen flames.     

Efficient Generation of Inflow Conditions for Large Eddy Simulation of Street-Scale Flows

Using a numerical weather forecasting code to provide the dynamic large-scale inlet boundary conditions for the computation of small-scale urban canopy flows requires a continuous specification of appropriate inlet turbulence. For such computations to be practical, a very efficient method of generating such turbulence is needed. Correlation functions of typical turbulent shear flows have forms not too dissimilar to decaying exponentials. A digital-filter-based generation of turbulent inflow conditions exploiting this fact is presented as a suitable technique for large eddy simulations computation of spatially developing flows. The artificially generated turbulent inflows satisfy the prescribed integral length scales and Reynolds-stress-tensor. The method is much more efficient than, for example, Klein’s (J Comp Phys 186:652–665, 2003) or Kempf et al.’s (Flow Turbulence Combust, 74:67–84, 2005) methods because at every time step only one set of two-dimensional (rather than three-dimensional) random data is filtered to generate a set of two-dimensional data with the appropriate spatial correlations. These data are correlated with the data from the previous time step by using an exponential function based on two weight factors. The method is validated by simulating plane channel flows with smooth walls and flows over arrays of staggered cubes (a generic urban-type flow). Mean velocities, the Reynolds-stress-tensor and spectra are all shown to be comparable with those obtained using classical inlet-outlet periodic boundary conditions. Confidence has been gained in using this method to couple weather scale flows and street scale computations.