Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “limitations”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Lower and upper threshold limit for artificial neural network based chilled and condenser water temperatures set-point control in a chilled water system

In this study, an ANN (artificial neural network) based real-time optimized control algorithm for a chilled water cooling system was developed and applied in an actual building to analyze its cooling energy saving effects through in-situ application and actual measurements. For this purpose, the cooling tower’s CndWT (condenser water temperature) and the chiller’s ChWT (chilled water temperature) were set as system control variables. To evaluate algorithm performance, the electric consumption and the COP (coefficient of performance) were compared and analyzed when ChWT and CndWTs were controlled conventionally and controlled based on the ANN. During the analysis, unexpected abnormal data was observed due to insufficient training data and limited consideration of OWBT (outdoor air wet-bulb temperature) when determining the CndWT set-point. Therefore, it is necessary to further build training data from a wider range of conditions and to set the lower limit of CndWT set-point to at least +3.6 °C above OWBT when the OWBT is higher than 23 °C, so that further energy savings can be achieved.

42 ENGINEERING↗
DOE OSTI2021

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Performance Analysis and Limiting Parameters of Cross-flow Membrane-based Liquid-desiccant Air Dehumidifiers

We report that to dehumidify a humid air stream, existing air conditioning (AC) systems substantially overcool the outdoor humid air below its dew point, thereby significantly reducing energy efficiency. Directly capturing humidity, membrane-based liquid-desiccant dehumidification systems separate sensible and latent cooling (SSLC) loads and thus offer a promising pathway for a high-performance AC solution. Design of an energy-efficient SSLC-AC system, however, rests largely on detailed understating of the dehumidification process. While some studies have identified the dehumidification process mainly depends on membrane characteristics, other studies have argued that the process is limited by desiccant liquid or alternatively air thermo-hydraulic physics for typical humid climate conditions. The present study examines performance and physics of the membrane-based liquid-desiccant dehumidification process over a wide range of climate conditions through a novel 3D, two-phase, multi-species CFD model. Decoupling the thermodynamic and hydraulic effects, the study reveals that the dehumidification rate is a linear function of the water vapor pressure potential ($J=α ΔP$) summarizing the system's thermodynamic state. The slope of the curve (i.e., α) depends on hydraulic transport characteristics of the membrane pores, air stream, and desiccant solution. More importantly, it was found that the air dehumidification process is mainly limited by the air-side transport physics for thin liquid-desiccant films and commonly used porous superhydrophobic membranes. Additionally, results show that, depending on ambient/desiccant conditions and physical dehumidifier characteristics, energy effectiveness and dehumidification rate vary from 13 to 34% and from 0.13 to 1.4 g m -2 s -1 , respectively. Therefore, the present study allows to efficiently design future SSLC-based AC systems exhibiting high performance energy metrics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Underlying limitations behind impedance rise and capacity fade of single crystalline Ni-rich cathodes synthesized via a molten-salt route

Layered oxide LiNi x Mn y Co z O 2 (NMC) cathodes are often synthesized as polycrystalline secondary particles. Due to intergranular fracture stemming from volume changes of randomly oriented primary particles during charge/discharge, the synthesis of larger single-crystalline cathodes is of high interest. In this work, molten salt assisted growth of micron-sized Ni-rich crystals is achieved with excellent crystallinity, low cation mixing, and negligible impurities. However, electrochemical performance is compromised by high surface reactivity resulting in decomposition of electrolyte and subsequent formation of a thick CEI layer. While intergranular fracture is eliminated, planar gliding and severe intragranular fracture along the (003) plane occurs in the high voltage region within the first few cycles and is associated primarily with H2 to H3 structural transitions. In addition, H2 to H3 transitions are highly irreversible with cyclic voltammograms revealing polarization growth within <5 cycles. Subsequently, the single-crystalline material exhibited markedly reduced available capacity and enhanced capacity fade from sharp impedance growth compared to its polycrystalline counterpart. Here, this work furthers a fundamental understanding into the limitations of single-crystalline Ni-rich cathodes, and the obstacles limiting the advantages offered by the single-crystalline morphology.

25 ENERGY STORAGE↗
DOE OSTI2022

Proper orthogonal decomposition based reduced-order modeling of flux-Limited gray thermal radiation

