Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient structure”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2026

Subseasonal Tropical Convection Characteristics in the Energy Exascale Earth System Model Version 2

Accurate simulation of subseasonal tropical moist convection remains a key challenge for Earth system models. The difficulties stem from the reliance of cumulus cloud processes on model parameterizations and the need to represent the multiscale nature of interactions among clouds, radiation, moisture, circulation, and surface energy fluxes. Equatorial convection drives circulation anomalies that can affect weather patterns and extremes globally, motivating efforts to better understand and simulate these tropical disturbances. Here, a detailed review of subseasonal tropical convective behavior as simulated in the Energy Exascale Earth System Model version 2 (E3SMv2) is presented, with comparison to its predecessor version 1 (E3SMv1) and reference data sets. Model structural changes to the deep convective trigger function and surface fluxes, along with parametric tuning of the cloud and microphysics schemes, together result in an improved depiction of organized tropical convection across scales. In particular, E3SMv2 exhibits a more realistic Madden‐Julian oscillation (MJO) and low‐frequency Kelvin waves—owing to a sharper time mean equatorial meridional moisture gradient and improved convection‐circulation coupling —as well as a better depiction of MJO Northern Hemisphere teleconnections. Despite these improvements, subseasonal precipitation variance continues to be strongly underestimated in E3SMv2. Use of a cloud plume model also reveals that the coupling between daily averaged tropical precipitation and lower tropospheric instability in E3SM is inconsistent with observations, a bias that could potentially impact the simulation of intraseasonal disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Machine-learning prediction for quasiparton distribution function matrix elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-x dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken x) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small- x reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_za)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements, especially in finer lattice units, increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning algorithms can make predictions of correlators to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different prediction accuracy and systematic errors. The predictions from smaller displacement z to larger ones work better than those for momentum p due to the higher correlation among the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Manipulating the diffusion energy barrier at the lithium metal electrolyte interface for dendrite-free long-life batteries

Abstract Constructing an artificial solid electrolyte interphase (SEI) on lithium metal electrodes is a promising approach to address the rampant growth of dangerous lithium morphologies (dendritic and dead Li 0 ) and low Coulombic efficiency that plague development of lithium metal batteries, but how Li + transport behavior in the SEI is coupled with mechanical properties remains unknown. We demonstrate here a facile and scalable solution-processed approach to form a Li 3 N-rich SEI with a phase-pure crystalline structure that minimizes the diffusion energy barrier of Li + across the SEI. Compared with a polycrystalline Li 3 N SEI obtained from conventional practice, the phase-pure/single crystalline Li 3 N-rich SEI constitutes an interphase of high mechanical strength and low Li + diffusion barrier. We elucidate the correlation among Li + transference number, diffusion behavior, concentration gradient, and the stability of the lithium metal electrode by integrating phase field simulations with experiments. We demonstrate improved reversibility and charge/discharge cycling behaviors for both symmetric cells and full lithium-metal batteries constructed with this Li 3 N-rich SEI. These studies may cast new insight into the design and engineering of an ideal artificial SEI for stable and high-performance lithium metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Autonomous phase mapping of gold nanoparticles synthesis with differentiable models of spectral shape

Autonomous experimentation–or self-driving labs–offers a systematic approach to accelerate materials discovery by integrating automated synthesis, characterization, and data-driven decision-making. We present a closed-loop workflow for the on-demand synthesis and structural characterization of colloidal gold nanoparticles, enabling direct mapping from composition to nanoscale structure. Our framework leverages differentiable models of spectral shape to address two central tasks in self-driving labs: (a) phase mapping, or identifying compositional regions with distinct structural behavior; and (b) material retrosynthesis, or optimizing compositions for target structure. Using functional data analysis, we develop a data-driven model with generative pre-training, active learning, and high-throughput experiments to predict spectral responses across composition space. We demonstrate the approach on seed-mediated growth of gold nanoparticles, showcasing its ability to extract design rules, reveal secondary interactions, and efficiently navigate morphology space. Gradient-based optimization of the models enables inverse design, making this a unified platform.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stochastic and deterministic processes shape bioenergy crop microbiomes along a vertical soil niche

