Signal-processing techniques for low signal-to-noise ratio laser Doppler velocimetry signals

A variety of methods have been developed to obtain acurate frequency estimates from laser Doppler velocimetry (LDV) signals. Rapid scanning and fiber optic LDV systems require robust methods for extracting accurate frequency estimates with computational efficiency from data with poor signal-to-noise ratios. These methods typically fall into two general categories, time domain parametric techniques and frequency domain techniques. The frequency domain approach is initiated by transforming the Doppler bursts into the frequency domain using the fast Fourier transform (FFT). From this basic transformation a variety of interpolation procedures (parabolic, Gaussian, and centroid fits) have been developed to optimize the frequency estimation accuracy. The time domain approaches are derived from the parametric form of a sinusoid. The estimation of constants in this relationship is performed to satisfy specific constraints, typically a minimization of a variance expression. A comparison of these techniques is presented using simulated signals and additive Gaussian and Poisson white noise. The statistical bias and random errors for each method are presented from 200 signal simulations at each condition. Frequency estimation via the FFT with zero-padding and a Gaussian interpolation scheme was found to produce the lowest bias and random errors.List of symbols A(z) eigenfilter or characteristic polynomial - a _{m + 1} eigenvector - f frequency, Hz - f normalized frequency f = f/f _s - _d Doppler frequency estimate - f _i frequency of FFT spectral bin - f ^s sampling frequency - N number of sample points in data set - P _i ith power spectral line from PSD - r _xx (i) autocorrelation coefficient for time lag i - R_{Mn + 1} autocorrelation matrix of order M+1 - T sampling period - f spectral resolution for FFT, f = 1/N t - t sampling interval     

Acoustic Lagrangian velocity measurement in a turbulent air jet

Measuring Lagrangian velocities in a turbulent flow is of a great interest for turbulence modeling. We report measurements made in an axisymmetric turbulent air jet at Reynolds number R _λ ≃ 320, using acoustical Doppler scattering. Helium-filled soap bubbles are used as Lagrangian tracers. We describe an experimental setup which allows the simultaneous measurement of the full three-component Lagrangian velocity and the longitudinal Eulerian one. Lagrangian velocity probability density functions (PDF) are found Gaussian, close to Eulerian ones. Velocity correlations are analysed as well as the statistical dependence between components.     

DNS Scrutiny of the <Emphasis Type="Italic">ζ</Emphasis>-<Emphasis Type="Italic">f</Emphasis> Elliptic-Relaxation Eddy-Viscosity Model in Channel Flows with a Moving Wall

In this paper we present a Direct Numerical Simulations (DNS) of channel flow with stationary and moving walls. Three cases, Poiseuille-type with U_W/U_b = 0.75, intermediate-type with U_W/U_b = 1.215, and Couette-type with U_W/U_b = 1.5 (U_W and U_b are the wall and the bulk velocity), were compared with the pure Poiseuille U_W/U_b = 0, at a bulk Reynolds number equal to 4,800 corresponding to Re _\uptau = 288_{\uptau} =288. The DNS results were used to scrutinize the capabilities of ζ-f eddy viscosity model (based on the elliptic relaxation concept) in reproducing the near-wall turbulence in non conventional flows where the shear stress structures are strongly different with respect to the cases used for models calibration. The ζ-f model (also in its basic formulation) demonstrated to have good prospects to reproduce the main phenomenology of such class of flows due to its built-in capabilities to account separately for the different (and opposite) near wall effects on turbulence: the damping due to viscosity and pressure reflection. The results of the computations demonstrated that standard ζ-f model can reasonably reproduce the phenomenology of these flows in terms of velocity and turbulent kinetic energy profiles and budgets.     

Joint Scalar versus Joint Velocity-Scalar PDF Simulations of Bluff-Body Stabilized Flames with REDIM

