Next Article in Journal
Turbulence Measurements Downstream of a Combustor Simulator Designed for Studies on the Combustor–Turbine Interaction
Previous Article in Journal
Comparison Study of the kkLω and γReθ Transition Model in the Open-Water Performance Prediction of a Rim-Driven Thruster
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Experimental Database for the Analysis of Bursting of a Laminar Separation Bubble †

Department of Mechanical, Energy, Management and Transport Engineering (DIME), University of Genova, 16145 Genova, Italy
*
Author to whom correspondence should be addressed.
This manuscript is an extended version of the ETC2023-176 meeting paper published in the Proceedings of the 15th European Turbomachinery Conference, Budapest, Hungary, 24–28 April 2023.
Int. J. Turbomach. Propuls. Power 2024, 9(1), 3; https://doi.org/10.3390/ijtpp9010003
Submission received: 5 July 2023 / Revised: 10 July 2023 / Accepted: 29 November 2023 / Published: 10 January 2024

Abstract

:
The bursting phenomenon consists in the switch of a laminar separation bubble from a short to a long configuration. In the former case, reduced effects on profile pressure distribution are typically observed with respect to the attached condition. On the contrary, long bubbles provoke significant variations in the loading coefficient upstream of the separation position, with increased risk of stall of the lifting surfaces. The present work presents an experimental database describing separated boundary layers evolving under different Reynolds numbers, adverse pressure gradients and free-stream turbulence levels. Overall, more than 80 flow conditions were tested concerning short and long bubbles for the characterization of separated flows under turbine-like conditions. Measurements were performed on a flat plate geometry using a fast-response Particle Image Velocimetry (PIV) system. For each flow case, two sets of 6000 flow records were acquired with an acquisition frequency equal to 300 and 1000 Hz. Based on existing criteria for the identification of the bursting phenomenon, the flow cases were clustered in terms of short and long bubble states. Additionally, the kind of instability (i.e., convective or absolute) developing into the separated boundary layer was identified based on flow statistics. The present data captures the existing link between the bursting of a laminar separation bubble and the onset of the absolute instability of the separated shear layer, with stationary vortices forming in the dead air region.

Graphical Abstract

1. Introduction

Under conditions of strong adverse pressure gradients (APG), low turbulence intensity ( T u ), and low Reynolds numbers, the boundary layer (BL) may experience separation, leading to increased losses. If turbulent reattachment of the boundary layer occurs, a laminar separation bubble (LSB) forms on the lifting surface. Various applications, such as unmanned aerial vehicles, low-pressure turbines, low-pressure compressors, helicopter blades, and wind turbines, operate under this specific flow configuration. Understanding the flow conditions that cause a relatively short bubble to transition into a long one is crucial, as it significantly alters the pressure distribution profile and heightens the risk of stall. This abrupt change in the bubble state is commonly referred to in the literature as bubble bursting (Marxen and Henningson [1]). Due to its detrimental nature, significant efforts have been dedicated over the past decades to comprehending bubble bursting. This includes the development of empirical criteria, particularly focused on its prediction (see, for instance, Von Doenhoff [2], Owen and Klanfer [3], Crabtree [4]). Gaster [5] conducted a comprehensive analysis of the bursting process, proposing a two-parameter criterion for predicting changes in the bubble state. The author conducted experiments on a flat plate with a variable Reynolds number and pressure gradient, revealing a correlation between the momentum thickness Reynolds number at separation and non-dimensional parameter P G = θ s 2 / ν ( Δ U / Δ X ) . Here, P G considers the momentum thickness at separation ( θ s ) and the velocity gradient across the bubble length ( Δ U / Δ X ) with respect to an attached flow condition. In more recent works, Diwan et al. [6] and Mitra and Ramesh [7] introduced a single-parameter criterion based on the newly defined quantity P D C R = h 2 / ν ( Δ U / Δ X ) a c t . Diwan et al. [6] incorporated the maximum bubble height (h) as an influencing parameter for identifying the bubble state, and ( Δ U / Δ X ) a c t represents the actual velocity gradient between the separation and reattachment locations. Bursting is expected to occur when P D C R < 28 . This criterion has been adopted in subsequent literature works focusing on external aerodynamics at very low turbulence levels (see, for example, Serna and Lázaro [8], Alferez et al. [9]). Alferez et al. [9] suggested that the bursting process may be triggered by modifications in shear layer characteristics due to increased distance from the wall rather than a shift in the instability nature from convective to absolute, contrary to the original asumption of Gaster [5]. Among other studies, Almutairi et al. [10] conducted numerical simulations of laminar separation bubbles (LSB) undergoing bursting at very low Reynolds numbers, observing more irregular bubble behavior compared to that in previous experimental works at higher Reynolds numbers. In the recent work of Eljack et al. [11], low-frequency oscillations of an LSB forming over an airfoil were identified as precursors to complete profile stall. Then, the above-mentioned works clearly highlight the complexity of the bursting phenomenon, which is still far to be completely understood. Indeed, the effects occurring due to flow parameter variation are mostly obtained by comparison to different works, concerning both experiments and numerical simulations.
This paper presented a recently obtained experimental database concerning laminar separation bubbles that form on a flat plate under various T u levels and pressure gradients. The experiments encompassed different Reynolds numbers for each combination of these two parameters, enabling the observation of both short and long bubble states. In total, approximately 80 distinct flow conditions were examined using Particle Image Velocimetry. The detailed measurements conducted in this study make the dataset well-suited for characterizing the statistical progression of separation-induced transitions in both short and long bubble states. Additionally, the dataset facilitates an examination of the predominant flow dynamics and their associated stability characteristics. Interested parties may request access to the data presented in this paper directly from the authors. This work is an extended version of the conference paper by Dellacasagrande et al. [12].

