Matching Items (18)
Description
This thesis focuses on studying the interaction between floating objects and an air-water flow system driven by gravity. The system consists of an inclined channel in which a gravity driven two phase flow carries a series of floating solid objects downstream. Numerical simulations of such a system requires the solution of not only the basic Navier-Stokes equation but also dynamic interaction between the solid body and the two-phase flow. In particular, this requires embedding of dynamic mesh within the two-phase flow. A computational fluid dynamics solver, ANSYS fluent, is used to solve this problem. Also, the individual components for these simulations are already available in the solver, few examples exist in which all are combined. A series of simulations are performed by varying the key parameters, including density of floating objects and mass flow rate at the inlet. The motion of the floating objects in those simulations are analyzed to determine the stability of the coupled flow-solid system. The simulations are successfully performed over a broad range of parametric values. The numerical framework developed in this study can potentially be used in applications, especially in assisting the design of similar gravity driven systems for transportation in manufacturing processes. In a small number of the simulations, two kinds of numerically instability are observed. One is characterized by a sudden vertical acceleration of the floating object due to a strong imbalance of the force acting on the body, which occurs when the mass flow of water is weak. The other is characterized by a sudden vertical movement of air-water interface, which occurs when two floating objects become too close together. These new types of numerical instability deserve future studies and clarifications. This study is performed only for a 2-D system. Extension of the numerical framework to a full 3-D setting is recommended as future work.
ContributorsMangavelli, Sai Chaitanya (Author) / Huang, Huei-Ping (Thesis advisor) / Kim, Jeonglae (Committee member) / Forzani, Erica (Committee member) / Arizona State University (Publisher)
Created2018
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
Particle Image Velocimetry (PIV) has become a cornerstone of modern experimental fluid mechanics due to its unique ability to resolve the entire instantaneous two-dimensional velocity field of an experimental flow. However, this methodology has historically been omitted from undergraduate curricula due to the significant cost of research-grade PIV systems and safety considerations stemming from the high-power Nd-YAG lasers typically implemented by PIV systems. In the following undergraduate thesis, a low-cost model of a PIV system is designed to be used within the context of an undergraduate fluid mechanics lab. The proposed system consists of a Hele-Shaw water tunnel, a high-power LED lighting source, and a modern smartphone camera. Additionally, a standalone application was developed to perform the necessary image processing as well as to perform Particle Streak Velocimetry (PSV) and PIV image analysis. Ultimately, the proposed system costs $229.33 and can replicate modern PIV techniques albeit for simple flow scenarios.
ContributorsZamora, Matthew Alan (Author) / Adrian, Ronald (Thesis director) / Kim, Jeonglae (Committee member) / Mechanical and Aerospace Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2021-05
Description
The current work aims to understand the influence of particles on scalar transport in particle-laden turbulent jets using point-particle direct numerical simulations (DNS). Such turbulence phenomena are observed in many applications, such as aircraft and rocket engines (e.g., engines operating in dusty environments and when close to the surface) and geophysical flows (sediment-laden rivers discharging nutrients into the oceans), etc.This thesis looks at systematically understanding the fundamental interplay between (1) fluid turbulence, (2) inertial particles, and (3) scalar transport. This work considers a temporal jet of Reynolds number of 5000 filled with the point-particles and the influence of Stokes number (St). Three Stokes numbers, St = 1, 7.5, and 20, were considered for the current work. The simulations were solved using the NGA solver, which solves the Navier-Stokes, advection-diffusion, and particle transport equations.
The statistical analysis of the mean and turbulence quantities, along with the Reynolds stresses, are estimated for the fluid and particle phases throughout the domain. The observations do not show a significant influence of St in the mean flow evolution of fluid, scalar, and particle phases. The scalar mixture fraction variance and the turbulent kinetic energy (TKE) increase slightly for the St = 1 case, compared to the particle-free and higher St cases, indicating that an optimal St exists for which the scalar variation increases. The current preliminary study establishes that the scalar variance is influenced by particles under the optimal particle St. Directions for future studies based on the current observations are presented.
ContributorsPaturu, Venkata Sai Sushant (Author) / Pathikonda, Gokul (Thesis advisor) / Kasbaoui, Mohamed Houssem (Committee member) / Kim, Jeonglae (Committee member) / Prabhakaran, Prasanth (Committee member) / Arizona State University (Publisher)
Created2023
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
Theoretical analyses of liquid atomization (bulk to droplet conversion) and turbulence have potential to advance the computability of these flows. Instead of relying on full computations or models, fundamental conservation equations can be manipulated to generate partial or full solutions. For example, integral form of the mass and energy for spray flows leads to an explicit relationship between the drop size and liquid velocities. This is an ideal form to integrate with existing computational fluid dynamic (CFD), which is well developed to solve for the liquid velocities, i.e., the momentum equation(s). Theoretical adaption to CFD has been performed for various injection geometries, with results that compare quite well with experimental data. Since the drop size is provided analytically, computational time/cost for simulating spray flows with liquid atomization is no more than single-phase flows. Some advances have also been made on turbulent flows, by using a new set of perspectives on transport, scaling and energy distributions. Conservation equations for turbulence momentum and kinetic energy have been derived in a coordinate frame moving with the local mean velocities, which produce the Reynolds stress components, without modeling. Scaling of the Reynolds stress is also found at the first- and second-gradient levels. Finally, maximum-entropy principle has been used to derive the energy spectra in turbulent flows.
ContributorsPark, Jung Eun (Author) / Lee, Taewoo (Thesis advisor) / Gardner, Carl (Committee member) / Huang, Huei-Ping (Committee member) / Kim, Jeonglae (Committee member) / Chen, Kangping (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
The study tested the parameterized neural ordinary differential equation (PNODE) framework with a physical system exhibiting only advective phenomenon. Existing deep learning methods have difficulty learning multiple dynamic, continuous time processes. PNODE encodes the input data and initial parameter into a set of reduced states within the latent space. Then the reduced states are fitted to a system of ordinary differential equations. The outputs from the model are then decoded back to the data space for a desired input parameter and time. The application of the PNODE formalism to different types of physical systems is important to test the methods robustness. The linear advection data was generated through a high-fidelity numerical tool for multiple velocity parameters. The PNODE code was modified for the advection dataset, whose temporal domain and spatial discretization varied from the original study configuration. The L2 norm between the reconstruction and surrogate model and the reconstruction plots were used to analyze the PNODE model performance. The model reconstructions presented mixed results. For a temporal domain of 20-time units, where multiple advection cycles were completed for each advection speed, the reconstructions did not agree with the surrogate model. For a reduced temporal domain of 5-time units, the reconstructions and surrogate models were in close agreement. Near the end of the temporal domain, deviations occurred likely resulting from the accumulation of numerical errors. Note, over the 5-time units, smaller advection speed parameters were unable to complete a cycle. The behavior for the 20-time units highlighted potential issues with imbalanced datasets and repeated features. The 5-time unit model illustrates PNODEs adaptability to this class of problems when the dataset is better posed.
ContributorsReithal, Richard Robert (Author) / Kim, Jeonglae (Thesis director) / Lee, Kookjin (Committee member) / Barrett, The Honors College (Contributor) / Mechanical and Aerospace Engineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor)
Created2022-12
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
The Tesla valve, originating from Nikola Tesla's "valvular conduit" patent in 1920, offers a unique solution to fluid control challenges by enabling unidirectional flow while impeding reverse flow. With applications ranging from fluid pumps to high-power engines, Tesla's design functions as a fluidic diode, inducing pressure drops across the valve to define its efficiency through diodicity. Through computational fluid dynamics (CFD) simulations using ANSYS Fluent, the impact of removing the bifurcated section on Tesla valve efficiency is explored. The T45-R, D-Valve, and GMF Valve designs are analyzed across a range of Reynolds numbers (Re). Results show that while the absence of bifurcation can lead to higher diodicity values due to increased flow divergence and vortex formation, efficiency varies depending on flow conditions. The T45-R valve exhibits linear diodicity increase with Reynolds number, plateauing at higher values due to reduced fluid inertia influence. Conversely, the D-Valve with bifurcation excels at lower Re values, while the non-bifurcated version proves more efficient at higher Re values. The GMF Valve with bifurcation demonstrates efficiency at lower Re values but decreases in effectiveness as Re rises, with the non-bifurcated version showing lower efficiency overall. Overall, this research provides insights into the fundamental physics and design considerations of Tesla valves, offering guidance for optimizing fluid control applications across diverse industries. The study underscores the importance of considering geometric variations and flow conditions when designing Tesla valves for specific applications, highlighting the intricate relationship between valve geometry, flow dynamics, and efficiency.
ContributorsWiley, Sean (Author) / Huang, Huei-Ping (Thesis director) / Kim, Jeonglae (Committee member) / Barrett, The Honors College (Contributor) / Mechanical and Aerospace Engineering Program (Contributor)
Created2024-05