Separating and reattaching flows are important in a large number of engineering configurations. The flow through a sudden expansion (backward-facing step) represents a conceptually simple case of this class of flows and hence has been the subject of numerous studies. The present study focuses on the effect of the expansion ratio (defined as the ratio of downstream channel height to upstream channel height) on the unsteady flow features in the reattachment region and further downstream. It is known that this flow demonstrates two different instabilities; the Kelvin-Helmholtz shear layer instability, which scales with the shear layer thickness, and the instability associated with the separation bubble, which scales with the step height and has similarities to K´arman vortex shedding behind a cylinder.In addition to these, there is also a possibility of the presence of the ‘preferred’ mode of the jet issuing from the inlet channel of the sudden expansion, especially at high expansion ratios, where the flow resembles a wall jet. The aim of this study is to investigate experimentally the changes in the instability of the separation bubble, as the expansion ratio is changed, and its possible interactions with the other instabilities in the flow.One might expect some changes in the flow with expansion ratio, as at low expansion ratios, the configuration represents a simple backward-facing step geometry, while at high expansion ratios, the geometry approaches that of a wall jet. A variable expansion ratio backward-facing step facility has been developed in an open circuit wind tunnel.This facility permits continuous variation of the expansion ratio from 1 to around 6. Attention is focused on the turbulent regime of the flow, where the flow structure has been found in previous studies to be relatively insensitive to the Reynolds number. The inlet conditions have been kept constant with a thin turbulent boundary layer at the step, the boundary layer thickness at separation being approximately 14 % of the inlet channel height. The Reynolds number based on the inlet channel height, H, is kept constant at Re=48,000 and the expansion ratio is varied by changing the channel height downstream of the step. Detailed hot wire measurements have been made to characterize the spatial variation of the dominant frequencies in the flow at different expansion ratios. The expansion ratio has been varied from a low value of 1.14 to a high value of 3.25, and detailed measurements are obtained for five expansion ratios of 1.14, 1.3, 1.5, 2.0 and 3.0. Further, to elucidate the dominant vortical structures in the flow, Particle Image Velocimetry measurements have been undertaken simultaneously with hot wire measurements for the case of expansion ratio 1.5, which have permitted the conditional averaging of vorticity fields.These investigations have brought forth some interesting features of the flow over a backward-facing step. Results for the time-mean properties of the flow indicate that the shear layer separating from the step deviates from a free mixing layer behaviour away from the step, possibly due to its interaction with the wall and the recirculation region underneath it. At any given streamwise location, the shear layer momentum thickness, θ, is seen to increase with the expansion ratio. Further, upto reattachment, the momentum thickness of the shear layer is seen to scale with the step height, h, independent of its initial thickness at separation, θo, as long as the boundary layer at separation is sufficiently thin as compared to the step height. Investigations for the unsteady flow features show that the frequency of the dominant peak in the velocity spectrum, supposed to represent the passage frequency (Strouhal number, S, based on the step height, h, and the inlet velocity, U) of the vortical structures, varies in the cross stream (y) direction, in addition to its expected variation in the streamwise (x) direction. The variation of the Strouhal number in the cross stream direction is seen to scale with the local momentum thickness of the shear layer, except for locations very close to the step. To characterize the development of the dominant frequency in the streamwise direction, the maximum value of the Strouhal number at a streamwise location is taken to be the representative value for that streamwise location. The Strouhal number is seen to decrease in the streamwise direction, from a very high value near the step, to a value of approximately 0.08 in the reattachment region, and remains constant further downstream. This value, supposed to represent the large scale structures shed from the reattachment region, is seen to remain very close to 0.08 for all Expansion ratios investigated. Conditional averaging of the vorticity fields in the reattachment region is done for an expansion ratio of 1.5, to get a detailed picture of the unsteady flow field. The hot wire signal at the outer edge of the shear layer in the reattachment region, which represents the non-dimensional structure passage frequency of S=0.08, is used as the conditioning signal. Results seem to indicate that the recirculation region, or the ‘bubble’ divides into two cells, and sheds the downstream cell quasi-periodically. The passage of these structures through the reattachment region seems to be concomitant With a local vertical motion of the shear layer. Further, the streamwise development of the local Strouhal number, Sθ, based on the local momentum thickness of the shear layer, and the local free stream velocity, Umax, indicates a possibility of a coupling between the shear layer and the structures shed from the reattachment region.

