Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface parameters”

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

Completely Multipolar Model for Many-Body Water–Ion and Ion–Ion Interactions

This work constructs an advanced force field, the Completely Multipolar Model (CMM), to quantitatively reproduce each term of an energy decomposition analysis (EDA) for aqueous solvated alkali metal cations and halide anions and their ion pairings. We find that all individual EDA terms remain well-approximated in the CMM for ion-water and ion-ion interactions, except for polarization, which shows errors due to the partial covalency of ion interactions near their equilibrium. We quantify the onset of the dative bonding regime by examining the change in molecular polarizability and Mayer bond indices as a function of distance, showing that partial covalency manifests by breaking the symmetry of atomic polarizabilities while strongly damping them at short-range. This motivates an environment-dependent atomic polarizability parameter that depends on the strength of the local electric field experienced by the ions to account for strong damping, with anisotropy introduced by atomic multipoles. The resulting CMM model for ions provides accurate dimer surfaces and three-body polarization and charge transfer compared to EDA, and shows excellent performance on various ion benchmarks including vibrational frequencies and cluster geometries.

Heindel, Joseph P↗
DOE OSTI2026

Core plasma fueling by fast inward particle transport after hydrogen pellet injection in Wendelstein 7-X

A large database of more than 1000 individual cryogenic hydrogen pellets injected into Wendelstein 7-X for plasma fueling was analyzed to improve the understanding of the three phases of the process: the ablation, deposition and transport of the pellet material. Kilohertz-sampled electron density and temperature measurements revealed a more complex drift behavior than predicted by numerical code simulation. It could be explained by the poloidal plasma E r x B- drift rotation, which plays a significant role in stellarators, but was not previously considered in pellet injection codes like HPI2. The drift results in a fast poloidal rotation of the pellet material around the plasma core, leading to an almost homogeneous deposition over the involved flux surfaces regardless of magnetic high and low field side injection geometry. Additionally, a novel fast inward directed transport mechanism (‘FIT-effect’) was observed. The effect occurs on timescales of tens of milliseconds and cannot be explained by neoclassical transport or diffusion. It might be linked to the turbulence pinch recently found in Wendelstein 7-X. When the FIT-effect occurs, the pellet particles are rapidly transferred from the deposition flux surfaces to the plasma core, causing the plasma density profile to peak, which is beneficial for confinement in Wendelstein 7-X. The large pellet injection database was statistical analyzed with regard to pellet and plasma parameters, which delivered some starting points towards developing an understanding of the physics behind the FIT-effect. The results indicate, that plasma core fueling via pellet injection is largely independent of the injection geometry in stellarators under certain conditions, reducing the technical complexity of the injection system.

Wendelstein 7-X↗
DOE OSTI2025

Pilot-Scale Validation of Distributed Optical Fiber Sensors for Underground Pipeline Monitoring

Monitoring parameters such as hoop strain, pressure, and acoustic vibrations is key to detecting potential leaks, intrusions, or structural issues. Distributed optical fiber sensor (DOFS) systems provide a compelling solution for continuous, real-time monitoring over long distances. This paper details the development and pilot-scale implementation of DOFS systems for underground pipeline monitoring, evolving from a proof-of-concept stage. Multiple custom-designed DOFS interrogator units—such as optical frequency-domain reflectometry (OFDR), Brillouin optical time-domain analysis (BOTDA), and multimodal interferometer-based fiber acoustic sensor systems were tested to measure the key parameters, such as hoop strain, pipe pressure, surrounding soil temperature, and acoustic vibrations. The underground product pipeline’s outer diameter is 30 inches, the wall thickness is 1.28 inches, and 3 feet deep from the surface. The fiber deployment strategies and sensing data acquisition methods for these systems are discussed. The results demonstrate the effectiveness of DOFS in detecting hoop strain, temperature changes, and acoustic vibrations, showcasing their potential for real-time monitoring and enhancing pipeline safety. These findings from pilot-scale testing offer valuable insights into advancing pipeline monitoring technologies and improving the reliability of underground pipeline systems.

fiber optic sensors↗
DOE OSTI2026

Planetary Boundary-Layer Height (PBLHT) Value-Added Product: Remote-Sensing Retrievals

