Measurement-Integrated simulations and Kalman filter applied to a co-flowing jet

This paper deals with the experimental evaluation of a flow analysis system based on the integration between an under-resolved Navier-Stokes simulation and experimental measurements with the mechanism of feedback (referred to as Measurement-Integrated simulation), applied to the case of a planar turbulent co-flowing jet. The experiments are performed with inner-to-outer-jet velocity ratio around 2 and the Reynolds number based on the inner-jet heights about 10000. The measurement system is a high-speed PIV, which provides time-resolved data of the flow-field, on a field of view which extends to 20 jet heights downstream the jet outlet. The experimental data can thus be used both for providing the feedback data for the simulations and for validation of the MI-simulations over a wide region. The effect of reduced data-rate and spatial extent of the feedback (i.e. measurements are not available at each simulation time-step or discretization point) was investigated. At first simulations were run with full information in order to obtain an upper limit of the MI-simulations performance. The results show the potential of this methodology of reproducing first and second order statistics of the turbulent flow with good accuracy. Then, to deal with the reduced data different feedback strategies were tested. It was found that for small data-rate reduction the results are basically equivalent to the case of full-information feedback but as the feedback data-rate is reduced further the error increases and tend to be localized in regions of high turbulent activity. Moreover, it is found that the spatial distribution of the error looks qualitatively different for different feedback strategies. Feedback gain distributions calculated by optimal control theory are presented and proposed as a mean to make it possible to perform MI-simulations based on localized measurements only. So far, we have not been able to low error between measurements and simulations by using these gain distributions.

Simulation of spherical turbulent premixed flames

Outwardly-propagating spherical turbulent premixed flames are studied using Reynolds-Averaged Navier-Stokes methodology. The reaction rate closure is based on recently developed strained flamelets, where the flamelets are parametrised using the scalar dissipation rate. It is shown that the leading edge speed for the spherical flame is higher than for the planar case, for a given turbulence and thermo-chemical conditions. In addition, it is shown that including the mean curvature effects in the reaction rate closure do not influence the spherical flame speeds when compared with a model that excludes these effects.

Improved jet noise modeling using a new time-scale.

To calculate the noise emanating from a turbulent flow using an acoustic analogy knowledge concerning the unsteady characteristics of the turbulence is required. Specifically, the form of the turbulent correlation tensor together with various time and length-scales are needed. However, if a Reynolds Averaged Navier-Stores calculation is used as the starting point then one can only obtain steady characteristics of the flow and it is necessary to model the unsteady behavior in some way. While there has been considerable attention given to the correct way to model the form of the correlation tensor less attention has been given to the underlying physics that dictate the proper choice of time-scale. In this paper the authors recognize that there are several time dependent processes occurring within a turbulent flow and propose a new way of obtaining the time-scale. Isothermal single-stream flow jets with Mach numbers 0.75 and 0.90 have been chosen for the present study. The Mani-Gliebe-Balsa-Khavaran method has been used for prediction of noise at different angles, and there is good agreement between the noise predictions and observations. Furthermore, the new time-scale has an inherent frequency dependency that arises naturally from the underlying physics, thus avoiding supplementary mathematical enhancements needed in previous modeling.

Effect of flow-geometry on turbulence-scalar interaction in premixed flames

Turbulent combustion of stoichiometric hydrogen-air mixture is simulated using direct numerical simulation methodology, employing complex chemical kinetics. Two flame configurations, freely propagating and V-flames stabilized behind a hot rod, are simulated. The results are analyzed to study the influence of flame configuration on the turbulence-scalar interaction, which is critical for the scalar gradient generation processes. The result suggests that this interaction process is not influenced by the flame configuration and the flame normal is found to align with the most extensive strain in the region of intense heat release. The combustion in the rod stabilized flame is found to be flamelet like in an average sense and the growth of flame-brush thickness with the downstream distance is represented well by Taylor theory of turbulent diffusion, when the flame-brushes are non-interacting. The thickness is observed to saturate when the flame-brushes interact, which is found to occur in the simulated rod stabilized flame with Taylor micro-scale Reynolds number of 97. © 2011 American Institute of Physics.

Experimental study of homogenisation in grid turbulence

Turbulence statistics have been measured immediately downstream of a regular grid made of round rods with rod spacing M. 2D-2C PIV was used to analyse a measurement area of 14M x 4M in the down and cross-stream directions respectively. The relevant Reynolds number span the range Re M = U ∞M/ν = 5 500 - 16 500. The Reynolds shear stresses recorded on two parallel measurement planes differently located relative to the grid exhibit significant discrepancies over the first 5M, but have completely homogenised in the cross-stream direction by x/M = 7. The downstream evolution of the two-point velocity correlation functions shows a progressive loss of coherence and a clear trend towards the expected isotropic behavior. The same conclusions apply to measurements taken in the wake of another regular grid made of square rods. Changes in the vortex shedding pattern from the grid were observed at the lowest Reynolds number, with two of the four rod wakes captured shedding in phase with each other but in anti-phase with a third one. The impact of this early flow coherence on the turbulence statistics did not persist due to the homogenisation process.

