Electric fields can be used to change and control flame properties, for example changing flame speed, enhancing flame stability, or reducing pollutant emission. The ions generated in flames are believed to play the primary role. Although experiments have been carried out to study electric field enhanced combustion, they are not sufficient to explain how the ions in a flame are affected by an electric field. It is therefore necessary to investigate the problem through numerical simulations. In the present work, the electric structure of stabilized CH4/air premixed flames at atmospheric pressure within a direct current field is studied using numerical simulations. This study consists of three parts. First, the transport equations are derived from the Boltzmann kinetic equation for each individual species. Second, a general method for computing the diffusivity and mobility of ions in a gas mixture is introduced. Third, the mechanisms for neutral and charged species are improved to give better predictions of the concentrations of charged species, based on experimental data. Following from this, comprehensive numerical results are presented, including the concentrations and fluxes of charged species, the distributions of the electric field and electric potential, and the electric current-voltage relation. Two new concepts introduced with the numerical results are the plasma sheath and dead zone in the premixed flame. A reactive plasma sheath and a Boltzmann relation sheath are discovered in the region near the electrodes. The plasma sheath penetrates into the flame gas when a voltage is applied, and penetrating further if the voltage is higher. The zone outside the region of sheath penetration is defined as the dead zone. With the two concepts, analytical solutions for the electric field, electric potential and current-voltage curve are derived. The solutions directly describe the electric structure of a premixed flame subject to a DC field. These analytical solutions, together with the discovery of the plasma sheath and dead zone in flames, are the novel contributions of this work.

The first objective of the dissertation is to develop a suitable reduced order model capable of investigating the nonlinear mechanical behavior of von-Karman plates under electrostatic actuation. The second objective is to investigate the nonlinear static and dynamic behavior of rectangular microplates under small and large actuating forces. In the first part, we present and compare various approaches to develop reduced order models for the nonlinear von-Karman rectangular microplates actuated by nonlinear electrostatic forces. The reduced-order models aim to investigate the static and dynamic behavior of the plate under small and large actuation forces. A fully clamped microplate is considered. Different types of basis functions are used in conjunction with the Galerkin method to discretize the governing equations. First we investigate the convergence with the number of modes retained in the model. Then for validation purpose, a comparison of the static results is made with the results calculated by a nonlinear finite element model. The linear eigenvalue problem for the plate under the electrostatic force is solved for a wide range of voltages up to pull-in. In the second part, we present an investigation of the static and dynamic behavior of a fully clamped microplate. We investigate the effect of different non-dimensional design parameters on the static response. The forced-vibration response of the plate is then investigated when the plate is excited by a harmonic AC load superimposed to a DC load. The dynamic behavior is examined near the primary and secondary (superharmonic and subharmonic) resonances. The microplate shows a strong hardening behavior due to the cubic nonlinearity of midplane stretching. However, the behavior switches to softening as the DC load is increased. Next, near-square plates are studied to understand the effect of geometric imperfections of microplates. In the final part of the dissertation, we investigate the mechanical behavior of initially curved microplates. Microplates often experience an initial curvature imperfection, due to the micro fabrication process, which affects significantly their mechanical behavior. In this case a clamped-free-clamped-free microplate is considered. We validate the reduced order model by comparing the calculated static behavior and the fundamental natural frequency with those computed by a finite element model. As case studies, we consider two commonly encountered profiles of the initial curvature imperfection and study their effects on both the static and dynamic responses of the microplates. Next, an initially curved microplate made of silicon nitride is studied. The static behavior of the microplate is investigated when applying a DC voltage. Then, the dynamic behavior of the microplate is examined under the application of a harmonic AC voltage, superimposed to a DC voltage. Simulation results calculated by the reduced order model are compared with experimental data for model validation purpose, which show good agreement.

