Matching Items (7)
Filtering by
- All Subjects: Spectral Element Method
- All Subjects: Wind power
- All Subjects: stereolithography
- Creators: Peet, Yulia
- Resource Type: Text
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
Description
With a ground-based Doppler lidar on the upwind side of a wind farm in the Tehachapi Pass of California, radial wind velocity measurements were collected for repeating sector sweeps, scanning up to 10 kilometers away. This region consisted of complex terrain, with the scans made between mountains. The dataset was utilized for techniques being studied for short-term forecasting of wind power by correlating changes in energy content and of turbulence intensity by tracking spatial variance, in the wind ahead of a wind farm. A ramp event was also captured and its propagation was tracked.
Orthogonal horizontal wind vectors were retrieved from the radial velocity using a sector Velocity Azimuth Display method. Streamlines were plotted to determine the potential sites for a correlation of upstream wind speed with wind speed at downstream locations near the wind farm. A "virtual wind turbine" was "placed" in locations along the streamline by using the time-series velocity data at the location as the input to a modeled wind turbine, to determine the extractable energy content at that location. The relationship between this time-dependent energy content upstream and near the wind farm was studied. By correlating the energy content with each upstream location based on a time shift estimated according to advection at the mean wind speed, several fits were evaluated. A prediction of the downstream energy content was produced by shifting the power output in time and applying the best-fit function. This method made predictions of the power near the wind farm several minutes in advance. Predictions were also made up to an hour in advance for a large ramp event. The Magnitude Absolute Error and Standard Deviation are presented for the predictions based on each selected upstream location.
Orthogonal horizontal wind vectors were retrieved from the radial velocity using a sector Velocity Azimuth Display method. Streamlines were plotted to determine the potential sites for a correlation of upstream wind speed with wind speed at downstream locations near the wind farm. A "virtual wind turbine" was "placed" in locations along the streamline by using the time-series velocity data at the location as the input to a modeled wind turbine, to determine the extractable energy content at that location. The relationship between this time-dependent energy content upstream and near the wind farm was studied. By correlating the energy content with each upstream location based on a time shift estimated according to advection at the mean wind speed, several fits were evaluated. A prediction of the downstream energy content was produced by shifting the power output in time and applying the best-fit function. This method made predictions of the power near the wind farm several minutes in advance. Predictions were also made up to an hour in advance for a large ramp event. The Magnitude Absolute Error and Standard Deviation are presented for the predictions based on each selected upstream location.
ContributorsMagerman, Beth (Author) / Calhoun, Ronald (Thesis advisor) / Peet, Yulia (Committee member) / Huang, Huei-Ping (Committee member) / Krishnamurthy, Raghavendra (Committee member) / Arizona State University (Publisher)
Created2014
Numerical simulation of environmental flow over urban landscape for applications to renewable energy
Description
Development of renewable energy solutions has become a major interest among environmental organizations and governments around the world due to an increase in energy consumption and global warming. One fast growing renewable energy solution is the application of wind energy in cities. To qualitative and quantitative predict wind turbine performance in urban areas, CFD simulation is performed on real-life urban geometry and wind velocity profiles are evaluated. Two geometries in Arizona is selected in this thesis to demonstrate the influence of building heights; one of the simulation models, ASU campus, is relatively low rise and without significant tall buildings; the other model, the downtown phoenix model, are high-rise and with greater building height difference. The content of this thesis focuses on using RANS computational fluid dynamics approach to simulate wind acceleration phenomenon in two complex geometries, ASU campus and Phoenix downtown model. Additionally, acceleration ratio and locations are predicted, the results are then used to calculate the best location for small wind turbine installments.
ContributorsYing, Xiaoyan (Author) / Huang, Huei-Ping (Thesis advisor) / Peet, Yulia (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2015
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 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.
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
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 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.
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
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