Abstract Sustainable biofuel cropping systems aim to address climate change while meeting energy needs. Understanding how soil and plant‐associated microbes respond to these different cropping systems is key to promoting agriculture sustainability and evaluating changes in ecosystem functions. Here, we leverage a long‐term biofuel cropping system field experiment to dissect soil and root microbiome changes across a soil‐depth gradient in poplar, restored prairie and switchgrass to understand their effects on the microbial communities. High throughput amplicon sequencing of the fungal internal transcribed spacer (ITS) and prokaryotic 16S DNA regions showed a common trend of root and soil microbial community richness decreasing and evenness increasing with depth. Ecological niche (root vs. soil) had the strongest effect on community structure, followed by depth, then crop. Stochastic processes dominated the structuring of fungal communities in deeper soil layers while operational taxonomic units (OTUs) in surface soil layers were more likely to co‐occur and to be enriched by plant hosts. Prokaryotic communities were dispersal limited at deeper depths. Microbial networks showed a higher density, connectedness, average degree and module size in deeper soils. We observed a decrease in fungal‐fungal links and an increase of bacteria–bacteria links with increasing depth in all crops, particularly in the root microbiome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Fungal community structure and function shifts with atmospheric nitrogen deposition

Fungal decomposition of soil organic matter depends on soil nitrogen (N) availability. This ecosystem process is being jeopardized by changes in N inputs that have resulted from a tripling of atmospheric N deposition in the last century. Soil fungi are impacted by atmospheric N deposition due to higher N availability, as soils are acidified, or as micronutrients become increasingly limiting. Fungal communities that persist with chronic N deposition may be enriched with traits that enable them to tolerate environmental stress, which may trade-off with traits enabling organic matter decomposition. We hypothesized that fungal communities would respond to N deposition by shifting community composition and functional gene abundances toward those that tolerate stress but are weak decomposers. We sampled soils at seven eastern US hardwood forests where ambient N deposition varied from 3.2 to 12.6 kg N ha -1 year -1 , five of which also have experimental plots where atmospheric N deposition was simulated through fertilizer application treatments (25–50 kg N ha -1 year -1 ). Fungal community and functional responses to fertilizer varied across the ambient N deposition gradient. Fungal biomass and richness increased with simulated N deposition at sites with low ambient deposition and decreased at sites with high ambient deposition. Fungal functional genes involved in hydrolysis of organic matter increased with ambient N deposition while genes involved in oxidation of organic matter decreased. One of four genes involved in generalized abiotic stress tolerance increased with ambient N deposition. In summary, we found that the divergent response to simulated N deposition depended on ambient N deposition levels. Fungal biomass, richness, and oxidative enzyme potential were reduced by N deposition where ambient N deposition was high suggesting fungal communities were pushed beyond an environmental stress threshold. Fungal community structure and function responses to N enrichment depended on ambient N deposition at a regional scale.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2021

Thermoelectric current in a graphene Cooper pair splitter

Generation of electric voltage in a conductor by applying a temperature gradient is a fundamental phenomenon called the Seebeck effect. This effect and its inverse is widely exploited in diverse applications ranging from thermoelectric power generators to temperature sensing. Recently, a possibility of thermoelectricity arising from the interplay of the non-local Cooper pair splitting and the elastic co-tunneling in the hybrid normal metal-superconductor-normal metal structures was predicted. Here, we report the observation of the non-local Seebeck effect in a graphene-based Cooper pair splitting device comprising two quantum dots connected to an aluminum superconductor and present a theoretical description of this phenomenon. The observed non-local Seebeck effect offers an efficient tool for producing entangled electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Bolometer tomography on Wendelstein 7-X for study of radiation asymmetry

