Matching Items (30)
Filtering by
- All Subjects: Computational Fluid Dynamics
- Creators: Herrmann, Marcus
- Creators: Harrington Bioengineering Program
- Creators: Kostelich, Eric
Description
A new method of adaptive mesh generation for the computation of fluid flows is investigated. The method utilizes gradients of the flow solution to adapt the size and stretching of elements or volumes in the computational mesh as is commonly done in the conventional Hessian approach. However, in the new method, higher-order gradients are used in place of the Hessian. The method is applied to the finite element solution of the incompressible Navier-Stokes equations on model problems. Results indicate that a significant efficiency benefit is realized.
ContributorsShortridge, Randall (Author) / Chen, Kang Ping (Thesis advisor) / Herrmann, Marcus (Thesis advisor) / Wells, Valana (Committee member) / Huang, Huei-Ping (Committee member) / Mittelmann, Hans (Committee member) / Arizona State University (Publisher)
Created2011
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
The purpose of this thesis was to explore how changes in the geometry of a bifurcating cerebral aneurysm will affect the hemodynamics in idealized models after stent treatment. This thesis explores the use of a computationally modeled Enterprise Vascular Reconstruction Device (Cordis, East Bridgewater, NJ), a high porosity and closed cell design. The models represent idealized cases of saccular aneurysms with dome sizes of either 4mm or 6mm and a dome to neck ratio of either 3:2 or 2:1. Two aneurysm contact angles are studied, one at 45 degrees and the other at 90 degrees. The stent was characterized and deployed with the use of Finite Element Analysis into each model. Computational Fluid Dynamic principles were applied in series of simulations on treated and untreated models. Data was gathered in the neck plane for the average velocity magnitude, root mean squared velocity, average flow vector angle of deflection, and the cross neck flow rate. Within the aneurysm, the average velocity magnitude, root mean squared velocity, and average pressure were calculated. Additionally, the mass flow rate at each outlet was recorded. The results of this study indicate that the Enterprise Stent was most effective in the sharper, 90 degree geometry of Model 3. Additionally, the stent had an adverse effect on the Models 1 and 4, which had the smallest neck sizes. Conclusions are that the Enterprise Stent, as a stand-alone treatment method is only reliable in situations that take advantage of its design.
ContributorsThomas, Kyle Andrew (Author) / Frakes, David (Thesis director) / LaBelle, Jeffrey (Committee member) / Babiker, Haithem (Committee member) / Barrett, The Honors College (Contributor) / Harrington Bioengineering Program (Contributor)
Created2013-05
Description
This thesis outlines the development of a vector retrieval technique, based on data assimilation, for a coherent Doppler LIDAR (Light Detection and Ranging). A detailed analysis of the Optimal Interpolation (OI) technique for vector retrieval is presented. Through several modifications to the OI technique, it is shown that the modified technique results in significant improvement in velocity retrieval accuracy. These modifications include changes to innovation covariance portioning, covariance binning, and analysis increment calculation. It is observed that the modified technique is able to make retrievals with better accuracy, preserves local information better, and compares well with tower measurements. In order to study the error of representativeness and vector retrieval error, a lidar simulator was constructed. Using the lidar simulator a thorough sensitivity analysis of the lidar measurement process and vector retrieval is carried out. The error of representativeness as a function of scales of motion and sensitivity of vector retrieval to look angle is quantified. Using the modified OI technique, study of nocturnal flow in Owens' Valley, CA was carried out to identify and understand uncharacteristic events on the night of March 27th 2006. Observations from 1030 UTC to 1230 UTC (0230 hr local time to 0430 hr local time) on March 27 2006 are presented. Lidar observations show complex and uncharacteristic flows such as sudden bursts of westerly cross-valley wind mixing with the dominant up-valley wind. Model results from Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS®) and other in-situ instrumentations are used to corroborate and complement these observations. The modified OI technique is used to identify uncharacteristic and extreme flow events at a wind development site. Estimates of turbulence and shear from this technique are compared to tower measurements. A formulation for equivalent wind speed in the presence of variations in wind speed and direction, combined with shear is developed and used to determine wind energy content in presence of turbulence.
ContributorsChoukulkar, Aditya (Author) / Calhoun, Ronald (Thesis advisor) / Mahalov, Alex (Committee member) / Kostelich, Eric (Committee member) / Huang, Huei-Ping (Committee member) / Phelan, Patrick (Committee member) / Arizona State University (Publisher)
Created2013
Description
Next generation gas turbines will be required to produce low concentrations of pollutants such as oxides of nitrogen (NOx), carbon monoxide (CO), and soot. In order to design gas turbines which produce lower emissions it is essential to have computational tools to help designers. Over the past few decades, computational fluid dynamics (CFD) has played a key role in the design of turbomachinary and will be heavily relied upon for the design of future components. In order to design components with the least amount of experimental rig testing, the ensemble of submodels used in simulations must be known to accurately predict the component's performance. The present work aims to validate a CFD model used for a reverse flow, rich-burn, quick quench, lean-burn combustor being developed at Honeywell. Initially, simulations are performed to establish a baseline which will help to assess impact to combustor performance made by changing CFD models. Rig test data from Honeywell is compared to these baseline simulation results. Reynolds averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) turbulence models are both used with the presumption that the LES turbulence model will better predict combustor performance. One specific model, the fuel spray model, is evaluated next. Experimental data of the fuel spray in an isolated environment is used to evaluate models for the fuel spray and a new, simpler approach for inputting the spray boundary conditions (BC) in the combustor is developed. The combustor is simulated once more to evaluate changes from the new fuel spray boundary conditions. This CFD model is then used in a predictive simulation of eight other combustor configurations. All computer simulations in this work were preformed with the commercial CFD software ANSYS FLUENT. NOx pollutant emissions are predicted reasonably well across the range of configurations tested using the RANS turbulence model. However, in LES, significant under predictions are seen. Causes of the under prediction in NOx concentrations are investigated. Temperature metrics at the exit of the combustor, however, are seen to be better predicted with LES.
ContributorsSpencer, A. Jeffrey (Author) / Herrmann, Marcus (Thesis advisor) / Chen, Kangping (Committee member) / Adrian, Ronald (Committee member) / Arizona State University (Publisher)
Created2012
Description
This study performs numerical modeling for the climate of semi-arid regions by running a high-resolution atmospheric model constrained by large-scale climatic boundary conditions, a practice commonly called climate downscaling. These investigations focus especially on precipitation and temperature, quantities that are critical to life in semi-arid regions. Using the Weather Research and Forecast (WRF) model, a non-hydrostatic geophysical fluid dynamical model with a full suite of physical parameterization, a series of numerical sensitivity experiments are conducted to test how the intensity and spatial/temporal distribution of precipitation change with grid resolution, time step size, the resolution of lower boundary topography and surface characteristics. Two regions, Arizona in U.S. and Aral Sea region in Central Asia, are chosen as the test-beds for the numerical experiments: The former for its complex terrain and the latter for the dramatic man-made changes in its lower boundary conditions (the shrinkage of Aral Sea). Sensitivity tests show that the parameterization schemes for rainfall are not resolution-independent, thus a refinement of resolution is no guarantee of a better result. But, simulations (at all resolutions) do capture the inter-annual variability of rainfall over Arizona. Nevertheless, temperature is simulated more accurately with refinement in resolution. Results show that both seasonal mean rainfall and frequency of extreme rainfall events increase with resolution. For Aral Sea, sensitivity tests indicate that while the shrinkage of Aral Sea has a dramatic impact on the precipitation over the confine of (former) Aral Sea itself, its effect on the precipitation over greater Central Asia is not necessarily greater than the inter-annual variability induced by the lateral boundary conditions in the model and large scale warming in the region. The numerical simulations in the study are cross validated with observations to address the realism of the regional climate model. The findings of this sensitivity study are useful for water resource management in semi-arid regions. Such high spatio-temporal resolution gridded-data can be used as an input for hydrological models for regions such as Arizona with complex terrain and sparse observations. Results from simulations of Aral Sea region are expected to contribute to ecosystems management for Central Asia.
ContributorsSharma, Ashish (Author) / Huang, Huei-Ping (Thesis advisor) / Adrian, Ronald (Committee member) / Herrmann, Marcus (Committee member) / Phelan, Patrick E. (Committee member) / Vivoni, Enrique (Committee member) / Arizona State University (Publisher)
Created2012
Description
This study investigates the application of Computational Fluid Dynamics (CFD) to the medical field. An overview of recent advances in computational simulation and modeling in medical applications is provided, with a particular emphasis on CFD. This study attempts to validate CFD and demonstrate the possibility for applying CFD to the clinical treatment and evaluation of atherosclerotic disease. Three different geometric configurations are investigated: one idealized bifurcation with a primary diameter of 8 mm, and two different patient-specific models of the bifurcation from the common femoral artery to the superficial and deep femoral arteries. CFD is compared against experimental measurements of steady state and pulsatile flow acquired with Particle Image Velocimetry (PIV). Steady state and pulsatile flow rates that are consistent with those observed in the femoral artery are used. In addition, pulsatile CFD simulations are analyzed in order to demonstrate meaningful clinical applications for studying and evaluating the treatment of atherosclerotic disease. CFD was successfully validated for steady state flow, with an average percent error of 6.991%. Potential for validation was also demonstrated for pulsatile flow, but methodological errors warrant further investigation to reformulate methods and analyze results. Quantities frequently associated with atherosclerotic disease and arterial bifurcations, such as large variations in wall shear stress and the presence of recirculation zones are demonstrated from the pulsatile CFD simulations. Further study is required in order to evaluate whether or not such phenomena are represented by CFD accurately. Further study must also be performed in order to evaluate the practicality and utility of CFD for the evaluation of atherosclerotic disease treatment.
ContributorsMortensen, Matthew James (Author) / VanAuker, Michael (Thesis director) / Frakes, David (Committee member) / Harrington Bioengineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
Dividing the plane in half leaves every border point of one region a border point of both regions. Can we divide up the plane into three or more regions such that any point on the boundary of at least one region is on the border of all the regions? In fact, it is possible to design a dynamical system for which the basins of attractions have this Wada property. In certain circumstances, both the Hénon map, a simple system, and the forced damped pendulum, a physical model, produce Wada basins.
ContributorsWhitehurst, Ryan David (Author) / Kostelich, Eric (Thesis director) / Jones, Donald (Committee member) / Armbruster, Dieter (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Department of Chemistry and Biochemistry (Contributor)
Created2013-05
Description
Morphine is a commonly used analgesic in pain management. Opioid administration to a patient after surgery, such as spinal decompression surgery, can lead to adverse side effects. To demonstrate these adverse side effects could be decreased we created a model of how morphine and its metabolites are transported and excreted from the body. Using the of morphine and a standard compartment approach this thesis aimed at projecting pharmacokinetics trends of morphine overtime. A Matlab compartment model predicting the transport of morphine through the body can contribute to a better understanding of the concentrations at the systemic level, specifically with respect to a CSF, and what happens when you compare an intravenous injection to a local delivery. Other studies and models commonly utilized patient data over small periods of time2,3,5. An extended period of time will provide information into morphine’s time course after surgery. This model focuses on a compartmentalization of the major organs and the use of a simple Mechalis-Menten enzyme kinetics for the metabolites in the liver. Our results show a CSF concentration of about 1.086×〖10〗^(-12) nmol/L in 6 weeks and 1.0097×〖10〗^(-12) nmol/L in 12 weeks. The concentration profiles in this model are similar to what was expected. The implications of this suggest that patients who reported effects of morphine paste, a locally administered opioid, weeks after the surgery were due to other reasons. In creating a model we can determine important variables and dosage information. This information allows for a greater understanding of what is happening in the body and how to improve surgical outcomes. We propose this study has implications in general research in the pharmacokinetics and dynamics of pharmacology through the body.
ContributorsJacobs, Danielle Renee (Author) / Caplan, Michael (Thesis director) / Giers, Morgan (Committee member) / Barrett, The Honors College (Contributor) / Harrington Bioengineering Program (Contributor)
Created2014-05