This thesis deals with receptivity mechanisms of three-dimensional, crossflow-dominated boundary layers. The receptivity of two model problems, a swept-flat-plate and a swept-wing boundary layer, is investigated by solving the parabolised stability equations (PSE) as well as by performing direct numerical simulations (DNS).Both flow cases are known to exhibit strong inflectional instabilities, the crossflow disturbances, whose excitation by external disturbances such as surface roughness or free-stream vorticity is studied. One focus is on worst-case scenarios. This involves the determination of optimal conditions, i.e. those disturbance environments yielding the largest possible response inside the boundary layer. A new method on the basis of the PSE is presented which allows to study optimal disturbances of swept-flat-plate boundary layers. These take the form of tilted streamwise vortices. While convected downstream they develop into streamwise streaks experiencing strong non-modal growth. Eventually, they turn into crossflow disturbances and undergo exponential growth. Non-modal growth is thus found to optimally excite crossflow disturbances and can be related to a receptivity mechanism of three-dimensional boundary layers. Evaluating effects of compressibility reveals that the potential for both non-modal and modal growth increases for higher Mach numbers. It is shown that wall cooling has diverse effects on disturbances of non-modal and modal nature. While destabilising the former it attenuates the growth of modal disturbances. Concave curvature on the other hand is found to be equally destabilising for both types of disturbances. The adjoint of the linearised Navier-Stokes equations is solved for a swept-wing boundary layer by means of DNS. The adjoint solution of a steady crossflow disturbance is computed in the boundary layer as well as in the free-stream upstream of the leading edge. This allows to determine receptivity to incoming free-stream disturbances and surface roughness as well as the corresponding worst-case scenarios. Upstream of a swept wing the optimal initial free-stream disturbance is found to be of streak-type which convects downstream towards the leading edge. It entrains the boundary layer a short distance downstream of the stagnation line. While minor streamwise vorticity is present the streak component is dominant all the way into the boundary layer where the optimal disturbance turns into a crossflow mode. Futher, the worst-case surface roughness is determined. It takes a wavy shape and is distributed in the chordwise direction. It is shown that, under such optimal conditions, the swept-wing boundary layer is more receptive to surface roughness than to free-stream disturbances. Another focus of this work has been the development and evaluation of tools for receptivity prediction. Both DNS and direct and adjoint solutions of the PSE are used to predict the receptivity of a swept-wing boundary layer to localised surface roughness. The configuration conforms to wind tunnel experiments performed by Saric and coworkers at the Arizona State University. Both the DNS and the PSE are found to predict receptivity amplitudes which are in excellent agreement with each other. Though the predicted disturbance amplitudes are slightly lower than experimental measurements the overall agreement with experimental results is very satisfactory. Finally, a DNS of the stabilisation of a transitional swept-wing boundary layer by means of discrete roughness elements is presented. This control approach is found to completely suppress transition to turbulence within the domain studied and confirms experimental results by Saric & coworkers. / QC 20111124

The correct tribological design will have a considerable effect on a gear’s service life and efficiency. The purpose of this thesis is to clarify the impact of variation in the gear tooth flank tribological system on the gear contact load capacity – to increase the understanding of how surface topography and lubricant interact. In this thesis the variation in surface topography inherent in the manufacturing method has been shown, by experimental work and computer simulations, to be an important factor for the contact condition in the early life of gears. Surface analysis revealed that the formation and composition of surface boundary layers depends strongly on the chemical composition of the lubricant, but also on pre-existing surface boundary layers. Additionally, surface boundary layers play a major role in frictional behaviour, wear and in allowing the lubricant to react properly with the surfaces. Paper A presents the current ISO 6336 calculation of surface durability. A robust design approach was used to investigate the extent to which the current standard for calculation of surface durability allows for manufacturing variations and the choice of lubricant. Paper B investigates the extent to which a logarithmical profile modification can increase gear contact pressure robustness compared to traditional lead profiles for gears. Paper C compares different gear manufacturing methods and their as-manufactured (fresh unworn) surface topographies, using measured surface topographies as input to a contact simulation program. Paper D examines surface boundary layer formation and the corresponding wear in relation to different anti-wear additives in an environmentally adapted base oil. Papers E and F make use of specimens with surface topographies imitating two gear manufacturing methods (grinding and superfinishing) to be used in a twin-disc and barrel-on-disc machine respectively. The contacts are analysed by friction measurements and simulations combined with methods for surface analysis. / <p>QC 20120925</p>

