Matching Items (88)
Filtering by
- Creators: Huang, Huei-Ping
- Member of: ASU Electronic Theses and Dissertations
- Resource Type: Text
Description
A new method of adaptive mesh generation for the computation of fluid flows is investigated. The method utilizes gradients of the flow solution to adapt the size and stretching of elements or volumes in the computational mesh as is commonly done in the conventional Hessian approach. However, in the new method, higher-order gradients are used in place of the Hessian. The method is applied to the finite element solution of the incompressible Navier-Stokes equations on model problems. Results indicate that a significant efficiency benefit is realized.
ContributorsShortridge, Randall (Author) / Chen, Kang Ping (Thesis advisor) / Herrmann, Marcus (Thesis advisor) / Wells, Valana (Committee member) / Huang, Huei-Ping (Committee member) / Mittelmann, Hans (Committee member) / Arizona State University (Publisher)
Created2011
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
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 heat and mass transfer phenomena in micro-scale for the mass transfer phenomena on drug in cylindrical matrix system, the simulation of oxygen/drug diffusion in a three dimensional capillary network, and a reduced chemical kinetic modeling of gas turbine combustion for Jet propellant-10 have been studied numerically. For the numerical analysis of the mass transfer phenomena on drug in cylindrical matrix system, the governing equations are derived from the cylindrical matrix systems, Krogh cylinder model, which modeling system is comprised of a capillary to a surrounding cylinder tissue along with the arterial distance to veins. ADI (Alternative Direction Implicit) scheme and Thomas algorithm are applied to solve the nonlinear partial differential equations (PDEs). This study shows that the important factors which have an effect on the drug penetration depth to the tissue are the mass diffusivity and the consumption of relevant species during the time allowed for diffusion to the brain tissue. Also, a computational fluid dynamics (CFD) model has been developed to simulate the blood flow and oxygen/drug diffusion in a three dimensional capillary network, which are satisfied in the physiological range of a typical capillary. A three dimensional geometry has been constructed to replicate the one studied by Secomb et al. (2000), and the computational framework features a non-Newtonian viscosity model for blood, the oxygen transport model including in oxygen-hemoglobin dissociation and wall flux due to tissue absorption, as well as an ability to study the diffusion of drugs and other materials in the capillary streams. Finally, a chemical kinetic mechanism of JP-10 has been compiled and validated for a wide range of combustion regimes, covering pressures of 1atm to 40atm with temperature ranges of 1,200 K - 1,700 K, which is being studied as a possible Jet propellant for the Pulse Detonation Engine (PDE) and other high-speed flight applications such as hypersonic missiles. The comprehensive skeletal mechanism consists of 58 species and 315 reactions including in CPD, Benzene formation process by the theory for polycyclic aromatic hydrocarbons (PAH) and soot formation process on the constant volume combustor, premixed flame characteristics.
ContributorsBae, Kang-Sik (Author) / Lee, Taewoo (Thesis advisor) / Huang, Huei-Ping (Committee member) / Calhoun, Ronald (Committee member) / Phelan, Patrick (Committee member) / Lopez, Juan (Committee member) / Arizona State University (Publisher)
Created2012
Description
This doctoral thesis investigates the predictability characteristics of floods and flash floods by coupling high resolution precipitation products to a distributed hydrologic model. The research hypotheses are tested at multiple watersheds in the Colorado Front Range (CFR) undergoing warm-season precipitation. Rainfall error structures are expected to propagate into hydrologic simulations with added uncertainties by model parameters and initial conditions. Specifically, the following science questions are addressed: (1) What is the utility of Quantitative Precipitation Estimates (QPE) for high resolution hydrologic forecasts in mountain watersheds of the CFR?, (2) How does the rainfall-reflectivity relation determine the magnitude of errors when radar observations are used for flood forecasts?, and (3) What are the spatiotemporal limits of flood forecasting in mountain basins when radar nowcasts are used into a distributed hydrological model?. The methodology consists of QPE evaluations at the site (i.e., rain gauge location), basin-average and regional scales, and Quantitative Precipitation Forecasts (QPF) assessment through regional grid-to-grid verification techniques and ensemble basin-averaged time series. The corresponding hydrologic responses that include outlet discharges, distributed runoff maps, and streamflow time series at internal channel locations, are used in light of observed and/or reference data to diagnose the suitability of fusing precipitation forecasts into a distributed model operating at multiple catchments. Results reveal that radar and multisensor QPEs lead to an improved hydrologic performance compared to simulations driven with rain gauge data only. In addition, hydrologic performances attained by satellite products preserve the fundamental properties of basin responses, including a simple scaling relation between the relative spatial variability of runoff and its magnitude. Overall, the spatial variations contained in gridded QPEs add value for warm-season flood forecasting in mountain basins, with sparse data even if those products contain some biases. These results are encouraging and open new avenues for forecasting in regions with limited access and sparse observations. Regional comparisons of different reflectivity -rainfall (Z-R) relations during three summer seasons, illustrated significant rainfall variability across the region. Consistently, hydrologic errors introduced by the distinct Z-R relations, are significant and proportional (in the log-log space) to errors in precipitation estimations and stream flow magnitude. The use of operational Z-R relations without prior calibration may lead to wrong estimation of precipitation, runoff magnitude and increased flood forecasting errors. This suggests that site-specific Z-R relations, prior to forecasting procedures, are desirable in complex terrain regions. Nowcasting experiments show the limits of flood forecasting and its dependence functions of lead time and basin scale. Across the majority of the basins, flood forecasting skill decays with lead time, but the functional relation depends on the interactions between watershed properties and rainfall characteristics. Both precipitation and flood forecasting skills are noticeably reduced for lead times greater than 30 minutes. Scale dependence of hydrologic forecasting errors demonstrates reduced predictability at intermediate-size basins, the typical scale of convective storm systems. Overall, the fusion of high resolution radar nowcasts and the convenient parallel capabilities of the distributed hydrologic model provide an efficient framework for generating accurate real-time flood forecasts suitable for operational environments.
ContributorsMoreno Ramirez, Hernan (Author) / Vivoni, Enrique R. (Thesis advisor) / Ruddell, Benjamin L. (Committee member) / Gochis, David J. (Committee member) / Mays, Larry W. (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2012
Description
This study performs numerical modeling for the climate of semi-arid regions by running a high-resolution atmospheric model constrained by large-scale climatic boundary conditions, a practice commonly called climate downscaling. These investigations focus especially on precipitation and temperature, quantities that are critical to life in semi-arid regions. Using the Weather Research and Forecast (WRF) model, a non-hydrostatic geophysical fluid dynamical model with a full suite of physical parameterization, a series of numerical sensitivity experiments are conducted to test how the intensity and spatial/temporal distribution of precipitation change with grid resolution, time step size, the resolution of lower boundary topography and surface characteristics. Two regions, Arizona in U.S. and Aral Sea region in Central Asia, are chosen as the test-beds for the numerical experiments: The former for its complex terrain and the latter for the dramatic man-made changes in its lower boundary conditions (the shrinkage of Aral Sea). Sensitivity tests show that the parameterization schemes for rainfall are not resolution-independent, thus a refinement of resolution is no guarantee of a better result. But, simulations (at all resolutions) do capture the inter-annual variability of rainfall over Arizona. Nevertheless, temperature is simulated more accurately with refinement in resolution. Results show that both seasonal mean rainfall and frequency of extreme rainfall events increase with resolution. For Aral Sea, sensitivity tests indicate that while the shrinkage of Aral Sea has a dramatic impact on the precipitation over the confine of (former) Aral Sea itself, its effect on the precipitation over greater Central Asia is not necessarily greater than the inter-annual variability induced by the lateral boundary conditions in the model and large scale warming in the region. The numerical simulations in the study are cross validated with observations to address the realism of the regional climate model. The findings of this sensitivity study are useful for water resource management in semi-arid regions. Such high spatio-temporal resolution gridded-data can be used as an input for hydrological models for regions such as Arizona with complex terrain and sparse observations. Results from simulations of Aral Sea region are expected to contribute to ecosystems management for Central Asia.
ContributorsSharma, Ashish (Author) / Huang, Huei-Ping (Thesis advisor) / Adrian, Ronald (Committee member) / Herrmann, Marcus (Committee member) / Phelan, Patrick E. (Committee member) / Vivoni, Enrique (Committee member) / Arizona State University (Publisher)
Created2012
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
Image processing in canals, rivers and other bodies of water has been a very important concern. This research using Image Processing was performed to obtain a photographic evidence of the data of the site which helps in monitoring the conditions of the water body and the surroundings. Images are captured using a digital camera and the images are stored onto a datalogger, these images are retrieved using a cellular/ satellite modem. A MATLAB program was designed to obtain the level of water by just entering the file name into to the program, a curve fit model was created to determine the contrast parameters. The contrast parameters were obtained using the data obtained from the gray scale image mainly the mean and variance of the intensity values. The enhanced images are used to determine the level of water by taking pixel intensity plots along the region of interest. The level of water obtained is accurate to less than 2% of the actual level of water observed from the image. High speed imaging in micro channels have various application in industrial field, medical field etc. In medical field it is tested by using blood samples. The experimental procedure proposed determines the flow duration and the defects observed in these channel using a fluid introduced into the micro channel the fluid being water based dye and whole milk. The viscosity of the fluid shows different types of flow patterns and defects in the micro channel. The defects observed vary from a small effect to the flow pattern to an extreme defect in the channel such as obstruction of flow or deformation in the channel. The sample needs to be further analyzed by SEM to get a better insight on the defects.
ContributorsShasedhara, Abhijeet Bangalore (Author) / Lee, Taewoo (Thesis advisor) / Huang, Huei-Ping (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2011
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 thesis focuses on studying the interaction between floating objects and an air-water flow system driven by gravity. The system consists of an inclined channel in which a gravity driven two phase flow carries a series of floating solid objects downstream. Numerical simulations of such a system requires the solution of not only the basic Navier-Stokes equation but also dynamic interaction between the solid body and the two-phase flow. In particular, this requires embedding of dynamic mesh within the two-phase flow. A computational fluid dynamics solver, ANSYS fluent, is used to solve this problem. Also, the individual components for these simulations are already available in the solver, few examples exist in which all are combined. A series of simulations are performed by varying the key parameters, including density of floating objects and mass flow rate at the inlet. The motion of the floating objects in those simulations are analyzed to determine the stability of the coupled flow-solid system. The simulations are successfully performed over a broad range of parametric values. The numerical framework developed in this study can potentially be used in applications, especially in assisting the design of similar gravity driven systems for transportation in manufacturing processes. In a small number of the simulations, two kinds of numerically instability are observed. One is characterized by a sudden vertical acceleration of the floating object due to a strong imbalance of the force acting on the body, which occurs when the mass flow of water is weak. The other is characterized by a sudden vertical movement of air-water interface, which occurs when two floating objects become too close together. These new types of numerical instability deserve future studies and clarifications. This study is performed only for a 2-D system. Extension of the numerical framework to a full 3-D setting is recommended as future work.
ContributorsMangavelli, Sai Chaitanya (Author) / Huang, Huei-Ping (Thesis advisor) / Kim, Jeonglae (Committee member) / Forzani, Erica (Committee member) / Arizona State University (Publisher)
Created2018