Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energy Transition Initiative”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Snow Albedo Evolution (SALVO) Campaign at the North Slope of Alaska Field Campaign Report

The springtime surface-albedo transition in the Alaskan Arctic, along with the forces that determine its duration and nature, was the focus of the Snow ALbedo eVOlution (SALVO) campaign. The SALVO team investigated the reasons and durations of the stages of spring melt, during which albedo values decrease from 0.8 to 0.1, signifying the year’s largest and most significant radiative energy change. SALVO II (2022, 2024) built on findings from the successful SALVO I (2019) melt season campaign and previous research on the melt transition period conducted by Grenfell and Perovich (2004). SALVO used the ARM NSA observatory, thus aligning with a core principle of the ARM Decadal Vision to “provide comprehensive and impactful field measurements to support scientific advancement of atmospheric process understanding.” SALVO fieldwork was conducted in the spring near the ARM NSA Central Facility in Utqiagvik, Alaska (Figure 1). Comprehensive data sets were collected in 2019, 2022, and 2024, which included spectral and broadband albedos, snow and meltwater depths, snow stratigraphy, snow grain size data, and aerial imagery. In April of each project year, before the snow began to melt, we established survey lines at three or four locations: inland tundra at NSA E12, coastal tundra at the NSA C1 site, Elson Lagoon, and offshore on the Chukchi Sea sea ice (2022 only). We set up a 200-meter-long line, marked every five meters, to enable repeated measurements of surface albedo and snow depth at each location. During most site visits, the SALVO II team dug at least one snow pit to assess the characteristics of the vertical snow layers, including the types and sizes of snow grains, snow density, and the distribution of liquid water. We evaluated the characteristics of flowing or pooled water at the bottom of the snowpack. Watching snow melt in the Arctic has been compared by some to being as exciting as watching paint dry, but it is far more exciting and dynamic than that. Initially, changes happen slowly while the snow cover remains above 80%, but then reductions in the albedo of the landscape, from about 0.8 to 0.4, combined with rising spring temperatures, begin to accelerate the melt. It happens so quickly that no matter how hard the field team works, they cannot keep up with documenting the changes. At first, liquid water is scarce—only a few wet layers of snow in the snowpack. But then, water becomes ubiquitous, transporting melt energy and creating an extremely heterogeneous albedo and melt landscape. At the end of each field season, the SALVO team is exhausted, relieved to see the snow gone, and joyful to play in the melt ponds on the sea ice under the midnight sun. Overall, the 2024 SALVO field campaign was a highlight, reflecting numerous lessons learned from 2019 and 2022. Daily measurements were taken more than 10 times at each site, resulting in a comprehensive time series of snow conditions and albedo. The team deployed using snowmachines and then on foot when snowmachine travel was no longer permitted on the melting tundra. Instrument mounts and packaging were optimized for quick deployment, regardless of the mode of transportation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modulating Surface Anionic Redox Chemistry toward Highly Stable Li-Rich Cathodes with Negligible Oxygen Loss

Low initial Coulombic efficiency and severe capacity/voltage fading during cycling caused by serious irreversible oxygen release, especially in the initial cycle, and resultantly induced unstable electrode/electrolyte interfacial chemistry, largely prohibit the commercial application of high-capacity Li-rich layered oxide cathodes (LLOs). Here, in this work, a dual reductive gas interface cotreatment strategy is applied to regulate the lattice oxygen redox activity and reversibility with a multiple defective structure design including Li/O/TM (TM = transition metal) vacancies and the intrinsic TM doping as well as a full-surface protective layer, which can suppress the irreversible TM migration and then undesirable phase transformation, resisting the corrosion of electrolyte during cycling effectively. Importantly, the introduced reversible SO 3 2- /SO 4 2- redox couple that provides extra capacity compensation could alleviate the distortion of oxygen-central octahedral structure and structural collapse caused by immoderate oxygen oxidation. Thus, the lattice oxygen redox chemistry is optimized, with negligible oxygen loss during the initial cycle. And the designed AS-LLO cathode with greatly enhanced structure stability shows high-capacity retentions of 99.2% at 0.3C after 100 cycles and 82.4% even after 1000 cycles at 5C. This work provides a guideline for manipulating the oxygen redox chemistry to achieve long-lifespan Li-rich layered oxide cathodes for high-energy-density lithium batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Time-Temperature-Transformation (TTT) Diagram for a Sludge Batch 9 Glass Composition Based on Coupled-Operation with the Salt Waste Processing Facility

