In particulate flows, the flow inertia impacts the motion and size distribution of the particles and this in turn, has a strong implication on global behavior of the emulsions such as their rheological properties. As such, the central goal of most of the investigations on dispersed multiphase flow, so far, has been to understand the phase distribution of particles and to correlate the global behavior of the system with this parameter. For pressure-driven particulate flows in a channel, it is known that the velocity gradient in the channel leads to a lateral force whose magnitude and direction depends on the viscosity and density ratios of the fluids and the drop deformation. This lateral (lift) force is the primary reason behind the various observed modes of phase distribution of the particles. Unfortunately, most of the studies conducted so far have been concerned with the solid particles and for flows at low to moderate Reynolds numbers. Little is known about the dynamics of deformable drops at high Reynolds numbers. The goal of this study is to bridge the gap by direct numerical simulations. A front tracking/finite difference technique is used to solve the Navier-Stokes equations in the fluids inside and outside of the drops. Initially, the drops are randomly distributed in the computational domain their evolutions are followed for a sufficiently long time so that the system reaches a quasi-steady state. The statistics about the flow then will be extracted. The flow inertia is increased incrementally by increasing the pressure gradient.

Potable water distribution systems are broken into major and minor distribution networks. Major water distribution networks refer to large-scale municipal pipe systems extending from the treatment plant to the upstream node of the water service line for buildings. Minor water distribution systems, also referred to as plumbing water distribution systems, run from the upstream node of the water service line to all interior plumbing fixtures and demand nodes associated with the building. Most texts and research papers focus on major systems, while only a small number of documents are available concerning the design and analysis of minor systems. In general, the available minor system documents are quite prescriptive in nature. This thesis presents a comprehensive evaluation of contemporary plumbing water distribution system design. All underlying theory is explained and advantages and drawbacks are discussed. Furthermore, contemporary methods for designing minor distribution systems have come under recent scrutiny. Issues have been raised regarding the accuracy of water demand estimation procedures for plumbing systems, namely, Hunter's method. Demand estimates are crucial for designing minor piping systems. The formulation and application of a pressure-driven design approach to replace Hunter-based design methods is presented. EPANET, a commonly used hydraulic modeling software package, is utilized to evaluate network behavior. Example applications are presented to illustrate the robustness of a pressure-driven approach, while also allowing the evaluation of plumbing water distribution system performance under worst-case loading conditions. / Master of Science

Microfluidic devices could find applications in many areas, such as BioMEMs, miniature fuel cells and microfluidic cooling of electronic circuitry. One of the important considerations of microfluidic device in analytical and bioanalytical chemistry is the dispersion of solute. In this study, we have developed an analytical solution, which considers the axial dispersion of a solute along the flow direction, to simulate convection and diffusion transport in a pressure driven creeping flow for a rectangular shape slit. During flow, the balance of competing effects of diffusion (especially cross-section diffusion) and convective diffusion in the flow direction are investigated. / Singapore-MIT Alliance (SMA)

Electroosmotic flow is a convenient mechanism for transporting polar fluid in a microfluidic device. The flow is generated through the application of an external electric field that acts on the free charges that exists in a thin Debye layer at the channel walls. The charge on the wall is due to the chemistry of the solid-fluid interface, and it can vary along the channel, e.g. due to modification of the wall. This investigation focuses on the simulation of the electroosmotic flow (EOF) profile in a cylindrical microchannel with step change in zeta potential. The modified Navier-Stoke equation governing the velocity field and a non-linear two-dimensional Poisson-Boltzmann equation governing the electrical double-layer (EDL) field distribution are solved numerically using finite control-volume method. Continuities of flow rate and electric current are enforced resulting in a non-uniform electrical field and pressure gradient distribution along the channel. The resulting parabolic velocity distribution at the junction of the step change in zeta potential, which is more typical of a pressure-driven velocity flow profile, is obtained. / Singapore-MIT Alliance (SMA)