A new large eddy simulation (LES) code is developed and used to investigate non-local features of turbulent planetary boundary layers (PBLs). The LES code is based on filtered Navier-Stokes equations, which describe motions of incompressible, Boussinesq fluid at high Reynolds numbers. The code computes directly large-scale, non-universal turbulence in the PBL whereas small-scale, universal turbulence is parameterized by a dynamic mixed subgrid closure. The LES code is thoroughly tested against high quality laboratory and field data. This study addresses non-local properties of turbulence which emphasis on the stable stratification. Its basic results are as follows. The flow stability in PBLs is generally caused by two mechanisms: the negative buoyancy force (in the stable density stratification) and the Coriolis force (in the rotating system). The latter stabilizes the flow if the earth’s vorticity and the turbulent vorticity are anti-parallel. The Coriolis force stability suppresses large-scale turbulence and makes large eddies asymmetric. The density stratification suppresses vertical scales of turbulence. Joint actions of the Coriolis and the buoyancy forces result in a more complex behavior of turbulence. Particularly, the layers of vigorous turbulence may appear in the course of development of low-level jets in baroclinic atmosphere. Non-local effects determine integral measures of PBLs, first of all the PBL depth. This study clearly demonstrates its pronounced dependences on the Coriolis parameter, the Kazanski-Monin internal stability parameter, and newly introduced imposed-stability and baroclinicity parameters. An LES database is created and used to validate an advanced PBL-depth formulation. LES support the idea that PBLs interact with the stably stratified free flow through the radiation of gravity waves, excited by large turbulent eddies at the interface.

The boundary layer is the part of the atmosphere that is in direct contact with the ground via turbulent motion. At mid-latitudes the boundary layer is usually one or a few kilometers deep, while in the Arctic it is much more shallow, typically a few hundred meters or less. The reason is that here the absolute temperature increases in the lowest kilometer, making the boundary layer semi-permanently stably stratified. The exchange of heat, momentum and tracers between the atmosphere, ocean and ground under stable stratification is discussed from an observational, modeling and climate-change point of view. A compilation of six observational datasets, ordered by the Richardson number (rather than the widely used Monin-Obukhov length) reveals new information about turbulence in the very stably stratified regime. An essentially new turbulence closure model, based on the total turbulent energy concept and these observational datasets, is developed and tested against large-eddy simulations with promising results. The role of mesoscale motion in the exchange between the atmosphere and surface is investigated both for observations and in idealized model simulations. Finally, it is found that the stably stratified boundary layer is more sensitive to external surface forcing than its neutral and convective counterparts. It is speculated that this could be part of the explanation for the observed Arctic amplification of climate change.

