Matching Items (9)
Filtering by
- All Subjects: Fluid Mechanics
- Creators: Adrian, Ronald
Description
A cerebral aneurysm is an abnormal ballooning of the blood vessel wall in the brain that occurs in approximately 6% of the general population. When a cerebral aneurysm ruptures, the subsequent damage is lethal damage in nearly 50% of cases. Over the past decade, endovascular treatment has emerged as an effective treatment option for cerebral aneurysms that is far less invasive than conventional surgical options. Nonetheless, the rate of successful treatment is as low as 50% for certain types of aneurysms. Treatment success has been correlated with favorable post-treatment hemodynamics. However, current understanding of the effects of endovascular treatment parameters on post-treatment hemodynamics is limited. This limitation is due in part to current challenges in in vivo flow measurement techniques. Improved understanding of post-treatment hemodynamics can lead to more effective treatments. However, the effects of treatment on hemodynamics may be patient-specific and thus, accurate tools that can predict hemodynamics on a case by case basis are also required for improving outcomes.Accordingly, the main objectives of this work were 1) to develop computational tools for predicting post-treatment hemodynamics and 2) to build a foundation of understanding on the effects of controllable treatment parameters on cerebral aneurysm hemodynamics. Experimental flow measurement techniques, using particle image velocimetry, were first developed for acquiring flow data in cerebral aneurysm models treated with an endovascular device. The experimental data were then used to guide the development of novel computational tools, which consider the physical properties, design specifications, and deployment mechanics of endovascular devices to simulate post-treatment hemodynamics. The effects of different endovascular treatment parameters on cerebral aneurysm hemodynamics were then characterized under controlled conditions. Lastly, application of the computational tools for interventional planning was demonstrated through the evaluation of two patient cases.
ContributorsBabiker, M. Haithem (Author) / Frakes, David H (Thesis advisor) / Adrian, Ronald (Committee member) / Caplan, Michael (Committee member) / Chong, Brian (Committee member) / Vernon, Brent (Committee member) / Arizona State University (Publisher)
Created2013
Description
The flow around a golf ball is studied using direct numerical simulation (DNS). An immersed boundary approach is adopted in which the incompressible Navier-Stokes equations are solved using a fractional step method on a structured, staggered grid in cylindrical coordinates. The boundary conditions on the surface are imposed using momentum forcing in the vicinity of the boundary. The flow solver is parallelized using a domain decomposition strategy and message passing interface (MPI), and exhibits linear scaling on as many as 500 processors. A laminar flow case is presented to verify the formal accuracy of the method. The immersed boundary approach is validated by comparison with computations of the flow over a smooth sphere. Simulations are performed at Reynolds numbers of 2.5 × 104 and 1.1 × 105 based on the diameter of the ball and the freestream speed and using grids comprised of more than 1.14 × 109 points. Flow visualizations reveal the location of separation, as well as the delay of complete detachment. Predictions of the aerodynamic forces at both Reynolds numbers are in reasonable agreement with measurements. Energy spectra of the velocity quantify the dominant frequencies of the flow near separation and in the wake. Time-averaged statistics reveal characteristic physical patterns in the flow as well as local trends within dimples. A mechanism of drag reduction due to the dimples is confirmed, and metrics for dimple optimization are proposed.
ContributorsSmith, Clinton E (Author) / Squires, Kyle D (Thesis advisor) / Balaras, Elias (Committee member) / Herrmann, Marcus (Committee member) / Adrian, Ronald (Committee member) / Stanzione, Daniel C (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2011
Description
The evolution of single hairpin vortices and multiple interacting hairpin vortices are studied in direct numerical simulations of channel flow at Re-tau=395. The purpose of this study is to observe the effects of increased Reynolds number and varying initial conditions on the growth of hairpins and the conditions under which single hairpins autogenerate hairpin packets. The hairpin vortices are believed to provide a unified picture of wall turbulence and play an important role in the production of Reynolds shear stress which is directly related to turbulent drag. The structures of the initial three-dimensional vortices are extracted from the two-point spatial correlation of the fully turbulent direct numerical simulation of the velocity field by linear stochastic estimation and embedded in a mean flow having the profile of the fully turbulent flow. The Reynolds number of the present simulation is more than twice that of the Re-tau=180 flow from earlier literature and the conditional events used to define the stochastically estimated single vortex initial conditions include a number of new types of events such as quasi-streamwise vorticity and Q4 events. The effects of parameters like strength, asymmetry and position are evaluated and compared with existing results in the literature. This study then attempts to answer questions concerning how vortex mergers produce larger scale structures, a process that may contribute to the growth of length scale with increasing distance from the wall in turbulent wall flows. Multiple vortex interactions are studied in detail.
ContributorsParthasarathy, Praveen Kumar (Author) / Adrian, Ronald (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2011
Description
This work helps to explain the drag reduction mechanisms at low and moderate turbulent Reynolds numbers in pipe flows. Through direct numerical simulation, the effects of wall oscillations are observed on the turbulence in both the near wall and the bulk region. Analysis of the average Reynolds Stresses at various phases of the flow is provided along with probability density functions of the fluctuating components of velocity and vorticity. The flow is also visualized to observe, qualitatively, changes in the total and fluctuating field of velocity and vorticity. Linear Stochastic Estimation is used to create a conditional eddy (associated with stress production) in the flow and visualize the effects of transverse wall oscillations on hairpin growth, auto-generation and structure.
ContributorsCoxe, Daniel (Author) / Peet, Yulia (Thesis advisor) / Adrian, Ronald (Thesis advisor) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2019
Description
Intracranial aneurysms, which form in the blood vessels of the brain, are particularly dangerous because of the importance and fragility of the human brain. When an intracranial aneurysm gets large it poses a significant risk of bursting and causing subarachnoid hemorrhaging (SAH), a possibly fatal condition. One possible treatment involves placing a stent in the vessel to act as a flow diverter. In this study we look at the hemodynamics of two geometries of idealized basilar tip aneurysms, at 2,3, and 4 ml/s pulsatile flow, at three different points in the cardiac cycle. The smaller model had neck and dome diameters of 2.67 mm and 4 mm respectively, while the larger aneurysm had neck and dome diameters of 3 mm and 6 mm respectively. Both diameters and the dome to neck ratio increased in the second model, representing growth over time. Flow was analyzed using stereoscopic particle image velocimetry (PIV) for both geometries in untreated models, as well as after treatment with a high porosity Enterprise stent (Codman and Shurtleff Inc.). Flow in the models was characterized by root mean square velocity in the aneurysm and neck plane, cross neck flow, max aneurysm vorticity, and total aneurysm kinetic energy. It was found that in the smaller aneurysm model (model 1), Enterprise stent treatment reduced all flow parameters substantially. The smallest reduction was in max vorticity, at 42.48%, and the largest in total kinetic energy, at 75.69%. In the larger model (model 2) there was a 52.18% reduction in cross neck flow, but a 167.28% increase in aneurysm vorticity. The other three parameters experienced little change. These results, along with observed velocity vector fields, indicate a noticeable diversion of flow away from the aneurysm in the stent treated model 1. Treatment in model 2 had a small flow diversion effect, but also altered flow in unpredictable ways, in some cases having a detrimental effect on aneurysm hemodynamics. The results of this study indicate that Enterprise stent treatment is only effective in small, relatively undeveloped aneurysm geometries, and waiting until an aneurysm has grown too large can eliminate this treatment option altogether.
ContributorsLindsay, James Bryan (Author) / Frakes, David (Thesis director) / LaBelle, Jeffrey (Committee member) / Nair, Priya (Committee member) / Barrett, The Honors College (Contributor) / School of Humanities, Arts, and Cultural Studies (Contributor)
Created2013-05
Description
This thesis focuses on the turbulent bluff body wakes in incompressible and compressible flows. An incompressible wake flow past an axisymmetric body of revolution at a diameter-based Reynolds number Re=5000 is investigated via a direct numerical simulation. It is followed by the development of a compressible solver using a split-form discontinuous Galerkin spectral element method framework with shock capturing. In the study on incompressible wake flows, three dominant coherent vortical motions are identified in the wake: the vortex shedding motion with the frequency of St=0.27, the bubble pumping motion with St=0.02, and the very-low-frequency (VLF) motion originated in the very near wake of the body with the frequencies St=0.002 and 0.005. The very-low-frequency motion is associated with a slow precession of the wake barycenter. The vortex shedding pattern is demonstrated to follow a reflectional symmetry breaking mode, with the detachment location rotating continuously and making a full circle over one vortex shedding period. The VLF radial motion with St=0.005 originates as m = 1 mode, but later transitions into m = 2 mode in the intermediate wake. Proper orthogonaldecomposition (POD) and dynamic mode decomposition (DMD) are further performed to analyze the spatial structure associated with the dominant coherent motions. Results of the POD and DMD analysis are consistent with the results of the azimuthal Fourier analysis. To extend the current incompressible code to be able to solve compressible flows, a computational methodology is developed using a high-order approximation for the compressible Navier-Stokes equations with discontinuities. The methodology is based on a split discretization framework with a summation-by-part operator. An entropy viscosity method and a subcell finite volume method are implemented to capture discontinuities. The developed high-order split-form with shock-capturing methodology is subject to a series of evaluation on cases from subsonic to hypersonic, from one-dimensional to three dimensional. The Taylor-Green vortex case and the supersonic sphere wake case show the capability to handle three-dimensional turbulent flows without and with the presence of shocks. It is also shown that higher-order approximations yield smaller errors than lower-order approximations, for the same number of total degrees of freedom.
ContributorsZhang, Fengrui (Author) / Peet, Yulia (Thesis advisor) / Kostelich, Eric (Committee member) / Kim, Jeonglae (Committee member) / Hermann, Marcus (Committee member) / Adrian, Ronald (Committee member) / Arizona State University (Publisher)
Created2022
Description
This work aims to address the design optimization of bio-inspired locomotive devices in collective swimming by developing a computational methodology which combines surrogate-based optimization with high fidelity fluid-structure interactions (FSI) simulations of thunniform swimmers. Three main phases highlight the contribution and novelty of the current work. The first phase includes the development and bench-marking of a constrained surrogate-based optimization algorithm which is appropriate to the current design problem. Additionally, new FSI techniques, such as a volume-conservation scheme, has been developed to enhance the accuracy and speed of the simulations. The second phase involves an investigation of the optimized hydrodynamics of a solitary accelerating self-propelled thunniform swimmer during start-up. The third phase extends the analysis to include the optimized hydrodynamics of accelerating swimmers in phalanx schools. Future work includes extending the analysis to the optimized hydrodynamics of steady-state and accelerating swimmers in a diamond-shaped school. The results of the first phase indicate that the proposed optimization algorithm maintains a competitive performance when compared to other gradient-based and gradient-free methods, in dealing with expensive simulations-based black-box optimization problems with constraints. In addition, the proposed optimization algorithm is capable of insuring strictly feasible candidates during the optimization procedure, which is a desirable property in applied engineering problems where design variables must remain feasible for simulations or experiments not to fail. The results of the second phase indicate that the optimized kinematic gait of a solitary accelerating swimmer generates the reverse Karman vortex street associated with high propulsive efficiency. Moreover, the efficiency of sub-optimum modes, in solitary swimming, is found to increase with both the tail amplitude and the effective flapping length of the swimmer, and a new scaling law is proposed to capture these trends. Results of the third phase indicate that the optimal midline kinematics in accelerating phalanx schools resemble those of accelerating solitary swimmers. The optimal separation distance in a phalanx school is shown to be around 2L (where L is the swimmer's total length). Furthermore, separation distance is shown to have a stronger effect, ceteris paribus, on the propulsion efficiency of a school when compared to phase synchronization.
ContributorsAbouhussein, Ahmed (Author) / Peet, Yulia (Thesis advisor) / Adrian, Ronald (Committee member) / Kim, Jeonglae (Committee member) / Kasbaoui, Mohamed (Committee member) / Mittelmann, Hans (Committee member) / Arizona State University (Publisher)
Created2022
Description
Realistic engineering, physical and biological systems are very complex in nature, and their response and performance are governed by multitude of interacting processes. In computational modeling of these systems, the interactive response is most often ignored, and simplifications are made to model one or a few relevant phenomena as opposed to a complete set of interacting processes due to a complexity of integrative analysis. In this thesis, I will develop new high-order computational approaches that reduce the amount of simplifications and model the full response of a complex system by accounting for the interaction between different physical processes as required for an accurate description of the global system behavior. Specifically, I will develop multi-physics coupling techniques based on spectral-element methods for the simulations of such systems. I focus on three specific applications: fluid-structure interaction, conjugate heat transfer, and modeling of acoustic wave propagation in non-uniform media. Fluid-structure interaction illustrates a complex system between a fluid and a solid, where a movable and deformable structure is surrounded by fluid flow, and its deformation caused by fluid affects the fluid flow interactively. To simulate this system, two coupling schemes are developed: 1) iterative implicit coupling, and 2) explicit coupling based on Robin-Neumann boundary conditions. A comprehensive verification strategy of the developed methodology is presented, including a comparison with benchmark flow solutions, h-, p- and temporal refinement studies. Simulation of a turbulent flow in a channel interacting with a compliant wall is attempted as well. Another problem I consider is when a solid is stationary, but a heat transfer occurs on the fluid-solid interface. To model this problem, a conjugate heat transfer framework is introduced. Validation of the framework, as well as studies of an interior thermal environment in a building regulated by an HVAC system with an on/off control model with precooling and multi-zone precooling strategies are presented. The final part of this thesis is devoted to modeling an interaction of acoustic waves with the fluid flow. The development of a spectral-element methodology for solution of Lighthill’s equation, and its application to a problem of leak detection in water pipes is presented.
ContributorsXu, Yiqin (Author) / Peet, Yulia (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Adrian, Ronald (Committee member) / Baer, Steven (Committee member) / Arizona State University (Publisher)
Created2021
Description
Four-Dimensional Emission Tomography (4DET) and Four-Dimensional Absorption Tomography (4DAT) are measurement techniques that utilize multiple 2D images (or projections) acquired via an optical device, such as a camera, to reconstruct scalar and velocity fields of a flow field being studied, using either emission- or absorption-based measurements, respectively. Turbulence is inherently three-dimensional, and thus research in the field benefits from a comprehensive understanding of coherent structures to fully explain the flow physics involved, for example, in the phenomena resulting from a turbulent jet. This thesis looks at the development, application and validity/practicality of emission tomography as an experimental approach to a obtaining a comprehensive understanding of coherent structures in turbulent flows. A pseudo test domain is decided upon, with a varying number of camera objects created to image the region of interest. Rays are then modelled as cylindrical volumes to build the weight matrix. Projection images are generated with Gaussian concentration defined as a spatial function of the domain to build the projection matrix. Finally, concentration within the domain, evaluated via the Least Squares method, is compared against original concentration values. The reconstruction algorithm is validated and checked for accuracy with DNS data of a steady turbulent jet. Reconstruction accuracy and a statistical analysis of the reconstructions are also presented.
ContributorsRodrigues, Cossack (Author) / Pathikonda, Gokul (Thesis advisor) / Grauer, Samuel (Committee member) / Adrian, Ronald (Committee member) / Kasbaoui, Mohamed (Committee member) / Kim, Jeonglae (Committee member) / Arizona State University (Publisher)
Created2023