SEARCH · Engineering Papers

Results for “Material balance time”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Nickel-rich Nickel Manganese Cobalt (NMC622) Cathode Lithiation Mechanism and Extended Cycling Effects Using Operando X-ray Absorption Spectroscopy

Ni-rich NMC materials are a particularly promising class of Li-ion cathodes for various applications. LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) offers a unique balance of thermal stability and energy density, thus attracting attention for electric vehicle implementation. However, upon extended cycling, capacity fade is prevalent due to structural degradation, which is a major drawback for layered oxide cathodes. Therefore, exploring the underlying phenomena that drive detrimental structural response can lead to future improvements. For the first time, operando X-ray absorption spectroscopy (XAS) was performed on NMC622 pouch cells at three different stages. An extensive description of the first cycle (de)lithiation mechanisms was achieved through X-ray absorption near-edge structure analyses and extended X-ray absorption fine structure modeling. Transition metal specific electrochemical participation and structural variation revealed that much of the delivered capacity and distortion is a result of Ni redox behavior, while the local structure of Co and Mn are impacted due to their interdependencies. Key mechanistic components were identified, as the local structural variation from redox processes and Ni 3+ Jahn–Teller distortion were decoupled. Further, operando XAS was used to investigate the structural response of NMC622 to extended cycling and was supported by X-ray diffraction, X-ray photoelectron spectroscopy, and electrochemical impedance spectroscopy analyses. In this work, reduced capacity was found in the cell after 100 cycles and is attributed to structural degradation and cathode–electrolyte interphase buildup from repeated Li (de)insertion processes, which limit the electrochemical reversibility of Ni and Co through increased polarization. These results expand the understanding of a Ni-rich NMC material under extended cycling, which is vital to the future design of electrode materials.

25 ENERGY STORAGE↗

DOE OSTI2023

Dynamic study of direct CO 2 capture from indoor air using poly(ethylenimine)-impregnated fiber sorbents

Supported amine adsorbents are promising materials for direct air capture (DAC) of CO 2 due to their high CO 2 capacity and relatively low energy requirement for regeneration. For a DAC process, it is essential to properly define operating parameters to achieve high sorbent productivity (amount of CO 2 captured per unit quantity of sorbent material over unit time). It is furthermore essential to understand the kinetic behavior of the process under the influence of various operating conditions such as the inlet air velocity, sorbent composition, and humidity to select an effective range of operating conditions to maximize sorbent productivity. Here, the dynamic behavior of a DAC process is probed using a fixed fiber sorbent contactor containing poly(ethylenimine) (PEI)-impregnated composite silica/cellulose acetate (CA) fibers. Throughout this study, experiments are conducted using both simulated air (398 ppm CO 2 balanced by N 2 ) and real indoor air (~400–500 ppm CO 2 ). The experimental behavior of the fibers using simulated air and indoor air is compared, and the influence of the inlet air velocity on the breakthrough behavior is assessed. By changing operating conditions, the impact on the fiber sorbent productivity (mmol CO 2 g fiber –1 h –1 ) is quantified to identify conditions that could favor high rates of CO 2 removal. The kinetics of steam-assisted CO 2 desorption are studied, identifying achievable desorption times. Productivities of 1.2 mmol CO 2 g fiber –1 h –1 are obtained using an inlet air velocity of 1.1 m s –1 . Performance trends show that further increasing the inlet air velocity will likely lead to even higher productivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗

DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗

DOE OSTI2022

Hydrogen Storage in Partially Exfoliated Magnesium Diboride Multilayers

