The separation of atomic and molecular gases from helium by cataphoresis

The cataphoretic purification of helium was investigated for binary mixtures of He with Ar, Ne, N₂, O₂, CO, and CO₂ in DC glow discharge. An experimental technique was developed to continuously measure the composition in the anode end-bulb without sample withdrawal. Discharge currents ranged from 10 ma to 100 ma. Total gas pressure ranged from 2 torr to 9 torr. Initial compositions of the minority component in He ranged from 1.2 mole percent to 7.5 mole percent.

The cataphoretic separation of Ar and Ne from He was found to be in agreement with previous investigators. The cataphoretic separation of N₂, O₂, and CO from He was found to be similar to noble gas systems in that the steady-state separation improved with (1) increasing discharge current, (2) increasing gas pressure, and (3) decreasing initial composition of the minority component. In the He-CO₂ mixture, the CO₂ dissociated to CO plus O₂. The fraction of CO₂ dissociated was directly proportional to the current and pressure and independent of initial composition.

The experimental results for the separation of Ar, Ne, N₂, O₂, and CO from He were interpreted in the framework of a recently proposed theoretical model involving an electrostatic Peclet number. In the model the electric field was assumed to be constant. This assumption was checked experimentally and the maximum variation in electric field was 35% in time and 30% in position. Consequently, the assumption of constant electric field introduced no more than 55% variation in the electrostatic Peclet number during a separation.

To aid in the design of new cataphoretic systems, the following design criteria were developed and tested in detail: (1) electric field independent of discharge current, (2) electric field directly proportional to total pressure, (3) ion fraction of impurity directly proportional to discharge current, and (4) ion fraction of impurity independent of total pressure. Although these assumptions are approximate, they enabled the steady-state concentration profile to be predicted to within 25% for 75% of the data. The theoretical model was also tested with respect to the characteristic time associated with transient cataphoresis. Over 80% of the data was within a factor of two of the calculated characteristic times.

The electrostatic Peclet number ranged in value from 0.13 to 4.33. Back-calculated ion fractions of the impurity component ranged in value from 4.8x10^-6 to 178x10^-6.

Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

Small scale turbulence in high Karlovitz number premixed flames

The purpose of this thesis is to characterize the behavior of the smallest turbulent scales in high Karlovitz number (Ka) premixed flames. These scales are particularly important in the two-way coupling between turbulence and chemistry and better understanding of these scales will support future modeling efforts using large eddy simulations (LES). The smallest turbulent scales are studied by considering the vorticity vector, ω, and its transport equation.

Due to the complexity of turbulent combustion introduced by the wide range of length and time scales, the two-dimensional vortex-flame interaction is first studied as a simplified test case. Numerical and analytical techniques are used to discern the dominate transport terms and their effects on vorticity based on the initial size and strength of the vortex. This description of the effects of the flame on a vortex provides a foundation for investigating vorticity in turbulent combustion.

Subsequently, enstrophy, ω² = ω • ω, and its transport equation are investigated in premixed turbulent combustion. For this purpose, a series of direct numerical simulations (DNS) of premixed n-heptane/air flames are performed, the conditions of which span a wide range of unburnt Karlovitz numbers and turbulent Reynolds numbers. Theoretical scaling analysis along with the DNS results support that, at high Karlovitz number, enstrophy transport is controlled by the viscous dissipation and vortex stretching/production terms. As a result, vorticity scales throughout the flame with the inverse of the Kolmogorov time scale, τ_η, just as in homogeneous isotropic turbulence. As τ_η is only a function of the viscosity and dissipation rate, this supports the validity of Kolmogorov’s first similarity hypothesis for sufficiently high Ka numbers (Ka ≳ 100). These conclusions are in contrast to low Karlovitz number behavior, where dilatation and baroclinic torque have a significant impact on vorticity within the flame. Results are unaffected by the transport model, chemical model, turbulent Reynolds number, and lastly the physical configuration.

