Engineering PapersSearch

SEARCH · Engineering Papers

Results for “ELECTRON IMPACT”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Electron-impact excitation data for W 2+ in support of tungsten spectroscopy and re-deposition measurements for magnetically-confined plasmas

Abstract To better understand plasma wall interactions involving tungsten, accurate atomic structure and electron-impact driven collisional processes for near-neutral ion stages of tungsten are required. Complementing existing work on neutral and singly ionised tungsten, atomic structure and collisional calculations for W 2+ electron-impact excitation have been completed. These excitation calculations are an important component of S/XB coefficients for near-neutral charge states, which may be used to spectroscopically infer re-deposition of tungsten at the plasma-solid boundary of fusion relevant devices. With W 2+ in particular having emission lines that can be observed at ultraviolet (UV) wavelengths, while higher charge states of tungsten are unlikely to have lines possible to observe outside of the vacuum UV range. The atomic structure was generated using the General-purpose Relativistic Atomic Structure Package (GRASP 0 ), implementing the Multi-configuration Dirac Fock approach. This structure was the basis for a subsequent Dirac R -matrix electron-impact excitation calculation to provide Maxwellian averaged rate coefficients. A synthetic spectrum was generated from this data using a collisional-radiative model to predict the strongest W III spectral lines and these lines were compared to emission from the Compact Toroidal Hybrid (CTH) plasma device. Several of the strongest W III lines are observed in CTH and agree well with the modelled line wavelengths and intensities, a table of these lines is provided that could be observed in other devices.

McCann, M. (ORCID:0000000215321240)

Electron-impact ionization from excited states method

We present a distorted-wave with exchange (DWE) formulation for computing electron-impact ionization (EII) cross sections within the recently developed excited states method (ESM) for plasmas. The ESM provides a self-consistent quantum-mechanical description of both bound and continuum electronic states in dense plasmas, incorporating finite-temperature and screening effects absent in isolated-atom approaches. Using this framework, we calculate EII cross sections for lithium across a wide range of temperatures and densities. We show that the ESM reproduces the isolated-atom limit at low densities while capturing strong plasma effects, such as pressure ionization and shape resonances, at solid density. Comparisons with isolated-atom DWE and average-atom (AA) calculations reveal that, although AA cross sections can approximate configuration average behavior, the ESM provides more accurate, state-resolved answers. In conclusion, these results demonstrate that the ESM allows one to obtain rate coefficients for collisional-radiative modeling of plasmas out of local thermodynamical equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Electron-impact excitation of W VII

an adf04 file for Tungsten 6+ (W VII) resulting from a Dirac R-Matrix Electron Impact Excitation calculation.