The goal of this thesis is the development of a fast and robust algorithm that is able to detect patterns in flow fields independent from their orientation and adequately visualize the results for a human user. This thesis is an interdisciplinary work in the field of vector field visualization and the field of pattern recognition. A vector field can be best imagined as an area or a volume containing a lot of arrows. The direction of the arrow describes the direction of a flow or force at the point where it starts and the length its velocity or strength. This builds a bridge to vector field visualization, because drawing these arrows is one of the fundamental techniques to illustrate a vector field. The main challenge of vector field visualization is to decide which of them should be drawn. If you do not draw enough arrows, you may miss the feature you are interested in. If you draw too many arrows, your image will be black all over. We assume that the user is interested in a certain feature of the vector field: a certain pattern. To prevent clutter and occlusion of the interesting parts, we first look for this pattern and then apply a visualization that emphasizes its occurrences. In general, the user wants to find all instances of the interesting pattern, no matter if they are smaller or bigger, weaker or stronger or oriented in some other direction than his reference input pattern. But looking for all these transformed versions would take far too long. That is why, we look for an algorithm that detects the occurrences of the pattern independent from these transformations. In the second part of this thesis, we work with moment invariants. Moments are the projections of a function to a function space basis. In order to compare the functions, it is sufficient to compare their moments. Normalization is the act of transforming a function into a predefined standard position. Moment invariants are characteristic numbers like fingerprints that are constructed from moments and do not change under certain transformations. They can be produced by normalization, because if all the functions are in one standard position, their prior position has no influence on their normalized moments. With this technique, we were able to solve the pattern detection task for 2D and 3D flow fields by mathematically proving the invariance of the moments with respect to translation, rotation, and scaling. In practical applications, this invariance is disturbed by the discretization. We applied our method to several analytic and real world data sets and showed that it works on discrete fields in a robust way.

Recent advances in nanotechnology have led to the production of materials with nanoscale dimensions (nm) and properties distinctly different from their bulk (>100 nm) counterparts. With increased use, it is inevitable that nanomaterials will accumulate in the environment and there is concern that the novel properties of nanomaterials could result in detrimental environmental and human health effects. In particular, there has been concern recently regarding the use of silver (Ag) based nanomaterials as antimicrobial agents in consumer and medical products. Current regulations dealing with the discharge of metals into the environment are based on total concentrations with no consideration for the form (e.g., ionic, nanoparticle, colloid) which can largely determine toxicity. Methods for the identification and characterization of nanoparticulates within complex matrices are lacking and the development of robust methods for this purpose are considered a high priority research area. This research focuses on the development and application of a novel method for characterizing Ag manufactured nanoparticles (MNPs) within terrestrial environments, in particular in soil pore water, with applications relevant to other metal MNPs as well. The method was then applied to understand the dynamics and behavior of Ag MNPs in soil and soil amended with sewage sludge biosolids.

Increasing environmental concerns as well as diminishing fossil fuel reserves call for a new generation of energy conversion technologies. Fuel cells, which convert the chemical energy of a fuel directly to electrical energy, have been identified as one of the leading alternative energy conversion technologies. Fuel cells are more efficient than conventional heat engines with minimal pollutant emissions and superior scalability. Proton Exchange Membrane Fuel Cells (PEMFCs) which produce electricity from hydrogen have been widely investigated for transportation and stationary applications. The focus of this study is on the Direct Propane Fuel Cell (DPFC), which belongs to the PEMFC family, but consumes propane instead of hydrogen as feedstock. A drawback associated with DPFCs is that the propane reaction rate is much slower than that of hydrogen. Two ideas were suggested to overcome this issue: (i) operating at high temperatures (150-230oC), and (ii) keeping the propane partial pressure at the maximum possible value. An electrolyte material composed of zirconium phosphate (ZrP) and polytetrafluoroethylene (PTFE) was suggested because it is an acceptable proton conductor at high temperatures. In order to keep the propane partial pressure at the maximum value, interdigitated flow-fields were chosen to distribute propane through the anode catalyst layer. In order to evaluate the performance of a DPFC which operates at high temperature and uses interdigitated flow-fields, a computational approach was chosen. Computational Fluid Dynamics (CFD) was used to create two 2-D mathematical models for DPFCs based on differential conservation equations. Two different approaches were investigated to model species transport in the electrolyte phase of the anode and cathode catalyst layers and the membrane layer. In the first approach, the migration phenomenon was assumed to be the only mechanism of proton transport. However, both migration and diffusion phenomena were considered as mechanisms of species transport in the second approach. Therefore, Ohm's law was used in the first approach and concentrated solution theory (Generalized Stefan-Maxwell equations) was used for the second one. Both models are isothermal. The models were solved numerically by implementing the partial differential equations and the boundary conditions in FreeFEM++ software which is based on Finite Element Methods. Programming in the C++ language was performed and the existing library of C++ classes and tools in FreeFEM++ were used. The final model contained 60 pages of original code, written specifically for this thesis. The models were used to predict the performance of a DPFC with different operating conditions and equipment design parameters. The results showed that using a specific combination of interdigitated flow-fields, ZrP-PTFE electrolyte having a proton conductivity of 0.05 S/cm, and operating at 230oC and 1 atm produced a performance (polarization curve) that was (a) far superior to anything in the DPFC published literature, and (b) competitive with the performance of direct methanol fuel cells. In addition, it was equivalent to that of hydrogen fuel cells at low current densities (30 mA/cm2).

