Junction flows originate from the interaction between a fluid moving over a wall with an obstacle mounted on the same surface. Understanding the physics of such flows is of great interest to engineers responsible for the design of systems consisting of wall-body junctions. From aerodynamics to turbomachinery and electronics to bridge hydraulics, a number of phenomena (drag, heat transfer, scouring) are driven by the behavior of the most prominent feature of junction flows: the horseshoe vortex system (HVS). Focusing on turbulent flows, the complex dynamics of the HVS is established through its unsteadiness and non-uniformity. The fundamentals of this dynamically-rich phenomenon have been described within the body of a rapidly-expanding literature. Nevertheless, important aspects remain inadequately understood and call for further scrutiny. This study emphasized three of them, by investigating the effects of: model scale, wall roughness, and bed geometry. High-resolution experiments were carried out using Particle Image Velocimetry (PIV). Statistical analyses, vortex identification schemes, and Proper Orthogonal decomposition were employed to extract additional information from the large PIV datasets. The time-averaged topology of junction flows developing over a smooth and impermeable wall was independent of the flow Reynolds number, Re (parameter that expresses the effects of scale). On the contrary, time-resolved analysis revealed a trend of increasing vorticity, momentum, and eruptions of near-wall fluid with Re. New insights on the modal dynamics of the HVS were also documented in a modified flow mechanism. Wall roughness (modeled with a permeable layer of crushed stones) diffused turbulence and vorticity throughout the domain. This effect manifested with high levels of intermittency and spatial irregularity for the HVS. Energetic flow structures were also identified away from the typical footprint of the HVS. Finally, a novel implementation of PIV allowed for unique velocity measurements over an erodible bed. It was demonstrated that, during the initial stages of scouring, the downflow at the face of the obstacle becomes the dominant flow characteristic in the absence of the HVS. Notwithstanding modeling limitations, the physical insight contributed here could be used to enhance the design of systems with similar flow and geometrical characteristics. / Ph. D.

The development of aerodynamic boundary layers on wind turbine blades is an important consideration in their performance. It can be quite challenging to replicate full scale conditions in laboratory experiments, and advanced diagnostics become valuable in providing data not available from traditional means. A new variant of Doppler global velocimetry (DGV) known as cross-correlation DGV is developed to measure boundary layer profiles on a wind turbine blade airfoil in the large scale Virginia Tech Stability Wind Tunnel. The instrument provides mean velocity vectors with reduced sensitivity to external conditions, a velocity measurement range from 0ms^-1 to over 3000ms^-1, and an absolute uncertainty. Monte Carlo simulations with synthetic signals reveal that the processing routine approaches the Cramér-Rao lower bound in optimized conditions. A custom probe-beam technique is implanted to eliminate laser flare for measuring boundary layer profiles on a DU96-W-180 wind turbine airfoil model. Agreement is seen with laser Doppler velocimetry data within the uncertainty estimated for the DGV profile. Lessons learned from the near-wall flow diagnostics development were applied to a novel benchmark model problem incorporating the relevant physical mechanisms of the high amplitude periodic turbulent flow experienced by turbine blades in the field. The model problem is developed for experimentally motivated computational model development. A circular cylinder generates a periodic turbulent wake, in which a NACA 63215b airfoil with a chord Reynolds number Re_c = 170, 000 is embedded for a reduced frequency k = (pi)fc/V = 1.53. Measurements are performed with particle image velocimetry on the airfoil suction side and in highly magnified planes within the boundary layer. Outside of the viscous region, the Reynolds stress profile is consistent with the prediction of Rapid Distortion Theory (RDT), confirming that the redistribution of normal stresses is an inviscid effect. The fluctuating component of the phase- averaged turbulent boundary layer profiles is described using the exact solution to laminar Stokes flow. A phase lag similar to that in laminar flow is observed with an additional constant phase layer in the buffer region. The phase lag is relevant for modeling the intermittent transition and separation expected at full scale. / Ph. D.