Interruption to water services and low water pressure conditions are commonly observed problems in water distribution systems (WDSs). Of particular concern are the unplanned events, such as pipe bursts. The current regulation in the UK requires water utilities to provide reliable water service to consumers resulting in as little as possible interruptions and of as short possible duration. All this pushes water utilities toward developing and using smarter responses to these events, based on advanced tools and solutions. All with the aim to change network management style from reactive to a proactive, and reduce water losses, optimize energy use and provide better services for consumers. This thesis presents a novel methodology for efficient and effective operational, short time response to an unplanned failure event (such as pipe burst) in a WDS. The proposed automated, near real-time operational response methodology consists of isolating the failure event followed by the recovery of the affected system area by restoring the flows and pressures to normal conditions. The isolation is typically achieved by manipulating the relevant on/off valves that are located closely to the event location. The recovery involves selecting an optimal combination of suitable operational network interventions. These are selected from a number of possible options with the aim to reduce the negative impact of the failure over a pre-specified time horizon. The intervention options considered here include isolation valve manipulations, changing the pressure reducing valve’s (PRV) outlet pressure and installation and use of temporary overland bypasses from a nearby hydrant(s) in an adjacent, unaffected part of the network. The optimal mix of interventions is identified by using a multi-objective optimization approach driven by the minimization of the negative impact on the consumers and the minimization of the corresponding number of operational interventions (which acts as a surrogate for operational costs). The negative impact of a failure event was quantified here as a volume of water undelivered to consumers and was estimated by using a newly developed pressure-driven model (PDM) based hydraulic solver. The PDM based hydraulic solver was validated on a number of benchmark and real-life networks under different flow conditions. The results obtained clearly demonstrate its advantages when compared to a number of existing methods. The key advantages include the simplicity of its implementation and the ability to predict network pressures and flows in a consistently accurate, numerically stable and computationally efficient manner under both pressure-deficient and normal-flow conditions and in both steady-state and extended period simulations. The new real-time operational response methodology was applied to a real world water distribution network of D-Town. The results obtained demonstrate the effectiveness of the proposed methodology in identifying the Pareto optimal network type intervention strategies that could be ultimately presented to the control room operator for making a suitable decision in near real-time.

For the past two decades, flow cytometry has been widely used as a powerful analysis tool for the diagnosis of many diseases due to its ability to count, characterize and sort cells. However, conventional flow cytometers are often bulky, expensive and complicated because sophisticated fluidic, electronic and optical systems are required to realize the functions of flow cytometry. The high cost and the complexity in operation and maintenance associated with flow cytometers as well as the large size have limited its use. In recent years, the rapid development of microfluidics-based lab-on-a-chip technology has created a new pathway for flow cytometry. Microfluidic devices allow for the integration of multiple liquid handling processes required in the diagnostic assays, such as pumping, metering, sampling, dispensing, sequential loading and washing. These lab-on-a-chip solutions have been recognized as an opportunity to bring portable, accurate and sensitive diagnostic tests to the flow cytometry. However, most current microfluidic flow cytometry devices are micro- only in the microfluidic chip, the rest of most apparatuses are still large and costly, usually involving tubes, microscopes, lasers and mechanical pumps. Therefore, the objective of this study is to develop a novel lab-on-a-chip system based on the electrokinetically-induced pressure-driven flow and dual-wavelength fluorescent detection, which lights a promising pathway for making a real portable, compact, low-cost microfluidic flow cytometry device. In this study, the core of this microfluidic system is the custom-designed PDMS (polydimethylsiloxane) microchip. A novel method was applied to generate the electrokinetically-induced pressure-driven flow in a T-shaped microchannel using parameters settings that had been optimized by numerical study. This method combined both the electrokinetic pumping force and the pressure pumping force to eliminate their shortcomings associated with the use of each force alone. This is the fundamental of my study. By using this microchip, the size of the fluidic control subsystem is reduced significantly. Furthermore, the dual-wavelength fluorescent detection strategy is proposed in this thesis. On the optical detection side, excitation lights of two different wavelengths are provided by a single LED (light-emitting diode) from one side of the microchannel. Then the two emission lights are captured individually by two photo-detectors placed on the top and the bottom of the microchip. Compared with other microfluidic detection devices reported in the literatures that use lasers or PMTs (Photomultiplier tubes), this design allows for a significant reduction of 90% in the volume and cost. As another important part of my thesis research, a novel flow focusing method that allows the hydrodynamic focusing in a T-shaped microchannel with two sheath flows is developed. This method solves the biggest obstacle which exists in current microfluidic flow cytometry devices. In this method, no external pumps, valves and tubing are involved in the system. Although substantial progress has been made in current microfluidic flow cytometry, there is still a need for a low-cost, compact, portable microfluidic devices, especially in low-resource settings as well as the developing world for POC (point-of-care) diagnosis and analysis. This thesis work has made a great achievement towards the final goal.