Soot is a common byproduct of hydrocarbon based combustion systems. It poses a risk to human and environmental health, and can negatively or positively affect combustor performance. As a result, there is significant interest in understanding soot formation in order to better control it. More recently, the need to study soot formation in engine relevant conditions has become apparent. One engine relevant parameter that has had little focus is the ambient pressure. This body of work focuses on the formation of soot in elevated pressure environments, and a number of investigations are carried out with this purpose. Laminar coflow diffusion flames are used as steady, simple soot producers. First, a commonly studied flame configuration is further characterized. Coflow flames are frequently used for fundamental flame studies, particularly at elevated pressures. However, they are more susceptible to buoyancy induced instabilities at elevated pressures. The velocity of the coflow is known to have an effect on flame stability and soot formation, though these have not been characterized at elevated pressures. A series of flames are investigated covering a range of flowrates, pressures, and nozzle diameters. The stability limits of coflow flames in this range is investigated. Additionally, an alternative strategy for scaling these flames to elevated pressures is proposed. Finally, the effect of coflow rate on soot formation is evaluated. Identification of fundamental flames for coordinated research can facilitate our understanding of soot formation. The next study of this work focuses on adding soot concentration and particle size information to an existing fundamental flame dataset for the purpose of numerical model validation. Soot volume fraction and average particle diameters are successfully measured in nitrogen-diluted ethylene-air laminar coflow flames at pressures of 4, 8, 12, and 16 atm. An increase in particle size with pressure is found up to 12 atm, where particle sizes plateau. Particle size in the annulus is more sensitive to pressure. Next, the development of an alternative particle size measuring technique is studied. Time Resolved Laser Induced Incandescence (TiRe-LII) is a commonly used technique to measure soot concentrations and particle size at atmospheric pressure. However, Laser Induced Incandescence (LII) models suffer from an incomplete understanding of the effects of elevated pressures on the absorption, annealing, and cooling of soot. The present study focuses on what affect the laser temporal pulse shape and duration may have on particle sizing. TiRe-LII in flames at 1 and 15 bar is carried out, using laser pulses with tophat or Gaussian temporal profiles of varying duration. Mono-disperse equivalent primary particle diameters are calculated using the KAUST LII model. Little difference in particle sizing is found for different laser pulses. However, this data will be useful for validating the KAUST LII model when absorption and poly-dispersion are accounted for. In an effort to move one step closer to logistical fuel studies, the sooting tendencies of a number of liquid fuels are studied at pressures up to 10. Of parallel relevance, a sooting index for surrogate development is evaluated for elevated pressure applications. The Yield Sooting Index (YSI) methodology is applied to 11 normal, cyclic, and branched alkanes. When referencing to two n-alkane fuels, the YSI of n-alkanes determined at atmospheric pressures accurately reflects the relative sooting tendencies of these fuels at elevated pressures. The relative sooting tendencies of cyclo- and methyl-alkanes have a lower pressure sensitivity than n-alkanes.