This thesis describes the theoretical bases, development and testing of an integrated moving object detection framework in 2D and 3D scenes. The detection problem is analyzed in stationary and non-stationary camera sequences and different algorithms are developed for each case. Two methods are proposed in stationary camera sequences: background extraction followed by differencing and thresholding, and motion detection using optical flow field calculated by &ldquo / Kanade-Lucas Feature Tracker&rdquo / . For non-stationary camera sequences, different algorithms are developed based on the scene structure and camera motion characteristics. In planar scenes where the scene is flat or distant from the camera and/or when camera makes rotations only, a method is proposed that uses 2D parametric registration based on affine parameters of the dominant plane for independently moving object detection. A modified version of the 2D parametric registration approach is used when the scene is not planar but consists of a few number of planes at different depths, and camera makes translational motion. Optical flow field segmentation and sequential registration are the key points for this case. For 3D scenes, where the depth variation within the scene is high, a parallax rigidity based approach is developed for moving object detection. All these algorithms are integrated to form a unified independently moving object detector that works in stationary and non-stationary camera sequences and with different scene and camera motion structures. Optical flow field estimation and segmentation is used for this purpose.

A circular cylinder is a geometrically simple bluff body that occurs in various practical applications. As with any bluff body, it exhibits large drag forces and a strong fluctuating lift force, both related to the strong shedding of vortices from the body, which is commonly referred to as the Karman Street. Rigid splitter plates in the wake of the cylinder are known to suppress shedding from the body, and thereby result in reduced drag and fluctuating lift forces, the latter being important to reduce flow-induced vibrations of the body. In the present work, the flow past a cylinder with a downstream flexible splitter plate/flap is studied, the length (L) and flexural rigidity (EI) of the flap being the main parameters besides the flow speed (U). Two flaps length to cylinder diameter ratios (L/D), namely, a short (L/D = 2) and a long (L/D = 5) flaps have been studied, the shorter one being smaller than the recirculation zone, while the larger is longer than the recirculation zone. In both these cases, the flexural rigidity (EI) and the flow speed are systematically varied. In all cases, the flaps motion are directly visualized, the lift and drag forces are measured with a force balance, and the wake velocity field is measured using PIV. In both the long and short flaps cases, the flexural rigidity (EI) of the flexible flap has been varied over a large range of values, and it has been found that the results for flaps tip motion and forces collapse well when plotted with a non-dimensional bending stiffness (K∗), which is defined as K∗ = EI/(1/2ρU2L3). This collapse occurs across flexible flaps with different values of EI, as long as Re > 5000. The collapse is not found to be good for Re < 5000. This difference appears to be related to the large reduction in fluctuating lift for a bare cylinder in the Re range between approximately 1600 and 5000 discussed by Norberg[41]. In the long flap case, the existence of two types of periodic modes is found within the range of K∗ values from 5 × 10−6 to 1 × 10−1 studied. The first one corresponds to a local peak in amplitude at K∗ ≈ 1.5 × 10−3 that is referred to as mode I, and the second that occurs at low values of K∗ (K∗ < 3 × 10−5) that is referred to as mode II. The fluctuating lift is found to be minimum for the mode I oscillation. The mean drag is also found to reach a broad minimum that starts at K∗ corresponding to mode I and continues to be at the same low level of approximately 65% of the bare cylinder drag for all higher K∗ values, representing an approximately 35% decrease in mean drag of the cylinder. The wake measurements also show significant changes with K∗. The formation length (lf /D) obtained from the closure point of the mean separation bubble is found to continuously increase with K∗, reaching values of approximately 2.6 at mode I and thereafter only small increases are seen as K∗ is increased to large values corresponding to the rigid splitter plate case, consistent with the observed variations in the mean drag. The stream wise and cross-stream turbulent intensities and the Reynolds shear stress are all found to be strikingly lower in the mode I case compared to the bare cylinder case, and more importantly, these values are even lower than the rigid splitter plate case. This is consistent with the shedding of weaker vortices and with the minimum in fluctuating lift found in the mode I case. The results for this flap length show that the mode I flap oscillation, corresponding to K∗ ≈ 1.5 × 10−3, may be useful to reduce lift, drag, velocity fluctuations in the wake and the strength of the shed vortices. In particular, the wake fluctuations corresponding to this mode are found to be significantly lower than the rigid splitter plate case. In the short flap case (L/D = 2), it is found that there exists a richer set of flapping modes compared to the long flap, with these modes being dependent on K∗. At low K∗ values, the flap exhibits large amplitude symmetric flap motion that is referred to as mode A, while clearly asymmetric flaps motion are seen at higher K∗ values corresponding to modes B and C. Mode B corresponds to asymmetric large amplitude flapping motion, while mode C is also asymmetric with the flap clearly deflected off to one side, but having small oscillation amplitudes. At even higher K∗ values, corresponding to mode D, symmetric flaps motion are again seen with the amplitudes being smaller than in mode A. Apart from the flap tip amplitude, the non-dimensional frequency of flap tip motion also changes as the flap changes modes. In this case, there is a minimum in the fluctuating lift corresponding to mode B and C oscillation. The mean drag is found to reach a minimum again corresponding to mode C, which corresponds to an approximately 35% decrease in mean drag of the cylinder. In this case, there is a large increase in fluctuating lift (approximately 150% of the bare cylinder case) at higher values of K∗ that appears to correspond to a “resonant” condition between the structural natural frequency of the flexible splitter plate/flap and the wake shedding frequency of the bare cylinder. The wake measurements show that the formation length (lf /D) is the largest for mode C (deflected flap state), which is consistent with the observed minimum in mean drag observed for this mode. The stream wise and cross-stream turbulent intensities and the Reynolds shear stress are all found to be strikingly lower in the mode C case compared to the bare cylinder case, with the values for the Reynolds shear stress being lower than the rigid splitter plate case. This is again consistent with the minimum in fluctuating lift found in the mode C case. The results for this flap length show that the mode C flap oscillation, corresponding to K∗ ≈ 5 × 10−2 that correspond to a deflected flap state with very small oscillation may be useful to reduce lift, drag, velocity fluctuations in the wake and the strength of the shed vortices. The results from the present study show that the flexible flap/splitter plate down-stream of the cylinder exhibits a variety of mode shapes depending on the effective bending rigidity of the flap K∗ for both the long and short flaps cases. The forces and the wake are also found to be strongly dependent on this parameter K∗ with the wake fluctuations, lift fluctuations and the drag being very effectively suppressed at an intermediate value of K∗ that is found to be dependent on the plate/flap length.

