W.J. Baars
Please Note
24 records found
1
This thesis first establishes an experimental reference dataset for a canonical vortex generator test case and combines it with fully-resolved simulations, which serve as a numerical benchmark for assessing the performance of a jBAY low-fidelity model. The major aspects of vortex generator behaviour are characterized through the integrated aerodynamic forces, the force distributions over the vortex generator surface, and the resulting downstream vortex development to allow direct force and indirect vorticity measures to be used for model analysis.
The jBAY model is evaluated at varying incidence angles. Across all cases, the model produces a more diffused downstream vortex than the corresponding fully resolved simulations and experiments. However, at the lower incidence angle, the discrepancy is more pronounced, with the vortex shape and structure not adequately captured. This indicates an angle-dependent mesh requirement, with finer downstream resolution required at lower incidence than at higher incidence angles or in the corresponding fully resolved simulations. Further sensitivity analyses show that obtaining a normal-force magnitude and source distribution closer to the experimental measurements does not necessarily improve the downstream flow prediction. Overall, the results highlight the importance of downstream mesh resolution for the jBAY model, particularly at lower incidence, and demonstrate that an accurate force response alone does not guarantee an accurate prediction of the resulting vortex. ...
This thesis first establishes an experimental reference dataset for a canonical vortex generator test case and combines it with fully-resolved simulations, which serve as a numerical benchmark for assessing the performance of a jBAY low-fidelity model. The major aspects of vortex generator behaviour are characterized through the integrated aerodynamic forces, the force distributions over the vortex generator surface, and the resulting downstream vortex development to allow direct force and indirect vorticity measures to be used for model analysis.
The jBAY model is evaluated at varying incidence angles. Across all cases, the model produces a more diffused downstream vortex than the corresponding fully resolved simulations and experiments. However, at the lower incidence angle, the discrepancy is more pronounced, with the vortex shape and structure not adequately captured. This indicates an angle-dependent mesh requirement, with finer downstream resolution required at lower incidence than at higher incidence angles or in the corresponding fully resolved simulations. Further sensitivity analyses show that obtaining a normal-force magnitude and source distribution closer to the experimental measurements does not necessarily improve the downstream flow prediction. Overall, the results highlight the importance of downstream mesh resolution for the jBAY model, particularly at lower incidence, and demonstrate that an accurate force response alone does not guarantee an accurate prediction of the resulting vortex.
Large-Scale Aerodynamic Wake Survey System
DSE Fall 2026, Group 24
This study experimentally investigates whether inclining the HR orifice against the flow can induce a phase shift between the orifice pressure and velocity. If achieved, this could lead to attenuation at the resonance frequency as well. In addition, since ambient noise is significant in practical applications and its effect on resonator-turbulence interaction remains unexplored, the influence of noise on a single HR in grazing flow is also investigated. The impact of noise and orifice orientation on sweep, ejection, and turbulence production is also studied. Experiments were conducted in the Delft University Boundary Layer Facility at Reτ ≈ 2500. Time-resolved particle image velocimetry, hot-wire anemometry, and microphone measurements were employed. Three HR orientations (Vertical, Flow-Opposed, and Flow Aligned) were tested under no-noise conditions, tonal excitation at the resonance frequency of the HR, and white noise at two levels. The hot-wire results show that the Flow-Opposed HR attenuates more at sub-resonance and amplifies less at resonance than the other configurations. This suggests that a phase shift is induced, but it is not sufficient to achieve attenuation at resonance. The effect of tonal noise is strongly orientation-dependent: Vertical HR generally shows attenuation of Reynolds stresses, while Flow-Opposed HR shows amplification. Quadrant analysis reveals that sweep and ejection events follow the same trend across all configurations and frequencies. ...
This study experimentally investigates whether inclining the HR orifice against the flow can induce a phase shift between the orifice pressure and velocity. If achieved, this could lead to attenuation at the resonance frequency as well. In addition, since ambient noise is significant in practical applications and its effect on resonator-turbulence interaction remains unexplored, the influence of noise on a single HR in grazing flow is also investigated. The impact of noise and orifice orientation on sweep, ejection, and turbulence production is also studied. Experiments were conducted in the Delft University Boundary Layer Facility at Reτ ≈ 2500. Time-resolved particle image velocimetry, hot-wire anemometry, and microphone measurements were employed. Three HR orientations (Vertical, Flow-Opposed, and Flow Aligned) were tested under no-noise conditions, tonal excitation at the resonance frequency of the HR, and white noise at two levels. The hot-wire results show that the Flow-Opposed HR attenuates more at sub-resonance and amplifies less at resonance than the other configurations. This suggests that a phase shift is induced, but it is not sufficient to achieve attenuation at resonance. The effect of tonal noise is strongly orientation-dependent: Vertical HR generally shows attenuation of Reynolds stresses, while Flow-Opposed HR shows amplification. Quadrant analysis reveals that sweep and ejection events follow the same trend across all configurations and frequencies.
This thesis concerns the estimation of the cooling drag caused by these installations. The developed methodology considers the internal resistances from parts of the ducting and the heat exchanger as well as the external drag caused by the installation, offering nuanced insights into the field of thermal management that help engineers accurately estimate cooling drag at an early stage and make substantiated design choices. ...
This thesis concerns the estimation of the cooling drag caused by these installations. The developed methodology considers the internal resistances from parts of the ducting and the heat exchanger as well as the external drag caused by the installation, offering nuanced insights into the field of thermal management that help engineers accurately estimate cooling drag at an early stage and make substantiated design choices.
Transverse Forcing by Acoustic Excitation
For friction drag reduction in a turbulent boundary layer
Interactions of Multi-Rotors with Surfaces
Aerodynamic Characterization and Performance Modelling
The high power-to-weight ratio and efficiency of electric motors provide an excellent platform to explore disruptive propulsion system configurations. An observable design trend here is the tight integration of distributed, fixed-pitched rotors with aerodynamic surfaces, aimed at generating beneficial aerodynamic coupling effects to increase the aircraft’s efficiency and reduce acoustic emissions.
An auspicious integrated tilt-wing propulsion system layout is known as Over-The-Wing propulsion, which promises favorable aerodynamic effects by the rotor-induced flow and reduced noise signature during fly-over by shielding. However, during the operation of Over-The-Wing propulsion in eVTOL flight conditions, several multi-rotor-surface interactions are encountered, resulting in unexplored aerodynamic and aeroacoustic effects. As a consequence, design guidelines to maximize aerodynamic performance and minimize noise signature for Over-The-Wing propulsion for eVTOL flight conditions are missing in the public domain. This thesis focuses on enabling Over-The-Wing propulsion for vertical flight through experimentation and low-order modeling of the fundamental aerodynamic interactions between distributed rotors and surfaces. ...
The high power-to-weight ratio and efficiency of electric motors provide an excellent platform to explore disruptive propulsion system configurations. An observable design trend here is the tight integration of distributed, fixed-pitched rotors with aerodynamic surfaces, aimed at generating beneficial aerodynamic coupling effects to increase the aircraft’s efficiency and reduce acoustic emissions.
An auspicious integrated tilt-wing propulsion system layout is known as Over-The-Wing propulsion, which promises favorable aerodynamic effects by the rotor-induced flow and reduced noise signature during fly-over by shielding. However, during the operation of Over-The-Wing propulsion in eVTOL flight conditions, several multi-rotor-surface interactions are encountered, resulting in unexplored aerodynamic and aeroacoustic effects. As a consequence, design guidelines to maximize aerodynamic performance and minimize noise signature for Over-The-Wing propulsion for eVTOL flight conditions are missing in the public domain. This thesis focuses on enabling Over-The-Wing propulsion for vertical flight through experimentation and low-order modeling of the fundamental aerodynamic interactions between distributed rotors and surfaces.
A stochastic spectral correlation analysis between wall-pressure fluctuations and velocity fluctuations in the logarithmic of a turbulent pipe flow reveal Reynolds-number– independence of the wall-pressure linear coherence spectrum. This is a first-of-its-kind result, hinting at the feasibility of scaling an input sensing strategy based on wall-pressure fluctuations from a low-Reynolds-number environment to operational engineering conditions.
An initial controller based on wall-shear stress fluctuations was developed to target drag-producing large-scale structures in the logarithmic region. The flow response was measured in terms of both the statistical (and spectral) response of the TBL flow to realtime control and in terms of the effect the control has not only on the friction coefficient, but also on the integral measures. This analysis revealed three main findings: (1) an attenuation of energy at streamwise wavelengths characteristic of large-scale motions, (2) a decrease in skin-friction and (3) an attenuation of the statistical integral measures of skin-friction (i.e. bulk production and FIK terms).
A similar control architecture was also implemented that employed wall-pressure (and wall-pressure–squared) as the input quantity. It was found that the wall-pressure–squared term improves the accuracy of an estimator enabling the prediction of off-the-wall velocity fluctuations from a wall-based position. Furthermore, its inclusion is essential, given that the linear term does not retain sufficient coherence over the relatively large streamwise extent separating input and actuation locations.
The final controller that was developed in the context of this dissertation is an adaptive one, relying on the Filtered-X Least Means Squares (Fx-LMS) algorithm. This strategy does not rely on a-priori system identification, as was the case for the previous two control strategies mentioned above. Instead, it automatically identifies the coefficients of the transfer functions relating input to output in the controller. For this study, this algorithm was deployed both to a flow case that was strongly modulated by cylinder vortex shedding and to a fully broadband turbulent boundary layer flow. In the former case, the controller readily identified the shedding frequency as the control target. For the latter, the controller converges to a situation where the large-scales were targeted and their intensity successfully attenuated.
...
A stochastic spectral correlation analysis between wall-pressure fluctuations and velocity fluctuations in the logarithmic of a turbulent pipe flow reveal Reynolds-number– independence of the wall-pressure linear coherence spectrum. This is a first-of-its-kind result, hinting at the feasibility of scaling an input sensing strategy based on wall-pressure fluctuations from a low-Reynolds-number environment to operational engineering conditions.
An initial controller based on wall-shear stress fluctuations was developed to target drag-producing large-scale structures in the logarithmic region. The flow response was measured in terms of both the statistical (and spectral) response of the TBL flow to realtime control and in terms of the effect the control has not only on the friction coefficient, but also on the integral measures. This analysis revealed three main findings: (1) an attenuation of energy at streamwise wavelengths characteristic of large-scale motions, (2) a decrease in skin-friction and (3) an attenuation of the statistical integral measures of skin-friction (i.e. bulk production and FIK terms).
A similar control architecture was also implemented that employed wall-pressure (and wall-pressure–squared) as the input quantity. It was found that the wall-pressure–squared term improves the accuracy of an estimator enabling the prediction of off-the-wall velocity fluctuations from a wall-based position. Furthermore, its inclusion is essential, given that the linear term does not retain sufficient coherence over the relatively large streamwise extent separating input and actuation locations.
The final controller that was developed in the context of this dissertation is an adaptive one, relying on the Filtered-X Least Means Squares (Fx-LMS) algorithm. This strategy does not rely on a-priori system identification, as was the case for the previous two control strategies mentioned above. Instead, it automatically identifies the coefficients of the transfer functions relating input to output in the controller. For this study, this algorithm was deployed both to a flow case that was strongly modulated by cylinder vortex shedding and to a fully broadband turbulent boundary layer flow. In the former case, the controller readily identified the shedding frequency as the control target. For the latter, the controller converges to a situation where the large-scales were targeted and their intensity successfully attenuated.
The numerical analysis involves the application of aerodynamics, thermal management and performance techniques to model and simulate the complex airflow within the cooling ducts. Parameters such as duct geometry, airflow properties, and heat transfer rates are systematically investigated to optimise the cooling process. The study also explores the interaction between the cooling ducts and the overall performance of the aircraft, considering the impact on drag and fuel efficiency.
...
The numerical analysis involves the application of aerodynamics, thermal management and performance techniques to model and simulate the complex airflow within the cooling ducts. Parameters such as duct geometry, airflow properties, and heat transfer rates are systematically investigated to optimise the cooling process. The study also explores the interaction between the cooling ducts and the overall performance of the aircraft, considering the impact on drag and fuel efficiency.
Pulsed jet arrays for turbulent separation control
An experimental study
Both steady blowing and periodic pulsing actuation strategies are tested and compared. Preliminary jet velocity and pulsing frequency sweeps are carried out to identify optimal actuator operating parameters, relying on wall static pressure measurements to evaluate control effectiveness. Select cases of interest are then investigated using two-dimensional two-component particle image velocimetry and compared against the uncontrolled baseline which is characterized using PIV and hot wire anemometry. Additional PIV-derived metrics are utilized to assess system performance.
For steady blowing, a jet-to-crossflow velocity ratio VR>1 was required to produce a separation delay, while diminishing improvements in control effect with increasing jet velocity started at VR=1.6 (actuation momentum ratio of Cμ=2.3%). This nominal velocity ratio was adopted for all further investigation. The actuator was found to produce alternating strong and weak downwash regions in the TBL resulting in an artificial sweep/ejection pattern at detachment. Periodic forcing with the same nominal velocity ratio was able to achieve better or comparable results to steady actuation, while requiring less input momentum (Cμ=1.2-1.8%). The optimum actuation frequency was determined to be the natural frequency of the uncontrolled bubble, with the performance of higher frequency actuation tending towards steady blowing levels. As shown by an analysis of flow dynamics based on phase-averaged PIV velocity fields, actuation at the bubble time scales produces significant flow oscillation in phase with actuation. This resonant behaviour results in transient high momentum sweeps between actuation pulses that boost actuator performance, achieving double the performance benefit afforded by steady actuation according to multiple metrics. In comparison, actuation at time scales multiple times shorter than that of the bubble produces a quasi-steady flow and performance comparable to that of steady actuation.
Additionally, a novel alternating actuation strategy is tested, in which the period of active blowing is composed of high frequency alternation between two inverted actuator rows. This aimed to produce a quasi-2D periodic control effect using 3D actuators, which Squire's theorem suggests could excite the separated shear layer instability more than conventional 3D perturbation. While high frequency alternation did achieve a quasi-2D effect, it also prevented the sweep/ejection pattern characteristic of 3D perturbation from forming, thus significantly limiting the actuator performance. ...
Both steady blowing and periodic pulsing actuation strategies are tested and compared. Preliminary jet velocity and pulsing frequency sweeps are carried out to identify optimal actuator operating parameters, relying on wall static pressure measurements to evaluate control effectiveness. Select cases of interest are then investigated using two-dimensional two-component particle image velocimetry and compared against the uncontrolled baseline which is characterized using PIV and hot wire anemometry. Additional PIV-derived metrics are utilized to assess system performance.
For steady blowing, a jet-to-crossflow velocity ratio VR>1 was required to produce a separation delay, while diminishing improvements in control effect with increasing jet velocity started at VR=1.6 (actuation momentum ratio of Cμ=2.3%). This nominal velocity ratio was adopted for all further investigation. The actuator was found to produce alternating strong and weak downwash regions in the TBL resulting in an artificial sweep/ejection pattern at detachment. Periodic forcing with the same nominal velocity ratio was able to achieve better or comparable results to steady actuation, while requiring less input momentum (Cμ=1.2-1.8%). The optimum actuation frequency was determined to be the natural frequency of the uncontrolled bubble, with the performance of higher frequency actuation tending towards steady blowing levels. As shown by an analysis of flow dynamics based on phase-averaged PIV velocity fields, actuation at the bubble time scales produces significant flow oscillation in phase with actuation. This resonant behaviour results in transient high momentum sweeps between actuation pulses that boost actuator performance, achieving double the performance benefit afforded by steady actuation according to multiple metrics. In comparison, actuation at time scales multiple times shorter than that of the bubble produces a quasi-steady flow and performance comparable to that of steady actuation.
Additionally, a novel alternating actuation strategy is tested, in which the period of active blowing is composed of high frequency alternation between two inverted actuator rows. This aimed to produce a quasi-2D periodic control effect using 3D actuators, which Squire's theorem suggests could excite the separated shear layer instability more than conventional 3D perturbation. While high frequency alternation did achieve a quasi-2D effect, it also prevented the sweep/ejection pattern characteristic of 3D perturbation from forming, thus significantly limiting the actuator performance.
Motivated by the importance of understanding how the interaction between the developing flow and the vibrating nozzle walls has an effect on supersonic noise generation and propagation, this work has studied the effect that wall compliancy has on the vibroacoustic loading of TOP contoured nozzles. This is demonstrated by means of cold flow tests carried out in the High Speed Laboratories of the Delft University of Technology on a stiff-walled aluminum nozzle, which serves as a baseline test case, and on a urethane-based compliant walled nozzle. Tests are conducted under comparable flow conditions and test parameters are measured by means of acoustic and optical techniques. Simultaneous recordings are performed and include the nozzle-wall deformation, by means of stereoscopic tracking of tracers on the nozzle lip, the imprint of the near-field acoustic signature, by means of arrays of pressure-microphones, and Schlieren imaging of the jet plume.
Measurement data allows for a Fourier decomposition of the nozzle lip displacement and of the acoustic pressure field in azimuth. Reconstruction of the instantaneous plume development enables the identification of the main flow structures responsible for noise generation.
Comparison of results between the two test articles highlights a different spectral content and directivity pattern. Correlation between the structural displacements and the acoustic signal, together with the use of DMD, quantitatively aids the investigation of how FSI has an impact on the generation of an aeroelastic tone at 180 Hz. Findings suggest that its production is the result of the periodic thickening and thinning of the shear layer owing to the heightened flapping motion of the nozzle lip preceding RSS transition, driven by an intensified shock foot instability. ...
Motivated by the importance of understanding how the interaction between the developing flow and the vibrating nozzle walls has an effect on supersonic noise generation and propagation, this work has studied the effect that wall compliancy has on the vibroacoustic loading of TOP contoured nozzles. This is demonstrated by means of cold flow tests carried out in the High Speed Laboratories of the Delft University of Technology on a stiff-walled aluminum nozzle, which serves as a baseline test case, and on a urethane-based compliant walled nozzle. Tests are conducted under comparable flow conditions and test parameters are measured by means of acoustic and optical techniques. Simultaneous recordings are performed and include the nozzle-wall deformation, by means of stereoscopic tracking of tracers on the nozzle lip, the imprint of the near-field acoustic signature, by means of arrays of pressure-microphones, and Schlieren imaging of the jet plume.
Measurement data allows for a Fourier decomposition of the nozzle lip displacement and of the acoustic pressure field in azimuth. Reconstruction of the instantaneous plume development enables the identification of the main flow structures responsible for noise generation.
Comparison of results between the two test articles highlights a different spectral content and directivity pattern. Correlation between the structural displacements and the acoustic signal, together with the use of DMD, quantitatively aids the investigation of how FSI has an impact on the generation of an aeroelastic tone at 180 Hz. Findings suggest that its production is the result of the periodic thickening and thinning of the shear layer owing to the heightened flapping motion of the nozzle lip preceding RSS transition, driven by an intensified shock foot instability.
Manipulation Of Large Scales Via A Spanwise Array Of Wall-Normal Jets In A Turbulent Boundary Layer
An Attempt At Large Scale Control For Drag Reduction
LSMs convecting in the logarithmic region (log region) of TBL play a critical role in the dynamics of wall-bounded turbulence, as they carry a significant portion of the turbulent kinetic energy (TKE) [Abbassi et al., 2017]. The active wall oscillation technique tuned to the frequency of the LSMs has proven to be energy-efficient, as drag reduction increases with the frictional Reynolds number (Reτ ) and the contribution of LSMs to the TKE and wall shear stress increases with the Reτ [Marusic et al., 2021]. In this thesis work, the receptivity of LSMs to active large-scale control strategy is assessed via a spanwise array of wall-normal jets in a TBL. Multiple wall-normal jets manipulate the flow over a flat plate by introducing a spanwise traveling wave, which aims to mimic wall oscillation tuned to the large scales in the log region of the TBL.
To assess the efficacy of the suggested control strategy with spatial and temporal tuning, particle image velocimetry (PIV) and hot-wire anemometry (HWA) are employed to quantify second-order turbulence statistics and analyze the organization of the coherent patterns introduced in the flow via two-point correlations. In addition, the effect of different control cases derived from different tuning and control strategies is analyzed at Reτ = 2227. The TKE production plots show that control cases targeting the LSMs effectively attenuate the TKE production near the wall. Furthermore, based on the analysis of the spectral energy plot of streamwise velocity fluctuations, it is observed that the energy content associated with larger length scales convecting in the middle of the log region decreases in the controlled case compared to the uncontrolled base case. This finding suggests that certain control cases notably impact the organization of LSM and second-order turbulence statistics far downstream of the actuation point. Additionally, the reduction in TKE production near the wall represents a reduction in wall shear stress and viscous drag. However, additional research must be conducted in the future based on the current findings to reach a definitive conclusion regarding the efficacy of individual control cases.
...
LSMs convecting in the logarithmic region (log region) of TBL play a critical role in the dynamics of wall-bounded turbulence, as they carry a significant portion of the turbulent kinetic energy (TKE) [Abbassi et al., 2017]. The active wall oscillation technique tuned to the frequency of the LSMs has proven to be energy-efficient, as drag reduction increases with the frictional Reynolds number (Reτ ) and the contribution of LSMs to the TKE and wall shear stress increases with the Reτ [Marusic et al., 2021]. In this thesis work, the receptivity of LSMs to active large-scale control strategy is assessed via a spanwise array of wall-normal jets in a TBL. Multiple wall-normal jets manipulate the flow over a flat plate by introducing a spanwise traveling wave, which aims to mimic wall oscillation tuned to the large scales in the log region of the TBL.
To assess the efficacy of the suggested control strategy with spatial and temporal tuning, particle image velocimetry (PIV) and hot-wire anemometry (HWA) are employed to quantify second-order turbulence statistics and analyze the organization of the coherent patterns introduced in the flow via two-point correlations. In addition, the effect of different control cases derived from different tuning and control strategies is analyzed at Reτ = 2227. The TKE production plots show that control cases targeting the LSMs effectively attenuate the TKE production near the wall. Furthermore, based on the analysis of the spectral energy plot of streamwise velocity fluctuations, it is observed that the energy content associated with larger length scales convecting in the middle of the log region decreases in the controlled case compared to the uncontrolled base case. This finding suggests that certain control cases notably impact the organization of LSM and second-order turbulence statistics far downstream of the actuation point. Additionally, the reduction in TKE production near the wall represents a reduction in wall shear stress and viscous drag. However, additional research must be conducted in the future based on the current findings to reach a definitive conclusion regarding the efficacy of individual control cases.
Study of high aspect-ratio dual intersecting jets and the installation effect on control authority
An experimental and numerical investigation
The intersecting twin jet has recently captured significant attention from the Aerothermal team for several reasons. Firstly, it provides the ability to actively control the flow angle from the HVAC unit into the vehicle cabin without requiring direct interaction with the vent. Additionally, it facilitates the flushing of the vent outlets with the A-class surfaces in the vehicle interior, helps achieve a simple and minimalist design for the instrument panel and other interior surfaces that align with Tesla's design language
The objective pursued in this thesis is to investigate the physics represented by the bulk flow characteristics of intersecting twin jets, and to scrutinize the installation effect on the merging jet, focusing on the influence of the installation surface offset height. This objective is achieved by characterizing the evolution of large-scale flow structures in High AR intersecting twin jets at moderate Reynolds numbers experimentally using planar Particle Image Velocimetry (PIV) in both installed and uninstalled conditions. Once accomplished, the experimental results are then used to help select, validate, and calibrate the appropriate RANS turbulence model.
The experimental results show that the intersecting jets converge before merging and forming a single jet. This jet behaves like a single jet emanating from a more recessed origin. The variation of the flow ratio between the two nozzle outlets results in varying the jet angle. The study reveals that the placement of the surface impacts the jet angle and velocity and introduces a circulation zone that affects the jet's behavior and direction. The surface's impact depends on the offset distance, length, and profile of the surface. The results also suggest that accurate computational prediction of the jet characteristics requires high mesh resolution and an appropriate selection of turbulence intensity. Among the four turbulence models assessed, the k-omega SST model is found to be adequate for the application, although it exhibits some numerical hysteresis. Overall, the findings provide insights into the dynamics of jet-surface interactions and can guide the design of systems involving such flows.
...
The intersecting twin jet has recently captured significant attention from the Aerothermal team for several reasons. Firstly, it provides the ability to actively control the flow angle from the HVAC unit into the vehicle cabin without requiring direct interaction with the vent. Additionally, it facilitates the flushing of the vent outlets with the A-class surfaces in the vehicle interior, helps achieve a simple and minimalist design for the instrument panel and other interior surfaces that align with Tesla's design language
The objective pursued in this thesis is to investigate the physics represented by the bulk flow characteristics of intersecting twin jets, and to scrutinize the installation effect on the merging jet, focusing on the influence of the installation surface offset height. This objective is achieved by characterizing the evolution of large-scale flow structures in High AR intersecting twin jets at moderate Reynolds numbers experimentally using planar Particle Image Velocimetry (PIV) in both installed and uninstalled conditions. Once accomplished, the experimental results are then used to help select, validate, and calibrate the appropriate RANS turbulence model.
The experimental results show that the intersecting jets converge before merging and forming a single jet. This jet behaves like a single jet emanating from a more recessed origin. The variation of the flow ratio between the two nozzle outlets results in varying the jet angle. The study reveals that the placement of the surface impacts the jet angle and velocity and introduces a circulation zone that affects the jet's behavior and direction. The surface's impact depends on the offset distance, length, and profile of the surface. The results also suggest that accurate computational prediction of the jet characteristics requires high mesh resolution and an appropriate selection of turbulence intensity. Among the four turbulence models assessed, the k-omega SST model is found to be adequate for the application, although it exhibits some numerical hysteresis. Overall, the findings provide insights into the dynamics of jet-surface interactions and can guide the design of systems involving such flows.
To continue building upon this previous work, subscale rocket nozzles were poured using various types of polyurethane. To cast these nozzles, multiple mould parts were manufactured. To guarantee a smooth surface finish for the inner nozzle wall, an inner mould was CNC machined. To enable the casting of multiple nozzles with different wall thicknesses, outer moulds were 3D-printed. After having determined the correct variant of polyurethane, Smooth-On's 'Smooth-Cast 66D', three different nozzles were cast on which tests were conducted. The nominal nozzle was designed with a wall thickness of 2mm, with two additional nozzles being cast with a wall thickness of 1.5mm and 2.5mm.
Having manufactured the nozzles, they were mounted to the test facility. The test facility consists of a stand connected to an air tank capable of containing air at pressures of up to 40 bar. Due to pressure losses throughout the valves, the maximum attainable nozzle pressure ratio, or NPR, amounted to slightly over 30. Various tests were performed during which the NPR was either altered during the test, done by manually opening or closing the control valve contained in the system. Further tests were conducted at a number of constant values of NPR at which excessive vibrations occurred to gain further insight into the flow phenomena present at these critical values of NPR. Lip tracking was done as well as schlieren photography. For the lip tracking, correction fluid was applied to the nozzle lip as well as spray painting the outer wall of the nozzle black. The area surrounding the nozzle was made dark such that the predominant feature photographed would be the circular nozzle lip. Once the images were taken, the nozzle lip was discretised in 180 points and was compared to an image of the nozzle taken when the wind tunnel was closed. The end result is a discretised function of 180 points displaying the lip deflection at each azimuthal location. With the help of the discrete Fourier transform, the dominant eigenmodes could be extracted from the nozzle vibrations in the form of their Fourier coefficients. In addition to lip tracking, schlieren images were taken simultaneously with the aim of visualising the flow downstream of the nozzle exit. The brightness of a set of pixels was tracked and their spectral make-up was determined.
Along the entire scale of NPR, two regions were found in which most of the vibrations occurred. The first region of activity, was purported to be linked to the transition of separation regime from free shock separation, FSS, to restricted shock separation, RSS. The second region of NPR, in which the deflection was often even more excessive than the first, is highly likely due to the presence of a recirculation bubble passing over the nozzle lip. Out of the first three eigenmodes, named the breathing, bending and ovalisation mode respectively, the ovalisation mode accounted for the largest proportion of the deflection by far. Based on these transient tests, further tests were conducted at a constant NPR of NPR = 21.7 and NPR = 27.7. These were chosen as all three nozzles showed large deflections at these NPR. For the nominal 2mm nozzle, two additional values were chosen of NPR = 22.7 and NPR = 27.5 to conduct further tests at. The energy present in each of the first three eigenmodes was determined by calculating the variance of the Fourier coefficients. The energy in the modes during the constant NPR runs were compared to the energy in the modes during the transient tests. These were obtained by determining the energy of an NPR-envelop spanning 1 unit of NPR centered around the NPR at which the constant run was performed. By and large, the energy storage in the first three modes remained similar when comparing transient tests with constant NPR runs. The bending mode showed the largest deviation from this trend, mainly when in RSS. Frequency analyses showed the nozzle to oscillate at rates predominantly within 130Hz-400Hz. Wavelet power spectra showed no clear trend between the NPR and the frequency with which the nozzles vibrated. Strouhal calculations further confirmed that NPR would not have a drastic effect on the unsteadiness frequency of the flow. The testing helped provide insight into the mechanical vibrations of nozzles of various wall thicknesses in different flow conditions. However, due to the absence of quantitative flow visualisation and measured flow unsteadiness frequency, no concrete link could be made to the aeroelastic coupling and the fluid's role in the forcing of the vibrations. ...
To continue building upon this previous work, subscale rocket nozzles were poured using various types of polyurethane. To cast these nozzles, multiple mould parts were manufactured. To guarantee a smooth surface finish for the inner nozzle wall, an inner mould was CNC machined. To enable the casting of multiple nozzles with different wall thicknesses, outer moulds were 3D-printed. After having determined the correct variant of polyurethane, Smooth-On's 'Smooth-Cast 66D', three different nozzles were cast on which tests were conducted. The nominal nozzle was designed with a wall thickness of 2mm, with two additional nozzles being cast with a wall thickness of 1.5mm and 2.5mm.
Having manufactured the nozzles, they were mounted to the test facility. The test facility consists of a stand connected to an air tank capable of containing air at pressures of up to 40 bar. Due to pressure losses throughout the valves, the maximum attainable nozzle pressure ratio, or NPR, amounted to slightly over 30. Various tests were performed during which the NPR was either altered during the test, done by manually opening or closing the control valve contained in the system. Further tests were conducted at a number of constant values of NPR at which excessive vibrations occurred to gain further insight into the flow phenomena present at these critical values of NPR. Lip tracking was done as well as schlieren photography. For the lip tracking, correction fluid was applied to the nozzle lip as well as spray painting the outer wall of the nozzle black. The area surrounding the nozzle was made dark such that the predominant feature photographed would be the circular nozzle lip. Once the images were taken, the nozzle lip was discretised in 180 points and was compared to an image of the nozzle taken when the wind tunnel was closed. The end result is a discretised function of 180 points displaying the lip deflection at each azimuthal location. With the help of the discrete Fourier transform, the dominant eigenmodes could be extracted from the nozzle vibrations in the form of their Fourier coefficients. In addition to lip tracking, schlieren images were taken simultaneously with the aim of visualising the flow downstream of the nozzle exit. The brightness of a set of pixels was tracked and their spectral make-up was determined.
Along the entire scale of NPR, two regions were found in which most of the vibrations occurred. The first region of activity, was purported to be linked to the transition of separation regime from free shock separation, FSS, to restricted shock separation, RSS. The second region of NPR, in which the deflection was often even more excessive than the first, is highly likely due to the presence of a recirculation bubble passing over the nozzle lip. Out of the first three eigenmodes, named the breathing, bending and ovalisation mode respectively, the ovalisation mode accounted for the largest proportion of the deflection by far. Based on these transient tests, further tests were conducted at a constant NPR of NPR = 21.7 and NPR = 27.7. These were chosen as all three nozzles showed large deflections at these NPR. For the nominal 2mm nozzle, two additional values were chosen of NPR = 22.7 and NPR = 27.5 to conduct further tests at. The energy present in each of the first three eigenmodes was determined by calculating the variance of the Fourier coefficients. The energy in the modes during the constant NPR runs were compared to the energy in the modes during the transient tests. These were obtained by determining the energy of an NPR-envelop spanning 1 unit of NPR centered around the NPR at which the constant run was performed. By and large, the energy storage in the first three modes remained similar when comparing transient tests with constant NPR runs. The bending mode showed the largest deviation from this trend, mainly when in RSS. Frequency analyses showed the nozzle to oscillate at rates predominantly within 130Hz-400Hz. Wavelet power spectra showed no clear trend between the NPR and the frequency with which the nozzles vibrated. Strouhal calculations further confirmed that NPR would not have a drastic effect on the unsteadiness frequency of the flow. The testing helped provide insight into the mechanical vibrations of nozzles of various wall thicknesses in different flow conditions. However, due to the absence of quantitative flow visualisation and measured flow unsteadiness frequency, no concrete link could be made to the aeroelastic coupling and the fluid's role in the forcing of the vibrations.
Chevron-shaped protrusions for turbulent drag reduction
Experimental investigation into the drag performance and flow mechanics
In this study, the experiments from the literature are replicated, and new array configurations are tested to characterise the effect of individual parameters on the drag performance. The test plates are manufactured by applying vinyl protrusions of roughly 100 μm in thickness to an aluminium base plate. This thickness corresponds to 5δν - 6δν for the design Reynolds number of Reτ= 1270. Direct force measurements are performed in the M-tunnel at a Reynolds number range of approximately 630 < Reτ < 1850 to determine the drag performance. Furthermore, the coherent structures are characterised by means of 2D-2C PIV of a wall-parallel plane at a minimum distance of 17δν from the wall.
The drag reduction reported in the literature could not be replicated, and the balance measurement results offer relevant insights into the drag performance of chevron-shaped protrusions. The results consistently show that the added roughness due to the presence of the protrusions is not the only parameter that determines the drag performance, confirming a meaningful interaction between the protrusions and the flow. Moreover, the results are found to be highly sensitive to the randomisation of the array. No substantial effect of the protrusions on the coherent structures close to the wall has been observed. In particular, no evidence has been observed that supports the working mechanism as proposed in the literature.
An analysis of possible causes to explain the discordance between this study and the literature is performed. Based on this analysis and the results from the aforementioned parametric study, an improved design is proposed, and recommendations for future research are postulated. This technology has inherent benefits for real-world implementations as it can easily be (retro)fitted to aircraft by means of a foil. Further research into this flow control technique is thereby deemed relevant due to the combination of the large drag reduction reported in the literature, the advantages in practical applications, and the novel opportunities for additional investigations. ...
In this study, the experiments from the literature are replicated, and new array configurations are tested to characterise the effect of individual parameters on the drag performance. The test plates are manufactured by applying vinyl protrusions of roughly 100 μm in thickness to an aluminium base plate. This thickness corresponds to 5δν - 6δν for the design Reynolds number of Reτ= 1270. Direct force measurements are performed in the M-tunnel at a Reynolds number range of approximately 630 < Reτ < 1850 to determine the drag performance. Furthermore, the coherent structures are characterised by means of 2D-2C PIV of a wall-parallel plane at a minimum distance of 17δν from the wall.
The drag reduction reported in the literature could not be replicated, and the balance measurement results offer relevant insights into the drag performance of chevron-shaped protrusions. The results consistently show that the added roughness due to the presence of the protrusions is not the only parameter that determines the drag performance, confirming a meaningful interaction between the protrusions and the flow. Moreover, the results are found to be highly sensitive to the randomisation of the array. No substantial effect of the protrusions on the coherent structures close to the wall has been observed. In particular, no evidence has been observed that supports the working mechanism as proposed in the literature.
An analysis of possible causes to explain the discordance between this study and the literature is performed. Based on this analysis and the results from the aforementioned parametric study, an improved design is proposed, and recommendations for future research are postulated. This technology has inherent benefits for real-world implementations as it can easily be (retro)fitted to aircraft by means of a foil. Further research into this flow control technique is thereby deemed relevant due to the combination of the large drag reduction reported in the literature, the advantages in practical applications, and the novel opportunities for additional investigations.
A parameter study on HRs under a fully developed grazing TBL at Reτ ≈ 2200 was performed. A strong pressure-velocity coupling was found between the HR and the grazing TBL when the HR design frequency either matched the frequency of the most energetic wall-normal velocity fluctuations or when the HR design frequency was below this frequency. The pressure-velocity coupling extends to y+ ≈ 25, in which phase-interlocking of the grazing TBL occurs. Clear inflow and outflow regions were identified with a streamwise width related to the HR resonance frequency. The inflow and outflow regions cause an increase and decrease in streamwise velocity, respectively. These velocity fluctuations appear as an increase in spectral energy around the wavelengths of the HR resonance and are accompanied by a reduction in spectral energy at higher wavelengths. While quadrant analysis indicates an increase in turbulence production from an increased relative strength of Q2 and Q4 events, variable-interval time averaging indicates no significant changes to the intensity and duration of the near-wall cycle turbulence events as a result of the HR resonance. To date, no significant changes to the mean boundary layer statistics were found that could directly be attributed to the achieved pressure-velocity coupling. Note that only a single HR has been the focus of the current study, with valuable information about HR scaling, resonance and domain-of influence on the TBL flow. Future arrays of HRs may be able to show a more pronounced global effect on the mean flow.
...
A parameter study on HRs under a fully developed grazing TBL at Reτ ≈ 2200 was performed. A strong pressure-velocity coupling was found between the HR and the grazing TBL when the HR design frequency either matched the frequency of the most energetic wall-normal velocity fluctuations or when the HR design frequency was below this frequency. The pressure-velocity coupling extends to y+ ≈ 25, in which phase-interlocking of the grazing TBL occurs. Clear inflow and outflow regions were identified with a streamwise width related to the HR resonance frequency. The inflow and outflow regions cause an increase and decrease in streamwise velocity, respectively. These velocity fluctuations appear as an increase in spectral energy around the wavelengths of the HR resonance and are accompanied by a reduction in spectral energy at higher wavelengths. While quadrant analysis indicates an increase in turbulence production from an increased relative strength of Q2 and Q4 events, variable-interval time averaging indicates no significant changes to the intensity and duration of the near-wall cycle turbulence events as a result of the HR resonance. To date, no significant changes to the mean boundary layer statistics were found that could directly be attributed to the achieved pressure-velocity coupling. Note that only a single HR has been the focus of the current study, with valuable information about HR scaling, resonance and domain-of influence on the TBL flow. Future arrays of HRs may be able to show a more pronounced global effect on the mean flow.