Matching Items (14)
Filtering by
- All Subjects: Fluid Mechanics
- All Subjects: Climate Change
- All Subjects: stereolithography
- Creators: Peet, Yulia
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
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
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
Description
This dissertation investigates the complex dynamics of semi-dilute inertial particles suspended in vortices using the Eulerian-Lagrangian method. The study explores the modulation of flow induced by inertial particles, focusing on the characteristics of a single vortex, instability analysis within particle-laden flows, and the merging process of co-rotating vortices. Simulations reveal a preferential concentration mechanism, where inertial particles cluster around a void fraction bubble, accelerating the decay of the vortex tube. Small-scale clusters, arising from particle-trajectory crossings, induce significant gradients in the fluid vorticity field, contributing to rapid vortex breakdown. Within a specific Stokes number range, increased particle inertia results in faster vortex decay and stronger inhomogeneity in the particle phase. The instability mechanism in particle-laden flows is explored using a Rankine vortex model. Two-way coupling triggers azimuthal perturbations, leading to the breakdown of the vortex structure. Linear Stability Analysis and Two-Fluid modeling demonstrate that the dusty vortex flow exhibits unstable modes, with growth rates increasing with wavenumber. Eulerian-Lagrangian simulations validate these results, showing excellent agreement between computed and predicted growth rates. The dissertation also delves into the co-rotating vortex merger in a semi-dilute dusty flow. For weak inertial effects, merger experiences a delay compared to particle-free vortices. Under moderate inertial conditions, the merger process exhibits repulsion, increased separation, and eventual convective merger stages. Highly inertial particles stretch the vortex core, initiating a merger with an outcome of a particle-free vortex core surrounded by a halo of concentrated particles. In conclusion, the feedback force from the dispersed phase induces instability and significantly influences the dynamics of vortices in particle-laden flows. The findings contribute to a deeper understanding of the intricate interactions between inertial particles and vortical structures.
ContributorsShuai, Shuai (Author) / Kasbaoui, Mohamed Houssem (Thesis advisor) / Herrmann, Marcus (Committee member) / Peet, Yulia (Committee member) / Huang, Huei-Ping (Committee member) / Wang, Zhihua (Committee member) / Arizona State University (Publisher)
Created2024
Description
This thesis presents a computational fluid dynamics (CFD) model of fluid flow driven by the motion of cilia, a cellular appendage found in organisms used to either move the fluid around them or to move themselves by propelling the fluid. Originating from an initial investigation to the flow patterns inside the third ventricle of a rat’s brain, this project expanded to improve the inadequacies of existing models of ciliary motion in fluid. This model was developed using the actuator line model to include the cilia motion to get an accurate representation of the cilia motion and its effect on the flow. This model not only provides exciting potential in various fields including soft robotics, biomedical research, environmental engineering, but also holds promise for drug delivery systems, and enhancing microfluidic designs. This thesis investigates the effect of the phase difference, the spacing and the frequency of the cilia motion on the fluid flow and the formation of the metachronal waves.
ContributorsShameer, Waseem (Author) / Peet, Yulia (Thesis advisor) / Herrmann, Marcus (Committee member) / Ma, Leixin (Committee member) / Arizona State University (Publisher)
Created2024
Description
Autonomic closure is a new general methodology for subgrid closures in large eddy simulations that circumvents the need to specify fixed closure models and instead allows a fully- adaptive self-optimizing closure. The closure is autonomic in the sense that the simulation itself determines the optimal relation at each point and time between any subgrid term and the variables in the simulation, through the solution of a local system identification problem. It is based on highly generalized representations of subgrid terms having degrees of freedom that are determined dynamically at each point and time in the simulation. This can be regarded as a very high-dimensional generalization of the dynamic approach used with some traditional prescribed closure models, or as a type of “data-driven” turbulence closure in which machine- learning methods are used with internal training data obtained at a test-filter scale at each point and time in the simulation to discover the local closure representation.
In this study, a priori tests were performed to develop accurate and efficient implementations of autonomic closure based on particular generalized representations and parameters associated with the local system identification of the turbulence state. These included the relative number of training points and bounding box size, which impact computational cost and generalizability of coefficients in the representation from the test scale to the LES scale. The focus was on studying impacts of these factors on the resulting accuracy and efficiency of autonomic closure for the subgrid stress. Particular attention was paid to the associated subgrid production field, including its structural features in which large forward and backward energy transfer are concentrated.
More than five orders of magnitude reduction in computational cost of autonomic closure was achieved in this study with essentially no loss of accuracy, primarily by using efficient frame-invariant forms for generalized representations that greatly reduce the number of degrees of freedom. The recommended form is a 28-coefficient representation that provides subgrid stress and production fields that are far more accurate in terms of structure and statistics than are traditional prescribed closure models.
In this study, a priori tests were performed to develop accurate and efficient implementations of autonomic closure based on particular generalized representations and parameters associated with the local system identification of the turbulence state. These included the relative number of training points and bounding box size, which impact computational cost and generalizability of coefficients in the representation from the test scale to the LES scale. The focus was on studying impacts of these factors on the resulting accuracy and efficiency of autonomic closure for the subgrid stress. Particular attention was paid to the associated subgrid production field, including its structural features in which large forward and backward energy transfer are concentrated.
More than five orders of magnitude reduction in computational cost of autonomic closure was achieved in this study with essentially no loss of accuracy, primarily by using efficient frame-invariant forms for generalized representations that greatly reduce the number of degrees of freedom. The recommended form is a 28-coefficient representation that provides subgrid stress and production fields that are far more accurate in terms of structure and statistics than are traditional prescribed closure models.
ContributorsKshitij, Abhinav (Author) / Dahm, Werner J.A. (Thesis advisor) / Herrmann, Marcus (Committee member) / Hamlington, Peter E (Committee member) / Peet, Yulia (Committee member) / Kim, Jeonglae (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
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
Climate change has been one of the major issues of global economic and social concerns in the past decade. To quantitatively predict global climate change, the Intergovernmental Panel on Climate Change (IPCC) of the United Nations have organized a multi-national effort to use global atmosphere-ocean models to project anthropogenically induced climate changes in the 21st century. The computer simulations performed with those models and archived by the Coupled Model Intercomparison Project - Phase 5 (CMIP5) form the most comprehensive quantitative basis for the prediction of global environmental changes on decadal-to-centennial time scales. While the CMIP5 archives have been widely used for policy making, the inherent biases in the models have not been systematically examined. The main objective of this study is to validate the CMIP5 simulations of the 20th century climate with observations to quantify the biases and uncertainties in state-of-the-art climate models. Specifically, this work focuses on three major features in the atmosphere: the jet streams over the North Pacific and Atlantic Oceans and the low level jet (LLJ) stream over central North America which affects the weather in the United States, and the near-surface wind field over North America which is relevant to energy applications. The errors in the model simulations of those features are systematically quantified and the uncertainties in future predictions are assessed for stakeholders to use in climate applications. Additional atmospheric model simulations are performed to determine the sources of the errors in climate models. The results reject a popular idea that the errors in the sea surface temperature due to an inaccurate ocean circulation contributes to the errors in major atmospheric jet streams.
ContributorsKulkarni, Sujay (Author) / Huang, Huei-Ping (Thesis advisor) / Calhoun, Ronald (Committee member) / Peet, Yulia (Committee member) / Arizona State University (Publisher)
Created2014
A level set approach for denoising and adaptively smoothing complex geometry stereolithography files
Description
Stereolithography files (STL) are widely used in diverse fields as a means of describing complex geometries through surface triangulations. The resulting stereolithography output is a result of either experimental measurements, or computer-aided design. Often times stereolithography outputs from experimental means are prone to noise, surface irregularities and holes in an otherwise closed surface.
A general method for denoising and adaptively smoothing these dirty stereolithography files is proposed. Unlike existing means, this approach aims to smoothen the dirty surface representation by utilizing the well established levelset method. The level of smoothing and denoising can be set depending on a per-requirement basis by means of input parameters. Once the surface representation is smoothened as desired, it can be extracted as a standard levelset scalar isosurface.
The approach presented in this thesis is also coupled to a fully unstructured Cartesian mesh generation library with built-in localized adaptive mesh refinement (AMR) capabilities, thereby ensuring lower computational cost while also providing sufficient resolution. Future work will focus on implementing tetrahedral cuts to the base hexahedral mesh structure in order to extract a fully unstructured hexahedra-dominant mesh describing the STL geometry, which can be used for fluid flow simulations.
A general method for denoising and adaptively smoothing these dirty stereolithography files is proposed. Unlike existing means, this approach aims to smoothen the dirty surface representation by utilizing the well established levelset method. The level of smoothing and denoising can be set depending on a per-requirement basis by means of input parameters. Once the surface representation is smoothened as desired, it can be extracted as a standard levelset scalar isosurface.
The approach presented in this thesis is also coupled to a fully unstructured Cartesian mesh generation library with built-in localized adaptive mesh refinement (AMR) capabilities, thereby ensuring lower computational cost while also providing sufficient resolution. Future work will focus on implementing tetrahedral cuts to the base hexahedral mesh structure in order to extract a fully unstructured hexahedra-dominant mesh describing the STL geometry, which can be used for fluid flow simulations.
ContributorsKannan, Karthik (Author) / Herrmann, Marcus (Thesis advisor) / Peet, Yulia (Committee member) / Frakes, David (Committee member) / Arizona State University (Publisher)
Created2014