Next, the isotropy of vorticity is assessed. It is found that given a sufficiently large value of the Karlovitz number (Ka ≳ 100) the vorticity is isotropic. At lower Karlovitz numbers, anisotropy develops due to the effects of the flame on the vortex stretching/production term. In this case, the local dynamics of vorticity in the strain-rate tensor, S, eigenframe are altered by the flame. At sufficiently high Karlovitz numbers, the dynamics of vorticity in this eigenframe resemble that of homogeneous isotropic turbulence.

Combined, the results of this thesis support that both the magnitude and orientation of vorticity resemble the behavior of homogeneous isotropic turbulence, given a sufficiently high Karlovitz number (Ka ≳ 100). This supports the validity of Kolmogorov’s first similarity hypothesis and the hypothesis of local isotropy under these condition. However, dramatically different behavior is found at lower Karlovitz numbers. These conclusions provides/suggests directions for modeling high Karlovitz number premixed flames using LES. With more accurate models, the design of aircraft combustors and other combustion based devices may better mitigate the detrimental effects of combustion, from reducing CO₂ and soot production to increasing engine efficiency.

I. Nickel-Exchanged Zincosilicate Catalysts for the Oligomerization of Propylene and II. Organic SDA-Free Catalysts for the Methanol-to-Olefins Reaction

Nickel-containing catalysts are developed to oligomerize light olefins. Two nickel-containing zincosilicates (Ni-CIT-6 and Ni-Zn-MCM-41) and two nickel-containing aluminosilicates (Ni-HiAl-BEA and Ni-USY) are synthesized as catalysts to oligomerize propylene into C_3n (C₆ and C₉) products. All catalysts oligomerize propylene, with the zincosilicates demonstrating higher average selectivities to C_3n products, likely due to the reduced acidity of the Zn heteroatom.

To test whether light alkanes can be incorporated into this oligomerization reaction, a supported homogeneous catalyst is combined with Ni-containing zincosilicates. The homogeneous catalyst is included to provide dehydrogenation/hydrogenation functions. When this tandem catalyst system is evaluated using a propylene/n-butane feed, no significant integration of alkanes are observed.

Ni-containing zincosilicates are reacted with 1-butene and an equimolar propylene/1-butene mixture to study other olefinic feeds. Further, other divalent metal cations such as Mn²⁺, Co²⁺, Cu²⁺, and Zn²⁺ are exchanged onto CIT-6 samples to investigate stability and potential use for other reactions. Co-CIT-6 oligomerizes propylene, albeit less effectively than Ni-CIT-6. The other M-CIT-6 samples, while not able to oligomerize light olefins, may be useful for other reactions, such as deNO_x.

Molecular sieves are synthesized, characterized, and used to catalyze the methanol-to-olefins (MTO) reaction. The Al concentration in SSZ-13 samples is varied to investigate the effect of Al number on MTO reactivity when compared to a SAPO-34 sample with only isolated Si Brønsted acid sites. These SSZ-13 samples display reduced transient selectivity behavior and extended reaction lifetimes as Si/Al increases; attributable to fewer paired Al sites. MTO reactivity for the higher Si/Al SSZ-13s resembles the SAPO-34 sample, suggesting that both catalysts owe their stable reaction behavior to isolated Brønsted acid sites.

Zeolites CHA and RHO are prepared without the use of organic structure-directing agents (OSDAs), dealuminated by steam treatments (500°C-800°C), and evaluated as catalysts for the MTO reaction. The effects of temperature and steam partial pressure during steaming are investigated. X-ray diffraction (XRD) and Ar physisorption show that steaming causes partial structural collapse of the zeolite, with degradation increasing with steaming temperature. ²⁷Al MAS NMR spectra of steamed materials reveal the presence of tetrahedral, pentacoordinate, and hexacoordinate aluminum.

Proton forms of as-synthesized CHA (Si/Al=2.4) and RHO (Si/Al=2.8) rapidly deactivate under MTO testing conditions (400°C, atmospheric pressure). CHA samples steamed at 600°C performed best among samples tested, showing increased olefin selectivities and catalyst lifetime. Acid washing these steamed samples further improved activity. Reaction results for RHO were similar to CHA, with the RHO sample steamed at 800°C producing the highest light olefin selectivities. Catalyst lifetime and C₂-C₃ olefin selectivities increase with increasing reaction temperature for both CHA-type and RHO-type steamed samples.