The amorphous structure of a glass waste form has the potential to rearrange into crystalline phases at temperatures between the liquidus temperature and the glass transition temperature (Tg). Certain phases that can form will be detrimental to the durability of the glass and it is important to know the conditions that promote devitrification. The canister-centerline-cooling (CCC) profile is used to replicate the area within the center of the Defense Waste Processing Facility (DWPF) canister with the slowest cooling during the initial cool down after pouring, which has the greatest potential for crystallization. Other time-temperature conditions that cause significant changes in either phase structure or phase composition are identified by a time-temperature-transformation (TTT) study. The phase stability of a waste form must be determined as a part of the Waste Acceptance Product Specifications (WAPS) if it is to eventually be stored in a geologic repository. This requires the creation of a TTT diagram and analysis of the Tg, as defined by the Department of Energy (DOE). The previous TTT study for a DWPF glass waste form was completed in 2010 prior to coupled operation with the Salt Waste Processing Facility (SWPF). SWPF transfers two high activity waste streams to DWPF for vitrification: a cesium-containing strip effluent and a stream containing monosodium titanate/sludge solids. These SWPF streams were first transferred to DWPF for vitrification during Sludge Batch 9 (SB9) in 2021. The impact of these SWPF streams on crystallization behavior was not determined in previous studies and the need for data to satisfy WAPS Specification 1.4 for SB9 was identified.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗
DOE OSTI2024

The path to 2060: Saudi Arabia's long-term pathway for GHG emission reduction

Saudi Arabia, as part of its Saudi Green Initiative, has announced its goal to achieve net zero green-house gas emissions by 2060. This ambitious target underscores the nation's dedication to address-ing climate change. However, there is a significant gap in comprehensive analysis regarding the long-term effects of Saudi Arabia's climate policies and their collective contribution towards the net-zero objective. This study endeavors to bridge this gap through a detailed examination using the GCAM-KSA, a specialized version of the Global Change Analysis Model tailored for Saudi Arabia, employing a multi-sectoral methodology that integrates economic, energy, and land use systems within a coherent framework to assess the impact of climate policies on GHG emissions. Our anal-ysis reveals that reaching net-zero GHG emissions by 2060 is a complex challenge requiring con-certed efforts across all sectors of the economy. While transitioning to low-carbon electricity and improving energy efficiency offer considerable emission reductions, fully decarbonizing the indus-trial and transportation sectors poses a significant hurdle. Our findings suggest that Saudi Arabia must triple its emission reduction commitments in its next Nationally Determined Contributions (NDCs) update to align with its 2060 net-zero goal. Early action and increased ambition could avoid the chances of getting locked into the high emission assets and give enough time to transform the energy system. Furthermore, the adoption and integration of Carbon Dioxide Removal (CDR) tech-nologies are identified as crucial for offsetting residual emissions, especially in sectors that might continue to rely on fossil fuels.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS2011

Insect Residue Contamination on Wing Leading Edge Surfaces: A Materials Investigation for Mitigation

Flight tests have shown that residue from insect strikes on aircraft wing leading edge surfaces may induce localized transition of laminar to turbulent flow. The highest density of insect populations have been observed between ground level and 153 m during light winds (2.6 -- 5.1 m/s), high humidity, and temperatures from 21 -- 29 C. At a critical residue height, dependent on the airfoil and Reynolds number, boundary layer transition from laminar to turbulent results in increased drag and fuel consumption. Although this represents a minimal increase in fuel burn for conventional transport aircraft, future aircraft designs will rely on maintaining laminar flow across a larger portion of wing surfaces to reduce fuel burn during cruise. Thus, insect residue adhesion mitigation is most critical during takeoff and initial climb to maintain laminar flow in fuel-efficient aircraft configurations. Several exterior treatments investigated to mitigate insect residue buildup (e.g., paper, scrapers, surfactants, flexible surfaces) have shown potential; however, implementation has proven to be impractical. Current research is focused on evaluation of wing leading edge surface coatings that may reduce insect residue adhesion. Initial work under NASA's Environmentally Responsible Aviation Program focused on evaluation of several commercially available products (commercial off-the-shelf, COTS), polymers, and substituted alkoxy silanes that were applied to aluminum (Al) substrates. Surface energies of these coatings were determined from contact angle data and were correlated to residual insect excrescence on coated aluminum substrates using a custom-built "bug gun." Quantification of insect excrescence surface coverage was evaluated by a series of digital photographic image processing techniques.

