Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

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 829 records · Page 46

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗

Enabling Efficient Surveillance, Control, and Automation of Geothermal Operations with Advanced Predictive Analytics

Automation and control of geothermal energy production and operations require reliable and efficient predictive tools. While physics-based simulation offers a comprehensive tool for predicting energy production performance in geothermal systems, predicting the behavior of geothermal reservoirs involves complex multi-physics processes with coupling effects, highly uncertain input parameters and subsurface descriptions. Moreover, building, running, and integrating simulation models into standard model calibration and optimization workflows entail significant technical and computational efforts. An emerging alternative to physics-based simulation is data-driven predictive analytics models that have gained popularity in energy industry. In this report, we develop novel predictive models for integration into real-time fault diagnosis and model predictive control algorithms to improve the efficiency of energy production operations in geothermal reservoirs. The report includes two major research Thrust Areas, that is, the surface power plant and the subsurface reservoir.

15 GEOTHERMAL ENERGY↗

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $ \$ $6B annually (at $ \$ $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year), and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE) missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

42 ENGINEERING↗

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $6B annually (at $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year),and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE)missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

99 GENERAL AND MISCELLANEOUS↗

Coupled Experimental–Theoretical Characterization of a Carbon Electrode in Vanadium Redox Flow Batteries using X-ray Absorption Spectroscopy

Vanadium redox flow batteries (VRFBs) have emerged as promising solutions for stationary grid energy storage due to their high efficiency, scalability, safety, near room-temperature operation conditions, and the ability to independently size power and energy capacities. The performance of VRFBs heavily relies on the redox couple reactions of V 2+ /V 3+ and VO 2+ /VO 2 + on carbon electrodes. Therefore, a thorough understanding of the surface functionality of carbon electrodes and their propensity for degradation during electrochemical cycles is crucial for designing VRFBs with extended lifespans. In this study, we present a coupled experimental–theoretical approach based on carbon K edge X-ray absorption spectroscopy (XAS) to characterize carbon electrodes prepared under different conditions and identify relevant functional groups that contribute to unique spectroscopic features. Atomic models were created to represent functional groups, such as hydroxyl, carboxyl, methyl, and aldehyde, bonded to carbon atoms in either sp 2 or sp 3 environments. The interactions between functionalized carbon and various solvated vanadium complexes were modeled using density functional theory. A library of carbon K-edge XAS spectra was generated for distinct carbon atoms in different functional groups, both before and after interacting with solvated vanadium complexes. Here we demonstrate how these simulated spectra can be used to deconvolve ex situ experimental spectra measured from carbon electrodes and to track changes in the electrode composition following immersion in different electrolytes or extended cycling within a functional VRFB. By doing so, we identify the active species present on the carbon electrodes, which play a crucial role in determining their electrochemical performance.

25 ENERGY STORAGE↗

Vegetation heterogeneity reflects soil thermal state and surface soil displacement in a thawing permafrost landscape

Thawing permafrost has the potential to dramatically alter the physical and ecological structure of northern landscapes. Warming of the Arctic and subsequent degradation of permafrost have created a need to assess the stability and movement of soils on hillslopes and the potential impacts on ecosystem structure. In this work, we explore the relationships among vegetation heterogeneity, soil temperature, and soil surface displacements observed from 2019 to 2022 in a watershed in the discontinuous permafrost region on the Seward Peninsula of Alaska. Vegetation heterogeneity was measured as the standard deviation (SD) of the normalized difference vegetation index (NDVI) from 3 m PlanetScope satellite imagery around each soil temperature and active layer thickness observation. Locations of observations were clustered into three soil thermal groups, warm, intermediate, and cold, based on soil temperature and active layer thickness. Average annual horizontal surface displacements were significantly lower for soils within the warm thermal group (median = 0.033 m yr −1 ) compared to soils within the cold thermal group (median = 0.090 m yr −1 ; p < 0.001). Conversely, vegetation heterogeneity was significantly higher in the warm (median = 0.014 SD NDVI; p = 0.002) and intermediate (median = 0.015 SD NDVI; p = 0.002) groups compared with the cold thermal group (median = 0.012 SD NDVI), suggesting a warming-induced shift in vegetation community complexity. Because of the observed associations of ground surface displacement rates and vegetation heterogeneity with soil thermal state, we hypothesize that warming soil conditions induce changes in the rates and patterns of hillslope erosion due to an increase in surface movement as near-surface permafrost thaws, followed by a decrease as the permafrost table deepens and excess ice content diminishes. The transition to warm soils promotes surface ecosystem transformation, shifting the dominant vegetation at the site, given the warming climatic conditions of the region. We integrated our observations of soil temperature, vegetation heterogeneity, and soil surface displacements into a conceptual model that describes the co-evolution of hillslopes and vegetation in warming permafrost environments, which is currently unrepresented in earth system models.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the City-Scale Impacts of Impervious Surfaces on Groundwater Recharge Potential: An Urban Application of WRF–Hydro