2. Setup of Experiments

The present experiments were performed in the wind tunnel of the Aerodynamic and Turbomachinery laboratory at the University of Genova. The wind tunnel operated at low speeds, with a Mach number well below 0.1. Although the experiments were conducted at subsonic regimes, the collected data are representative of laminar separation bubbles evolving at high speeds, provided that shock-induced separation does not occur, and the same effective adverse pressure gradients are considered. Indeed, shock-induced boundary layer separation differs significantly from that induced by a continuous APG distribution. In the conducted experiments, a flat plate with a 4:1 elliptic leading edge was positioned between two adjustable endwalls, allowing for variation in the streamwise pressure gradient imposed on the boundary layer (see Figure 1). The flow was accelerated from the plate leading edge to the channel throat using a converging fixed geometry. The pressure gradient could be adjusted in the rear part of the plate. The opening angle of the endwalls ( α ) was set to 7, 9, and 12 degrees, corresponding to acceleration parameter ( A P ) values of −0.18, −0.27, and −0.41, respectively, with A P = L U , 0 Δ U Δ x . This parameter is defined based on plate length (L), the free-stream velocity at the channel throat ( U , 0 ), and the average velocity gradient along the rear part of the plate. The plate dimensions were 300 mm in length and width, providing a 2D time-mean flow at the mid-section of the channel, where Particle Image Velocimetry (PIV) measurements were conducted.
Free-stream turbulence intensity was controlled using turbulence-generating grids located 500 mm upstream of the plate leading edge. The root mean square of velocity fluctuations was measured by Laser Doppler Velocimetry (LDV) at the inlet section of the test section. Velocity data were acquired in the free-stream region at a mean sampling rate of 10 kHz over a sampling period of 120 s. Table 1 provides geometric parameters for the grids used, corresponding T u levels, and a streamwise integral length-scale ( L x ), computed using the auto-correlation function of LDV data under Taylor’s hypothesis of frozen turbulence. For each combination of free-stream turbulence and pressure gradient, approximately 9 Reynolds numbers ( R e L ) were tested. This parameter is defined based on the plate length and the free-stream velocity at the plate leading edge. The acquired dataset includes various flow conditions, as listed in Table 2, presenting the tested values of influencing parameters.
PIV measurements were performed in this work in a wall-normal plane aligned with the meridional section of the plate. To reduce laser reflections, the plate was painted black and the laser incidence angle was minimized by positioning the laser optics far behind the plate trailing edge instead of illuminating the plate from the top. Two cameras were used simultaneously to acquire the evolution of the boundary layer from the channel throat to the end of the plate. This guaranteed that the entire evolution of the laminar separation bubbles forming over the plate surface was recorded for the overall combinations of the above-mentioned parameters. Moreover, using two cameras allowed us extension of the field of view while keeping the spatial resolution sufficiently high to solve, with great detail, the dominant coherent structures within the flow. With the aim of providing a qualitative view of the flow dynamics captured in the present data base, Figure 2 reports exemplary flow visualizations of the operating PIV system for a very low Reynolds number case. The dark region in the pictures corresponds to the separated shear layer evolving above the plate surface. Depending on the flow Reynolds number, wall-normal oscillations of the separated shear layer were observed, characterized by significantly different central frequency (Boutilier and Yarusevych [13]). As it is further discussed in this work, vortical structures were seen to form within the recirculating flow region at a low Reynolds number. Such structures interact with the overlaying shear layer, likely provoking the shedding of larger vortices near the position of maximum shear, similar to what is shown in Figure 2. On the contrary, the classic inviscid convective Kelvin–Helmholtz instability was found to dominate at a high Reynolds number, for which short laminar separation bubbles occur. This provided evidence of the usability of the present experimental database for the characterization of separation-induced transition with reference to different states of a laminar separation bubble.
The PIV system employed in this study was a dual-cavity ND:YLF pulsed laser (Litron LDY 300) with a pulse energy of 30 mJ at a repetition rate of 1000 Hz and a wavelength of 527 nm. Two high-sensitive SpeedSense M340 digital cameras equipped with a cooled 2560 × 1600 pixels CMOS matrix were used to capture the overall velocity field. The magnification factor for these experiments was set to approximately 0.16. For the PIV image analysis, a multi-grid algorithm was implemented to compute adaptive cross-correlation. Initially, a 32 × 32 pixel interrogation area was used, followed by a successive refinement to 16 × 16 pixels with a 50% overlap. This resulted in a vector grid spacing of 0.41 mm. A peak validation method was employed to distinguish between valid and invalid vectors. Additionally, a Gaussian fitting procedure was applied to achieve sub-pixel recognition accuracy of particle displacement equal to 0.1 pixel. The uncertainty in instantaneous velocity measurements was estimated to be less than 3% in the free-stream region and increased to 6% in the boundary layer. This estimation was based on the peak ratio method proposed by Charonko and Vlachos [14], which correlates particle displacement uncertainty with the ratio between the dominant correlation peak and the second highest one. Commercial software DynamicStudio 7.3, used for these measurements, implemented this method. Instantaneous maps of uncertainty for particle displacement were generated, from which velocity uncertainty values were computed (see also Canepa et al. [15,16]). Complementary analysis regarding the impact of near-wall velocity gradients was conducted following the approach outlined in the work of Wilson and Smith [17]. This comprehensive assessment ensures a thorough understanding of the uncertainties associated with the measured velocity field.
In order to achieve a detailed characterization of the flow dynamics and also of the statistics of short and long laminar separation bubbles, two sets of 6000 snapshots were collected at 1000 and 300 Hz for each combination of the influencing parameters. The statistical reliability of the so-acquired database was assessed based on the temporal auto-correlation functions of the fluctuating velocity field for different cases. Figure 3 shows exemplary auto-correlation functions for low (left) and high (right) Reynolds numbers at the maximum pressure gradient and the lowest T u level, i.e., α = 12 and T u = 1.5%. This represents the most burdensome combination of APG and the T u level for the occurrence of bubble bursting. The colored curves reported in Figure 3 refer to different numbers of PIV snapshots adopted for the computation of the auto-correlation function. Data reported in Figure 3 clearly indicate that the characteristic period of the main coherent fluctuations is well captured using only 3000 snapshots (see the position of the first negative peak in the auto-correlation function). It should be noted that the time shift ( Δ t ) is normalized by the overall acquisition period ( T a c q ). Accordingly, about 100 to 500 shedding cycles were acquired from the lowest to the highest Reynolds numbers at a sampling rate of 1 kHz. For the highest Reynolds numbers, the bubble time-scale was discretized with about 20 sample points, which increased to about 70 for the lowest R e L values. Further details on the adopted test section and measuring techniques can be found in the previous works by Verdoya et al. [18] and Canepa et al. [19].