The electrocatalytic and stoichiometric oxidation chemistry of binuclear ruthenium complexes incorporating the anionic tripod ligand {(η5-C5H5)Co[P(OCH3)2(=O)]3}-

The anionic tripod ligand NaLoMe (L_(oMe) - = [(η^5-C_5H_5)Co{P(O)(OCH_3)_2}_3]^-) reacts with RuO_4 in a biphasic reaction mixture of 1% H_2SO_4 and CCI_4 to afford [(L_(oMe) (HO)Ru^(IV) (µ-O)_2Ru ^(IV)(OH)(L_(oMe)] (1), which is treated with aqueous CF_3S0_3H to generate [(L_(oMe)(H_2O)Ru^(IV) (µ-O)_2R^(IV) (OH_2)(L_(oMe)][CF_3SO_3]_2 ([H_21][CF_3SO_3]_2). Addition of iodosobenzene to an acetonitrile solution of this salt yields [(L_(oMe)(O)Ru^v(µ-0)2Ru^v-(O)(_(LoMe)] (2). The dimer 1 can be reduced chemically or electrochemically to the Ru^(III)- Ru^(III) dimers [(L_(oMe)(H_20)Ru^(III) (µ-OH)_2Ru^(III) (OH_2)(L_(oMe)) ]^2+ and [(L_(oMe)) ^(III) (µ-0Hh(µ-0H2)Ru^(III) (L_(oMe)]^2+ which interconvert in aqueous media. Two electron processes dominate both the bulk chemistry and the electrochemistry of 1. Among these processes are the quasi-reversible Ru^(IV) - Ru^(IV)/Ru^(III)- Ru^(III) and Ru^(III)- Ru^(III)/ Ru^(II)- Ru^(II) reductions and a largely irreversible Ru^(V) - Ru^(V)/ Ru^(IV) - Ru^(IV)/oxidation. The dioxo dimer 2 oxidizes alcohols and aldehydes in organic media to afford 1 and the corresponding aldehydes and acids. Analogously, the Ru^(V) - Ru^(V)/ Ru^(IV)- Ru^(IV) redox wave mediates the electrooxidation of alcohols and aldehydes in aqueous buffer. In this system, substrates can be oxidized completely to CO_2. The kinetic behavior of these oxidations was examined by UV-vis and chronoamperometry, respectively, and the chemistry is typical of metal-oxo complexes, indicating that electronic coupling between two metal centers does not dramatically affect the metal-oxo chemistry. Dimer [H_21]^(2+) also reacts with alcohols, aldehydes, and triphenylphosphine in CH_3CN to afford Ru^(III)- Ru^(III) products including [(L_(oMe))CH_3CN) Ru^(III) (µ-OH)_2 Ru^(III) (NCCH_3)( L_(oMe))][CF_3SO_3]2 (characterized by X-ray crystallography) and the corresponding organic products. Reaction of 1 with formaldehyde in aqueous buffer quantitatively affords the triply bridged dimer [(L_(oMe)Ru^(III) (µ-OH)2- (µ-HCOO) Ru^(III) (L_(oMe)][CF_3SO_3] (characterized by X-ray crystallography). This reaction evidently proceeds by two parallel inner-sphere pathways, one of which is autocatalytic. Neither pathway exhibits a primary isotope effect suggesting the rate determining process could be the formation of an intermediate, perhaps a Ru^(IV) - Ru^(IV) formate adduct. The Ru^(III)- Ru^(III)formate adduct is easily oxidized to the Ru^(IV) - Ru^(IV) analog [(L_(oMe)Ru^(IV)(µ-OH)_2-(µ-HCOO) Ru^(IV) (L_(oMe)][CF_3SO_3], which, after isolation, reacts slowly with aqueous formaldehyde to generate free formate and the Ru^(III)- Ru^(III) formate adduct. These dimers function as catalysts for the electrooxidation of formaldehyde at low anodic potentials (+0.0 V versus SCE in aqueous buffer, pH 8.5) and enhance the activity of Nafion treated palladium/carbon heterogeneous fuel cell catalysts.