Separation of the turbulent boundary layer (BL) on a flat plate under adverse pressure gradient was studied experimentally using Time-Resolved PIV technique. The results of spatio-temporal analysis of flow-field in the separation zone are presented. For this purpose, the POD (Proper Orthogonal Decomposition) and its extension BOD (Bi-Orthogonal Decomposition) techniques are applied as well as dynamical approach based on POPs (Principal Oscillation Patterns) method.
The study contributes to understanding physical mechanisms of a boundary layer separation process. The acquired information could be used to improve strategies of a boundary layer separation control.V článku jsou prezentovány výsledky prostorově-časové analýzy experimentů popisujících odtržení turbulentní mezní vrstvy na rovinné stěně při nepříznivém tlakovém gradientu. Pro měření byla použita rychlá metoda PIV, pro analýzu potom vlastní ortogonální dekompozice (POD) a její zobecnění biortogonální dekompozice (BOD). Dále pro výzkum dynamických vlastností byla aplikována metoda hlavních oscilačních struktur (POPs).
Studie přispívá k pochopení fyzikálních mechanismů, které se uplatňují v procesu odtržení mezní vrstvy. Získané poznatky mohou být uplatněny při vylepšování strategie řízení odtržení mezní vrstvy

Uruba, Václav

Transactions of the VŠB - Technical University of Ostrava Mechanical Series

DSpace at VSB Technical University of Ostrava