The entrainment process in two-phase buoyant jets differs significantly from their singlephase counterparts, and is not well understood. Entrainment models developed for singlephase flow are often used in two-phase jetting simulations, albeit with limited success. In this work, Particle Image Velocimetry (PIV) and shadowgraph flow visualization experiments have been conducted on submerged round gas jets of varying speeds and nozzle diameters with the goal of improving our understanding of the entrainment process in a two-phase (gas-liquid) jet. The total entrainment estimated using the PIV measurements is higher than the respective values suggested by a common empirical model developed for singlephase buoyant jets. A two-phase theoretical entrainment model used for comparison shows an overestimation of entrainment, but predicts the increase in the rate of entrainment with axial distance from the jet nozzle seen in the PIV results. This thesis also presents advances in PIV processing methodology that were developed concurrently with the entrainment research. The novel Spectral Phase Correlation (SPC) allows for particle displacement to be determined directly from phase information in the Fourier domain. Some of the potential benefits of the SPC explored here include (1) avoidance of errors introduced by spatial peak-finding routines; (2) use of a modal analysis that can be used to provide information such as correlation quality; and (3) introduction of a means of incorporating information from multiple image windows. At low image noise levels, the method performs as well as an advanced CC-based method. However, difficulties unwrapping the aliased phase information cause the SPC's performance to degrade at high noise levels. / Master of Science

Heart disease is one of the primary causes of morbidity and mortality for the adult population over the age of 65. Furthermore, ailments such as hypertension can affect as many as 50% of the adult population over the age of 45. If left untreated, these ailments eventually precipitate the onset of diastolic dysfunction and heart failure. Diastolic dysfunction is the alteration or impairment of performance in either the left or right ventricle of the heart. Although there has been a marked increase in study of this disease, there is still an apparent difficulty to diagnose patients. Flow visualization techniques have been commonly employed to study the development of these diseases as they relate to the filling process of the ventricles. One method, Echo Particle Image Velocimetry (Echo-PIV) is a relatively new method for cardiac flow chamber visualization, with the potential to provide physicians with a cost-effective and safe method for obtaining high temporal resolution recordings for extending knowledge on the filling processes in cardiac chamber flow. This work presents a new approach to extending the capabilities of Echo-PIV for more accurate measurement of cardiac flows for patients with poor quality recordings. Currently, much of the literature notes that temporal resolution and poor acoustic windows results in exclusion from study. These recordings are more representative of the contrast-enhancement studies used by physicians to better identify chamber walls. When applying standard PIV cross-correlation techniques, measurements tend to fail due to image noise and artifacts. By implementing a Moving Ensemble (MWE) with Product of Correlation (PoC) processing scheme, measurement accuracy, reliability, and robustness can be obtained for measurement in left ventricular filling assessment. / Master of Science

Aerodynamic measurements were taken over models of the Boeing 777 high fidelity isolated landing gear in the 6- by 6-foot Virginia Tech Stability Wind Tunnel (VT-SWT) at a free-stream Mach number of 0.16. Noise control devices (NCD) were developed at Virginia Tech [9] to reduce noise by shielding gear components, reducing wake interactions and by streamlining the flow around certain landing gear components. Aerodynamic measurements were performed to understand the flow over the landing gear and also changes in the flow between "Baseline" and "NCD" configurations (without and with Noise Control Devices respectively). Hot-film, Pitot-static measurements and flow visualization using tufts were performed over an isolated 26% scale-model high fidelity landing gear for the "Baseline" and "NCD" configurations. Contours of turbulence intensity, normalized wake velocity and normalized total pressure loss for both configurations are compared. The "Baseline" configuration was also compared with the NASA Ames study conducted by Horne et al [7]. Hot-film measurements are also compared to Microphone Phased Array results which were acquired at Virginia Tech by Ravetta [8]. A novel technique for processing hot-film measurements by breaking turbulence into octave bands as acoustic measurements is presented. Particle Image Velocimetry (PIV) measurements were taken at six different locations over the 13% scale-model landing gear with no door and at a truck angle of zero degrees. Results are compared to PIV measurements taken over the wheels of a four-wheel landing gear by Lazos [10-12]. PIV results such as average velocity contours and vectors, streamlines and instantaneous velocity contours and vectors are presented. Results presented from PIV and flow visualization are in good agreement with results from Lazos [10-12]. / Master of Science