Metal boride nanostructures have shown significant promise for hydrogen storage applications. However, the synthesis of nanoscale metal boride particles is challenging because of their high surface energy, strong inter- and intraplanar bonding, and difficult-to-control surface termination. Here, it is demonstrated that mechanochemical exfoliation of magnesium diboride in zirconia produces 3-4 nm ultrathin MgB 2 nanosheets (multilayers) in high yield. High-pressure hydrogenation of these multilayers at 70 MPa and 330 °C followed by dehydrogenation at 390 °C reveals a hydrogen capacity of 5.1 wt%, which is ≈50 times larger than the capacity of bulk MgB 2 under the same conditions. This enhancement is attributed to the creation of defective sites by ball-milling and incomplete Mg surface coverage in MgB 2 multilayers, which disrupts the stable boron-boron ring structure. The density functional theory calculations indicate that the balance of Mg on the MgB 2 nanosheet surface changes as the material hydrogenates, as it is energetically favorable to trade a small number of Mg vacancies in Mg(BH 4 ) 2 for greater Mg coverage on the MgB 2 surface. To conclude, the exfoliation and creation of ultrathin layers is a promising new direction for 2D metal boride/borohydride research with the potential to achieve high-capacity reversible hydrogen storage at more moderate pressures and temperatures.

08 HYDROGEN↗

DOE OSTI2025

Simulation-driven design optimization of reaction injection molding (RIM) process for polydicyclopentadiene (pDCPD): Minimizing cycle time, defects, and warpage

Replacing metal components in trucks, trailers, and buses with lightweight polymer composites is challenging due to high temperatures and complex manufacturing. The Reaction Injection Molding (RIM) process using Dicyclopentadiene (DCPD) resin offers a solution by producing robust parts with excellent stiffness, impact strength, and resistance properties. Simulations are essential for optimizing this process, predicting defects, and improving quality. However, most commercial software is tailored for thermoplastics, requiring thermoset users to generate their own datasets. In this study, a material data card for DCPD was developed to perform RIM simulations. Design of Experiments (DOE) was used to identify key factors affecting filling, curing, and warpage, aiming to minimize cycle time and defects. The simulations explored varying injection gate parameters (size, location, number) and process conditions (mold/resin temperature, injection/curing pressure). Results showed that gate design significantly impacts filling behavior and defects. A single central gate provided balanced flow with fewer defects, while two corner gates led to more defects. Additionally, lower injection pressure increased filling time, while higher mold temperature accelerated curing but led to more warpage. In conclusion, this optimization framework aims to enhance DCPD part performance and promote sustainable manufacturing by reducing waste and energy consumption.

42 ENGINEERING↗

DOE OSTI2025

Dynamic Competition between Hubbard and Superexchange Interactions Selectively Localizes Electrons and Holes through Polarons

Controlling the effects of photoexcited polarons in transition metal oxides can enable the long time-scale charge separation necessary for renewable energy applications and controlling new quantum phases through dynamically tunable electron–phonon coupling. In previously studied transition metal oxides, polaron formation is facilitated by a photoexcited ligand-to-metal charge transfer (LMCT). When the polaron is formed, oxygen atoms move away from iron centers, which increases carrier localization at the metal center and decreases charge hopping. Studies of yttrium iron garnet and erbium iron oxide have suggested that strong electron and spin correlations can modulate photoexcited polaron formation. To understand the interplay between strong spin and electronic correlations in highly polar materials, we studied gadolinium iron oxide (GdFeO 3 ), which selectively forms photoexcited polarons through an Fe–O–Fe superexchange interaction. Excitation-wavelength-dependent transient extreme ultraviolet (XUV) spectroscopy selectively excites LMCT and metal-to-metal charge transfer (MMCT) transitions. Here, the LMCT transition suppresses photoexcited polaron formation due to the balance between superexchange and Hubbard interactions, while MMCT transitions result in photoexcited polaron formation within 250 ± 40 fs. Ab initio theory demonstrates that electron and hole polarons localize on iron centers following MMCT. In addition to understanding how strong electronic and spin correlations can control strong electron–phonon coupling, these experiments separately measure electron and hole polaron interactions on neighboring metal centers for the first time, providing insight into a large range of charge-transfer and Mott–Hubbard insulators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Water limiting current measurements in anion exchange membrane fuel cells (AEMFCs); part 1: Water limiting current method development

