Matching Items (8)
Filtering by
- All Subjects: Mathematics
- Status: Published
Description
The numerical climate models have provided scientists, policy makers and the general public, crucial information for climate projections since mid-20th century. An international effort to compare and validate the simulations of all major climate models is organized by the Coupled Model Intercomparison Project (CMIP), which has gone through several phases since 1995 with CMIP5 being the state of the art. In parallel, an organized effort to consolidate all observational data in the past century culminates in the creation of several "reanalysis" datasets that are considered the closest representation of the true observation. This study compared the climate variability and trend in the climate model simulations and observations on the timescales ranging from interannual to centennial. The analysis focused on the dynamic climate quantity of zonal-mean zonal wind and global atmospheric angular momentum (AAM), and incorporated multiple datasets from reanalysis and the most recent CMIP3 and CMIP5 archives. For the observation, the validation of AAM by the length-of-day (LOD) and the intercomparison of AAM revealed a good agreement among reanalyses on the interannual and the decadal-to-interdecadal timescales, respectively. But the most significant discrepancies among them are in the long-term mean and long-term trend. For the simulations, the CMIP5 models produced a significantly smaller bias and a narrower ensemble spread of the climatology and trend in the 20th century for AAM compared to CMIP3, while CMIP3 and CMIP5 simulations consistently produced a positive trend for the 20th and 21st century. Both CMIP3 and CMIP5 models produced a wide range of the magnitudes of decadal and interdecadal variability of wind component of AAM (MR) compared to observation. The ensemble means of CMIP3 and CMIP5 are not statistically distinguishable for either the 20th- or 21st-century runs. The in-house atmospheric general circulation model (AGCM) simulations forced by the sea surface temperature (SST) taken from the CMIP5 simulations as lower boundary conditions were carried out. The zonal wind and MR in the CMIP5 simulations are well simulated in the AGCM simulations. This confirmed SST as an important mediator in regulating the global atmospheric changes due to GHG effect.
ContributorsPaek, Houk (Author) / Huang, Huei-Ping (Thesis advisor) / Adrian, Ronald (Committee member) / Wang, Zhihua (Committee member) / Anderson, James (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2013
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
The dynamics of a fluid flow inside 2D square and 3D cubic cavities
under various configurations were simulated and analyzed using a
spectral code I developed.
This code was validated against known studies in the 3D lid-driven
cavity. It was then used to explore the various dynamical behaviors
close to the onset of instability of the steady-state flow, and explain
in the process the mechanism underlying an intermittent bursting
previously observed. A fairly complete bifurcation picture emerged,
using a combination of computational tools such as selective
frequency damping, edge-state tracking and subspace restriction.
The code was then used to investigate the flow in a 2D square cavity
under stable temperature stratification, an idealized version of a lake
with warmer water at the surface compared to the bottom. The governing
equations are the Navier-Stokes equations under the Boussinesq approximation.
Simulations were done over a wide range of parameters of the problem quantifying
the driving velocity at the top (e.g. wind) and the strength of the stratification.
Particular attention was paid to the mechanisms associated with the onset of
instability of the base steady state, and the complex nontrivial dynamics
occurring beyond onset, where the presence of multiple states leads to a
rich spectrum of states, including homoclinic and heteroclinic chaos.
A third configuration investigates the flow dynamics of a fluid in a rapidly
rotating cube subjected to small amplitude modulations. The responses were
quantified by the global helicity and energy measures, and various peak
responses associated to resonances with intrinsic eigenmodes of the cavity
and/or internal retracing beams were clearly identified for the first time.
A novel approach to compute the eigenmodes is also described, making accessible
a whole catalog of these with various properties and dynamics. When the small
amplitude modulation does not align with the rotation axis (precession) we show
that a new set of eigenmodes are primarily excited as the angular velocity
increases, while triadic resonances may occur once the nonlinear regime kicks in.
under various configurations were simulated and analyzed using a
spectral code I developed.
This code was validated against known studies in the 3D lid-driven
cavity. It was then used to explore the various dynamical behaviors
close to the onset of instability of the steady-state flow, and explain
in the process the mechanism underlying an intermittent bursting
previously observed. A fairly complete bifurcation picture emerged,
using a combination of computational tools such as selective
frequency damping, edge-state tracking and subspace restriction.
The code was then used to investigate the flow in a 2D square cavity
under stable temperature stratification, an idealized version of a lake
with warmer water at the surface compared to the bottom. The governing
equations are the Navier-Stokes equations under the Boussinesq approximation.
Simulations were done over a wide range of parameters of the problem quantifying
the driving velocity at the top (e.g. wind) and the strength of the stratification.
Particular attention was paid to the mechanisms associated with the onset of
instability of the base steady state, and the complex nontrivial dynamics
occurring beyond onset, where the presence of multiple states leads to a
rich spectrum of states, including homoclinic and heteroclinic chaos.
A third configuration investigates the flow dynamics of a fluid in a rapidly
rotating cube subjected to small amplitude modulations. The responses were
quantified by the global helicity and energy measures, and various peak
responses associated to resonances with intrinsic eigenmodes of the cavity
and/or internal retracing beams were clearly identified for the first time.
A novel approach to compute the eigenmodes is also described, making accessible
a whole catalog of these with various properties and dynamics. When the small
amplitude modulation does not align with the rotation axis (precession) we show
that a new set of eigenmodes are primarily excited as the angular velocity
increases, while triadic resonances may occur once the nonlinear regime kicks in.
ContributorsWu, Ke (Author) / Lopez, Juan (Thesis advisor) / Welfert, Bruno (Thesis advisor) / Tang, Wenbo (Committee member) / Platte, Rodrigo (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2019
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 formation and stability of a slowly evolving zonal jet in 2-D flow with beta effect is analyzed using the model developed by Manfroi and Young in which the final governing equation was derived by means of a perturbation analysis of a barotropic vorticity equation with sinusoidal meridional mean flow. However in the original study the term β0, that represents the effect of large-scale Rossby waves, was dropped and was proceeded on a path of finding solutions for a simplified 1-D flow. The idea of this study is to understand the effects of the dropped term on the overall dynamics of the zonal jet evolution. For this purpose the system that is entirely deterministic with no additional forcing is solved by means of a standard finite difference scheme. The Numerical solutions are found for varying β0 and μ values where μ represents the bottom drag. In addition to this the criteria for the formation of zonal jets developed originally for the 1-D system is verified for the 2-D system as well. The study reveals the similarity in some of the results of the 1-D and the 2-D system like the merging of jets in the absence of bottom drag, formation of steady jets in presence of a non-zero bottom drag and the adherence to the boundary criteria for the formation of zonal jets. But when it comes to the formation of steady jets, a finite β0 value is required above which the solution is similar to the 1-D system. Also the jets formed under the presence of non-zero bottom drag seem wavy in nature which is different from the steady horizontal jets produced in the 1-D system.
ContributorsRaghunathan, Girish Nigamanth (Author) / Huang, Huei-Ping (Thesis advisor) / Herrmann, Marcus (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2017
Description
This study uses the Weather Research and Forecasting (WRF) model to simulate and predict the changes in local climate attributed to the urbanization for five desert cities. The simulations are performed in the fashion of climate downscaling, constrained by the surface boundary conditions generated from high resolution land-use maps. For each city, the land-use maps of 1985 and 2010 from Landsat satellite observation, and a projected land-use map for 2030, are used to represent the past, present, and future. An additional set of simulations for Las Vegas, the largest of the five cities, uses the NLCD 1992 and 2006 land-use maps and an idealized historical land-use map with no urban coverage for 1900.
The study finds that urbanization in Las Vegas produces a classic urban heat island (UHI) at night but a minor cooling during the day. A further analysis of the surface energy balance shows that the decrease in surface Albedo and increase effective emissivity play an important role in shaping the local climate change over urban areas. The emerging urban structures slow down the diurnal wind circulation over the city due to an increased effective surface roughness. This leads to a secondary modification of temperature due to the interaction between the mechanical and thermodynamic effects of urbanization.
The simulations for the five desert cities for 1985 and 2010 further confirm a common pattern of the climatic effect of urbanization with significant nighttime warming and moderate daytime cooling. This effect is confined to the urban area and is not sensitive to the size of the city or the detail of land cover in the surrounding areas. The pattern of nighttime warming and daytime cooling remains robust in the simulations for the future climate of the five cities using the projected 2030 land-use maps. Inter-city differences among the five urban areas are discussed.
The study finds that urbanization in Las Vegas produces a classic urban heat island (UHI) at night but a minor cooling during the day. A further analysis of the surface energy balance shows that the decrease in surface Albedo and increase effective emissivity play an important role in shaping the local climate change over urban areas. The emerging urban structures slow down the diurnal wind circulation over the city due to an increased effective surface roughness. This leads to a secondary modification of temperature due to the interaction between the mechanical and thermodynamic effects of urbanization.
The simulations for the five desert cities for 1985 and 2010 further confirm a common pattern of the climatic effect of urbanization with significant nighttime warming and moderate daytime cooling. This effect is confined to the urban area and is not sensitive to the size of the city or the detail of land cover in the surrounding areas. The pattern of nighttime warming and daytime cooling remains robust in the simulations for the future climate of the five cities using the projected 2030 land-use maps. Inter-city differences among the five urban areas are discussed.
ContributorsKamal, Samy (Author) / Huang, Huei-Ping (Thesis advisor) / Anderson, James (Thesis advisor) / Herrmann, Marcus (Committee member) / Calhoun, Ronald (Committee member) / Myint, Soe (Committee member) / Arizona State University (Publisher)
Created2015
Description
In this study, the stereolithography (SLA) 3D printing method is used to manufacture honeycomb-shaped flat sorbents that can capture CO2 from the air. The 3D-printed sorbents were synthesized using polyvinyl alcohol (PVA), propylene glycol, photopolymer resin, and an ion exchange resin (IER). The one-factor-at-a-time (OFAT) design-of-experiment approach was employed to determine the best combination ratio of materials to achieve high moisture swing and a good turnout of printed sorbents. The maximum load limit of the liquid photopolymer resin to enable printability of sorbents was found to be 44%. A series of moisture swing experiments was conducted to investigate the adsorption and desorption performance of the 3D-printed sorbents and compare them with the performance of IER samples prepared by a conventional approach. Results from these experiments conducted indicate that the printed sorbents showed less CO2 adsorptive characteristics compared to the conventional IER sample. It is proposed for future research that a liquid photopolymer resin made up of an IER be synthesized in order to improve the CO2-capturing ability of manufactured sorbents.
ContributorsObeng-Ampomah, Terry (Author) / Phelan, Patrick (Thesis advisor) / Lackner, Klaus (Committee member) / Shuaib, Abdelrahman (Committee member) / Arizona State University (Publisher)
Created2020
Description
This honors thesis explores and models the flow of air around a cylindrical arrow that is rotating as it moves through the air. This model represents the airflow around an archery arrow after it is released from the bow and rotates while it flies through the air. This situation is important in archery because an understanding of the airflow allows archers to predict the flight of the arrow. As a result, archers can improve their accuracy and ability to hit targets. However, not many computational fluid dynamic simulations modeling the airflow around a rotating archery arrow exist. This thesis attempts to further the understanding of the airflow around a rotating archery arrow by creating a mathematical model to numerically simulate the airflow around the arrow in the presence of this rotation. This thesis uses a linearized approximation of the Navier Stokes equations to model the airflow around the arrow and explains the reasoning for using this simplification of the fully nonlinear Navier Stokes equations. This thesis continues to describe the discretization of these linearized equations using the finite difference method and the boundary conditions used for these equations. A MATLAB code solves the resulting system of equations in order to obtain a numerical simulation of this airflow around the rotating arrow. The results of the simulation for each velocity component and the pressure distribution are displayed. This thesis then discusses the results of the simulation, and the MATLAB code is analyzed to verify the convergence of the solution. Appendix A includes the full MATLAB code used for the flow simulation. Finally, this thesis explains potential future research topics, ideas, and improvements to the code that can help further the understanding and create more realistic simulations of the airflow around a flying archery arrow.
ContributorsCholinski, Christopher John (Author) / Tang, Wenbo (Thesis director) / Herrmann, Marcus (Committee member) / Mechanical and Aerospace Engineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2019-05