3. Time-Mean Flow Fields and Boundary Layer Statistical Parameters

In order to highlight the response of the laminar separation bubble to flow parameter variation, this section reports some statistical quantities of the separated boundary layer for selected flow cases among those tested.
Figure 4 shows the contour plots of the normalized streamwise velocity component ( u / U e ) at fixed T u = 1.5% and α = 12 for a variable Reynolds number. This latter is the parameter showing the greatest influence on the definition of the bubble state (i.e., short and long one). Similar time-mean velocity distributions were indeed observed for the other combinations of the pressure gradient and the T u level, not shown here for brevity. Results reported in Figure 4 clearly highlight the increase in the bubble size with the reduction in the Reynolds number. Particularly, moving from the bottom to the top plot, the acquired flow field switches from a very short bubble to a fully stalled case in which the turbulent reattachment does not take place before the end of the plate. The boundary layer separation occurs for all cases in the range of 0.39 < x / L < 0.44 ; thus, the growth of the bubble length occurs mostly due to the downstream shift of the bubble maximum displacement and reattachment positions. Interestingly, for R e L < 40,600, the growth rate of the bubble due to Reynolds number reduction significantly increases (the Reynolds number was almost uniformly sampled). This sudden increment of the bubble size can be directly ascribed to the occurrence of the bursting process, as further discussed in the following sections; thus, a critical Reynolds number for bursting was identified for the different flow cases.
Figure 5 shows the contour plots of the root mean square (rms) of the streamwise velocity fluctuations ( u r m s / U e ) for the same cases reported in Figure 4. The region where the maximum rms values occur is seen to move downstream as the Reynolds number is reduced, according to the increment of the bubble length. Additionally, the saturation level of the rms of fluctuations becomes increasingly higher as R e L is reduced. This is due to the occurrence of significantly larger vortices which are responsible for higher levels of fluctuations. As is shown in the following section, the range of the Reynolds number variation tested here was found to cause a substantial modification of the leading structures observed near the bubble maximum displacement, where the maximum rms values occur. Moreover, the growth of fluctuations in the fore part of the bubble was shown to reduce at low R e L values, which is the cause for delayed transition and reattachment.
To better illustrate the impact of variations in flow parameters on the time-mean shape of the laminar separation bubble and its stability characteristics, Figure 6 displays streamwise distributions of non-dimensional boundary layer (BL) displacement thickness ( δ * / L , top plots), recirculating height ( h r / L , center plots), and local maxima of root mean square of normal-to-the-wall fluctuations ( v r m s / U 0 , bottom plots) for cases with α = 12 T u = 1.5 % , α = 12 T u = 2.5 % , and α = 9 T u = 2.5 % (from left to right). The influence of Reynolds number variation is shown within the same range of R e L values in each plot. The displacement thickness distributions (top plots) exhibit a consistent trend in the initial part of the measurement domain when varying the Reynolds number at a constant adverse pressure gradient and turbulence intensity. Higher values of δ * / L near the inlet at lower R e L result from boundary layer thickening, as expected. Comparing cases with different turbulence intensities ( α = 12 T u = 1.5 % and α = 12 T u = 2.5 % ), an increase in T u leads to a reduced growth of δ * / L upstream of the separation position before notable divergence occurs. A reduction in streamwise pressure gradient (comparing α = 12 T u = 2.5 % and α = 9 T u = 2.5 % ) produces a similar effect. Elevated free-stream turbulence enhances the boundary layer resistance to adverse pressure gradients. Conversely, a decrease in channel opening angle reduces velocity gradients, retarding boundary layer growth and delaying detachment according to previous separation criteria (e.g., Curle and Skan [20]). Moving towards the plate end, significant divergence in the curves for different combinations of T u and APG reveals the effects of Reynolds number variation on maximum bubble height and its location. The peak of δ * / L shifts upstream and diminishes notably as the Reynolds number increases, consistent with findings in prior literature works (e.g., Istvan et al. [21], Dellacasagrande et al. [22]).
Distributions of h r / L (mid-plots of Figure 6) further demonstrate the effects, occurring due to flow parameter variation, on the separation and reattachment positions. At the lowest T u level (left plot), the separation position remains almost unchanged at different R e L values, according to Figure 4. A delay in separation position is instead observed when increasing the T u level, as well as when reducing the APG, and the role of R e L variation becomes more relevant in this sense. Nevertheless, the largest changes still concern the reattachment position and the bubble height. Notably, a marked growth of the bubble dimension occurs below a certain R e L value, which is explained by bursting occurrence. A similar behavior was indeed observed for the other combinations of the flow parameters tested here, which further indicates the presence of different bubble states within the current database.
Finally, the distributions of v r m s / U 0 (bottom plots of Figure 6) highlight the changes in the dynamic properties of the tested laminar separation bubbles. For all APG and T u levels, an earlier growth of v r m s / U 0 occurs at higher Reynolds numbers, which highlights the anticipation of the laminar to turbulent transition. This is also confirmed by the upstream shift of the streamwise position at which v r m s / U 0 saturates for the different cases, which occurs near the bubble maximum displacement position. Moreover, the higher the R e L , the higher the spatial growth of v r m s / U 0 , which denotes a modification of stability characteristics of the separated shear layer (see, e.g., Simoni et al. [23]). Again, there exists an R e L value below which a marked modification, i.e., reduction, in v r m s / U 0 growth is observed. Data reported in Figure 6 make clear that this behavior is linked to the increment of the bubble size (see h r / L distributions). The growth rate of fluctuations is found to also reduce at lower APG, while the T u level has minor effects in this sense.