The non-local nature of structure functions

Kolmogorov's two-thirds, ((Δv) 2) ∼ e 2/ 3r 2/ 3, and five-thirds, E ∼ e 2/ 3k -5/ 3, laws are formally equivalent in the limit of vanishing viscosity, v → 0. However, for most Reynolds numbers encountered in laboratory scale experiments, or numerical simulations, it is invariably easier to observe the five-thirds law. By creating artificial fields of isotropic turbulence composed of a random sea of Gaussian eddies whose size and energy distribution can be controlled, we show why this is the case. The energy of eddies of scale, s, is shown to vary as s 2/ 3, in accordance with Kolmogorov's 1941 law, and we vary the range of scales, γ = s max/s min, in any one realisation from γ = 25 to γ = 800. This is equivalent to varying the Reynolds number in an experiment from R λ = 60 to R λ = 600. While there is some evidence of a five-thirds law for g > 50 (R λ > 100), the two-thirds law only starts to become apparent when g approaches 200 (R λ ∼ 240). The reason for this discrepancy is that the second-order structure function is a poor filter, mixing information about energy and enstrophy, and from scales larger and smaller than r. In particular, in the inertial range, ((Δv) 2) takes the form of a mixed power-law, a 1+a 2r 2+a 3r 2/ 3, where a 2r 2 tracks the variation in enstrophy and a 3r 2/ 3 the variation in energy. These findings are shown to be consistent with experimental data where the polution of the r 2/ 3 law by the enstrophy contribution, a 2r 2, is clearly evident. We show that higherorder structure functions (of even order) suffer from a similar deficiency.

Velocity profiles in a cylindrical liquid jet by reconstructed velocimetry

An experimental setup and a simple reconstruction method are presented to measure velocity fields inside slightly tapering cylindrical liquid jets traveling through still air. Particle image velocimetry algorithms are used to calculate velocity fields from high speed images of jets of transparent liquid containing seed particles. An inner central plane is illuminated by a laser sheet pointed at the center of the jet and visualized through the jet by a high speed camera. Optical distortions produced by the shape of the jet and the difference between the refractive index of the fluid and the surrounding air are corrected by using a ray tracing method. The effect of the jet speed on the velocity fields is investigated at four jet speeds. The relaxation rate for the velocity profile downstream of the nozzle exit is reasonably consistent with theoretical expectations for the low Reynolds numbers and the fluid used, although the velocity profiles are considerably flatter than expected. © 2012 American Society of Mechanical Engineers.

Turbulent flow structure in experimental laboratory wind-generated gravity waves

This paper is the third part of a report on systematic measurements and analyses of wind-generated water waves in a laboratory environment. The results of the measurements of the turbulent flow on the water side are presented here, the details of which include the turbulence structure, the correlation functions, and the length and velocity scales. It shows that the mean turbulent velocity profiles are logarithmic, and the flows are hydraulically rough. The friction velocity in the water boundary layer is an order of magnitude smaller than that in the wind boundary layer. The level of turbulence is enhanced immediately beneath the water surface due to micro-breaking, which reflects that the Reynolds shear stress is of the order u *w 2. The vertical velocities of the turbulence are related to the relevant velocity scale at the still-water level. The autocorrelation function in the vertical direction shows features of typical anisotropic turbulence comprising a large range of wavelengths. The ratio between the microscale and macroscale can be expressed as λ/Λ=a Re Λ n, with the exponent n slightly different from -1/2, which is the value when turbulence production and dissipation are in balance. On the basis of the wavelength and turbulent velocity, the free-surface flows in the present experiments fall into the wavy free-surface flow regime. The integral turbulent scale on the water side alone underestimates the degree of disturbance at the free surface. © 2012 Elsevier B.V.

Breaking of symmetry in microfluidic propulsion driven by artificial cilia.

In this Brief Report we investigate biomimetic fluid propulsion due to an array of periodically beating artificial cilia. A generic model system is defined in which the effects of inertial fluid forces and the spatial, temporal, and orientational asymmetries of the ciliary motion can be individually controlled. We demonstrate that the so-far unexplored orientational asymmetry plays an important role in generating flow and that the flow increases sharply with Reynolds number and eventually becomes unidirectional. We introduce the concept of configurational symmetry that unifies the spatial, temporal, and orientational symmetries. The breaking of configurational symmetry leads to fluid propulsion in microfluidic channels.

Assessment of computational fluid dynamics and experimental data for shock boundary-layer interactions

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop, numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric, and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and, in general, it was difficult to discern clear trends in the data. For the Reynolds-averaged Navier-Stokes (RANS) methods, the choice of turbulence model appeared to be the largest factor in solution accuracy. Scale-resolving methods, such as large-eddy simulation (LES), hybrid RANS/LES, and direct numerical simulation, produced error levels similar to RANS methods but provided superior predictions of normal stresses. Copyright © 2012 by Daniella E. Raveh and Michael Iovnovich.

2D PIV measurements in the near field of grid turbulence using stitched fields from multiple cameras