Lorenzi, Tyler M.↗
DOE OSTI2020

Initial Implementation Of Machine Learning System for SRF Cavity Fault Classification at CEBAF

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Laboratory is high power CW electron accelerator. It uses a mixture of two types of SRF modules, older, lower energy C20/C50 modules and newer higher energy C100 modules, arrayed in two anti-parallel linear accelerators. Accurately classifying the type of cavity faults is essential to maintaining and improving accelerator performance. Each C100 cryomodule contains eight 7-cell cavities. When a cavity fault occurs within a cryomodule, all eight cavities generate 17 waveforms each containing 8192 points. This data is exported from EPICS and saved for review. Analysis of these waveforms is time intensive and requires subject matter expertise (SME). SMEs examine the data from each event and label it according to one of several known cavity fault types. Multiple machine learning models have been developed on this labeled dataset with sufficient performance to warrant the creation of a limited machine learning software system for use by accelerator operations staff. This paper discusses the transition from model development to implementation of a prototype system.

Carpenter, A.↗
NASA NTRS2004

Advanced Life Support Project: Crop Experiments at Kennedy Space Center

Crop production systems provide bioregenerative technologies to complement human crew life support requirements on long duration space missions. Kennedy Space Center has lead NASA's research on crop production systems that produce high value fresh foods, provide atmospheric regeneration, and perform water processing. As the emphasis on early missions to Mars has developed, our research focused on modular, scalable systems for transit missions, which can be developed into larger autonomous, bioregenerative systems for subsequent surface missions. Components of these scalable systems will include development of efficient light generating or collecting technologies, low mass plant growth chambers, and capability to operate in the high energy background radiation and reduced atmospheric pressures of space. These systems will be integrated with air, water, and thermal subsystems in an operational system. Extensive crop testing has been done for both staple and salad crops, but limited data is available on specific cultivar selection and breadboard testing to meet nominal Mars mission profiles of a 500-600 day surface mission. The recent research emphasis at Kennedy Space Center has shifted from staple crops, such as wheat, soybean and rice, toward short cycle salad crops such as lettuce, onion, radish, tomato, pepper, and strawberry. This paper will review the results of crop experiments to support the Exploration Initiative and the ongoing development of supporting technologies, and give an overview of capabilities of the newly opened Space Life Science (SLS) Lab at Kennedy Space Center. The 9662 square m (104,000 square ft) SLS Lab was built by the State of Florida and supports all NASA research that had been performed in Hanger-L. In addition to NASA research, the SLS Lab houses the Florida Space Research Institute (FSRI), responsible for co-managing the facility, and the University of Florida (UF) has established the Space Agriculture and Biotechnology Research and Education (SABRE) Center with several faculty.

Sager, John C.↗
DOE OSTI2020

Closing the Loops on Solar Photovoltaics Modules: An Agent-Based Modeling Approach for the Study of Circular Economy Strategies

Solar photovoltaics (PV) installed capacity have grown exponentially since the early 2000s (average annual growth rate of 50%). With 4,700 GW projected installed capacity by 2050, the volume of PV panels waste is also expected to become substantial. Though renewables are a sine qua non to the establishment of a truly circular economy (CE), the issue arising from their end-of-life management needs to be resolved. The management of PV end-of-life also represent a singular opportunity to create value, from recovered valuable, rare or critical materials (e.g. silver, tellurium, indium). Moreover, the PV case illustrates some of the current barriers to CE; for instance, the need for a common definition of waste (PV are defined as e-waste in the European Union but as general waste in the United States) and potential loss of innovations' advantages (PV average efficiency has continuously grown the past 10 years). In this study, an agent-based modeling (ABM) approach is proposed to simulate PV end-of-life management in the United States. The model explores how the decisions of the various actors involved in handling PV waste affect the quantities of PV that are reused, recycled, or land-filled. In the model, manufacturers, residential, and nonresidential PV owners, installers, and recyclers are represented by agents while governmental policies and regulations are treated as exogenous variables. PV-market data are used to define agents' characteristics (e.g., installed PV capacities or producing costs). Agents' decision rules draw on the literature related to industrial symbiosis (IS) and peoples' waste behaviors as they appear as the most widely used model to implement CE principles at a meso level. Specifically, the concepts of mutual trust between IS actors and knowledge about the IS philosophy are yielded to define the agents' decision process. The primary outputs of the ABM are the costs and volume of waste associated with each end-of-life pathway. Preliminaries results indicate recycling of residential PV is affected by the cost, perceived difficulty of recycling behaviors, and social norms. Moreover, the type of network connecting agents and the initial recycling rate strongly influence results. The model also highlights the crucial role of recycling behavior adoption: in case of early failure of PV, recycled volumes increase by about 40%. This represents an additional 5.7 billion USD from potential recovered silver. Finally, although results for the PV case are presented here, the developed ABM aims at being a general tool for the study of CE strategies and the CE transition. Furthermore, steps of this research include extending the model to encompass circular economy strategies (e.g., design for recycling or lifetime extension) and validating its general architecture from various case studies.