Transactions of the VŠB – Technical University of Ostrava, Mechanical Series No. 3, 2010, vol. LVI article No. 1830 Václav URUBA* DECOMPOSITION METHODS FOR A PIV DATA ANALYSIS WITH APPLICATION TO A BOUNDARY LAYER SEPARATION DYNAMICS DEKOMPOZICE DAT ZÍSKANÝCH METODOU PIV S APLIKACÍ NA VÝZKUM DYNAMIKY ODTRŽENÍ MEZNÍ VRSTVY Abstract Separation of the turbulent boundary layer (BL) on a flat plate under adverse pressure gradient was studied experimentally using Time-Resolved PIV technique. The results of spatio-temporal analysis of flow-field in the separation zone are presented. For this purpose, the POD (Proper Orthogonal Decomposition) and its extension BOD (Bi-Orthogonal Decomposition) techniques are applied as well as dynamical approach based on POPs (Principal Oscillation Patterns) method. The study contributes to understanding physical mechanisms of a boundary layer separation process. The acquired information could be used to improve strategies of a boundary layer separation control. Abstrakt V článku jsou prezentovány výsledky prostorově-časové analýzy experimentů popisujících odtržení turbulentní mezní vrstvy na rovinné stěně při nepříznivém tlakovém gradientu. Pro měření byla použita rychlá metoda PIV, pro analýzu potom vlastní ortogonální dekompozice (POD) a její zobecnění biortogonální dekompozice (BOD). Dále pro výzkum dynamických vlastností byla aplikována metoda hlavních oscilačních struktur (POPs). Studie přispívá k pochopení fyzikálních mechanismů, které se uplatňují v procesu odtržení mezní vrstvy. Získané poznatky mohou být uplatněny při vylepšování strategie řízení odtržení mezní vrstvy.  1 INTRODUCTION Classical concept of a BL separation starts from the Prandtl’s idea of a BL and is built for 2D cases. The necessary condition for BL separation is the increasing pressure in the streamwise direction, i.e. positive (or adverse) pressure gradient along the flow path. The necessary second determining factor is presence of viscous effects in the BL no matter being of laminar or turbulent in nature. Above a critical Reynolds number, the flow field is shown to undergo self-sustained 2D low-frequency fluctuations in the upstream region of the separation, resulting into aperiodical vortex shedding further downstream. Low-frequency fluctuations, also called “flapping”, have been shown to be a characteristic feature of separated layers in general ([3]). It has been argued that they are due to a global instability manifested in the reattachment region ([4]), triggered by topological flow changes generating secondary recirculation zones ([2]).                                                                                                                                                                    * doc., Ing., CSc., Institute of Thermomechanics AS CR, v.v.i., Fluid Dynamics Department, Dolejškova 5, Praha 8, phone (+420) 286 588 547, e-mail uruba@it.cas.cz  157 2 APPROACHES Apparently the phenomenon in question is substantially 3D and highly dynamical in nature in the same time. To analyze such a phenomenon, appropriate methods of analysis should be applied on experimental data. Actually no universal methods for such analysis are available. The main problem is a huge number of degrees of freedom of the dynamical system, which should be reduced substantially. For this purpose the decomposition methods based on energetic (POD) and stability (POPs) analysis are to be used. Reduction of the dynamical system dimension is to be done using the BOD method. This approach has been proposed and tested in [11].  2.1 Methods Based on POD The Proper Orthogonal Decomposition (POD) method has applications in almost any scientific field where extended dynamical systems are involved. Recently, the POD has been widely used in studies of turbulence. Historically, it was introduced in the context of turbulence by Lumley [6] as an objective definition of what was previously called big eddies and which is now widely known as coherent structures. Details on this classical method could be found in [7]. The next step is the Bi-Orthogonal Decomposition (BOD) representing itself an extension of the POD. While POD analyses data in spatial domain only, the BOD performs spatiotemporal decomposition. Aubry [1] presented the BOD as a deterministic analysis tool for complex spatiotemporal signals. First, a complete two-dimensional decomposition was performed. These decompositions were based on two-point temporal and spatial velocity correlations. A set of orthogonal spatial (“topos”) and temporal (“chronos”) eigenmodes are to be computed to allow the expansion of the velocity field.  2.2 Stability Analysis – POPs The Principal Oscillation Patterns (POPs) decomposition is based on modal structures based on temporal or spatial linear evolution dynamics. Each mode is characterized by a complex frequency involving information on frequency, phase and growth/decay information. On the other hand, the POPs analysis is unable to resolve standing oscillations. The basis of the POPs analysis was formulated by Hasselmann [5] for discrete Markov processes in linearized dynamical systems driven by white noise with application in climatology. The POPs theory is a special case of more general Principal Interaction Patterns (PIPs) method for nonlinear dynamical systems. In the POPs approach the fluctuating part of Navier-Stokes equation is modeled by Langevin equation for the linear Markov process. The method provides the complex modes characterized with frequency and so called e-fold time e which characterizes the delay for the mode amplitude reduction in ratio 1 : e. The least stable modes are pointed out. Detailed description of this method is given in [9].  3 EXPERIMENT The existing blow-down wind tunnel was used with one-sided diffuser angle 15.5 deg. Mean flow velocity in the diffuser inlet was 7.8 m/s, turbulent boundary layer on the bottom wall. In Fig. 1 the schema of experimental setup is shown. The coordinate system was defined with the x-axis in the input flow direction and y-axis is perpendicular to the wall. Origin of the coordinates is located in the beginning of the diffuser; the cross-section here is 100 x 100 mm2. Downstream of this section, the upper wall is inclined, while the bottom plane-wall is used to study the boundary layer separation. To prevent separation from the upper wall this is permeable and sucked out. The bottom-wall boundary layer was of turbulent nature, about 5 mm thick. The suction velocity along the upper wall could be estimated to 3 m/s.  158 Fig. 1 Experimental setup Standard configuration of time resolved PIV was used for measurements, particles were generated using fog SAFEX generator, introduced to the inlet of the blow-down facility. Measurements take place in the plane of symmetry of the channel, two velocity components in this plane are evaluated in grid 79 x 18 points, 5108 subsequent complete vector fields are evaluated with frequency 1545 Hz representing approx. 3.2 s. Details on used measuring techniques are given in Uruba, Knob [8].  4 RESULTS A few preliminary results have been already obtained and presented in [10]. Here classical evaluation of “separation point” position is shown together with the first information on flow dynamics in the separation zone. The evaluation of separation point as a place of zero mean wall shear stress, or zero mean streamwise velocity component close to the wall is shown in Fig. 2. The line  represents locus of zero mean streamwise velocity. 0lxy350 400 450 500 55002040l0D Fig. 2 Mean streamwise velocity distribution The velocity fluctuation could be recognized in the whole separation region, but maximum fluctuation activity have been detected in the region above the  line - see Fig. 3. 0lxy350 400 450 500 55002040D: 0.1 0.3 0.5 0.7 0.9 1.1 1.3 1.5 1.7 1.9 2.1 Fig. 3 Velocity fluctuation energy distribution The first results of POD and POPs analysis will be shown. As an example in Fig. 4 the first 4 most energetic toposes are presented. The toposes consist of vortical structures, lower modes with  159high energy of one or two of them, while higher order modes 10 and 25 of 5 or even 10 small structures respectively. xy350 400 450 500 55002040Top 1xy350 400 450 500 55002040Top 2xy350 400 450 500 55002040Top 3xy350 400 450 500 55002040Top 4  Fig. 4 Toposes No. 1, 2, 3, 4, 10 and 25 Then, the POPs analysis is to be carried out after the system truncation. The system dynamics is modeled by 30 most energetic POD modes representing more than 82 % of the total fluctuating energy. Re β [rad/s]Imβ[rad/s]-80 -60 -40 -20 0-400-2000200400 Fig. 5 The POPs eigenvalues The POPs analysis shows the less stable modes. Fig. 5 shows complex eigenvalues of the 30 evaluated POPs modes. The eigenvalues are either real or complex conjugate. The real part of a given eigenvalue   is connected with stability of the mode and value of the e-fold time e , while imaginary part determines the frequency f :  160   Im1 ,Re 2ef     .  All modes are stable, as the real parts of eigenvalues are negative and thus the e-fold time is positive. Then, the modes are ordered according to theirs stability, less stable first (i.e.   closer to 0 and e  bigger). Tab. 1 First 10 modes Mode No. f [Hz] e [ms] n = e/T 1 0 330 02 0.8 138 0.123 1.8 132 0.244 3.2 132 0.425 9.1  65 0.59 6 10.6  47 0.50 7 0  42 0 8 12.2  41 0.50 9 6.3  39 0.25 10 2.4 37 0.12 In Tab. 1 there are first 10 modes with corresponding frequencies f , e-fold times e and ratios . The ratio  shows relationship of e-fold time n n e and the period 1T f  and could be expressed as  /e en f T   . Its physical meaning is number of periods while the mode amplitude decays by factor e . The given mode could be considered as periodic only when the decay during one period is not very big, let say not more than factor 10. This situation corresponds to 0.43n  , for  bigger then this value the mode behaves periodically, because we could positively register more than one period during the mode single appearance. For  we have nearly aperiodical mode behavior as it disappears before the first period is accomplished. In the case n0.43n 0f n   we have completely aperiodical behavior with proportional decay, when 0 0.n 43   than the behavior is nearly aperiodical. In our case we have proportionally decaying modes 1 and 7, nearly aperiodical modes 2, 3, 4, 9, 10 and periodic modes 5, 6, 8. xy350 400 450 500 5500204060POPs 1real Fig. 6 Real part of the least stable POPs 1 mode Now we show topology of some typical modes. The first mode No. 1 is proportionally decaying with the imaginary part vanishing, its real part is shown in Fig. 6, isolines map the vorticity (negative values by dashed line) showing position of vortices. The vortices are distributed more or  161less equally within the whole separation region. During the mode occurrence its structure does not change, it only decays. xy350 400 450 500 5500204060POPs 2realxy350 400 450 500 5500204060POPs 2imaginary Fig. 7 Real and imaginary parts of the nearly aperiodical POPs mode No. 2 The nearly aperiodical mode No. 2 is shown in Fig. 7. The vortical structures are distributed along the  line. They move slightly during the mode decay. 0lxy350 400 450 500 5500204060POPs 5realxy350 400 450 500 5500204060POPs 5imaginary Fig. 8 Real and imaginary parts of the periodic POPs mode No. 5 The example of periodic mode is in Fig. 8, where the mode No. 5 is shown. This mode exhibits maximal vorticity concentration along the  line, where system of traveling contra-rotational vortex pairs could be identified. 0lThe other evaluated modes manifest qualitatively the same type of behavior however the detailed topology is different, of course.  5 CONCLUSIONS The stability modes have been studied on the reduced dynamical model of a boundary layer separation region. The pulsating and traveling structures associated with flapping process have been identified. Traveling structures are concentrated close to the zero streamwise mean velocity component line, while pulsating structures are distributed in the whole separation region. ACKNOWLEDGEMENTS The author gratefully acknowledges financial support of the Grant Agency of the Czech Republic, projects No. 101/08/1112 and P101/10/1230.         162REFERENCES [1] AUBRY, N., GUYONNET, R., LIMA, R., 1991, “Spatiotemporal Analysis of Complex Signals: Theory and Applications”, Journal of Statistical Physics, vol.64, Nos.2/3, pp.683-739. [2] DALLMANN, U., HERBERG, TH., GEBING, H., SU, W.-H. & ZHANG, H.-Q., 1995, “Flow field diagnostics: topological flow changes and spatio-temporal flow structure”. AIAA Paper 95-0791. [3] CHERRY, N. J., HILLIER, R. & LATOUR, M.P., 1984, “Unsteady measurements in a separated and reattaching flow”. J. Fluid Mech., vol.144, pp.13–46. [4] HAGGMARK, C.P., BAKCHINOV, A.A. & ALFREDSSON, P.H., 2000, “Experiments on a two-dimensional laminar separation bubble”. Phil. Trans. R. Soc. Lond. A, vol. 358, 3193–3205. [5] HASSELMANN, K., 1976, “Stochastic climate models”, Tellus, vol.28, pp.473-485. [6] LUMLEY, J.L., 1967, “The structure of inhomogeneous turbulent flows”, Atm.Turb. and Radio Wave Prop., Yaglom and Tatarsky eds., Nauka, Moskva, 1967, pp.166-178. [7] TROPEA, C., YARIN, A.L., FOSS, J.F. eds., 2007, “Handbook of Experimental Fluid Mechanics”, Springer. [8] URUBA, V., KNOB, M., 2007, “Dynamics of Controlled Boundary Layer Separation”, Colloquium Fluid Dynamics 2007, Praha, ÚT AV ČR, v.v.i., (Eds.: Jonáš P.; Uruba V.), pp.95-96. [9] URUBA, V., 2009, “Methods of Dynamical Analysis of Spatio-Temporal Data”, XXIII Symposium on Anemometry, Holany-Litice, Institute of Hydrodynamics AS CR, v.v.i., (Eds. Chára, Z.; Klaboch, L.), pp.110-115. [10] URUBA, V., 2009, “Dynamics of a Boundary Layer Separation”, Conference on Modelling Fluid Flow (CMFF’09), The 14th International Conference on Fluid Flow Technologies, Budapest, Hungary, September 9-12, 2009, pp.269-275. [11] URUBA, V., 2010, "A New Method for Vector Fields Dynamics Analysis", Engineering Mechanics 2010, (Ed. Zolotarev, I.), Svratka, May 10 - 13, pp.  163 164  

Dekompozice dat získaných metodou PIV s aplikací na výzkum dynamiky odtržení mezní vrstvy

http://dspace.vsb.cz/bitstream/10084/89181/2/1830.pdf

Dekompozice dat získaných metodou PIV s aplikací na výzkum dynamiky odtržení mezní vrstvy

Abstract

Similar works

Full text

Available Versions

DSpace at VSB Technical University of Ostrava