McElroy, Niall [Queen's University Belfast] (ORCID

Characterization of runaway electron impact on instrumented sacrificial limiters on DIII-D

Instrumented sacrificial limiter heads, both domed (proud) and flat (flush) are used in DIII-D runaway electron (RE) wall strikes to study the wall impact dynamics with high spatial and time resolution. The approximate structure of the RE wetted area and heating depth on the domed limiter heads were predicted qualitatively using orbit-tracking simulations, although a strong left–right asymmetry (about the magnetic field direction) was not captured well by the simulations. It is hypothesized that this difference is perhaps due to the local 3D magnetic field perturbation of the dome limiter head. The average kinetic energy K and pitch angle θ of REs striking the limiter head were estimated from the spatial distribution of local HXR emission and were estimated to be roughly K ≈ 4 MeV and θ ≈ 0.2 rad. These values are roughly consistent with in-plasma values estimated before the loss event, indicating that RE kinetic energy and pitch angle are not drastically altered when transporting to the wall. Large shot–shot variations (1–10 kJ) in energy deposition into the limiter head were observed and were explained by shot–shot variations in locked magneto-hydrodynamics mode toroidal phase. For the largest deposited energies (10 kJ), graphite material failure and explosive dust release was observed, and the depth of material failure at higher energy deposition was successfully reproduced using modelling of volumetric energy deposition and brittle failure. The presence of energetic (keV) level ion impact during the RE wall strike was confirmed by three different surface analysis techniques. The ratio of energetic ion to RE flux appears to be larger on flat surfaces within the RE wetted area, although the energetic ion flux and total energy flux due to energetic ions have not yet been quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Understanding cold electron impact on parallel-propagating whistler chorus waves via moment-based quasilinear theory

Earth's magnetosphere hosts a wide range of collisionless particle populations that interact through various wave-particle processes. Among these, cold electrons, with energies below 100 eV, often dominate the plasma density but remain poorly characterized due to measurement challenges such as spacecraft charging and photoelectron contamination. Understanding the contribution of these cold populations to wave–particle interaction is of significant interest. Recent kinetic simulations identified a secondary drift-driven instability, in which parallel-propagating whistler-mode chorus waves excite oblique electrostatic whistler waves near the resonance cone and Bernstein-mode turbulence. These secondary modes enable a new channel of energy transfer from the parallel-propagating whistler wave to the cold electrons. In this work, we develop a moment-based quasilinear theory of the secondary instabilities to quantify such energy exchange. Our results show that these secondary instabilities persist for a wide range of parameters and, in many cases, lead to nearly complete damping of the primary wave. Such secondary instability might limit the amplitude of parallel-propagating whistler waves in Earth's magnetosphere and might explain why high-amplitude oblique whistler or electron Bernstein waves are rarely observed simultaneously with high-amplitude field-aligned whistler waves in the inner magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Unexpected Hydricity of a Bis-Carbene Iridium Complex Provides Insight into Electronic Structure Impacts on Hydride Donor Ability

The hydricity (ΔG° H– ), or hydride donor ability, of a transition metal complex is a thermodynamic parameter which can aid in the design and interpretation of various catalytic reactions that involve hydride transfer as a key step. In an attempt to generate a strong hydride donor, the bis-carbene ligand 3,3′-methylenebis(1-methyl-imidazol-2-ylidene) (“bis-mim”) was installed in an iridium hydride complex, [Cp*Ir(bis-mim)H] + . Experimental and computational studies show that [Cp*Ir(bis-mim)H] + is actually a relatively weak hydride donor, however. To understand why the complex is an unexpectedly weak hydride donor, experimental and computational studies probing the steric and electronic effects on hydricity were conducted. Steric factors had a minimal impact on thermodynamic hydricity but dampened kinetic hydricity. In conclusion, the poor thermodynamic hydride donor ability can be attributed to an electronic structure that results in relatively long Cp*−Ir bonds, an unusually high Ir–H BDFE and pK a values that have an outsized influence on thermochemical cycles for hydricity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS

Balmer emission measurements in JET-ILW hydrogen, deuterium, tritium and deuterium-tritium low-confinement mode plasmas

Measurements of Balmer-alpha (B-α) and Balmer-gamma (B-γ), where B = H, D, or T isotopes of hydrogen, emission across the low-field side (LFS) divertor in JET ITER-like Wall (ILW), low-confinement mode (L-mode) plasma indicate that the divertor plasma was more strongly detached in deuterium, tritium and deuterium–tritium than in hydrogen plasmas. These measurements are consistent with the lower degree of detachment in hydrogen versus tritium plasmas inferred from Langmuir probes (Groth et al., 2023). The ratio of B-γ to B-α is a suitable indicator of high density/low temperature conditions because B-α emission is dominated by electron impact excitation with a small electron–ion recombination component regardless of electron temperature and density (for conditions in tokamak plasmas) whilst for B-γ the recombination component increases and the excitation component decreases as the electron temperature decreases and the electron density increases. Within the uncertainty of the measurements, the electron density at the low-field side separatrix and the power crossing the separatrix are similar to permit these comparisons. The B-©/B-α ratio is well correlated with the onset of detachment as measured by Langmuir probes (McCracken et al., 1998; Groth et al., 2023). The ratios of B-γ and B-α varied between 0.005 and 0.035 for low-recycling, high-recycling and partially detached conditions, and were 60–70 % higher in partially detached hydrogen than in tritium plasmas. In contrast, the difference between deuterium, tritium and 40–60 % deuterium–tritium is within 10 % and thus within the uncertainties of the measurements. Profiles of the B-α and B-γ emission was measured with two filtered CCDs connected to the visible lost light leg of JET’s mirror link spectroscopy system (Meigs et al., 2010; Lomanowski et al., 2014). Spatially (radial) narrower B-α profiles were measured with increasing isotope mass, while similarly broadened B-γ profiles were measured.

Balmer

An upper pressure limit for low- Z benign termination of runaway electron beams in TCV

Abstract We present a model for the particle balance in the post-disruption runaway electron plateau phase of a tokamak discharge. The model is constructed with the help of, and applied to, experimental data from TCV discharges investigating the so-called ‘low- Z benign termination’ runaway electron mitigation scheme. In the benign termination scheme, the free electron density is first reduced in order for a subsequently induced MHD instability to grow rapidly and spread the runaway electrons widely across the wall. We show that the observed non-monotonic dependence of the free electron density with the measured neutral pressure is due to plasma re-ionization induced by runaway electron impact ionization. At higher neutral pressures, more target particles are present in the plasma for runaway electrons to collide with and ionize. Parameter scans are conducted to clarify the role of the runaway electron density and energy on the free electron density, and it is found that only the runaway electron density has a noticeable impact. While the free electron density is shown to be related to the spread of heat fluxes at termination, the exact cause for the upper neutral pressure limit remains undetermined and an object for further study.

Hoppe, M. (ORCID:0000000339948977)

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis

Characterization of uncertainties in electron-argon collision cross sections

Abstract The predictive capability of a plasma discharge model depends on accurate representations of electron-impact collision cross sections, which determine the corresponding reaction rates and electron transport properties. The values of cross sections can be known only approximately either through experiments or simulations and are thus subject to uncertainties. Quantifying the uncertainties in plasma simulations allows us to assess the reliability of simulations and to provide a basis for interpreting discrepancies between simulations and experiments. For such uncertainty quantification of plasma simulations, it is essential to quantify the uncertainties of the underlying cross sections. Although much effort has been committed to calibrate the cross section values, their uncertainties are not well investigated. We characterize uncertainties in electron-argon atom collision cross sections using a Bayesian framework. Six collision processes—elastic momentum transfer, ionization, and four excitations—are characterized with semi-empirical models, which effectively capture the features important to the macroscopic properties of the plasma. A probability model for the uncertain parameters of these semi-empirical models is developed. Specifically, a Gaussian-process likelihood model is proposed to capture discrepancies among data sets, as well as the model-form inadequacies of the semi-empirical models. Two other likelihood models are compared with the proposed Gaussian-process model, to illustrate the importance of the choice of the likelihood model. The cross section models are calibrated using the electron-beam experiments and ab-inito quantum simulations. The resulting calibrated uncertainties capture well the scattering among the data sets. The calibrated cross section models are further validated against swarm-parameter experiments and zero-dimensional Boltzmann equation simulations of widely used cross section datasets.

Chung, Seung Whan (ORCID:0000000302501549)

Benchmark of the Fe xvv 𝓡 ratio in photoionized plasma during eclipse of Centaurus X-3 with XRISM/Resolve

The $\mathcal {R}$ ratio is a useful diagnostic of the X-ray emitting astrophysical plasmas and is defined as the intensity ratio of the forbidden over the inter-combination lines in the K$\alpha$ line complex of He-like ions. The value is altered by excitation processes (electron impact or UV photoexcitation) from the metastable upper level of the forbidden line, thereby constraining the electron density or UV field intensity. The diagnostic has been applied mostly in electron density constraints in collisionally ionized plasmas using low-Z elements, as was originally proposed for the Sun (Gabriel & Jordan, 1969a, MNRAS, 145, 241), but it can also be used in photoionized plasmas. To make use of this diagnostic, we need to know its value in the limit of no excitation of metastables ($\mathcal {R}_{0}$), which depends on the element, how the plasmas are formed, how the lines are propagated, and the spectral resolution affecting line blending principally with satellite lines from Li-like ions. We benchmark $\mathcal {R}_0$ for photoionized plasmas by comparing calculations using radiative transfer codes and observation data taken with the Resolve X-ray microcalorimeter onboard XRISM. We use the Fe xxv He$\alpha$ line complex of the photo-ionized plasma in Centaurus X-3 observed during eclipse, in which the plasma is expected to be in the limit of no metastable excitation. The measured $\mathcal {R} = 0.65 \pm 0.08$ is consistent with the value calculated using xstar for the plasma parameters derived from other line ratios of the spectrum. We conclude that the $\mathcal {R}$ ratio diagnostic can be used for high-Z elements such as Fe in photoionized plasmas, which has wide applications in plasmas around compact objects at various scales.

X-rays: binaries

Experimental validation of a collision-radiation dataset for molecular hydrogen in plasmas

Quantitative spectroscopy of molecular hydrogen has generated substantial demand, leading to the accumulation of diverse elementary process data encompassing radiative transitions, electron-impact transitions, predissociations, and quenching. However, their rates currently available are still sparse, and there are inconsistencies among those proposed by different authors. In this study, we demonstrate an experimental validation of such a molecular dataset by composing a collisional-radiative model (CRM) for molecular hydrogen and comparing experimentally obtained vibronic populations across multiple levels. From the population kinetics of molecular hydrogen, the importance of each elementary process in various parameter space is studied. In low-density plasmas (electron density ne≲1017 m−3) the excitation rates from the ground states and radiative decay rates, both of which have been reported previously, determine the excited state population. The inconsistency in the excitation rates affects the population distribution the most significantly in this parameter space. However, in higher density plasmas (ne≳1018 m−3), the excitation rates from excited states become important, which have never been reported in the literature, and may need to be approximated in some way. In order to validate these molecular datasets and approximated rates, we carried out experimental observations for two different hydrogen plasmas; a low-density radio frequency heated plasma (ne≈1016 m−3) and the Large Helical Device (LHD) divertor plasma (ne≳1018 m−3). The visible emission lines from EF1Σg+, HH¯1Σg+, D1Πu±, GK1Σg+, I1Πg±, J1Δg±, h3Σg+, e3Σu+, d3Πu±,g3Σg+, i3Πg±, and j3Δg± states were observed simultaneously and their population distributions were obtained from their intensities. We compared the observed population distributions with the CRM prediction, in particular the CRM with the rates compiled by Janev et al., Miles et al., and those calculated with the molecular convergent close-coupling (MCCC) method. The MCCC prediction gives the best agreement with the experiment, particularly for the emission from the low-density plasma. However, the population distribution in the LHD divertor shows a worse agreement with the CRM than those from low-density plasma, indicating the necessity of the precise excitation rates from excited states. We also found that the rates for the electron attachment is inconsistent with experimental results. This requires further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Nonlocal, diamagnetic electromagnetic effects in magnetically insulated transmission lines

We identify the time-dependent physics responsible for the critical reduction of current losses in magnetically insulated transmission lines (MITLs) due to uninsulated space charge-limited currents of electrons emitted by field stress. A drive current of sufficiently short pulse length introduces a strong enough time dependence that steady-state results alone become inadequate for the complete understanding of current losses. The time-dependent physics can be described as a nonlocal, diamagnetic electromagnetic response of space charge limited currents. As the pulse length is increased or equivalently, the MITL length reduced, these time-dependent effects diminish and current losses converge to those predicted by the well-known Child–Langmuir law in the external (vacuum) fields. We present a simple one-dimensional (1D) model that encapsulates the essence of this physics. We find excellent agreement with 2D particle-in-cell simulations for two MITL geometries, Cartesian parallel plate and azimuthally symmetric straight coaxial. Based on the 1D model, we explore various scaling dependencies of MITL losses with relevant parameters, e.g., peak current, pulse length, geometrical dimensions, etc. We propose an improved physics model of magnetic insulation in the form of a Hull curve, which could also help improve predictions of current losses by common circuit element codes, such as BERTHA. Finally, we describe how to calculate the temperature rise due to electron impact within the 1D model.

Computer simulation

Spatial profile of argon (1s 5 ) metastables in an electron beam generated plasma

Electron beams with an applied magnetic field generate a secondary cold plasma with a selective chemical composition, featuring low-energy ions and metastable species in the discharge periphery, ideal for low-damage plasma treatment of material substrates. In this work, we studied the plasma generated by an e-beam using a 4 kV voltage in a pure argon gas environment under a magnetic field of 150 G and in the pressure range of 25–90 mTorr. We measured the absolute spatial density profile of argon (1 s 5 ) metastables in an electron beam generated plasma by laser-induced fluorescence and found it to be of the order of 10 16 m −3 . The electron temperature and the electron density measured by a Langmuir probe were of the order of 10 16 m −3 and less than an eV respectively. Electron-impact quenching was identified as a significant loss mechanism for the Ar(1s 5 ) state, leading to the saturation of the metastable density at higher pressures. Outside the primary ionization region, the spatial distribution of argon metastables followed a linear diffusion profile, indicating negligible additional production in those regions.

EEDF

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN

Mechanistic Insights for Plasma-Catalytic CO 2 Reduction over TiO 2 in a Dielectric Barrier Discharge Reactor

Reaction kinetics experiments coupled with phenomenological kinetic modeling and parameter estimation are used to elicit insights into the mechanism and active sites for the plasma-catalytic dissociation of CO 2 on TiO 2 . Experimental and model insights showed that gas-phase reactions contribute at least two-thirds of the overall product formation at explored conditions; weak temperature dependence, strong sensitivity to specific energy input (SEI), apparent first order in CO 2 , and positive influence of cofed argon (Ar) and oxygen (O 2 ) for the gas-phase contributions all suggest that expected plasma reaction steps such as electron-impact and high-energy collisions are the dominant modes for CO 2 dissociation. The Arrhenius-like expression for gas contributions resulted in a preexponential of 4.40 × 10 –3 s –1 , an E SEI,g of 7.90 × 10 –4 mol/kJ, and an E a,g of 1.00 × 10 –3 J/mol. For surface contributions, the small apparent barrier of 16.3 kJ/mol, relatively weaker dependence on SEI, first-order dependence on CO 2 , and insensitivity to cofed Ar and O 2 all point to CO 2 dissociation on TiO 2 surface facets without vacancies and aided by plasma (leading to vibrationally excited CO 2 and/or a reactive surface with significant surface charge accumulation). The Arrhenius-like expression resulted in a preexponential of 7.81 × 10 –2 s –1 , an E SEI,s of 1.90 × 10 –3 mol/kJ, and an E a,s of 1.63 × 10 4 J/mol. The derived kinetic model further enabled a systematic evaluation of the effect of inputs (plasma power, flow rate, CO 2 inlet concentration, and temperature) to identify process trends and optimal operating conditions.

catalyst

Experimental studies of high-temperature thermal dissociation of iso-propanol

iso-Propanol is an important biofuel for transportation and industrial applications, as well as a canonical compound for understanding alcohol reaction chemistry. There are few studies of the thermal decomposition of iso-propanol, hence the objective of this work was to investigate the mechanism of iso-propanol pyrolysis at high-temperature conditions. Shock tube studies were used to investigate the thermal decomposition of iso-propanol and identify the dominant reaction pathways and the intermediate species present in the temperature range of 1395–2053 K and at pressures of 0.3–4 bar. Time-resolved mass spectra were obtained from pyrolysis experiments in shock tube/mass spectrometer apparatuses using electron impact ionization or synchrotron vacuum ultraviolet photoionization. The relative concentrations of the intermediate species were used to identify the dominant reaction channels. iso-Propanol dissociated by water loss to propene and by elimination of a methyl radical yielding a 1-hydroxyethyl radical at the conditions studied. The time-dependent mass spectra allowed secondary products to be identified along with the orders of appearance. Notably, peaks at m/z = 50 and 52 were observed providing insight into secondary reactions. These features have not been previously reported in propanol-pyrolysis literature.

Kim, John H. [Univ. of Michigan, Ann Arbor, MI (Un