Physical experiments are used to investigate the influence of the Coriolis forces on flow structures in channelized turbidity currents, and their implication for the evolution of straight and sinuous submarine channels. Initial tests were used to determine whether or not saline density currents are a good surrogate for particle-laden currents. Results imply that this assumption is valid when turbidity currents are weakly-depositional and have similar velocity and turbulence structures to saline density currents. Second, the controls of Coriolis forces on flow structures in straight channel sections are compared with two mathematical models: Ekman boundary layer dynamics and the theory of Komar [1969]. Ekman boundary layer dynamics prove to be a more suitable description of flow structures in rotating turbidity currents and should be used to derive flow parameters from submarine channels systems that are subjected to Coriolis forces. The significance of Coriolis forces for submarine channel systems were determined by evaluating the dimensionless Rossby number RoW. The Rossby number is defined as the ratio of the flow velocity, U, of a turbidity current to the channel width, W, and the rotation rate of the Earth represented by the Coriolis parameter, f. Coriolis forces are very significant for channel systems with RoW ≤ O(1). Third, the effect of Coriolis forces on the internal flow structure in sinuous submarine channels is considered. Since previous studies have only considered pressure gradient and centrifugal forces, the Coriolis force provides a crucial contribution to the lateral momentum balance in channel bends. In a curved channel, both the Rossby number RoW and the ratio of the channel curvature radius R to the channel width W, determine whether Coriolis forces affect the internal flow structure. The results demonstrate that Coriolis forces can cause a significant shift of the density interface and the downstream velocity core of channelized turbidity currents. The sediment transport regime in high-latitude channel systems, which have RoW << R/W, is therefore strongly influenced by Coriolis forces. Finally, these findings are incorporated into a conceptual model describing the evolution of submarine channels at different latitudes. For instance, the Northern Hemisphere channels have a distinctly higher right levee system and migrate predominantly to the left side and generally exhibit a low sinuosity. In contrast, low latitude channel systems have RoW >> R/W so that centrifugal forces are more dominant. This results in more sinuous submarine channel systems with varying levee asymmetries in subsequent channel bends. In conclusion, Coriolis forces are negligible around the equator but should be considered in high latitude systems, particularly when RoW ~ O(1) and RoW << R/W.

Physical experiments are used to investigate the influence of the Coriolis forces on flow structures in channelized turbidity currents, and their implication for the evolution of straight and sinuous submarine channels. Initial tests were used to determine whether or not saline density currents are a good surrogate for particle-laden currents. Results imply that this assumption is valid when turbidity currents are weakly-depositional and have similar velocity and turbulence structures to saline density currents. Second, the controls of Coriolis forces on flow structures in straight channel sections are compared with two mathematical models: Ekman boundary layer dynamics and the theory of Komar [1969]. Ekman boundary layer dynamics prove to be a more suitable description of flow structures in rotating turbidity currents and should be used to derive flow parameters from submarine channels systems that are subjected to Coriolis forces. The significance of Coriolis forces for submarine channel systems were determined by evaluating the dimensionless Rossby number RoW. The Rossby number is defined as the ratio of the flow velocity, U, of a turbidity current to the channel width, W, and the rotation rate of the Earth represented by the Coriolis parameter, f. Coriolis forces are very significant for channel systems with RoW ≤ O(1). Third, the effect of Coriolis forces on the internal flow structure in sinuous submarine channels is considered. Since previous studies have only considered pressure gradient and centrifugal forces, the Coriolis force provides a crucial contribution to the lateral momentum balance in channel bends. In a curved channel, both the Rossby number RoW and the ratio of the channel curvature radius R to the channel width W, determine whether Coriolis forces affect the internal flow structure. The results demonstrate that Coriolis forces can cause a significant shift of the density interface and the downstream velocity core of channelized turbidity currents. The sediment transport regime in high-latitude channel systems, which have RoW << R/W, is therefore strongly influenced by Coriolis forces. Finally, these findings are incorporated into a conceptual model describing the evolution of submarine channels at different latitudes. For instance, the Northern Hemisphere channels have a distinctly higher right levee system and migrate predominantly to the left side and generally exhibit a low sinuosity. In contrast, low latitude channel systems have RoW >> R/W so that centrifugal forces are more dominant. This results in more sinuous submarine channel systems with varying levee asymmetries in subsequent channel bends. In conclusion, Coriolis forces are negligible around the equator but should be considered in high latitude systems, particularly when RoW ~ O(1) and RoW << R/W.

