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.

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.

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.

The project is mainly aimed at detecting the gas flow rate in Biosensors and medical health applications by means of an acoustic method using whistle based device. Considering the challenges involved in maintaining particular flow rate and back pressure for detecting certain analytes in breath analysis the proposed system along with a cell phone provides a suitable way to maintain the flow rate without any additional battery driven device. To achieve this, a system-level approach is implemented which involves development of a closed end whistle which is placed inside a tightly fitted constant back pressure tube. By means of experimentation pressure vs. flowrate curve is initially obtained and used for the development of the particular whistle. Finally, by means of an FFT code in a cell phone the flow rate vs. frequency characteristic curve is obtained. When a person respires through the device a whistle sound is generated which is captured by the cellphone microphone and a FFT analysis is performed to determine the frequency and hence the flow rate from the characteristic curve. This approach can be used to detect flow rate as low as low as 1L/min. The concept has been applied for the first time in this work to the development and optimization of a breath analyzer.

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.

Hydraulic fracturing is an effective technique used in well stimulation to increase petroleum well production. A combination of multi-stage hydraulic fracturing and horizontal drilling has led to the recent boom in shale gas production which has changed the energy landscape of North America.

During the fracking process, highly pressurized mixture of water and proppants (sand and chemicals) is injected into to a crack, which fractures the surrounding rock structure and proppants help in keeping the fracture open. Over a longer period, however, these fractures tend to close due to the difference between the compressive stress exerted by the reservoir on the fracture and the fluid pressure inside the fracture. During production, fluid pressure inside the fracture is reduced further which can accelerate the closure of a fracture.

In this thesis, we study the stress distribution around a hydraulic fracture caused by fluid production. It is shown that fluid flow can induce a very high hoop stress near the fracture tip. As the pressure gradient increases stress concentration increases. If a fracture is very thin, the flow induced stress along the fracture decreases, but the stress concentration at the fracture tip increases and become unbounded for an infinitely thin fracture.

The result from the present study can be used for studying the fracture closure problem, and ultimately this in turn can lead to the development of better proppants so that prolific well production can be sustained for a long period of time.

Multi-touch tablets and smart phones are now widely used in both workplace and consumer settings. Interacting with these devices requires hand and arm movements that are potentially complex and poorly understood. Experimental studies have revealed differences in performance that could potentially be associated with injury risk. However, underlying causes for performance differences are often difficult to identify. For example, many patterns of muscle activity can potentially result in similar behavioral output. Muscle activity is one factor contributing to forces in tissues that could contribute to injury. However, experimental measurements of muscle activity and force for humans are extremely challenging. Models of the musculoskeletal system can be used to make specific estimates of neuromuscular coordination and musculoskeletal forces. However, existing models cannot easily be used to describe complex, multi-finger gestures such as those used for multi-touch human computer interaction (HCI) tasks. We therefore seek to develop a dynamic musculoskeletal simulation capable of estimating internal musculoskeletal loading during multi-touch tasks involving multi digits of the hand, and use the simulation to better understand complex multi-touch and gestural movements, and potentially guide the design of technologies the reduce injury risk. To accomplish these, we focused on three specific tasks. First, we aimed at determining the optimal index finger muscle attachment points within the context of the established, validated OpenSim arm model using measured moment arm data taken from the literature. Second, we aimed at deriving moment arm values from experimentally-measured muscle attachments and using these values to determine muscle-tendon paths for both extrinsic and intrinsic muscles of middle, ring and little fingers. Finally, we aimed at exploring differences in hand muscle activation patterns during zooming and rotating tasks on the tablet computer in twelve subjects. Towards this end, our musculoskeletal hand model will help better address the neuromuscular coordination, safe gesture performance and internal loadings for multi-touch applications.

Multi-pulse particle tracking velocimetry (multi-pulse PTV) is a recently proposed flow measurement technique aiming to improve the performance of conventional PTV/ PIV. In this work, multi-pulse PTV is assessed based on PTV simulations in terms of spatial resolution, velocity measurement accuracy and the capability of acceleration measurement. The errors of locating particles, velocity measurement and acceleration measurement are analytically calculated and compared among quadruple-pulse, triple-pulse and dual-pulse PTV. The optimizations of triple-pulse and quadruple-pulse PTV are discussed, and criteria are developed to minimize the combined error in position, velocity and acceleration. Experimentally, the velocity and acceleration fields of a round impinging air jet are measured to test the triple-pulse technique. A high speed beam-splitting camera and a custom 8-pulsed laser system are utilized to achieve good timing flexibility and temporal resolution. A new method to correct the registration error between CCDs is also presented. Consequently, the velocity field shows good consistency between triple-pulse and dual-pulse measurements. The mean acceleration profile along the centerline of the jet is used as the ground truth for the verification of the triple-pulse PIV measurements of the acceleration fields. The instantaneous acceleration field of the jet is directly measured by triple-pulse PIV and presented. Accelerations up to 1,000 g's are measured in these experiments.

The role of environmental factors that influence atmospheric propagation of sound originating from freeway noise sources is studied with a combination of field experiments and numerical simulations. Acoustic propagation models are developed and adapted for refractive index depending upon meteorological conditions. A high-resolution multi-nested environmental forecasting model forced by coarse global analysis is applied to predict real meteorological profiles at fine scales. These profiles are then used as input for the acoustic models. Numerical methods for producing higher resolution acoustic refractive index fields are proposed. These include spatial and temporal nested meteorological simulations with vertical grid refinement. It is shown that vertical nesting can improve the prediction of finer structures in near-ground temperature and velocity profiles, such as morning temperature inversions and low level jet-like features. Accurate representation of these features is shown to be important for modeling sound refraction phenomena and for enabling accurate noise assessment. Comparisons are made using the acoustic model for predictions with profiles derived from meteorological simulations and from field experiment observations in Phoenix, Arizona. The challenges faced in simulating accurate meteorological profiles at high resolution for sound propagation applications are highlighted and areas for possible improvement are discussed.

A detailed evaluation of the environmental forecast is conducted by investigating the Surface Energy Balance (SEB) obtained from observations made with an eddy-covariance flux tower compared with SEB from simulations using several physical parameterizations of urban effects and planetary boundary layer schemes. Diurnal variation in SEB constituent fluxes are examined in relation to surface layer stability and modeled diagnostic variables. Improvement is found when adapting parameterizations for Phoenix with reduced errors in the SEB components. Finer model resolution (to 333 m) is seen to have insignificant ($<1\sigma$) influence on mean absolute percent difference of 30-minute diurnal mean SEB terms. A new method of representing inhomogeneous urban development density derived from observations of impervious surfaces with sub-grid scale resolution is then proposed for mesoscale applications. This method was implemented and evaluated within the environmental modeling framework. Finally, a new semi-implicit scheme based on Leapfrog and a fourth-order implicit time-filter is developed.

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.