4. Instantaneous Perturbation Velocity Fields

The present section reports exemplary PIV snapshots to highlight the main structures driving the transition process of an LSB as well as changes in the shear layer stability properties due to bursting occurrence. Since the bubble state has been found to be mostly affected by the flow Reynolds number variation, Figure 7, Figure 8 and Figure 9 compare high and low Reynolds number fluctuating velocity fields at different T u levels and pressure gradients.
In the top plot of Figure 7 (case T u = 1.5 % α = 12 with high R e L ), counter-rotating vortices are evident near the maximum displacement of the bubble (indicated by blue circles). These vortices form in the region of maximum shear and subsequently grow as they move in the main flow direction, indicating convective amplification. Further downstream, smaller structures for x / L > 0.75 suggest that vortex breakup and transition occur. Such spanwise vortices have been shown to constitute a significant source of losses due to the associated high levels of turbulence production (see e.g., Dellacasagrande et al. [24]). What is shown in the top plot in Figure 7 illustrates the typical pattern associated with the Kelvin–Helmholtz (K-H) instability process, driving the convective growth of disturbances in short laminar separation bubbles and the subsequent turbulent reattachment of the boundary layer. In contrast, the bottom plot in Figure 7 depicts a different scenario for the low R e L case. As the Reynolds number decreases below a certain threshold, the train of counter-rotating vortices characteristic of the K-H instability no longer occurs. Instead, significantly larger vortices start to form within the recirculating flow region (indicated by blue circles in the first snapshot). These spanwise vortices are observed to grow over time at nearly fixed positions (indicated by blue arrows in the plots) rather than being convectively amplified. This indicates a substantial modification of the stability characteristics of the laminar separation bubble, as also evident in v r m s / U 0 distributions shown in Figure 6. Specifically, for the low R e L cases, the occurrence of absolute instability is expected, as further discussed in the following section. This results in a modification of the growth rate of fluctuations and the delayed transition and reattachment of the laminar separation bubble.
Similar effects due to changes in the Reynolds number are observed in Figure 8 and Figure 9 at different T u levels and pressure gradients: K-H instability dominates at high R e L , regardless of the free-stream turbulence and the average velocity gradient, whereas the dynamics of an LSB are seen to change significantly at low R e L values. It is worth noting that all the low Reynolds number cases reported herein suggest that interaction may occur between vortices forming within the recirculating region, which are likely due to the absolute instability of disturbances, and those originating in the region of maximum shear. More precisely, as spanwise vortex-like structures grow over time within the recirculating region, they appear to provoke wall-normal oscillation of the overlaying shear layer, which then exhibits marked sinuous motions leading to the generation of larger-scale vortical structures near the bubble maximum displacement (see the bottom plots of Figure 8 and Figure 9). Namely, coherent fluctuations moving at a low speed within the recirculating region act to shrink adjacent vorticity nuclei (blue circles) and the separated shear layer, provoking the successive shedding of convective structures moving at much higher speed (see the inclination of the blue arrows depicted in the plots). The present results indicate that in the long bubble regime, convective and absolute instability may participate in provoking the shear layer instability and its turbulent transition. The present experimental database is therefore usable for the characterization of the flow physics driving the transition of laminar separation bubbles forming over airfoils and turbomachinery bladings.
Based on previous studies (see e.g., Simoni et al. [25]), the change in stability properties of an LSB may be quantified considering the wave number (k) of the largest coherent structures shed by the bubble once scaled via the shear layer thickness at the separation position (l). To overcome the limitation of the visual inspection of PIV data, the wavelength and frequency of the main coherent structures shed by the bubble were computed for each case by means of Proper Orthogonal Decomposition (POD): POD modes provided the dominant wavelength ( λ ) for each case, while the associated central frequency was obtained by means of Fourier transform of the corresponding temporal coefficients. A non-dimensional wave-number k l was therefore defined as k l = 2 π l / λ . This non-dimensional parameter was equal to about 0.81–0.92 in the short bubble regime, regardless of R e L , the T u level and APG. This indicates a self-similarity of the main vortical structures shed by the bubble. Instead, k l significantly deviated from this value when the Reynolds number reduced below a certain threshold. For the cases reported in Figure 7, Figure 8 and Figure 9, the scaling parameter k l reduced to about 0.55–0.68 when passing from high to low R e L conditions. Similarly, the group velocity of the dominant coherent vortical structures was computed and scaled with the local external velocity. The non-dimensional traveling speed of the vortical structures formed in the LSB reduced from about 0.5 in the short regime to null and even negative values for the structures originating within the recirculating region of the bubble. This is linked to the dominance of temporal growth of fluctuations and the onset of absolute instability of disturbances. Again, this indicates a substantial modification of the instability processes driving turbulence generation in the separated shear layer.

