Matching Items (20)
154540-Thumbnail Image.png
Description
A moving overlapping mesh methodology that achieves spectral accuracy in space and up to second-order accuracy in time is developed for solution of unsteady incompressible flow equations in three-dimensional domains. The targeted applications are in aerospace and mechanical engineering domains and involve problems in turbomachinery, rotary aircrafts, wind turbines and

A moving overlapping mesh methodology that achieves spectral accuracy in space and up to second-order accuracy in time is developed for solution of unsteady incompressible flow equations in three-dimensional domains. The targeted applications are in aerospace and mechanical engineering domains and involve problems in turbomachinery, rotary aircrafts, wind turbines and others. The methodology is built within the dual-session communication framework initially developed for stationary overlapping meshes. The methodology employs semi-implicit spectral element discretization of equations in each subdomain and explicit treatment of subdomain interfaces with spectrally-accurate spatial interpolation and high-order accurate temporal extrapolation, and requires few, if any, iterations, yet maintains the global accuracy and stability of the underlying flow solver. Mesh movement is enabled through the Arbitrary Lagrangian-Eulerian formulation of the governing equations, which allows for prescription of arbitrary velocity values at discrete mesh points.

The stationary and moving overlapping mesh methodologies are thoroughly validated using two- and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal global convergence, for both methods, is documented and is in agreement with the nominal order of accuracy of the underlying solver.

Stationary overlapping mesh methodology was validated to assess the influence of long integration times and inflow-outflow global boundary conditions on the performance. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics are validated against the available data.

Moving overlapping mesh simulations are validated on the problems of two-dimensional oscillating cylinder and a three-dimensional rotating sphere. The aerodynamic forces acting on these moving rigid bodies are determined, and all results are compared with published data. Scaling tests, with both methodologies, show near linear strong scaling, even for moderately large processor counts.