The planetary boundary layer (PBL) is fundamental to numerous atmospheric processes, including aerosol mixing and transport, cloud evolution, and precipitation formation. A critical parameter in these studies is the PBL height (PBLHT). This vertical depth is essential for characterizing PBL structures in numerical simulations and serves as a primary metric for estimating flux exchanges between the Earth’s surface and the atmosphere. Radiosonde (SONDE) observations provide high-vertical-resolution measurements of temperature and moisture profiles and are widely used to estimate PBLHT (Liu and Liang 2010, Seidel et al. 2010). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s PBLHT value-added product (VAP) for radiosonde measurements, known as PBLHTSONDE, applies three commonly used methods—the Heffter (1980) method, the Liu and Liang (2010) method, and the bulk Richardson number approach (Seibert et al. 2000)—to derive PBLHT. The PBLHTSONDE VAP operates routinely at ARM observatories and mobile facilities, with data available from the ARM Data Center shortly after sounding observations are collected (Sivaraman et al. 2013). However, radiosonde observations are limited by their low temporal resolution. Most stations launch soundings only twice daily, which constrains the ability to investigate and characterize the temporal evolution of the PBL using radiosonde data alone. The use of continuous remote-sensing observations provides high temporal resolution of PBLHT estimates. These observations include aerosol lidars (Dang et al. 2019, Su et al. 2020), Doppler lidar (DL; Tucker et al. 2009, Krishnamurthy et al. 2021), and water vapor and/or temperature lidars and radiometers (Turner et al. 2014). These observations provide valuable data on the PBL’s thermodynamic properties (e.g., water vapor and/or temperature lidars and radiometers), dynamic properties (e.g., DL), and distribution of tracer substances (e.g., aerosol lidars), all of which can be used to estimate PBLHT. ARM developed PBLHT estimates from the micropulse lidar (MPL; PBLHTMPL), Doppler lidar (PBLHTDL), and combined Raman lidar (RL)/atmospheric emitted radiance interferometer (AERI) thermodynamic profiles (PBLHTTHERMO). Each estimate captures different physical characteristics of the boundary layer—aerosol tracers, vertical velocity turbulence, and thermodynamic structure—and exhibits distinct strengths and limitations depending on the PBL regime and time of day. In addition, the ARM ceilometer (CEIL) provides three potential PBLHT candidates derived from the vendor's built-in algorithm. Building on these individual retrievals, ARM developed the PBLHTBEML VAP, which combines the four remote-sensing-based estimates with ancillary meteorological variables using the machine learning approach of Zhang et al. (2025) to produce a best-estimate PBLHT at 10-minute resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Structural and Electrocatalytic Studies of Pulsed Laser Deposited Epitaxial RuO 2 Thin Films

Two sets of high-quality epitaxial ruthenium oxide (RuO 2 ) thin films with different thicknesses were synthesized in situ on cost-competitive sapphire substrates by using a pulsed laser deposition technique. The first set of films, with a thickness of 40 nm and a sheet resistance of 15.7 Ω/□, was prepared using 2100 laser pulses, while the second set of films, with a thickness of 87 nm and a sheet resistance of 6.7 Ω/□, was prepared using 4800 laser pulses. All other deposition parameters were kept the same. The post-deposition structural and morphological measurements showed that both sets of films grew at the same growth rate, had the same crystallinity, similar grain boundary density, and slightly different surface roughness. The thicker RuO 2 films achieved an overpotential of 280 mV for the oxygen evolution reaction at a current density of 100 μA/cm 2 , comparable to or exceeding the performance of films grown on more expensive substrates. Comparatively, the thinner RuO 2 films, which have a significantly higher charge transfer resistance (250 Ω versus 100 Ω for thick films), display a higher overpotential of 320 mV. These results indicate that the lower electrical resistance of thicker films promotes charge transfer through the film body, leading to superior electrocatalytic properties.

charge-transfer↗
DOE OSTI2025

Generation of strong fields with subcritical density plasmas to study the phase transitions of magnetized warm dense matter

Warm dense matter (WDM) is a regime where Fermi degenerate electrons play an important role in the macroscopic properties of a material. Recent experiments have brought us closer to understanding unmagnetized processes in WDM, but magnetized WDM remains unexplored because kilotesla magnetic fields are required. Although there are examples of field compression generating such fields by imploding pre-magnetized targets, these existing methods give no independent control over the parameters of the magnetized plasma and result in limited laser access for sample creation and diagnosis. In this paper, numerical simulations show that kilotesla magnetic fields can be obtained by shining laser beams onto the inner surface of a cylindrical target, rather than on the outer surface. This approach relies on field compression by a low-density, high-temperature plasma, rather than a high-density, low-temperature plasma, used in the more conventional approach. With this novel configuration, the region of peak magnetic field is mostly free of plasma, hence, other beams can reach a sample placed in the region of the peak field to form WDM and diagnose it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Molecular Dynamics Study of the Effect of Grafting Density on Ion Diffusivity in a MARTINI Coarse-Grained Strong Polyelectrolyte Brush

Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separation technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact the PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and the counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI force field. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N = 150 monomers; M W = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm 2 . Using 5 μs simulations, we investigate the effects of grafting density on the PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. Finally, the interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning Framework for Characterizing Processing–Structure Relationship in Block Copolymer Thin Films