14 SOLAR ENERGY↗
DOE OSTI2023

Redox Mechanisms and Migration Tendencies in Earth-Abundant 0.7Li 2 MnO 3 ·0.3LiFeO 2 Cathodes: Coupling Spin-Resolved X-ray Absorption Near Edge and X-ray Absorption Fine Structure Spectroscopies

We report the use of iron 1s3p resonant X-ray emission processes to conduct spin-selective, high-energy resolution fluorescence detected X-ray absorption near-edge spectroscopy (HERFD-XANES) on an iron-containing, lithium- and manganese-rich, fully earth-abundant cathode material, Li 1.3 Mn 0.5 Fe 0.2 O 2 (0.7Li 2 MnO 3 ·0.3LiFeO 2 ). Coupling this technique with conventional Mn K-edge XANES and detailed extended X-ray absorption fine structure (EXAFS) analysis from both the Mn and Fe vantage points, we gain fundamental insights into the redox processes and migration tendencies of transition metals in this cathode material at the bulk level. We show that during the first charge, Fe 3+ undergoes oxidation to form Fe 4+ prior to the activation plateau. Toward the end of activation, a significant fraction of the iron is present as tetrahedral Fe 3+ . This observation reveals that iron migration from octahedral to tetrahedral sites and iron reduction are initiated during activation. Upon first discharge from the activated state, a continuous and overlapping reduction of both Fe and Mn is observed, with Fe largely restored back as an octahedrally coordinated Fe 3+ . The manganese local environment gradually changes to a distorted cooperative Jahn–Teller Mn 3+ structure during discharge, with the clear presence of two Mn–O as well as two Mn–Mn correlation distances at 2.0 V. The significant reduction of manganese in the very first discharge is distinctly different from that seen in typical nickel-based lithium-manganese-rich materials but is similar to that observed for pure Li 2 MnO 3 . In conclusion, these findings shed light on key structure–property correlations in the cathode material and point to a causative relationship between the redox mechanisms as well as structural changes endured by the material and relatively poor performance during extended electrochemical cycling.

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

EMPOWERED Distributed Energy Resources Permit Accelerator Pilot