The tissue microenvironment is a complex living system containing heterogeneous mechanical and biophysical cues. Cellular components are surrounded by extracellular matrix and interstitial fluid, while transport of nutrients and biochemical factors is achieved via the vasculature. Each constituent of the tissue microenvironment can play a significant role in its ability to function normally. Many diseases including cancer have been linked with dysfunction in the tissue microenvironment; therefore an improved understanding of interaction between components of this complex system is needed. In vitro platforms mimicking the tissue microenvironment appear to provide the most promising avenue for studies of cell-cell and cell-matrix interactions as well as elucidation of the mechanisms leading to disease phenomena such as tumor metastasis. However, successful recapitulation of all three primary components of the tissue microenvironment in three dimensions has remained challenging. In particular, matching mechanical cues and biochemical transport to in vivo conditions is difficult because of lack of quantitative characterization of the physical properties and parameters of such platforms. In this work, extensive characterization of collagen I hydrogels, popular for use as extracellular matrix mimics, was performed in order to enable tuning to specific in vivo conditions. Additionally, perfusion of blood in a 3D tissue microenvironment platform fabricated using collagen hydrogels was characterized to enable future advances in in vitro modeling of the in vivo microenvironment. Finally, the tissue microenvironment platform is modified to enable biochemical gradients within the hydrogel and used to examine directed migration (chemotaxis) of human breast cancer cells in response to gradients in growth factor combined with varied stiffness and pore diameter of the extracellular matrix. / Ph. D.

Time varying flow structures are involved in a large percentage of fluid flows although there is still much unknown regarding their behavior. With the development of high spatiotemporal resolution measurement systems it is becoming more feasible to measure these complex flow structures, which in turn will lead to a better understanding of their impact. One method that has been developed for studying these flow structures is finite time Lyapunov exponents (FTLEs). These exponents can reveal regions in the fluid, referred to as Lagragnian coherent structures (LCSs), where fluid elements diverge or attract. Better knowledge of how these time varying structures behave can greatly impact a wide range of applications, from aircraft design and performance, to an improved understanding of mixing and transport in the human body. This work provides the development of new methodologies for measuring and studying three-dimensional time varying structures. Provided herein is a method to improve replacement of erroneous measurements in particle image velocimetry data, which leads to increased accuracy in the data. Also, a method for directly measuring the finite time Lyapunov exponents from particle images is developed, as well as an experimental demonstration in a three-dimensional flow field. This method takes advantage of the information inherently contained in these images to improve accuracy and reduce computational requirements. Lastly, this work provides an in depth look at the flow field for developing wall jets across a wide range of Reynolds numbers investigating the mechanisms that contribute to their development. / Ph. D.

Chrysopelea paradisi, also known as the flying snake, possesses one of the most unique forms of aerial locomotion found in nature, using its entire body as a dynamic lifting surface without the use of wings or membranes. Unlike other airborne creatures, this species lacks appendages to aid in controlling its flight trajectory and producing lift. The snake exhibits exception gliding and maneuvering capabilities compared with other species of gliders despite this lack of appendages. While gliding, C. paradisi morphs its body by expanding its ribs, essentially doubling its width and utilizing its entire length as a reconfigurable wing. Its cross-sectional shape transforms into a thick, airfoil shape with a concave ventral surface, outwards protruding lips at the leading and trailing edges, a somewhat triangular dorsal surface with a round apex, and fore-aft symmetry. This study investigated the aerodynamic performance of this unique shape by simulating a single, static segment of the snake's body over a wide range of Reynolds numbers (3,000 to 15,000) and angles of attack (-10 to 60o) to simulate the full range of the snake's flight kinematics. This is the first study on an anatomically accurate snake model, and few aerodynamic studies have been performed in this low Reynolds number regime. Load cell measurements and time-resolved digital particle image velocimetry (TRDPIV) were performed on a 2D anatomically accurate model to determine the lift and drag coefficients, wake dynamics, and vortex shedding characteristics. This geometry produced a maximum lift coefficient of 1.9 and maximum lift to drag ratio of 2.7, and maintained increases in lift up to 35o. Overall, this geometry demonstrated robust aerodynamic behavior by maintain significant lift production and near maximum lift to drag ratios over a wide range of test parameters. These aerodynamic characteristics may enable the flying snake to glide at steep angles and over a wide range of angles of attack, often encountered in gliding trajectories. This geometry also produced larger maximum lift coefficients than many other bluff bodies and airfoils in this low Reynolds number regime. This thesis is organized as follows. The first section contains a broad introduction on gliding flight and C. paradisi's unique mode of gliding. The following section is a manuscript that will be submitted to a journal and contains the experimental analysis on the snake's cross-sectional shape. Several appendices attached to the end of this thesis contain additional analysis and work performed throughout the duration of this project and unique Matlab algorithms developed during this research. / Master of Science

