Theses and Dissertations
This collection includes both ASU Theses and Dissertations, submitted by graduate students, and the Barrett, Honors College theses submitted by undergraduate students.
Displaying 1 - 10 of 17
Filtering by
- All Subjects: Computational Fluid Dynamics
- Creators: Herrmann, Marcus
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
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
Atomization of fluids inside combustion chamber has been a very complex and long-lasting subject that is still researched into for maximum efficiency in mixing oxidizer and fuel. This thesis focuses on an injector called the Liquid-Liquid Swirl Coaxial Injector (LLSC) to be used in a small-scale rocket engine due to its high efficiency in spray angles and low pressure drops. Injectors are the elements that exist as a connection in between the plumbing and the combustion chamber of the rocket engine. The performance of injectors can greatly affect the stability and efficiency of the engine. Injectors proportionally help breakup the fluid into small droplets that help in the efficiency of vaporization of fluids while combusting. Helios Rocketry, Arizona State University’s student-led engineering organization, is working to design and successfully launch a small-scale bi-propellant liquid rocket engine to a 100 km (Karman Line) in space as part of the Base11 challenge. For this task a highly efficient injector element needed to be designed that can achieve high amounts of atomization with a large spray angle, to help with combustion in a relatively small sized chamber. The purpose of this thesis is to explore a specific type of injector element called a LLSC injector element. This is performed by simulating it through an LES model in computational fluid dynamics using a Voronoi based meshing scheme, by using codes from Cascade Technologies. In the end a 35-injector element design was used for an injector plate. This helped minimize the pressure drop and keep the wall stress below the safety limit.
ContributorsDave, Himanshu Hitendra (Author) / Herrmann, Marcus (Thesis director) / Adrian, Ronald (Committee member) / Mechanical and Aerospace Engineering Program (Contributor, Contributor) / Barrett, The Honors College (Contributor)
Created2019-05
Description
Advancements to a dual scale Large Eddy Simulation (LES) modeling approach for immiscible turbulent phase interfaces are presented. In the dual scale LES approach, a high resolution auxiliary grid, used to capture a fully resolved interface geometry realization, is linked to an LES grid that solves the filtered Navier-Stokes equations. Exact closure of the sub-filter interface terms is provided by explicitly filtering the fully resolved quantities from the auxiliary grid. Reconstructing a fully resolved velocity field to advance the phase interface requires modeling several sub-filter effects, including shear and accelerational instabilities and phase change. Two sub-filter models were developed to generate these sub-filter hydrodynamic instabilities: an Orr-Sommerfeld model and a Volume-of-Fluid (VoF) vortex sheet method. The Orr-Sommerfeld sub-filter model was found to be incompatible with the dual scale approach, since it is unable to generate interface rollup and a process to separate filtered and sub-filter scales could not be established. A novel VoF vortex sheet method was therefore proposed, since prior vortex methods have demonstrated interface rollup and following the LES methodology, the vortex sheet strength could be decomposed into its filtered and sub-filter components. In the development of the VoF vortex sheet method, it was tested with a variety of classical hydrodynamic instability problems, compared against prior work and linear theory, and verified using Direct Numerical Simulations (DNS). An LES consistent approach to coupling the VoF vortex sheet with the LES filtered equations is presented and compared against DNS. Finally, a sub-filter phase change model is proposed and assessed in the dual scale LES framework with an evaporating interface subjected to decaying homogeneous isotropic turbulence. Results are compared against DNS and the interplay between surface tension forces and evaporation are discussed.
ContributorsGoodrich, Austin Chase (Author) / Herrmann, Marcus (Thesis advisor) / Dahm, Werner (Committee member) / Kim, Jeonglae (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2023
Description
An interface reconstruction algorithm for the Volume of Fluid (VOF) method is required for two-phase flow problems for advection of phase interface. The primary method for interface reconstruction has been through piecewise linear interface calculation (PLIC) reconstruction. However, while PLIC reconstruction is highly accurate at representing small curvature interfaces by approximating planes across multiple grid cells, accuracy problems arise when the size of the mesh is too coarse to accurately approximate a large curvature without resorting to refining the mesh. An elliptic interface reconstructing algorithm is explored for two-phase flow problems in 2D to determine the viability of a higher-order interface reconstruction algorithm. This requires first developing an area overlap function between an arbitrary triangle and ellipse, which is then extended to calculate the area fraction field of an ellipse within a mesh. Then, the "reverse" problem of elliptic interface reconstruction given an area fraction field is examined. A study is conducted to determine the presence of any local minimums when varying the ellipse parameters. In the future, a multi-dimensional root-finding solver using Newton's Method will be developed to properly reconstruct the elliptic interface given the area fraction field.
ContributorsLee, Chase (Author) / Herrmann, Marcus (Thesis director) / Kasbaoui, Mohamed (Committee member) / Wells, Valana (Committee member) / Barrett, The Honors College (Contributor) / Mechanical and Aerospace Engineering Program (Contributor)
Created2022-05
Description
This study identifies the influence that leading-edge shape has on the aerodynamic characteristics of a wing using surface far-field and near-field analysis. It examines if a wake survey is the appropriate means for measuring profile drag and induced drag. The paper unveils the differences between sharp leading-edge and blunt leading-edge wings with the tools of pressure loop, chordwise pressure distribution, span load plots and with wake integral computations. The analysis was performed using Computational Fluid Dynamics (CFD), vortex lattice potential flow code (VORLAX), and a few wind-tunnels runs to acquire data for comparison. This study found that sharp leading-edge wings have less leading-edge suction and higher drag than blunt leading-edge wings.
The blunt leading-edge wings have less drag because the normal vector of the surface in the front section of the airfoil develops forces at opposed skin friction. The shape of the leading edge, in conjunction with the effect of viscosity, slightly alter the span load; both the magnitude of the lift and the transverse distribution. Another goal in this study is to verify the veracity of wake survey theory; the two different leading-edge shapes reveals the shortcoming of Mclean’s equation which is only applicable to blunt leading-edge wings.
The blunt leading-edge wings have less drag because the normal vector of the surface in the front section of the airfoil develops forces at opposed skin friction. The shape of the leading edge, in conjunction with the effect of viscosity, slightly alter the span load; both the magnitude of the lift and the transverse distribution. Another goal in this study is to verify the veracity of wake survey theory; the two different leading-edge shapes reveals the shortcoming of Mclean’s equation which is only applicable to blunt leading-edge wings.
ContributorsOu, Che Wei (Author) / Takahashi, Timothy (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
Description
This dissertation introduces FARCOM (Fortran Adaptive Refiner for Cartesian Orthogonal Meshes), a new general library for adaptive mesh refinement (AMR) based on an unstructured hexahedral mesh framework. As a result of the underlying unstructured formulation, the refinement and coarsening operators of the library operate on a single-cell basis and perform in-situ replacement of old mesh elements. This approach allows for h-refinement without the memory and computational expense of calculating masked coarse grid cells, as is done in traditional patch-based AMR approaches, and enables unstructured flow solvers to have access to the automated domain generation capabilities usually only found in tree AMR formulations.
The library is written to let the user determine where to refine and coarsen through custom refinement selector functions for static mesh generation and dynamic mesh refinement, and can handle smooth fields (such as level sets) or localized markers (e.g. density gradients). The library was parallelized with the use of the Zoltan graph-partitioning library, which provides interfaces to both a graph partitioner (PT-Scotch) and a partitioner based on Hilbert space-filling curves. The partitioned adjacency graph, mesh data, and solution variable data is then packed and distributed across all MPI ranks in the simulation, which then regenerate the mesh, generate domain decomposition ghost cells, and create communication caches.
Scalability runs were performed using a Leveque wave propagation scheme for solving the Euler equations. The results of simulations on up to 1536 cores indicate that the parallel performance is highly dependent on the graph partitioner being used, and differences between the partitioners were analyzed. FARCOM is found to have better performance if each MPI rank has more than 60,000 cells.
The library is written to let the user determine where to refine and coarsen through custom refinement selector functions for static mesh generation and dynamic mesh refinement, and can handle smooth fields (such as level sets) or localized markers (e.g. density gradients). The library was parallelized with the use of the Zoltan graph-partitioning library, which provides interfaces to both a graph partitioner (PT-Scotch) and a partitioner based on Hilbert space-filling curves. The partitioned adjacency graph, mesh data, and solution variable data is then packed and distributed across all MPI ranks in the simulation, which then regenerate the mesh, generate domain decomposition ghost cells, and create communication caches.
Scalability runs were performed using a Leveque wave propagation scheme for solving the Euler equations. The results of simulations on up to 1536 cores indicate that the parallel performance is highly dependent on the graph partitioner being used, and differences between the partitioners were analyzed. FARCOM is found to have better performance if each MPI rank has more than 60,000 cells.
ContributorsBallesteros, Carlos Alberto (Author) / Herrmann, Marcus (Thesis advisor) / Adrian, Ronald (Committee member) / Chen, Kangping (Committee member) / Huang, Huei-Ping (Committee member) / Lopez, Juan (Committee member) / Arizona State University (Publisher)
Created2019
Description
The goal of this paper was to do an analysis of two-dimensional unsplit mass and momentum conserving Finite Volume Methods for Advection for Volume of Fluid Fields with interfaces and validating their rates of convergence. Specifically three unsplit transport methods and one split transport method were amalgamated individually with four Piece-wise Linear Reconstruction Schemes (PLIC) i.e. Unsplit Eulerian Advection (UEA) by Owkes and Desjardins (2014), Unsplit Lagrangian Advection (ULA) by Yang et al. (2010), Split Lagrangian Advection (SLA) by Scardovelli and Zaleski (2003) and Unsplit Averaged Eulerian-Lagrangian Advection (UAELA) with two Finite Difference Methods by Parker and Youngs (1992) and two Error Minimization Methods by Pilliod Jr and Puckett (2004). The observed order of accuracy was first order in all cases except when unsplit methods and error minimization methods were used consecutively in each iteration, which resulted in second-order accuracy on the shape error convergence. The Averaged Unsplit Eulerian-Lagrangian Advection (AUELA) did produce first-order accuracy but that was due to a temporal error in the numerical setup. The main unsplit methods, Unsplit Eulerian Advection (UEA) and Unsplit Lagrangian Advection (ULA), preserve mass and momentum and require geometric clipping to solve two-phase fluid flows. The Unsplit Lagrangian Advection (ULA) can allow for small divergence in the velocity field perhaps saving time on the iterative solver of the variable coefficient Poisson System.
ContributorsAnsari, Adil (M.S.) (Author) / Herrmann, Marcus (Thesis advisor) / Peet, Yulia (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019