This report summarizes the EMPOWERED Distributed Energy Resources (DER) Permit Accelerator Pilot project. Its primary objective was to streamline the design, installation, permitting, and inspections of targeted Distributed Energy Resource (DER) solutions. Streamlining these processes is critical to strengthening and accelerating the adoption of DER solutions nationwide. Moreover, streamlined processes help expedite the clean energy transition and expand energy resilience by improving understanding and capability among key stakeholders like property owners, code officials, and the general workforce. The project team accomplished this by developing and deploying a set of design and permitting guides for simple DER solutions, which they rolled out via a permitting pilot program. The project team developed nine permitting and inspection guides for simple DER solutions, as follows: 1. Single-Family and Duplex: Electric Vehicle Service Equipment (EVSE), and Storage, Solar + Storage (each in two code cycles) 2. Multi-family and Office: EVSE, Storage, Solar + Storage (each in two code cycles) 3. A permitting process guide to support the implementation of the permitting guides. The guides cover key code requirements via plan review and field inspection checklists for code officials and building owners to follow to streamline the permitting and inspections process for DER solutions. Following the development of the initial set of guides, the project team formed two cohorts, each made up of four jurisdictions, to participate in the Permit Accelerator Pilot. Each cohort received training and technical assistance to support them in incorporating the guides into use. The project successfully engaged jurisdictions in two disparate US regions: Chelsea, Somerville, Natick, and Norwood in metro Boston; and Pima County, Town of Gilbert, Flagstaff, and Sedona in Arizona. Feedback gathered from the project was valuable in refining the guides to better address jurisdictional needs. After the pilot period, the project team developed final versions of the guides based on extensive feedback from the cohorts and from external technical peer review. Feedback included needs and barriers related to implementation, as well as detailed technical improvements from external peer review. The final guides were presented in a series of webinars with a national reach, with 192 attendees from 35 states in the live session of the final webinar, and the final guides are posted online and available for use nationwide. Key findings: • The most receptive jurisdictions were those adjacent to other jurisdictions that had already adopted similar initiatives. The least receptive jurisdictions were those in the midst of adopting new code cycles. • The most common barrier to jurisdictional adoption of new guidelines were internal administrative and process delays, more so than technical barriers. • The main technical barrier to adoption proved to be the concern of fire hazard risks posed by energy storage systems. There is some apprehension in the code enforcement community about the safety of the batteries. Education coordinated with fire and safety services will help alleviate concerns and increase confidence in the acceptance of innovative technologies. • Opportunities to further leverage and advance the guides and related resources include outreach, education, and jurisdictional support; updates to the resources to align with code cycle changes and local requirements; and advancements supporting emerging technologies.

14 SOLAR ENERGY↗
DOE OSTI2022

Clean energy technology pathways from research to commercialization: Policy and practice case studies

Clean energy research and development (R and D) leading to commercial technologies is vital to economic development, technology competitiveness, and reduced environmental impact. Over the past 30 years, such efforts have advanced technology performance and reduced cost by leveraging network effects and economies of scale. After demonstrating promise in applied R and D, successful clean energy and energy efficiency technologies are incorporated into an initial product sold by the private sector. Despite its importance, processes by which first commercialization occurs are difficult to generalize while capturing specific insights from practitioners in markets and technologies. This paper presents a policy-focused qualitative assessment of the first commercialization of four diverse energy technologies: thin film photovoltaics, wind turbine blades, dual-stage refrigeration evaporators, and fuel cells for material handling equipment. Each technology presents distinct value propositions, markets, and regulatory drivers. The case studies indicate three common characteristics of successful first commercialization for new energy technologies: 1) good fit between the technology, R&D infrastructure, and public-private partnership models; 2) high degree of alignment of government regulations and R&D priorities with market forces; and 3) compatibility between time scales required for R&D, product development, and opportunities. These findings may inform energy investment decision-making, maximize benefits from R&D, and advance the transition to a low-emission future.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Modeling the Effect of Surface Platinum–Tin Alloys on Propane Dehydrogenation on Platinum–Tin Catalysts

Uncertainty analysis, reported experimental literature data, and density functional theory were synthesized to model the effect of surface tin coverage on platinum-based catalysts for nonoxidative propane dehydrogenation to propylene. Here, this study tests four different platinum–tin skin surface models as potential catalytic sites, Pt 3 Sn/Pt(100), PtSn/Pt(100), Pt 3 Sn/Pt(111), and Pt 2 Sn/Pt(211), and compares them to the corresponding pure Pt surface sites using an uncertainty analysis methodology that uses BEEF-vdW with its ensembles (BMwE) to generate the uncertainty for the energies of the intermediates and transition states. One experimental data set with two experimental observations, selectivity to propylene and turnover frequency of propylene, was used as a calibration data set to evaluate the impact of the experimental data on informing the models. This study finds that the prior model for Pt 3 Sn/Pt(100) is the most active and Pt 2 Sn/Pt(211) is the most selective toward propylene. Active sites on the (100) facet have the highest probability of being responsible for C 1 and C 2 product formations (C–C bond cleavage). Increasing the Sn coverage on the (100) surface facet to a PtSn/Pt(100) active site leads to a significantly reduced rate and might explain the experimentally observed higher selectivity of Sn-doped catalysts relative to pure Pt catalysts. Next, this study finds that for all surfaces, except PtSn/Pt(100), the rate-controlling steps are the initial dehydrogenation steps alongside some partially rate-controlling second dehydrogenation steps. For PtSn/Pt(100), only the initial terminal dehydrogenation step to CH 3 CH 2 CH 2 * and second dehydrogenation steps are rate-controlling. Next, the calibrated models for all surfaces were found to be selective toward propylene production and model the reported turnover frequency successfully. Nevertheless, Pt 2 Sn/Pt(211) emerges as the active site with some (minor) evidence as the main active site based on Jeffreys’ scale interpretation of Bayes factors. This observation agrees with prior studies that also found step sites to be most likely the most relevant active sites for pure Pt catalysts. Overall, the results indicate that tin, in addition to affecting the binding strength of the adsorbed species, prevents deeper dehydrogenation (reducing coking) and cracking reactions through increasing activation barriers for unwanted side reactions.

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