The algorithm for bolometer tomography at Wendelstein 7-X (W7-X) has been recently improved using a novel regularization functional, based on relative gradient smoothing (RGS) of the sought radiation profile. It has been validated using radiation patterns provided by 3D modeling under real plasma conditions as phantoms and then applied to bolometer measurements performed during the first divertor operation phase of W7-X. The following results are presented: 1) edge-localized 2D radiation patterns with clearly resolved magnetic island radiation structures, 2) an up-down asymmetry in the impurity radiation that is not captured by the 3D edge plasma transport modeling, 3) reversal of the asymmetry with reversed magnetic field direction. Further analysis reveals a poloidal variation of the emissivity in the outer confined plasma region with a field-direction dependent asymmetry, also supported by the soft X-ray measurements. This asymmetry is considered to be related to asymmetric impurity distributions, driven by the pronounced ion-impurity friction force at the plasma edge where the collisionality of the W7-X plasma is sufficiently high to develop impurity asymmetry as predicted by neoclassical theory of parallel impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Association between soil organic carbon and calcium in acidic grassland soils from Point Reyes National Seashore, CA

All data is collected from across 3 soil cores spanning a pH gradient at Point Reyes National Seashore, as part of a recent publication. The bulk data csv includes hygroscopic moisture content, soil pH, total C and N values, particle size distributions, cation exchange capacity and extractable cations. The STXM data package includes raw spectra data from scanning transmission X-ray microscopy C near-edge X-ray absorption fine structure spectroscopy analysis, extracted from STXM Image Reader analysis and then subsequently normalised in Athena (see details in paper). The data file also contains a bulk chemistry and mineralogy dataset, including total and trace element contents from 3 depths and mineral compositions as measured at BL 11-3 SSRL and semi-quantified in HighScore. Finally the final data file is calcium K-edge X-ray absorption near-edge structure spectroscopy (BL 4-3) and micro-XANES and micro-X-ray fluorescence (BL 14-3b) data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗
DOE OSTI2023

Advancements in High-Q development for novel medium-velocity 644 MHz 5-cell elliptical superconducting RF cavities for continuous-wave operation in heavy-ion linacs

The Facility for Rare Isotope Beams is a world-leading center for experimental nuclear physics research, and relies on a first-in-kind superconducting RF driver linac to supply 400 kW beams of a uniquely large range of particles, from protons through the heaviest uranium ions. The FRIB400 project proposes to double the end energy of the current superconducting FRIB linac for the heaviest uranium ions from 200 MeV per nucleon (MeV/u) to 400 MeV/u, which equates to approximately 1 GeV for protons. The increased rare isotope production from higher-energy drive beams would deliver new and exciting capabilities to FRIB users, which include extending the facility's reach along the neutron drip-line, greatly increasing the facility's rare isotope yield, and significantly improving the precision of parameter measurements for the nuclear matter equation of state, which is a particularly timely complement to the recent advent of multi-messenger astronomy, and the study of neutron star mergers. The unique set of design parameters defined by the 400 MeV/u (uranium) energy goal, the physical size, continuous-wave operation, and cryogenic capacity of the extant facilities at FRIB, requires development of a first-in-kind beta = v/c = 0.65, 644 MHz 5-cell elliptical superconducting RF (SRF) cavity, with field-leading efficiency, measured by the cavity's intrinsic quality factor, Q_0. The field of SRF has recently experienced a very active period of development, which demonstrated the potential to raise the Q_0 of 1.3 GHz cavities by factors of two to three or more. In addition to uncovering fascinating new properties of the fundamental physics of superconducting RF phenomena, these improvements suggest substantial improvements could be made over the minimum Q_0 of the proposed FRIB400 upgrade SRF cavities, if these new techniques can be adapted to the beta = 0.65, 644 MHz regime. This work presents the first validation of the the novel beta = 0.65, 644 MHz 5-cell elliptical superconducting RF cavity design, and the first results of both conventional and advanced RF surface preparation techniques. In particular, this work focuses on two recently-developed high-Q_0 recipes: N-doping, and furnace baking, setting the current world-record Q_0 for this type of resonant cavity at 3.8 x 10^10 at the FRIB400 operating gradient of 17.5 MV/m. We also find evidence that these results may be further improved upon, and make the case for further RF surface preparation recipe optimization. We then pair these high-power RF investigations with techniques borrowed from material science to try to better understand the ways in which specific properties of the niobium material, such as grain structure and superconducting flux pinning force, can be related to its superconducting RF performance. These experiments better elucidate the effects of high-temperature annealing on niobium samples, which suggest methods of further optimizing cavity flux expulsion and thus SRF performance.Read

McGee, Kellen E.↗
DOE OSTI2022