The aim of this work was to understand the physical requirements as well as to develop methodology required to employ Time Resolved Digital Particle Image Velocimetry (TRDPIV) for measuring high speed, high magnification, near wall flow fields. Previous attempts to perform measurements such as this have been unsuccessful because of both limitations in equipment as well as proper methodology for processing of the data. This work addresses those issues and successfully demonstrates a test inside of a transonic turbine cascade as well as a high speed high magnification wall jet. From previous studies it was established that flow tracer delivery is not a trivial task in a high speed high back pressure environment. Any TRDPIV measurement requires uniform spatial seeding density, but time-resolved measurements require uniform temporal seeding density as well. To this end, a high pressure particle generator was developed. This advancement enhanced current capability beyond what was previously attainable. Unfortunately, this was not sufficient to resolve the issue of seeding all together, and an advanced data reconstruction methodology was developed to reconstruct areas of the flow field that where lost do to inhomogeneous seeding. This reconstruction methodology, based on Proper Orthogonal Decomposition (POD), has been shown to produce errors in corrected velocities below tradition spatial techniques alone. The combination of both particle generator and reconstruction methodology was instrumental for successfully acquiring TRDPIV measurements in a high speed high pressure environment such as a transonic wind tunnel facility. This work also investigates the development of a turbulent wall jet. This experiment helped in demonstrating the capability of taking high speed high magnification TRDPIV measurements. This experiment was very unique in that it is one of only a few experiments that studied the developing region of these jets. The Reynolds number ranged for this experiment from 150 – 10,000 which corresponded to velocities of 1 - 80 m/s. The results from this experiment showed good agreement with currently published time averaged data. Using scaling laws for fully developed jets a new scaling law was found for the developing region of the jet that could be applied to all Reynolds numbers in this study. A temporal investigation was also carried out using the temporal coefficients from POD. A vortex identification scheme was also applied to all of the Reynolds numbers showing clear trends as Reynolds number increased. / Master of Science

A stator cascade was developed to simulate the flow conditions within a close-stage-spacing transonic axial compressor. Experiments were conducted in a linear transonic blowdown cascade wind tunnel with an inlet Mach number of 0.65. The bow shock from the downstream rotor was simulated by a single moving normal shock generated with a shock tube. First, steady pressure data were gathered to ensure that the stator cascade operated properly without the presence of the shock. Next, the effects of the passing shock on the stator flow field were investigated using shadowgraph photography and Digital Particle Image Velocimetry (DPIV). Measurements were taken for three different shock strengths. In every case studied, a vortex formed near the stator trailing edge as the shock impacted the blade. The size of this vortex was shown to be directly related to the shock strength, and the vortex remained present in the trailing edge flow field throughout the cycle duration. Analysis of the DPIV data showed that the vortex acts as a flow blockage, with the extent of this blockage ranging from 2.9% of the passage for the weakest shock, to 14.3% of the passage for the strongest shock. The vortex was also shown to cause flow deviation up to 75° for the case with the strongest shock. Further analysis estimated that the total pressure losses due to shock-induced vorticity ranged from 46% to 113% of the steady wake losses. Finally, the total pressure loss purely due to the upstream-propagating normal shock was estimated to be roughly 0.22%. / Master of Science