We report water management is critical in optimizing anion exchange membrane fuel cell (AEMFCs) performance and durability. The role of membranes on determining cell water balance in AEMFCs has received limited study. A more thorough understanding of membrane water transport parameters is necessary for thorough cell optimization. In this study, a novel method for probing limiting current based on water flux across an AEMFC is presented for the first time. The water limiting current has been investigated as a function of relative humidity of anode and cathode streams, cathode backpressure and cell operating temperature. From these studies, it is clearly shown that water diffusion through the membrane from the anode is the most critical mechanism for achieving high current density performance and durability. The developed method can act as a novel screening technique that focuses on the impacts of water transport parameters and helps identify promising membrane materials based on their water transport properties. A critical but underappreciated facet of AEMFC performance and durability.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2021

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗

DOE OSTI2022

Pareto Optimization of Oligomer Polarizability and Dipole Moment Using a Genetic Algorithm

High-performance electronic components are highly sought after in order to produce increasingly smaller and cheaper electronic devices. Drawing inspiration from inorganic dielectric materials, in which both polarizability and polarization contribute, organic materials can also maximize both. For a large set of small molecules drawn from PubChem, a Pareto-like front appears between the polarizability and dipole moment, indicating the presence of an apparent trade-off between these two properties. We tested this balance in π-conjugated materials by searching for novel conjugated hexamers with simultaneously large polar- izabilities and dipole moments with potential use for dielectric materials. Using a genetic algorithm (GA) screening technique in conjunction with an approximate density functional tight-binding method for property calculations, we were able to efficiently search chemical space for optimal hexamers. Given the scope of chemical space, using the GA technique saves considerable time and resources by speeding up molecular searches compared to a systematic search. Here, we also explored the underlying structure–function relationships, including sequence and monomer properties, that characterize large polarizability and dipole moment regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Three rate-determining protein roles in photosynthetic O 2 -evolution addressed by time-resolved experiments on genetically modified photosystems

Light-driven water splitting by plants, algae and cyanobacteria is pivotal for global bioenergetics and biomass formation. A manganese cluster bound to the photosystem II proteins catalyzes the complex reaction at high rate, but the rate-determining factors are insufficiently understood. Here we trace the oxygen-evolution transition by time-resolved polarography and infrared spectroscopy for cyanobacterial photosystems genetically modified at two strategic sites, complemented by computational chemistry. Our results highlight three rate-determining roles of the protein environment of the metal cluster: acceleration of proton-coupled electron transfer, acceleration of substrate-water insertion after O 2 -formation, and balancing of rate-determining enthalpic and entropic contributions. Whereas in general the substrate-water insertion step may be unresolvable in time-resolved experiments, here it likely becomes traceable because of deceleration by genetic modification. Our results may stimulate new time-resolved experiments on substrate-water insertion in photosynthesis, clarification of enthalpy-entropy compensation in enzyme catalysis, and knowledge-guided development of inorganic catalyst materials.

Bioenergetics↗

DOE OSTI2020

Catalytic Bio-crude Production in a Novel, Short-Contact Time Reactor