Here, in this work, a proper orthogonal decomposition (POD) based reduced-order model (ROM) is developed to solve gray, flux-limited thermal radiation diffusion. We focus on the variable opacity radiation penetration benchmark posed by Olson, Auer, and Hall. The T -3 relationship for opacity in conjunction with high-temperature radiation penetrating an initially cold material produces a strong thermal radiation shock. This class of problems is particularly challenging for standard POD-based reduced-order modeling due to the nonlinearities presented by 1) the T 4 source term, and 2) flux-limited diffusion operator. To address these challenges and develop a cost competitive ROM, we employ a “hyper-reduction” technique through discrete empirical interpolation (DEIM) and allow for adaptive reduced-order projections through principal interval decomposition (PID). Performance of the proposed methodology is quantified by comparing the cost savings and accuracy relative to a full-order computation. Reference solutions and snapshot data are obtained through a full-order calculation performed by the University of Chicago maintained astrophysics code, FLASH. For consistency and potential extensibility, the developed ROM is also implemented in FLASH. We find that in the initialization regime, where the thermal radiation wave is initially created by the warming of the material, this class of problems is highly reducible and suitable for POD-based ROMs. However, the strong convective nature of the wave propagation regime is less reducible and more challenging to create an efficient ROM.

42 ENGINEERING↗
DOE OSTI2023

Encoding nonlinear and unsteady aerodynamics of limit cycle oscillations using nonlinear sparse Bayesian learning

This article investigates the applicability of a recently proposed, nonlinear sparse Bayesian learning (NSBL) algorithm to identify and estimate the complex aerodynamics of limit cycle oscillations. NSBL provides a semi-analytical framework for determining the data-optimal sparse model nested within a (potentially) over-parameterized model. This is particularly relevant to nonlinear dynamical systems where modelling approaches involve the use of physics-based and data-driven components. In such cases, the data-driven components, where analytical descriptions of the physical processes are not readily available, are often prone to overfitting, meaning that the empirical aspects of these models will often involve the calibration of an unnecessarily large number of parameters. While an overparameterized model may fit the observed data well, such models may be inadequate for making predictions in regimes that are different from those wherein the data were recorded. In view of this, it is desirable to not only calibrate the model parameters, but also identify the optimal compromise between data fit and model complexity. In this article, we exhibit the optimal model discovery for an aeroelastic system wherein the structural dynamics are well-known and described by a differential equation model, coupled with a semi-empirical aerodynamic model for laminar separation flutter, resulting in low-amplitude limit cycle oscillations (LCO). To illustrate the performance of the algorithm, in this article, we use synthetic data and demonstrate the ability of the algorithm to correctly rediscover the optimal model and model parameters, given a known data-generating model. The synthetic data are generated from a forward simulation of a known differential equation model with parameters selected so as to mimic the dynamics observed in wind-tunnel experiments. Subsequently, we demonstrate the performance of the algorithm for model selection using noisy LCO data from wind tunnel experiments. As there is no ground truth available for the experimental data case, we provide a comparison between NSBL and Bayesian model selection to validate the results, and demonstrate the use of NSBL as an efficient alternative to traditional methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Limited-constraint WAAM fixture for hybrid manufacturing

This paper describes a limited-constraint alternative to the traditional over-constrained build plate clamping in metal additive manufacturing. A fixture design is presented that enables build plate thermal growth within its plane while restricting deformation perpendicular to that plane. A prototype fixture is designed, fabricated, and tested for wire arc additive manufacturing (WAAM). Experimental results are presented for depositing 1.2 mm diameter 5356 aluminum wire on a 5052 aluminum build plate. The build plate deformation is measured for various clamping torques and a 77.3% reduction in distortion is achieved using the limited-constraint approach with identical WAAM parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mg vacancy and impurity-limited MgO single crystal thermal conductivity

Magnesium oxide (MgO) exhibits one of the highest thermal conductivities among oxides and is widely used as a dielectric material and substrate in semiconductor devices, in refractory applications, and as a promising filler in thermal interface materials for electronics. Its high thermal conductivity may be sensitive to impurity and defects, yet this influence is still uncertain. Here, in this study, the impact of the common impurities, i.e., Al, Ca, Ti, V, Fe, Si, B, Nb, Zr, Na, and K, as well as Mg and O vacancies on phonon scattering and thermal conductivity of MgO is studied using a fully first-principles T-matrix framework. It is found that B, Nb, and Zr impurities, along with Mg vacancies, lead to exceptionally strong reductions in thermal conductivity. By contrast, O vacancies and other impurities have modest to minimal impacts. Leveraging the T-matrix results, we reassess the perturbative, mass-only formalism whose use is pervasive in the literature and show that neglecting bond disorder does not necessarily lead to underestimation: for all transition-metal impurities studied, bond perturbations partially cancel mass disorder, causing the traditional perturbative model to overestimate scattering. We propose a simple modified perturbative expression that incorporates both mass and bond disorder and closely reproduces the T-matrix trends. Our predicted low-temperature trends by including phonon-impurity and phonon-boundary scattering match reasonably well with experiments. This work provides an in-depth study of impurity- and vacancy-limited thermal conductivity of MgO and suggests that reported “high-purity” MgO values have likely not yet reached the intrinsic upper limit, which may be substantially higher.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning for reactor power monitoring with limited labeled data

