This thesis focused on verifying previous literature and research that has been conducted on different spherical objects. Mainly, verifying literature that examines both how surface roughness contributes to the overall drag and how wake turbulence is affected by different surface roughness. The goal of this project is to be able to capture data that shows that the flow transition from laminar to turbulent occurs at lower Reynolds numbers for a rough spherical object rather than a perfectly smooth sphere. In order to achieve this goal, both force balance testing and hot-wire testing were conducted in the Aero-lab complex in USE170. The force balance was mounted and used in the larger wind tunnel while the hot-wire probe was mounted and used in the smaller wind tunnel. Both of the wind tunnels utilized LABVIEW software in order to collect and convert the qualitative values provided by the testing probes and equipment. The two main types of testing equipment that were used in this project were the force balance and the hot-wire probe. The overall results from the experiment were inconclusive based on the limitations of both the testing probes and the testing facility itself. Overall, the experiment yielded very limited results due to these limitations.

The liquid rocket engine, more specifically, the bi-propellant liquid rocket engine, is a popular type of chemical propulsion system within the propulsion industry due to its relatively high specific impulse and high thrust levels compared to the other chemical propulsion choices. For the purposes of this thesis, a bi-propellant liquid rocket engine system consists of a rocket engine, a set of tanks for the storage and supply of liquid propellants, and everything required in between for thrust-producing operation. Among the hardware in this "in between" necessary for a liquid rocket engine to produce thrust exists an injector, or an assembly of injector elements, whose purpose is to introduce and meter the flow of the fuel and oxidizer of the liquid rocket engine into the combustion chamber. To do this the injector or injector assembly, upon injection into the combustion chamber, must cause the two liquids to break up into small droplets, proportionally and uniformly distribute and mix the liquid into a spray pattern within the combustion chamber, and allow for engine combustion to occur as efficiently as possible. Daedalus Astronautics @ ASU, one of Arizona State University's engineering student organizations, has been working to design, construct, and successfully test a bi-propellant liquid rocket engine of its own. In doing so, Daedalus Astronautics has designed a bi-propellant liquid rocket engine injector assembly consisting of a forward bulkhead and an injector plate. The purpose of this thesis is to experimentally verify the flow of liquid through this injector assembly modeled using computational fluid dynamics methods. During the two semester time line allowed for this thesis project, a mesh was created for a single orifice geometry injector plate and combustion chamber assembly in ANSYS ICEM CFD and an experiment was designed for imaging the spray pattern from the injector plate and forward bulkhead assembly, from which several things about the injector geometry design were discovered.

We studied left ventricular flow patterns for a range of rotational orientations of a bileaflet mechanical heart valve (MHV) implanted in the mitral position of an elastic model of a beating left ventricle (LV). The valve was rotated through 3 angular positions (0, 45, and 90 degrees) about the LV long axis. Ultrasound scans of the elastic LV were obtained in four apical 2-dimensional (2D) imaging projections, each with 45 degrees of separation. Particle imaging velocimetry was performed during the diastolic period to quantify the in-plane velocity field obtained by computer tracking of diluted microbubbles in the acquired ultrasound projections. The resulting velocity field, vorticity, and shear stresses were statistically significantly altered by angular positioning of the mechanical valve, although the results did not show any specific trend with the valve angular position and were highly dependent on the orientation of the imaging plane with respect to the valve. We conclude that bileaflet MHV orientation influences hemodynamics of LV filling. However, determination of ‘optimal’ valve orientation cannot be made without measurement techniques that account for the highly 3-dimensional (3D) intraventricular flow.

The central purpose of this work is to investigate the large-scale, coherent structures that exist in turbulent Rayleigh-Bénard convection (RBC) when the domain is large enough for the classical ”wind of turbulence” to break down. The study exclusively focuses on the structures that from when the RBC geometry is a cylinder. A series of visualization studies, Fourier analysis and proper orthogonal decomposition are employed to qualitatively and quantitatively inspect the large-scale structures’ length and time scales, spatial organization, and dynamic properties. The data in this study is generated by direct numerical simulation to resolve all the scales of turbulence in a 6.3 aspect-ratio cylinder at a Rayleigh number of 9.6 × 107 and Prandtl number of 6.7. Single and double point statistics are compared against experiments and several resolution criteria are examined to verify that the simulation has enough spatial and temporal resolution to adequately represent the physical system.