This thesis presents the results and analyses of computational studies on certain problems of combustion instability in solid propellants. Specifically, effects of relaxing certain assumptions made in previous models of unsteady burning of solid propellants are investigated. Knowledge of unsteady burning of solid propellants is essential in studying the phenomenon of combustion instability in solid propellant rocket motors. In Chapter 1, an introduction to different types of unsteady combustion investigated in this thesis, such as 1) intrinsic instability, 2) pressure-driven dynamic burning, 3) extinction by depressurization, and 4) L* -instability, is given. Also, a review of previous experimental and theoretical studies of these phenomena is presented. From this review it is concluded that all the previous studies, which investigated the unsteady combustion of solid propellants, made one or more of the following assumptions: 1) quasi-steady gas-phase (QSG), 2) quasi-steady condensed phase reaction zone (QSC), 3) small perturbations, and 4) unity Lewis number. These assumptions limit the validity of the results obtained with such models to: 1) relatively low frequencies (< 1 kHz) of pressure oscillations and 2) small deviations in pressure from its steady state or mean values. The objectives of the present thesis are formulated based on the above conclusions. These are: 1) to develop a nonlinear numerical model of unsteady solid propellant combustion, 2) to relax the assumptions of QSG and QSC, 3) to study the consequent effects on the intrinsic instability and pressure-driven dynamic burning, and 4) to investigate the L* -instability in solid propellant rocket motors. In Chapter 2, a nonlinear numerical model, which relaxes the QSG and QSC assumptions, is set up. The transformation and nondimensionalization of the governing equations are presented. The numerical technique based on the method of operator-splitting, used to solve the governing equations is described. In Chapter 3, the effect of relaxing the QSG assumption on the intrinsic instability is investigated. The stable and unstable solutions are obtained for parameters corresponding to a typical composite propellant. The stability boundary, in terms of the nondimensional parameters identified by Denison and Baum (1961), is predicted using the present model. This is compared with the stability boundary obtained by previous linear stability theories, based on activation energy asymptotics in the gas-phase, which employed QSC and/or QSG assumptions. It is found that in the limit of large activation energy and low frequencies, present result approaches the previous theoretical results. This serves as a validation of the present method of solution. It is confirmed that relaxing the QSG assumption widens the stable region. However, it is found that a distributed reaction in the gas-phase destabilizes the burning. The effect of non-unity Lewis number on the stability boundary is also investigated. It is found that at parametric values corresponding to low pressures and large flame stand-off distances, small amplitude, high frequency (at frequencies near the characteristic frequency of the gas-phase) oscillations in burning rate appear when the Lewis number is greater than one. In Chapter 4, the effect of relaxing the QSG assumption is further investigated with respect to the pressure-driven dynamic burning. Comparison of the pressure-driven frequency response function, Rp, obtained with the present model, both in the QSG and non-QSG framework, with those obtained with previous linear stability theories invoking QSG and QSC assumptions are made. As the frequency of pressure oscillations approaches zero, |RP| predicted using present models approached the value obtained by previous theoretical studies. Also, it is confirmed that the effect of relaxing QSG is to decrease the |Rp| at frequencies near the first resonant frequency. Moreover, relaxing QSG assumption produces a second resonant peak in |Rp| at a frequency near the characteristic frequency of the gas-phase. Further, Rp calculated using the present model is compared with that obtained by a previous linear theory which relaxed the QSG assumption. The two models predicted the same resonant frequencies in the limit of small amplitudes of pressure oscillations. Finally, it is found that the effect of large amplitude of pressure oscillations is to introduce higher harmonics in the burning rate and to reduce the mean burning rate. In Chapter 5, first the present non-QSC model is validated by comparing its results with that of a previous non-QSC model for radiation-driven burning. The model is further validated for steady burning results by comparing with experimental data for a double base propellant (DBP). Then, the effect of relaxing the QSC assumption on steady state solution is investigated. It is found that, even in the presence of a strong gas-phase heat feedback, QSC assumption is valid for moderately large values of condensed phase Zel'dovich number, as far as steady state solution is concerned. However, for pressure-driven dynamic burning, relaxing the QSC assumption is found to increase |RP| at all frequencies. The error due to QSC assumption is found to become significant, either when |Rp| is large or as the driving frequency approaches the characteristic frequency of the condensed phase reaction zone. The predicted real part of the response function is quantitatively compared with experimental data for DBP. The comparison seems to be better with a value of condensed phase activation energy higher than that suggested by Zenin (1992). In Chapter 6, burning rate transients for a DBP during exponential depressurization are computed using non-QSG and non-QSC models. Salient features of extinction and combustion recovery are predicted. The predicted critical initial depressurization rate, (dp/dt)i, is found to decrease markedly when the QSC assumption is relaxed. The effect of initial pressure level on critical (dp/dt)i is studied. It is found that the critical (dp/dt)i decreases with the initial pressure. Also, the overshoot of burning rate during combustion recovery is found to be relatively large with low initial pressures. However as the initial pressure approached the final pressure, the reduction in initial pressure causes a large increase in the critical (dp/dt)i. No extinction is found to occur when the initial pressure is very close to the final pressure. In Chapter 7, a numerical model is developed to simulate the L* -instability in solid propellant motors. This model includes a) the propellant burning model that takes into account nonlinear pressure oscillations and that takes into account an unsteady gas- and condensed phase, and b) a combustor model that allows pressure and temperature oscillations of finite amplitude. Various regimes of L* -burning of a motor, with a typical composite propellant, namely 1) steady burning, 2) oscillatory burning leading to steady state, 3) oscillatory burning leading to extinction, 4) reignition and 5) chuffing are predicted. The predicted dependence of frequency of L* -oscillations on mean pressure is compared with one set of available experimental data. It is found that proper modeling of the radiation heat flux from the chamber walls to the burning surface may be important to predict the re-ignition. In Chapter 8, the main conclusions of the present study are summarized. Certain suggestions for possible future studies to enhance the understanding of dynamic combustion of solid propellants are also given.

