Many physical phenomena and industrial applications involve multiphase fluid flows and hence it is of high importance to be able to simulate various aspects of these flows accurately. The Dynamic Contact Angles (DCA) and the contact lines at the wall boundaries are a couple of such important aspects. In the past few decades, many mathematical models were developed for predicting the contact angles of the inter-face with the wall boundary under various flow conditions. These models are used to incorporate the physics of DCA and contact line motion in numerical simulations using various interface capturing/tracking techniques. In the current thesis, a simple approach to incorporate the static and dynamic contact angle boundary conditions using the level set method is developed and implemented in multiphase CFD codes, LIT (Level set Interface Tracking) (Herrmann (2008)) and NGA (flow solver) (Desjardins et al (2008)). Various DCA models and associated boundary conditions are reviewed. In addition, numerical aspects such as the occurrence of a stress singularity at the contact lines and grid convergence of macroscopic interface shape are dealt with in the context of the level set approach.

This dissertation describes a process for interface capturing via an arbitrary-order, nearly quadrature free, discontinuous Galerkin (DG) scheme for the conservative level set method (Olsson et al., 2005, 2008). The DG numerical method is utilized to solve both advection and reinitialization, and executed on a refined level set grid (Herrmann, 2008) for effective use of processing power. Computation is executed in parallel utilizing both CPU and GPU architectures to make the method feasible at high order. Finally, a sparse data structure is implemented to take full advantage of parallelism on the GPU, where performance relies on well-managed memory operations.

With solution variables projected into a kth order polynomial basis, a k+1 order convergence rate is found for both advection and reinitialization tests using the method of manufactured solutions. Other standard test cases, such as Zalesak's disk and deformation of columns and spheres in periodic vortices are also performed, showing several orders of magnitude improvement over traditional WENO level set methods. These tests also show the impact of reinitialization, which often increases shape and volume errors as a result of level set scalar trapping by normal vectors calculated from the local level set field.

Accelerating advection via GPU hardware is found to provide a 30x speedup factor comparing a 2.0GHz Intel Xeon E5-2620 CPU in serial vs. a Nvidia Tesla K20 GPU, with speedup factors increasing with polynomial degree until shared memory is filled. A similar algorithm is implemented for reinitialization, which relies on heavier use of shared and global memory and as a result fills them more quickly and produces smaller speedups of 18x.

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.

Conventional fluid dynamics models such as the Navier-Stokes equations are derived for prediction of fluid motion at or near equilibrium, classic examples being the motion of fluids for which inter-molecular collisions are dominant. Flows at equilibrium permit simplifications such as the introduction of viscosity and also lead to solutions that are single-valued. However, many other regimes of interest include "fluids"' far from equilibrium; for example, rarefied gases or particle-laden flows in which the dispersed phase can be comprised of granular solids, droplets, or bubbles. Particle motion in these flows is not typically dominated by collisions and may exhibit significant memory effects; therefore, is often poorly described using continuum, field-based (Eulerian) approaches. Non-equilibrium flows generally lack a straightforward counterpart to viscosity and their multi-valued solutions cannot be represented by most Eulerian methods. This strongly motivates different strategies to address current shortcomings and the novel approach adopted in this work is based on the Conditional Quadrature Method of Moments (CQMOM). In CQMOM, moment equations are derived from the Boltzmann equation using a quadrature approximation of the velocity probability density function (PDF). CQMOM circumvents the drawbacks of current methods and leads to multivariate and multidimensional solutions in an Eulerian frame of reference. In the present work, the discretized PDF is resolved using an adaptive two-point quadrature in three-dimensional velocity space. The method is applied to computation of a series of non-equilibrium flows, ranging from simple two-dimensional test cases to fully-turbulent three-dimensional wall-bounded particle-laden flows. The primary contribution of the present effort is on development, application, and assessment of CQMOM for predicting the key features of dilute particle-laden flows. Statistical descriptors such as mean concentration and mean velocity are in good agreement with previous results, for both collision-less and collisional flows at varying particle Stokes numbers. Turbulent statistics and measures of local accumulation agree less favorably with prior results and identify areas for improvement in the modeling strategy.

The Volume-of-Fluid method is a popular method for interface tracking in Multiphase applications within Computational Fluid Dynamics. To date there exists several algorithms for reconstruction of a geometric interface surface. Of these are the Finite Difference algorithm, Least Squares Volume-of-Fluid Interface Reconstruction Algorithm, LVIRA, and the Efficient Least Squares Volume-of-Fluid Interface Reconstruction Algorithm, ELVIRA. Along with these geometric interface reconstruction algorithms, there exist several volume-of-fluid transportation algorithms. This paper will discuss two operator-splitting advection algorithms and an unsplit advection algorithm. Using these three interface reconstruction algorithms, and three advection algorithms, a comparison will be drawn to see how different combinations of these algorithms perform with respect to accuracy as well as computational expense.

Cavitation erosion is a significant cause of wear in marine components, such as impellers, propellers or rudders. While the erosion process has been widely studied on metals, the effect of cavitation on polymers is not well-understood. The stress response in metals differs greatly from that of polymers, e.g. rate and temperature effects are far more important, thus damage and wear mechanisms of polymers under cavitating flows are significantly different. In this work, heat-driven failure caused by viscous dissipation and void nucleation resulting from tensile stresses arising from stress wave reflections are investigated as two possible material failure mechanisms.

