ANALYSIS OF HYPORHEIC FLOW INDUCED BY A BAR IN A GRAVEL STREAM: AN EXPERIMENTAL STUDY

Pipelines are one of the safest means to transport crude oil, but are not spill-free. This is of concern in North America, due to the large volumes of crude oil shipped by Canadian producers and the lengthy network of pipelines. Each pipeline crosses many rivers, supporting a wide variety of human activities, and rich aquatic life. However, there is a knowledge gap on the risks of contamination of river beds due to oil spills. This thesis addresses this knowledge gap by focussing on mechanisms that transport water (and contaminants) from the free surface flow to the bed sediments, and vice-versa. The work focuses on gravel rivers, in which bed sediments are sufficiently permeable that pressure gradients caused by the interactions of flow with topographic elements (gravel bars), or changes in direction induce exchanges of water between the free surface flow and the bed, known as hyporheic flows. The objectives of the thesis are: to present a new method to visualize and quantify hyporheic flows in laboratory experiments; to conduct a novel series of experiments on hyporheic flow induced by a gravel bar under different free surface flows. The new method to quantify hyporheic flows rests on injections of a solution of dye and water. The method yielded accurate flow lines, and reasonable estimates of the hyporheic flow velocities. The present series of experiments was carried out in a 11 m long, 0.39 m wide, and 0.41 m deep tilting flume. The gravel had a mean particle size of 7.7 mm. Different free surface flows were imposed by changing the flume slope and flow depth. Measured hyporheic flows were turbulent. Smaller free surface flow depths resulted in stronger hyporheic flows (higher velocities, and deeper dye penetration into the sediment). A significant finding is that different free surface flows (different velocities, Reynolds number, etc.) produce similar hyporheic flows as long as the downstream hydraulic gradients are similar. This suggests, that for a specified bar geometry, the characteristics of the hyporheic flows depend on the downstream hydraulic gradients, and not or only minimally on the internal dynamics of the free surface flow.

Laboratory Investigation of Diluted Bitumen Trapping and Dissolution in Gravel

The Canadian economy is largely dependent on the distribution of large volumes of oil to domestic and international markets by a long network of pipelines. Unfortunately, accidents occur, and oil can leak or spill from these pipelines before it reaches its destination. Of particular concern are the long-term consequences of oil spills in freshwater, which include sinking of oil in water and the contamination of sensitive areas, such as where fish (e.g., salmon) deposit their eggs in gravel-dominated river sediments. There is a knowledge gap regarding the fate and behaviour of oil in river sediment. To this end, this study aimed at finding the potential for diluted bitumen (dilbit) oil to become trapped in gravel and to transfer hydrocarbons into water by dissolution, which are harmful to aquatic life. Two sets of laboratory experiments were conducted to simulate conditions of an oil spill on an exposed bank or in shallow water. In the first set, by conducting capillary pressure-saturation (Pc-Sw) experiments it was found that dilbit can enter gravel pores without much resistance and approximately 14% of the pore volume can be occupied by discontinuous single or multipore blobs of dilbit following imbibition of water. Air-water Pc-Sw experiments done in laboratory 1-D columns required gravity correction and could be successfully scaled to predict dilbit-water Pc-Sw curves, except for the trapped amount of dilbit. Trapped dilbit constituents can be dissolved into the water flowing through gravel pores (hyporheic flow) at different velocities. In the second set, dissolution experiments suggested that for the duration of the test, hydrocarbons that cause acute toxicity dissolve rapidly, likely resulting in a decrease in their effective solubility. However, dilbit saturation changed only <2% within that time. Chronically toxic PAH compounds were also detected in the effluent water. The total concentration of all detected PAHs and alkylPAHs exceeded the threshold literature value to protect early-life stage fish. Observations of decreased concentrations with increased aqueous velocities as well as less than equilibrium concentrations indicated that the mass transfer was rate-limited. A correlation was developed for the mass transfer rate coefficient to understand the mass transfer behaviour beyond the conditions used in the experiments, which had a Reynolds number exponent similar to the studies of NAPL dissolution in groundwater.

Experimental and RANS-CFD Study of Nacelle Drag

The drag on a nacelle model was investigated experimentally and computationally to provide guidance and insight into the capabilities of RANS-based CFD. The research goal was to determine whether industry constrained CFD could participate in the aerodynamic design of nacelle bodies. Grid refinement level, turbulence model and near wall treatment settings, to predict drag to the highest accuracy, were key deliverables. Cold flow low-speed wind tunnel experiments were conducted at a Reynolds number of 6∙〖10〗^5, 293 K and a Mach number of 0.1. Total drag force was measured by a six-component force balance. Detailed wake analysis, using a seven-hole pressure probe traverse, allowed for drag decomposition via the far-field method. Drag decomposition was performed through a range of angles of attack between 0o and 45o. Both methods agreed on total drag within their respective uncertainties. Reversed flow at the measurement plane and saturation of the load cell caused discrepancies at high angles of attack. A parallel CFD study was conducted using commercial software, ICEM 15.0 and FLUENT 15.0. Simulating a similar nacelle geometry operating under inlet boundary conditions obtained through wind tunnel characterization allowed for direct comparisons with experiment. It was determined that the Realizable k-ϵ was best suited for drag prediction of this geometry. This model predicted the axial momentum loss and secondary flow in the wake, as well as the integrated surface forces, within experimental error up to 20o angle of attack. SST k-ω required additional surface grid resolution on the nacelle suction side, resulting in 15% more elements, due to separation point prediction sensitivity. It was further recommended to apply enhanced wall treatment to more accurately capture the viscous drag and separated flow structures. Overall, total drag was predicted within 5% at 0o angle of attack and 10% at 20o, each within experimental uncertainty. What is more, the form and induced drag predicted by CFD and measured by the wake traverse shared good agreement. Which indicated CFD captured the key flow features accurately despite simplification of the nacelle interior geometry.