Hysteresis between gas breakdown and plasma discharge

In direct-current (DC) discharge, it is well known that hysteresis is observed between the Townsend (gas breakdown) and glow regimes. Forward and backward voltage sweep is performed using a one-dimensional particle-in-cell Monte Carlo collision (PIC-MCC) model considering a ballast resistor. When increasing the applied voltage after reaching the breakdown voltage (V b ), transition from Townsend to glow discharges is observed. When decreasing the applied voltage from the glow regime, the discharge voltage (V d ) between the anode–cathode gap can be smaller than the breakdown voltage, resulting in a hysteresis, which is consistent with experimental observations. Next, the PIC-MCC model is used to investigate the self-sustaining voltage (V s ) in the presence of finite initial plasma densities between the anode and cathode gap. It is observed that the self-sustaining voltage coincides with the discharge voltage obtained from the backward voltage sweep. In addition, the self-sustaining voltage decreases with increased initial plasma density and saturates above a certain initial plasma density, which indicates a change in plasma resistivity. The decrease in self-sustaining voltage is associated with the electron heat loss at the anode for the low pd (rarefied) regime. In the high pd (collisional) regime, the ion energy loss toward the cathode due to the cathode fall and the inelastic collision loss of electrons in the bulk discharge balance out. Finally, it is demonstrated that the self-sustaining voltage collapses to a singular value, despite the presence of a initial plasma, for microgaps when field emission is dominant, which is also consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Development of a Numerical Model of Hypervelocity Impact into a Pressurized Composite Overwrapped Pressure Vessel