The moving overlapping mesh methodology is utilized to investigate the effect of an upstream turbulent wake on a three-dimensional oscillating NACA0012 extruded airfoil. A direct numerical simulation (DNS) at Reynolds Number 44,000 is performed for steady inflow incident upon the airfoil oscillating between angle of attack 5.6 and 25 degrees with reduced frequency k=0.16. Results are contrasted with subsequent DNS of the same oscillating airfoil in a turbulent wake generated by a stationary upstream cylinder.
ContributorsMerrill, Brandon Earl (Author) / Peet, Yulia (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2016
155243-Thumbnail Image.png
Description
A numerical study of chemotaxis in 3D turbulence is presented here. Direct Numerical

Simulation were used to calculate the nutrient uptake for both motile and non-motile bacterial

species and by applying the dynamical systems theory the effect of flow topology on the

variability of chemotaxis is analyzed. It is done

A numerical study of chemotaxis in 3D turbulence is presented here. Direct Numerical

Simulation were used to calculate the nutrient uptake for both motile and non-motile bacterial

species and by applying the dynamical systems theory the effect of flow topology on the

variability of chemotaxis is analyzed. It is done by injecting a highly localized patch of nutrient

in the turbulent flow, and analyzing the evolution of reaction associated with the observed

high and low stretching regions. The Gaussian nutrient patch is released at different locations

and the corresponding nutrient uptake is obtained. The variable stretching characteristics of

the flow is depicted by Lagrangian Coherent Structures and the roles they play in affecting the

uptake are analyzed. The Lagrangian Coherent Structures are quantified by the Finite Time

Lyapunov Exponents which is a measure of the average stretching experienced by the flow in

finite time. It is found that in high stretching regions, the motile bacteria are attracted to the

nutrient patch very quickly, but also dispersed quickly; whereas in low stretching regions the

bacteria respond slower towards the nutrient patch. However the total uptake is intricately

determined by stretching history. These reaction characteristics are reflected in the several

realizations of simulations. This helps in understanding turbulence intensity and how it affects

the uptake of the nutrient.
ContributorsGeorge, Jino (Author) / Tang, Wenbo (Thesis advisor) / Peet, Yulia (Thesis advisor) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2017
155810-Thumbnail Image.png
Description
The effect of reduced frequency on dynamic stall behavior of a pitching NACA0012 airfoil in a turbulent wake using Direct Numerical Simulations is presented in the current study. Upstream turbulence with dynamically oscillating blades and airfoils is associated with ambient flow unsteadiness and is encountered in many operating conditions. Wake

The effect of reduced frequency on dynamic stall behavior of a pitching NACA0012 airfoil in a turbulent wake using Direct Numerical Simulations is presented in the current study. Upstream turbulence with dynamically oscillating blades and airfoils is associated with ambient flow unsteadiness and is encountered in many operating conditions. Wake turbulence, a more realistic scenario for airfoils in operation, is generated using a small solid cylinder placed upstream, the vortices shed from which interact with the pitching airfoil affecting dynamic stall behavior.

A recently developed moving overlapping grid approach is used using a high-order Spectral Element Method (SEM) for spatial discretization combined with a dynamic time-stepping procedure allowing for up to third order temporal discretization. Two cases of reduced frequency (k = 0:16 and 0:25) for airfoil oscillation are investigated and the change in dynamic stall behavior with change in reduced frequency is studied and documented using flow-fields and aerodynamic coefficients (Drag, Lift and Pitching Moment) with a focus on understanding vortex system dynamics (including formation of secondary vortices) for different reduced frequencies and it’s affect on airfoil aerodynamic characteristics and fatigue life. Transition of the flow over the surface of an airfoil for both undisturbed and disturbed flow cases will also be discussed using Pressure coefficient and Skin Friction coefficient data for a given cycle combined with a wavelet analysis using Morse wavelets in MATLAB.
ContributorsGandhi, Anurag (Author) / Peet, Yulia (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2017
155305-Thumbnail Image.png
Description
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

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.
ContributorsSakievich, Philip Sakievich (Author) / Peet, Yulia (Thesis advisor) / Adrian, Ronald (Committee member) / Squires, Kyle (Committee member) / Herrmann, Marcus (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2017
168774-Thumbnail Image.png
Description
Modern aircraft propulsion systems such as the ultra high bypass ratio turbofan impose constraints on engine installation below the wing, causing jet–wing interactions. Similar interactions are encountered when a jet-powered aircraft takes off on airport runway or aircraft carrier deck. High-speed jet flow near a solid surface shows markedly different

Modern aircraft propulsion systems such as the ultra high bypass ratio turbofan impose constraints on engine installation below the wing, causing jet–wing interactions. Similar interactions are encountered when a jet-powered aircraft takes off on airport runway or aircraft carrier deck. High-speed jet flow near a solid surface shows markedly different turbulence characteristics compared with free jet, including attached turbulent jet and development of non-equilibrium boundary layer down- stream. Wall pressure fluctuations tend to be more unsteady and stronger, leading to increased vibration affecting aircraft cabin noise and modified jet noise radiation. Large-eddy simulation (LES) is useful to characterize turbulent jet flows over a solid surface as well as wall pressure distribution to promote physical understanding and modeling studies. In this study, LES is performed for an installed setup of a Mach 0.7 turbulent jet where the jet–plate distance is fixed at 2D where D is the nozzle-exit diameter. Unstructured-grid LES is used to validate the corresponding experiment (from literature). In addition, a high-fidelity numerical database is built for further analysis and modeling. Turbulence statistics and energy spectra show that agreement with the experimental measurement for the installed case is encouraging, paving a way for future analysis and modeling.
ContributorsTamhane, Nikhil (Author) / Kim, Jeongale (Thesis advisor) / Peet, Yulia (Thesis advisor) / Jeun, Jinah (Committee member) / Arizona State University (Publisher)
Created2022
158194-Thumbnail Image.png
Description
Compressible fluid flows involving multiple physical states of matter occur in both nature and technical applications such as underwater explosions and implosions, cavitation-induced bubble collapse in naval applications and Richtmyer-Meshkov type instabilities in inertial confinement fusion. Of particular interest is the atomization of fuels that enable shock-induced mixing of fuel

Compressible fluid flows involving multiple physical states of matter occur in both nature and technical applications such as underwater explosions and implosions, cavitation-induced bubble collapse in naval applications and Richtmyer-Meshkov type instabilities in inertial confinement fusion. Of particular interest is the atomization of fuels that enable shock-induced mixing of fuel and oxidizer in supersonic combustors. Due to low residence times and varying length scales, providing insight through physical experiments is both technically challenging and sometimes unfeasible. Numerical simulations can help provide detailed insight and aid in the engineering design of devices that can harness these physical phenomena.

In this research, computational methods were developed to accurately simulate phase interfaces in compressible fluid flows with a focus on targeting primary atomization. Novel numerical methods which treat the phase interface as a discontinuity, and as a smeared region were developed using low-dissipation, high-order schemes. The resulting methods account for the effects of compressibility, surface tension and viscosity. To aid with the varying length scales and high-resolution requirements found in atomization applications, an adaptive mesh refinement (AMR) framework is used to provide high-resolution only in regions of interest. The developed methods were verified with test cases involving strong shocks, high density ratios, surface tension effects and jumps in the equations of state, in one-, two- and three dimensions, obtaining good agreement with theoretical and experimental results. An application case of the primary atomization of a liquid jet injected into a Mach 2 supersonic crossflow of air is performed with the methods developed.
ContributorsKannan, Karthik (Author) / Herrmann, Marcus (Thesis advisor) / Huang, Huei-Ping (Committee member) / Lopez, Juan (Committee member) / Peet, Yulia (Committee member) / Wang, Liping (Committee member) / Arizona State University (Publisher)
Created2020
161883-Thumbnail Image.png
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

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
171849-Thumbnail Image.png
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

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
171461-Thumbnail Image.png
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

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
156800-Thumbnail Image.png
Description
Owing to the surge in development of endovascular devices such as coils and flow diverter stents, doctors are inclined to approach surgical cases non-invasively more often than before. Treating brain aneurysms as a bulging of a weakened area of a blood vessel is no exception. Therefore, promoting techniques that can

Owing to the surge in development of endovascular devices such as coils and flow diverter stents, doctors are inclined to approach surgical cases non-invasively more often than before. Treating brain aneurysms as a bulging of a weakened area of a blood vessel is no exception. Therefore, promoting techniques that can help surgeons have a better idea of treatment outcomes are of invaluable importance.

In order to investigate the effects of these devices on intra-aneurysmal hemodynamics, the conventional computational fluid dynamics (CFD) approach uses the explicit geometry of the device within an aneurysm and discretizes the fluid domain to solve the Navier-Stokes equations. However, since the devices are made of small struts, the number of mesh elements in the boundary layer region would be considerable. This cumbersome task led to the implementation of the porous medium assumption. In this approach, the explicit geometry of the device is eliminated, and relevant porous medium assumptions are applied. Unfortunately, as it will be shown in this research, some of the porous medium approaches used in the literature are over-simplified. For example, considering the porous domain to be homogeneous is one major drawback which leads to significant errors in capturing the intra-aneurysmal flow features. Specifically, since the devices must comply with the complex geometry of an aneurysm, the homogeneity assumption is not valid.

In this research, a novel heterogeneous porous medium approach is introduced. This results in a substantial reduction in the total number of mesh elements required to discretize the flow domain while not sacrificing the accuracy of the method by over-simplifying the utilized assumptions.
ContributorsYadollahi Farsani, Hooman (Author) / Herrmann, Marcus (Thesis advisor) / Frakes, David (Thesis advisor) / Chong, Brian (Committee member) / Peet, Yulia (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2018