Dynamics of settling charged particles in turbulence : theory and experiments

It has been proposed that inertial clustering may lead to an increased collision rate of water droplets in clouds. Atmospheric clouds and electrosprays contain electrically charged particles embedded in turbulent flows, often under the influence of an externally imposed, approximately uniform gravitational or electric force. In this thesis, we present the investigation of charged inertial particles embedded in turbulence. We have developed a theoretical description for the dynamics of such systems of charged, sedimenting particles in turbulence, allowing radial distribution functions to be predicted for both monodisperse and bidisperse particle size distributions. The governing parameters are the particle Stokes number (particle inertial time scale relative to turbulence dissipation time scale), the Coulomb-turbulence parameter (ratio of Coulomb ’terminalar speed to turbulence dissipation velocity scale), and the settling parameter (the ratio of the gravitational terminal speed to turbulence dissipation velocity scale). For the monodispersion particles, The peak in the radial distribution function is well predicted by the balance between the particle terminal velocity under Coulomb repulsion and a time-averaged ’drift’ velocity obtained from the nonuniform sampling of fluid strain and rotation due to finite particle inertia. The theory is compared to measured radial distribution functions for water particles in homogeneous, isotropic air turbulence. The radial distribution functions are obtained from particle positions measured in three dimensions using digital holography. The measurements support the general theoretical expression, consisting of a power law increase in particle clustering due to particle response to dissipative turbulent eddies, modulated by an exponential electrostatic interaction term. Both terms are modified as a result of the gravitational diffusion-like term, and the role of ’gravity’ is explored by imposing a macroscopic uniform electric field to create an enhanced, effective gravity. The relation between the radial distribution functions and inward mean radial relative velocity is established for charged particles.

Comprehensive review and application of particle image velocimetry

For a fluid dynamics experimental flow measurement technique, particle image velocimetry (PIV) provides significant advantages over other measurement techniques in its field. In contrast to temperature and pressure based probe measurements or other laser diagnostic techniques including laser Doppler velocimetry (LDV) and phase Doppler particle analysis (PDPA), PIV is unique due to its whole field measurement capability, non-intrusive nature, and ability to collect a vast amount of experimental data in a short time frame providing both quantitative and qualitative insight. These properties make PIV a desirable measurement technique for studies encompassing a broad range of fluid dynamics applications. However, as an optical measurement technique, PIV also requires a substantial technical understanding and application experience to acquire consistent, reliable results. Both a technical understanding of particle image velocimetry and practical application experience are gained by applying a planar PIV system at Michigan Technological University’s Combustion Science Exploration Laboratory (CSEL) and Alternative Fuels Combustion Laboratory (AFCL). Here a PIV system was applied to non-reacting and reacting gaseous environments to make two component planar PIV as well as three component stereographic PIV flow field velocity measurements in conjunction with chemiluminescence imaging in the case of reacting flows. This thesis outlines near surface flow field characteristics in a tumble strip lined channel, three component velocity profiles of non-reacting and reacting swirled flow in a swirl stabilized lean condition premixed/prevaporized-fuel model gas turbine combustor operating on methane at 5-7 kW, and two component planar PIV measurements characterizing the AFCL’s 1.1 liter closed combustion chamber under dual fan driven turbulent mixing flow.

QUADRUPOLE LEVITATION OF PARTICLES IN A THERMODYNAMICALLY REALISTIC CLOUD ENVIRONMENT