As the outlook for space exploration becomes more ambitious and spacecraft travel deeper into space than ever before, it is increasingly important that propulsion systems perform reliably within the space environment. The increased reliability compels designers to increase design margin at the expense of system mass, which contrasts with the need to limit vehicle mass to maximize payload. Such are the factors that motivate the integration of high specific strength composite materials in the construction of pressure vessels commonly referred to as composite overwrapped pressure vessels (COPV). The COPV consists of a metallic liner for the inner shell of the COPV that is stiff, negates fluid permeation and serves as the anchor for composite laminates or filaments, but the liner itself cannot contain the stresses from the pressurant it contains. The compo-site-fiber reinforced polymer (CFRP) is wound around the liner using a combination of hoop (circumferential) and helical orientations. Careful consideration of wrap orientation allows the composite to evenly bear structural loading and creates the COPV's characteristic high strength to weight ratio. As the CFRP overwrap carries most of the stresses induced by pressurization, damage to the overwrap can affect mission duration, mission success and potentially cause loss-of-vehicle/loss-of-crew. For this reason, it is critical to establish a fundamental understanding of the mechanisms involved in the failure of a stressed composite such as that of the COPV. One of the greatest external threats to the integrity of a spacecraft's COPV is an impact from the meteoroid and orbital debris environments (MMOD). These impacts, even from submillimeter particles, generate extremely high stress states in the CFRP that can damage numerous fibers. As a result of this possibility, initial assumptions in survivability analysis for some human-rated NASA space-craft have assumed that any alteration of the vessel due to impact is considered a catastrophic failure. This assumption is conservative and made due to lack of knowledge on the level of allow-able damage to the composite overwrap that can be sustained and still allow successful completion of the mission. To quantify the allowable damage level to the composite overwrap involves assessing stress redistribution following damage as well as evaluating possible time-dependent mechanisms involved in the COPV response to an impact event. Limited published work in this subject has shown that COPV can withstand at least some level of damage due to high energy impacts. These observations have been confirmed and expanded upon in recent experimental research performed by NASA. This research has demonstrated that there is not only robustness in a COPV to compensate for CFRP damage, but has also identified two significant failure modes for pressurized COPV. The lowest threshold failure mode involves the perforation of the vessel, and the highest threshold failure mode is the catastrophic rupture. While both of these failure modes mean a loss of the COPV, system robustness affords some tolerance to the venting as opposed to the more catastrophic rupture. As a consequence, it is necessary to understand the conditions that result in the transition between these failure modes. The aforementioned experimental research has been performed in both the unpressurized and pressurized condition to identify the damage level that triggered the failure thresh-old. This COPV test program was sponsored by the NASA Engineering and Safety Center (NESC), and tests were performed at NASA White Sands Test Facility (WSTF). Planning and coordination were provided by NASA JSC Hypervelocity Impact Technology (HVIT) group, and the COPVs were provided by the ISS Program. Unpressurized testing has been conducted at the pressure of the vacuum test chamber, while, the pressurized testing has been conducted at 290 +/- 10 bar (4,200  100 psi) using nitrogen as the pressurizing gas, which corresponds to the design pressure for the target COPV. In this research, spherical aluminum projectiles with varying diameter has been chosen as the impactor. For the unpressurized COPV, the dependence of penetration up to the dependence of hole size in the liner has been obtained as a function of impact conditions. For the pressurized research, the dependence of penetration up to rupture has been obtained as a function of im-pact conditions. Two representative post-test photographs of the failed COPV's from a nor-mal impact into the COPV surface are shown in Fig. 1. These images display the dramatic difference between failure modes, venting (Fig. 1a) and rupture (Fig. 1b). For venting, liner perforation, severed composite fibers/tows and ply delamination are commonly observed damage characteristics of this COPV failure mode. In the case of rupture, the COPV typically experienced a separation of its domed regions and severe break-up of the cylindrical region. Fully understanding the transition from venting to rupture experimentally is costly and potentially unachievable for conditions that cannot be generated in the laboratory. These shortcomings have motivated the performance of three-dimensional numerical simulations to expand the existing experimental database. These simulations have been carried out with the nonlinear-structural-dynamics, analysis-tool, CTH. A typical pressure contour plot from an impact simulation of an entire COPV is shown in Fig. 2. To generate the COPV stress state without initiating a shock wave, the pressure in the simulated COPV is ramped up to the final pressure over a millisecond prior to impact of the projectile with nitrogen gas. Figure 2a shows the system in this initial condition. After one millisecond, a projectile is initiated into the simulation and impacts the COPV. Figure 2b shows the system after this impact. In the figure, the onset of venting is represented as the change in pressure (μbar), red to green, at the perforation site. Also seen in the figure is the eroded projectile that had passed into the COPV vessel with the generated shock wave in the pressurant propagating just ahead of the material. In this paper, pertinent experimental details and the development of the material constitutive models necessary for this work along with the efforts to validate their use are dis-cussed. The simulation results are presented and compared with the NASA experimental observations. While work is on-going from this effort, early observations pertinent to the failure threshold are presented.

Garcia, M. A.↗
DOE OSTI2024

Evaluating Pressure‐dependent Discharge Behavior of Foil Versus In situ Plated Lithium Metal Anodes in Solid‐State Batteries

Abstract Anode‐free manufacturing of solid‐state batteries (SSBs) shows promise to maximize energy density by eliminating excess lithium (Li) and simplifying battery production. However, high reversibility during discharge (stripping of Li) is necessary for long‐lifetime SSBs with a limited Li reservoir. Further, the plastic flow of Li changes depending on the Li thickness, leading to possible differences in discharge performance under stack pressure. This work investigates the pressure‐dependent discharge performance of anode‐free manufactured SSBs with in situ plated Li and compares the performance to that of conventional thick Li foil cells. Distinct stripping behavior is observed at low pressures (0–1 MPa), where Li diffusivity and initial interfacial contact may control accessible capacity, compared to high pressures (3–10 MPa) where mechanical deformation of Li likely governs stripping behavior. Analysis of impedance spectra collected during stripping shows that additional stack pressure delays the formation of deep, as opposed to lateral, voids in the Li anode. These results provide insights to guide the transition from thick Li foil anodes to anode‐free manufactured SSBs.

Haslam, Catherine G.↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