Decades of urbanization have created sprawling, complex, and vulnerable cities, half of which are located in water-scarce areas. With the looming effects of climate change, including increasing droughts and water shortages, there is an urgent need to better understand how urbanization impacts the water cycle at city scale. Impervious surfaces disrupt the natural flow of water, affecting groundwater recharge in water-scarce cities, such as Los Angeles, looking to harness local water resources. In the face of growing water demand, informing on opportunities to maximize potential groundwater recharge can help increase cities’ resilience. WRF–Hydro, a physics-based hydrological modeling system, capable of resolving atmospheric, land surface, and hydrological processes at city scale, is adapted to represent urban impervious surfaces. The modified model is used to assess the hydrological implications of historical urbanization. Pre- and post-urban scenarios are used to quantify the impacts of impervious surfaces on the local water budget. Our results show that urbanization in LA has vastly decreased the potential for groundwater recharge, with up to half of the water inflow being redirected from infiltration in highly urbanized watersheds, while doubling surface runoff’s share of the city’s water budget, from ~15% to 30%. This study not only sheds light on the role of imperviousness on groundwater recharge in water-scarce cities, but also offers a robust and transferable tool for the management of urban land and water resources.

54 ENVIRONMENTAL SCIENCES↗

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗

DESI DR2 Galaxy luminosity functions

We present galaxy luminosity functions (LFs) for the Dark Energy Spectroscopic Instrument (DESI) DR2 Bright Galaxy Survey (BGS) in the g, r, z, and $w1$ bands over $0.002\lt z\lt 0.6$. Our analysis uses updated k-corrections and evolutionary corrections, including new polynomial kcorrection fits derived from BGS Year 1 data that supersede earlier GAMA-based prescriptions. Exploiting the statistical power of DESI, we measure LFs to very faint magnitudes, reaching $^{0.1}M_r-5\log h\sim -10$. Independent measurements from the North and South survey regions show excellent agreement around the LF knee, but the very small statistical uncertainties reveal that simple analytic forms fail to capture the full LF shape. The bright end departs from a pure exponential decline, while the faint end exhibits complex, non-powerlaw behaviour, including a pronounced upturn at $^{0.1}M_r-5\log h\gtrsim -15$, which is stronger for red galaxies than for blue. We show that our LFs are largely complete for galaxies with surface brightness $\mu _{50}\lt 25$, and that an apparent steepening fainter than $-13$ is driven primarily by local overdensity and fragmentation of large galaxies. A systematic North–South offset at the brightest magnitudes is traced to red galaxies and may reflect shallower North photometry underestimating extended earlytype profiles, although this remains inconclusive. We therefore also provide LFs based on model Petrosian magnitudes. Redshift splitting reveals small but significant residuals, indicating limitations of a simple global evolutionary model. Using the redshift limits of J. Loveday et al. (2012), we find excellent agreement with GAMA, with substantially reduced statistical errors. These measurements provide a precise reference for studies of environmental and population-dependent LFs and for testing galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗

Air Blast Mesh Sensitivity and Pressure Mapping Study

Nose cone structural and aerodynamic characteristics are essential for intelligent design of aircraft, spacecraft, and ballistic systems. Finite element analysis can be used to help understand the structural integrity and flight characteristics of different nose cones. A mesh sensitivity study was undertaken for a particular nose cone geometry that was used in tests at LANL facilities in order to confirm the integrity of the meshed geometry. A simple cone that best matched closed-form theoretical solutions was modeled, and received good correlation to the theory. Complexity was then added back to the nose cone. Parameters applied to the simple cone were then implemented in the nose cone geometry giving assurance of accuracy after the geometry was changed. Nose cone results averaged 6.3% error for radial displacement when compared with the theoretical. Hoop stress averaged 6.0% error and meridional stress averaged 5.7% error at the finest mesh level. Meshes showed signs of convergence when compared to all three theoretical solutions. Finally, pressure time-history data from LANL computational fluid dynamics simulations was applied to the surface of the final nose cone geometry. The pressure data was interpolated from pressure gauge locations onto nearby meshed elements, which allowed for FEA software to run simulations on the cone with the pressure data as a loading condition. The pressure mapping resulted in the ability to understand the nose cone’s rigid body motion that in turn can inform design of future nose cones.

42 ENGINEERING↗

CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance

Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.

42 ENGINEERING↗