The impact of a spherical object onto a surface of a liquid, solid or granular material, is a configuration which occurs in numerous industrial and natural phenomena. The resulting dynamics can produce complex outcomes and often occur on very short time-scales. Their study thereby requires high-speed video imaging, as is done herein. This three-part dissertation investigates widely disparate but kindred impact configurations, where the impacting object is a solid steel sphere, or a molten metal droplet. The substrate, on the other hand, is either granular material, a liquid, or solid ice. Therefore both fluid mechanics and thermodynamics play a key role in some of these dynamics. Part I, investigates the penetration depth of a steel sphere which impacts onto a granular bed containing a mixture of grains of two different sizes. The addition of smaller grains within a bed of larger grains can promote a “lubrication” effect and deeper penetration of the sphere. However, there needs to be enough mass fraction of the smaller grains so that they get lodged between the larger grains and are not simply like isolated rattlers inside the voids between the larger grains. This lubrication occurs even though the addition of the small grains increases the overall packing fraction of the bed. We compare the enhanced penetration for the mixtures to a simple interpolative model based on the results for monodispersed media of the constitutive sizes. The strongest lubrication is observed for large irregular shaped Ottawa sand grains, which are seeded with small spherical glass beads. Part II, tackles the topic of a molten metal drop impacting onto a pool of water. When the drop temperature is far above the boiling temperature of water, a continuous vapor layer can form at the interface between the metal and water, in what is called the Leidenfrost phenomenon. This vapor layer can become unstable forming what is called a vapor explosion, which can break up the molten metal drop. We study the details of these explosions and characterize the metal debris. We contrast the results for two different metals, i.e. tin and a special metal alloy called Field’s metal. For tin the drop solidifies and forms a porous foam-like solid, whereas the Field’s metal breaks up into a multitude of spherical beads, with a range of sizes as small as a few microns. We attribute this difference to the much lower melting point of the Field’s metal, which is only 60oC, compared to 230oC for the tin. This allows more fragmentation of the Field’s metal drop before it solidifies. When the temperature of the impacting metal is increased, high-speed imaging reveals a sequence of up to three vapor explosions, each of increasing intensity. We measure the acceleration of the vapor interface and compare the size-distribution of the microbeads to the fastest growing instability mode of the corresponding Rayleigh-Taylor instability. Part III, investigates the coefficient of restitution when a steel sphere impacts on an ice surface. As observed in earlier studies the restitution coefficient is largest for the smallest impact velocities, where the surface is not greatly fragmented. Our focus is on greatly heating the sphere up to 400oC to investigate how the thermal load affects the short term interaction of the sphere with the ice. We see a clear trend where hotter spheres rebound less than cold spheres. We also track the speed of ice-fragments ejected during the earliest stages of the impact.

Future thrusts for gasoline engine development can be broadly summarized into two categories: (i) efficiency improvements in conventional spark ignition engines, and (ii) development of advance compression ignition (ACI) concepts. Efficiency improvements in conventional spark ignition engines requires downsizing (and turbocharging) which may be achieved by using high octane gasolines, whereas, low octane gasolines fuels are anticipated for ACI concepts. The current work provides the essential combustion kinetic data, targeting both thrusts, that is needed to develop high fidelity gasoline surrogate mechanisms and surrogate complexity guidelines. Ignition delay times of a wide range of certified gasolines and surrogates are reported here. These measurements were performed in shock tubes and rapid compression machines over a wide range of experimental conditions (650 – 1250 K, 10 – 40 bar) relevant to internal combustion engines. Using the measured the data and chemical kinetic analyses, the surrogate complexity requirements for these gasolines in homogeneous environments are specified. For the discussions presented here, gasolines are classified into three categories: (i) Low octane gasolines including Saudi Aramco’s light naphtha fuel (anti-knock index, AKI = (RON + MON)/2 = 64; Sensitivity (S) = RON – MON = 1), certified FACE (Fuels for Advanced Combustion Engines) gasoline I and J (AKI ~ 70, S = 0.7 and 3 respectively), and their Primary Reference Fuels (PRF, mixtures of n-heptane and iso-octane) and multi-component surrogates. (ii) Mid octane gasolines including FACE A and C (AKI ~ 84, S ~ 0 and 1 respectively) and their PRF surrogates. Laser absorption measurements of intermediate and product species formed during gasoline/surrogate oxidation are also reported. (iii) A wide range of n-heptane/iso-octane/toluene (TPRF) blends to adequately represent the octane and sensitivity requirements of high octane gasolines including FACE gasoline F and G (AKI ~ 91, S = 5.6 and 11 respectively) and certified Haltermann (AKI ~ 87, S = 7.6) and Coryton (AKI ~ 92, S = 10.9) gasolines. To assess conditions where shock tubes may not be ideal devices for ignition delay measurements, this work also presents a detailed discussion on shock tube pre-ignition affected ignition data and the ignition regimes in homogeneous environments. The shock tube studies on pre-ignition and associated bulk ignition advance may help engines research community understand and control super-knock events.