Two transported PDF strategies, joint velocity-scalar PDF (JVSPDF) and joint scalar PDF (JSPDF), are investigated for bluff-body stabilized jet-type turbulent diffusion flames with a variable degree of turbulence–chemistry interaction. Chemistry is modeled by means of the novel reaction-diffusion manifold (REDIM) technique. A detailed chemistry mechanism is reduced, including diffusion effects, with N ₂ and CO ₂ mass fractions as reduced coordinates. The second-moment closure RANS turbulence model and the modified Curl’s micro-mixing model are not varied. Radiative heat loss effects are ignored. The results for mean velocity and velocity fluctuations in physical space are very similar for both PDF methods. They agree well with experimental data up to the neck zone. Each of the two PDF approaches implies a different closure for the velocity-scalar correlation. This leads to differences in the radial profiles in physical space of mean scalars and mixture fraction variance, due to different scalar flux modeling. Differences are visible in mean mixture fraction and mean temperature, as well as in mixture fraction variance. In principle, the JVSPDF simulations can be closer to physical reality, as a differential model is implied for the scalar fluxes, whereas the gradient diffusion hypothesis is implied in JSPDF simulations. Yet, in JSPDF simulations, turbulent diffusion can be tuned by means of the turbulent Schmidt number. In the neck zone, where the turbulent flow field results deteriorate, the joint scalar PDF results are in somewhat better agreement with experimental data, for the test cases considered. In composition space, where results are reported as scatter plots, differences between the two PDF strategies are small in the calculations at hand, with a little more local extinction in the joint scalar PDF results.     

K-ε-T model of dense liquid-solid two-phase turbulent flow and its application to the pipe flow

To predict the characteristics of dense liquid-solid two-phase flow, K-ε-T model is established, in which the turbulent flow of fluid phase is described with fluid turbulent kinetic energy K_f and its dissipation rate ε_f, and the particles random motion is described with particle turbulent energy K_p and its dissipation rate ε_p and pseudothermal temperature T_p. The governing equations are also derived. With K-ε-T model, numerical study of dense liquid-solid two-phase turbulent up-flow in a pipe is performed. The calculated results are in good agreement with experimental data of Alajbegovic et al. (1994), and some flow features are captured.     

On model equations for particle dispersion in inhomogeneous turbulence

Comparisons are made between the Advection–Diffusion Equation (ADE) approach for particle transport and the two-fluid model approach based on the PDF method. In principle, the ADE approach offers a much simpler way of calculating the inertial deposition of particles in a turbulent boundary layer than that based on the PDF approach. However the ADE equations that have recently been used are only strictly valid for a simple Gaussian process when particle inertia is small. Using a prescribed, but in general non-Gaussian random particle velocity field, it is shown that the net particle mass flux contains a drift term in addition to that from the mean velocity of the particle velocity field, associated with the compressibility of the velocity field. Furthermore the diffusive flux in general depends not only upon the gradient of the mean concentration (true only for a Gaussian random flow field) but also upon higher order derivatives whose relative contribution depends on diffusion coefficients D_ijk… etc. These coefficients depend upon the statistical moments associated with random displacements and compressibility of the particle flow field along particle trajectories which in turn depend upon particle inertia. In contrast the PDF approach offers the advantage of using a simple gradient (Gaussian) approximation in particle phase space which can lead to a non-Gaussian spatial dispersion process when particle inertia is important. Conditions based on the particle mean free path are derived for which a simple ADE is appropriate. Some of the features of particle transport in an inhomogeneous turbulent flow are illustrated by examining particle dispersion in a random flow field composed of pairs of counter rotating vortices which has an rms velocity which increase linearly from a stagnation point.     

Flow Structure of Lifted Swirling Jet Flames

Swirling combustion is widely applied in various applications such as gas turbines, utility boilersor waste incinerators. This article contributes to the ongoing research by providing experimentaldata that are gathered in the mixing zone of a lifted swirling premixed natural gas flame. Theobjective of this paper is fivefold: (1) to introduce the lifted swirling flame featuring lowNO _x emissions (2) to provide experimental data such as major species distributions, temperature and streamlines of the flow pattern, (3) to report on velocity bias in probability density function (PDF) distributions and to present PDF sequences of velocities in medium scale swirling flows, (4) to make an assessment on the local small-scale turbulence that is present in the swirling mixinglayer and (5) to provide new experimental data for model verification and development.The PDFs are corrected in order to compensate for the velocity bias phenomenon, which is typicalfor randomly sampled LDA data. Sequences of axial PDF data are presented and measurement locationsof interest are selected to look at the PDF characteristics of the internal and externalrecirculation zones, the mixing layer and the onset of the reacting flow into detail. The mixinglayer PDFs covered a wide velocity range and revealed bimodality; even the concept ofmulti-modality is suggested and explored. Analysis showed that a sum of two Gaussian distributionscan accurately envelop the experimental PDFs. The reason for this broadband turbulence behavior isto be found in combination of precessing and flapping motion of the flow structures, and also incombustion generated instabilities of the lifted flame. As a result, the flame brush is wide (largescale motion) and the mixing (small-scale turbulence) flattens any high temperatures in thecombustion process.The multi-scale turbulence concept is subsequently used to make anassessment of the local turbulence characteristics in the mixing layer.The idea is that the PDFs capture both contributions of the flow-inherent fine grain turbulence (u _l ) which is superposed on slowlarge scale fluctuating structures. It is this u _l that will be of interest in continued research on the classification of the lifted flame into acombustion regime diagram (e.g. Borghi diagram). Finally, the bimodalitycharacter in reacting flows and the prediction of large-scale structuresmay be a challenge for LES researchers.     