As a first step in developing a fundamental understanding of the cavitation erosion process on polymer surfaces, simulations are performed of the collapse of individual bubbles against a compliant surface e.g. metallic substrates with polyurea coatings. The surface response of collapse-driven impact loads is represented by a idealized, time-dependent, Gaussian pressure distribution on the surface. A two-dimensional distribution of load radii and durations is considered corresponding to characteristic of cavitating flows accelerated erosion experiments. Finite element simulations are performed to fit a response curve that relates the loading parameters to the energy dissipated in the coating and integrated with collapse statistics to generate an expected heat input into the coating.

The impulsive pressure, which is generated due to bubble collapse, impacts the material and generates intense shock waves. The stress waves within the material reflects by interaction with the substrate. A transient region of high tensile stress is produced by the interaction of these waves. Simulations suggests that maximum hydrostatic tension which cause failure of polyurea layer is observed in thick coating. Also, the dissipated viscous energy and corresponding temperature rise in a polyurea is calculated, and it is concluded that temperature has influence on deformation.

ABSTRACT Electronics especially mobile electronics such as smart phones, tablet PCs, notebooks and digital cameras are undergoing rapid development nowadays and have thoroughly changed our lives. With the requirement of more transistors, higher power, smaller size, lighter weight and even bendability, thermal management of these devices became one of the key challenges. Compared to active heat management system, heat pipe, which is a passive fluidic system, is considered promising to solve this problem. However, traditional heat pipes have size, weight and capillary limitation. Thus new type of heat pipe with smaller size, lighter weight and higher capillary pressure is needed. Nanofiber has been proved with superior properties and has been applied in multiple areas. This study discussed the possibility of applying nanofiber in heat pipe as new wick structure. In this study, a needleless electrospinning device with high productivity rate was built onsite to systematically investigate the effect of processing parameters on fiber properties as well as to generate nanofiber mat to evaluate its capability in electronics cooling. Polyethylene oxide (PEO) and Polyvinyl Alcohol (PVA) nanofibers were generated. Tensiometer was used for wettability measurement. The results show that independent parameters including spinneret type, working distance, solution concentration and polymer type are strongly correlated with fiber morphology compared to other parameters. The results also show that the fabricated nanofiber mat has high capillary pressure.

Production from a high pressure gas well at a high production-rate encounters the risk of operating near the choking condition for a compressible flow in porous media. The unbounded gas pressure gradient near the point of choking, which is located near the wellbore, generates an effective tensile stress on the porous rock frame. This tensile stress almost always exceeds the tensile strength of the rock and it causes a tensile failure of the rock, leading to wellbore instability. In a porous rock, not all pores are choked at the same flow rate, and when just one pore is choked, the flow through the entire porous medium should be considered choked as the gas pressure gradient at the point of choking becomes singular. This thesis investigates the choking condition for compressible gas flow in a single microscopic pore. Quasi-one-dimensional analysis and axisymmetric numerical simulations of compressible gas flow in a pore scale varicose tube with a number of bumps are carried out, and the local Mach number and pressure along the tube are computed for the flow near choking condition. The effects of tube length, inlet-to-outlet pressure ratio, the number of bumps and the amplitude of the bumps on the choking condition are obtained. These critical values provide guidance for avoiding the choking condition in practice.

The atomization of a liquid jet by a high speed cross-flowing gas has many applications such as gas turbines and augmentors. The mechanisms by which the liquid jet initially breaks up, however, are not well understood. Experimental studies suggest the dependence of spray properties on operating conditions and nozzle geom- etry. Detailed numerical simulations can offer better understanding of the underlying physical mechanisms that lead to the breakup of the injected liquid jet. In this work, detailed numerical simulation results of turbulent liquid jets injected into turbulent gaseous cross flows for different density ratios is presented. A finite volume, balanced force fractional step flow solver to solve the Navier-Stokes equations is employed and coupled to a Refined Level Set Grid method to follow the phase interface. To enable the simulation of atomization of high density ratio fluids, we ensure discrete consistency between the solution of the conservative momentum equation and the level set based continuity equation by employing the Consistent Rescaled Momentum Transport (CRMT) method. The impact of different inflow jet boundary conditions on different jet properties including jet penetration is analyzed and results are compared to those obtained experimentally by Brown & McDonell(2006). In addition, instability analysis is performed to find the most dominant insta- bility mechanism that causes the liquid jet to breakup. Linear instability analysis is achieved using linear theories for Rayleigh-Taylor and Kelvin- Helmholtz instabilities and non-linear analysis is performed using our flow solver with different inflow jet boundary conditions.

Derived from the necessity to increase testing capabilities of hybrid rocket motor (HRM) propulsion systems for Daedalus Astronautics at Arizona State University, a small-scale motor and test stand were designed and developed to characterize all components of the system. The motor is designed for simple integration and setup, such that both the forward-end enclosure and end cap can be easily removed for rapid integration of components during testing. Each of the components of the motor is removable allowing for a broad range of testing capabilities. While examining injectors and their potential it is thought ideal to obtain the highest regression rates and overall motor performance possible. The oxidizer and fuel are N2O and hydroxyl-terminated polybutadiene (HTPB), respectively, due to previous experience and simplicity. The injector designs, selected for the same reasons, are designed such that they vary only in the swirl angle. This system provides the platform for characterizing the effects of varying said swirl angle on HRM performance.