Large-scale structures are found to organize as roll-cells aligned along the cell’s side walls, with rays of vorticity pointing toward the core of the cell. Two different large- scale organizations are observed and these patterns are well described spatially and energetically by azimuthal Fourier modes with frequencies of 2 and 3. These Fourier modes are shown to be dominant throughout the entire domain, and are found to be the primary source for radial inhomogeneity by inspection of the energy spectra. The precision with which the azimuthal Fourier modes describe these large-scale structures shows that these structures influence a large range of length scales. Conversely, the smaller scale structures are found to be more sensitive to radial position within the Fourier modes showing a strong dependence on physical length scales.

Dynamics in the large-scale structures are observed including a transition in the global pattern followed by a net rotation about the central axis. The transition takes place over 10 eddy-turnover times and the subsequent rotation occurs at a rate of approximately 1.1 degrees per eddy-turnover. These time-scales are of the same order of magnitude as those seen in lower aspect-ratio RBC for similar events and suggests a similarity in dynamic events across different aspect-ratios.

Cerebral aneurysms are pathological balloonings of blood vessels in the brain, commonly found in the arterial network at the base of the brain. Cerebral aneurysm rupture can lead to a dangerous medical condition, subarachnoid hemorrhage, that is associated with high rates of morbidity and mortality. Effective evaluation and management of cerebral aneurysms is therefore essential to public health. The goal of treating an aneurysm is to isolate the aneurysm from its surrounding circulation, thereby preventing further growth and rupture. Endovascular treatment for cerebral aneurysms has gained popularity over traditional surgical techniques due to its minimally invasive nature and shorter associated recovery time. The hemodynamic modifications that the treatment effects can promote thrombus formation within the aneurysm leading to eventual isolation. However, different treatment devices can effect very different hemodynamic outcomes in aneurysms with different geometries.

Currently, cerebral aneurysm risk evaluation and treatment planning in clinical practice is largely based on geometric features of the aneurysm including the dome size, dome-to-neck ratio, and parent vessel geometry. Hemodynamics, on the other hand, although known to be deeply involved in cerebral aneurysm initiation and progression, are considered to a lesser degree. Previous work in the field of biofluid mechanics has demonstrated that geometry is a driving factor behind aneurysmal hemodynamics.

The goal of this research is to develop a more combined geometric/hemodynamic basis for informing clinical decisions. Geometric main effects were analyzed to quantify contributions made by geometric factors that describe cerebral aneurysms (i.e., dome size, dome-to-neck ratio, and inflow angle) to clinically relevant hemodynamic responses (i.e., wall shear stress, root mean square velocity magnitude and cross-neck flow). Computational templates of idealized bifurcation and sidewall aneurysms were created to satisfy a two-level full factorial design, and examined using computational fluid dynamics. A subset of the computational bifurcation templates was also translated into physical models for experimental validation using particle image velocimetry. The effects of geometry on treatment were analyzed by virtually treating the aneurysm templates with endovascular devices. The statistical relationships between geometry, treatment, and flow that emerged have the potential to play a valuable role in clinical practice.

Image segmentation is of great importance and value in many applications. In computer vision, image segmentation is the tool and process of locating objects and boundaries within images. The segmentation result may provide more meaningful image data. Generally, there are two fundamental image segmentation algorithms: discontinuity and similarity. The idea behind discontinuity is locating the abrupt changes in intensity of images, as are often seen in edges or boundaries. Similarity subdivides an image into regions that fit the pre-defined criteria. The algorithm utilized in this thesis is the second category.

This study addresses the problem of particle image segmentation by measuring the similarity between a sampled region and an adjacent region, based on Bhattacharyya distance and an image feature extraction technique that uses distribution of local binary patterns and pattern contrasts. A boundary smoothing process is developed to improve the accuracy of the segmentation. The novel particle image segmentation algorithm is tested using four different cases of particle image velocimetry (PIV) images. The obtained experimental results of segmentations provide partitioning of the objects within 10 percent error rate. Ground-truth segmentation data, which are manually segmented image from each case, are used to calculate the error rate of the segmentations.