5. Data Clustering Based on the Bubble State and BL Instability Characteristics

In this section, statistical parameters defined in previous literature works are used to cluster the current data in terms of short and long bubble types, as well as with reference to the instability mechanisms driving the disturbance growth in the boundary layer. This allows discussion of the eventual link between the occurrence of the bubble bursting and the change in the stability characteristics of the separated shear layer that was highlighted via the inspection of instantaneous flow records. The bubble state is identified based on the P D C R parameter previously defined by Mitra and Ramesh [7]. Then, two additional parameters are used for the characterization of the instability process for the different cases examined. The first one is the ratio between the maximum height of the recirculating flow region and the local boundary layer displacement thickness, i.e., h r / δ * . According to the work of Maucher et al. [26], when this quantity is higher than 0.6, the absolute instability of fluctuations is expected to occur. Similarly, the switch from the convective to the absolute instability is expected to take place when the maximum reversed flow is larger than 15% of the local external velocity, that is, u m i n / U e < 0.15 (see Maucher et al. [26], Pauley et al. [27]). These two different criteria for the occurrence of the absolute instability of fluctuations are combined in Figure 10 with the P D C R -based criterion for the occurrence of bubble bursting. The blue plane in the plot indicates the bursting threshold corresponding to P D C R = 28 , as suggested in the work of Mitra and Ramesh [7]. Based on the adopted bursting criterion, most of the acquired flow cases are representative of a short bubble state (i.e., P D C R > 28 ), whereas part of them is classified as long ones. These latter conditions are characterized by h r / δ * and u m i n / U e values that are higher than the proposed thresholds for the occurrence of the absolute instability in the laminar separation bubble (see the dashed lines in the plot). Otherwise, all the short bubbles fall within the range typical of convective one, i.e., that of the K-H instability process. What is shown in Figure 10 clearly indicates the link between the change in the bubble state and the modification of its stability characteristics, which can be further studied in future works by means of the present experimental database.

6. Conclusions

The current study introduced an extensive experimental database focusing on laminar separation bubbles forming on a flat plate, exploring variable Reynolds numbers, turbulence intensity levels, and pressure gradients. Through the analysis of time-mean flow distributions and the application of a bursting criterion from the existing literature, the study clearly identified the presence of both short and long bubble states across different combinations of flow parameters. This dataset is well suited for characterizing such phenomena, potentially offering engineering correlations for predicting bursting events. Furthermore, the stability characteristics of the acquired laminar separation bubbles were found to be significantly influenced by changes in external flow conditions. The scale of vortical structures driving transition and the mechanisms promoting disturbance growth in the boundary layer exhibited notable variations with changes in the flow Reynolds number. Notably, the bursting process seemed closely linked to the onset of absolute instability in the separated shear layer. This experimental database is intended to be a valuable resource for studying separation-induced transition across a wide range of flow conditions and gaining insights into the underlying mechanisms governing the transition between short and long bubble states. To this end, the present data base is intended to be open and available to the scientific community.

Author Contributions

Conceptualization, M.D., D.S. and D.L.; methodology, M.D., D.S. and D.L.; software, M.D.; validation, M.D. and D.L.; writing—original draft preparation, M.D. and M.U.; writing—review and editing, M.D., D.L. and M.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available within the article. Additional data are available upon request to the authors.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
R e L Reynolds number based on plate length
T u Free-stream turbulence intensity
APGAdverse pressure gradient
APAcceleration parameter
LPlate length
LSBLaminar separation bubble
R u u Streamwise auto-correlation function
h r Recirculating region height
δ * Boundary layer displacement thickness
θ Boundary layer momentum thickness
α Opening angle of the top-wall of the channel
ν Kinematic viscosity

