13.11.2014 Views

Integrating CFD and Experiment in Aerodynamics - CFD4Aircraft

Integrating CFD and Experiment in Aerodynamics - CFD4Aircraft

Integrating CFD and Experiment in Aerodynamics - CFD4Aircraft

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

2.2 Data Aquired<br />

The Embedded Laser Doppler Velocimeter (ELDV) has an optical head mounted on a support<strong>in</strong>g turntable l<strong>in</strong>ked to<br />

the oscillat<strong>in</strong>g frame as sketched <strong>in</strong> figures 1 <strong>and</strong> 2. The optical head is equipped with a beam-exp<strong>and</strong>er <strong>in</strong> order to<br />

<strong>in</strong>crease the focal distance up to 400 mm. The laser beams focus through a 45 deg mirror at a given position <strong>in</strong> chord<br />

<strong>and</strong> span. The support<strong>in</strong>g turntable is l<strong>in</strong>ked with the oscillat<strong>in</strong>g frame, so that U <strong>and</strong> V velocity components can be<br />

directly measured <strong>in</strong> the same reference frame of the motion. Due to the periodicity of the motion, each period is<br />

considered as a specific sample of the same phenomenon. So, each velocity component is recorded at each phase<br />

angle ωt between 0 deg <strong>and</strong> 360 deg by steps of 0.703 deg over a large number of periods. Data are then statistically<br />

analyzed at prescribed values of the period, e.g. the <strong>in</strong>stantaneous <strong>in</strong>cidence, with an uncerta<strong>in</strong>ty of:<br />

δα=4∆α/360=24/360=0.066deg (1)<br />

Data acquisitions are made on a microcomputer from two Burst Spectrum Analyzer (BSA) deliver<strong>in</strong>g the values of the<br />

two components U <strong>and</strong> V together with the arrival validation time for each frequency measurement. The software used for<br />

acquisition <strong>and</strong> data reduction (COMBSA) was developed at LABM under the Apple LABVIEW system.<br />

The unsteady data reduction is performed us<strong>in</strong>g an ensemble average procedure suited for periodic flows <strong>in</strong>vestigation<br />

as described below. Dur<strong>in</strong>g a boundary-layer survey (s/C fixed) <strong>and</strong> for each normal distance y above the wall, 38 blocks of<br />

512 values are acquired which correspond to 19456 measurements of U <strong>and</strong> t u (arrival validation time of U) <strong>and</strong> the same<br />

number for V <strong>and</strong> t v. The two BSA work <strong>in</strong> a master-master mode : the (U,V) acquisition is <strong>in</strong>dependently performed with the<br />

same time orig<strong>in</strong>. Moreover, this particular acquisition mode is coupled with a «dead time» function which makes sure that<br />

only one particle can be validated dur<strong>in</strong>g a fixed time <strong>in</strong>terval (here ∆t = 10 -4 s). This function avoids the simultaneous<br />

validation of a burst of particles. Based on such options, the oscillation cycle number varies <strong>in</strong> a range between 150 <strong>and</strong> 300<br />

depend<strong>in</strong>g on the flow configuration (with or without separation). F<strong>in</strong>ally, one data file conta<strong>in</strong><strong>in</strong>g the 19456 data l<strong>in</strong>es (t U, U,<br />

t V, V) is created for each normal distance y <strong>and</strong> the analysis software, provides the phase averaged (,) velocity<br />

components <strong>and</strong> their associated RMS quantities 1 . Indeed, due to the periodicity of the flow, each period is considered as a<br />

specific sample of the same phenomenon, so that each velocity component can be obta<strong>in</strong>ed at each phase angle t as the<br />

averaged value of the velocity samples recorded at the same given phase angle <strong>and</strong> over a large number of oscillation cycles<br />

(greater than 150).<br />

The negative arrival times due to the <strong>in</strong>ternal clock reset of the BSA are cut out <strong>in</strong> the first step. In the second step, the<br />

oscillation cycles are counted to keep only the common Nc cycles correspond<strong>in</strong>g to the U <strong>and</strong> V simultaneous<br />

measurements. In the third step, the Nc cycles are pooled to obta<strong>in</strong> only one fictitious cycle. This cycle is split up <strong>in</strong>to different<br />

divisions (512 <strong>in</strong> the present study) where velocities are averaged <strong>in</strong> order to provide discrete equidistributed representations<br />

of U <strong>and</strong> V velocities, written as <strong>and</strong> . The discrete representation obta<strong>in</strong>ed (512 po<strong>in</strong>ts) results for each division,<br />

from an average of measurements from one or different cycles. Indeed, the division length (which represents 2.10 -3 second<br />

for a typical frequency of 1Hz), appears to be larger than the «dead time» <strong>in</strong>terval fixed to ∆t = 10 -4 s.<br />

The next step is to produce a Fourier series (256 harmonics) represent<strong>in</strong>g the two mean functions of velocity<br />

components (written as U ~ <strong>and</strong> V ~ ) <strong>in</strong> pitch<strong>in</strong>g motion case. The fluctuat<strong>in</strong>g quantities u’ <strong>and</strong> v’ can be then calculated by :<br />

u' n (t) = U n (t) -U(ωt), v' n (t) = V n (t) -V(ωt) .<br />

Then the RMS <strong>in</strong>tensities are def<strong>in</strong>ed by:<br />

σ u(<br />

t)<br />

=<br />

σ v(<br />

t)<br />

=<br />

1<br />

Nu<br />

1<br />

Nv<br />

Nu<br />

[ un' ( t)<br />

]<br />

n=<br />

1<br />

Nv<br />

[ vn' ( t)<br />

]<br />

n=<br />

1<br />

2<br />

2<br />

(2a)<br />

(2b)<br />

Where Nu <strong>and</strong> Nv are respectively the number of U <strong>and</strong> V velocity component values <strong>in</strong>cluded <strong>in</strong> one division.<br />

All these fluctuat<strong>in</strong>g quantities are presented accord<strong>in</strong>g to the phase averaged procedure, by represent<strong>in</strong>g the , ,<br />

, , , <strong>and</strong> <<br />

u 'v'<br />

> quantities as a function of the normal distance y to the surface.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!