This dissertation describes a process for interface capturing via an arbitrary-order, nearly quadrature free, discontinuous Galerkin (DG) scheme for the conservative level set method (Olsson et al., 2005, 2008). The DG numerical method is utilized to solve both advection and reinitialization, and executed on a refined level set grid (Herrmann, 2008) for effective use of processing power. Computation is executed in parallel utilizing both CPU and GPU architectures to make the method feasible at high order. Finally, a sparse data structure is implemented to take full advantage of parallelism on the GPU, where performance relies on well-managed memory operations.

With solution variables projected into a kth order polynomial basis, a k+1 order convergence rate is found for both advection and reinitialization tests using the method of manufactured solutions. Other standard test cases, such as Zalesak's disk and deformation of columns and spheres in periodic vortices are also performed, showing several orders of magnitude improvement over traditional WENO level set methods. These tests also show the impact of reinitialization, which often increases shape and volume errors as a result of level set scalar trapping by normal vectors calculated from the local level set field.

Accelerating advection via GPU hardware is found to provide a 30x speedup factor comparing a 2.0GHz Intel Xeon E5-2620 CPU in serial vs. a Nvidia Tesla K20 GPU, with speedup factors increasing with polynomial degree until shared memory is filled. A similar algorithm is implemented for reinitialization, which relies on heavier use of shared and global memory and as a result fills them more quickly and produces smaller speedups of 18x.

The application of novel visualization and modeling methods to the study of cardiovascular disease is vital to the development of innovative diagnostic techniques, including those that may aid in the early detection and prevention of cardiovascular disorders. This dissertation focuses on the application of particle image velocimetry (PIV) to the study of intracardiac hemodynamics. This is accomplished primarily though the use of ultrasound based PIV, which allows for in vivo visualization of intracardiac flow without the requirement for optical access, as is required with traditional camera-based PIV methods.

The fundamentals of ultrasound PIV are introduced, including experimental methods for its implementation as well as a discussion on estimating and mitigating measurement error. Ultrasound PIV is then compared to optical PIV; this is a highly developed technique with proven accuracy; through rigorous examination it has become the “gold standard” of two-dimensional flow visualization. Results show good agreement between the two methods.

Using a mechanical left heart model, a multi-plane ultrasound PIV technique is introduced and applied to quantify a complex, three-dimensional flow that is analogous to the left intraventricular flow. Changes in ventricular flow dynamics due to the rotational orientation of mechanical heart valves are studied; the results demonstrate the importance of multi-plane imaging techniques when trying to assess the strongly three-dimensional intraventricular flow.

The potential use of ultrasound PIV as an early diagnosis technique is demonstrated through the development of a novel elasticity estimation technique. A finite element analysis routine is couple with an ensemble Kalman filter to allow for the estimation of material elasticity using forcing and displacement data derived from PIV. Results demonstrate that it is possible to estimate elasticity using forcing data derived from a PIV vector field, provided vector density is sufficient.

Wall-bounded turbulence manifests itself in a broad range of applications, not least of which in hydraulic systems. Here we briefly review the significant advances over the past few decades in the fundamental study of wall turbulence over smooth and rough surfaces, with an emphasis on coherent structures and their role at high Reynolds numbers. We attempt to relate these findings to parallel efforts in the hydraulic engineering community and discuss the implications of coherent structures in important hydraulic phenomena.

Vortex organization in the outer layer of a turbulent boundary layer overlying sparse, hemispherical roughness elements is explored with two-component particle-image velocimetry (PIV) in multiple streamwise-wall-normal measurement planes downstream and between elements. The presence of sparse roughness elements causes a shortening of the streamwise length scale in the near-wall region. These measurements confirm that vortex packets exist in the outer layer of flow over rough walls, but that their organization is altered, and this is interpreted as the underlying cause of the length-scale reduction. In particular, the elements shed vortices which appear to align in the near-wall region, but are distinct from the packets. Further, it is observed that ejection events triggered in the element wakes are more intense compared to the ejection events in smooth wall. We speculate that this may initiate a self-sustaining mechanism leading to the formation of hairpin packets as a much more effective instability compared to those typical of smooth-wall turbulence.