References

  1. Marxen, O.; Henningson, D.S. The Effect of Small-Amplitude Convective Disturbances on the Size and Bursting of a Laminar Separation Bubble. J. Fluid Mech. 2011, 671, 1–33. [Google Scholar] [CrossRef]
  2. Von Doenhoff, A.E. A Preliminary Investigation of Boundary-Layer Transition along a Flat Plate with Adverse Pressure Gradient; NTRS-NASA Technical Report; 300 Hidden Figures Way SW: Washington, DC, USA, 1938. [Google Scholar]
  3. Owen, P.; Klanfer, L. On the Laminar Boundary Layer Separation from the Leading Edge of a Thin Aerofoilv; Technical Report; Aeronautical Research Council: London, UK, 1953. [Google Scholar]
  4. Crabtree, L. The Formation of Regions of Separated Flow on Wing Surfaces; Citeseer: London, UK, 1959. [Google Scholar]
  5. Gaster, M. The Structure and Behaviour of Laminar Separation Bubbles; Technical Report; Aeronautical Research Council: London, UK, 1967. [Google Scholar]
  6. Diwan, S.S.; Chetan, S.; Ramesh, O. On the bursting criterion for laminar separation bubbles. In Proceedings of the IUTAM Symposium on Laminar-Turbulent Transition, Bangalore, India, 2004; Springer: Dordrecht, The Netherlands, 2006; pp. 401–407. [Google Scholar]
  7. Mitra, A.; Ramesh, O. New correlation for the prediction of bursting of a laminar separation bubble. AIAA J. 2019, 57, 1400–1408. [Google Scholar] [CrossRef]
  8. Serna, J.; Lázaro, B. On the bursting condition for transitional separation bubbles. Aerosp. Sci. Technol. 2015, 44, 43–50. [Google Scholar] [CrossRef]
  9. Alferez, N.; Mary, I.; Lamballais, E. Study of stall development around an airfoil by means of high fidelity large eddy simulation. Flow Turbul. Combust. 2013, 91, 623–641. [Google Scholar] [CrossRef]
  10. Almutairi, J.H.; Jones, L.E.; Sandham, N.D. Intermittent bursting of a laminar separation bubble on an airfoil. AIAA J. 2010, 48, 414–426. [Google Scholar] [CrossRef]
  11. Eljack, E.; Soria, J.; Elawad, Y.; Ohtake, T. Simulation and characterization of the laminar separation bubble over a NACA-0012 airfoil as a function of angle of attack. Phys. Rev. Fluids 2021, 6, 034701. [Google Scholar] [CrossRef]
  12. Dellacasagrande, M.; Lengani, D.; Simoni, D.; Ubaldi, M. An Experimental Database for the Analysis of Bursting of a Laminar Separation Bubble. In Proceedings of the 15th European Turbomachinery Confeerence, Budapest, Hungary, 24–28 April 2023; Paper n. ETC2023-176. Available online: https://www.euroturbo.eu/publications/conference-proceedings-repository/ (accessed on 5 July 2023).
  13. Boutilier, M.S.; Yarusevych, S. Parametric study of separation and transition characteristics over an airfoil at low Reynolds numbers. Exp. Fluids 2012, 52, 1491–1506. [Google Scholar] [CrossRef]
  14. Charonko, J.J.; Vlachos, P.P. Estimation of uncertainty bounds for individual particle image velocimetry measurements from cross-correlation peak ratio. Meas. Sci. Technol. 2013, 24, 065301. [Google Scholar] [CrossRef]
  15. Canepa, E.; Cattanei, A.; Jafelice, F.; Zecchin, F.M.; Parodi, D. Effect of rotor deformation and blade loading on the leakage noise in low-speed axial fans. J. Sound Vib. 2018, 433, 99–123. [Google Scholar] [CrossRef]
  16. Canepa, E.; Cattanei, A.; Mazzocut Zecchin, F. Leakage noise and related flow pattern in a low-speed axial fan with rotating shroud. Int. J. Turbomach. Propuls. Power 2019, 4, 17. [Google Scholar] [CrossRef]
  17. Wilson, B.M.; Smith, B.L. Uncertainty on PIV mean and fluctuating velocity due to bias and random errors. Meas. Sci. Technol. 2013, 24, 035302. [Google Scholar] [CrossRef]
  18. Verdoya, J.; Dellacasagrande, M.; Lengani, D.; Simoni, D.; Ubaldi, M. Inspection of structures interaction in laminar separation bubbles with extended proper orthogonal decomposition applied to multi-plane particle image velocimetry data. Phys. Fluids 2021, 33, 043607. [Google Scholar] [CrossRef]
  19. Canepa, E.; Cattanei, A.; Moradi, M.; Nilberto, A. Experimental Study of the Leakage Flow in an Axial-Flow Fan at Variable Loading. Int. J. Turbomach. Propuls. Power 2021, 6, 40. [Google Scholar] [CrossRef]
  20. Curle, N.; Skan, S. Approximate methods for predicting separation properties of laminar boundary layers. Aeronaut. Q. 1957, 8, 257–268. [Google Scholar] [CrossRef]
  21. Istvan, M.S.; Kurelek, J.W.; Yarusevych, S. Turbulence Intensity Effects on Laminar Separation Bubbles Formed over an Airfoil. AIAA J. 2018, 56, 1335–1347. [Google Scholar] [CrossRef]
  22. Dellacasagrande, M.; Barsi, D.; Lengani, D.; Simoni, D.; Verdoya, J. Response of a flat plate laminar separation bubble to Reynolds number, free-stream turbulence and adverse pressure gradient variation. Exp. Fluids 2020, 61, 128. [Google Scholar] [CrossRef]
  23. Simoni, D.; Lengani, D.; Ubaldi, M.; Zunino, P.; Dellacasagrande, M. Inspection of the dynamic properties of laminar separation bubbles: Free-stream turbulence intensity effects for different Reynolds numbers. Exp. Fluids 2017, 58, 66. [Google Scholar] [CrossRef]
  24. Dellacasagrande, M.; Lengani, D.; Simoni, D.; Ubaldi, M.; Zunino, P. Experimental Investigation on the Loss Production Mechanisms in Transitional Boundary Layers. In Proceedings of the Turbo Expo: Power for Land, Sea, and Air. American Society of Mechanical Engineers, Virtual, 21–25 September 2020; Volume 84102, p. V02ET41A019. [Google Scholar]
  25. Simoni, D.; Ubaldi, M.; Zunino, P. Experimental investigation of flow instabilities in a laminar separation bubble. J. Therm. Sci. 2014, 23, 203–214. [Google Scholar] [CrossRef]
  26. Maucher, U.; Rist, U.; Wagner, S. Refined interaction method for direct numerical simulation of transition in separation bubbles. AIAA J. 2000, 38, 1385–1393. [Google Scholar] [CrossRef]
  27. Pauley, L.L.; Moin, P.; Reynolds, W.C. The structure of two-dimensional separation. J. Fluid Mech. 1990, 220, 397–411. [Google Scholar] [CrossRef]