Real-time reactor power monitoring is critical for a variety of nuclear applications, spanning safety, security, operations, and maintenance. While machine learning methods have shown promise in monitoring reactor power levels, there is limited research on their efficacy in label-starved environments. The goal of this work is to assess the feasibility of classifying nuclear reactor power level using multisource data in scenarios with limited labels. Data were collected using low-resolution multisensors at four nuclear reactor facilities: two large research reactors and two TRIGA reactors. Within each pair, one reactor dataset served as the source and the other as the target in a transfer learning paradigm. Twenty-three supervised models were trained on labeled sequences of magnetic field and acceleration data from each of the target sites. Self-learning and transfer learning methods were applied to the top performing models to assess their classification performance with increasing amounts of labeled data. While reactor power level classification was achieved with a Matthews Correlation Coefficient of up to 0.739 ± 0.003 and 0.622 ± 0.009 with only 400 sequences per power state for the large research reactor and TRIGA target sites, respectively, self-learning and transfer learning leveraging source site data did not improve target classification performance. These findings suggest that alternative methods, such as higher sensitivity sensors, digital twins, or the use of physics-informed models, are required to enable high-performance classification in machine learning approaches to reactor monitoring with a dearth of target ground truth.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Switching limits of top-gated carbon nanotube field-effect transistors

The performance limits of carbon nanotube field-effect transistors (CNFETs) based on a recently reported process are studied by computational techniques at temperatures between 300 K and 4 K. Here, the impact of band-to-band tunneling (BTBT) and source-to-drain tunneling (SDT) is examined for devices with varying gate length through the use of simulations based on the non-equilibrium Green’s function (NEGF) formalism, and calibrated to measurements. Additionally, the case of junctionless chemical doping profiles is analyzed in contrast to electrostatic doping recently reported for test structures. The switching limits of CNFETs are further explored for devices based on carbon nanotubes (CNTs) with more favorable electronic structures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Enhancement of disposal efficiency for deep geological repositories based on three design factors - Decay heat optimization, increased thermal limit of the buffer and double-layer concept

This study investigates the enhancement of disposal efficiency for deep geological repositories (DGRs) based on three design factors: decay heat optimization, increased thermal limit of the buffer, and double-layer concept using coupled thermo-hydro-mechanical (THM) numerical simulations. Decay heat optimization is achieved by iteratively emplacing spent nuclear fuels having the maximum and minimum decay heat in a canister. Disposal areas can be reduced by 20 % to 40 % compared to the current reference disposal system in Korea (KRS+) in accordance with the combinations of the three design factors, alleviating challenges in site selection for the DGR. This study additionally identifies an optimal layer spacing of 500 m for the double-layer concept in the viewpoint of the buffer temperature, where thermal interaction between the upper and lower layers nearly disappears. However, determining the ultimate disposal and layer spacing requires engineering judgement, considering not only the thermal performance of the DGR but also various factors such as cost and difficulties of the construction and rock mass stability. DGRs designed with an increased thermal limit of the buffer poses a greater probability of rock mass failure around disposal tunnels and deposition holes due to elevated thermal stresses. Densely arranged heat sources for the DGRs with enhanced disposal efficiency lead to larger temperature increase even at the far-field scale, raising a possibility of thermally driven fracture shear activation with associated hydraulic, mechanical, and seismic changes.

58 GEOSCIENCES↗
DOE OSTI2023

Thermodynamic and Kinetic Barriers Limiting Solid-State Reactions Resolved through In Situ Synchrotron Studies of Lithium Halide Salts