Increasing environmental concerns as well as diminishing fossil fuel reserves call for a new generation of energy conversion technologies. Fuel cells, which convert the chemical energy of a fuel directly to electrical energy, have been identified as one of the leading alternative energy conversion technologies. Fuel cells are more efficient than conventional heat engines with minimal pollutant emissions and superior scalability. Proton Exchange Membrane Fuel Cells (PEMFCs) which produce electricity from hydrogen have been widely investigated for transportation and stationary applications. The focus of this study is on the Direct Propane Fuel Cell (DPFC), which belongs to the PEMFC family, but consumes propane instead of hydrogen as feedstock. A drawback associated with DPFCs is that the propane reaction rate is much slower than that of hydrogen. Two ideas were suggested to overcome this issue: (i) operating at high temperatures (150-230oC), and (ii) keeping the propane partial pressure at the maximum possible value. An electrolyte material composed of zirconium phosphate (ZrP) and polytetrafluoroethylene (PTFE) was suggested because it is an acceptable proton conductor at high temperatures. In order to keep the propane partial pressure at the maximum value, interdigitated flow-fields were chosen to distribute propane through the anode catalyst layer. In order to evaluate the performance of a DPFC which operates at high temperature and uses interdigitated flow-fields, a computational approach was chosen. Computational Fluid Dynamics (CFD) was used to create two 2-D mathematical models for DPFCs based on differential conservation equations. Two different approaches were investigated to model species transport in the electrolyte phase of the anode and cathode catalyst layers and the membrane layer. In the first approach, the migration phenomenon was assumed to be the only mechanism of proton transport. However, both migration and diffusion phenomena were considered as mechanisms of species transport in the second approach. Therefore, Ohm's law was used in the first approach and concentrated solution theory (Generalized Stefan-Maxwell equations) was used for the second one. Both models are isothermal. The models were solved numerically by implementing the partial differential equations and the boundary conditions in FreeFEM++ software which is based on Finite Element Methods. Programming in the C++ language was performed and the existing library of C++ classes and tools in FreeFEM++ were used. The final model contained 60 pages of original code, written specifically for this thesis. The models were used to predict the performance of a DPFC with different operating conditions and equipment design parameters. The results showed that using a specific combination of interdigitated flow-fields, ZrP-PTFE electrolyte having a proton conductivity of 0.05 S/cm, and operating at 230oC and 1 atm produced a performance (polarization curve) that was (a) far superior to anything in the DPFC published literature, and (b) competitive with the performance of direct methanol fuel cells. In addition, it was equivalent to that of hydrogen fuel cells at low current densities (30 mA/cm2).