Figure 1. Test section and PIV instrumentation layout. PIV measuring planes are highlighted with green boxes.
Figure 1. Test section and PIV instrumentation layout. PIV measuring planes are highlighted with green boxes.
Ijtpp 09 00003 g001
Figure 2. Flow visualization of shear layer instability at low Reynolds number. The shear layer region corresponds to the shaded portion of the frames.
Figure 2. Flow visualization of shear layer instability at low Reynolds number. The shear layer region corresponds to the shaded portion of the frames.
Ijtpp 09 00003 g002
Figure 3. Temporal auto-correlation function of velocity fluctuations for cases (left) R e L = 21,017 − T u = 1.5 % α = 12 and (right) R e L = 66,173 − T u = 1.5 % α = 12 . The auto-correlation functions were obtained from the dataset sampled at 1 kHz over an acquisition period of 6 s. The time shift is scaled with the overall acquisition period. The number of snapshots adopted is indicated in the figure legend.
Figure 3. Temporal auto-correlation function of velocity fluctuations for cases (left) R e L = 21,017 − T u = 1.5 % α = 12 and (right) R e L = 66,173 − T u = 1.5 % α = 12 . The auto-correlation functions were obtained from the dataset sampled at 1 kHz over an acquisition period of 6 s. The time shift is scaled with the overall acquisition period. The number of snapshots adopted is indicated in the figure legend.
Ijtpp 09 00003 g003
Figure 4. Normalized time-mean velocity ( u / U e ) for a variable Reynolds number at a fixed T u = 1.5 % and α = 12 .
Figure 4. Normalized time-mean velocity ( u / U e ) for a variable Reynolds number at a fixed T u = 1.5 % and α = 12 .
Ijtpp 09 00003 g004
Figure 5. Normalized root mean square of velocity fluctuations ( u r m s / U e ) for a variable Reynolds number at a fixed T u = 1.5 % and α = 12 .
Figure 5. Normalized root mean square of velocity fluctuations ( u r m s / U e ) for a variable Reynolds number at a fixed T u = 1.5 % and α = 12 .
Ijtpp 09 00003 g005
Figure 6. Streamwise variation in the displacement thickness scaled with the plate length (top, δ * / L ), of the recirculating height scaled with the plate length (center, h r / L ) and of the rms of normal to the wall fluctuations scaled with the inlet free-stream velocity (bottom, v r m s / U 0 ). From left to right: α = 12 T u = 1.5 % , α = 12 T u = 2.5 % and α = 9 T u = 2.5 % . The Reynolds variation effect is highlighted in each plot.
Figure 6. Streamwise variation in the displacement thickness scaled with the plate length (top, δ * / L ), of the recirculating height scaled with the plate length (center, h r / L ) and of the rms of normal to the wall fluctuations scaled with the inlet free-stream velocity (bottom, v r m s / U 0 ). From left to right: α = 12 T u = 1.5 % , α = 12 T u = 2.5 % and α = 9 T u = 2.5 % . The Reynolds variation effect is highlighted in each plot.
Ijtpp 09 00003 g006
Figure 7. Instantaneous vector maps of fluctuating velocity field for cases (top) R e L = 66,200 − T u = 1.5 % α = 12 and (bottom) R e L = 21,000 − T u = 1.5 % α = 12 . Iso-contour lines of u / U e = 0 , 0.3 and 0.9 are depicted with red color. Blue arrows are adopted to track vortices at the different time instants. One in every two measured points is depicted.
Figure 7. Instantaneous vector maps of fluctuating velocity field for cases (top) R e L = 66,200 − T u = 1.5 % α = 12 and (bottom) R e L = 21,000 − T u = 1.5 % α = 12 . Iso-contour lines of u / U e = 0 , 0.3 and 0.9 are depicted with red color. Blue arrows are adopted to track vortices at the different time instants. One in every two measured points is depicted.
Ijtpp 09 00003 g007
Figure 8. Instantaneous vector maps of fluctuating velocity field for cases (top) R e L = 66,900 − T u = 1.5 % α = 9 and (bottom) R e L = 22,500 − T u = 1.5 % α = 9 . Iso-contour lines of u / U e = 0 , 0.3 and 0.9 are depicted with red color. Blue arrows are adopted to track vortices at the different time instants. One in every two measured points is depicted.
Figure 8. Instantaneous vector maps of fluctuating velocity field for cases (top) R e L = 66,900 − T u = 1.5 % α = 9 and (bottom) R e L = 22,500 − T u = 1.5 % α = 9 . Iso-contour lines of u / U e = 0 , 0.3 and 0.9 are depicted with red color. Blue arrows are adopted to track vortices at the different time instants. One in every two measured points is depicted.
Ijtpp 09 00003 g008
Figure 9. Instantaneous vector maps of fluctuating velocity field for cases (top) R e L = 57,400 − T u = 2.5 % α = 9 and (bottom) R e L = 20,700 − T u = 2.5 % α = 9 . Iso-contour lines of u / U e = 0 , 0.3 and 0.9 are depicted with red color. Blue arrows are adopted to track vortices at the different time instants. One in every two measured points is depicted.
Figure 9. Instantaneous vector maps of fluctuating velocity field for cases (top) R e L = 57,400 − T u = 2.5 % α = 9 and (bottom) R e L = 20,700 − T u = 2.5 % α = 9 . Iso-contour lines of u / U e = 0 , 0.3 and 0.9 are depicted with red color. Blue arrows are adopted to track vortices at the different time instants. One in every two measured points is depicted.
Ijtpp 09 00003 g009
Figure 10. Three-dimensional plot of flow cases in terms of the maximum reversed flow velocity ( u m i n / U e ), the non-dimensional height of the recirculating region ( h r / δ * ) and the bursting parameter P D C R . Purple and red lines indicate the h r / δ * and u m i n / U e limits, respectively, for absolute instability. Only part of the entire database is shown to increase plot readability.
Figure 10. Three-dimensional plot of flow cases in terms of the maximum reversed flow velocity ( u m i n / U e ), the non-dimensional height of the recirculating region ( h r / δ * ) and the bursting parameter P D C R . Purple and red lines indicate the h r / δ * and u m i n / U e limits, respectively, for absolute instability. Only part of the entire database is shown to increase plot readability.
Ijtpp 09 00003 g010
Table 1. Turbulence-generating grid characterization: free-stream turbulence intensity ( T u ), bar width (d), mesh size (M), porosity parameter P = ( 1 d / M ) 2 and streamwise integral length scale L x .
Table 1. Turbulence-generating grid characterization: free-stream turbulence intensity ( T u ), bar width (d), mesh size (M), porosity parameter P = ( 1 d / M ) 2 and streamwise integral length scale L x .
Case Tu (%)d (mm)M (mm)P L x (mm)
No-Grid1.5%---6.2
Low-Tu2.5%280.5611.5
High-Tu3.5%480.2513.2
Table 2. Flow cases investigated: the Reynolds number, acceleration parameter (AP) and the T u level are reported.
Table 2. Flow cases investigated: the Reynolds number, acceleration parameter (AP) and the T u level are reported.
AP−0.18−0.18−0.18−0.27−0.27−0.27−0.41−0.41−0.41
Tu1.5%2.5%3.5%1.5%2.5%3.5%1.5%2.5%3.5%
R e L 29,60028,80024,80021,00016,10018,20015,30015,90017,200
R e L 34,70033,20027,50022,50020,70023,80021,00020,60021,400
R e L 39,10038,10033,60029,10028,20025,30025,10023,80025,500
R e L 41,00043,10036,70037,50032,60029,20031,10028,00030,600
R e L 44,10045,40040,90046,80037,50035,10040,60032,00039,000
R e L 49,00050,90049,90057,80042,10040,40050,80035,20043,000
R e L 53,20056,80053,20060,10047,20044,90056,00040,60051,600
R e L 71,50063,60056,10066,90052,20051,40060,60046,70054,800
R e L 76,10068,80057,00069,70057,40052,60066,20056,60060,300
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dellacasagrande, M.; Lengani, D.; Simoni, D.; Ubaldi, M. An Experimental Database for the Analysis of Bursting of a Laminar Separation Bubble. Int. J. Turbomach. Propuls. Power 2024, 9, 3. https://doi.org/10.3390/ijtpp9010003

AMA Style

Dellacasagrande M, Lengani D, Simoni D, Ubaldi M. An Experimental Database for the Analysis of Bursting of a Laminar Separation Bubble. International Journal of Turbomachinery, Propulsion and Power. 2024; 9(1):3. https://doi.org/10.3390/ijtpp9010003

Chicago/Turabian Style

Dellacasagrande, Matteo, Davide Lengani, Daniele Simoni, and Marina Ubaldi. 2024. "An Experimental Database for the Analysis of Bursting of a Laminar Separation Bubble" International Journal of Turbomachinery, Propulsion and Power 9, no. 1: 3. https://doi.org/10.3390/ijtpp9010003

APA Style

Dellacasagrande, M., Lengani, D., Simoni, D., & Ubaldi, M. (2024). An Experimental Database for the Analysis of Bursting of a Laminar Separation Bubble. International Journal of Turbomachinery, Propulsion and Power, 9(1), 3. https://doi.org/10.3390/ijtpp9010003

Article Metrics

Back to TopTop