Matching Items (9)
Filtering by
- All Subjects: Aerospace Engineering
- Creators: Adrian, Ronald
- Resource Type: Text
Description
Multiphase flows are an important part of many natural and technological phe- nomena such as ocean-air coupling (which is important for climate modeling) and the atomization of liquid fuel jets in combustion engines. The unique challenges of multiphase flow often make analytical solutions to the governing equations impos- sible and experimental investigations very difficult. Thus, high-fidelity numerical simulations can play a pivotal role in understanding these systems. This disserta- tion describes numerical methods developed for complex multiphase flows and the simulations performed using these methods. First, the issue of multiphase code verification is addressed. Code verification answers the question "Is this code solving the equations correctly?" The method of manufactured solutions (MMS) is a procedure for generating exact benchmark solutions which can test the most general capabilities of a code. The chief obstacle to applying MMS to multiphase flow lies in the discontinuous nature of the material properties at the interface. An extension of the MMS procedure to multiphase flow is presented, using an adaptive marching tetrahedron style algorithm to compute the source terms near the interface. Guidelines for the use of the MMS to help locate coding mistakes are also detailed. Three multiphase systems are then investigated: (1) the thermocapillary motion of three-dimensional and axisymmetric drops in a confined apparatus, (2) the flow of two immiscible fluids completely filling an enclosed cylinder and driven by the rotation of the bottom endwall, and (3) the atomization of a single drop subjected to a high shear turbulent flow. The systems are simulated numerically by solving the full multiphase Navier- Stokes equations coupled to the various equations of state and a level set interface tracking scheme based on the refined level set grid method. The codes have been parallelized using MPI in order to take advantage of today's very large parallel computational architectures. In the first system, the code's ability to handle surface tension and large tem- perature gradients is established. In the second system, the code's ability to sim- ulate simple interface geometries with strong shear is demonstrated. In the third system, the ability to handle extremely complex geometries and topology changes with strong shear is shown.
ContributorsBrady, Peter, Ph.D (Author) / Herrmann, Marcus (Thesis advisor) / Lopez, Juan (Thesis advisor) / Adrian, Ronald (Committee member) / Calhoun, Ronald (Committee member) / Chen, Kangping (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 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 thesis encompasses research performed in the focus area of structural health monitoring. More specifically, this research focuses on high velocity impact testing of carbon fiber reinforced structures, especially plates, and evaluating the damage post-impact. To this end, various non-destructive evaluation techniques such as ultrasonic C-scan testing and flash thermography were utilized for post-impact analysis. MATLAB algorithms were written and refined for the localization and quantification of damage in plates using data from sensors such as piezoelectric and fiber Bragg gratings sensors. Throughout the thesis, the general plate theory and laminate plate theory, the operations and optimization of the gas gun, and the theory used for the damage localization algorithms will be discussed. Additional quantifiable results are to come in future semesters of experimentation, but this thesis outlines the framework upon which all the research will continue to advance.
ContributorsMccrea, John Patrick (Author) / Chattopadhyay, Aditi (Thesis director) / Borkowski, Luke (Committee member) / Barrett, The Honors College (Contributor) / Mechanical and Aerospace Engineering Program (Contributor) / Department of Military Science (Contributor)
Created2015-05
Description
Advanced aerospace materials, including fiber reinforced polymer and ceramic matrix composites, are increasingly being used in critical and demanding applications, challenging the current damage prediction, detection, and quantification methodologies. Multiscale computational models offer key advantages over traditional analysis techniques and can provide the necessary capabilities for the development of a comprehensive virtual structural health monitoring (SHM) framework. Virtual SHM has the potential to drastically improve the design and analysis of aerospace components through coupling the complementary capabilities of models able to predict the initiation and propagation of damage under a wide range of loading and environmental scenarios, simulate interrogation methods for damage detection and quantification, and assess the health of a structure. A major component of the virtual SHM framework involves having micromechanics-based multiscale composite models that can provide the elastic, inelastic, and damage behavior of composite material systems under mechanical and thermal loading conditions and in the presence of microstructural complexity and variability. Quantification of the role geometric and architectural variability in the composite microstructure plays in the local and global composite behavior is essential to the development of appropriate scale-dependent unit cells and boundary conditions for the multiscale model. Once the composite behavior is predicted and variability effects assessed, wave-based SHM simulation models serve to provide knowledge on the probability of detection and characterization accuracy of damage present in the composite. The research presented in this dissertation provides the foundation for a comprehensive SHM framework for advanced aerospace materials. The developed models enhance the prediction of damage formation as a result of ceramic matrix composite processing, improve the understanding of the effects of architectural and geometric variability in polymer matrix composites, and provide an accurate and computational efficient modeling scheme for simulating guided wave excitation, propagation, interaction with damage, and sensing in a range of materials. The methodologies presented in this research represent substantial progress toward the development of an accurate and generalized virtual SHM framework.
ContributorsBorkowski, Luke (Author) / Chattopadhyay, Aditi (Thesis advisor) / Liu, Yongming (Committee member) / Mignolet, Marc (Committee member) / Papandreou-Suppappola, Antonia (Committee member) / Rajadas, John (Committee member) / Arizona State University (Publisher)
Created2015
Description
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.
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.
ContributorsJibben, Zechariah J (Author) / Herrmann, Marcus (Thesis advisor) / Squires, Kyle (Committee member) / Adrian, Ronald (Committee member) / Chen, Kangping (Committee member) / Treacy, Michael (Committee member) / Arizona State University (Publisher)
Created2015
Description
Advanced fibrous composite materials exhibit outstanding thermomechanical performance under extreme environments, which make them ideal for structural components that are used in a wide range of aerospace, nuclear, and defense applications. The integrity and residual useful life of these components, however, are strongly influenced by their inherent material flaws and defects resulting from the complex fabrication processes. These defects exist across multiple length scales and govern several scale-dependent inelastic deformation mechanisms of each of the constituents as well as their composite damage anisotropy. Tailoring structural components for optimal performance requires addressing the knowledge gap regarding the microstructural material morphology that governs the structural scale damage and failure response. Therefore, there is a need for a high-fidelity multiscale modeling framework and scale-specific in-situ experimental characterization that can capture complex inelastic mechanisms, including damage initiation and propagation across multiple length scales. This dissertation presents a novel multiscale computational framework that accounts for experimental information pertinent to microstructure morphology and architectural variabilities to investigate the response of ceramic matrix composites (CMCs) with manufacturing-induced defects. First, a three-dimensional orthotropic viscoplasticity creep formulation is developed to capture the complex temperature- and time-dependent constituent load transfer mechanisms in different CMC material systems. The framework also accounts for a reformulated fracture mechanics-informed matrix damage model and the Curtin progressive fiber damage model to capture the complex scale-dependent damage and failure mechanisms through crack kinetics and porosity growth. Next, in-situ experiments using digital image correlation (DIC) are performed to capture the damage and failure mechanisms in CMCs and to validate the high-fidelity modeling results.
The dissertation also presents an exhaustive experimental investigation into the effects of temperature and manufacturing-induced defects on toughened epoxy adhesives and hybrid composite-metallic bonded joints. Nondestructive evaluation techniques are utilized to characterize the inherent defects morphology of the bulk adhesives and bonded interface. This is followed by quasi-static tensile tests conducted at extreme hot and cold temperature conditions. The damage mechanisms and failure modes are investigated using in-situ DIC and a high-resolution camera. The information from the morphology characterization studies is used to reconstruct high-fidelity geometries of the test specimens for finite element analysis.
ContributorsKhafagy, Khaled Hassan Abdo (Author) / Chattopadhyay, Aditi (Thesis advisor) / Fard, Masoud Y. (Committee member) / Milcarek, Ryan (Committee member) / Stoumbos, Tom (Committee member) / Borkowski, Luke (Committee member) / Arizona State University (Publisher)
Created2022
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
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