Adaptive wall functions for the v2‐f turbulence model

Adaptive wall functions for the v²‐f turbulence model have been derived for the flow over a flat plate at zero pressure gradient. These wall functions were implemented via tables for the turbulence quantities and the friction velocity u_τ. A special treatment for the ε and f boundary conditions is proposed. On fine grids (y⁺<1) this approach yields results consistent with the wall integration solution. Detailed numerical results are presented for a zero pressure gradient boundary layer and separated flow over a ramp. Copyright © 2006 John Wiley & Sons, Ltd.     

Piezo rotary and axial vibrator (PRAV) characterization of a fresh coating during its drying

Using a piezo rotary and axial vibrator (PRAV) the viscoelastic properties of a fresh Al₂O₃ coating that is deposited on a well-chosen plate can be measured simultaneously during drying. There are three types of vibrations available to be excited and detected in their resonance modes, by measuring resonance frequencies f _k and half-widths h _k, before and after coating as a function of time: axial vibration, bending vibration to get the evolving Young’s modulus and rotary vibration to follow the viscosity increase during drying of the thin coating. This information is contained in the complex frequency shifts Ω_k = 2Δf _k /f _k + iΔh _k /f _k of the three vibration modes caused by the coating on the plate. A derivation of the relationships, their validation on Newtonian liquids and experimental applications carried out using the PRAV are given in the paper.     

Presumed Mapping Functions for Eulerian Modelling of Turbulent Mixing

The mapping closure of Chen et al. [Phys. Rev. Lett., 63, 1989] is a transported probability density function (PDF) method that has proven very efficient for modelling of turbulent mixing in homogeneous turbulence. By utilizing a Gaussian reference field, the solution to the mapping function (in homogeneous turbulence) can be found analytically for a range of initial conditions common for turbulent combustion applications, e.g. for binary or trinary mixing. The purpose of this paper is to investigate the possibility of making this solution a presumed mapping function (PMF) for inhomogeneous flows. The PMF in turn will imply a presumed mixture fraction PDF that can be used for a wide range of models in turbulent combustion, e.g. flamelet models, the conditional moment closure (CMC) or large eddy simulations. The true novelty of the paper, though, is in the derivation of highly efficient, closed algebraic expressions for several existing models of conditional statistics, e.g. for the conditional scalar dissipation/diffusion rate or the conditional mean velocity. The closed form expressions nearly eliminates the overhead computational cost that usually is associated with nonlinear models for conditional statistics. In this respect it is argued that the PMF is particularly well suited for CMC that relies heavily on manipulations of the PDF for consistency. The accuracy of the PMF approach is shown with comparison to DNS of a single scalar mixing layer to be better than for the β-PDF. Not only in the shape of the PDF itself, but also for all conditional statistics models computed from the PDF.     

Equilibrium states in homogeneous turbulence with baroclinic instability

