This thesis is a study of the effects of distributed roughness elements on a two-dimensional turbulent boundary layer. Measurements were taken on a total of ten rough wall configurations: four involving Gaussian spikes, and six with circular cylindrical posts. Rough wall flows are particularly suited to study with Laser Doppler Velocimetry (LDV) due to the fact that measurements are required near a solid surface, as well has in highly turbulent fluid. The LDV system used in this study is a fine resolution (~50 micron), three-component, fiber optic system. All mean velocities, Reynolds stresses, and triple products are measured. This study is unique in the range and variety of roughness cases for which data was taken. The data show that the flow over a rough wall is characterized by high levels of turbulence near the roughness element peaks at the interface between low-speed, near-wall fluid and the higher speed fluid above. Behind an element, high-momentum fluid sweeps toward the wall, and there is a small region of ejection of low-momentum fluid. Cylindrical elements typically have larger magnitudes of turbulent stresses at their peaks compared to Gaussian elements. Trends in mean velocity profile parameters such as displacement height, roughness effect, and wake parameter are examined with respect to roughness element geometry and spacing. / Master of Science

Turbulent flow and transition are some of the most important phenomena of fluid mechanics and aerodynamics and represent a challenging engineering problem for aircraft manufacturers looking to improve aerodynamic efficiency. Laminar flow technology has the potential to provide a significant reduction to aircraft drag by manipulating the instabilities within the laminar boundary layer to achieve a delay in transition to turbulence. Currently prediction and simulation of laminar-turbulent transition is con- ducted using either a low-fidelity approach involving the stability equations or via a full Direct Numerical Simulation (DNS). The work in this thesis uses an alternative high-fidelity simulation method that aims to bridge the gap between the two simulation streams. The methodology uses an LES approach with a low-computational cost sub-grid scale model (WALE) that has inherent ability to reduce its turbulent viscosity contribution to zero in laminar regions. With careful grid spacing the laminar regions can be explicitly modelled as an unsteady Navier-Stokes simulation while the turbulent and transitional regions are simulated using LES. The methodology has been labelled as an unsteady Navier-Stokes/Large Eddy Simulation (UNS/LES) approach. Two test cases were developed to test the applicability of the method to simulate and control the crossflow instability. The first test case replicated the setup from an experiment that ran at a chord-based Reynolds number of 390, 000. Two methods were used to generate the initial disturbance for the crossflow vortices, firstly using a continuous suction hole and secondly an isolated roughness element. The results for this test case showed that the approach was capable of modelling the full transition process, from explicitly modelling the growth of the initial amplitude of the disturbances to final breakdown to turbulence. Results matched well with the available experimental data. The second test case replicated an experimental setup using a custom- designed aerofoil run at a chord-based Reynolds number of 2.4 million. The test case used Distributed Roughness Elements (DRE) to induce crossflow vortices at both a critical and a control wavelength. By forcing the crossflow vortices at a stable (control) wavelength a delay in laminar-turbulent transition can be achieved. The results showed that the UNS/LES approach was capable of capturing the initial disturbance amplitudes due to the roughness elements and their growth rates matched well with experimental data. Finally, downstream a transitional region was assessed with low-freestream turbulence provided using a modified Synthetic Eddy Method (SEM). The full laminar-turbulent transition pro- cess was simulated and results showed significant promise. In conclusion, the method employed in this thesis showed promising results and demonstrated a possible route to high-fidelity transition simulation run at more realistic flow conditions and geometries than DNS. Further work and validation is required to test the secondary instability region and the final breakdown to turbulence.

Modeling the effects of surface roughness is an area of concern in many practical engineering applications. Many current roughness models to this point have involved the use of empirical 'constants' and equivalent sand grain roughness. These underdeveloped concepts have little direct relationship to realistic roughness and cannot predict accurately and consistently the flow characteristics for different roughness shapes. In order to aid in the development of turbulence models, the present research is centered around the experimental investigation of seven various shaped single roughness elements and their effects on turbulence quantities in a two-dimensional turbulent boundary layer. The elements under scrutiny are as follows: cone, cone with spatial variations equal to the smallest sublayer structure length scale, cone with spatial variations equal to 2.5 times the smallest sublayer structure length scale, Gaussian-shaped element, hemisphere, cube aligned perpendicular to the flow (cube at 90°), and a cube rotated 45° relative to the flow. The roughness element heights, k+, non-dimensionalized by the friction velocity (U_tau) of the approaching turbulent boundary layer, are 145, 145, 145, 145, 80, 98, and 98 respectively. Analysis of a three-dimensional fetch of the same Gaussian-shaped elements described previously was also undertaken. In order to analyze the complex flow fields, detailed measurements were obtained using a fine-measurement-volume (50 micron diameter) three-velocity component laser-Doppler velocimetry (LDV) system. The data reveals the formation of a horseshoe vortex in front of the element, which induces the downwash of higher momentum fluid toward the wall. This 'sweep' motion not only creates high Reynolds stresses (v^2, w^2, -uv) downstream of the element, but also leads to higher skin-friction drag. Triple products were also found to be very significant near the height of the element. These parameters are important in regards to the contribution of the production and diffusion of the turbulent kinetic energy in the flow. The 'peakiness' of the roughness element was found to have a direct correlation to the production of circulation, whereas the spatial smoothing does not have an immense effect on this parameter. The peaked elements were found to have a similar trend in the decay of circulation in the streamwise direction. These elements tend to show a decay proportional to (x/d)^-1.12, whereas the cube elements and the hemisphere do not have a common trend. A model equation is proposed for a drag correlation common to all roughness elements. This equation takes into account the viscous drag and pressure drag terms in the calculation of the actual drag due to the roughness elements presence in the boundary layer. The size, shape, frontal and wetted surface areas of the roughness elements are related to one another via this model equation. Flow drawings related to each element are presented which gives rise to a deeper understanding of the physics of the flow associated with each roughness element. / Master of Science