The morphology of block copolymers (BCPs) critically influences material properties and applications. This work introduces a machine learning (ML)-enabled, high-throughput framework for analyzing grazing incidence small-angle X-ray scattering (GISAXS) data and atomic force microscopy (AFM) images to characterize BCP thin film morphology. A convolutional neural network was trained to classify AFM images by surface features, achieving 97% testing accuracy. Classified images were then analyzed to extract 2D grain size measurements from the samples in a high-throughput manner. ML models were trained to predict domain orientation based on processing parameters such as solvent ratio, additive type, and additive ratio. GISAXS-based properties were predicted with strong performances (R 2 > 0.75), while AFM-based property predictions were less accurate (R 2 < 0.60), likely due to the localized nature of AFM measurements compared to the bulk information captured by GISAXS. Beyond model performance, interpretability was addressed using SHapley Additive exPlanations (SHAP). SHAP analysis revealed that the additive ratio had the largest impact on morphological predictions, where additive provides the BCP chains with increased volume to rearrange into thermodynamically favorable morphologies. This interpretability helps validate model predictions and offers insight into parameter importance. Altogether, the presented framework combining high-throughput characterization and interpretable ML offers an approach to exploring and optimizing BCP thin film morphology across a broad processing landscape.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Characterization of volcanic tuff pores pre- and post-underground nuclear detonation using ultra-small and small angle neutron scattering

The ability to accurately model the subsurface transport of radionuclides is fundamental to the remote detection and characterization of underground nuclear explosion (UNE) events. Developing more sophisticated transport models presents a significant opportunity to enhance monitoring capabilities, particularly in the reliable prediction of signature migration. Experimentally determined characterization of geologic materials associated with transport properties is the pertinent base information for such robust model development and calibration. Here, we report results from an unprecedented study demonstrating changes to the pore and fracture network structures in geological materials in response to UNEs over nanometer to micrometer length scales. Volcanic tuffs of five different lithological formations from pre- and post-UNE environments were collected from the Nevada National Security Site. Combined ultra-small and small-angle neutron scattering techniques were used to characterize the tuff pore structure. The results demonstrate measurable differences in the specific surface area and porosity of samples pre- and post-shot from texturally similar lithological formations, indicating that pore properties can serve as a direct physical signature of a UNE. The results also provide experimentally determined transport parameters in support of advanced model development through the integration of gas migration, hydrodynamic simulations, and geologic framework models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Predictive numerical modeling of plasma-induced surface roughness and wettability evolution in LM-PAEK/CF tape

Plasma surface modification effectively enhances adhesion in thermoplastic composites, yet its impacts on high-performance polymers like low-melting polyaryletherketone (LM-PAEK) remain inadequately quantified. Here, this study integrates experimental analysis and numerical modeling to characterize surface roughness and wettability changes in LM-PAEK/carbon fiber composites treated with atmospheric plasma. Atomic Force Microscopy quantified surface topography (n = 10 per condition for contact angles), while static contact-angle assessments measured wettability. Roughness rapidly increased from ∼0.2 nm to 1.6 nm, and contact angle reduced from ∼90° to 24°, both stabilizing after 25–30 s of exposure. A semi-empirical, physics-informed framework was calibrated to these data, coupling surface chemistry via the Owens–Wendt decomposition with topography via the Wenzel roughness factor, and evaluated using out-of-sample (cross-validated) tests, while static contact angle assessments measured wettability. Numerical predictions matched experimental results closely (RMSE <5%, R 2 > 0.95). Incorporating material-specific parameters, the calibrated model supports plasma-treatment optimization and provides quantitative guidance for improving interfacial adhesion in thermoplastic composite manufacturing.

Atomic Force Microscopy↗
DOE OSTI2024

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Effect of Ni deposition on CrI 3 monolayer: First-principles study