Capillary-Enhanced Two-Phase Micro-Cooler Using Copper-Inverse-Opal Wick with Silicon Microchannel Manifold for High-Heat-Flux Cooling Application

In this work, we demonstrate a two-phase capillary-fed boiling micro-cooler that consists of a ~ 25-..mu..m-thick copper inverse opal (CIO) porous wicking structure for high-heat-flux boiling and a silicon 3D-manifold for distributed liquid delivery and vapor extraction across a 0.5 cm x 0.5 cm heated area. At low inlet water mass flow rates of 1.5 to 1.9 g(min)-1, the micro-cooler displays nearly two-phase boiling with exit vapor quality ~ 1 and a high critical heat flux (CHF) of 253 to 320 W cm-2 with low superheat of ~ 10 degrees C resulting in a thermal resistance of boiling ~ 0.025 cm2 degrees C W-1 or heat transfer coefficient of 0.4 MW m-2 degrees C-1. For higher flow rates of 5, 10, and 15 g(min)-1, the micro-cooler exhibits a hybrid single-phase and two-phase cooling regime where the contribution of the sensible heat (single-phase) cooling is linearly added to that of the two-phase cooling. For the highest flow rate of 15 g(min)-1, the CHF is increased to ~ 500 W cm-2 resulting in an overall thermal resistance of ~ 0.18 cm2 degrees C W-1. However, the two-phase heat transfer effectiveness, which estimates the utilization level of the inlet mass flow rate for two-phase boiling, is reduced to ~ 0.11. To achieve the best cooling system performances, the micro-cooler must operate entirely within the two-phase boiling regime (exit vapor quality or two-phase heat transfer effectiveness ~ 1). Ideally, the "coolant" should be delivered near its saturation temperature (~ 100 degrees C for water), which provides significant advantages for the energy-efficient operation of data centers and power electronics. We present detailed analysis with Infrared and high speed camera images at various inlet flow rates and heat fluxes to understand complex heat transfer in the micro-cooler. Furthermore, a conjugate thermofluidic simulation model, which incorporates the physics of capillary-fed boiling in a porous copper layer, agrees well with the experimental data.

capillary-enhanced boiling↗

Volumetric analysis and mesh generation of real and artificial microstructural geometries

Producing a viable finite element mesh of realistic microstructural structural geometry is a critical step in analyzing the thermo-mechanical behavior of complex multi-material composites. Advancements in imaging technology such as micro computed tomography have allowed modelers to access high resolution mesoscale geometries for direct numerical simulation. However, converting from voxel based 3D images to usable finite element meshes has been challenging. A robust method including algorithms and software scripts for generating finite element meshes from 3D imaged microstructures is presented in this paper. It includes a routine for inserting cohesive elements around material interfaces to enable modeling of interface properties including delamination and damage. The algorithms and procedures presented in this method leverage currently available software packages for processing surface based geometry into volume based meshes. In addition to converting real geometry from physical imaging systems, algorithms for producing numerically generated and statistically equivalent microstructural geometry are also included. These artificial microstructures can be a valuable resource for modelers when physical specimens do not exist or are limited in quantity: Method establishes a workflow from voxel data to viable finite element mesh including interface information, and Includes a method for synthetic geometry generation based on metrics from real microstructures.

36 MATERIALS SCIENCE↗

Dynamics of a Persistent Insulator-to-Metal Transition in Strained Manganite Films

Transition metal oxides possess complex free-energy surfaces with competing degrees of freedom. Photoexcitation allows shaping of such rich energy landscapes. In epitaxially strained La 0.67 Ca 0.33 MnO 3 , optical excitation with a sub-100-fs pulse above 2 mJ/cm 2 leads to a persistent metallic phase below 100 K. Here, using single-shot optical and terahertz spectroscopy, we show that this phase transition is a multistep process. We conclude that the phase transition is driven by partial charge-order melting, followed by growth of the persistent metallic phase on longer timescales. A time-dependent Ginzburg-Landau model can describe the fast dynamics of the reflectivity, followed by longer timescale in-growth of the metallic phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Efficient Hydrogen Delivery for Microbial Electrosynthesis via 3D-Printed Cathodes