The sorption of humic acid (HA) and HA size fractions onto magnetite has been studied. There is considerable irreversibility in the interaction of the humic with the magnetite surface, but the presence of Eu3+ ions has no effect on the sorption of humic onto magnetite. The magnitude of the sorption to magnetite increases with HA fraction size for all ionic strengths between 0.01 and 3 mol dm-3. Increasing ionic strength also increases sorption. Asymmetric Flow Field Flow Fractionation analysis of HA sorption to magnetite after 1 day revealed preferential sorption of lower molecular weight material. Eu3+ sorption onto magnetite was studied as a function of Eu concentration, which showed an increase in relative sorption as Eu concentration decreased. The behaviour of Eu3+ in ternary (HA/Eu3+/magnetite) systems is heavily influenced by HA, and from the data there is direct evidence for ternary complex formation. Larger HA size fractions retain more Eu3+ in solution than the smaller fractions. The binding strengths of HA size fractions were determined through ion exchange resin experiments: generally the larger fractions (> 10 kDa) showed stronger binding than the smaller components, but the unfractionated sample showed the strongest binding.First order dissociation rate constants have been determined for the whole HA and HA size fractions. The dissociation rate constants are independent of HA fraction size, but the larger species bind more Eu non-exchangeably. Time series ultrafiltration of Eu3+/whole humic mixtures has shown a shift in the distribution of metal ions to larger size fractions after a few days. Two ternary system kinetic speciation models have been developed to predict the behaviour of HA and Eu3+ in ternary systems. The two differ in their description of the multi-component behaviour of the binary HA-mineral interaction. The first assumes a single HA species and two surface binding sites and was found to perform better overall than the second, which has a single surface sorption site and two HA species in solution. The exchangeable binding strengths for the different HA samples calculated from both models showed similarities to those measured experimentally.

In this contribution we use wavefront shaping approaches for image correlation based flow-field measurements for the first time. Aberrations introduced by a single phase boundary in the detection beam path were explored. Variations of the optical path-length result in strong errors in position allocation and thus to an enhancement of the measurement uncertainty of the velocity. Our results show that the usage of wavefront shaping enables to reduce these errors and to strongly improve the quality of image correlation based flow-field measurements. First experimental and simulated results underline the importance of these approaches.

We propose the usage of wavefront shaping approaches for image correlation based flow-field measurements. Aberrations introduced by a single phase boundary in the detection beam path were explored. Variations of the optical path-length result in strong errors in position allocation and thus to an enhancement of the measurement uncertainty of the velocity. Our results show that the usage of wavefront shaping enables to reduce these errors and to strongly improve the quality of image correlation based flow-field measurements. First experimental and simulated results underline the importance of these approaches