In our recent experimental work, we demonstrated that Ni thin film deposition can be used to tune the surface magnetic states in bulk CrI 3 . Here, we present the results of a computational study of Ni deposition on a CrI 3 monolayer. In particular, we compare structural, electronic, and magnetic properties of three different systems: a pristine CrI 3 monolayer, a CrI 3 monolayer with partial Ni coverage, and a CrI 3 monolayer with full Ni coverage. The corresponding cells used in our calculations are Cr 2 I 6 , Cr 2 I 6 Ni 1 , and Cr 2 I 6 Ni 3 . In Cr 2 I 6 Ni 1 , the Ni atom is positioned over the hollow site, which corresponds to the lowest energy configuration. In Cr 2 I 6 Ni 3 , the three Ni atoms are positioned over the single hollow site, and over the two Cr sites. Our calculations indicate that the partial coverage of CrI 3 monolayer with Ni retains the semiconducting nature of the parent compound, but results in a strong reduction of the energy band gap. At the same time, the full coverage of CrI 3 monolayer with Ni results in a metallic transition, with a highly spin-polarized (possibly half-metallic) electronic structure. The magnetic alignment of the CrI 3 monolayer with the deposited Ni is ferrimagnetic, due to the anti-aligned magnetic moments of Cr and Ni atoms. While the partial coverage of the CrI 3 monolayer with Ni has a moderate impact on the lattice parameter of the parent compound, full coverage by Ni results in a significant increase of the lattice constant. For both partial and full coverages, the lowest energy configuration corresponds to the Ni atoms being approximately aligned with the out-of-plane coordinate of the surface iodine layer. In addition, in the full coverage mode, the Cr layer exhibits a significant shift away from the Ni atoms, effectively “detaching” the Ni / I atomic layer from the rest of the cell.

2D magnetism↗
DOE OSTI2025

Dynamic Carrier Modulation via Nonlinear Acoustoelectric Transport in van der Waals Heterostructures

Dynamically manipulating carriers in van der Waals heterostructures could enable solid-state quantum simulators with tunable lattice parameters. A key requirement is the formation of deep potential wells to reliably trap excitations. Here, we report the observation of nonlinear acoustoelectric transport and dynamic carrier modulation in boron nitride-encapsulated graphene devices coupled to intense surface acoustic waves (SAWs) on LiNbO 3 substrates. SAWs generate strong acoustoelectric current densities ( J AE ), transitioning from linear to nonlinear regimes with increasing SAW intensity. In the nonlinear regime, periodic carrier (electrons, holes, or their mixtures) stripes emerge. Using counter-propagating SAWs, we create standing SAWs (SSAWs) to dynamically manipulate charge distributions without static gates. The saturation of J AE , attenuation transitions, and tunable resistance peaks confirms strong carrier localization. Finally, these results establish SAWs as a powerful tool for controlling carrier dynamics in two-dimensional (2D) materials, paving the way for the development of time-dependent quantum systems and acoustic lattices for quantum simulation.

2D materials↗
DOE OSTI2025

Quantum Sensing in Micro-Architected Scaffolds

Quantum sensing with nitrogen-vacancy centers in diamond has emerged as a powerful tool for measuring diverse physical parameters, yet the versatility of these measurement approaches is often limited by the achievable layout and dimensionality of bulk-crystal platforms. Here, in this study, we demonstrate a versatile approach to creating designer quantum sensors by surface-functionalizing multiphoton lithography microstructures with NV-containing nanodiamonds. We showcase this capability by fabricating a 150 μm × 150 μm × 150 μm triply periodic minimal surface gyroid structure with millions of attached nanodiamonds. We demonstrate a means to volumetrically image these structures using a refractive index matching confocal imaging technique and extract ODMR spectra from 1.86 μm × 1.86 μm areas of highly concentrated nanodiamonds across a cross-section of the gyroid. Furthermore, the high density of sensing elements enables ensemble temperature measurements with a sensitivity of 0.548 ± 0.084 K/√Hz at 5 mW excitation power. This approach to creating quantum-enabled microarchitectures opens new possibilities for multimodal sensing in complex three-dimensional environments.

Materials science↗
DOE OSTI2026

Quantitative Analysis of the Semiconductor–Electrolyte Interface Using Cyclic Voltammetry Measurements

