LIBRARY Michigan State University PLACE IN RETURN BOX to remove this checkout from your record. To AVOID FINES return on or before date due. MAY BE RECALLED with earlier due date if requested. DATE DUE DATE DUE DATE DUE 6/01 c:/ClRC/DateOuo.p65-p.15 __——— Stochastic Estimation of the Flow Structure Downstream of a Separating/Reattaching Flow Region Using Wall-Pressure Array Measurements By Mohamed Ibrahim Daoud A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Mechanical Engineering 2004 ABSTRACT Stochastic Estimation of the Flow Structure Downstream of a Separating/ Reattaching Flow Region Using Wall-Pressure Array Measurements By Mohamed Ibrahim Daoud This study examines the spatio-temporal characteristics of the surface-pressure fluctuations and associated flow structures in the developing flow downstream of the reattachment point of a fence-with-splitter-plate flow. The investigation focuses on understanding the wall-pressure field characteristics, and the flow sources responsible for its generation in the non-equilibrium boundary layer originating from the separating/reattaching shear layer associated with the flow over the fence, using a wall- pressure database that was simultaneously acquired with X-hotwire time series. This is motivated by guiding efforts to predict and/or control flow-induced noise and vibration in applications involving flows downstream of appendages and surface protrusions. Characterization of the wall-pressure data alone showed that the wall-pressure fluctuations were dominated by large-scale downstream-traveling disturbances that were generated upstream in the separated shear layer. Notwithstanding this dominance, the p' signature of these structures decayed with increasing downstream distance as the vortices underwent a relaxation process while the contribution of eddies, associated with the development of a "sub-boundary layer", became more significant with increasing downstream distance. In addition, wavenumber-frequency-spectrum results showed that pressure signatures of all wavenumbers and frequencies were associated with flow disturbances that travel downstream with the same convection velocity. Finally, multi-point Linear Stochastic Estimation of the flow field based on instantaneous wall-pressure information confirmed the dominance of wall-pressure generation by the passage of the outer-shear layer vortical structures and their mutual interaction. Examination of the linear source term in Poisson's equation of the pressure in conjunction with the stochastically-estimated velocity field revealed two mechanisms for p' generation associated with the quasi-periodic vortex passage. One mechanism was related to sources localized at the height of the vortex centers in the outer-shear layer, which represented direct generation of p' caused by the strong vortex-induced disturbances in the outer part of the flow. The other was located near the wall (y/6 < 0.2, where 5 is the boundary layer thickness) and resulted from interaction of the weak near- wall disturbances generated by the vortex passage with the strong mean wall shear. It is important to realize that while it is likely that the former mechanism could be modeled by the wall-pressure field associated with quasi-periodic passing of vortices embedded in uniform inviscid flow, the latter mechanism requires proper account of viscous effects near the wall and associated mean shear. DEDICATED TO MY MOTHER AND THE MEMORY OF MY FATHER iv ACKNOWLEDGMENT I would like to deeply acknowledge my family, advisor, committee members and colleagues, because this thesis work could not have been done without their support and/or contribution. I owe a great debt of gratitude and appreciation to my little family, Nermeen, Mariam and Omar, who have had selflessly sacrificed to stand by and assist me by any means to overcome the toughest moments when I became stuck between a rock and a hard place. My lovely family, whose age is the same as that of this thesis work, has made my whole life much easier due to their understanding and support. Also, I would like to express my love and respect to my mom who also has greatly suffered because of my stay away from her. I would like also to thank my sisters and brother, my brothers in law, and my nieces and nephews for their love and support. I would like to express my special appreciation to my advisor, Dr. Ahmed Naguib, whose guidance and directions have been greatly supportive of my work. I wish also to emphasize that his warm friendship and brotherhood-like relationship with his students and me myself has been the thrust of my Ph.D. rocket. I would like also to thank my committee members Dr. Giles Brereton, Dr. Dennis Gilliland and Dr. Manooch Koochesfahani whose advice and recommendations were greatly helpful throughout my thesis work. Finally, I wish to thank my colleagues whom I have worked with at the Flow Physics and Control Lab. This includes Laura Hudy, Yongxiang Li, Antonious Aditjandra, and Chad Stimson. This work has been partially supported by NSF grant CTS-0116907. TABLE OF CONTENTS LIST OF FIGURES ............................................................................ viii LIST OF TABLES ............................................................................. xv NOMENCALTURE ............................................................................ xvi 1. INTRODUCTION ........................................................................... 1 1.1. Literature Review ................................................................ 3 1.1.1. Separating, Reattaching and Developing Flows .................. 3 1.1.2. Characteristics of the Turbulent Wall Pressure ................... 13 1.1.3. Stochastic Estimation ................................................. 18 1.2. Motivation ......................................................................... 22 1.3. Objectives ......................................................................... 23 2. EXPERIMENTAL SET-UP AND METHODOLOGY ................................ 25 2.1. Experimental Set-up ............................................................. 25 2.1.1. Wind Tunnel Facility ................................................. 25 2.1.2. Test Model ............................................................. 28 2.1.3. Instrumentation ........................................................ 31 2.2. Experimental Methodology ..................................................... 40 2.2.1. Static Pressure System ................................................ 40 2.2.2. Microphones ........................................................... 42 2.2.3. Tufts ..................................................................... 50 2.2.4. Hotwire Sensors ........................................................ 55 3. WALL-PRESSURE RESULTS ........................................................... 76 3.1. Mean—pressure Distribution and the Reattachment Length ................. 77 3.2. Fluctuating Pressure Distribution .............................................. 79 3.3. Autocorrelation ................................................................... 82 3.4. Power Spectra ..................................................................... 87 3.5. Cross Correlation ................................................................. 96 3.6. Wavenumber—Frequency Spectrum ............................................ 103 4. VELOCITY-PRESSURE ANALYSIS AND STOCHASTIC ESTIMATION ............................................................................... 107 4.1. Velocity Characteristics ......................................................... 107 4.1.1. Characteristics of the Boundary-Layer Mean and Turbulent Velocity Profiles ........................................... 107 4.1.2. Velocity Spectra ....................................................... 123 4.2. Velocity-Pressure Characteristics .............................................. 134 4.2.1. Velocity/Wall-Pressure Cross-Spectra ............................. 134 4.2.2. Velocity/Wall-Pressure Cross-Correlations ....................... 145 4.2.3. Conditionally-Averaged-Data Analysis ............................ 156 4.3. Stochastic Estimation ............................................................ 165 vi 4.3.1. Principle and Equations .............................................. 167 4.3.2. Comparison between ML, SL and SQ .............................. 174 4.3.3. Stochastic Estimation Results and Wall-Pressure Sources ...... 184 5. CONCLUSIONS AND RECOMMENDATION FOR FUTURE WORK ........... 207 5.1 . Conclusions ........................................................................ 207 5.2. Future Work ....................................................................... 211 A. Derivation of Equation (2.6) ............................................................... 213 B. Data Acquisition Settings .................................................................. 214 REFERENCES .................................................................................. 216 vii Figure 1.1. Figure 2.1. I.) I) , 2' '. ( 1‘4; c. -1 (1; (I I.) U) ._.' I. r. a ._., (b Figure 2.4. J! Figure 2.-. Figure 2.6.. Figure 2. Figure 2.3. Figure 2.9. Figure 2.1. Figure 2.1 Figure 2. ii Figure 2.1 Figure 2.1J FightelH FiEllie 7 1. Figure 1.1. Figure 2.1. Figure 2.2. Figure 2.3. Figure 2.4. Figure 2.5. Figure 2.6. Figure 2.7. Figure 2.8. Figure 2.9. Figure 2.10. Figure 2.11. Figure 2.12. Figure 2.13. Figure 2.14. Figure 2.15. Figure 2.16. Figure 2.17. Figure 2.18. Figure 2.19. Figure 2.20. Figure 2.21. LIST OF FIGURES A schematic of an ideal two-dimensional flow over a splitter-plate-with—fence ....................................................... 2 Schematic of the wind facility (dimensions in meters) .................... 26 A schematic of the test model ................................................ 28 A frontal picture of the model showing the blockage effect of the fence ....................................................................... 30 A schematic of the instrument and middle plates depicting the locations of the pressure taps and microphones ........................ 30 A picture of the fence-with-splitter-plate model inside the wind tunnel .................................................................. 31 An image of a tufi attached to the I-plate .................................... 33 An image ofone ofthe microphones used in the array... . 35 A schematic of the microphone driving circuit .............................. 35 A schematic of the X-probe .................................................... 36 Images showing top and side views of the X-probe ........................ 37 Images of the traversing mechanism and the LVDT: overall (left) and close-up (right) views ...................................... 39 LVDT calibration ............................................................... 39 Mean-pressure-coefficient distribution for the top and bottom sides of the splitter plate .......................................... 42 Microphone installation in the Instrument and Middle plates .............................................................................. 43 A picture of the microphone-calibration setup .............................. 45 PWT calibration results: pressure ratio (top) and phase shift (bottom) ............................................................. 47 Frequency response of a microphone/hole assembly ....................... 48 Mean sensitivity (Km) of themicrophones (top) and Corresponding time-delay (bottom). Microphone #1 is the most upstream microphone ............................................ 51 Cross-talk check of Microphone # 13 (Mic. #1 is the most upstream one) ...................................................... 52 A sample wide-aperture tufl image ........................................... 54 FFP of the reattaching flow versus the x-location of the tuft . . . . 55 viii figm“ hguel‘ Figure 2.22. Figure 2.23. Figure 2.24. Figure 2.25. Figure 2.26. Figure 2.27. Figure 2.28. Figure 2.29. Figure 2.30. Figure 2.31. Figure 2.32. Figure 2.33. Figure 2.34. Figure 2.35. Figure 3.1. Figure 3.2. Figure 3.3. Figure 3.4. Figure 3.5. Figure 3.6. A schematic of the X-wire probe and angles definition .................... 57 An image of the yaw calibration setup ....................................... 59 Schematic of the yaw calibration setup ...................................... 62 Yaw calibration: voltage results ............................................... 62 Sample calibration data for the X-probe hotwires .......................... 63 Yaw calibration data after conversion of voltage to velocity ............. 63 Yaw-calibration data and corresponding curve fits for hotwire #1 (top plot) and hotwire #2 (bottom plot) ................................... 64 A photograph of the X-wire probe and its reflection in the splitter plate ...................................................................... 66 A close—up image of the X sensor ............................................. 67 Leading edge for X- and single-wire in-situ calibration ................... 69 Boundary-layer streamwise mean-velocity profiles from single- and X-wire measurements at x/Xr = 3.05 ............................ 71 Boundary-layer RMS streamwise-velocity profiles from single- and X-wire measurements at x/Xr = 3.05 ............................ 72 Pressure spectra of the microphone immediately upstream (x/XIr = 3.0) of the X—wire for three y locations of the wire ................ 74 Pressure spectra of the microphone immediately downstream (x/Xr = 3.11) ofthe X-wire for three y locations ofthe wire 75 Streamwise distribution of the mean pressure coefficient from the present study compared to Hudy et al. (2003) .............................. 78 Streamwise distribution of the coefficient of the RMS pressure fluctuation ........................................................................ 80 A comparison between the va streamwise distribution of Hudy et al. (2003) and the present work with that of Farabee and Casarella (1986) ................................................................. 81 Autocorrelation coefficient at five different streamwise locations downstream of the reattachment region for Re = 7600 and 15700 ............................................................ 84 Contour maps of the autocorrelation coefficient for all 16 Microphones and the two Reynolds numbers; Re = 7600 and 15700 ........................................................................ 85 A comparison of the auto-correlation results at zero time delay and five different streamwise locations downstream of the reattachment region for Re = 7600 .................................... 86 ix Figure 3.1 Figure 3.1 Figure .1_ ”we 4. iERre 4 Figure 3.7. Figure 3.8. Figure 3.9. Figure 3.10. Figure 3.11. Figure 3.12. Figure 3.13. Figure 3.14. Figure 3.15. Figure 3.16. Figure 3.17. Figure 3.18. Figure 4.1. Figure 4.2. Figure 4.3. Figure 4.4. Figure 4.5. Figure 4.6 Figure 4.7. Figure 4.8. Figure 4.9. A full contour map of the p’ autocorrelation coefficient for the separating/reattaching (Hudy et al. 2003) and the present recovering flows ............................................................. Spectra of the wall-pressure fluctuation for Re = 7600: logarithmic (top) and semi-logarithmic (bottom) plots Spectra of the wall-pressure fluctuation for Re = 15700: logarithmic (top) and semi-logarithmic (bottom) plots . Definition of the frequency used in the splitting procedure Of p' spectrum ............................................................... The ratio between hi gh- and low-frequency pressure fluctuations energy .......................................................... Cross-correlation results at five different locations for Re = 7600 and 15700 ........................................................ Cross-correlation results for all 16 microphones and Re = 7600 (tOp); and 15700 (bottom) ..................................... Plot for extraction of the convection velocity for the two Reynolds numbers of 7600 and 15700 ................................... Contour map of the cross-correlation coefficient for all 16 microphones and Re = 7600 ........................................... Contour map of the cross-correlation coefficient for all 16 microphones and Re = 15700 .......................................... Wavenumber-frequency spectrum for Re = 7600 ....................... Wavenumber-frequency spectrum from the separating/ reattaching flow of Hudy et al. (2003) ................................... Wall scaling of boundary-layer mean-velocity profile ................. Shear-layer scaling of the boundary-layer mean-velocity profile ........ 115 ... 117 119 ... 120 121 Wall scaling of the boundary-layer ulrms profile .......................... Shear-layer scaling of the boundary-layer urms profile .................. Wall scaling of the boundary-layer vrms profile ......................... Shear-layer scaling of the boundary-layer Vrms profile .................. Wall scaling of the boundary-layer Reynolds stress profile ........... Shear-layer scaling of the boundary-layer Reynolds stress profile ......................................................................... Streamwise-velocity spectra at y/8 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 and g = 2.05 ................................................. 87 89 90 94 95 98 99 100 102 102 104 106 110 114 123 125 Figure 4.1 Figure 4.1 Figure 4.1 Figure 4.2 Figure 4.". Figure 4.1 Figure 4.1 Figure 4.1 Figure 4.1 Figure 41 Flgllre 4.: FigUre 4.: Figure 4.10. Figure 4.11. Figure 4.12. Figure 4.13. Figure 4.14. Figure 4.15. Figure 4.16. Figure 4.17. Figure 4.18. Figure 4.19. Figure 4.20. Figure 4.21. Figure 4.22. Figure 4.23. Semi-log plots of the streamwise-velocity spectra at y/6 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 and E = 2.05 ............................. 128 Normal-velocity spectra at y/5 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 and E = 2.05 .......................................................... 128 Semi-log normal-velocity spectra at y/6 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 at E = 2.05 .............................................. 129 Boundary-layer velocity cross spectra at y/8 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 and E = 2.05 ...................................... 131 Semi-log plot of the boundary layer velocity cross spectra at y/6 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 and E = 2.05 132 Cross-spectra between the streamwise velocity (at E = 2.05) and wall-pressure (at E = 1.33) for y/5 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ............................................................ 137 Semi-log plot of the cross-spectra between the streamwise velocity (at E = 2.05) and wall-pressure (at E: 1.33) for y/5= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ................................ 137 Cross-spectra between the streamwise velocity (at E = 2.05) and wall-pressure (at E = 2.0) for y/S = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ............................................................ 138 Semi-log plot of the cross-spectra between the streamwise velocity (at E = 2.05) and wall-pressure (at E = 2.0) for y/5= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ................................ 138 Cross-spectra between the streamwise velocity (at E = 2.05) and wall-pressure (at E = 2.33) for y/8 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ............................................................ 139 Semi-log plot of the cross-spectra between the streamwise velocity (at E = 2.05) and wall-pressure (at E = 2.33) for y/5= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ................................ 139 Cross-spectra between the normal velocity (at E = 2.05) and wall-pressure (at E = 1.33) for y/5 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ..................................................................... 141 Semi-log plot of the cross-spectra between the normal velocity (at E = 2.05) and wall-pressure (at E = 1.33) at for y/8 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ............................................. 141 Cross-spectra between the normal velocity (at E= 2.05) and wall-pressure (at E = 2.0) for y/6= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ..................................................................... 143 xi Figure 4.: Figure 4.2 Figure 4.2 Figure 4.: Figure 4.: Figure 4.2 Figure 4.3 Figure 4.3 FiElite 4.3 Figure 4.24. Figure 4.25. Figure 4.26. Figure 4.27. Figure 4.28. Figure 4.29. Figure 4.30. Figure 4.31. Figure 4.32. Figure 4.33. Figure 4.34. Figure 4.35. Semi-log plot of the cross-spectra between the normal velocity (at E= 2.05) and wall-pressure (at E= 2.0) for y/5= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ..................................................... 143 Cross-spectra between the normal velocity (at E: 2.05) and wall-pressure (at E= 2.33) for y/8= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ..................................................................... 144 Semi-log plot of the cross-spectra between the normal velocity (at E= 2.05) and wall-pressure (at E= 2.33) at for y/6= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ............................................. 144 Contour plots of the cross—correlation coefficient at zero time shift between the streamwise (top plot) and normal (bottom plot) velocity and the wall-pressure at the location of the ten downstream most microphones (fi'om E = 1.33 to 2.33) 147 Maps of the cross-correlation coefficient between the streamwise velocity at E = 2.05 and wall-pressure at E= 1.33 (top plot), E= 2.0 (middle plot) and E: 2.33 (bottom plot) ....................................... 151 Maps of the cross-correlation coefficient between the normal velocity at E = 2.05 and wall-pressure at E= 1.33 (top plot), E= 2.0 (middle plot) and E= 2.33 (bottom plot) .............................. 152 The averaged cross-correlation between the normal-velocity and wall-pressure in the range of y/f) = 0.375 to 0.625 for E=1.67,1.89 and2.11 .......................................................... 155 Plot for extraction of the convection velocity from Rvp results ........... 155 Contour plots of the cross-correlation coefficient at zero- time shift between the low-pass filtered streamwise (top plot) and normal (bottom plot) velocity and the wall- pressure at the ten downstream most microphones (from E= 1.33 to 2.33) .......................................................... 158 Maps of the cross-correlation coefficient between the low- pass filtered streamwise velocity at E = 2.05 and wall- pressure at E = 1.33 (top plot), E= 2.0 (middle plot) and E= 2.33 (bottom plot) ............................................................ 159 Maps of the cross-correlation coefficient between the low- pass filtered normal velocity at E = 2.05 and wall-pressure at E= 1.33 (top plot), E: 2.0 (middle plot) and E = 2.33 (bottom plot) ..................................................................... 160 Wall-pressure,
, and conditionally-averaged mean-removed velocity-vector field associated with positive (top plot) and negative (bottom plot) wall-pressure peaks at E = 2.05 .................... 161 xii Figure 4.3 Figure 4.4 Figure 4.; Hgm44 HE‘ure 4.4 Figllre 4.; Figure 44 503644 (I Figure 4.36. Figure 4.37. Figure 4.38. Figure 4.39. Figure 4.40. Figure 4.41. Figure 4.42. Figure 4.43 Figure 4.44. Figure 4.45. Figure 4.46. Comparison between the mean streamwise-velocity profile and the conditionally-averaged streamwise-velocity profiles associated with positive and negative wall-pressure peaks ............................. 164 Wall-pressure,
, and conditionally-averaged velocity-vector
Field associated with positive (top plot) and negative (bottom
plot) wall-pressure peaks, viewed in a frame of reference
moving with 0.8 ano ............................................................ 166
Correlation coefficient between the measured and stochastically-
estimated streamwise- (top plot) and normal- (bottom plot)
velocity using ML, SL1, SL2 and SQ2 ........................................ 176
Comparison between the measured and stochastically-estimated
streamwise-velocity spectra using ML, SL1, SL2 and SQ2
for y/5 = 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625.. .................................. 179
Comparison between the measured and stochastically-estimated
normal-velocity spectra using ML, SL1, SL2 and SQ2
for y/6= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ............................. 180
Comparison of the conditionally-averaged and stochastically-
estimated streamwise-velocity using SL2 and SQ2 for y/8 =
0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ...................................... 182
Comparison of the conditionally-averaged and stochastically-
estimated normal-velocity using SL2 and SQ2 for y/5 =
0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 ....................................... 183
Wall-pressure time series at E = 2.0 (bottom of each plot) and
the associated stochastically-estimated velocity vector field (top
of each plot) using ML for three consecutive time windows,
viewed in a frame of reference translating with velocity of
0.81Uco (largest vector = 0.25Uw) ............................................. 187
The stochastically—estimated velocity vector field using
ML viewed in frames of reference translating with three
different velocities of 0.76U00 (top plot), 0.81Uco (middle
plot) and 0.86U..o (bottom plot); largest vector = 0.25Uco
(note that the "x" and "0" locations are fixed for all three
plots to help identify the variability in the vortex and
saddle-point locations, respectively) .......................................... 189
Vorticity and associated stochastically—estimated velocity
vector field at three consecutive time windows, viewed in
a frame of reference translating with 0.81Uno ................................. 193
Pseudo-instantaneous linear pressure source (q L) flooded
contour maps and associated wall-pressure for three
consecutive time windows ...................................................... 198
xiii
Figure 4.-
Figure 4.-
Figure 4.-
Figure 4.7”
Figure 4.47.
Figure 4.48.
Figure 4.49.
Figure 4.50.
The conditionally-averaged normal velocity and its streamwise
gradient associated with negative wall-pressure peak for
y/5= 0.05, 0.125, 0.25, 0.375, 0.5 and 0.625 .............................
Pseudo-instantaneous linear pressure source ((1%) flooded
contour maps and associated wall-pressure for three consecutive
time windows .................................................................
RMS profile of the linear pressure source calculated from
stochastically-estimated and measured velocity time series . ..
RMS profile of the linear weighted pressure source calculated
from stochastically-estimated and measured velocity time series
xiv
201
203
.. 205
206
Table B. 1
LIST OF TABLES
Table B. 1. Data acquisition settings ......................................................... 215
XV
1".
I .‘G‘V’lm
-.- .
.1"
..mw rV‘r'»
I . V‘ "V. 4.9
Abbrei'i:
AD
AC
CCD
FFP
FPaC L
HPF
H117
l-Plate
1.4
l d
LSE
LVDT
ML
M‘Piale
PDF
Ply
PM
031:
Abbreviation
A/D
AC
CCD
DC
FFP
FFT
F PaCL
HPF
HW
I-Plate
kx-f
1/d
LSE
LVDT
ML
M-Plate
MS
PDF
PIV
PWT
QSE
NOMENCLATURE
Analog to Digital
Alternating Current
Charge Coupled Device
Direct Current
Forward Flow Probability
Fast Fourier Transform
Flow Physics and Control Laboratory
High Pass Filter
Hotwire
Instrument Plate
Wavenumber-Frequency
length-to-diameter ratio of hotwire
Linear Stochastic Estimation
Linear Variable Differential Transformer
Multi-point Linear Stochastic Estimation
Middle Plate
Mean-shear pressure-source term
Probability Density Function
Particle Image Velocimetry
Plane Wave Tube
Quadratic Stochastic Estimation
xvi
1015
SAR
51.
SP1.
Q
S
N.
I'll.
A: in: & :5
ALE-qua 0;
.eg,
B.
RMS
SL
SPL
SQ
TM
TT
TTL
Au, lin & Av, lin
Au, quad & Av, quad
Bu&Bv
Cp, max & Cp, min
Root Mean Square
Signal Attenuation Ratio
Single-point Linear Stochastic Estimation
Sound Pressure Level
Single-point Quadratic Stochastic Estimation
Turbulence-Manipulation section of wind tunnel
Turbulence-Turbulence pressure-source term
Transistor—Transistor Logic
Half of the plane-wave-tube width
Stochastic Estimation coefficient of the linear term
LSE coefficients of the streamwise and normal velocities,
respectively
QSE coefficients of the linear term in the estimation of the
streamwise and normal velocities, respectively
Stochastic Estimation coefficient of the quadratic term
QSE coefficient of the quadratic terms in the estimation of the
streamwise and normal velocities, respectively
High-pass-filter capacitance
Speed of sound
Mean-pressure coefficient
Maximum and minimum mean-wall-pressure coefficient values,
respectively
xvii
0&C.
K&Km
kick:
Dt
K&Km
k1&k2
patm
Fluctuating wall-pressure coefficient
Cross-correlation coefficient of the wall-pressure
Cross-correlation coefficient between the measured and
stochastically—estimated streamwise and normal velocities,
respectively
Half of the wind tunnel test-section height
Voltage output of hotwires
Voltage output of Linear Variable Differential Transformer
frequency
Yaw-response functions of hotwire #1 and #2, respectively, of the
X-wire sensor
Half of the total fence height
Fence height above the splitter plate
Frequency-dependent and mean microphone sensitivity,
respectively
Yaw-response constants of hotwire #1 and #2, respectively, of the
X-wire sensor
Streamwise wavenumber
Correlation time delay expressed in number of data points
Total number of data points in a time series
Fluctuating wall pressure
Root mean square of the wall-pressure fluctuations
Atmospheric pressure
xviii
pi; & P?
. ...Ivfl‘hbwi: "E,
RJE&R
w§&(
po & pi
pr
ps
R2
Re
Reference and variable microphone locations, respectively
Reference mean static-pressure upstream of the model
mean static-pressure at the pressure taps
Hi gh-pass-filter resistance
Reynolds number based on the step height
Reynolds number based on local momentum thickness
Acoustic RMS pressure ratio
Auto-correlation coefficient of the wall-pressure
Position vector and its magnitude, respectively, of the pressure-
observation point
Temporal cross-correlation coefficient between the streamwise-
and normal-velocity, respectively, and wall-pressure
Spatial cross-correlation coefficient at zero time shift between the
streamwise- and normal—velocity, respectively, and wall-pressure
Time
Time reference of the stochastically-estimated velocity vector field
Total flow velocity for X-wire calibration
Mean streamwise velocity
Conditionally-averaged mean-removed streamwise and normal
velocity, respectively
Wall-scaled mean streamwise velocity, u/ uI
Stochastically-estimated velocity
xix
UyVS
um"
f I
uSaVs
Uel & U62
uref
urrns
X,y
x0, yo and z()
Stochastically—estimated streamwise and normal velocity,
respectively
Freestrearn velocity
Convection velocity
Effective velocity of hotwire
Effective velocities of hotwire #1 and #2, respectively, of the X-
wire sensor
Mean streamwise velocity at the location of the Reynolds stress
peak
Root mean square of streamwise-velocity fluctuations
Voltage output of microphones
Root mean square of normal-velocity fluctuations
Streamwise and wall-normal coordinates, respectively
Coordinates of the pressure-observation point
Mean reattachment length
Coordinates of the pressure source
Wall-scaled y location, yut/v
y location of the Reynolds stress peak
BOX/Um)
Yaw-angle offset of hotwire
Dimensionless wall-pressure spectrum, (Dpvpr/(l/Z p U002)2
Boundary-layer thickness
Ratio between outer- and inner-boundary-layer scaling, SUI/V
XX
It
“P1 ’ “P2 ,... uPro
$1.11.], q)VV & q)UV
VP: ’ V132 "" V1310
Ko&7ti
<(DZ, s>
Wall-pressure spectrum
Cross-spectrum of the wall-pressure with the streamwise and
normal velocity, respectively
Cross-spectra between the streamwise velocity and wall-pressure
at the ten downstream-most microphone locations
Streamwise- and normal-velocity auto— and cross—spectra,
respectively
Cross-spectra between normal velocity and wall-pressure at the ten
downstream-most microphone locations
Characteristic length of the dominant eddies
Acoustic wavelength
The characteristics size of the pressure signature associated with
the dominant flow structures in the outer-shear and sub-boundary
layers, respectively
Kinematic viscosity of air
Rotational-traverse yaw angle
Density of air
Time delay
Time-delay of microphones
Time—shift of cross-correlation peaks
Angular frequency, 21rf
Spanwise vorticity of the conditionally averaged velocity field
xxi
Downstream distance from reattachment point, normalized by X,,
(x-x,yxr
xxii
1.I.\'T
(A
large vi:
excitatio
excitatio
surfaee's
and und:
domains
pressure
turbulen:
generatic
understa:
T
means (‘
remath:
Simultan
purpose .
IO ark-am.
SiSnail”L
Single-p(.
I
Chara“?
01E” app
1. INTRODUCTION
Separating/reattaching flows contain very energetic structures, which generate
large wall-pressure fluctuations. These fluctuations are a direct representation of the
excitation forces produced by the turbulent flow on the underlying surface. If such
excitation takes place at frequencies and wavenumbers of one or more of the underlying
surface's resonant modes, unwanted vibrations and noise will be generated. Investigating
and understanding the wall-pressure-field characteristics in both the spatial and temporal
domains is important to predict and/or control such undesired effects. Moreover, wall-
pressure measurements can be used as a non-intrusive technique for capturing the
turbulent flow activities above the surface that are responsible for the wall-pressure
generation. In this context, the wall-pressure signature can also be used to gain better
understanding of the turbulence processes that occur in wall-bounded flows.
The present investigation examines the surface pressure spatiotemporally by
means of a wall-microphone array in the developing flow downstream of the
reattachment zone of the flow over a fence-with-splitter-plate. Furthermore,
simultaneous measurements of the velocity field and wall pressure are conducted for the
purpose of investigating their relationship. The simultaneous data are also used as a tool
to examine the flow structures associated with the generation of various wall-pressure
signatures using Linear and Quadratic Stochastic Estimation methods based on multi- and
single—point wall-pressure information.
The flow geometry of a fence-with-splitter-plate is illustrated in Figure 1.1. The
characteristic features of such a complex flow field may be described in terms of five
overlapping flow zones. In zone I, a freestream flow approaches the fence. At the tip of
'_:n. ..
."
.L‘.'~
Figure l
Smaiifcc
The "11ml
n
”Maxed
60.99
\ Freestream
U00 .- _;_\ ________ 7 _______ 1
T . 7 : Outer—shear layer .
:. . . 7 ,,,,,, 4|
: I ' ,‘ 2 ,_...’ —————— Sub-boundarylayer[ |
> . ..
1 I X
1' ..r:' i 1
Fence Splitter-plate : l
, J 11 ‘ 111 1 2 IV v
Figure 1.1. A schematic of an ideal two—dimensional flow over a splitter-plate-with-fence
the fence, the flow is forced to separate forming a free shear layer and recirculating flows,
zone H. Due to entrainment, the shear layer grows in the downstream direction till it
reattaches on the surface of the splitter-plate forming an energized reattachment zone
(zone III). At the reattachment point, a portion of the flow goes upstream into the
recirculation region, while the other goes with the downstream flow. As a result of the
impingement of the shear layer on the splitter-plate’s surface, the vortical structures that
have been energized in the free shear layer create high-pressure fluctuations. As these
structures travel downstream, they are exposed to a continually weakening mean-flow
shear and, hence, undergo a relaxation process, progressively loosing their energy in zone
IV. In parallel, the newly created strong mean shear near the wall produces energetic
small-scale turbulence similar to that found near the wall of turbulent boundary layers.
The "border" between this small-scale near—wall sustained turbulence, on one hand, and
the continually-decaying, yet still energetic, large-scale vortices that were "born" in the
separated shear layer upstream splits the flow in zone IV into two main regions: sub-
boundar
through
direeuot
sufiieier
develop.
zone IV
1.1. Lit
reattaeh:
understu
boundar;
PIOV’ideg
SIUdieS I
1
of a blur'
boundary layer and outer-shear layer, respectively. This terminology will be used
throughout this document to designate these two different regions that "stack" in the y
direction to form a non-equilibrium, or developing, boundary layer region. Given
sufficient streamwise length, this boundary layer develops into an equilibrium, or fully-
developed, turbulent boundary layer in zone V. The focus of the present measurements is
zone IV.
1.1. Literature Review
Because the present study investigates a flow field that is transitioning from a
reattaching free-shear-layer state to an equilibrium-boundary—layer one, reviewing and
understanding the main physics of typical separating/reattaching as well as equilibrium-
boundary-layer flows are important. These flows have been extensively investigated in
the literature, and in order to review them comprehensively the discussion will be
prohibitively lengthy. Therefore, given the specific scope of the present work, the review
provided herein focuses mainly on investigations that address the relationship between
the velocity-field and associated wall-pressure signature. In addition, a few selected
studies dealing with fundamental flow physics pertaining to both flows are also
considered.
1.1.1. Separating, Reattaching and Developing Flows
Many researchers have investigated separating/reattaching flows, in order to
characterize the main flow features in the recirculating-flow and reattaching-free-shear-
layer zones. Cherry et al. (1984) conducted simultaneous wall-pressure and velocity
measurements in addition to smoke flow visualization in the separating/reattaching flow
of a blunt-face splitter plate. Pressure was measured using two pressure transducers,
while V
dominar
attribute
zone be
simultar
reattaeh.
with the
3550(13IC
between
shear-la;
dimensi.
SUTJCIUIL
~—4
been 1'6;
1 1987) 5
backu'ar
flapping
IOCaIiOn;
employ.
lite Step
dEIachm' I
separallo
while velocity was measured using a conventional hotwire. Their data showed a
dominant low frequency in the wall-pressure fluctuations near separation, which they
attributed to the successive processes of expansion and contraction of the recirculation
zone beneath the free shear layer, or the so called flapping motion. Furthermore, their
simultaneous smoke flow visualization and pressure measurements downstream of the
reattachment point revealed that negative wall—pressure-peak events were synchronized
with the passage of what appeared to be the cores of vortical structures while those
associated with positive pressure peaks occurred in the vicinity of inward-flow regions
between the vortices. In addition, their Spanwise measurements showed that the free-
shear-layer structures became three dimensional soon after separation, but this three
dimensionality did not seem to be influenced by the reattachment process as the
structures "impinged"/interacted with the wall.
The flapping motion of the free shear layer in separating/reattaching flows has
been reported by several researchers in addition to Cherry et al. (1984). Driver et al.
(1987) studied the flapping motion observed in a reattaching free shear layer of a
backward-facing step in order to determine the frequency and spatial extent of the
flapping motion. To this end, they used thermal-tuft measurements at different x
locations downstream of the step, in addition to velocity and wall-pressure measurements
employing a wall-flush pressure transducer mounted at 5.5 step heights downstream of
the step. Driver et al. explained the reasoning for the flapping motion as a disorder of the
shear layer that arose when a vortical structure escaped the reattachment zone. Such a
detachment process of a vortical structure reduced the engulfed reverse-flow by the
separation bubble causing it to collapse momentarily, which increased the angle of
'2 :rrr‘xrcgg-HH
impingt
thatal:
toexpm
be less
of relati
explana‘
zone IHl
moxed
Structure
nzecau
point.
imbalanr
reattach,
FCCII C 111;;
Caused 1:
(.2001) s
“1335 Ur e?
bellt'een
hiporhe,
0n the I.
c0000]
impingement and, consequently, the curvature of the shear layer. Therefore, they stated
that a larger streamwise pressure gradient was created at reattachment, forcing the bubble
to expand back again. Driver et al. also estimated the amplitude of the flapping motion to
be less than 20% of the shear layer thickness. However, they stated that the motion was
of relatively low energy and might be ignored. Lee and Sung (2002) suggested a similar
explanation of the free-shear-layer flapping motion. They stated that the recirculation
zone increased linearly in size as the shear-layer vortical structures grew in size and
moved closer to wall. The expansion of the recirculation zone continued till a vortical
structure left the free shear layer, after which the recirculation zone abruptly shrank in
size causing a sawtooth-like behavior of the instantaneous location of the reattachment
point.
Eaton and Johnston (1982), and later Kiya and Sasaki (1985), suggested the
imbalance between rates of flow entrainment by the shear layer and reinjection at the
reattachment point was the main source of the enlargement and shrinkage of the
recirculation bubble. Furthermore, Kiya and Sasaki stated that such an imbalance was
caused by the breakdown of the Spanwise vortices in the shear layer. Spazzini et al.
(2001) studied the unsteady behavior of a backward facing step flow using skin friction
measurements and flow visualization. Their results exhibited a strong correlation
between the growth and successive breakdown of the secondary re-circulation bubble (at
the comer of the step) and the flapping motion of the free shear layer. This led them to
hypothesize that the flapping motion was linked to the behavior of the secondary bubble.
On the other hand, Heenan and Morrison (1998) conducted experiments to passively
control the low-frequency, buffeting, wall-pressure fluctuation downstream of a
pm... -———
back“?
found t
with It?
in the s
Spatiote:
micropn
backward-facing step using a reattachment surface with different permeability. They
found that the flapping motion totally vanished when a permeable reattachment surface
with length extending from the step to 0.56 of the mean reattachment length, X,, was
employed. The permeable surface apparently inhibited the recirculating flow and
associated upstream convection of disturbances produced at reattachment.
Lee and Sung (2001) made laboratory measurements of wall-pressure fluctuations
in the separating/reattaching flow over a backward-facing step. They investigated the
spatiotemporal statistical properties of the wall-pressure fluctuations using a 32 electret-
microphone array in both the streamwise and Spanwise directions. Based on the
wavenumber-frequency spectra of their data, they suggested that the shear-layer vortical
structures were modulated by the flapping motion of the free shear layer and moved
downstream with a convection velocity of 60% of the freestream velocity. In a more
recent work, Lee and Sung (2002) introduced a new spatial filtering technique, which
they referred to as Spatial Box Filtering (SBF). This approach was basically utilized to
"isolate" the wall-pressure signature of certain streamwise wavenumbers. Using the
zeroth and second modes of the SBF, they could adequately isolate the fluctuations
corresponding to the flapping motion from those corresponding to the passage of the
vortical structures generated in the free shear layer. The data showed good agreement
between the wavelength of the second mode of the SBF and the vortical structures
streamwise spacing, which was approximately half of the mean reattachment length.
The separating/reattaching flow over a fence-with-splitter—plate model was first
studied by Arie and Rouse (1956) and further investigated by many authors. For example,
Castro and Haque (1987) reported detailed measurements within the
Fl
separz
motiv;
lamina
the str
the sta
layer 1
would
splitter
investig
X.) C a
separate
the Shea
being in
“‘35 3131
163113: ,5}
51315 (la
axis-mm
of Rem
[UFbUien-
initially
Spread r;
separating/reattaching flow of a fence-with—splitter-plate model. They stated that their
motivation for selecting this flow geometry was the attractiveness of having a thin
laminar boundary layer at the separation point. The thin boundary layer is created due to
the strong favorable pressure gradient towards the edge of the fence that is produced by
the stagnation flow on the upstream face of the fence. Thus, the effect of the boundary-
layer thickness at separation on the flow was negligible and the reattachment length
would be only affected by the blockage ratio, h/Dt (h being the fence height above the
splitter plate and Dr is the wind tunnel test-section’s half width. Smits (1982) extensively
investigated the relation between the blockage ratio and the mean reattachment length,
Xr). Castro and Haque’s measurements demonstrated that the turbulent structures of the
separated shear layer and the plane-mixing layer were quite different. They reported that
the shear-layer grth rate was neither linear nor equal to that of the plane-mixing-layer,
being initially rather higher but reducing gradually as reattachment was approached. This
was also associated with a continuous increase in turbulence energy all the way to
reattachment, followed by a relatively rapid fall thereafter.
Hussain and Zedan (1978) also investigated the effect of the initial boundary layer
state (laminar or turbulent) and Reynolds number on the flow characteristics of an
axisymmetric free shear layer. They could show that the flow features were independent
of Reynolds number, but dependent on whether the initial boundary layer is laminar or
turbulent. In particular, they observed that while the separating boundary-layer was
initially either laminar or turbulent, its momentum thickness showed independence of the
spread rates, similarity parameters, and evolution of the shear layer. In contrast, those
values
mrbule
pressur
pressur:
They cr
length (
scale d‘
i WhiCh \
k‘ i UpSirear
second
fOUnd It
SthIurt
(100715“-
the free
middle .
OSCIllalti.
Collllacj ~
boundar
[hey Ira
understa
One far
values showed significant dependency on whether the initial boundary layer was
turbulent or laminar.
Recently, Hudy (2001) and Hudy et al. (2003) compiled a database of wall-
pressure-array measurements for studying the spatiotemporal character of the surface
pressure within the separating/reattaching flow region in a fence-with-splitter-plate flow.
They could distinguish two regions, which are defined based on the mean reattachment
length of the separated shear layer. In the first region, from the fence to 0.25X,, large-
scale disturbances dominated the signature of the wall-pressure. These disturbances,
which were associated with the shear layer flapping, were found to convect in both the
upstream and the downstream directions with a convection velocity of 0.21Uw. In the
second region, which was located beyond 0.25X,, smaller time-scale structures were
found to be responsible for the generation of the wall-pressure fluctuations. These
structures corresponded to the free-shear-layer vortices and traveled with an average
downstream convection velocity of 0.57Um. Hudy et al. could also relate the flapping of
the free shear layer to an absolute instability zone, or self-sustained oscillator, near the
middle of the recirculation region. They suggested that due to such a self-sustained
oscillator, the separation bubble continuously underwent processes of expansion and
contraction leading to the free-shear-layer flapping motion.
As mentioned above, the flow downstream of reattachment is a non-equilibrium
boundary-layer flow, in which the shear-layer vortices undergo a relaxation process as
they travel downstream. Several authors have investigated this flow seeking better
understanding of how the flow relaxes from its shear-layer—like state to a boundary-layer
one far downstream. For example, Bradshaw and Wong (1972) used existing
expenn
their 0\
nature :
druid
categor
byarn
which C
the stro
unny
OVem'h
m anot
Ckuufi;
t110 ove
free 3}]:
the latte
such as
shearlai
l
the ma
Ci}:
hޤu4%
anequfli
me“ah,
Li“- B:
experimental results on the low-speed flow downstream of steps and fences in addition to
their own measurements downstream of a backward-facing step to demonstrate the flow
nature in the separation/reattachment and relaxation (developing) zones. They defined
three degrees of perturbations; i.e., weak, strong and overwhelming perturbation to
categorize different separating/reattaching flows. The weak perturbation is that caused
by a minor disturbance such as a mild change in pressure gradient or wall roughness,
which does not significantly change the flow velocity or length scale. On the other hand,
the strong perturbation resembles the flow of a boundary layer over a very small notch or
cavity that significantly alters the turbulent structures in the flow, while the
overwhelming perturbation is one that changes, all together, the flow from one "species"
to another; e.g., a shear layer changing to a boundary layer. According to the
classification of Bradshaw and Wong, the typical backward-facing-step flow involves
two overwhelming perturbations: one when the flow changes from a boundary layer to a
free shear layer at separation, and the other, at reattachment, when the flow switches from
the latter to the former state back again. However in the case of a thin boundary layer,
such as in the fence flow, the first perturbation may be ignored assuming that the free
shear layer starts developing at the separation point.
Bradshaw and Wong (1972) conducted their experiments employing a 0.13 step-
height-thin laminar boundary layer at separation. They observed a marked deviation in
the mean-velocity profile of the boundary layer downstream of reattachment from that of
an equilibrium boundary layer. This deviation manifested itself as an overshoot, closer to
the wall, and undershoot, away from the wall, of the profile relative to the universal log-
law. Bradshaw and Wong suggested that the main reason for the deviation was the
displ
u all
the It
the st
irmer
obser
up u)
produc
measui
thin in
I“Dillard.
disproportionality between the heights of the dominant turbulent structures above the
wall and their length scales. This disprOportionality (i.e., length scale ~ y), upon which.
the log-law is based, was attributed to the practically constant length scale (of the order of
the step height) of the free shear layer vortical structures at reattachment, except in the
inner layer very near the wall. Such a deviation from the equilibrium boundary layer was
observed as far as 52 step heights (approximately 8.6X,) downstream of the step.
Chandrsuda and Bradshaw (1981) made hotwire measurements, which extended
up to two reattachment lengths downstream of the step, of the second- and third-order
products of the turbulent velocity components behind a backward-facing step. Their
measurements showed that the free shear layer in the separation/reattachment zone with a
thin initial laminar layer was not sensitive to its initial conditions. Unlike Castro and
Haque (1987), Chandrsuda and Bradshaw reported that the free shear layer was not
greatly different from a plane mixing layer with a uniform external stream. Moreover,
they observed that the shear layer underwent a rapid change near the reattachment zone,
which was related to the confinement effect of the wall on the vortical structures. This
wall effect also led to an attenuation of the normal component of the flow velocity and
the transport of turbulence energy towards the wall.
F arabee (1986) and Farabee and Casarella(l986) conducted measurements of the
wall-pressure field underneath two separated/reattached boundary-layer flows: over a
forward- and backward-facing step. They found that the process of separation and
reattachment of a turbulent boundary layer produces very large wall-pressure fluctuations.
In the backward-facing step, the wall-pressure fluctuations were higher than those of the
forward-facing step and equilibrium boundary layer by factors of five and ten,
10
" .. ‘3 Irwin-aw!
afllic
hnghu
meob:
bound;
layer 0'
layer d1
bhuu fl
1
that thc‘
also ob
lhan the
that. n ,
dex'elori
the H hoi
i
FECOV‘er)
Similar 5
vortices
I:arther '
reia\allc
Ref-hold
respectively. Although the energetic flow structures responsible for the generation of the
wall-pressure decayed as they were convected downstream by the mean flow, they were
still identifiable in Farabee and Casarella’s measurements at locations up to 72 step
heights downstream of the backward-facing step. The latter observation, which sustains
the observations of Bradshaw and Wong (1972), indicates that the upstream history of a
boundary layer could be significant on the aero/hydro-acoustic features of the boundary
layer over substantial streamwise distances.
Castro and Epik (1998) also conducted measurements in a developing boundary
layer downstream of a separating/reattaching flow downstream of the leading edge of a
blunt flat plate (the same model used in Cherry et al. 1984). Castro and Epik reported
that the log-law did not exist immediately beyond the mean reattachment point. They
also observed a slow development process in the inner sublayer; although it was faster
than the development of the outer part of the flow. Nevertheless, Castro and Epik argued
that, notwithstanding the faster rate of development of the inner region, the rate of
development of the outer part of the flow determined the overall rate of development of
the whole flow.
Recently, Song and Eaton (2002) investigated the separation, reattachment and
recovery regions of a boundary layer flow over a curved ramp. Their flow exhibited
similar structural features to those of the backward-facing step, including the roll-up of
vortices downstream of separation followed by their partial distortion at reattachment.
Farther downstream, a non-equilibrium boundary layer existed and underwent a similar
relaxation process to that discussed above. Song and Eaton (2002) also investigated
Reynolds number effects, varying the Reynolds number by changing the air density and
11
freestr.
They al
the out
ICCOVE‘T
110w C
RudenT
model.
pulsed-
geomet
dou‘nst
dongai
separati
mec‘iSUrt
Relholt
l
dismbu.
e‘lLIilibr
Farabee
311d Bra
EVidel’lc:
heights 1
found In
freestream velocity. They found evidence of the flapping motion of the free shear layer.
They also reported that the remnants of the separated flow vortical structures dominated
the outer layer in the recovery region. Moreover, this outer portion of the flow seemed to
recover to the equilibrium state at much slower rate than that of the near-wall flow.
On the other hand, a few investigations studied the developing boundary layer
flow downstream of the separating/reattaching flow over a fence-with-splitter—plate.
Ruderich and Femholz (1986) investigated the flow over a fence—with-splitter—plate
model, and carried out mean and fluctuating velocity measurements using hotwire and
pulsed-wire anemometry. Similar to Castro and Haque (1987), they selected this flow
geometry because of the negligible effect of the upstream boundary layer on the flow
downstream of separation. Another attractive feature of this flow field was the
elongation of the recirculation region due to the steep angle of the velocity vector at
separation, which improved the measurement resolution. Ruderich and F ernholz
measurements were used to estimate the mean-velocity, Reynolds shear stress and
Reynolds normal stress distributions in the wall-normal direction. The mean-velocity
distribution of their flow downstream of reattachment exhibited a deviation from the
equilibrium boundary layer similar to that reported by Bradshaw and Wong (1972) and
Farabee (1986) in the back-step flow. Also, consistent with Farabee and Casarella (1986)
and Bradshaw and Wong (1972), Ruderich and Femholz (1986) could not observe any
evidence of the equilibrium boundary layer even at the end of the splitter plate (68 step
heights). Interestingly, Ruderich and Femholz's study is one of the very few studies that
found no evidence of the free-shear-layer flapping motion.
12
surface
problet
physic
betyyee
forces
probler
related
by taki:
leg. SC
Where i
apply:
imPOI‘Ia
FIESSUn
“here
imbue?
r espeeti'
the Gov.
1M 3 ). o
1.1.2. Characteristics of the Turbulent Wall Pressure
In the past few decades, there have been many investigations of the fluctuating
surface-pressure field in turbulent flows. Although the solution of the mathematical
problem of turbulence (the closure problem) depends on the ability to understand the
physics of, and successfully model, the Reynolds stress terms, understanding the relation
between the pressure and velocity fields is important for predicting the unsteady flow
forces acting on a surface for devising solutions for flow-induced noise and vibration
problems and controlling the flow field. Basically, turbulent pressure fluctuations are
related to the velocity field of incompressible flows through Poisson’s equation, obtained
by taking the divergence of the Navier-Stokes equations. Poisson’s equation is given by
(e.g., see Willmarth 1975 and Kim 1989):
1
—V2p=-—ui‘j up, (1.1)
p
where p is the pressure, p is the fluid density and u, is the total velocity vector. By
applying the Reynolds decomposition for a two-dimensional mean flow with only one
important mean-shear component, and considering only the turbulent component of the
pressure, equation (1 .1) simplifies to:
(1.2)
where u is the mean streamwise velocity, the prime denotes the mean-removed, or
turbulent, quantities, and x and y are the streamwise and wall-normal coordinates,
respectively. Equation (1.2) shows how the wall-pressure fluctuations are a function of
the flow "sources" in the turbulent flow field. The source terms consist of the mean shear
(MS), or linear term, and turbulence-turbulence (TT), or nonlinear terms (e.g., Chang et
13
111.
the
char
using
.\
“mu—“k
H0“ 6
‘l‘__ v.
the 110
been I
teehnic
researc
their re.
method
3133110113
lhe OCCl
l'aiUeS.
al. 1999). Because the MS term represents the interaction between the mean shear and
the fluctuating-velocity gradients (first term on the right hand side of equation 1.2), it
changes instantaneously in response to changes in mean-flow conditions and is called the
rapid term. On the other hand, because the TT term represents the interactions between
the fluctuating-velocity gradients, it is affected only once the turbulence has had a chance
to adjust to the new mean-flow conditions, and therefore this it is called the slow term
(Kim 1989).
Ideally, it is desired to measure the pressure fluctuations inside a turbulent flow
using a non-intrusive technique to avoid introducing any error into the measurements.
However, measuring the pressure fluctuations in a turbulent flow without interfering with
the flow field is not possible to date. Therefore, pressure fluctuations measurements have
been limited to the wall in the case of wall-bounded flows. Such a measurement
technique is non-intrusive, which made wall-pressure measurements attractive to many
researchers who aim to gain better understanding of both the turbulent structures and
their relationship to the wall—pressure signature.
Typically turbulent wall-pressure fluctuations are characterized using statistical
methods, either conditional or long-time averaged. The former are used to study the
spatiotemporal features of the wall-pressure using its ensemble-averaged data relative to
the occurrence of strong positive or negative pressure peaks, while the latter include RMS
values, spectra, probability density function (PDF), and higher-order moments of the
wall-pressure. Below, an account is given of the wall-pressure characteristics in
separating/reattaching, developing-boundary-layer and equilibrium-boundary-layer flows
based on existing literature.
14
1. Sep‘.
(Farah
haye r;
shear-l
SlrUClt.‘
0n the
shear-)
dOWnSt
b)? ChL
rePorte
Smittu:
110W er
slou- re
I. Separating/Reattaching Flows
As mentioned earlier in 1.1.1, several investigations have studied the wall-
pressure characteristics in separating/reattaching (e.g., Cherry et al. 1984, Driver er al.
1987, Lee and Sung 2001 & 2002, and Hudy et al. 2003) and developing-boundary-layer
(Farabee 1986, and Farabee and Casarella 1986) flows. The former set of investigations
have reported that the wall-pressure beneath the recirculation bubble is dominated by the
shear-layer flapping motion very close to separation. Farther downstream, the vortical
structures in the free shear layer grow in size and approach the surface till they impinge
on the wall at reattachment. As a result, the wall-pressure signature is dominated by the
shear-layer vortices, and the wall-pressure fluctuation level increases with increasing
downstream distance till it reaches a peak slightly upstream of reattachment, as discussed
by Cherry et al. (1984) and Hudy et al. (2003). Similarly, Farabee and Casarella (1986)
reported that the wall-pressure signature was dominated by the free-shear—layer vortical
structures and their level remained higher than those of the equilibrium boundary layer
flow even at 72 step heights (z12xr) downstream of the back step. They attributed that to
slow relaxation of the vortical structures.
The physical relation between the large-scale vortices in the reattaching free shear
layer and farther downstream in the non-equilibrium boundary layer and the wall-
pressure signature has been an interest for many researchers (e. g. Cherry et al. 1984, Lee
and Sung 2001, Lee and Sung 2002, Kiya and Sasaki 1985). Their investigations have
shown that negative peaks in wall-pressure fluctuations are associated with the passage of
large-scale vortex cores. On the other hand, positive peaks were found to occur beneath
the downward inrush of freestream fluid inbetween the vortical structures.
15
:
. I! “’7‘
ofPoh
doyxns‘
Hmanq
fluctua
(1mm 5
H.Ful
boundd
l
at the \t
l
Bun 1g
Offlou'
‘ 1
hfihfip
that the:
53f“ i
Which
A few studies have investigated the spectral characteristics of wall-pressure in
separating/reattaching flows, revealing that as the flow approaches reattachment its wall-
pressure signature becomes more dominated by the free-shear-layer vortical structures
(e.g., Lee and Sung 2001, Hudy et al. 2003, and Farabee and Casarella 1986).
Additionally, Lee and Sung (2001) observed a build up of a frequency range in the wall-
pressure spectra with a slope of —7/3 as the flow progressed towards and downstream of
reattachment. This agrees with the study of George et al. (1984), who reported —11/3 and
—7/3 spectral slopes associated with the MS and TT terms, respectively, in their solution
of Poisson’s equation in a free shear flow. Lee and Sung (2001) suggested that farther
downstream of reattachment, the TT interaction term became more prominent due to the
enhanced small-scale activity. Moreover, they concluded that the wall-pressure
fluctuations are largely attributable to the free-shear-layer vortical structures even
downstream of reattachment.
II. Fully-Developed Turbulent Boundary Layers
Most conceptual models of near-wall coherent structures in turbulent wall-
bounded flows are based on the theme of a horseshoe (hairpin) vortex that forms locally
at the wall, with the legs of the horseshoe trailing upstream of the arch (e.g., Thomas and
Bull 1983 and Lueptow 1997). The horseshoe vortex structure is associated with ejection
of slow speed fluid away from the wall (typically referred to as Q2 event) and inrush of
high-speed fluid towards the wall (or Q4 event). Thomas and Bull (1983) found evidence
that the wall-pressure positive peak was associated with a sudden step-like rise pattern,
even in the outer part of the boundary layer, in the streamwise velocity fluctuations,
which reflected a large-scale strong sweep event. On the other hand, Farabee and
16
50 UFC t‘
lilal Ft
bound
Casarella (1991) established the scaling law for the low-, mid-, and high-frequency
regions of the wall-pressure spectra. Their results showed the existence of two
"wavenumber" groups: a high wavenumber group that was associated with turbulent
sources in the logarithmic region of the boundary layer, and a low wavenumber group
that represented the large-scale turbulence contribution from the outer region of 'the
boundary layer.
Johansson et al. (1987) used conditional-average data in turbulent boundary
layers to show that the buffer region structures are responsible for the generation of large
positive wall-pressure peaks, which indicated a link between wall-pressure peaks and the
turbulence-producing mechanisms. Johansson et al.’s data indicated that the pressure-
peak amplitude was found to scale linearly with the conditional-velocity amplitude
indicating the dominance of the MS source term in generating the wall pressure.
However, Kim (1989) and Chang et al. (1999) found that in htrbulent channel flow the
slow and the rapid pressure fluctuations are of equal importance very near the wall.
Over several decades, many researchers have used scaling arguments to show that
the wall-pressure spectrum ((Dpvpv) for a certain frequency range should obey a power-law
type behavior; i.e., (Dp'pr ~ 00", where n is a function of the frequency range of the power
spectrum. Bradshaw (1967), Panton and Linbarger (1974), Farabee and Casarella (1991),
Gravante et al. (1998), and Chang et al. (1999) have shown that the wall-pressure power
spectrum of a turbulent boundary layer should exhibit a fall-off rate of a)“ in the middle
range of frequencies. However, Gravante et al. (1998) showed that the frequency extent
of this region decreased or even disappeared with decreasing Reynolds number, which
17
0
’-
o
o
. -
..
1:.
-..
-"
. ..
..
..
O
.\.i
so"
I," 1
1?,
:‘HD
'7':
I 4
”0‘
1‘- -
'ki
.‘ m
.1
."
Jit
’Lr
'".
‘
i‘
a
in "n feat-“aa-
ryimi-‘u in? an .2.
lug-3:: "111v: ..-
manure;
..h 31'1“}? 13"
J1"
LM'td 'I
'Azi’vs
.i'.
353.
J:
t: ‘o
. 17,33);
' 3" WM“ "‘P’ “fi—Zli
.1473. ‘
"3
MW ‘
pressu
could-
also it
this it
freque
b€h3\ j
\‘ari ab
literat;
11999;
1.1.3.5
here a
SIOChg
COndh:
“W St
3h Le“
Adfiar
and 1h.
agreed with Kim (1989), who reported the nonexistence of the —1 region in their low
Reynolds number channel flow.
On the other hand, Monin and Yaglom (1975) theorized a spectral fall-off rate of
(0’73 in the inertial subrange of the wall-pressure power spectrum, corresponding to the
pressure signature of locally isotropic turbulence. The results of Gravante et al. (1998)
could not reveal a substantial —7/3 spectral range. Earlier, Farabee and Casarella (1991)
also noted the absence of the -7/3 range from their spectral measurements and attributed
this to the spatial-resolution limitation of their pressure sensors. Finally, in the high-
frequency region of the turbulent wall-pressure spectrum, Blake (1986) theorized an 00‘5
behavior that was independent of Reynolds number when scaled using wall, or viscous,
variables. The decay rate of 03—5 has been also reported in several investigations in the
literature of the turbulent boundary layer; e. g., Gravante et al. ( 1998) and Chang et al.
(1999).
1.1.3. Stochastic Estimation
As will be outlined in the objectives of this work, stochastic estimation is used
here as a tool to estimate the flow velocity field from its wall-pressure signature.
Stochastic estimation was first used by Adrian (1977 & 1979) to characterize the
conditional eddies of isotropic turbulence. He examined the existence of the conditional
flow structures by computing the estimated velocity u,(x+r, t) from a known velocity u(x,
t); i.e., , which was referred to as “conditional eddies”. Basically,
Adrian proposed that the estimated velocity (u,) could be expanded in a Taylor series of u,
and then truncated at a certain order. When only the first term is included in the series,
18
:LIJL-
the es
lClTllS
t0 e51:
panic
isotro;
qualit.
W35 5
and .\T
a horr.
tensor
genera
Were
the estimation is referred to as Linear Stochastic Estimation (LSE) while if the first two
terms are contained in the expansion, Quadratic Stochastic Estimation (QSE) is obtained.
In the 19805 and 19905 many researchers employed stochastic estimation as a tool
to estimate the turbulent velocity field from known velocity and/or wall-shear stress at a
particular point in space. Tung and Adrian (1980) examined the conditional eddies in
isotropic turbulence by estimating the velocity field around a point where the velocity
was specified. They concluded that LSE was a satisfactory method for studying the
qualitative large-scale features of the conditional eddy field. On the other hand, Adrian
and Moin (1988) could characterize quantitatively the large—scale organized structures of
a homogenous turbulent shear flow using LSE based on the velocity and the deformation
tensor at an arbitrary point. They applied the linear estimate to the turbulent field data
generated from direct numerical simulation, DNS. They found that using the Reynolds
shear stress as the event in the estimation, the largest contribution to the Reynolds stress
were associated with second- and fourth-quadrant events.
Guezennec (1989) applied the stochastic estimation technique to a fully-
developed turbulent boundary layer. His study depicted the advantage of the technique in
obtaining conditional averages from the unconditional statistics to estimate the
spatiotemporal characteristics of the second- and fourth-quadrant events in the boundary
layer. He also extended the usage of the stochastic estimation technique, based on a wall-
shear-stress condition, to single- and two-point conditional averages. Although, the
results showed that the difference between the linear and quadratic estimation was
minimal, it appeared that the quadratic estimation included additional information.
19
in the
comp?
as im:
preser
time-t
that t?
proyig~
appro.
ICITTIS.
recons
bound.
the o:
(LSE)
identi l;
GlaUSq
Guezennec also found that using higher estimation orders resulted in higher noise
in the estimation, which indicated that the quadratic estimation represented a “good
compromise” between the convergence of the series expansion and practical limitations,
as indicated earlier by Tung and Adrian (1980). On the other hand, Brereton (1992)
presented a procedure to assess the accuracy in the stochastic estimation model of the
time-delayed conditional-averaged velocity in a turbulent boundary layer. He showed
that the inclusion of higher-order terms in the estimation biases the stochastic model to
provide better representation of rarer events. Therefore, Brereton introduced a new
approach, which was based on tailoring the estimation model to include negative-order
terms, seeking better accuracy in estimating more-frequent events.
In addition, Choi and Guezennec (1990) used stochastic estimation as a tool to
reconstruct the conditional structure and to examine their asymmetry in a turbulent
boundary layer for various levels of Spanwise velocity perturbations near the wall. On
the other hand, Bonnet et al. (1994) introduced a linear stochastic estimation
(LSE)/Proper Orthogonal Decomposition (POD) complementary technique for
identifying structures in an axisymmetric jet and a 2-D mixing layer. Furthermore,
Glauser et al. (1999) used the LSE/POD as a low-dimensional model to estimate the
large—scale unsteadiness of the reattachment region in an axisymmetric sudden expansion.
Estimating the flow field features using their wall-pressure signatures has been
also sought by a few investigators. Naguib et al. (2001) were the first to seek to
understand the interrelation between the flow field structures and the turbulent wall-
pressure through stochastic estimation. They assessed the ability of stochastic estimation
to capture the conditionally averaged flow field associated with negative and positive
20
wall-pressure peaks by conducting simultaneous measurements of the wall-pressure at a
single-point and velocity in a turbulent boundary layer. The comparison between the
conditional and stochastic estimation (linear and quadratic) results showed that the linear
stochastic estimation based on the wall-pressure did not converge to the conditional
average, and that it was necessary to include the quadratic term. The latter was
inconsistent with the findings of earlier investigations employing stochastic estimation
e.g., Adrian et al. (1987) and Guezennec et al. (1987). However, the discrepancy was
attributed by Naguib et al. (2001) to the use of the wall-pressure as the estimation
condition. Specifically, Naguib et al. demonstrated that the need for the inclusion of the
quadratic term in the estimation was attributed to the effect of the TT pressure source
term.
Recently, Murray and Ukeiley (2003) have demonstrated the potential for multi-
point stochastic estimation, based on wall pressure, to accurately capture the
instantaneous flow field in a resonating cavity flow. The success of the estimation in the
cavity flow is believed to be due to the highly organized nature of this kind of flow and
the well-defined phase relationship between the pressure and the vortical structures in the
shear layer, which dominate the flow. In particular, feedback of the pressure disturbances
due to periodic vortex impingement on the downstream lip of the cavity to the upstream
lip sets the entire flow field in a highly organized resonant state. In a latter work, Murray
and Ukeiley (2004) also employed the same estimation procedure in a non-resonating
cavity flow. Although no direct comparison with the instantaneous flow field was
provided in this case, the authors demonstrated that the estimation captured the essential
features of the instantaneous flow structure.
21
ilk'r—SLF 1”“ ‘; " on. -«- .~ru
. - ..
estin
Thes
teehr
flow
separa‘
l 1986 '1
Separa
Only 11
the \'e:
and \‘e
Two other example studies that employ surface measurements for flow field
estimation are those by Taylor and Glauser (2002) and Schmit and Glauser (2004).
These studies employ a combined LSE/Proper Orthogonal Decomposition (POD)
technique for estimating the flow field above the wall. In the first work, wall-pressure-
array measurements were successfully employed to capture low-frequency fluctuations of
the instantaneous velocity field associated with the flow over a backward-facing ramp.
On the other hand, Schmit and Glauser (2004) employed an array of dynamic strain
gauges mounted on the wing of a Micro Air Vehicle to estimate the instantaneous flow
field in the immediate wake behind the wing.
1.2. Motivation
This work was motivated by the interest to understand the relation between the
flow field and wall-pressure in the non-equilibrium flow downstream of a
separating/reattaching flow. In particular, with the exception of Farabee and Casarella
(1986), no other study of the wall-pressure beneath the flow downstream of a
separating/reattaching flow was found. Furthermore, Farabee and Casarella employed
only two-point measurements in their work. In the current investigation, simultaneous
measurements of the wall-pressure field p'(x,t), using a sixteen-microphone array, and
the velocity field, employing an X-hotwire sensor, will be undertaken. The coupled array
and velocity measurements will allow:
1. Determination of the space-time characteristics, including the fi'equency-
wavenumber spectrum, of the wall pressure field. The documentation and
understanding of such characteristics is important to flow-induced noise and
vibration problems since it is the combined spatial and temporal characteristics
22
1.3. C
2.
SJ
that determine if the wall-pressure field excitation of the wall will lead to
substantial vibration and noise generation.
Multi-point Linear Stochastic Estimation (LSE) of the velocity field based on wall
measurements. The estimation is motivated by examining the flow structure and
wall-pressure sources associated with various types of wall-pressure signature.
Ultimately, such examination may aid in developing structure-based models of
the flow pressure sources, allowing simplified computations of the wall-pressure
field in engineering calculations of flow-induced noise and vibration problems
downstream of appendages and other separated flows in engineering devices.
1.3. Objectives
1.
The specific objectives of this work may be listed as follows:
Construct a 16-point fluctuating-wall-pressure-sensor array and develop an
appropriate calibration procedure. Integrate the sensor array into the fence-with-
splitter-plate setup of Hudy (2001) beneath the non-equilibrium boundary layer
downstream of reattachment.
Characterize the wall-pressure field through one- and two-point statistics as well
as the wavenumber-frequency spectrum.
Conduct simultaneous measurements of the wall-pressure and flow velocity using
the microphone array and X-hotwire sensor, respectively.
Examine the velocity-field characteristics and its relation to the wall-pressure
field utilizing various conventional and conditional statistical analyses.
Estimate the "pseudo-instantaneous" flow field associated with typical
instantaneous spatial wall-pressure patterns using multi-point wall-pressure based
23
' a ' 3."
descri
preset
bound
result:
FCCOII‘.
LSE. The estimation is expected to not only provide some sense of the variability
in the characteristics of the flow structures associated with typical wall-pressure
signatures, but also to lead to better understanding of the physical nature and
relative importance of the flow pressure sources. It is noted here that such a
detailed study into the nature and significance of the wall-pressure-generating
flow structures have not been conducted to date in the flow field considered here.
The remainder of this work consists of four main chapters. Chapter 2 provides a
description of the experimental set-up and methodology, Chapter 3 contains a
presentation of the wall—pressure results, while Chapter 4 includes a discussion of the
boundary layer characteristics, pressure-velocity correlations and stochastic estimation
results. Lastly, Chapter 5 highlights the main conclusions of the present work and
recommended future work.
24
F
[mi-.- -:9 z ‘ war
F." \I ‘J'. UHF-)1} "74
2. E2
in the
ofin:
hotyy
2. EXPERIMENTAL SET-UP AND METHODOLOGY
This chapter describes the experimental set-up, techniques and methodology used
in the present work. Four measurement techniques were used to investigate the flow field
of interest: static-pressure taps, flow visualization using tufts, microphones and single/X-
hotwire probes. Two different sets of data-acquisition hardware were used to acquire the
experimental data resulting from application of these various techniques. A full
description of the set-up, methodology and data-acquisition hardware is provided in
subsequent sections.
2.1. Experimental Set-up
The experimental set-up consisted of a fence-with-splitter-plate model, an open-
circuit wind-tunnel facility, the sensing instrumentation and data-acquisition systems.
Each of these constituents is described in details in this section.
2.1.1 Wind Tunnel Facility
The present experiments were conducted in a low-speed open-circuit indraft wind
tunnel1 facility in the F low Physics and Control Laboratory (FPaCL) at Michigan State
University. The wind tunnel consisted of four main sections (see Figure 2.1): the
turbulence-manipulation section, contraction section, test section and diffuser and blower
section. The overall dimensions of the wind tunnel were approximately 7.5 m long by
2.0 m high and 1.2 m wide.
The turbulence-manipulation (TM) section directed and conditioned the air flow
to the contraction section. The main purpose of this section is to produce uniform mean
velocity profile and to break up the large-scale eddies in the flow into smaller ones so that
' The wind tunnel was originally designed and built at Northwestern University by Steve Snarski and
Richard Lueptow, and was moved to MSU in 1999.
25
CFOSS
they can dissipate prior to entering the contraction section. The overall dimensions of the
TM unit were 1.15m x 1.15 m in cross section by 1.06 m in length. It consisted of 10
separable wooden frames. Each of the first five frames supported an l8-mesh aluminum
window screen while each of the last five ones supported 30—mesh stainless steel screen.
A 25,000-cell aluminum honeycomb element (6.0 mm cell size, 50.0 mm length) was
also housed in the filth section with the 18-mesh screen.
The contraction section accelerated the flow through a contraction ratio of 10.821.
The section had a 1.15 m-wide square inlet cross-section and a 0.35 m-wide square outlet
cross-section. The total length of the contraction section was about 1.5 m and it was built
of plywood and fiberglass with a Formica inner surface.
Turbulence Manipulation Section Contraction Test Diffuser/Blower
(Screens and Honeycombs) Section Section Section
0508 -
5 x 0.125 u
5 x 0.0875 ~ < ‘ 0966 ' 7 2'77 if- «048* 1.35 >
Figure 2.1. Schematic2 of the wind facility (dimensions in meters)
The test section was 2.77 m long with inlet cross-section matching that at the exit
of the contraction. The walls of the test section had a divergence angle of 0.13 degrees
with respect to the centerline. The angle was selected to compensate for the displacement
thickness of the boundary layer growing on the test-section walls in order to minimize the
2 Done by Chad Stimson, Michigan State University, East Lansing, MI 48824
26
streamwise pressure gradient for most of the wind tunnel operating velocity range. The
test section was made of plywood with a Formica inner surface supported by an
aluminum external frame. One of the test section walls contained three windows that
were made out of clear plexiglass. Each window was 0.9 m long, 0.3 m high and 12.5
mm thick. The windows were mounted such that their inner faces were flush with the
inner test-section wall.
The fourth, and the last, component was the diffuser/blower section. The section
was used to decelerate and turn the flow horizontally by 90 degrees towards the blower.
The section was made out of 14 gage steel exterior case with perforated 22 gage
galvanized steel interior flow surface and fiberglass acoustical wool inserted between the
two. The relatively high contraction ratio (10.8) combined with the 90-degree bend and
the perforated steel/acoustical lining were responsible for reducing the frequency and the
strength of the blower noise, which is a potential source of contamination of the wall-
pressure measurements. The high contraction ratio helped to reduce the blower rotation
speed for a given wind tunnel velocity, and hence reduce the frequency of the noise
generated due to the passage of the blower blades. On the other hand, the geometry and
sound-absorption lining of the diffuser helped to attenuate the upstream—propagating
noise resulting from the blower-blade passage.
The air flow in the tunnel is induced by an airfoil-vane centrifugal blower, type 27
SQA from Chicago Blower Corporation. The blower was driven by a dual belt drive and
Minarik 180 Volt, 3hp permanent magnet DC motor, model number 504-36-043A with
speed range of 30-1800 rpm. The speed control was provided by General Electrical
rheostat controlled AC motor-DC generator set with variable frequency controller of the
27
bloyy
indic
Pitot
’Q
.—
Id
1 see
symr
were
plate
1W0
regar
blower speed. The controller has a remote-operating keypad with a digital frequency
indicator. The frequency indication was calibrated against flow velocity measured by a
Pitot tube in the test section to facilitate setting of the freestream velocity.
2.1.2. Test Model
The fence-with-splitter-plate model of Hudy (2001) was used in the present work
(see Figure 2.2). The design of the model, which was constructed from aluminum, was
symmetric with respect to top and bottom. The total length and the width of the model
were 160h and 44h, respectively; where h = 8 mm is the fence height above the splitter
plate (or "step" height). In order to improve the two dimensionality of the mean flow,
two endplates were attached to the two sides of the splitter-plate (for more details
regarding the model design, see Hudy 2001).
Extensions ‘ Fence
‘ \
Flow \ - Instrument Plate
Middle Plate
Endplates \
T " \i\‘\ I I);
‘y ' - \‘a Tail Plate
Z/ \x
28 Pressure Taps
l
. E
16 Microphones a
(’1
c5
1
_ <0 f
_, 0:115
/
Figure 2.2. A schematic of the test model
28
rahc
expl.
tneas
than
pres;
Inods
Skelc
F1911?
ificthi
The total fence height (2H = 2h + splitter plate thickness) was 35 mm. This fence
height was selected to produce reasonable flow blockage while maintaining a long
reattachment length to "stretch" the flow field and facilitate spatial resolution of the
measurements in the x direction. Smits (1982) defined a blockage ratio (h/D.)
representing the ratio between the step height (h) and half of the test section height (D1).
He found that the blockage ratio resulted in flow acceleration around the fence, which
decreased the mean reattachment length as the blockage ratio increased. The blockage
ratio of the present model was 4.5% (see Figure 2.3 for visual illustration). As will be
explained later, the reattachment length (X,) of the present work was estimated using tuft
measurements, which produced an Xr value of 180 mm (22.5h). This is slightly shorter
than Xr of Hudy et al. (2003) because their blockage ratio (2%) was lower than the
present study.
For the current study, an array of 16 microphones and 56 pressure taps (28 on
each side of the model) were mounted in the splitter plate. The splitter plate consisted of
a half-inch skeleton sandwiched between 3.175 mm-thickness aluminum plates. These
plates included an instrument plate (I-plate) and a middle plate (M-plate) on each side.
The I-plate was 51h (406 mm) long while the M-plate was 76h (610 mm) long. The tail
plate (T-plate) was used to reduce plate-wake disturbances by gradually reducing the
model thickness to zero. The T-plate length was 32h (254 mm). The volume within the
skeleton was used for microphones insertion, as well as wires and tubes connections.
Figure 2.4 shows schematic drawing of the “top” view of the model. A corresponding
picture of the model installed in the wind tunnel is shown in Figure 2.5. Note that,
29
. ' V ‘ 'v' " a: ‘Vn77'7'fiirm'.
Flt
PICS I
Figure 2.3. A frontal picture of the model showing the blockage effect of the fence
Support Test section
Members top wall Endplates
.................. . \ ._._.._ .\
........ \ ,\\ \.\..
q 7 1600 mm
‘ 406.4 mm __
20 mm _ _‘
28 Pressure taps
Flow 4 ~ ' a»,
, I ‘ 300 mm > 16 Microphones
Fence
I 1
Instrument plate Middle Plate
Figure 2.4. A schematic of the instrument and middle plates depicting the locations of the
pressure taps and microphones
30
images in this thesis are presented in color. The 28 pressure taps and six microphones
were installed in the l-plate while the middle plate contained the rest (10) of the
microphones. The pressure taps were 9.5 mm center-to-center apart while the
microphones were 20 mm apart with the microphone array stretching from 300 mm
(3 7.5h) to 600 mm (75h) downstream of the fence.
....
n
M-plate
-1
i ,
l
1
i
.1
Figure 2.5. A picture of the fence-with-splitter-plate model inside the wind tunnel
2.1.3. Instrumentation
1. Static pressure
Static pressure measurements were used to align the model with the flow direction
in order to achieve flow symmetry on both sides of the splitter plate. As mentioned
earlier, there were 56 pressure taps (28 on each side) mounted flush with the l-pate
31
ill I" -P JIM
' A" (F.
surt'i
appr
pnet
Ufifll
scar
“11h
08g;
pres:
refer
0P6:
PICS:
Piesl
PIES:
\
surface. The inner diameter of each tap was 1 mm and the outer diameter was
approximately 1.5 mm. Each of the taps was connected with a long Urethane tube to a
pneumatic connector (48D9M-l/2) of a pressure scanner. The scanner, designed and
manufactured by Scanivalve Corporation (48D9-1/2 Scanner Oiless Design with 48
ports), had a 100-psi range and was driven by a rotary solenoid (48D9M-1/2) coupled
with a homemade solenoid-actuation circuit (for more detail regarding the circuit design
the reader is referred to Hudy 2001). The solenoid was made to "step" from one pressure
port to the next by feeding a 5V control signal to the actuation circuit from an analog-
output channel of a PC-based D/A converter. The stepping and pressure-data-acquisition
algorithms were synchronized by containing them within the same data acquisition
program.
To measure the static pressure relative to some reference, the output port of the
scanivalve was connected to a 0-1.33 kPa Baratron pressure transducer (model 223BD)
with an output range of 0-1 Volt. The output port of the scanivalve was connected to the
negative side of the pressure transducer while the positive side was connected to a
pressure tap at the entrance of the test section, upstream of the model, to provide a
reference pressure (p,) for the measurements. The first port of the scanivalve was left
open to the atmosphere while the following ones were connected to the tubes from the
pressure taps. In this manner, the “home position” of the scanivalve gave the differential
pressure between the atmospheric and reference pressure (pm - pr) while the following
positions yielded the differential pressure between that acting at the corresponding
pressure tap and the reference pressure (ps - pr). Forty of the 56 pressure taps (28 taps on
the top side and 12 on bottom) were used to align the model with the flow direction.
32
1
E‘
,.
. a
ll.T‘
esuh
dime
that}
in d:
yam
thetH
caPlU:
II. Tufts Visualization
Tufts were used to visualize the flow direction close to the wall in order to
estimate the reattachment length, which is a relevant length scale for non-
dirnensionalizing the investigated flow quantities (e.g., Ruderich and Femholz 1986). A
black-yarn tuft was attached to the surface of the I-plate at different streamwise locations
in the reattachment zone. The tuft was constructed from a 10 mm-long piece of black
yarn attached to a short thread that was in turn attached to the I-plate using a tape, freeing
the tuft to follow the flow. An image of the tuft is shown in Figure 2.6.
Figure 2.6. An image of a tuft attached to the I-plate
A CCD camera (Sony XC-75) coupled with a flame grabber system was used to
capture the tuft image. The camera had a '/2- inch CCD image sensor with a resolution of
33
nun
pert
the
inta;
alin
lll..
the
Pan'
the
in th
Sensi
dElf
ofth
T6813
then
Emd I
768 x 494 pixels. The frame grabber was a National Instruments IMAQ PC1-1408 board
(for more information regarding the board specifications, the reader is referred to the
subsection V below in the current section, 2.1.3). The idea was to capture enough
number of images in order to statistically determine the dominant flow direction from the
percent of the time that the tuft is imaged in the up or downstream direction. Because of
the high reflectivity of aluminum, black yarn was used to increase the contrast of the
images and make it easier to determine the tuft location by looking for a dark area against
a bright background.
111. Microphones
A microphone array consisting of 16 Panasonic microphones was used to measure
the wall pressure fluctuations downstream of the separating/reattaching flow. The
Panasonic microphone (WM-6OAY) was an omnidirectional back electret condenser
microphone cartridge. Figure 2.7 shows an image that illustrates the geometry of one of
the microphones used in the array. The microphone was 6.0 mm in diameter and 5.0 mm
in thickness. The sensing diaphragm was exposed to the flow through a 2 mm round
sensing hole. The nominal sensitivity was specified by the manufacturer to be —42 i 3
dB for a bandwidth of 20 — 20,000 Hz.
Figure 2.8 shows a schematic of a homemade circuit that was used to drive each
of the microphones. The circuit was powered by 9 V DC power supply connected with a
resistance of 2.2 k!) in parallel with the microphone. The output of the microphone was
then high-pass filtered to remove any DC component in the output (V0). The capacitor (C)
and the resistor (R2) of the high-pass filter were selected for a cut-off frequency of 0.16
34
Hz, which was sufficiently low to block the DC without removing any significant
fluctuating-pressure energy.
Figure 2.7. An image of one of the microphones used in the array
I F=0.luF
R,=2.2 kn
V0
+
R2=10 M9 9 v Mic.
2.]
_I_
Figure 2.8. A schematic of the microphone driving circuit
35
IV.
inye
COlTl
TU
IV. Hotwire
Single- and X-hotwire probes were used to measure the velocity in the
investigated flow. The single hotwire was used to measure the streamwise (u) velocity
component for estimating the boundary-layer mean- and turbulent-velocity profiles. On
the other hand, the X-hotwire was used to measure the streamwise and wall-normal
velocity components simultaneously with the microphone-array data. The probe design
was similar to that of a straight X-probe except the four prongs were bent by an angle of
30°, as shown in Figure 2.9. Bending the prongs in this manner enabled positioning them
as close as 0.5 mm (or less) to the surface of the splitter plate. Furthermore, by making
the bent part lower than the tip of the probe support, the significance of any disturbances
that might be created by the probe support is reduced. Photographs of the probe can be
seen in Figure 2.10.
. .7 * ___i ,7 7 7,. l/x/T
’ Dimensions in mm ./
.. ,./"'/i
Q ’1,’ /
x???“ g “ /
e / '/ / ./
f2 / ,,// /,
° x// if \ ,/ 1’ I ,//
y-level of the tip of 31 /,/ /,/
the probe support / X / \\
\ ,/ ,/ .
/ , /. \\\
Probe Support
Figure 2.9. A schematic of the X-probe
36
norm
consi
brack
frams
hOfiZt
Kid in
C0510}:
“Hire :-
Rater.
“615 1';
SEFPC
motor
Top View
Figure 2.10. Images showing top and side views of the X-probe
The hotwire probes were operated by a Constant Temperature Anemometer, CTA,
from DANTEC (Model 54T30) with an Over Heat Ratio (OHR) of 1.6. The output
signal was low pass filtered using a built-in filter with a 10 kHz cut-off frequency.
A traversing mechanism was used to traverse the hotwire probes in the wall
normal (y) direction. Figure 2.1] shows an image of this traversing mechanism, which
consists of a toothed rod geared with a stepper motor and assembled to an aluminum
bracket. The figure illustrates how the traversing mechanism was fixed to a supporting
fi'ame that was built out of uni-strut structural elements next to the test section. A
horizontal slot was machined in the Plexiglass window to allow protrusion of the toothed
rod into the test section. The hotwire probe was clamped at the end of this rod using a
custom adapter. The position of the probe was then controlled by feeding TTL square
wave pulses to the stepper motor from a digital I/O controller interfaced with a PC. The
traversing mechanism was calibrated against a dial-gauge to determine its resolution. It
was found that the probe translated a distance of 58.8 mm for each revolution of the
stepper motor. This corresponded to a traversing resolution of approximately 74 um per
motor step (since 800 pulses, or steps were required to complete one full motor
resolution).
37
.;.::~:
- u
. 11": 1‘- WI— _
WI
p05
yer
OUI
one
nnc
daui
6024
PYOe.
SUD;
Chan
I0 th
VEiCig
The toothed rod was coupled with a Linear Variable Differential Transformer
(LVDT) sensor, as shown in Figure 2.11. The sensor was used to feedback the hotwire
position to verify the actual movement of the probe. The LVDT calibration was verified
versus that of the stepper motor. Figure 2.12 depicts the calibration results, where the
output voltage of the LVDT (e) is plotted versus the displacement (y). The figure
reveals the linear character of the sensor over a range of about 50 mm. This range was
large enough to accommodate the full movement of the hotwire probe in the y direction.
The slope of the calibration line represents the LVDT sensitivity, which was 0.35 V/mm.
V. Data Acquisition Hardware
Three different data-acquisition systems were used in the present work. The first
one was used throughout the experimental procedure of the wall-pressure measurements,
including aligning the model, calibrating the microphones and acquiring data fiom the 16
microphones. This was accomplished using a National Instruments AT-MIO-16E-10
data acquisition A/D card with 16-single-ended analog-input channels and maximum
sampling rate of 100,000 samples/s (corresponding to 6,250 samples/s per channel). This
sampling rate resulted in 10 us inter-channel time delay.
The second data acquisition card was also a National Instruments A/D card (NI
6024E-ADC), which was employed for driving the stepper motor, hotwire calibration
procedure, measurement of the boundary-layer velocity profiles, and acquiring the
simultaneous velocity and pressure data. The card also had 16—single-ended analog-input
channels with a maximum sampling frequency of 200,000 Hz. It was necessary to switch
to the second A/D board because the stepper motor had to be synchronized with the
velocity measurements for automated operation of the traverse/acquire procedure. This
38
I
1
I-l
III-III
F1
up
was not possible with the first board, which did not have the hardware provisions
necessary to control the motor.
'I
o
o
..
V!
F
0
O
s
'I‘t-
0.
ll.
..
‘I
I
I
.
o
l
...
s
-
Figure 2.11. Images of the traversing mechanism and the LVDT: overall (left) and close-
up (right) views
10 i T j 7'1 fl j ‘V T7 77 If I fl fl fir 7 f T ' ri-T 7 I 7—T "T 7 ‘l—" 7 ifi
g e =9.813—0.35y
. o r
5_ ................ . . . .......
i 0»
6
cu
_5 g ..........
. b
r J : i ,
_lOr_‘4l',L‘hi iJri. -774. ,L Val—‘1‘. ,1 274‘___1L. 2 ,l 4L_,.li A #.__l¥.J . J—; l 4;." a
0 10 20 30 40 50
y [mm]
Figure 2.12. LVDT calibration
39
1.1
j in!" am“ -.- 3"
tuft
11101
yidt
flour
cha.
lite
syst
anal
2.2.
1116;:
Subs
5"
Ix)
”181‘;
mO-tf
desc'
Para
“Er;
The third and last computer-based-acquisition board was used for grabbing the
tuft images. The board was an NI IMAQ PC1-1408 board, which was a high-accuracy,
monochrome, PCI-based IMAQ board that supported RS-170, CCIR, NTSC, and PAL
video standards. The board was connected to a BNC break-out terminal, which included
four video-input channels and four corresponding I/O triggering channels. The triggering
channels are used to start/stop sampling of images or indicate the status of acquisition.
The PC1-1408 acquired image frames in real time and transferred them directly to the
system memory after converting their data into digital formats using an 8-bit flash
analog-to-digital converter (ADC).
2.2. Experimental Methodology
A full description of the procedures used in implementation of the various
measurement techniques in the experiments conducted here is provided in the following
subsections.
2.2.1. Static pressure system
Static pressure was only measured in the separating/reattaching flow region. As
mentioned earlier, static pressure measurements were used for the purpose of aligning the
model parallel to the freestream inside the test section. This subsection gives a
description of the alignment procedure.
1. Acquisition Settings
The acquisitions settings were selected after considering the flow features and the
parameters of the static-pressure measuring system. The main settings of the acquisition
were the sampling time to converge to the average static pressure, and the time needed
after each step of the scanivalve for the pressure inside the pressure-measuring system
40
' “'3‘
I.“ , Him" .
“E
W e
531‘.
err
11.
niot
r Fl
pres
ofti
floyt
usek
Pics
disw
“as
(11811:
(connecting tubing, scanivalve and transducer) to stabilize. Tests conducted by Hudy
(2001) using the same setup had shown that 10 seconds and one second, respectively,
were more than sufficient for these parameters. Therefore, 1000 samples of pressure
were acquired at a sampling rate of 100 Hz before triggering the scanivalve to step. After
the trigger, the program waited for one full second before commencing with data
sampling. The procedure was then repeated for all static-pressure ports to complete the
experiment.
11. Model Alignment
Model alignment was done in two steps: a coarse alignment and a fine one. The
model was coarsely aligned by connecting the tubes of pressure tap # 18 (x/h = 21.4 and
x/Xr = 0.95) on opposite sides of the model to the positive and negative ports of the
pressure transducer. The "angle of attack" of the model was then adjusted till the readout
of the pressure transducer became zero, indicating rough alignment of the model with the
flow direction.
For fine alignment, 40 of the 56 pressure taps (28 on top and 12 on bottom) were
used in the alignment procedure. The taps were used to obtain a “fuller” picture of the
pressure distribution on each side of the splitter plate. By matching the static pressure
distributions on both sides of the model, the fine alignment of the model was achieved.
Figure 2.13 shows the top and bottom static-pressure distributions after the model
was aligned. It is evident that both the top and bottom mean-pressure-coefficient (Cp)
distributions match well, where:
_ P5 ‘Pr
Cp—KPUE.’ (2.1)
41
pr:
fre
bet
dei
Fig
Split
‘ ail
4-
«It
,6
“a1?
p, is the surface pressure at location x, pr is a reference pressure measured with a static
pressure tap located at the exit of the contraction, p is the fluid (air) density, and U.G is the
freestream velocity upstream of the model. The results in Figure 2.13 indicate that the
model is satisfactorily aligned with the direction of the flow. In fact the largest deviation
between the two curves is 3.5% of the total Cp variation (Cp, max - Cp, min), with the RMS
deviation being 2.5%.
'0.2.'" ’ T " 2 6 i ,2 7 fl 7 7 . + v 2?
l 0 “Top . O I C O
1 ‘ ' I ‘
I Bottom .'
T i ' o
-0.4" - , f , , - f
, O
'
DO- ‘06 —._.
I
O
-0.8' -.. Q,
j o
' 0 o I . . or
' a" o
_1‘. 2 , ‘. BAA- ,1_7_ .‘ __ ; _ -.., . ,
0 5 IO 15 20 25 3O 35
x/h
Figure 2.13. Mean-pressure-coefficient distribution for the top and bottom sides of the
splitter plate
2.2.2. Microphones
As mentioned in section 2.1.3, Panasonic microphones were used to measure the
wall-pressure fluctuations. The microphones installation, their calibration, time-delay
check and cross-talk check are discussed here.
42
1.1
and
mi;
the
The
the
1. Installation
Figure 2.14 shows the installation of one of the microphones. To mount the
microphones, counter-bored through holes were made on the back face of the instrument
and middle plates. The smaller diameter of the counter-bored holes matched the
microphone sensing diameter (2 mm), while the bigger diameter matched the diameter of
the microphone casing. The idea was to make the sensing holes completely flush with
the surface, which was achieved by inserting the microphone from the backside of the
instrument and the middle plates.
Panasonic
FIOW I Microphone
E 02 0 '
5 ' mm >‘ < ' Aluminum
" Plate
1 . l . ._ . , I .
-‘ _ " \ \ x x a E‘
r . _ e
‘ '\. O-
E ‘ . '"y
E . . c
tn
5?. £3 Conductive
M epoxy
‘0 6.0 mm;
Figure 2.14. Microphone installation in the Instrument and Middle plates
The microphones were held in place using conductive epoxy (see Figure 2.14).
The epoxy was applied around the perimeter of the microphone casing, electrically
connecting the latter in the process to the model, which acted as a ground plane. Using
the conductive epoxy guaranteed that all microphones had a common ground. All
43
pa
ael
microphone wiring was fed through the space between the top and the bottom plates of
the splitter-plate.
II. Calibration
Although the Panasonic microphones had a nominal sensitivity of —42 i 3 dB
for a bandwidth of 20 — 20,000 Hz, they had to be calibrated to determine the frequency
response of each microphone individually after their installation in the I- and M-plate.
The calibration yielded two important pieces of information. The first was the
determination of the mean sensitivity and phase response of each microphone. The
second was to check that each microphone was mounted properly without a gap between
the microphone face and the counter-bored hole, which could result in a Helmholtz-type
resonance during measurements.
A Plane Wave Tube (PWT) was used to determine the frequency response of the
Panasonic microphones. Basically, a PWT produces plane sound waves that travel
parallel to the axis of a tube or a duct. This one-dimensional planar wave propagation is
achievable if the width of the duct is small in comparison to the acoustic wavelength, 2,.
Specifically, for a square duct with rigid walls and side length of 2a, planar wave
propagation is accomplished when it, > 4a, or f < c/4a (where f is the sound frequency
and c is the speed of sound); e.g., Kinsler et al. (1982). Thus, f = c/4a establishes the
upper limit, or cut-off, frequency of the PWT. In the current investigation, the PWT
width selected is 12.7 mm, which gives an upper frequency limit of approximately 13.5
kHz for a speed of sound of 350 m/s.
A picture of the acoustic—wave-guide setup may be seen in Figure 2.15. The tube
was made of 1.8 m long aluminum square duct with a cross-section of 12.7 x 12.7 mm2
44
A‘
(D
and a wall thickness of 3.175 mm. The centers of the microphones were located along
the centerline of the PWT. Opposite to each microphone, a plug for mounting a reference
Larson & Davis (L&D) l/4-inch microphone (with known response) was made in the
upper wall of the PWT (see Figure 2.15). A sound generating system (function generator,
amplifier and speaker) was used to calibrate the Panasonic microphones. The system was
used to excite the microphones over a broad range of frequencies using white noise
containing frequencies up to the highest frequency of interest in the experiments (z 3.5
l
Spun 1\(‘l‘
l—plule
\lirltlle plate
Figure 2.15. A picture of the microphone-calibration setup
Before using the PWT for calibration purposes, it was necessary to check the
accuracy of the calibration procedure. The test was conducted using two L&D
microphones, where the first microphone was mounted at the upper wall of the PWT,
45
\K
6X1
fro
the
des
while the second was mounted in a plug in the side of the PWT. In this manner the
pressure ratio (R) and the phase shift (0) between the acoustic pressures acting on the top
and sidewalls of the PWT could be determined, where:
—
i
(p'nm).0p and (p'rms),,-de are the RMS values of the pressure fluctuations on the top and side
R (2.2)
wall of the PWT, respectively.
The results provided in Figure 2.16 are the average for two locations along the
PWT, corresponding to the beginning and end positions of the microphone array. Figure
2.16 (top plot) shows the amplitude uniformity of the sound waves over the cross section
of the PWT, as demonstrated by the Rp value of one (within 10%) for frequency range
extending up to 10 kHz. The planar character of the sound wave can also be confirmed
from the phase shift, 0 (bottom plot). The plot shows that 0 stays within a few degrees in
the same frequency range. Collectively, the amplitude and phase results show that the
desired planar sound waves have been achieved inside the PWT.
Using the above calibration procedure, the frequency response of the
microphone/hole assemblies could be determined. Figure 2.17 depicts the frequency
response of one of the assemblies. The top plot in the figure shows that the sensitivity of
the microphone (K) was flat and fell within :2 dB of the estimated mean sensitivity (Km,
which is obtained by averaging over all frequencies within the calibration range) over the
frequency range 50 < f < 5000 Hz. In addition, the plot of the phase shift between the
Panasonic and reference (L&D) microphone depicts a negligible phase delay. The results
46
01
01
~—
TWULMWdi C
verify the appropriateness of the frequency response of the microphone after embedding
it in the counter-bored holes of the I- and M-plates.
Figure 2.18 (top plot) shows the mean sensitivities of all microphones. The range
of sensitivities extends from 8.5 to 12.5 mV/Pa (—41.4 to —38 dB relative to a sensitivity
of 1.0 V/Pa; the nominal value reported by the manufacturer is —42 i 3 dB).
L5
4
H
trii—FiitTttti—rtnrr
...........................................................................................
(15
ALI—LlidLiiJlllLll
200
100 L
I I I I
0 [degree]
I I I I
-100 '—
| | T
-200
0 2000 4000 6000 8000 110
f [Hz]
Figure 2.16. PWT calibration results: pressure ratio (top) and phase shift (bottom)
47
:2.— .4—\.4_ 7:222‘; C
K/Km [dB]
0 [degree]
20 1
15
-10
-15
................
............................................................................................
l . 4
. . 1 .
1 . I l
,_ ............................. . ............................................ , ........................
>— . r . a
. l I
l . .
1 n
1 l
l l a
. I 9 I
r u I l
_5_. .............. . ...... . ......... . .....................................................................
t . .
t . y r
............................................................................................
Figure 2.17. Frequency response of a microphone/hole assembly
48
.....................................................................
##LIA._L_A~L 1214.1
III. Time-Delay Check
The time-delay of the microphone output voltage relative to the measured
pressure was deduced from the phase shift information (see Figure 2.17; bottom plot)
using the following equation:
“39’ 23
Td-dw (-)
where, Td is the time-delay in seconds, a) = 21tf (f is the frequency in Hz) and 6 is the
phase shifi angle in radians.
The deduction of the time-delay was based on estimating the average slope (if!)
in the plot of the phase shift versus frequency. The estimation indicated that Id was z 10
us, which was more than an order of magnitude (z 66 times) smaller than the average
convection time between two successive microphones (z 667 us), calculated based on the
highest U0o of 30 m/s. Figure 2.18 shows the time delays of all microphones.
IV. Cross-Talk Check
Because of the relatively large number of signals acquired simultaneously in this
study, checking the cross-talk among different acquisition channels was necessary before
sampling data. The cross-talk check was conducted by exciting only one of the
microphones by a sound source and determining the square of the RMS pressure
measured by all other microphones relative to the square of the RMS pressure measured
by the excited channel. The procedure, which utilized a Larson and Davis (L&D) hand-
held calibrator (model CAL 200), was repeated for all microphones in the array. The
calibrator was capable of generating sound waves at a frequency of 1 kHz and one of two
49
3!
TE
Sound Pressure Levels (SPL) selectable by a switch: 94 and 114 dB (relative to a
reference pressure of 20 uPa). To ensure that sound was only applied to a single
microphone, a short, flexible Urethane tube was connected to the cavity of the L&D
calibrator on one end, while the other end was pressed against the top of the excited
microphone only. Note that the test conducted in this manner enabled assessment of the
cross-talk of the whole system end to end; i.e., including the microphones while.
embedded in the I— and M-plates, the driving/conditioning circuitry, the data acquisition
terminal block, and the A/D card.
Figure 2.19 shows the cross-talk test results for the case when microphone # 13
was excited. This case is presented here because it had the most (i.e., worst) cross-talk
amongst all 16 microphones. The data used in the plot was obtained by exciting
microphone # 13 using the calibrator, as mentioned earlier, while simultaneously
acquiring data from all 16 channels. A signal attenuation ratio (SAR) was then calculated
by dividing the square of the RMS pressure measured at each channel by that of channel
13. The plot shows that the ratio is zero dB (decibels) at microphone # 13, which is
consistent with the fact that microphone # 13 is the excited microphone. The strength of
the signals measured by the other channels was at least two orders of magnitude (SAR < -
40 dB) less than the source-signal. That indicated that the cross-talk among channels was
negligible.
2.2.3. Tufts
Tufts were used to visualize the near-wall flow in order to estimate the
reattachment length (X,). As mentioned earlier, Xr is considered an important length
50
scale for characterizing separating/reattaching flows. The procedure employed for using
tufts to evaluate Xr is presented here.
15% T I I I W I I I I 7 f I I T
LIII.IOIiOOIOUOOIdEIIIIIII-é- é -£—=-IIIII:‘IIIIIII.: 3’ 3;--- 2 ' ‘-Ji‘ ““““ I .‘IIOII
L. . I - I - I -
10‘“ fi ---------- O ------------ ------------------------------------------------------
o . ’ f 9
m I . . Q .
& :Illounnq nnnnnn one nnnnnnnn ‘po-o-n-r: ; T---=-q' ----- ‘ -',_ ---u: nnnnnnnn Q nnnnnnnn Ill-IIII‘: IIIIII coo-nun..-
> : . I
E t .‘ '
ME ; Nomlnal range 5
F of sensitivity i
5L— .............................................................
. I
r
E 1
~ I
>— . ‘ 4
0; I I L I I I I I I I I I I I l
101 l 7 l l I I I I I I I I I f
I E i i '0 i E I
r : ' ' 3 . .
8,...n-..-....-, ................................ J ..... .-.-.-W.-., .......... . ..... L ...... ..... _'1
Z c r , 2 ¢ 26
7 F 1 . . I
I; 6%... .. ...... ................................ _,
:9 P . ' _
—' I. O 5 I
45’ , ; 9 o -
a) ' . : .
g 4“,... ................. . .. ................................... ‘l
E--‘ ' : 3 I
_ . l
t : . I
2_ ..................................... - ....... I ...................................................... j
; I
0 I I I I I I l I I I I I I I l
1 2 3 4 5 6 7 8 9 10 ll 12 13 14 15 16
Mic#
Figure 2.18. Mean sensitivity (Km) of the microphones (top) and corresponding time-
delay (bottom). Microphone #1 is the most upstream microphone
51
SAR [dB]
-30;
'40:” , . . . .
-60 - ~ - 7.1 .,, - -t v: .* Li W . 2 fig H . a z, A .
2 7
Mic #
Figure 2.19. Cross-talk check of Microphone # 13 (Mic. # 1 is the most upstream one)
1. Procedure
A 10 mm-long black-yam tufi was attached to the surface of the I-plate using tape
at different streamwise locations in the reattachment zone. The streamwise location of
the point of attachment (root) of the tuft was determined using the image shown in Figure
2.6. The root point was defined as the point at which the free part of the thread
intersected with the edge of the tape. The corresponding x location of that point could be
determined from knowing the x location of the closest pressure taps captured in the image
(see Figure 2.6).
The required number of tufi images was determined by testing the effect of the
number of the acquired images on the estimated Forward Flow Probability, FFP (to be
explained later). The test results showed that 5000 images were enough for the FFP to
converge within 5% percent, which was determined by acquiring 10,000 images. During
recording, the tuft was illuminated using a fiber-optic strobe light synchronized with the
52
frame grabber system, while the shutter of the CCD camera was set to "always-open"
mode. The synchronization was accomplished by triggering the strobe light using a
signal generated from the image-grabbing LabView program through the triggering
channel of the image acquisition board. This enabled “freezing” of the tuft location in
each image.
A sample image of the tufi is shown in Figure 2.20. The image was captured after
enlarging the aperture of the CCD camera lens relative to that used to capture the image
shown in Figure 2.6. The idea was to force the background of the image (the aluminum
surface) to be completely bright, leaving the black tuft highly visible in the foreground.
In processing the images, the main idea was to determine to which side the tuft was
located, relative to a "splitting line", which passes through the root of the tufi, as
exemplified in Figure 2.20. If the tuft was to the left of the line, this would indicate a
backward (upstream) flow, and vice versa. By using a MatLab program, each image was
split into right and left "halves" based on the splitting line. Then, the average intensity of
pixels in each half was calculated. The program was then able to indicate on which side
the tuft was located by detecting the image half with lower mean intensity. The
information from all 5000 images was used to determine the mean flow direction.
II. Forward Flow Probability (FFP)
Forward Flow Probability (FFP) is defined as the probability of; i.e., the fraction
of the time that, the flow is in the forward (downstream) direction. If FFP equals 1.0, the
flow is always forward. A zero FFP, on the other hand, reflects a flow direction that is
always in the upstream direction. Accordingly, the mean point of reattachment is located
by finding the x location where FFP equals 0.5.
53
§\Splitting
Line
Figure 2.20. A sample wide-aperture tuft image
The F FP information was obtained for eight different x locations and the results
are shown in Figure 2.21. The FFP value varies from around 0.1 to 0.9 at x locations
extending from 161 mm to 199 mm downstream of the fence. The data for all eight
locations were curve-fitted using a third-order polynomial function (the solid line in
Figure 2.21) and the location of FFP = 0.5 on the fitting curve was found at
approximately 180 mm downstream of the fence. For comparison purposes, the
reattachment length was also estimated from the data presented by Smits (1982) for X,
versus blockage ratio for a similar flow geometry. This yielded an X, value of 168 mm
(which is 7% less than that obtained from the tuft data). Additionally, similar data from a
shorter S-mm tuft yielded an X, value of 183 mm, which is within 2% of the value found
from Figure 2.21. Considering all three values, it appears that the uncertainly in the
estimated X, value is of the order of 5-10%.
54
—" *_ ~m--.~A.?_L -- _ ......hh
] llllIllVrfTrlellfiTfilljlfiIfifT—TIITjTTYTT—T
l
l
i 1
fl 0 5 _ l
LL. ' - I
I i
I _
I I
O P- 1 1 L; l l l l [J J l l 1 l l i 1 1 l 4 l J 1 A L l l l l l l l l j l l l J
160 165 170 175 180 185 190 195 200
x [mm]
Figure 2.21. F FP of the reattaching flow versus the x-location of the tuft
2.2.4. Hotwire Sensors
Both a single- and X-hotwire probes were used to measure the flow velocity in the
developing boundary layer at a single x location downstream of the reattachment point.
The procedures for the magnitude and yaw calibration as well as for determination of the
probe height above the wall are presented below.
I. Procedure
A 5 um-diameter annealed tungsten wire was used to construct the hotwire
sensors. First, a short length of the wire (z 50 mm) was electroplated with copper in a
copper sulfate solution, except for a short bare length in the middle that formed the
sensing length of the wire. This sensing portion was 1 mm long, resulting in a length-to-
diameter ratio (l/d) of 200. After plating, the wire was directly soldered to the probe
prongs, resulting in a typical sensor resistance of around 5 ohms.
55
The procedure adopted to measure the flow velocity using the hotwire probes may
be summarized as follows (note: for the single wire, only procedure steps 2 through 5 are
pertinent):
1. Conducting yaw calibration of the X-wire probe to determine the angle a1+a2
between the two wires forming the X configuration (see Figure 2.22 for definition
of the angles), and identify the relation between the effective cooling velocity (U,)
for each of the wires and the velocity components u and v.
2. Positioning the probe at the desired y location above the wall and determining the
angles of the wires of the X-probe relative to the x axis
3. Calibrating the probe in-situ at the specific streamwise location where
measurements are to be conducted (pre-calibration)
4. Conducting velocity measurements
5. Calibrating the probe after measurements (post-calibration)
II. Yaw Calibration
Yaw calibration means determining the relation between the effective wire
cooling velocity magnitude, on one hand, and the magnitude of the velocity vector and its
angle relative to the wire, on the other. Such a calibration can be accomplished by
changing the yaw angle of the X-wire probe over a certain range in a known-velocity
(direction and magnitude) flow and recording the output of both wires at every angle.
Using the resulting voltage versus angle information, the sought relation can be
determined. For a single-wire probe, the wire is aligned perpendicular to the flow
direction, experiencing the most cooling influence and resulting in the maximum output
voltage for the given velocity magnitude. In this case, the angle between the velocity
56
vector and the normal to the wire (in a plane containing both the wire and velocity vector)
is zero. If there is an angle between the flow direction and the normal to the wire, the
wire primarily 'feels' the cooling effect of the component normal to the wire, with very
little cooling resulting from the component parallel to the wire axis. The velocity
corresponding to the net cooling influence is referred to as the effective cooling velocity
(Ue). Ideally, for a wire with an infinite length-to-diameter (l/d) ratio, the cooling effect
of the velocity component parallel to the wire is zero and the relation between U, and the
actual velocity is a “cosine law”. For example, under ideal conditions, wire #1 of the X-
probe shown in Figure 2.22 experiences a cooling velocity of u cos(90-0t1) associated
with the x component of the velocity (u).
U / q
' ‘ v
n uFC _,
Wire #1 Wire #2
Figure 2.22. A schematic of the X-wire probe and angles definition
57
In reality, the length-to-diameter (l/d) ratio is finite and the cooling influence of
the velocity component parallel to the wire causes a deviation from the cosine law. In
this case, the effective velocity should be related to the actual velocity by a more general
yaw response function (F2(90-a); the square exponent of F indicates positive definiteness)
instead:
U. = U F2(90-oc) (2.4)
Where, U is the flow velocity magnitude and a is the angle between the wire axis and the
velocity direction.
The form of the yaw-response function used in the present study is that proposed
by Champagne (1965):
F2(90-a) = [0032(90-0t) + k2 Sin2(90-a)]“2 (2.5)
where, k is a constant that is determined from yaw calibration. Note that when k equals
zero, the cosine law is recovered.
Figure 2.23 shows an image of the yaw-calibration setup that was used in the
present study. The setup included a rotating stage with an angle indication of one degree
resolution. The rotating stage, which was used to control the yaw angle of the probe, was
installed above the upper wall of the wind tunnel test section, and was connected to the
probe holding fixture in the wind tunnel via a 1/4" threaded rod that passed through a slot
in the test section ceiling. The fixture holding the probe mainly consisted of a horizontal
aluminum bar and a vertical cylindrical rod that were almost of the same length of 0.15 m.
Together with the X-wire probe, these elements formed a U shape, as evident in Figure
2.23. The main purpose of the U-shape design was to keep the X-wire sensing part on
axis with the rotating stage, thus allowing a change in the yaw angle while keeping the
58
center of the X-wire at the same location. This is demonstrated more clearly in Figure
2.24.
- was” mm.
Figure 2.23. An image of the yaw calibration setup
Figure 2.24 provides a schematic of the setup for yaw calibration. The angle 9
shown in the figure is that between the flow direction and a line drawn from the center of
the X-wire to the center of the vertical rod. 6 is used here to indicate the rotation of the
wire relative to the freestream. Note that because of the bending angle of the sensor
prongs, discussed earlier, when 9 = 0, the axis of the cylindrical body of the probe makes
a very shallow angle relative to the flow direction. Also note that 9 does not represent
the yaw angle (which is that between the normal to the wire and the flow velocity vector),
but 9 increments equal yaw angle increments for both wires, with the largest wire output
voltage at a given freestream velocity corresponding to a yaw angle of zero. In
59
implementing the yaw calibration procedure, the setup was initially adjusted such that the
readout of the rotating stage was zero when the horizontal rod was visually aligned with
the x direction (flow direction). The stage was then rotated through a large negative
angle (—110 degrees) to start the calibration. At this angle, the corresponding hotwire
output voltage was recorded. The procedure was then repeated for a total of 55 angles
obtained by incrementing 0 in the positive direction in increments of 5 degrees up to
+160 degrees.
Figure 2.25 shows the primary output of the yaw calibration. The figure depicts a
plot of the output voltage versus 0 for both wires of the X-probe. The results were used
to determine the zero-yaw angle for each wire; i.e., where the sensor output reaches its
peak, which was approximately 0 = —15 and 75 degrees for wires #1 and #2, respectively.
The difference between these two values equals to 180-(11—(12 (see Figure 2.22),
showing that the two wires are in fact very close to being perpendicular to each other as
intended in the construction of the probe (note that the angle between the two wires is
a1+oc2).
To determine the value of the constant k in equation 2.5, each wire had to be
calibrated individually at zero yaw angle (i.e., 9 = —15° for wire #1 and 75° for wire #2,
as discussed above). The hotwire calibration was conducted by recording the mean wire
voltage (E) at different known values of the flow velocity (U). The data, a sample of
which is shown in Figure 2.26, were fitted using a fourth-order polynomial fit of the form:
u = a0+a1E+... (further details of the calibration process is given below). Using these
polynomials, it was then possible to convert the voltages shown in Figure 2.25 to velocity
60
and replot the results as U,,/U.O (or effective to actual velocity ratio, where U.0 is the
freestream velocity in this case) versus 9 in Figure 2.27.
Figure 2.28 shows the data points from Figure 2.27 replotted as Ue/U.o versus
(90—0L) for wires #1 (top plot) and #2 (bottom plot). This is accomplished by subtracting
the 0 value corresponding to the peak wire output from all 0 values. Plotted in this
format, the data in Figure 2.28 can then be fitted with an equation of the form given by
equation 2.5. The equation has a single adjustable parameter: k. A second parameter, ore,
representing an angle offset error in determination of the 0 value corresponding to the
peak wire output (zero yaw) was also added in order to determine the zero-yaw angle
with resolution better than the 5-degree increment used in the yaw calibration. The best
values of these two parameters were found by varying each of them independently over a
predetermined range, and calculating the corresponding average absolute error between
the data points and the fitting curve. The pair of values corresponding to minimum error
for each of the wires was then selected for the best fit. For reference, the resulting fits are
displayed in Figure 2.28 with the data. The corresponding parameters are: k1 = 0.04, kg
= 0.0, and on+ocz = 93.
Once k1 and kg were determined, the following equations were used to calculate
the streamwise and wall-normal velocity components from the effective cooling
velocities of the two wires, respectively (for derivation of the equations see Appendix A):
u 2 U... F301.) + U. ma.)
F.2(90— a.) F5612) + F301.) F§(90— a.)
(2.6)
v _ U... F12(90-aI)-U.IF§(90-az)
F12(90- “1)F22(a2)+ 1212(0'1) F2200 “ a2)
61
E [V]
Vertical rod
Probe ‘ \
\ ‘ \7 ‘ .
.< 7 a ,3 93 m >— 'I
Center of rotation
, +
Figure 2.24. Schematic of the yaw calibration setup
1.9 ' ' ‘ ' V I V
I 180 — a, - a2 : —o— Hot-Wire #1
H—> '
1 35 . 7 —I— Hot-Wire #2
50 A 100 150 200
9 [degree]
Figure 2.25. Yaw calibration: voltage results
62
U [m/s]
Ue/Uao
30;
25:"
20;
1.2'
0.8‘
0.6 '
0.4 " "
0.2 ‘
+HW#1
Figure 2.26.
-150 -lOO
1.7 1.8 1.9 2
E [V]
Sample calibration data for the X-probe hotwires
180-(11-(12 . _HOt-WlfC#l
—l— Hot-Wire #2
JL
-50 0 50 100 150 200
0 [degree]
Figure 2.27. Yaw calibration data after conversion of voltage to velocity
63
80
80
90—0t
Figure 2.28. Yaw-calibration data and corresponding curve fits for hotwire #1 (top
plot) and hotwire #2 (bottom plot)
64
where, U61 and U62 are the effective cooling velocities measured by wires #1 and #2,
respectively, and F12 and F22 are given by:
F12(a1)=1cos2(a1)+ k3- sin2(0c1)]m (2.7)
F22(0t2) = [0082012) + k2 sin2(0t2)]”2 (2.8)
The method used to determine on and or; will be discussed below.
[11. Probe Installation
The main purpose of the probe installation procedure was to position the probe at
the desired location (between microphones #13 and 14, or 550 mm (3.05 X,), downstream
of the fence), and to adjust/determine the angle (in the x-y plane) between the wires of
the X-probe and the x direction. The alignment of the probe prongs parallel to the splitter
plate at that location was achieved with the aid of a small piece (approximately 25 mm by
25 mm) of a 23 um-thick polymer foil with a thin layer of aluminum deposited on one
side. The polymer side of the foil was adhered to the splitter plate surface directly below
the probe to provide a mirror image of the prongs. Figure 2.29 shows a sample of the
bent probe prongs and their reflected image. By adjusting the probe till the real prongs
and their reflected image were parallel, the probe alignment was achieved.
After the alignment of the prongs, their image was also used to position the probe
at the desired starting y location. The main idea was to measure the distance between one
of the bottom prongs (the one closest to the wall) and its reflected image. Hence, the
corresponding y location of the prong was the difference between half the measured
distance and the thickness of the foil (23 um). To measure the distance between the
prong and its image, each image had to be calibrated (estimating how many pixels in the
image per mm in real dimensions). The image was calibrated by causing the traversing
65
mechanism of the hotwire to move a known distance (typically 1 mm) and measuring the
corresponding prong displacement in pixels, yielding an image scale factor in pixels per
mm. This scale factor was then combined with the measured distance between the prong
and its reflected image in pixels to position the bottom prongs 0.5 mm above the splitter
plate surface. Given the geometry of the X-wire probe, the corresponding y location of
the center of the sensing parts of the wires was 2.0 mm (0.5 + 1.5 mm). For the single
wire probe, the same procedure was used to position the single wire at a starting y
location of 0.5 mm. Because the wire of the single hotwire was parallel to the splitter
plate, it was used as a reference to determine the y location of the probe.
Figure 2.29. A photograph of the X-wire probe and its reflection in the splitter plate
66
Figure 2.30 shows a zoomed-in image of the X-wire sensor. The image was used
to measure or] and or; (required for evaluation of equations 2.6 and 2.7). Once the probe
was aligned as explained in the previous paragraph, the prongs of the probe became
parallel to the x-direction. Thus, it was possible to determine on] and 012 by measuring the
angle between the wires and the prongs as depicted in Figure 2.30. The figure
demonstrates how the angles were calculated by determining the pixel coordinates of the
three points (1, 2 and 3 for or], and 1, 4 and 3 for (12) for each of the angles. The
measured values of or] and 012 were 47° and 46°, respectively. Note how the sum of these
two geometrically determined angels equals to 93 degrees, which is in very good
agreement with the value determined earlier of the angle between the two wires from yaw
calibration.
Prongs
1 mm 6) j *
Figure 2.30. A close-up image of the X sensor
67
IV. Hotwire Calibration
In every test, the hotwire probes had to be calibrated before (pre-calibration) and
after (post calibration) velocity measurements. It was desired to accomplish this in-situ
by positioning the probe in the freestream above the re-developing boundary layer, where
the flow is steady and laminar and the velocity is known. However, the thickness of the
boundary layer (which was about 80 mm, as will be explained later in section 2.2.5) was
larger than the farthest possible position the traversing mechanism can reach away from
the wall. Thus, it was not possible for the probe to reach the freestream with the existing
experimental setup. This problem was overcome by replacing the fence with a sharp (V-
shaped) leading edge during calibration to reduce the boundary layer thickness by
eliminating the bluff-body effect of the fence.
An actual image, along with a schematic diagram, of the “v-attachment” used in
the calibration may be seen in Figure 2.31. The image provides a side view while the
schematic yields a frontal view of the attachment. The modification of the leading edge
resulted in diminishing the thickness of the boundary layer to approximately 25.0 mm. In
order to measure the freestream velocity (which is used as a reference in the hotwire
calibration) at the specific x location where the hotwire measurements were conducted, a
pitot-static probe located in the freestream at this x location was calibrated against a
static-pressure tap in the sidewall of the test section upstream of the model prior to
installation of the hotwire. During the actual calibration procedure, the static pressure tap
was then used to infer the freestream velocity “seen” by the wire.
68
V-attachment I-plate
Flow I \ i a
Figure 2.31. Leading edge for X- and single-wire in-situ calibration
The procedure for the hotwire calibration was as follows:
1. The wire mean output voltage was acquired at eight different wind tunnel
freestream velocities. The corresponding velocity magnitude was measured using the
output pressure of the reference pressure tap.
69
2. To deduce the effective velocities of both wires (U61 and U62) at each calibration
point, U.o together with the known values of (11 and (12, were plugged in:
Uel = U... F12(90-oc1) (2.8)
U,2 = U... F22(9O-a2) (2.9)
3. The effective velocity of each wire was plotted versus the output voltage and the
resulting data were curve-fitted using a 4th order polynomial to provide an analytic
calibration equation for converting the wire output voltage to effective cooling velocity.
The calibration procedure of the single hotwire only included steps 1 and 3 of the above
procedure.
V. Boundary Layer Profiles
A single hotwire was used to obtain the boundary-layer velocity profiles (mean
and RMS) at the selected streamwise location of x = 550 mm (x/XIr = 3.05 and x/h =
68.75) downstream of the fence. The velocity was measured at 182 wall-normal (y)
locations from y = 0.5 to 91 mm in increments of 0.5 m. 10,000 samples were acquired
at each y location with a sampling frequency of 1000 sample/s to construct the boundary
layer profiles as shown in Figures 2.32 and 2.33.
Figure 2.32 depicts the mean-velocity profile of the boundary layer. The ordinate
is the mean velocity in the x direction (11) normalized by the local freestream velocity (U0o
= 16.2 m/s), while the abscissa is the physical y locations of the hotwire normalized by
the boundary-layer thickness based on 0.99U.o (5 = 80 mm). The shape of the profile
generally agrees with the typical shape of a developing boundary layer downstream of a
reattaching shear layer (a detailed discussion of the profile and comparison to the
literature may be found in section 4.1.1). For verification of the X-wire measurement
7O
procedure, which is substantially more elaborate than the single wire, data from the X-
wire are also plotted using open circles in Figure 2.32. The X-wire was traversed from a
y location of 2 mm to 51.5 mm (y/5 = 0.025 to 0.64) in steps of 0.5 mm (Ay/B = 0.00625).
The figure depicts a generally good agreement between the two profiles with a maximum
error of 3% of Ugo.
0.4- .....
o.3~ ..... -
0.2: ......
01 - smgfimé
, w , 1 ..X-Yzirez .
00 ‘02 A04 T 0.6 T 0.8 i ’ 1.2
Figure 2.32. Boundary-layer mean-velocity profiles from single- and X-wire
measurements at x/Xr = 3.05
Figure 2.33 shows the RMS values of the flow velocity in the x direction. The
ordinate represents the RMS velocity normalized by U30 and the abscissa is the same as
Figure 2.32. The RMS profile shows a plateau of high RMS value of approximately
71
° Single wire
I f: X-Wire
8 0.08
rms
o
o
:3 0.06 ‘ - . . .415"- ,
0.04z’ * 7 7 7 ,, 4.5%;
0.02
(i i 0.2 0.4 H 8 0:6 0.8" i 1 1.2
Figure 2.33. Boundary-layer RMS streamwise-velocity profiles from single- and x-wire
measurements at x/Xr = 3.05
0.12Uoo. Beyond a certain height, the RMS value decays as the probe moves away from
the wall. Very close to the wall, within a narrow region of y/8 < 0.025 or so, there is a
rise in the RMS value above that given by the plateau. The RMS values calculated from
the X-wire data are also plotted in Figure 2.33. The figure depicts that the two profiles
agree fairly well. The estimated maximum deviation between the two profiles is less than
1% of Ugo. This deviation, which depicts the X-wire data to be consistently lower than
that of the single-wire data, is not surprising since the X-wire averages the measurements
over a spatial sensing volume of approximately 1 m3, while the single wire averages the
measurements over a 1 mm length in the Spanwise direction (the wire dimensions in x
72
and y directions may be neglected relative to the flow scales). The larger sensor volume
of the X-wire, particularly in the y direction results in attenuation of the small-scale
energy. This is most evident near the wall, where the X-wire data decreases in value as y
approaches the wall, while the single-wire RMS data rise.
VI. X—Probe Disturbance Check
Because of the simultaneous wall-pressure and velocity measurements, it was
necessary to check if any significant disturbances caused by the presence of the X-probe
in the flow contaminated the wall-pressure data. The check was done by examining the
microphones output when the X-wire probe resided at different y locations. Figures 2.34
depicts the power spectra of the microphone voltage output for three different heights of
the X-probe (y = 4, 10 and 50 mm). It should be noted that these three locations indicate
the height of the center of the probe (the sensing part), while the corresponding locations
of the closest prong to the splitter plate were 2.5, 8.5 and 48.5 mm, respectively. Figure
2.34 shows the spectra of the pressure measured by the microphone closest to the hotwire
location (x/Xr = 3.0) on the upstream side. In interpreting the data, the farthest y location
was considered as the no disturbance case. The figure depicts good agreement among
the three spectra, which indicates that the microphone measurements are insensitive to
the hotwire location, suggesting that the presence of the probe in the flow does not create
any undesired contamination of the microphone data.
73
10
-1
.— 10 3'
N ;~
«3 T.
E. ii
I“ *-
9“ -2
'9‘ 10 I:
-3
10 .
'4: . . ... 3 3 3... , .
10 43 - ~* 3 ~ 4* .3 3 . 3 , - +- -3
10 10 10 10
f [HZ]
Figure 2.34. Pressure spectra of the microphone immediately upstream (x/Xr = 3.0) of the
X-wire for three y locations of the wire
Moreover, the same test was done for a microphone immediately downstream of
the hotwire. The test results are shown in Figure 2.35. The data confirm that the hotwire
probe effects on the pressure data can be neglected for the range of frequency of interest.
Note that the same test was conducted for all of the microphones, but only the data for
microphones #13 and #14 were presented because they were the closest microphones to
the hotwire probe, and hence they should suffer the most from any interfering effects of
the X-probe.
74
0
10
-10
F, 10 z:
N
a .
E:
3:.
CL -2
6‘ 10 "f
-3
10
'4: . ; . , , . .. . . _ .
10 3.1 , ,z , - , . .m.2 as i .i. m -A3 W -, _. . W. ,. 4
10 10 10 10
f [Hz]
Figure 2.35. Pressure spectra of the microphone immediately downstream (x/Xr = 3.11)
of the X-wire for three y locations of the wire
75
3. WALL-PRESSURE RESULTS
This chapter focuses on the results of the wall-pressure measurements only. This
focus is important for two reasons: (1) one of the primary goals of this work is to
understand the nature of the wall-pressure field beneath the recovering boundary layer
and compare this to the fully-developed boundary layer, on one hand, and the
separating/reattaching flow, on the other; (2) because in the following chapter the
velocity field will be estimated using Linear Stochastic Estimation (LSE) based on wall-
pressure data, understanding of the characteristics of the latter is an important precursor
to proper interpretation of the estimated velocity field. Of course, an understanding of
the wall-pressure/velocity field relationship is also important for interpretation of the LSE
results. This relationship will be explored in detail in Chapter 4 prior to implementation
of the LSE procedure.
As explained in Chapter 2, the mean wall-pressure was measured using 28
pressure taps located on the center of the splitter plate, while the fluctuating wall-pressure
signature underneath the flow was captured using a l6-microphone array downstream of
the reattachment point. It should be noted that due to the lower limit of the microphones'
operating frequency range (20Hz), all microphone data are high-pass filtered at 20Hz.
This removes any ambiguity concerning the magnitude and phase of p' in the analyzed
data. As will be demonstrated herein, the wall-pressure spectral peak is in fact captured
within the analyzed band of frequencies. Moreover, there is negligible amount of
fluctuating—pressure energy lost because of the filtering process. In the following
subsection, data-processing details and results will be presented concerning the mean-
76
and fluctuating-pressure streamwise variation, autocorrelation, power spectra, cross
correlation, and wavenumber-frequency spectrum.
3.]. Mean-pressure distribution and the reattachment length
Figure 3.1 shows a comparison between the mean-wall-pressure coefficient (Cp,
see equation 2.1) distribution from the present work and that from Hudy et al (2003).
The ordinate represents the pressure coefficient while the abscissa represents the
streamwise location of the pressure taps (x) in mm with respect to the fence. The test
model used in Hudy et al. (2003) was the same as the current one except that their
blockage ratio was 1.94% versus 4.5% in the present work. This blockage ratio is
defined based on the fence height (h). The actual blockage in the present flow (based on
the total fence height, 2H) was 9.9% (compared to 4.26% in Hudy et al.). Because of the
higher blockage ratio and consequent flow acceleration, the distribution of C], from the
present data shows a more pronounced negative valley downstream of the fence (note
that in both studies the reference pressure in Cp is taken at the beginning of the test
section, or upstream of the model). Furthermore, the pressure recovery is less than that of
Hudy et al. (2003). This pressure recovery is faster in the present study, resulting in a
shorter reattachment length (180mm, 22.5h, as found using tufts, compared to 205mm,
25.5h, in Hudy et al. 2003).
The observation that the mean reattachment length of the present work is shorter
than that of Hudy et al. (2003) agrees well with the results of Smits (1982), who
investigated the effect of the blockage ratio on the mean reattachment length. Smits’
results showed that the reattachment length decreased for larger blockage ratio.
Furthermore, when the blockage ratio of the present model is used in conjunction with
77
Smits’ data to estimate the mean reattachment length, an XIr value of approximately 21h
is found, which is shorter than, but falls within 7% of, the measured reattachment length.
3T 3 PfieSent T T T3
-0.1 '4 3- 'I-Iudy et al. (2093) 3
-0,2.--
03-, ,,,,., ,3,3i7fi177,.ln
-o.3 ~
.0.4 ~
on -0.5 ~
-O.6 r
I :.‘1T.TT..
‘9 J .....
' l
I ‘7
9
if 2‘0TT ST 25 T ”T 3.0 TT T 355
Figure 3.1. Streamwise distribution of the mean pressure coefficient from the present
study compared to Hudy et al. (2003)
Ruderich and Femholz (1982) also used a fence-with-splitter—plate model with
two fence heights above the splitter plate (h) of 11 and 22mm. The cross-section height
of the wind tunnel used with their experiments was 0.5m, while the splitter plate
thickness was 6.3mm, which resulted in total fence heights (2H) of 28.3 and 50.3mm.
Consequently, the blockage ratios based on h = 11 and 22 mm were approximately 4.4%
and 8.8%, while the actual blockage ratios based on 2H were approximately 5.7% and
10%, and the corresponding mean reattachment lengths were 22.6h and 17.2h,
respectively. Therefore, Ruderich and F emholz’s results depict an agreement with that of
Smits that the mean reattachment length decreases for higher blockage ratio.
78
Furthermore, the mean reattachment length of the present work agrees well with that of
Ruderich and Femholz (1986) for the case of h = 22mm, which approximately have the
same blockage ratio based on h as the current work. Additionally, Castro and Haque
(1987) reported a mean reattachment length of 19.2h with a 6.2% blockage ratio based on
b.
Note that due to the flow acceleration caused by the blockage effect of the model,
the local (at the location of hotwire measurements) freestream velocity will be used as U.o
for the remainder of this thesis. This local value is approximately 14% higher than the
approaching freestream velocity.
3.2. Fluctuating pressure distribution
Figure 3.2 shows the streamwise distribution of the fluctuating-pressure
I
p [1115
coefficient (C3. = 1 U2 , where p'rms is the root mean square (RMS) of the fluctuating
2 co
pressure). The figure depicts the profiles of Cpl versus g, where Q = (x-X,)/X,, for
Reynolds numbers 7600 and 15700 based on step height. The distributions depict a
decay in the pressure fluctuation with increasing distance from the reattachment point.
This suggests that the flow structures dominating the wall-pressure generation are
continually loosing energy as they travel downstream of the reattachment point. This is
consistent with the relaxation process of the shear layer vortices in zone VI (see Figure
1.1) referred to in the Chapter 1. More specifically, the high mean shear stress energizing
those vortices during their earlier development in the free shear layer (before
reattachment) slowly gives way to a much weaker mean shear as the attached shear layer
gets thicker with increasing x. Thus, the source of turbulent energy received by these
vortices becomes weaker and weaker and so does their energy as well as their wall-
79
pressure signature. Furthermore, the plot in Figure 3.2 shows that at the higher Reynolds
number p' is more energetic and its rate of decay is slightly higher than that at the lower
Reynolds number. This decay presumably continues monotonically till p' reaches a level
corresponding to that of an equilibrium turbulent boundary layer.
01 T ' T ' T T T' ' TT T T T TTITNW ,T‘ 3;T : , '7
I Re= 7600
009-- 4 » y .3 Re=15700
0.08
0.07
0.06
0‘“- 0.05-
0.04~
0.03
0.02- - ~ ~— —- - —
,0. 292211001129...“ ........... .
boundary layer
012%. H ,A i - .- L, ..__ .__f__l____. _ ,-__ ..___ 1,2 , 7 -j
0.5 1 1.5 2 2.5
5.
Figure 3.2. Streamwise distribution of the coefficient of the RMS pressure fluctuation
Figure 3.3 depicts the full streamwise distribution of va in both the separating/
reattaching and redeveloping flow zones at Reynolds number of 7600. The solid squares
represent data obtained in the separating/reattaching flow by Hudy et a1. (2003), while
the solid circles show those from the present work. A polynomial fit is also added to the
plot for the purpose of showing the general trend in the streamwise profile of Cp'. The
profile shows the typical behavior of p' in a separating/reattaching flow before and after
reattachment, which is represented by a rise in the wall-pressure fluctuation till it reaches
80
a peak slightly upstream of reattachment, followed by a decay during the relaxation
process downstream of reattachment.
02 T T ' T T T,if’__"',_, T 1’ T T: T i ,
I Present (Re = 7600)
018- j , I Hudy et al. 2003 (Re = 8000) »-
. 3 — Curve-fit
0.16 - - — - :1: [arabesand QaSér¢11M1986> :
0.14- 7
0.121 *
0‘“ 01
0.08
-1 -05 0 0.5 T 1 T 1.5 2 2.15
Figure 3.3. A comparison between the va streamwise distribution of Hudy et al. (2003)
and the present work with that of Farabee and Casarella (1986)
Figure 3.3 also depicts a comparison between the Cp' streamwise distribution of
the fence-with-splitter-plate model used in the present work with that of a backward-
facing step flow from Farabee and Casarella (1986). Both distributions exhibit the same
qualitative behavior. However, the present flow exhibits much more energetic wall—
pressure fluctuations than that of the typical backward-facing step. This is believed to be
caused by the much stronger shear across the very thin, laminar boundary layer at
separation in the fence flow, in comparison to that experienced across the much thicker
turbulent boundary layer of F arabee and Casarella (1986). Considering that the stronger
81
the shear (ii-E) the higher the turbulent energy production term MR7 $5, the strong
Y Y
shear in the investigated flow results in more energetic vortices resulting from the roll-up
of the shear layer and their associated wall-pressure signature.
3.3. Autocorrelation
In order to analyze the time scales of the flow structures producing the most
energetic pressure fluctuation, the autocorrelation of p' time records was evaluated. The
autocorrelation, Rp'p', is defined as:
p' 0.0 p' (m - r)
(p'...(x))2 ‘3'”
where the overbar denotes time averaging and r is a time delay. The discrete form of
Rp.p.(t , x) =
equation (3.1) that is implemented in the processing of the discrete-time data is given by:
N—m-l
Zp'1x.n)p'(x.n — m)
R,.,.(m,x)= "=0 N (3.2)
N — lml '2
(X)
“2:313
where N is the total number of data points in the time series and m is the time delay in
data samples. The corresponding time delay in seconds is estimated from r = m/fs, where
fS is the sampling frequency in samples/sec. To calculate the autocorrelation, each
pressure data series, is split into records of 4096 samples each. Since N = 223 = 8,388,608
and 220 = 1,048,576 samples for Re = 7600 and 15700, respectively, 2048 and 256 data
records are obtained for the low and high Reynolds numbers, respectively. The
autocorrelation is finally calculated by averaging the inverse Fast Fourier Transforms
(FFT) of the product of the F FT and its conjugate for all records.
82
Figure 3.4 shows the autocorrelation coefficient (Rp'p') plotted versus the non-
dimensional time delay (tUoo/Xr) for the two Reynolds numbers examined here at five
sparse locations covering the streamwise range of the measurements: E, = 0.67, 1.11, 1.56,
2.0 and 2.33. Generally speaking, there is very little change in Rp'p' for all x positions.
This is more evident in the color contour maps in Figure 3.5. These maps
represent R W values for all x positions using color contours. The constant-shade
contours are almost parallel to the 1 axis, revealing a practically “frozen” auto-correlation
function.
In Figure 3.4, a “preferred” time scale corresponding to the negative peak in Rp'p'
is suggested. This time scale represents a quasi-periodic disturbance with a period of
tUm
X
r
z 1.3 (peak-to-peak time delay). The dominance of this disturbance decays
somewhat with increasing 8,. This is evident from the small decay in the largest negative
Rp'p. value with downstream distance. On the other hand, Figure 3.6 depicts a plot
enlarging the region around the autocorrelation peak for the lower Reynolds number and
all five x locations. The figure shows an increasing curvature of the auto-correlation at
zero time delay with increasing 5,, which implies a decrease in the Taylor microscale.
This suggests that the smaller scale turbulence is increasingly contributing to the wall-
pressure fluctuation with increasing downstream distance from X,. Finally, it is noted
here that Rp'p' results seem to be affected very little by the Reynolds number (at least for
the small Re range covered here).
83
'-2T T‘-T1 0 1 T T 2 '-2TTTT-1 T 0 1T * T2
‘ ...----:.'R.e= 1579(1‘
TU/X
00 r
Figure 3.4. Autocorrelation coefficient at five different streamwise locations downstream
of the reattachment region for Re = 7600 and 15700
84
Re=7600
Re=15700
IU/X
Figure 3.5. Contour maps of the autocorrelation coefficient for all 16 microphones and
the two Reynolds numbers; Re = 7600 and 15700
85
Figure 3.7 depicts the color contour map of the autocorrelation coefficient for the
flow upstream and downstream of reattachment. The results show combination of those
from the separating/reattaching flow region obtained by Hudy et al. (2003), g < 0.3, at Re
= 8000, and the present results. The map shows that the auto correlation width, and
hence the dominant time scale of the wall-pressure-generating structures, remain
practically unchanged with increasing x downstream of i z —0.5. This suggests that the
flow structures dominating p' generation within the development zone trace their origin
upstream to the middle of the separation bubble: an observation that is consistent with the
idea proposed earlier that the vortices originating in the separating shear layer do in fact
dominate the wall-pressure signature in the non-equilibrium boundary layer within the
streamwise extent investigated here. The persistence of the shear-layer vortices
downstream of reattachment has been reported by Bradshaw and Wong (1972), Farabee
and Casarella (1986) and Ruderich and Femholz (1986).
9 g: 0.67 T
3 3 g=1.11 1
0.95 - V 5:156 *1
‘ ‘ 552.0
g=2.33
0.9:
h.
b.
m
0.35
0.8
0.75 - ‘7 , - 3 -
-0.1 -0.08 -0.06 -0.04 -0.02 0 0.02 0.04 0.06 0.08 0.1
tU/X
m I’
Figure 3.6. A comparison of the auto-correlation results at zero time delay and five
different streamwise locations downstream of the reattachment region for Re = 7600
86
z UGO/X
-0.5 0 0.5 l 1.5 2
5.
Figure 3.7. A full contour map of the p' autocorrelation coefficient for the separating
reattaching (Hudy et al. 2003) and the present recovering flows
Upstream of g = —0.5, the dominant p' disturbances posses a substantially larger
time scale, as indicated by the significantly larger width of the auto-correlation. In
addition to Hudy et al. (2003), the larger-time-scale or lower-frequency has been reported
in different investigations and attributed to shear-layer flapping; e.g., Castro and Haque
(1987), F arabee and Casarella (1986) and Lee and Sung (2002); see Chapter 1.
3.4. Power Spectra
The wall-pressure power spectrum, (Dp'pl, is used for characterizing the frequency
content of the pressure signature at the wall. CD”: is calculated by multiplying the FFT
of the wall-pressure signal by its complex conjugate and then dividing the product by the
square of the number of samples in the spectrum. To reduce the random uncertainty error,
(DWI is calculated as the average of the power spectra of individual data records obtained
87
from partitioning the pressure time series. Each record contained 4096 samples, resulting
in a spectral resolution of 12.2 Hz. The total number of records was 2048 and 256,
resulting in random uncertainty error of 2.2% and 6.3%, for the low and high Reynolds
numbers, respectively. The specific equation that is used to calculate (Dprpc
¢prp'(0= NTLZ PINE-X”) 1 (3.3)
2
i=1 N
J
where P and P. are the FFT of the wall-pressure data record and its complex conjugate,
respectively, f is the associated frequency, fs is the sampling frequency, N is the total
number of samples in each record, NJ- is the total number of records, and j is the record
index.
Non-dimensional power spectra for the low and high Reynolds numbers, at the
same five locations of the autocorrelation results in Figure 3.4, are plotted in two
different forms in Figures 3.8 and 3.9, respectively. The top graph in both figures depicts
the non-dimensional power spectrum of the wall-pressure, B (where [3 = CDp'pl/(l/Z p
Uw2)2), versus the non-dimensional frequency (f Xr/Uw) using logarithmic scale for both
the ordinate and abscissa. For both Reynolds numbers, the top spectrum shows
concentration of the pressure fluctuations at low frequency as depicted from the broad
spectral peak close to fUX' = 0.7. This frequency value is consistent with that
documented in the literature to correspond to the passage of the vortical structures
generated in the separated shear layer upstream of reattachment. For example, the top
graph in Figure 3.8 shows a comparison between the spectrum of the present data and
88
a.~§‘5 £5.50! . .a -. ..
I
>1 ha» 12.1-
10 ;
5 o
1 _ “N -7/3 slope
-4 T ~ 11.3%” \
10 ‘ ‘3‘ \\ i
T ‘I‘AI“ \\ 4
~ TTlT“;‘§ \\ 3
5 “‘ ‘41.»33‘3‘ \\
10 1 x - {nix -~\ ~ - ~ 3
‘~ 1“ 5‘ ‘ a
3 TT‘x “if“ \\ T5 Slope i
-6 \s~~‘ 3 . 'fl 33533“? “ \\\ 3
ca. 10 “x “ “7, .‘y‘ * "1
.» ‘\ i_1.‘-\ 3
3 I 3 33 3 3 A 3 s 1"134339 I.“ 3
10-7: 1 §=0 (Hudy et al. 2003) 13 3 ‘13:} 35‘ l
E . a: 0.67 1 $93 3
E" = 1.11 1 .\ 3
-8 T i . §_ 1 i \\ l
10 3:1 £4.56 3 ~ ~ ~ ~ ~ * ~ ~ ~~ 13 ‘~~.«.\* - 3
1 a: 2.0 3 ; ‘ ‘s,
_ 3 a=z_-33 e , 1 * ii“ “ i
10 =53 ~ ~~ . 3 ~ ~ ~ - . L3* eeeeee .. ~ ~43
10 10 10 10
.4
x 10
2.5.~~~ ~ 3 3 3.. ereeeee+++j
* 1
1
2. - .-1
1 1
‘ 1
‘ 1
1.51~ ‘3
1
5 ‘ 1
1 1
l3 ‘
1 1
05+ ~ «
-1 7 7 W. , k, 3 _ 1 , ,, . :0 i. .1. . 1 1...- ; l l_ _, - . ’ -‘ i i,,,.. L. i 2
10 10 10 10
f Xr/U00
Figure 3.8. Spectra of the wall-pressure fluctuation for Re = 7600: logarithmic (top) and
semi-logarithmic (bottom) plots
89
-3
10 3
-4 I
10 1" 3 ,
l 3 3
T \
. \\
10's: ‘ \ _, , 1
' \\-5 slope 1
6 \\
©- 10 \\ ‘ 1
‘ \
' '- \
’ ~11 « ' ,\
-7 T ‘IT‘T‘T"‘iTTi.\\
T 1 7 7 , T .‘~'.-:\ , ,
E ‘ §= 1.11 ; * 1 .3 ..
'8 ' 1.» 1
10 1_ E1: 156 7 1 7 7 7 7. 3 _, 7 3, 133333333 _ _ 3
é: 2'0 I 1 \\
39~ §= 2.33 3 \:.‘\.J
10 _3 - -3 1 , ~10 7 ,. 77 +1 ..3 77 -2 L- a, a2
10 10 10 10
74
x 10
2.51 - 7 ~ 3 7 m 3 g g g 3
2 .
1
’ 1
1.53 3
C3 1 3
1 ;
0.5
1
.l . ,n, . 1 1 l l . i 0 . - L n . . ‘ . . 1 TT 1 3“, “‘ -' —- . 2 1 1 . 2
10 10 10 10
f Xr/Uw
Figure 3.9. Spectra of the wall-pressure fluctuation for Re = 15700: logarithmic (top) and
semi-logarithmic (bottom) plots
90
that of Hudy et al. (2003) at reattachment. The comparison shows that the pressure-
fluctuation—peak frequency found at reattachment is similar to that found farther
downstream beneath the redeveloping boundary layer, which is consistent with the auto-
correlation results discussed earlier. Similar values of the dimensionless frequency of the
wall-pressure spectral peak were also reported by a number of studies focused on the
separation/reattachment zone e.g., Cherry et al. (1984), Farabee and Casarella (1986) and
Lee and Sung (2001).
As the vortices are advected downstream by the mean flow, they appear to decay
as evident from the attenuation in the spectrum peak with increasing 5,. In contrast, the
spectral level at the high-frequency end is maintained by the smaller-scale turbulence.
Blake (1986) showed that the high-frequency end of the wall-pressure spectrum beneath
turbulent boundary layers is characterized by a slope of —5 on a log-log plot. The slope
of -5 (within approximately 5% uncertainty) is also found here at the high-frequency end
of the spectra shown in the top plot of Figures 3.8 and 3.9, respectively. Based on
Blake's analysis, the flow structures responsible for the -5 spectral region are those
residing in the buffer layer of a fully developed turbulent boundary layer. Thus, it is
possible that the existence of a -5 region in the spectra measured here is a manifestation
of the redevelopment of the inner region of the boundary layer downstream of
reattachment; i.e., that associated with the sub-boundary layer.
On the other hand, George et al. (1984) found that the turbulent-turbulent (TT)
source terms of pressure result in the establishment of a power-low-type spectral
characteristics of p' with exponent of -7/3; i.e., corresponding to a frequency range with
slope of -7/3 on a log-log plot of the spectrum. The top graph in both Figures 3.8 and 3.9
91
depict that the slope of the wall-pressure spectra is slightly shallower than -7/3 for 1 <
fX
' < 10 beneath the redeveloping boundary layer. However, the wall-pressure
w
spectrum of Hudy et al. (2003) at reattachment does posses a slope of -7/3 within the
same range of frequencies. Furthermore, Lee and Sung (2001) observed a slope of -7/3
in the wall-pressure spectra near reattachment, and a shallower slope downstream of
reattachment (E, z 0.35). Thus, it appears that within the separation/reattachment zone,
the wall-pressure signature reflects spectral characteristics similar to those found in a
free—shear layer. As the shear layer reattaches, and the vortical structures continue to
decay with increasing x, these characteristics gradually change, ultimately disappearing
and giving way to boundary-layer—like characteristics as x approaches infinity.
The bottom graph in both Figures 3.8 and 3.9 shows a plot of the spectral
information plotted with the ordinate on represented on linear scale (a = f B X,/ U00). The
rationale for plotting the spectra in this manner may be seen when integrating the
spectrum to obtain the fluctuating wall-pressure energy:
7 co m
p = 1d>,.,.df=1f<1>,,.dnog(m (3.4)
0 0
Equation (3.4) shows that when using a logarithmic frequency axis, the geometrical area
under the spectrum curve corresponds to the pressure fluctuations energy only if (1);,er is
multiplied by f and plotted on a linear scale. The corresponding non-dimensional
quantity is then or instead of B. The abscissa still represents the non-dimensional
f X . . . . .
frequency U ' . The figures confirm that the maxrmum contribution to p' occurs Within
w
a non-dimensional frequency band centered around O.7-O.9, which corresponds to the
92
passage frequency of the separated shear layer vortices obtained in previous literature
within the separation/reattachment zone as discussed above. The center of this range,
er/Uoo = 0.8, corresponds to a non-dimensional time scale of 1.25, which agrees quite
well with the dominant time scale of 1.3 identified from the auto-correlation results
earlier. Moreover, visual extrapolation of the data towards lower frequencies than those
shown in the plot suggests that only very small fraction of the overall p' energy is not
captured because of the low cut—off frequency of the microphones (20 Hz).
The semi-log plots also verify the growth of the boundary-layer fractional
contribution to p' energy. This may be examined by considering the ratio of the area
under the spectrum curve at high frequencies to that at low frequencies (p'h/p'1)2. This
ratio is estimated by splitting the spectrum of the wall-pressure fluctuation into two
portions. One is for low frequencies (due to the large vortical structures originating in the
separated shear layer) and the second is for boundary-layer contribution at high
frequencies (or small turbulent eddies). This idea of the p' spectrum beneath the non-
equilibrium flow being a composite of two spectra separated in frequency, one associated
with the outer-shear layer and the other with the sub-boundary layer, may be examined
further in Figure 3.10. The figure depicts a comparison of typical p' spectra in
separating/reattaching (Hudy et al. 2003) and equilibrium boundary layer (Gravante et al.
1995) flows. To superpose the spectrum from Gravante et a1. (1995) and the present
spectrum at g = 2.33 on the same plot, the start of the -5 region of the two spectra were
made to overlap3. On the other hand, the spectrum of Hudy et al. (2003) was re-scaled
3 As a side note: the equilibrium-boundary-layer literature shows the -5 region to start at f v/u,2 z 0.13,
where v is the kinematic viscosity and n1 is the friction velocity. Using this value suggests that n1 for the
current flow is 0.5 m/s
93
such that its peak coincides with the present spectrum to facilitate comparison of the
spectral shapes.
11 — geiflkeattEchfi = 0,1Hud_y et al.) 3
11 —- Equilibrium B.L. (Gravante et al.)
-9 3 — Present Data (2'; = 2.33) 1
10 1.. M- *4 :1 17? Jr 1 7 V ,, ;‘ f
71
0
10 10
1
1
1
1.
1
— 1
.—l
O
‘ I I .1
—O
O
H 1-
O
fX/U
1' 00
Figure 3.10. Definition of the frequency used in the splitting procedure of p' spectrum
Based on the comparison in Figure 3.10, the splitting of the spectral energy into
out-shear—layer and sub-boundary-layer components was based on the frequency at which
the spectrum changes its slope, which is indicated by fn in Figure 3.10. The choice of fn
is somewhat heuristic since the demarcation between the two frequency ranges is likely
to be not sharply defined, and one would expect a gradual switch from one spectrum to
the other with the middle range of frequencies receiving contribution from both the
smallest scales of the outer-shear layer and the largest of the sub-boundary layer.
94
. leI In Ila? or”. J... .. P1VWWU
P5055 .. .
. .
However, fn does coincide with the point at which the spectrum slope switches from the
value typical of the middle frequency range in turbulent boundary layers (z —1 , see Blake,
1986) to a more negative value that has never been observed in boundary layers at low
frequencies.
The streamwise variation of the ratio (p’h/p'l) 2 is provided in Figure 3.11 for the
two Reynolds numbers examined here. The results verify that the boundary layer
contribution to the wall-pressure fluctuation becomes increasingly important with
increasing x. Moreover, Figure 3.11 shows that the initial contribution of the small-scale
turbulence to p' energy increases as Reynolds number increases, which is consistent with
known turbulence physics. Note that the two solid curves in the figure are plotted to aid
in visualizing the trend of the data.
0.4 :1
6 Re =7600 ’1 '
0.351 ,',,' Re :‘5700 ;
0.3
0.251
0.21
(rah/p.)2
0.15-
0.1 "
0.05 ~ 1 ~ ~ 71
11:5 ‘ 1 " 1:5 ' i s 7.5
Figure 3.11. The ratio between high- and low-frequency pressure fluctuations energy
95
3.5. Cross Correlation
The cross-correlation coefficient (va3,!) is defined as:
C : p' , and conditionally-averaged mean-removed velocity-
2.05
vector field associated with positive (top plot) and negative (bottom plot) wall-pressure
Figure 4.35. Wall
peaks at g
values at I)?” z 0.3 and -0.5. In the flow over blunt flat plate, Kiya et al. (1982) and
f
Kiya & Sasaki (1985) reported positive u' fluctuation near the wall between successive
vortices. They attributed the instantaneous excess in streamwise momentum near the
wall to the penetration of the irrotational freestream fluid close to the wall in between the
vortical structures. They also associated this region of the flow with the formation of
saddle point, when viewed in a frame of reference that is translating with the convection
velocity of the vortices.
The bottom plot in Figure 4.35 shows the conditional-velocity-vector field
associated with negative wall-pressure peaks. As seen in the plot, the negative wall-
pressure fluctuations seem to coincide with negative u' fluctuation. This was also
observed by Kiya et al. (1982) and Kiya & Sasaki (1985) who also indicated that the
negative pressure coincided with the passage of a spanwise-vortex center, which is a
significant source of low pressure.
The above discussion leads to the following simple scenario for the successive
generation of positive and negative wall-pressure peaks. As the large-scale vortical
structures of the outer-shear layer travel downstream, they entrain irrotational freestream
fluid in between each pair of successive vortices. Accordingly, the accompanying
interface between the irrotational flow in the freestream and the vortical structures moves
up and down in the direction normal to the wall in synchronization with the passage of
the vortices above a certain streamwise location. Specifically, when a vortex passes over
a point of observation, it thickens the boundary layer leading to a streamwise momentum
deficit across the boundary layer and associated negative wall pressure, as found above.
On the other hand, the entrained irrotational flow between the vortices reduces the
162
boundary layer thickness and provides streamwise momentum excess. This coincides
with the saddle points up and downstream of each vortical structure and is associated
with the generation of positive wall pressure.
The x-momentum deficit and excess accompanying the passage of the vortex
center and saddle point, respectively, are examined in Figure 4.36. The figure depicts a
comparison between the mean streamwise-velocity profile and the conditionally averaged
streamwise-velocity profiles at the instant of occurrence of the positive and negative p'
peaks. The plot shows that the conditionally averaged streamwise-velocity profile
associated with the positive wall-pressure peak has higher velocity at all heights included
in the present measurements, which confirms the fact that positive pressure coincides
with x-momentum excess across the boundary layer. On the other hand, the conditionally
averaged profile synchronized with the negative wall-pressure peaks drops below the
mean-velocity profile at all heights. The above results clarify some of the features of the
spatial correlation results in Figure 4.32. In particular, the top plot in the figure shows a
vertical positive Rup correlation lobe that extends across the entire y range at a g location
of 2.05. Based on the analysis of the conditional streamwise-velocity profiles it is now
clear that positive p' is associated with positive u' across the whole y range and vice
versa, thus resulting in the vertical positive correlation lobe seen in Figure 4.32. The
corresponding feature in the Rvp correlation (bottom of Figure 4.32) is a strong
streamwise gradient of the correlation value that is associated with a switch in the
correlation sign from positive to negative. This would be consistent with the scenario
identified above that associates the pressure generation with vortex centers and saddle
points. Both of these features are characterized with a switch in the direction of the v'
163
component from the upstream side of the source to the downstream one. In the
following, the association of vortex centers and saddle points with wall-pressure peaks
will be established more convincingly.
0.7' 7 ., a, 7:7 7 7I7 7 *2i ’7' "':7.0
A Positivepeak ‘ ‘ 7':
V NegativePeak; v ‘
0.6179 rMeanprofilei ; .~ ~ . v'.? 7 7 16.0
1 V .A
1 V CA
0.51 ~ ' 1 , 7 7.7-..‘77 ' ~15.0
.V .A .
V .A
‘ V .A 1
04.7 :- 0.; 114.0
Q A 5
>1 '7 ..A >\
0.31 *V O A , w , 13.0
V O A
V O A
V O A ‘
1 ‘ 'V.O AA
0.1" 7 , ”V ‘7‘. 111.0
x r ‘ 1 1 1
V
(61.7 7 7 7 .1 7 7 . 1K. #7. -77 ..‘0
.5 0.6 0.7 0.8 0.9 1
u/U0°
Figure 4.36. Comparison between the mean streamwise-velocity profile and the
conditionally-averaged streamwise-velocity profiles associated with positive and negative
wall-pressure peaks
Although Figure 4.35 illustrates the main features of the conditional velocity field
and associated wall-pressure signature, it does not show the real physical picture of the
flow structures. To obtain such a picture, the velocity vector field has to be viewed from
a coordinate frame that moves with the convection velocity of the flow structures (Fiedler
1988). Fielder (1988) investigated this issue extensively and referred to the vector-field
representation in Figrre 4.35 as triple-decomposition, which is obtained by looking at the
mean-removed velocity vector field in laboratory coordinates. However, he argued that
164
in order to obtain the velocity vector field that represents the real physical picture of the
flow structures, one has to do a double-decomposition of the flow by viewing the total
velocity vector field (mean + fluctuation) relative to an observer moving with the
convection velocity of the flow structures. The top and the bottom plots of Figure 4.37
show the same results as in Figure 4.35 after adding the mean velocity and viewing the
result in a frame of reference that is moving with a convection speed of 0.81Uoo. The top
plot of the figure shows clear evidence of a saddle point that is synchronized with the
positive wall-pressure peak directly below it, while the bottom plot provides a clear
physical picture of a vortical structure with a center located directly above the negative
wall-pressure peak. This provides a strong support for the central role played by the
vortex structures and saddle points in the scenario of wall-pressure generation outlined
earlier. Note that strictly speaking, because the abscissa in Figure 4.37 represents time
rather than space, the results are indicative of the passage of the vortex structures and
saddle points past the point of observation of the pressure, rather than the actual spatial
structure of these features. However, as discussed in Chapter 3, it is believed that Taylor
hypothesis of frozen turbulence holds fairly well for the large-scale vortices, and hence
the results in Figure 4.37 are likely to represent well the actual spatial structure. A more
quantitative assessment of this statement is provided in 4.3.
4.3. Stochastic Estimation
The above analysis has led to the identification of the quasi-periodic passage of
the outer-shear layer vortices and saddle points resulting from their mutual interaction as
the primary source of wall-pressure generation. Generally, the identification is based on
165
'— 6.75
77
77.7.7. 7 m
7 7- 7 .
77,1 . 7 _
7,, . 77
7,7,7”, §
77/47” \3
7,
.=‘ €/
:4
I
7‘
.17"—
{/1
"'
_" I
73
Ebb?
735 7.
0.4
0.2
7 0.1
L011 7 -03 7—012
-0.5
_27 7.
0.1
tU/X
0.3 0.4 0.5
0.2
—o.2’ '
—o'.3
—o.4
71.5
and conditionally-averaged velocity-vector field
Figure 4.37. Wall-pressure, ,
-pressure peaks, viewed
associated with positive (top plot) and negative (bottom plot) wall
in a frame of reference moving with 0.81Um
166
time-averaged statistical information, which does not reveal an in-depth picture of the
variation in the characteristics of the individual flow structures or the nature of the
mechanism(s) leading to the wall-pressure generation. Hence, in this section an analysis
is conducted using the stochastic estimation of the flow field based on the instantaneous
p' signature across the ten most downstream microphones. As explained below, the
procedure yields an estimate of the "quasi-instantaneous" velocity over the full y range of
the measurements. Subsequent evolution in time of this estimate yields a picture of the
quasi-instantaneous flow structures traveling past the streamwise location of the hotwire,
which can be examined to reveal some of the variability in the individual flow structures
and associated wall-pressure-producing mechanisms.
4.3.1. Principle and Equations
Application of the stochastic estimation technique for identification of
organized, or coherent, motion in turbulent flows was first proposed by Adrian (1975)
and then fiirther refined by him and his coworkers (e.g. Tung and Adrian 1980). Since
then, applications of stochastic estimation have expanded to range from the identification
of the conditional eddies of isotropic turbulence, Adrian (1979) and Adrian & Moin
(1988), to the extraction of the spatio-temporal flow field characteristics of turbulent
boundary layer structures, Guezennec (1989) and Choi & Guezennec (1990).
Classically, the stochastic estimation technique estimates any variable at one location in
the flow from a known variable at some other location. However, and as pointed out in
the objectives section, stochastic estimation is used in the present work to estimate the
velocity field above the wall from its wall-pressure signature. In particular, the technique
seeks to get the estimated velocity ufi'sfij, + 'r‘,t + r) from a known wall-pressure p'(‘r‘o ,t);
167
where f is the position vector of the pressure-observation point, the components of
0
which are x0, 0 and 20, in the streamwise, wall-normal and Spanwise direction,
respectively, f (x-xo, y, z-zo) and r are the spatial-offset vector and time—shift between the
estimated velocity and wall-pressure event, respectively, subscript i is the tensor index
notation specifying the velocity component and subscript s is used to denote the
stochastic estimate. Note that “is“; +"r’,t+r) and p'(fo,t) are the mean-removed
values. The former is obtained from a Taylor-series expansion in terms of the latter as
follows:
ug,$(fo + f, t + r) = Ai('r‘o + it) p'(?0,t) + Bi('r‘o + it) p'2('fo,t) + - -- (4.5)
where, A, (f0 + f, r) and BIGo + 'f, r) are the estimation coefficients for the linear and
quadratic terms in the expansion, respectively. Note that both coefficients depend on the
location of observation of the event (to) but not its time of occurrence (t). This is
because the flow field considered is generally inhomogeneous and hence statistical
properties of the flow in the vicinity of a particular event location (1;) will differ from
those at another location. On the other hand, the flow is stationary in time, and therefore
its statistical properties are independent of the time of occurrence of the event.
The primary goal of the stochastic estimation procedure is to obtain the
coefficients A, and B,, so that the velocity field can be estimated from equation 4.5. If
only the linear term in the equation is included in the estimation, it is known as Linear
Stochastic Estimation (LSE). However, if the first two terms are included, the estimation
is a Quadratic Stochastic Estimation (QSE). Note that both LSE and QSE can be
implemented using wall-pressure conditions from multiple points in space
168
simultaneously, yielding a "multi-point" estimation. The multi-point estimate is
generally more effective in capturing the "true" nature of the flow structures, and it is the
one of primary focus here given the availability of multi-point information from the
pressure-array data. The multi-point estimate will only be of the LSE type because of the
complexity involved in obtaining a multi-point QSE. However, the results will be
compared to those from single-point LSE and QSE. Derivation of the equations used to
obtain the estimation coefficients for all of these cases is given below.
I. Single-Point LSE
As mentioned earlier, the single-point LSE (denoted by SL hereafter), provides a
linear estimation of the conditional flow filed, u;(?0 +'f,t + r) from the wall-pressure,
p'(fo,t). For brevity, the following analysis will only consider the x~component of the
velocity, but the resulting equations are the same for the y-component. The estimated u'
is calculated from the following equation:
u'.,(fo + it + T) : Au,lin(-fo + it) p'(ifi) (4-6)
where, Au, 1m is the coefficient of the linear estimate of u'. This is determined such that
the long-time mean squared error between the velocity and its estimate, e("fO +f) , is
minimized. This is the maximum likehood estimate, assuming a Gaussian distribution of
the error about the estimate. Explicitly, at a specific location within the flow domain the
mean squared error is:
e(f0 +r) =(u',(ro +f,t+r)—u'('fo +f,t+r))2
(4.7)
: (Au. lin p'(ro ) tC) _ U'Cfo 'l' f, t ‘l" 1))2
where, the overbar denotes time averaging. The value of Aqun that minimizes the error
may be found then by requiring:
169
de d ' ,_ ,1 _
dA =dA KAN.p(r.,t)-U(r.+r,t+t))’l=0 (4.8)
u. lin u. lin
Exchanging the order of the derivative and time integration operations:
d
dA
[(AuJin FIG), 0 ‘ “'60 + it + T))2l
u. lin
.—. 205...“. p'2(fo, t) — u'(f0 + i‘,t + r) x p'(Fo,t)) (4.9)
= (Au.lin p'2 (to, t) — u'('r; + 'r’, t + t)p'(i';,, t))= 0
Equation 4.9 is satisfied if
_ u'(fo + 131 + r) p'(fo,t) _ u'(ro + ‘r‘,t + r) p'(ro,t)
Au.lin — ,2 _ :2 -
p (r..t) pm..(r.,t)
(4.10)
Inspection of Equation 4.10 shows that Au,“ equals the ratio between the cross-
correlation of u' and p', at a time delay corresponding to the delay between the estimate
and the event, and the square of the root-mean-square, or second PDF moment, of the
wall-pressure. Equations 4.6 and 4.10 yield the final form of the equations employed for
SL estimation:
u;=["(r°+r.’.t+..t)p(r°’°]p' & v;=Iv(r.+r',2t+_r)p(r..t)Ip. (4.11)
pm(1;,,t) prms(ro’t)
II. Single-Point QSE
Similar to the single-point LSE, SL, the single-point QSE will be denoted by SQ.
Based on equation 4.5, the main equation for obtaining the SQ estimate is given by:
u'S(fo + 'r',t + r) = Au,quad ('r‘o + 'r‘, r) p'(‘r’o, t) + Bu(‘r‘o + 'r‘, r) p'2 (1;, t) (4.12)
where, Amquad and Bu are the unknown coefficients of the first- and second-order terms
for the SQ method. Similar to SL, those coefficients are determined by a least square
error minimization process:
170
e (to + r) = (AWad p'(ro, t) + Bup'2(ro,t) — u'(i"o + r, t + o)’ (4.13)
By taking the derivative of Equation 4.13 with respect to Amquad and Bu individually and
equating the result of each to zero, the error is minimized in a least-square sense:
dACrlr:uad = dquuad [(Au quad p’(izo ’ t) + Bup’2 (f0, 0 — u’(-f0 + f’t + 1))2]
= d [(All quad p'(‘g, t) + Bup'2 (to, t) — u'(f0 + f, t + 10):] (4.14)
dAu.quad ‘
= ZlA...q..d 13' +3.1)” -U')>< p' = (AM... 137’— +3.? -U'—p’)=0
548—6: : 21%[(Au‘quad p'(fo9t) + Bup'2(i:o’t) - u'(i:o +f,t+t))z:l
d y - l "’ ' ‘ -.
236-— (AILQuadp(roat)+Bup2(ro’t)—u(r0+r’t+1))2] (4'15)
= 2(1chl p'(fo,t) + Bup'2(r;,,t) — u'('fo + 'f,t+ 1))x p'2(ro,t)
= (Au.quad p'3('r;,,t) + Bup"’(fo,t) — u'(fo + ?,t+'c) p'2(‘r‘°,t))= 0
Solving Equations 4.14 and 4.15 for Au, quad and B, gives:
B _ p'2(?o, t) u'(f0 + it + 1:) p'2 (f0, t) — p'3('r'o, t) u'('r; + it + r)p'(fo, t)
u (4.16)
(Memo)? —p'“ &