Although halide salts such as LiCl and LiBr are routinely used as a source of Li ions during ion exchange reactions, a detailed understanding of the processes controlling the rates of these reactions is presently lacking. Recently, it has been discovered that the rate limiting barriers for ion exchange are commonly associated with these salts rather than the ceramic target of ion exchange, making it important to quantitatively understand of salt processes. Here it is demonstrated that in situ synchrotron studies of ion exchange reactions can be used to precisely quantify the thermodynamic activation energies associated with these solid state reactions in a manner that can be directly compared with predictions from density functional theory (DFT). While the temperature dependence of the LiCl reaction rate is found to be set by a barrier associated with ion hopping, it was discovered that for LiBr it is also affected by the defect formation energy – an energy found to be substantially lower than predicted by DFT. Furthermore, it is shown that when reaction rates for different relative amounts of reactants are varied, it is possible to identify the rate-limiting reagent and to elucidate an overall scaling relationship that controls the concentration-dependence of the reaction rate. Also, it is demonstrated that global fits across doped and undoped salts can be used to probe both intrinsic and extrinsic vacancy concentrations. This improved understanding of ion exchange mechanisms can be used to accelerate ion exchange reaction rates by orders of magnitude. In conclusion, the techniques demonstrated here can be broadly applied to probe the kinetics and thermodynamics of solid state reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solid-State, Single-Ion Conducting, Polymer Blend Electrolytes with Enhanced Li + Conduction, Electrochemical Stability, and Limiting Current Density

The development of solid-state polymer electrolytes with high lithium conductivity is crucial to improve lithium-ion battery performance and ameliorate the safety challenges associated with current solvent-based electrolytes. Unfortunately, sluggish polymer segmental dynamics are known to constrain conductivity enhancements in solid-state polymer electrolyte systems, limiting overall performance. In this work, a glassy single-ion-conducting polymer, poly[lithium sulfonyl(trifluoromethane sulfonyl)imide methacrylate] (PLiMTFSI), was blended with a flexible polymer, poly(oligo-oxyethylene methyl ether methacrylate) (POEM), and the impact of PLiMTFSI molecular weight and ion concentration on the thermal and ion-conducting behavior of blend electrolytes was investigated. High ionic conductivities approaching 1 × 10 -2 S/cm at 150 °C were realized in this polymer blend electrolyte system as a result of decoupling Li + transport from polymer segmental dynamics. The decoupled ion transport was attributed to the packing frustration of the glassy PLiMTFSI – sufficient percolating free volume was generated to produce effective ion diffusion pathways. This decoupling was tunable as the ion transport could be altered from being closely coupled to the polymer segmental dynamics (Vogel–Tammann–Fulcher-like) to hopping (Arrhenius-like) by increasing the PLiMTFSI molecular weight and ion concentration. Moreover, the immobilized TFSI anion resulted in high Li + selectivity (Li + transference number = 0.9), high electrochemical stability (up to 4.7 V against Li + / Li), and limiting current density of 1.8 mA/cm 2 (electrolyte thickness = 0.05 cm). These features suggest that this single-ion-conducting, polymer blend electrolyte might be a promising alternative to a benchmark system – salt-doped polyethylene oxide. Moreover, the above characteristics can support the battery operation at higher voltages using energy-dense Li metal anodes, with faster charging rates and enhanced energy/power densities. Altogether, the results suggest that polymer chain packing frustration can be exploited to overcome the constraints of slow polymer segmental relaxations to achieve rapid and highly selective ion transport and enhanced performance in solid-state polymer electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Designing Particle Morphologies for Materials with Solid Transport Limitations: A Case Study of Lithium and Manganese Rich Cathode Oxides

A lithium and manganese rich nickel-manganese-cobalt oxide (LMR-NMC) cathode is a promising candidate for next-generation batteries due to its high specific capacity, low cost, and low cobalt content. However, the material suffers from poor rate capability due to the diffusion limitations of lithium in the cathode particles. Understanding the material performance requires careful control of the morphology of the cathode particles, taking into account the primary and agglomerated diffusion pathways and the presence of pores, some of which could be closed from electrolyte infiltration. Here, in this study, we use a microstructure-based mathematical model combined with experimental data to understand the role of the complex cathode particle morphology in the rate performance of the material. Scanning electron microscopy images of cathodes made under different synthesis conditions, which results in different agglomerate morphologies, serve as the input into the mathematical model. The model is then compared to rate data to understand the controlling parameters. The presence of intra-agglomerate closed pores results in a large agglomerate diffusion length in comparison to the ideal condition, where the primary particles are agglomerated in an open and dispersed manner such that the entire interfacial area is available for electrochemical reaction. Smaller primary and agglomerate diffusion lengths result in better electrochemical performance. This points us toward designing the morphology of the cathode particles to compensate for the diffusion limitation of LMR-NMC while maximizing the density.

Tewari, Deepti↗
DOE OSTI2026

Interference-Limited Absorption in Dense Molecular Nanolayers Near Reflecting Surfaces