Ionic transport in electrostatic Janus Membranes. An explicit solvent molecular dynamic simulation.

Janus –or two-sided, charged membranes offer promise as ionic current rectifiers. In such systems, pores consisting of two regions of opposite charge can be used to generate a current from a gradient in salinity. The efficiency of nanoscale Janus pores increases dramatically as their diameter becomes smaller. However, little is known about the underlying transport processes, particularly under experimentally accessible conditions. In this work, we examine the molecular basis for rectification in Janus nanopores using an applied electric field. Molecular simulations with explicit water and ions are used to examine the structure and dynamics of all molecular species in aqueous electrolyte solutions. For several macroscopic observables, the results of such simulations are consistent with experimental observations on asymmetric membranes. Our analysis reveals a number of previously unknown features, including a pronounced local re-orientation of water molecules in the pores, and a segregation of ionic species that has not been anticipated by previously reported continuum analyses of Janus pores. Using these insights, a model is proposed for ionic current rectification in which electric leakage at pore entrance controls net transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Turbulence and transport in mirror geometries in the Large Plasma Device

Thanks to advances in plasma science and enabling technology, mirror machines are being reconsidered for fusion power plants and as possible fusion volumetric neutron sources. However, cross-field transport and turbulence in mirrors remains relatively understudied compared with toroidal devices. Turbulence and transport in mirror configurations were studied utilizing the flexible magnetic geometry of the Large Plasma Device (LAPD). Multiple mirror ratios from $M=1$ to $M=2.68$ and three mirror-cell lengths from $L=3.51$ to $L=10.86$ m were examined. Langmuir and magnetic probes were used to measure profiles of density, temperature, potential and magnetic field. The electric field-fluctuation-driven ${\tilde {\boldsymbol{E}}} \times {\boldsymbol{B}}$ particle flux, where $\boldsymbol{B}$ is the background field, was calculated from these quantities. Two probe correlation techniques were used to infer wavenumbers and two-dimensional structure. Cross-field particle flux and density fluctuation power decreased with increased mirror ratio. Core density and temperatures remain similar with mirror ratio, but radial line-integrated density increased. The physical expansion of the plasma in the mirror cell by using a higher field in the source region may have led to reduced density fluctuation power through the increased gradient scale length. This increased scale length reduced the growth rate and saturation level of rotational interchange and drift-like instabilities. Despite the introduction of magnetic curvature, no evidence of mirror-driven instabilities – interchange, velocity space or otherwise – were observed. For curvature-induced interchange, many possible stabilization mechanisms were present, suppressing the visibility of the instability.

Travis, Phil (ORCID:0000000343170432)↗
DOE OSTI2020

Generalized Grad–Shafranov equation for non-axisymmetric MHD equilibria

The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad–Shafranov equation, which is a single elliptic partial differential equation in two space dimensions. By revealing a hidden symmetry, here we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad–Shafranov equation. In contrast to solutions of the classical Grad–Shafranov equation, solutions of the generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad–Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enabling On-Demand Aerospace Component Manufacturing: Topology Optimization of GE Engine Bracket and Fabrication Using Metal FFF

Additive Manufacturing (AM) offers advantages over conventional manufacturing processes, particularly by reducing the number of parts produced through multistage combined technologies, but these often result in low manufacturing yields or require post-processing. AM facilitates the production of complex geometries with fine features, overhangs, and lattice structures. For instance, Laser Powder Bed Fusion (LPBF) technology enables the fabrication of intricate parts that can be easily post-processed by removing residual powder. Laser powder bed AM technologies are widely discussed in the literature due to their design freedom in creating complex geometries, with and without the need for support generation. However, rapid solidification due to a thermal gradient in the build direction, which leads to the formation of columnar grains and warpage, is one of the challenges. To address this challenge, we propose layer-by-layer metal FFF technology, followed by the debinding and sintering process, as an alternative to powder- and laser-based approaches. Furthermore, design for additive manufacturing (DfAM) principles are discussed to minimize the need for support generation and enable easy post-processing, thereby improving surface finish to meet high tolerances in fabricating components for aerospace and healthcare applications.

Singh, Abhishek [University of Michigan, Ann Arbor↗