RTI has developed a catalytic fast pyrolysis (CFP) process to produce biomass-based hydrocarbon replacement fuels. From a technology perspective, this advanced biofuels process produces liquid transportation fuels that can leverage capital expenditures in the existing petroleum distribution infrastructure. From a business perspective, this technology does not face potential market limitations caused by oxygenate blend limitations for light duty vehicle applications. A promising catalyst has been developed with excellent potential for deoxygenation to produce a hydrocarbon-rich biocrude that is more suitable for upgrading with traditional hydroprocessing technology. The role of the catalyst is to control the chemistry during biomass pyrolysis to produce a biocrude in a single step that has lower oxygen content and is more thermally stable than conventional biomass fast pyrolysis oil. Catalyst properties are optimized to minimize gas and coke production and improve catalytic deoxygenation and biocrude yields. Results from this project have demonstrated that biocrudes from catalytic fast pyrolysis contain <20 wt% oxygen, with the organic fraction as low as 12% oxygen. These oxygen numbers are > 50% lower compared to traditional fast pyrolysis oils. The biocrude produced from our CFP process can be upgraded to gasoline range hydrocarbons. A very conservative approach for upgrading included a preliminary stabilization step followed by a hydroprocessing step. The oxygen content in the upgraded product was decreased by 97% (to less than 3wt %) from the oxygen content of the biocrude sample. Additionally, the acid content was reduced by 99% from the acid content of the biocrude sample. In an optimized integrated process, the CFP step and the hydroprocessing (upgrading) step need to be carefully balanced to maximize the input biomass energy that is recovered in the finished biofuels. Hydrogen utilization and carbon recovery during both catalytic biomass pyrolysis and biocrude upgrading need to be optimized to maximize energy recovery in the finished biofuel. A nominal 1 ton per day (TPD) catalytic biomass pyrolysis system was designed based on the laboratory experimental data. The system was designed to achieve the short residence times (0.5–2 sec) and high heat transfer rates required for maximum liquid biocrude yields while optimizing process integration to maintain catalyst activity by continuous regeneration. The catalytic biomass pyrolysis reactor is a continuously circulating single-loop transport reactor design that is flexible enough to allow sensitivity studies around temperature, residence time, biomass feed rate, and catalyst-to-biomass ratio (i.e., catalyst circulation rate) for process optimization. For a nominal 100 lb/hr feed rate of dry biomass, the products include ~25 lb/hr char and coke, ~20 lb/hr permanent gases, and ~55 lb/hr condensable vapors. Pyrolysis temperatures are 350–600°C, and the system is operated close to ambient pressure. Production of the catalyst developed within the project was scaled up to provide enough material (200kg) for bench-scale testing. The 1TPD bench-scale unit was fabricated, installed and commissioned with this project. Biocrude samples from wood were produced in this unit with the CFP catalyst with material balance > 90%. The intermediate biocrude yields and chemical composition compare well with previous laboratory results and a hydroprocessing strategy for upgrading the resulting biocrude was developed. The experimental results obtained within this project were used to develop a process model for a 2000 TPD commercial facility. Detailed flow sheets were developed for all major unit operations in the process and provided a basis for an economic analysis to develop capital and operating costs. Estimates of the biofuel production costs were used to develop internal rates of return for the yields and compositions measured for the current state of technology at the end of this project and for the yields anticipated for an optimized (nth plant) scenario as additional process improvements are made.

09 BIOMASS FUELS↗

DOE OSTI2020

Task 12 PV Sustainability - Life Cycle Inventories and Life Cycle Assessments of Photovoltaic Systems