We present measurements of grid turbulence using 2D particle image velocimetry taken immediately downstream from the grid at a Reynolds number of Re M = 16500 where M is the rod spacing. A long field of view of 14M x 4M in the down- and cross-stream directions was achieved by stitching multiple cameras together. Two uniform biplanar grids were selected to have the same M and pressure drop but different rod diameter D and crosssection. A large data set (10 4 vector fields) was obtained to ensure good convergence of second-order statistics. Estimations of the dissipation rate ε of turbulent kinetic energy (TKE) were found to be sensitive to the number of meansquared velocity gradient terms included and not whether the turbulence was assumed to adhere to isotropy or axisymmetry. The resolution dependency of different turbulence statistics was assessed with a procedure that does not rely on the dissipation scale η. The streamwise evolution of the TKE components and ε was found to collapse across grids when the rod diameter was included in the normalisation. We argue that this should be the case between all regular grids when the other relevant dimensionless quantities are matched and the flow has become homogeneous across the stream. Two-point space correlation functions at x/M = 1 show evidence of complex wake interactions which exhibit a strong Reynolds number dependence. However, these changes in initial conditions disappear indicating rapid cross-stream homogenisation. On the other hand, isotropy was, as expected, not found to be established by x/M = 12 for any case studied. © Springer-Verlag 2012.

Oscillatory sensitivity patterns for global modes in wakes

Globally unstable wakes with co-flow at intermediate Reynolds numbers are studied, to quantify important spatial regions for the development and control of the global instability. One region of high structural sensitivity is found close to the inlet for all wakes, in agreement with previous findings for cylinder wakes. A second, elongated region of high structural sensitivity is seen downstream of the first one for unconfined wakes at Re = 400. When base flow modifications are considered, a spatially oscillating sensitivity pattern is found inside the downstream high structural sensitivity region. This implies that the same change in the base flow can either destabilize or stabilize the flow, depending on the exact position where it is applied. It is shown that the sensitivity pattern remains unchanged for different choices of streamwise boundary conditions and numerical resolution. Actual base flow modifications are applied in selected configurations, and the linear global modes recomputed. It is confirmed that the linear global eigenvalues move according to the predicted sensitivity pattern for small amplitude base flow modifications, for which the theory applies. We also look at the implications of a small control cylinder on the flow. Only the upstream high sensitivity region proves to be robust in terms of control, but one should be careful not to disturb the flow in the downstream high sensitivity region, in order to achieve control. The findings can have direct implications on the numerical resolution requirements for wakes at higher Reynolds numbers. Furthermore, they provide one more possible explanation to why confined wakes have a more narrow frequency spectrum than unconfined wakes.

Global nonlinear dynamics of thin aerofoil wakes

In the present investigation of thin aerofoil wakes we compare the global nonlinear dynamics, obtained by direct numerical simulations, to the associated local instability features, derived from linear stability analyses. A given configuration depends on two control parameters: the Reynolds number Re and the adverse pressure gradient m (with m < 0) prevailing at the aerofoil trailing edge. Global instability is found to occur for large enough Re and |m|; the naturally selected frequency is determined by the local absolute frequency prevailing at the trailing edge. © 2010 Springer Science+Business Media B.V.

Study of the turbulence in the air-side and water-side boundary layers in experimental laboratory wind induced surface waves

This study detailed the structure of turbulence in the air-side and water-side boundary layers in wind-induced surface waves. Inside the air boundary layer, the kurtosis is always greater than 3 (the value for normal distribution) for both horizontal and vertical velocity fluctuations. The skewness for the horizontal velocity is negative, but the skewness for the vertical velocity is always positive. On the water side, the kurtosis is always greater than 3, and the skewness is slightly negative for the horizontal velocity and slightly positive for the vertical velocity. The statistics of the angle between the instantaneous vertical fluctuation and the instantaneous horizontal velocity in the air is similar to those obtained over solid walls. Measurements in water show a large variance, and the peak is biased towards negative angles. In the quadrant analysis, the contribution of quadrants Q2 and Q4 is dominant on both the air side and the water side. The non-dimensional relative contributions and the concentration match fairly well near the interface. Sweeps in the air side (belonging to quadrant Q4) act directly on the interface and exert pressure fluctuations, which, in addition to the tangential stress and form drag, lead to the growth of the waves. The water drops detached from the crest and accelerated by the wind can play a major role in transferring momentum and in enhancing the turbulence level in the water side.On the air side, the Reynolds stress tensor's principal axes are not collinear with the strain rate tensor, and show an angle α σ≈=-20°to-25°. On the water side, the angle is α σ≈=-40°to-45°. The ratio between the maximum and the minimum principal stresses is σ a/σ b=3to4 on the air side, and σ a/σ b=1.5to3 on the water side. In this respect, the air-side flow behaves like a classical boundary layer on a solid wall, while the water-side flow resembles a wake. The frequency of bursting on the water side increases significantly along the flow, which can be attributed to micro-breaking effects - expected to be more frequent at larger fetches. © 2012 Elsevier B.V.