With global emission regulations becoming stringent, development of new combustion technologies that meet future emission regulations is essential. In this vein, this dissertation presents the application of sensitive diagnostic tools to validate and improve chemical kinetic mechanisms that play a fundamental role in the design of new combustion technologies. First, a novel high sensitivity laser-based sensor with a wide frequency tuning range (900 – 1000 cm-1) was developed utilizing pulsed cavity ringdown spectroscopy (CRDS) technique. The novel laser-based sensor was illustrated by measuring trace amounts of multiple combustion intermediates, namely ethylene, propene, allene, and 1-butene in a static cell at ambient conditions. Subsequently, pulsed CRDS technique was utilized to develop an ultra-fast, high sensitivity diagnostic to monitor trace concentrations of ethylene in shock tube pyrolysis experiments. This diagnostic represented the first ever successful application of CRDS technique to transient species measurements in a shock tube. The high sensitivity and fast time response (10μs) diagnostic may be utilized for measuring other key neutrals and radicals which are crucial in the oxidation chemistry of practical fuels. Secondly, a quadrupole mass spectrometer (QMS) was employed to measure relative cation mole fractions in atmospheric and low-pressure (30 Torr) flames of methane/oxygen diluted in argon. Lean, stoichiometric and rich flames were 4 examined to evaluate the dependence of ion chemistry on flame stoichiometry. Spatial distribution of cations was compared with predictions of an existing ion chemistry model. Based on the extensive measurements carried out in this work, modifications were suggested to improve the ion chemistry model to enhance the fidelity of such mechanisms. In-depth understanding of flame ion chemistry is vital to model the interaction of flames with electric fields and thereby pave the way to enable active combustion control for increased efficiency and reduced emissions. Finally, a compact fast time-response time-of-flight mass spectrometer (TOFMS) was coupled to the shock tube through a pin-hole end-wall to enable timeresolved species concentration measurements. This diagnostic tool was demonstrated by investigating the decomposition of 1,3,5-trioxane over a wide range of shock conditions. Reaction rate coefficients were extracted by the best fit to the experimentally measured species time-histories. TOF-MS coupled to the shock tube is an ideal diagnostic tool for developing kinetic mechanisms for future fuels due to its ability to simultaneously measure several species during fuel pyrolysis/oxidation processes.

The free-surface impact of spheres is important to several applications in the military, industry and sports such as the water-entry of torpedoes, dip-coating procedures and slamming of boats. This two-part thesis attempts to explore this field by investigating cavity formation during the impact of spheres with deep liquid pools and cavitation in thin ultra-viscous films. Part I reports results from an experimental study on the formation of stable- streamlined and helical cavity wakes following the free-surface impact of heated Leidenfrost spheres. The Leidenfrost effect encapsulates the sphere by a vapor layer to prevent any physical contact with the surrounding liquid. This phenomenon is essential for the pacification of acoustic rippling along the cavity interface to result in a stable-streamlined cavity wake. Such a streamlined configuration experiences drag coefficients an order of magnitude lower than those acting on room temperature spheres. A striking observation is the formation of helical cavities which occur for impact Reynolds numbers Re0 ≳ 1.4 × 105 and are characterized by multiple interfacial ridges, stemming from and rotating synchronously about an evident contact line around the sphere equator. This helical configuration has 40-55% smaller overall force coefficients than those obtained in the formation of stable cavity wakes. Part II of this thesis investigates the inception of cavitation and resulting structures when a sphere collides with a solid surface covered with a layer of non-Newtonian liquid having kinematic viscosities of up to v0 = 20,000,000 cSt. The existence of shear-stress- induced cavitation during sphere approach towards the base wall (i.e. the pressurization stage) in ultra-viscous films is shown using a synchronized dual-view high-speed imaging system. In addition, cavitation by depressurization is noted for a new class of non-contact cases whereby the sphere rebounds without any prior contact with the solid wall. Horizontal shear rates calculated using particle image velocimtery (PIV) measurements reveal the apparent fluid viscosity to vary substantially as the sphere approaches and rebounds away from the base wall. A theoretical model based on the lubrication assumption is also solved for the squeeze flow in the regime identified for shear-induced cavity events to investigate the criterion for cavity inception in further detail.