Life Cycle Assessment (LCA) is a structured, comprehensive method of quantifying material- and energy-flows and their associated impacts in the life cycles of products (i.e., goods and services). One of the major goals of IEA PVPS Task 12 is to provide guidance on assuring consistency, balance, transparency and quality of LCA to enhance the credibility and reliability of the results. The current report presents the latest consensus life cycle inventories among the authors, PV LCA experts in North America, Europe, Asia and Australia. At this time consensus is limited to four technologies for which there are well-established and up-to-date life cycle inventory (LCI) data (mono- and multi-crystalline Si, CdTe, CIGS, as well as one emerging technology (perovskite silicon tandem). LCIs are necessary for LCA and the availability of such data is often the greatest barrier for conducting LCA. The Task 12 LCA experts have put great efforts in gathering and compiling the LCI data presented in this report. These include detailed inputs and outputs during manufacturing of cell, wafer, module, and balance-of-system (i.e., structural and electrical components) that were estimated from actual production and operation facilities. In addition, data are presented to enable analyses of various types of PV installations; these include operational data of rooftop and ground-mount PV systems and country-specific PV-mixes. The LCI datasets presented in this report are the latest that are available to the public describing the status in 2018 for crystalline Si (some manufacturing data from 2011 were not updated), 2015 and 2017-2018 for CdTe, 2010 for CIGS, 2010 for HCPV, and 2017 for perovskite silicon tandem technology.

14 SOLAR ENERGY↗

DOE OSTI2022

Radiation hard gallium oxide scintillators for high count rate radiation detection

Despite decades of materials research, the availability of appropriate crystalline radiation-hard scintillators with ultrafast decay times (~10ns) and high light yields (>2000 ph/MeV) are still limited. In this study, we demonstrate the potential of gallium oxide (β-Ga 2 O 3 ) based scintillators for high count rate applications. The low-cost rugged β-Ga 2 O 3 scintillators were grown using the optical float zone (FZ) technique. Several dopants and growth atmospheres were used to demonstrate the balance between ultrafast primary decay time and the light yield of the scintillators. Light yields as high as 6446 ± 716 ph/MeV were obtained with 18.3 ns decay time for 662 keV gamma rays. Decay times as fast as 12ns were obtained with a 3212 ± 357 ph/MeV light yield. An excellent 662 keV gamma energy resolution of 7% was obtained using the β-Ga 2 O 3 :Ce crystals. Even with these excellent scintillation properties, the gamma radiation hardness of β-Ga 2 O 3 crystals was better than the leading radiation-hard lead tungstate (PbWO 4 ) crystals. Twenty-six scintillators fabricated from the grown β-Ga 2 O 3 crystals were tested for radiation hardness with a 60 Co gamma irradiation rate of 14 krad/h with a total radiation dose of 2 Mrad. The scintillation properties of these crystals remained unaffected. To date, no other scintillator has demonstrated such high radiation resistance. Furthermore, these FZ-grown low-cost, non-hygroscopic, and high-density β-Ga 2 O 3 scintillators can replace the existing scintillators in many applications ranging from nuclear and high energy physics experiments to nuclear security.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗

DOE OSTI2021

Humate's Role in Terrestrial Greenhouse Gas and Water Reduction

Carbon dioxide (CO 2 ) is considered the sole culprit for global warming; however, nitrous oxide (N 2 O), a greenhouse gas (GHG) with approximately 300 times more global warming potential than CO 2 , accounts for 6% of the GHG emissions in the United States. Seventy five percent of N 2 O emissions come from synthetic nitrogen (N) fertilizer usage in the agriculture sector primarily due to excess fertilization. Numerous studies have shown that changes in soil management practices, specifically optimizing N fertilizer use and amending soil with organic and humate materials can reverse soil damage and improve a farmer's or land reclamation company's balance sheet. Soil restoration is internationally recognized as one of the lowest cost GHG abatement opportunities available. Profitability improves in two ways: (1) lower operating costs resulting from lower input costs (water and fertilizer); and (2) increased revenue by participation in emerging GHG offsets markets, and water quality trading markets.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modeling of ExB effects on tungsten re-deposition and transport in the DIII-D divertor

Mixed-material DIVIMP-WallDYN modelling, now incorporating ExB drifts, is presented that simultaneously reproduces tungsten (W) erosion and deposition patterns observed during the DIII-D Metal Rings Campaign, in which a toroidally symmetric set of W-coated tiles were installed in the carbon (C) DIII-D divertor. Since most reactor plasma facing component (PFC) designs call for mixed-material environments, including ITER’s W/Be enviroment, the divertor targets will quickly evolve into reconstitued surfaces of multiple elements. This work identifies controlling physics that affects material migration patterns in the divertor, which impact PFC lifetimes and impurity leakage from the divertor to the core. These simulations indicate that radial and poloidal ExB transport dominates over parallel force balance for high-Z impurities such as W in the divertor region of DIII-D. It is demonstrated that ExB drifts are required to reproduce the experimental observation of non-local W and C co-accumulation in a band ~7-9 cm outboard of the outer-strike-point (OSP) W source, for attached Lmode conditions in the unfavorable ion grad-B drift direction. In addition, W gross erosion is localized to the region outboard of the OSP, as the formation of C co-deposits suppresses W erosion at the strike point. Time-dependent simulations with scaled ExB impurity drifts (60% of the OEDGE-calculated drift velocity) and W re-erosion quantitatively reproduce these features, including depth-resolved W/C ratios, within a factor of 2 over ~115 seconds of accumulated plasma exposure. The location of co-deposition regions is shown to be well represented by an analytic leakage model, driven largely by poloidal ExB drifts. Qualitative agreement is also found between campaign-integrated W deposition measurements and simulations for the favorable ion grad-B drift direction, the standard mode of operation for most tokamaks. Furthermore, these results imply that a longterm inward radial migration of material from the outer divertor through the private flux region may occur in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