The evolution of energies and fluxes in homogeneous turbulence with baroclinic instability is analyzed using the linear theory. The mean flow corresponds to a vertical shear having a uniform mean velocity gradient, ?U _i /?x _j  = S δ _i1 δ _j3, a system rotation about the vertical axis with rate Ω, Ω _i  = Ωδ _i3, and uniform buoyancy gradients in the spanwise ${(\partial B{/}\partial x_2\,{=}\, N_h^2\,{=}\,-2\Omega S)}The evolution of energies and fluxes in homogeneous turbulence with baroclinic instability is analyzed using the linear theory. The mean flow corresponds to a vertical shear having a uniform mean velocity gradient, ∂U _i /∂x _j  = S δ _i1 δ _j3, a system rotation about the vertical axis with rate Ω, Ω _i  = Ωδ _i3, and uniform buoyancy gradients in the spanwise (?B/?x₂ = N_h² = -2WS){(\partial B{/}\partial x_2\,{=}\, N_h^2\,{=}\,-2\Omega S)} and vertical (?B/?x₃ = N_v²){(\partial B{/}\partial x_3\,{=}\,N_v^2)} directions. Computations based on the rapid distortion theory (RDT) are performed for several values of the rotation number R = 2Ω/S and the Richardson number R_i = N_v²/S² < 1{R_i\,{=}\,N_v^2/S^2 <1 }. It is shown that, during an initial phase, the energies and the buoyancy fluxes are sensitive to the effects of pressure and viscosity. At large time, the ratios of energies, as well as the normalized fluxes, evolve to an asymptotically constant value, while the pressure–strain correlation scaled with the product of the turbulent kinetic energy by the shear rate approaches zero. Accordingly, an analytical parametric study based on the “pressure-less” approach (PLA) is also presented. The analytical study indicates that, when R _i  < 1, there is an exponential instability and equilibrium states of turbulence, in agreement with RDT. The energies and the buoyancy fluxes grow exponentially for large times with the same rate (γ in St units). The asymptotic value of the ratios of energies yielded by RDT is well described by its PLA counterpart derived analytically. At R _i  = 0, the asymptotic value of γ increases with increasing R approaching 2 for high rotation rates. At low rotation rates, an important contribution to the kinetic energy comes from the streamwise kinetic energy, whereas, at high rotation rates, the contribution of the vertical kinetic energy is dominant. When 0 < R _i  < 1 and R 1 0{R\ne 0}, the asymptotic value of γ decreases as R _i increases so as it becomes zero at R _i  = 1.     

Range-pair-counting fatigue gage

In this paper, the surface roughening phenomenon of an aluminum foil during fatigue process is utilized for monitoring of a variable amplitude stress. It is found that constant mean stress does not affect surface roughness. When the mean stress varies, however, it causes surface roughening. In such a case, the equivalent stress defined by (Σσ _i ^α N _i /ΣN _i )^1/α dominates surface roughness of an aluminum foil, where σ _i andN _i are stress amplitude and the number of cycles counted by the rangepair method, respectively, and α is a gage factor which is peculiar to the foil used. It is concluded that the aluminum foil can be used as the range-pair counting fatigue gage with high accuracy.     

Asymptotic Variational Wave Equations

We investigate the equation (u _t +(f(u)) _x ) _x =f ^{′ ′}(u) (u _x )²/2 where f(u) is a given smooth function. Typically f(u)=u ²/2 or u ³/3. This equation models unidirectional and weakly nonlinear waves for the variational wave equation u _tt − c(u) (c(u)u _x ) _x =0 which models some liquid crystals with a natural sinusoidal c. The equation itself is also the Euler–Lagrange equation of a variational problem. Two natural classes of solutions can be associated with this equation. A conservative solution will preserve its energy in time, while a dissipative weak solution loses energy at the time when singularities appear. Conservative solutions are globally defined, forward and backward in time, and preserve interesting geometric features, such as the Hamiltonian structure. On the other hand, dissipative solutions appear to be more natural from the physical point of view.We establish the well-posedness of the Cauchy problem within the class of conservative solutions, for initial data having finite energy and assuming that the flux function f has a Lipschitz continuous second-order derivative. In the case where f is convex, the Cauchy problem is well posed also within the class of dissipative solutions. However, when f is not convex, we show that the dissipative solutions do not depend continuously on the initial data.     

Study of airborne particles produced by normal impact of millimetric droplets onto a liquid film