Autoignition characteristics of pre-vaporized iso-octane, primary reference fuels, gasolines, and dimethyl ether (DME) have been investigated experimentally in a coflow with elevated temperature of air. With the coflow air at relatively low initial temperatures below autoignition temperature Tauto, an external ignition source was required to stabilize the flame. Non-autoignited lifted flames had tribrachial edge structures and their liftoff heights correlated well with the jet velocity scaled by the stoichiometric laminar burning velocity, indicating the importance of the edge propagation speed on flame stabilization balanced with local flow velocity. At high initial temperatures over Tauto, the autoignited flames were stabilized without requiring an external ignition source. The autoignited lifted flames exhibited either tribrachial edge structures or Mild combustion behaviors depending on the level of fuel dilution. For the iso-octane and n-heptane fuels, two distinct transition behaviors were observed in the autoignition regime from a nozzle-attached flame to a lifted tribrachial-edge flame and then a sudden transition to lifted Mild combustion as the jet velocity increased at a certain fuel dilution level. The liftoff data of the autoignited flames with tribrachial edges were analyzed based on calculated ignition delay times for the pre-vaporized fuels. Analysis of the experimental data suggested that ignition delay time may be much less sensitive to initial temperature under atmospheric pressure conditions as compared with predictions. For the gasoline fuels for advanced combustion engines (FACEs), and primary reference fuels (PRFs), autoignited liftoff data were correlated with Research Octane Number and Cetane Number. For the DME fuel, planar laser-induced fluorescence (PLIF) of formaldehyde (CH2O) and CH* chemiluminescence were visualized qualitatively. In the autoignition regime for both tribrachial structure and mild combustion, formaldehyde were found mainly between the fuel nozzle and the lifted flame edge. On the other hand, they were formed just prior to the flame edge for the non-autoignited lifted flames. The effect of fuel pyrolysis and partial oxidation were found to be important in explaining autoignited liftoff heights, especially in the Mild combustion regime. Flame structures of autoignited flames were investigated numerically for syngas (CO/H2) and methane fuels. The simulations of syngas fuel accounting for the differential diffusion have been performed by adopting several kinetic mechanisms to test the models ability in predicting the flame behaviors observed previously. The results agreed well with the observed nozzle-attached flame characteristics in case of non-autoignited flames. For autoignited lifted flames in high temperature regime, a unique autoignition behavior can be predicted having HO2 and H2O2 radicals in a broad region between the nozzle and stabilized lifted flame edge. Autoignition characteristics of laminar nonpremixed methane jet flames in high- temperature coflow air were studied numerically. Several flame configurations were investigated by varying the initial temperature and fuel mole fraction. Characteristics of chemical kinetics structures for autoignited lifted flames were discussed based on the kinetic structures of homogeneous autoignition and flame propagation of premixed mixtures. Results showed that for autoignited lifted flame with tribrachial structure, a transition from autoignition to flame propagation modes occurs for reasonably stoichiometric mixtures. Characteristics of Mild combustion can be treated as an autoignited lean premixed lifted flame. Transition behavior from Mild combustion to a nozzle-attached flame was also investigated by increasing the fuel mole fraction.