The flow profiles of pressure-driven and electro-kinetic driven flows were compared for a microfluidic chip. It was found that the pressure-driven flow had a parabolic profile while the electro-kinetic flow had a plug shaped flow profile. The measured velocities were similar to those determined by the Poiseuille flow model and the Helmholtz-Smoltchowski equation. Flow uniformity is very important for control in microfluidic mixers. Parabolic flow profiles lead to inconsistent reactions while the more uniform plug shape flow allow for a more steady reaction across the channel. Previous work had been performed to measure the flow of a solution of fluorescent polystyrene beads in PDMS channels using a laser confocal microscope. This showed that particles easily stuck to the channel making it difficult to measure over time. In addition, bubble formation in the channel made measuring velocities difficult. Current work used a LabSmith Video Synchronized microscope with software to measure the flow rates at different areas of the channel. Solutions of fluorescent polystyrene beads were used to visually observe the flow within a channel under a microscope. Four different channels were used for the pressure-driven flows of varying dimensions and materials. The channel with the best measured profile was also measured under electro-kinetic flow. A LabSmith High Voltage Sequencer was used to apply a voltage across the channel for electro-kinetic measurements. This research confirmed the different flow profiles under pressure-driven and electro-kinetic driven flow. Future work can be done to determine how this effects mixing in the channels.

This thesis is motivated by problems studied as part of the Diavik Waste Rock Pile Project. Located at the Diavik Diamond Mine in the Northwest Territories, with academic support from the University of Waterloo, the University of Alberta, and the University of British Columbia, this project focuses on constructing mine waste rock piles and studying their physical and chemical properties and the transport processes within them. One of the main reasons for this investigation is to determine the effect of environmental factors on acid mine drainage (AMD) due to sulfide oxidation and the potential environmental impact of AMD. This research is concerned with modeling pressure driven flow through waste rock piles. Unfortunately, because of the irregular shape of the piles, very little data for fluid flow about such an obstacle exists, and the numerical techniques available to work with this domain are limited. Since this restricts the study of the mathematics behind the flow, this thesis focuses on a cylindrical domain, since flow past a solid cylinder has been subjected to many years of study. The cylindrical domain also facilitates the implementation of a pseudo-spectral method. This thesis examines a pressure driven flow through a cylinder of variable permeability subject to turbulent forcing. An equation for the steady flow of an incompressible fluid through a variable permeability porous medium is derived based on Darcy's law, and a pseudo-spectral model is designed to solve the problem. An unsteady time-dependent model for a slightly compressible fluid is then presented, and the unsteady flow through a constant permeability cylinder is examined. The steady results are compared with a finite element model on a trapezoidal domain, which provides a better depiction of a waste rock pile cross section.

This thesis is motivated by problems studied as part of the Diavik Waste Rock Pile Project. Located at the Diavik Diamond Mine in the Northwest Territories, with academic support from the University of Waterloo, the University of Alberta, and the University of British Columbia, this project focuses on constructing mine waste rock piles and studying their physical and chemical properties and the transport processes within them. One of the main reasons for this investigation is to determine the effect of environmental factors on acid mine drainage (AMD) due to sulfide oxidation and the potential environmental impact of AMD. This research is concerned with modeling pressure driven flow through waste rock piles. Unfortunately, because of the irregular shape of the piles, very little data for fluid flow about such an obstacle exists, and the numerical techniques available to work with this domain are limited. Since this restricts the study of the mathematics behind the flow, this thesis focuses on a cylindrical domain, since flow past a solid cylinder has been subjected to many years of study. The cylindrical domain also facilitates the implementation of a pseudo-spectral method. This thesis examines a pressure driven flow through a cylinder of variable permeability subject to turbulent forcing. An equation for the steady flow of an incompressible fluid through a variable permeability porous medium is derived based on Darcy's law, and a pseudo-spectral model is designed to solve the problem. An unsteady time-dependent model for a slightly compressible fluid is then presented, and the unsteady flow through a constant permeability cylinder is examined. The steady results are compared with a finite element model on a trapezoidal domain, which provides a better depiction of a waste rock pile cross section.