The efficient delivery of electrochemically in situ produced H 2 can be a key advantage of microbial electrosynthesis over traditional gas fermentation. However, the technical details of how to supply large amounts of electric current per volume in a biocompatible manner remain unresolved. Here, we explored for the first time the flexibility of complex 3D-printed custom electrodes to fine tune H 2 delivery during microbial electrosynthesis. Using a model system for H 2 -mediated electromethanogenesis comprised of 3D fabricated carbon aerogel cathodes plated with nickel-molybdenum and Methanococcus maripaludis, we showed that novel 3D-printed cathodes facilitated sustained and efficient electromethanogenesis from electricity and CO 2 at an unprecedented volumetric production rate of 2.2 L CH4 /L catholyte /day and at a coulombic efficiency of 99%. Importantly, our experiments revealed that the efficiency of this process strongly depends on the current density. At identical total current supplied, larger surface area cathodes enabled higher methane production and minimized escape of H 2 . Specifically, low current density (<1 mA/cm 2 ) enabled by high surface area cathodes was found to be critical for fast start-up times of the microbial culture, stable steady state performance, and high coulombic efficiencies. Our data demonstrate that 3D-printing of electrodes presents a promising design tool to mitigate effects of bubble formation and local pH gradients within the boundary layer and, thus, resolve key critical limitations for in situ electron delivery in microbial electrosynthesis.

03 NATURAL GAS↗

First lithium experiments in HIDRA and evidence of helium retention during quasi-steady-state stellarator plasma operations

Recent experiments in Hybrid Illinois Device for Research and Applications (HIDRA) have had operational discharges between t discharge = 60 and 1000 s using electron cyclotron resonant heating (ECRH) of the plasma. This means that quasi-steady-state plasma discharges reach conditions to study long-pulse plasma material interactions (PMIs). The newly commissioned HIDRA-Material Analysis Test-stand PMI diagnostic is used to place a drop of lithium onto a heated tungsten surface, transfer the sample in-vacuo and expose it in a helium plasma. Helium is of interest as there is an open question to whether lithium will be able to remove helium ash in real fusion devices. The introduction of the W-Li sample in HIDRA resulted in evaporation of lithium into the helium plasma during a 600 s pulse and caused a reduction of over 90% in neutral pressure during the discharge. It was also observed that the plasma density and temperature increased by over 2.5 times. Using spectroscopy and a helium collisional radiative model, the peak temperature and density of the helium plasma can be monitored during the discharge. During lithium evaporation, as significant lithium ionization occurs, there is a 85% drop in the HIDRA vessel neutral pressure, despite a constant flow rate of He gas. This reduction in neutral pressure is supported by spectroscopy data with corresponding reductions in He I line intensities (587 nm, 667 nm, 706 nm, and 728 nm), as well as those of other impurities. At one point in the discharge a lithium plasma is created, as indicated by an increase in Li + emission and a complete reduction in He + emission, but the electron density jumps from ne = 3 × 10 18 m –3 to over ne = 8 × 10 18 m –3 while the core temperature stays relatively constant between Te = 16 eV and 20 eV. Once lithium has completely evaporated from the sample and the majority of the ionized lithium has diffused from the plasma to the vessel walls, pressure and spectroscopy data paired with He collisional radiative model calculations shows a re-establishment of a helium plasma in a low recycling regime. In this regime, the density drops down to ne = 2 × 10 18 m –3 and the electron temperature increases from T e = 20 eV to over T e = 50 eV indicating an increase in helium heating efficiency. This is also indicated by the He+ emission re-establishing and having a higher intensity. Here, we show the results from the first lithium campaign in HIDRA. In the presence of lithium, and in particular when lithium ions are present, the helium disappears from the plasma via an as of yet unknown complex relationship that needs to be further studied. The most likely explanation is that the lithium ions are distributed around the vessel and able to trap helium to the surface turning HIDRA into a large gettering surface. These results have potential implications on future plasma facing component design using liquid lithium for impurity and recycling control using limiters and divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Methods for Color Center Preserving Hydrogen‐Termination of Diamond

Abstract Chemical functionalization of diamond surfaces by hydrogen is an important method for controlling the charge state of near‐surface fluorescent color centers, an essential process in fabricating devices such as diamond field‐effect transistors and chemical sensors, and a required first step for realizing families of more complex terminations through subsequent chemical processing. In all these cases, termination is typically achieved using hydrogen plasma sources that can etch or damage the diamond, as well as deposited materials or embedded color centers. This work explores alternative methods for lower‐damage hydrogenation of diamond surfaces, specifically the annealing of diamond samples in high‐purity, non‐explosive mixtures of nitrogen and hydrogen gas, and the exposure of samples to microwave hydrogen plasmas in the absence of intentional stage heating. The effectiveness of these methods are characterized by x‐ray photoelectron spectroscopy (XPS), and comparison of the results to density‐functional modelling of the surface hydrogenation energetics implicates surface oxygen ligands as the primary factor limiting the termination quality of annealed samples. Finally, photoluminescence (PL) spectroscopy is used to verify that both the annealing and reduced sample temperature plasma methods are non‐destructive to near‐surface ensembles of nitrogen‐vacancy (NV) centers, in stark contrast to plasma treatments that use heated sample stages.

36 MATERIALS SCIENCE↗