Electrically assisted combustion has been studied in order to control or improve flame characteristics, and emphasizing efficiency and emission regulation. Many phenomenological observations have been reported on the positive impact of electric fields on flame, however there is a lack of detailed physical mechanisms for interpreting these. To clarify the effects of electric fields on flame, I have investigated flame structure, soot formation, and flow field with ionic wind electrical current responses in nonpremixed counterflow flames. The effects of direct current (DC) electric field on flame movement and flow field was also demonstrated in premixed Bunsen flames. When a DC electric field was applied to a lower nozzle, the flames moved toward the cathode side due to Lorentz force action on the positive ions, soot particles simultaneously disappeared completely and laser diagnostics was used to identify the results from the soot particles. To understand the effects of an electric field on flames, flow visualization was performed by Mie scattering to check the ionic wind effect, which is considered to play an important role in electric field assisted combustion. Results showed a bidirectional ionic wind, with a double-stagnant flow configuration, which blew from the flame (ionic source) toward both the cathode and the anode. This implies that the electric field affects strain rate and the axial location of stoichiometry, important factors in maintaining nonpremixed counterflow flames; thus, soot formation of the counterflow flame can also be affected by the electric field. In a test of premixed Bunsen flames having parallel electrodes, flame movement toward the cathode and bidirectional ionic wind were observed. Using PIV measurement it was found that a created radial velocity caused by positive ions (i.e. toward a cathode), was much faster than the velocity toward the anode. Even in a study of alternating current (AC) electric fields, bidirectional ionic wind could be observed, regardless of applied frequencies. Therefore, the effect of ionic wind cannot be considered negligible under both DC and AC electric fields. Detailed explanations for electrical current, flame behavior, and flow characteristics under various conditions are discussed herein.

This work investigates the influence of one or more OH groups present on the fuel molecule and the resultant formation of NOX emissions. Combustion of oxygenated fuels has been increasing globally and such fuels offer significant potential in the reduction of pollutant emissions. One such emission class is the oxides of nitrogen, which typically form through a combination of two regimes: the thermal and non-thermal mechanisms. While thermal NO formation can be reduced by lowering the combustion temperature, non-thermal NO formation is coupled to the fuel chemistry. An experimental and computational investigation of NOX formation in three different burner configurations and under a range of equivalence ratios and temperature regimes explored the differences in NO formation. Measurements of temperature profiles and in-flame species concentrations, utilizing both probed and non-intrusive laser based techniques, allowed for the investigation of NO formation through non-thermal pathways and the differences that exist between fuels with varying numbers of OH groups. The first burner configuration was composed of a high swirl liquid spray burner with insulted combustion chamber walls designed specifically for the combustion of low energy density fuels. In this system the combustion of alcohols and glycerol (the largest by-product of biodiesel production), along with other fuels with multiple hydroxyl groups, was studied. Measurements of the mean flame temperature and exhaust gas measurements of NOX showed significant reductions in non-thermal NO concentrations with increasing numbers of OH groups. An accompanying modeling study and detailed reaction path analysis showed that fuel decomposition pathways through formaldehyde were shown a preference due to the presence of the OH groups which resulted in reduced contributions to the hydrocarbon radical pools subsequent reductions to the Prompt NO mechanism. Two burner configurations with reduced dimensionality facilitated measurements in premixed flames for temperature and species in high and low temperature flames. These measurements included probed thermocouple temperature measurements, extractive gas sampling for NO and intermediate hydrocarbon species, and planar Laser Induced Fluorescence (LIF) measurements for 2OH-LIF thermometry, semiquantitative CH2O LIF, and quantitative NO LIF. Additionally, the simplified nature of the burner geometries allowed for the modeling of the flames incorporating detailed reaction kinetics for fuel decomposition and NOX formation. Significant reductions in NO formation were observed in comparisons of alcohol and alkane flames, resulting in up to 50% reductions in the pollutant. Computational analyses and nitrogen flux accounting allowed for the identification of the reduction in NO formation through all the known NOX formation pathways. It was observed that all of the known pathways exhibited reductions in contributions to NO formation in the presence of OH functional groups, indicating a complex coupling of fuel and NOX chemistry.