Small changes in the chemical potential at a semiconductor interface can result in dramatic changes to the space-charge layer that underpins applications in the electronic and photovoltaic industries as well as in photoelectrochemical cells for fuel production. There has hence been great interest in techniques that directly probe the space-charge layer, yet many fail at the semiconductor–electrolyte interface due to the potential drop in the electric double-layer region of the electrolyte. This article demonstrates that photovoltages, obtained from straightforward cyclic voltammetry measurements, provide an experimental and quantitative approach for characterizing the semiconductor–electrolyte interface. Key parameters accessible through this approach include the flat-band potential ( E fb ), the fraction of the total potential that drops across the space-charge layer (γ sc ) and the electric double layer, as well as the surface recombination lifetime (τ s ). Here, we report photovoltage measurements for p -type Si(111) photoelectrodes in contact with electrolytes containing redox-active species with a range of known reduction potentials that exceed the 1.1 eV bandgap. In tetrabutylammonium [NBu 4 ] + electrolyte, the flat-band potential determined for hydrogen-terminated ( p -Si–H), methyl-terminated ( p -Si–CH 3 ), and chemically oxidized ( p -Si–cSiO x ) surfaces were −0.02, −0.31, and 0.30 V vs Fc +/0 , respectively, agreeing well with expected shifts arising from surface dipole modifications. The quantitative analysis also reveals that 67% of the applied bias drops across the space-charge layer for p -Si–H, 73% for p -Si–CH 3 , and only 44% for p -Si–cSiO x . The remaining potential drop is attributed to the interfacial surface layer, which consists of a molecular dipole or oxide overlayer, and the Helmholtz layer within the electrolyte. When the larger [NBu 4 ] + electrolyte was replaced with Li + , the flat-band position showed minimal changes, but the fraction of the potential drop across the space-charge layer increased significantly, consistent with the small cation altering the structure of the electric double layer.

electrolytes↗
DOE OSTI2026

Impact of n > 1 neoclassical tearing modes on fast ions, plasma rotation, and the onset of disruptive tearing modes in DIII-D ITER baseline scenario

We present integrated TRANSP-kick analyses quantitatively characterizing how n > 1 magnetic islands redistribute fast ions (FIs) and modify momentum balance in low-torque DIII-D H-mode scenarios operated with the ITER normalized parameter set and shape (ITER baseline scenario, or IBS). In this plasma scenario, disruptive neoclassical tearing modes (NTMs) are most commonly seeded by nonlinear 3-wave coupling when the differential rotation between the q = 1 and q = 2 surfaces ($δf_{1, 2}$) approaches zero. Analysis of a DIII-D IBS database of shots unstable to the 2/1 NTM demonstrates that the flattening rotation profile is correlated with n > 1 NTM amplitude growth, with a strong drop in $δf_{1, 2}$ occurring when the radial magnetic perturbation amplitude at the rational surface ($δB$ r s ) is estimated to be about 20 G. The interpretive TRANSP-kick simulations show that the experimentally observed reduction in core rotation is due to the fast-ion kicks. Nonlinearity in momentum losses versus NTM amplitude is found to be correlated with overlap of FI island chains, occurring when $δB$ r s is about 35 G, in rough agreement with the experimental results. This agreement suggests the primary mechanism for the conductive momentum loss in the core is due to chaotic FI orbits and FI losses developing from overlapping resonant FI island chains.

ITER baseline scenario↗
DOE OSTI2026

Time-, space-, and molecule-resolved multiple plasma reactive species in complex settings using an innovative approach of cavity ringdown spectroscopy (Project Final Report)

The primary research goal of this 2-year project is to demonstrate the feasibility of an innovative approach to measure reactive plasma species in the most challenging settings: right on the surface of water and under water. We propose to use the Brewster angle-cavity ringdown spextroscopy (BA-CRDS). Research objectives include: 1) Measurements of OH and NO in the plasma gas-phase zone (Z1) using different types of LTP sources including continuous and pulsed plasma jets, sheets, torches, etc. 2) Measurements of OH and NO in the multiphase reaction zones, right on the surface of water (Z2) and under water (Z3). 3) RS formation and loss mechanisms and reaction dynamics in the three zones. 4) Generalization of the proposed system to be configurable with other selected innovative techniques for other LTP parameters; invitation of selected PIs to form a strong collaborative team to pursue the follow-up effort (Category-2).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

FIREFLY: heat load and particle exhaust approximations for rapid evaluation of divertor designs

The divertor in a magnetic confinement fusion reactor is an essential component for power dissipation and particle removal. The FIREFLY package for rapid evaluation of divertor designs is presented as an extension of the FLARE code for field line reconstruction from a flux tube mesh. First, divertor loads are approximated with a simplified heat transport model. Neutralized particles are then sampled from the resulting load distribution, and the EIRENE code is used to track molecules and atoms in a plasma background while accounting for dissociation, charge exchange and ionization. Particles are removed on pumping surfaces in order to estimate the exhaust efficiency for a given divertor geometry. Optimization of the divertor geometry for more efficient particle exhaust is explored by using W7-X as an example, and the sensitivity to model parameters for the plasma background in the proxy calculations is evaluated.

mesh generation, magnetic field lines, scrape-off ↗