Matching Items (14)
Filtering by
- All Subjects: Atmospheric science
- All Subjects: Physics
- Genre: Academic theses
- Creators: Huang, Huei-Ping
- Member of: Theses and Dissertations
Description
The evolution of single hairpin vortices and multiple interacting hairpin vortices are studied in direct numerical simulations of channel flow at Re-tau=395. The purpose of this study is to observe the effects of increased Reynolds number and varying initial conditions on the growth of hairpins and the conditions under which single hairpins autogenerate hairpin packets. The hairpin vortices are believed to provide a unified picture of wall turbulence and play an important role in the production of Reynolds shear stress which is directly related to turbulent drag. The structures of the initial three-dimensional vortices are extracted from the two-point spatial correlation of the fully turbulent direct numerical simulation of the velocity field by linear stochastic estimation and embedded in a mean flow having the profile of the fully turbulent flow. The Reynolds number of the present simulation is more than twice that of the Re-tau=180 flow from earlier literature and the conditional events used to define the stochastically estimated single vortex initial conditions include a number of new types of events such as quasi-streamwise vorticity and Q4 events. The effects of parameters like strength, asymmetry and position are evaluated and compared with existing results in the literature. This study then attempts to answer questions concerning how vortex mergers produce larger scale structures, a process that may contribute to the growth of length scale with increasing distance from the wall in turbulent wall flows. Multiple vortex interactions are studied in detail.
ContributorsParthasarathy, Praveen Kumar (Author) / Adrian, Ronald (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2011
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
Description
The partitioning of available solar energy into different fluxes at the Earth's surface is important in determining different physical processes, such as turbulent transport, subsurface hydrology, land-atmospheric interactions, etc. Direct measurements of these turbulent fluxes were carried out using eddy-covariance (EC) towers. However, the distribution of EC towers is sparse due to relatively high cost and practical difficulties in logistics and deployment. As a result, data is temporally and spatially limited and is inadequate to be used for researches at large scales, such as regional and global climate modeling. Besides field measurements, an alternative way is to estimate turbulent fluxes based on the intrinsic relations between surface energy budget components, largely through thermodynamic equilibrium. These relations, referred as relative efficiency, have been included in several models to estimate the magnitude of turbulent fluxes in surface energy budgets such as latent heat and sensible heat. In this study, three theoretical models based on the lumped heat transfer model, the linear stability analysis and the maximum entropy principle respectively, were investigated. Model predictions of relative efficiencies were compared with turbulent flux data over different land covers, viz. lake, grassland and suburban surfaces. Similar results were observed over lake and suburban surface but significant deviation is found over vegetation surface. The relative efficiency of outgoing longwave radiation is found to be orders of magnitude deviated from theoretic predictions. Meanwhile, results show that energy partitioning process is influenced by the surface water availability to a great extent. The study provides insight into what property is determining energy partitioning process over different land covers and gives suggestion for future models.
ContributorsYang, Jiachuan (Author) / Wang, Zhihua (Thesis advisor) / Huang, Huei-Ping (Committee member) / Vivoni, Enrique (Committee member) / Mays, Larry (Committee member) / Arizona State University (Publisher)
Created2012
Description
This thesis outlines the development of a vector retrieval technique, based on data assimilation, for a coherent Doppler LIDAR (Light Detection and Ranging). A detailed analysis of the Optimal Interpolation (OI) technique for vector retrieval is presented. Through several modifications to the OI technique, it is shown that the modified technique results in significant improvement in velocity retrieval accuracy. These modifications include changes to innovation covariance portioning, covariance binning, and analysis increment calculation. It is observed that the modified technique is able to make retrievals with better accuracy, preserves local information better, and compares well with tower measurements. In order to study the error of representativeness and vector retrieval error, a lidar simulator was constructed. Using the lidar simulator a thorough sensitivity analysis of the lidar measurement process and vector retrieval is carried out. The error of representativeness as a function of scales of motion and sensitivity of vector retrieval to look angle is quantified. Using the modified OI technique, study of nocturnal flow in Owens' Valley, CA was carried out to identify and understand uncharacteristic events on the night of March 27th 2006. Observations from 1030 UTC to 1230 UTC (0230 hr local time to 0430 hr local time) on March 27 2006 are presented. Lidar observations show complex and uncharacteristic flows such as sudden bursts of westerly cross-valley wind mixing with the dominant up-valley wind. Model results from Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS®) and other in-situ instrumentations are used to corroborate and complement these observations. The modified OI technique is used to identify uncharacteristic and extreme flow events at a wind development site. Estimates of turbulence and shear from this technique are compared to tower measurements. A formulation for equivalent wind speed in the presence of variations in wind speed and direction, combined with shear is developed and used to determine wind energy content in presence of turbulence.
ContributorsChoukulkar, Aditya (Author) / Calhoun, Ronald (Thesis advisor) / Mahalov, Alex (Committee member) / Kostelich, Eric (Committee member) / Huang, Huei-Ping (Committee member) / Phelan, Patrick (Committee member) / Arizona State University (Publisher)
Created2013
Description
The implications of a changing climate have a profound impact on human life, society, and policy making. The need for accurate climate prediction becomes increasingly important as we better understand these implications. Currently, the most widely used climate prediction relies on the synthesis of climate model simulations organized by the Coupled Model Intercomparison Project (CMIP); these simulations are ensemble-averaged to construct projections for the 21st century climate. However, a significant degree of bias and variability in the model simulations for the 20th century climate is well-known at both global and regional scales. Based on that insight, this study provides an alternative approach for constructing climate projections that incorporates knowledge of model bias. This approach is demonstrated to be a viable alternative which can be easily implemented by water resource managers for potentially more accurate projections. Tests of the new approach are provided on a global scale with an emphasis on semiarid regional studies for their particular vulnerability to water resource changes, using both the former CMIP Phase 3 (CMIP3) and current Phase 5 (CMIP5) model archives. This investigation is accompanied by a detailed analysis of the dynamical processes and water budget to understand the behaviors and sources of model biases. Sensitivity studies of selected CMIP5 models are also performed with an atmospheric component model by testing the relationship between climate change forcings and model simulated response. The information derived from each study is used to determine the progressive quality of coupled climate models in simulating the global water cycle by rigorously investigating sources of model bias related to the moisture budget. As such, the conclusions of this project are highly relevant to model development and potentially may be used to further improve climate projections.
ContributorsBaker, Noel C (Author) / Huang, Huei-Ping (Thesis advisor) / Trimble, Steve (Committee member) / Anderson, James (Committee member) / Clarke, Amanda (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2013
Description
This study considered the impact of grid resolution on wind velocity simulated by the Weather Research and Forecasting (WRF) model. The period simulated spanned November 2009 through January 2010, for which, multi-resolution nested domains were examined. Basic analysis was performed utilizing the data assimilation tools of NCEP/NCAR (National Center for Environmental Prediction/National Center for Atmospheric Research) to determine the ideal location to examine during the simulation was the Pacific Northwest portion of the United States, specifically the border between California and Oregon. The simulated mutli-resolution nested domains in this region indicated an increase in apparent wind speed as the resolution for the domain was increased. These findings were confirmed by statistical analysis which identified a positive bias for wind speed with respect to increased resolution as well as a correlation coefficient indicating the existence of a positive change in wind speed with increased resolution. An analysis of temperature change was performed in order to test the validity of the findings of the WRF simulation model. The statistical analysis performed on temperature change throughout the increased grid resolution did not indicate any change in temperature. In fact the correlation coefficient values between the domains were found in the 0.90 range, indicating the non-sensitivity of temperature across the increased resolutions. These results validate the findings of the WRF simulation: increased wind velocity can be observed at higher grid resolution. The study then considered the difference between wind velocity observed over the entire domains and the wind velocity observed solely over offshore locations. Wind velocity was observed to be significantly higher (an increase of 68.4%) in the offshore locations. The findings of this study suggest simulation tools should be utilized to examine domains at a higher resolution in order to identify potential locations for wind farms. The results go further to suggest the ideal location for these potential wind farms will be at offshore locations.
ContributorsBouey, Michael (Author) / Huang, Huei-Ping (Thesis advisor) / Trimble, Steve (Committee member) / Ronald, Ronald (Committee member) / Arizona State University (Publisher)
Created2012
Description
Rapid expansion of dense beds of fine, spherical particles subjected to rapid depressurization is studied in a vertical shock tube. As the particle bed is unloaded, a high-speed video camera captures the dramatic evolution of the particle bed structure. Pressure transducers are used to measure the dynamic pressure changes during the particle bed expansion process. Image processing, signal processing, and Particle Image Velocimetry techniques, are used to examine the relationships between particle size, initial bed height, bed expansion rate, and gas velocities.
The gas-particle interface and the particle bed as a whole expand and evolve in stages. First, the bed swells nearly homogeneously for a very brief period of time (< 2ms). Shortly afterward, the interface begins to develop instabilities as it continues to rise, with particles nearest the wall rising more quickly. Meanwhile, the bed fractures into layers and then breaks down further into cellular-like structures. The rate at which the structural evolution occurs is shown to be dependent on particle size. Additionally, the rate of the overall bed expansion is shown to be dependent on particle size and initial bed height.
Taller particle beds and beds composed of smaller-diameter particles are found to be associated with faster bed-expansion rates, as measured by the velocity of the gas-particle interface. However, the expansion wave travels more slowly through these same beds. It was also found that higher gas velocities above the the gas-particle interface measured \textit{via} Particle Image Velocimetry or PIV, were associated with particle beds composed of larger-diameter particles. The gas dilation between the shocktube diaphragm and the particle bed interface is more dramatic when the distance between the gas-particle interface and the diaphragm is decreased-as is the case for taller beds.
To further elucidate the complexities of this multiphase compressible flow, simple OpenFOAM (Weller, 1998) simulations of the shocktube experiment were performed and compared to bed expansion rates, pressure fluctuations, and gas velocities. In all cases, the trends and relationships between bed height, particle diameter, with expansion rates, pressure fluctuations and gas velocities matched well between experiments and simulations. In most cases, the experimentally-measured bed rise rates and the simulated bed rise rates matched reasonably well in early times. The trends and overall values of the pressure fluctuations and gas velocities matched well between the experiments and simulations; shedding light on the effects each parameter has on the overall flow.
The gas-particle interface and the particle bed as a whole expand and evolve in stages. First, the bed swells nearly homogeneously for a very brief period of time (< 2ms). Shortly afterward, the interface begins to develop instabilities as it continues to rise, with particles nearest the wall rising more quickly. Meanwhile, the bed fractures into layers and then breaks down further into cellular-like structures. The rate at which the structural evolution occurs is shown to be dependent on particle size. Additionally, the rate of the overall bed expansion is shown to be dependent on particle size and initial bed height.
Taller particle beds and beds composed of smaller-diameter particles are found to be associated with faster bed-expansion rates, as measured by the velocity of the gas-particle interface. However, the expansion wave travels more slowly through these same beds. It was also found that higher gas velocities above the the gas-particle interface measured \textit{via} Particle Image Velocimetry or PIV, were associated with particle beds composed of larger-diameter particles. The gas dilation between the shocktube diaphragm and the particle bed interface is more dramatic when the distance between the gas-particle interface and the diaphragm is decreased-as is the case for taller beds.
To further elucidate the complexities of this multiphase compressible flow, simple OpenFOAM (Weller, 1998) simulations of the shocktube experiment were performed and compared to bed expansion rates, pressure fluctuations, and gas velocities. In all cases, the trends and relationships between bed height, particle diameter, with expansion rates, pressure fluctuations and gas velocities matched well between experiments and simulations. In most cases, the experimentally-measured bed rise rates and the simulated bed rise rates matched reasonably well in early times. The trends and overall values of the pressure fluctuations and gas velocities matched well between the experiments and simulations; shedding light on the effects each parameter has on the overall flow.
ContributorsZunino, Heather (Author) / Adrian, Ronald J (Thesis advisor) / Clarke, Amanda (Committee member) / Chen, Kangping (Committee member) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
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
The study of deflagration to detonation transition (DDT) in explosives is of prime importance with regards to insensitive munitions (IM). Critical damage owing to thermal or shock stimuli could translate to significant loss of life and material. The present study models detonation and deflagration of a commonly used granular explosive: cyclotetramethylene-tetranitramine, HMX. A robust literature review is followed by computational modeling of gas gun and DDT tube test data using the Sandia National Lab three-dimensional multi-material Eulerian hydrocode CTH. This dissertation proposes new computational practices and models that aid in predicting shock stimulus IM response. CTH was first used to model experimental data sets of DDT tubes from both Naval Surface Weapons Center and Los Alamos National Laboratory which were initiated by pyrogenic material and a piston, respectively. Analytical verification was performed, where possible, for detonation via empirical based equations at the Chapman Jouguet state with errors below 2.1%, and deflagration via pressure dependent burn rate equations. CTH simulations include inert, history variable reactive burn and Arrhenius models. The results are in excellent agreement with published HMX detonation velocities. Novel additions include accurate simulation of the pyrogenic material BKNO3 and the inclusion of porosity in energetic materials. The treatment of compaction is especially important in modeling precursory hotspots, caused by hydrodynamic collapse of void regions or grain interactions, prior to DDT of granular explosives. The CTH compaction model of HMX was verified within 11% error via a five pronged validation approach using gas gun data and employed use of a newly generated set of P-α parameters for granular HMX in a Mie-Gruneisen Equation of State. Next, the additions of compaction were extended to a volumetric surface burning model of HMX and compare well to a set of empirical burn rates. Lastly, the compendium of detonation and deflagration models was applied to the aforementioned DDT tubes and demonstrate working functionalities of all models, albeit at the expense of significant computational resources. A robust hydrocode methodology is proposed to make use of the deflagration, compaction and detonation models as a means to predict IM response to shock stimulus of granular explosive materials.
ContributorsMahon, Kelly Susan (Author) / Lee, Taewoo (Thesis advisor) / Herrmann, Marcus (Committee member) / Chen, Kangping (Committee member) / Jiao, Yang (Committee member) / Huang, Huei-Ping (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