This thesis reports an experimental investigation of low Reynolds number particle-laden turbulent flows in a horizontal plane channel. Experiments were conducted over a smooth wall and over two rough surfaces made from sand grain and gravel of relative roughness k/h ≈ 0.08 and 0.25, respectively, where k is the roughness height and h is the channel half-height. The flow was loaded with small solid particles with diameters less than 1/10 of the length scale of the energy-containing eddies, and whose concentrations decreased with time due to sedimentation. A novel particle image velocimetry (PIV) method that employed colour filtering for phase discrimination was used to measure the velocities of the fluid and solid particles. Over the smooth wall, the particles mean velocity, turbulence intensities and Reynolds shear stress matched those of the unladen flow very well. There were substantial differences between particle and fluid profiles over the rough wall, which include more rapid reduction in the particle mean velocity and significantly larger turbulence intensities and Reynolds shear stress compared to the unladen flow values. Stratification of the particle concentration led to attenuation of the fluid wall-normal turbulence intensity. This effect was nullified by the roughness perturbation leading to collapse of the wall-normal turbulence intensities over the rough wall. The streamwise turbulence intensity also collapsed over the rough wall but it was found that particles augmented the fluid Reynolds shear stress due to enhanced correlation between the rough wall streamwise and wall-normal velocity fluctuations. A quadrant decomposition of the fluid Reynolds shear stress also revealed corresponding enhancements in ejections and sweeps, the dominant contributors to the Reynolds shear stress, over the rough wall. Based on two-point correlations between the velocity fluctuations and between the velocity fluctuations and swirling strength, it was concluded that both wall roughness and particles modified the turbulence structure by increasing the size of the larger-scale structures. The idea of eddies growing from the wall, thereby enhancing communication between the inner layer and outer parts of the flow, has implications for wall-layer models that assume that the outer layer is detached from the turbulence in the inner region.

The spatial and temporal structure of the flow, vorticity and stress over equilibrium orbital-scale sand ripples are investigated at turbulence-resolving scales with a wide-band coherent Doppler profiler (MFDop) in an oscillating tray apparatus. The oscillation period and horizontal excursion were 10 s and 0.5 m. Velocity profiles were acquired with 3 mm vertical resolution and at a 42 Hz sampling rate. Ripple wavelength and amplitude were 25 cm and 2.2 cm. The MFDop measurements are used to investigate the development of the lee vortex as a function of phase, and the co-evolution of turbulent kinetic energy, Reynolds stress and turbulence production. Shear stress is determined from the vertically-integrated vorticity equation and using the double-averaging approach. Friction factors obtained from the two methods are comparable and range from 0.1 to 0.2.

Large-eddy simulations are carried out to study the combined effects of roughness and favourable pressure gradient in boundary layer flows, where the high acceleration (on smooth walls) may cause flow reversion to the quasi-laminar state. A sand-grain roughness model is used, with the no-slip boundary condition modeled by an immersed boundary method. The properties and accuracies of the scheme are studied, the roughness model is validated, and the spatial-resolution requirements are determined. The roughness model is applied to boundary layers subject to mild or strong acceleration, with simulations carried out underlining the effects of three parameters: the acceleration parameter, the roughness height, and the inlet Reynolds number. The roughness effects are limited to the roughness sublayer; the outer layer is affected indirectly only, through the changes that roughness causes in the relaminarization and retransition processes. The roughness significantly affects the inner-layer quantities like the friction velocity and the friction coefficient, while the local Reynolds number, the outer-layer mean velocity, as well as the Reynolds stresses beyond the roughness sublayer, are not sensitive to the roughness. The acceleration decreases the Reynolds stresses in the overlap region and promotes a laminar-like velocity profile. The acceleration leads to stabilization of near-wall structures and causes one-dimensional turbulence. The roughness generates small-scale structures at the bottom wall, which disturb the larger structures originally stabilized by the pressure gradient, leading to a decrease in the Reynolds-stress anisotropy. Roughness increases the Reynolds stresses in the roughness sublayer and tends to restore the fully turbulence flow early. The inlet Reynolds number affects the flow stability by determining the viscous length scale compared to the roughness length scales, and by determining how far the roughness effect extents into the boundary layer. / Thesis (Master, Mechanical and Materials Engineering) -- Queen's University, 2011-10-17 11:19:08.063

Thesis (Ph.D.); Submitted to the Department of Physics and Astronomy, Rice University, Houston, TX (US); 1 Sep 2004. / Published through the Information Bridge: DOE Scientific and Technical Information. "LA-14149-T" J. Koller. US DOE (US) 09/01/2004. Report is also available in paper and microfiche from NTIS.