Understanding clouds and their role in climate depends in part on our ability to understand how individual cloud particles respond to environmental conditions. Keeping this objective in mind, a quadrupole trap with thermodynamic control has been designed and constructed in order to create an environment conducive to studying clouds in the laboratory. The quadrupole trap allows a single cloud particle to be suspended for long times. The temperature and water vapor saturation ratio near the trapped particle is controlled by the flow of saturated air through a tube with a discontinuous wall temperature. The design has the unique aspect that the quadrupole electrodes are submerged in heat transfer fluid, completely isolated from the cylindrical levitation volume. This fluid is used in the thermodynamic system to cool the chamber to realistic cloud temperatures, and a heated section of the tube provides for the temperature discontinuity. Thus far, charged water droplets, ranging from about 30-70 microns in diameter have been levitated. In addition, the thermodynamic system has been shown to create the necessary thermal conditions that will create supersaturated conditions in subsequent experiments. These advances will help lead to the next generation of ice nucleation experiments, moving from hemispherical droplets on a substrate to a spherical droplet that is not in contact with any surface.

LATTICE BOLTZMANN METHOD AND CELLULAR AUTOMATA SIMULATION OF PARTICLE MOTION AND DEPOSITION IN 2-D CASE

This technical report discusses the application of the Lattice Boltzmann Method (LBM) and Cellular Automata (CA) simulation in fluid flow and particle deposition. The current work focuses on incompressible flow simulation passing cylinders, in which we incorporate the LBM D2Q9 and CA techniques to simulate the fluid flow and particle loading respectively. For the LBM part, the theories of boundary conditions are studied and verified using the Poiseuille flow test. For the CA part, several models regarding simulation of particles are explained. And a new Digital Differential Analyzer (DDA) algorithm is introduced to simulate particle motion in the Boolean model. The numerical results are compared with a previous probability velocity model by Masselot [Masselot 2000], which shows a satisfactory result.

Study of particle movement in the nearshore region of Lake Superior with radioisotope tracers

Biogeochemical processes in the coastal region, including the coastal area of the Great Lakes, are of great importance due to the complex physical, chemical and biological characteristics that differ from those on either the adjoining land or open water systems. Particle-reactive radioisotopes, both naturally occurring (210Pb, 210Po and 7Be) and man-made (137Cs), have proven to be useful tracers for these processes in many systems. However, a systematic isotope study on the northwest coast of the Keweenaw Peninsula in Lake Superior has not yet been performed. In this dissertation research, field sampling, laboratory measurements and numerical modeling were conducted to understand the biogeochemistry of the radioisotope tracers and some particulate-related coastal processes. In the first part of the dissertation, radioisotope activities of 210Po and 210Pb in a variability of samples (dissolved, suspended particle, sediment trap materials, surficial sediment) were measured. A completed picture of the distribution and disequilibrium of this pair of isotopes was drawn. The application of a simple box model utilizing these field observations reveals short isotope residence times in the water column and a significant contribution of sediment resuspension (for both particles and isotopes). The results imply a highly dynamic coastal region. In the second part of this dissertation, this conclusion is examined further. Based on intensive sediment coring, the spatial distribution of isotope inventories (mainly 210Pb, 137Cs and 7Be) in the nearshore region was determined. Isotope-based focusing factors categorized most of the sampling sites as non- or temporary depositional zones. A twodimensional steady-state box-in-series model was developed and applied to individual transects with the 210Pb inventories as model input. The modeling framework included both water column and upper sediments down to the depth of unsupported 210Pb penetration. The model was used to predict isotope residence times and cross-margin fluxes of sediments and isotopes at different locations along each transect. The time scale for sediment focusing from the nearshore to offshore regions of the transect was on the order of 10 years. The possibility of sediment longshore movement was indicated by high inventory ratios of 137Cs: 210Pb. Local deposition of fine particles, including fresh organic carbon, may explain the observed distribution of benthic organisms such as Diporeia. In the last part of this dissertation, isotope tracers, 210Pb and 210Po, were coupled into a hydrodynamic model for Lake Superior. The model was modified from an existing 2-D finite difference physical-biological model which has previously been successfully applied on Lake Superior. Using the field results from part one of this dissertation as initial conditions, the model was used to predict the isotope distribution in the water column; reasonable results were achieved. The modeling experiments demonstrated the potential for using a hydrodynamic model to study radioisotope biogeochemistry in the lake, although further refinements are necessary.