Matching Items (96)
Filtering by
- Creators: Huang, Huei-Ping
- Member of: ASU Electronic Theses and Dissertations
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
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
For the increasing concerns of influence on environment by fossil-electricity generation, application of renewable energy becomes one of the most focused issues in society. Based on the limitation on urban environment, wind turbines, which can be mounted on rooftop or between buildings, are regarded as a feasible way for wind energy generation. This study presents wind flow simulations in a large-scale environment with certain dimension buildings. Different inlet velocity boundary conditions are tested firstly, and the non-uniform inlet boundary condition shows better agreement with realistic situation. Turbulence intensity is set to be 10% for comparison consistency. The k-epsilon turbulence model is regarded as a better simulation for this certain condition. After that, three different structures, which include single building, pristine double building and modified circular gap double building systems, are tested in this environment condition. The result shows 18.8% velocity increasing on the top of single building system. Pristine double building systems are tested with 4 different gap distances, and building with 10 meters gap achieved the best velocity condition, which 32.8% velocity increasing and 11.8% improvement comparing to single building system, respectively. But the location of maximum velocity moves to the gap and the maximum velocity on the rooftop of double building system is approximately 5.1% lower than single building system. Based on previous study, modified circular gap double building system is created with 10 meters gap. Comparing result with single building system, modified circular gap system achieves higher improvement for wind flow, whose improvement of velocity increasing in the gap and on the rooftop of building are 47.1% and 3.0%, respectively. As a result, the modified circular gap double building can be regarded as a high efficiency system of environmental wind flow over buildings for renewable energy system.
ContributorsLi, Guoyi (Author) / Huang, Huei-Ping (Thesis advisor) / Lee, Taewoo (Committee member) / Forzani, Erica (Committee member) / Arizona State University (Publisher)
Created2015
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
First, in a large-scale structure, a 3-D CFD model was built to simulate flow and temperature distributions. The flow patterns and temperature distributions are characterized and validated through spot measurements. The detailed understanding of them then allows for optimization of the HVAC configuration because identification of the problematic flow patterns and temperature mis-distributions leads to some corrective measures. Second, an appropriate form of the viscous dissipation term in the integral form of the conservation equation was considered, and the effects of momentum terms on the computed drop size in pressure-atomized sprays were examined. The Sauter mean diameter (SMD) calculated in this manner agrees well with experimental data of the drop velocities and sizes. Using the suggested equation with the revised treatment of liquid momentum setup, injection parameters can be directly input to the system of equations. Thus, this approach is capable of incorporating the effects of injection parameters for further considerations of the drop and velocity distributions under a wide range of spray geometry and injection conditions. Lastly, groundwater level estimation was investigated using compressed sensing (CS). To satisfy a general property of CS, a random measurement matrix was used, the groundwater network was constructed, and finally the l-1 optimization was run. Through several validation tests, correct estimation of groundwater level by CS was shown. Using this setup, decreasing trends in groundwater level in the southwestern US was shown. The suggested method is effective in that the total measurements of registered wells can be reduced down by approximately 42 %, sparse data can be visualized and a possible approach for groundwater management during extreme weather changes, e.g. in California, was demonstrated.
ContributorsLee, Joon Young (Author) / Lee, Taewoo (Thesis advisor) / Huang, Huei-Ping (Committee member) / Lopez, Juan (Committee member) / Phelan, Patrick (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
Description
The subject of this thesis is concerned with the amount of cooling air assigned to seal high pressure turbine rim cavities which is critical for performance as well as component life. Insufficient air leads to excessive hot annulus gas ingestion and its penetration deep into the cavity compromising disc life. Excessive purge air, adversely affects performance. Experiments on a rotating turbine stage rig which included a rotor-stator forward disc cavity were performed at Arizona State University. The turbine rig has 22 vanes and 28 blades, while the rim cavity is composed of a single-tooth rim lab seal and a rim platform overlap seal. Time-averaged static pressures were measured in the gas path and the cavity, while mainstream gas ingestion into the cavity was determined by measuring the concentration distribution of tracer gas (carbon dioxide). Additionally, particle image velocimetry (PIV) was used to measure fluid velocity inside the rim cavity between the lab seal and the overlap. The data from the experiments were compared to an 360-degree unsteady RANS (URANS) CFD simulations. Although not able to match the time-averaged test data satisfactorily, the CFD simulations brought to light the unsteadiness present in the flow during the experiment which the slower response data did not fully capture. To interrogate the validity of URANS simulations in capturing complex rotating flow physics, the scope of this work also included to validating the CFD tool by comparing its predictions against experimental LDV data in a closed rotor-stator cavity. The enclosed cavity has a stationary shroud, a rotating hub, and mass flow does not enter or exit the system. A full 360 degree numerical simulation was performed comparing Fluent LES, with URANS turbulence models. Results from these investigations point to URANS state of art under-predicting closed cavity tangential velocity by 32% to 43%, and open rim cavity effectiveness by 50% compared to test data. The goal of this thesis is to assess the validity of URANS turbulence models in more complex rotating flows, compare accuracy with LES simulations, suggest CFD settings to better simulate turbine stage mainstream/disc cavity interaction with ingestion, and recommend experimentation techniques.
ContributorsKanjiyani, Shezan (Author) / Lee, Taewoo (Thesis advisor) / Mirzamoghadam, Alexander (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2016
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
A numerical study of incremental spin-up and spin-up from rest of a thermally- stratified fluid enclosed within a right circular cylinder with rigid bottom and side walls and stress-free upper surface is presented. Thermally stratified spin-up is a typical example of baroclinity, which is initiated by a sudden increase in rotation rate and the tilting of isotherms gives rise to baroclinic source of vorticity. Research by (Smirnov et al. [2010a]) showed the differences in evolution of instabilities when Dirichlet and Neumann thermal boundary conditions were applied at top and bottom walls. Study of parametric variations carried out in this dissertation confirmed the instability patterns observed by them for given aspect ratio and Rossby number values greater than 0.5. Also results reveal that flow maintained axisymmetry and stability for short aspect ratio containers independent of amount of rotational increment imparted. Investigation on vorticity components provides framework for baroclinic vorticity feedback mechanism which plays important role in delayed rise of instabilities when Dirichlet thermal Boundary Conditions are applied.
ContributorsKher, Aditya Deepak (Author) / Chen, Kangping (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2011
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
Description
With the ever-increasing demand for high-end services, technological companies have been forced to operate on high performance servers. In addition to the customer services, the company's internal need to store and manage huge amounts of data has also increased their need to invest in High Density Data Centers. As a result, the performance to size of the data center has increased tremendously. Most of the consumed power by the servers is emitted as heat. In a High Density Data Center, the power per floor space area is higher compared to the regular data center. Hence the thermal management of this type of data center is relatively complicated.
Because of the very high power emission in a smaller containment, improper maintenance can result in failure of the data center operation in a shorter period. Hence the response time of the cooler to the temperature rise of the servers is very critical. Any delay in response will constantly lead to increased temperature and hence the server's failure.
In this paper, the significance of this delay time is understood by performing CFD simulation on different variants of High Density Modules using ANSYS Fluent. It was found out that the delay was becoming longer as the size of the data center increases. But the overload temperature, ie. the temperature rise beyond the set-point became lower with the increase in data center size. The results were common for both the single-row and the double-row model. The causes of the increased delay are accounted and explained in detail manner in this paper.
Because of the very high power emission in a smaller containment, improper maintenance can result in failure of the data center operation in a shorter period. Hence the response time of the cooler to the temperature rise of the servers is very critical. Any delay in response will constantly lead to increased temperature and hence the server's failure.
In this paper, the significance of this delay time is understood by performing CFD simulation on different variants of High Density Modules using ANSYS Fluent. It was found out that the delay was becoming longer as the size of the data center increases. But the overload temperature, ie. the temperature rise beyond the set-point became lower with the increase in data center size. The results were common for both the single-row and the double-row model. The causes of the increased delay are accounted and explained in detail manner in this paper.
ContributorsRamaraj, Dinesh Balaji (Author) / Gupta, Sandeep (Thesis advisor) / Hermann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2015