The aim of this work is to carry out an experimental investigation into the generation of airborne microparticles when millimetric droplets of aqueous solutions impact onto a liquid film. Impact experiments using 3.9 mm diameter droplets were carried out for Weber numbers between 159 and 808, with a fixed Ohnesorge number of 2 × 10⁻³ and film parameters S _f (the ratio between the thickness of the liquid film h _film and the diameter of the impacting droplet d _i) between 0.3 and 1. Observed results show that the deposition/splashing threshold is independent of the parameter S _f in agreement with the data in the literature. The aerosol measurement results demonstrate the production of solid particles from the evaporation of secondary microdroplets with diameters less than 30 μm formed when splash occurs. The median diameter of these microdroplets is around 20 μm, corresponding to a value of d ₅₀/d _i = 5 × 10⁻³. Taken together, the results show that the mass and the number of particles emitted increase as the Weber number increases. Moreover, at a Weber number of 808, the results show that the mass and number of particles emitted increases as the parameter S _f decreases. In this case, the mean number of microdroplets emitted per impact is equal to 14 for S _f = 1 and equal to 76 for S _f = 0.3.     

Reynolds' analogy for the thermal convection driven by nonisothermal stretching surfaces

In the present paper the steady boundary-layer flows induced by permeable stretching surfaces with variable temperature distribution are investigated under the aspect of Reynolds' analogy r = St _x /C _f(x). It is shown that for certain stretching velocities and wall temperature distributions, “Reynolds' function”r, i.e. the ratio of the local Stanton number St _x and the skin friction coefficient C _f(x) equals −1/2 for any value of the Prandtl number Pr and of the dimensionless suction/injection velocity f _w. In all of these cases, the dimensionless temperature field ϑ is connected to the dimensionless downstream velocity f ^′ by the simple relationship ϑ=(f ^′)^Pr. It is also shown that in the general case, Reynolds' function r may possess several singularities in f _w. The largest of them represents a critical value, so that for f _w<f _w,crit the solutions of the energy equation (although they still satisfy all the boundary conditions) become nonphysical.     

Turbulent boundary-layer analysis using finite elements

The use of finite element methods for turbulent boundary-layer flow is relatively recent and of limited extent.¹ In the present study, we extend the group variable approach of Fletcher and Fleer^2,3 to treat turbulent boundary layer flows with heat transfer using a two-equation turbulence model. The main concepts in the formulations include a Dorodnitsyn-type transformation which uses a velocity component as the transverse variable, a ‘variational’ formulation for the transformed equations using special test functions and development of a two-equation turbulence model in terms of the turbulent kinetic energy and turbulence dissipation rate as additional field variables. Several numerical test cases have been examined comparing the results with finite difference calculations and comparing the two-equation turbulence model with an algebraic turbulence model.     

A photoelastic study of friction at multipoint contacts

A new method of measuring the normal and sliding loads associated with multiple-point contact is introduced. A multiple-point contact is modeled with a steel die with a profile that simulates a rough surface. A very large scale factor is used in modeling this surface. The steel die is placed in contact with a photoelastic model of a half plane and is subjected to a normal load. This normal load is partitioned over the multiple points of contact producing an isochromatic fringe pattern that describes the stress distribution in the local neighborhood of the contact points. A sliding load is then imposed on the model which destroys the symmetry of this fringe pattern. The fringe data in this pattern are sufficient to determine the local loadsP _i andQ _i and the local coefficient of frictionf _i =Q _i /P _i at each contact point. An overdeterministic method is introduced which gives the solution forP _i ,Q _i andf _i using many data points taken from the isochromatic pattern in the local neighborhood of the contacts.Paper was presented at the 1991 SEM Spring Conference on Experimental Mechanics held in Milwaukee, WJ on June 9–13.     

Sensitivity of probabilistic fracture mechanics analysis to variations in statistical parameters

Probability of failure (p_f) of a structure is usually calculated for a specified set of statistical parameters (mean, standard deviation, and probability distribution) that characterize random variables. This approach may not be efficient in cases where one would like to know the effect of variations in statistical parameters on the probability of failure. A method based on generating and analyzing randomly selected statistical parameters is proposed. The method consists of generating databases of mean and coefficient of variation (COV = mean/standard deviation) values of relevant fracture mechanics variables through a random process. The method was applied to surface cracks in a flat wide plate loaded under elastic conditions. Probability of failure was calculated for each database record using the first-order reliability method (FORM). Multiple linear regression analyses of the database records were performed with p_f as dependent variable and statistical parameters as independent variables. The predicted p_f values were in very good agreement with the directly calculated p_f values for the specified variations of statistical parameters (±10%, ±15%, and ±20%), except those for fracture toughness and tensile stress, where variations should be limited to ±10% and ±15% ranges.