We investigate linear resonant absorption by a dense ensemble of molecules confined to a sub-wavelength layer in two geometries: (i) a free-standing film in a homogeneous space and (ii) the same film placed at a controlled distance from a reflecting surface. In both cases, increasing the effective light–matter coupling (via molecular density/oscillator strength) produces a nonmonotonic response: absorption rises to an optimum and then decreases as the film becomes increasingly radiatively bright and reflective. Finite-difference time-domain simulations and analytical transfer-matrix calculations agree quantitatively and yield compact ridge conditions for the optimum. We interpret the trends using a scattering/port picture: the isolated film is a symmetric two-port system (reflection and transmission), which bounds single-sided resonant absorption to ≤50% in the ultrathin limit (reflecting transition saturation), whereas adding a mirror suppresses transmission and converts the structure into an effectively one-port absorber. In the mirror-backed geometry, interference can cancel reflection, and unity absorption is obtained at critical coupling, when radiative leakage is balanced by intrinsic molecular loss. These results clarify fundamental limits and design rules for collective absorption in dense molecular layers near dielectric or metallic boundaries.

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

Ammonolysis Under NH 3 –Limiting Conditions as a Pathway to Improved LaTiO 2 N Water Splitting Photoanodes

LaTiO 2 N is a promising intermediate band gap semiconductor for the water splitting reaction, a pathway to hydrogen fuel from solar energy. However, the photoelectrochemical (PEC) activity of the material is hindered by defects, particularly Ti(III) species, which promote photocarrier recombination. These defects are formed during the high-temperature ammonolysis reaction. Here we show that improved LaTiO 2 N materials can be synthesized under NH 3 -limiting conditions by introducing N 2 to lower the NH 3 partial pressure to0.13atm.This reduces the Ti(III) defect density in the material from 6.06 × 10 16 to ∼4.61 × 10 15 cm −3 , by a factor of 13, based on electron paramagnetic resonance (EPR) spectroscopy. Any remaining Ti(III) defects are localized at the LaTiO 2 N surface, according to X-ray photoelectron spectroscopy (XPS), due to the formation of a depletion layer in the semico. Optical absorption spectra of the improved LaTiO 2 N reveal a blue-shifted band gap absorption edge and a suppressed sub-band gap absorption. Defect removal also reduces a sub-band gap surface photovoltage feature visible in the 1.0 atm reference material. The improved LaTiO 2 N supports a 1.57 mA cm −2 water oxidation photocurrent at 1.23 V RHE under simulated sunlight conditions, and an enhanced quantum efficiency of 4.5% (400 nm) for photocatalytic oxygen evolution from aqueous silver nitrate solution. Stable PEC operation is observed for over 55 min. This confirms that ammonolysis under NH 3 -limiting conditions improves the solar energy conversion properties of LaTiO 2 N. The ability to control metal ion defects in oxynitrides by varying the ammonia partial pressure during ammonolysis might be generally useful for the preparation of metal nitrides and oxynitrides.

defects↗
DOE OSTI2023

Limits of Detection for EXAFS Characterization of Heterogeneous Single-Atom Catalysts

Single-atom catalysts (SACs), consisting of individual metal atoms dispersed on a support, attract attention due to their unique reactivity, efficient use of precious metals, and precise chemical tunability. Characterization of the metal species is crucial to substantiate structure–function relationships. Authors often use—and referees often require—X-ray absorption spectroscopy (XAS) data to prove the absence of clustered metal (or metal oxide) structures after pre-treatment and under in situ or operando conditions. However, there has been no critical assessment of the limitations of XAS in substantiating such conclusive statements, which is particularly important given the potential outsized influence of minority catalyst structures in dictating catalytic activity. Here, in this article, we quantitatively assess the detection limits of XAS to identify metal (or metal oxide) clusters in samples containing predominantly single atoms by modeling the extended X-ray absorption fine structure (EXAFS) of mixtures of structures. We identified that a significant fraction of clusters can coexist with SAC active sites (e.g., ~10% metallic Pt or ~40% oxidized Pt clusters in Pt/CeO 2 SACs), while eluding detection via EXAFS with any statistical significance. To generalize these conclusions, a descriptor-based screening of bulk metal oxides using a continuous Cauchy wavelet transform was proposed that suggests certain materials for which differentiating atomically dispersed metal species and metal oxide clusters would be infeasible by EXAFS (e.g., ReO x ). Based on this analysis, we suggest best practices for the study of SACs using EXAFS and provide recommendations to ensure that conclusions do not outpace the evidence used to support them. In this rapidly expanding research area, rigorous characterization will lead to greater understanding of the behavior of SACs and ultimately improved catalytic materials.

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

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