Search NASA⌕ Search

SEARCH · Search NASA

Results for “inversion”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part I-Algorithm and Morphology

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in characterizing the global 3D distribution and variability of ionospheric electron density (N e ). In this study, we developed an optimal estimation (OE) method to retrieve N e profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles (ε < 0°). Although both OE and onion-peeling (OP) methods use the Abel weighting function in the N e inversion, they are significantly different in terms of performance in the lower ionosphere. The new OE results can overcome the large N e oscillations, sometimes negative values, seen in the OP retrievals in the E-region ionosphere. In the companion paper in this Special Issue, the HmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006–2019), COSMIC-2 (2019–present), and Spire (2019–present) constellations, showing a consistent ionospheric N e morphology. The unprecedented spatiotemporal sampling of the ionosphere from these constellations now allows a detailed analysis of the frequency–wavenumber spectra for the N e variability at different heights. In the lower ionosphere (~150 km), we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2, and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4, and SW4. The co-existence of eastward- and westward-propagating wave4 components implies the presence of a stationary wave4 (SPW4), as suggested by other earlier studies. Further improvements to the OE method are proposed, including a tomographic inversion technique that leverages the asymmetric sampling about the tangent point associated with GNSS-LEO links.

GNSS constellation↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2024

Comparison of machine learning and electrical resistivity arrays to inverse modeling for locating and characterizing subsurface targets

Here, this study evaluates the performance of multiple machine learning (ML) algorithms and electrical resistivity (ER) arrays for inversion with comparison to a conventional Gauss-Newton numerical inversion method. Four different ML models and four arrays were used for the estimation of only six variables for locating and characterizing hypothetical subsurface targets. The combination of dipole-dipole with Multilayer Perceptron Neural Network (MLP-NN) had the highest accuracy. Evaluation showed that both MLP-NN and Gauss-Newton methods performed well for estimating the matrix resistivity while target resistivity accuracy was lower, and MLP-NN produced sharper contrast at target boundaries for the field and hypothetical data. Both methods exhibited comparable target characterization performance, whereas MLP-NN had increased accuracy compared to Gauss-Newton in prediction of target width and height, which was attributed to numerical smoothing present in the Gauss-Newton approach. MLP-NN was also applied to a field dataset acquired at U.S. DOE Hanford site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Insights into the Surface Electronic Structure and Catalytic Activity of InO x /Au(111) Inverse Catalysts for CO 2 Hydrogenation to Methanol

In this article, the direct conversion of carbon dioxide (CO 2 ) into methanol (CH 3 OH) via low-temperature hydrogenation is crucial for recycling anthropogenic CO 2 emissions and producing fuels or high value chemicals. Nevertheless, it continues to be a great challenge due to the trade-off between selectivity and catalytic activity. For CO 2 hydrogenation, In 2 O 3 catalysts are known for their high CH 3 OH selectivity. Subsequent studies explored depositing metals on In 2 O 3 to enhance CO 2 conversion. Despite extensive research on metal (M) supported In 2 O 3 catalysts, the role of In-M alloys and M/In 2 O 3 interfaces in CO 2 activation and CH 3 OH selectivity remains unclear. In this work, we have examined the behavior of In/Au(111) alloys and InO x /Au(111) inverse systems during CO 2 hydrogenation using synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS) and catalytic tests in a batch reactor. Indium forms alloys with Au(111) after deposition. The In-Au(111) alloys display high reactivity towards CO 2 and can dissociate the molecule at room temperature to generate InO x nanostructures. At very low coverages of In (≤ 0.05 ML), the InO x nanostructures are not stable under CO 2 hydrogenation conditions and the active In-Au(111) alloys produces mainly CO and little methanol. An increase in indium coverage to 0.3 ML led to stable InOx nanostructures under CO 2 hydrogenation conditions. These InO x /Au(111) catalysts displayed a high selectivity (~ 80 %) towards CH 3 OH production and an activity for CO 2 conversion that was at least 10 times larger than that of plain In 2 O 3 or Cu(111) and Cu/ZnO(000$\overline{1)}$ benchmark catalysts. The results of AP-XPS show that InO x /Au(111) produces methanol via methoxy intermediates. Inverse oxide/metal catalysts containing InOx open up a possibility for improving CO 2 → CH 3 OH conversion in processes associated with the control of environmental pollution and the production of high value chemicals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Local Inversion Symmetry Breaking and Thermodynamic Evidence for Ferrimagnetism in Fe 3 GaTe 2

The layered compound Fe 3 GaTe 2 is attracting attention due to its high Curie temperature, low dimensionality, and the presence of topological spin textures above room temperature, making Fe 3 GaTe 2 a good candidate for applications in spintronics. Here, in this study, we show, through transmission electron microscopy (TEM) techniques, that Fe 3 GaTe 2 single crystals break local inversion symmetry while maintaining global inversion symmetry according to X-ray diffraction. Coupled to the observation of Néel skyrmions via Lorentz-TEM, our structural analysis provides a convincing explanation for their presence in centrosymmetric materials. Magnetization measurements as a function of the temperature displays a sharp first-order thermodynamic phase-transition leading to a reduction in the magnetic moment. This implies that the ground state of Fe 3 GaTe 2 is globally ferrimagnetic and not a glassy magnetic state composed of ferrimagnetic, and ferromagnetic domains as previously claimed. Neutron diffraction studies indicate that the ferromagnetic to ferrimagnetic transition upon reducing the external magnetic field is associated with a change in the magnetic configuration/coupling between Fe1 and Fe2 moments. We observe a clear correlation between the hysteresis observed in both the skyrmion density and the magnetization of Fe 3 GaTe 2 . This indicates that its topological spin textures are affected by the development of ferrimagnetism upon cooling. Observation, via magnetic force microscopy, of magnetic bubbles at the magnetic phase boundary suggests skyrmions stabilized by the competition among magnetic phases and distinct exchange interactions. Our study provides an explanation for the observation of Néel skyrmions in centrosymmetric systems, while exposing a correlation between the distinct magnetic phases of Fe 3 GaTe 2 and topological spin textures.

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

Electron Inversion and Tunneling at Silicon Thermal Oxide Interfaces for Solar-Driven Molecular Catalysis to Syngas

Semiconductor photoelectrodes are regularly coupled to solid-state heterogeneous catalysts to perform solar-driven reduction of CO 2 . Less frequently, molecular catalysts are employed to better control the reactivity toward desired products, yet the development of robust semiconductor/molecule interfaces has proven challenging. Here, we demonstrate that a 2–3 nm thermal oxide layer on Si exhibits stability in aqueous solution, high photovoltage, and a photocurrent density of ∼10 mA/cm 2 for the solar-driven photoelectrochemical reduction of a homogeneous molecular catalyst, producing syngas with an ∼2:1 H 2 to CO ratio. Because of a low defect density, the oxide interface forms an electron inversion layer with metal-like electron density at cathodic potentials. This inversion layer facilitates electron transfer to redox-active molecules via tunneling even if the molecule’s reduction potential is beyond the semiconductor’s conduction band edge. Using an electrolyte solution composed of a homogeneous cobalt bis(terpyridine) catalyst in a water/organic solvent mixture, stable photoelectrochemistry was observed under 1-sun illumination, exhibiting an ∼30% Faradaic efficiency for CO that was similar to a glassy carbon electrode under comparable conditions. Furthermore, the results demonstrate that an ultrathin thermal oxide interface is a robust platform for development of aqueous-stable, molecule-driven photoelectrocatalysis.

Catalysts↗
DOE OSTI2025

Catalytic Activity of an Ensemble of Sites for CO 2 Hydrogenation to Methanol on a ZrO 2 -on-Cu Inverse Catalyst

The significant increase in CO 2 emissions from heavy fossil fuel utilization has raised serious concerns, highlighting the need for effective methods to convert CO 2 into value-added chemicals. Here, in this work, we report a computational investigation on the catalytic activity of ZrO 2 -on-Cu inverse catalysts for CO 2 hydrogenation to methanol, considering highly dispersed ZrO 2 trimers on Cu (111). Such clusters present a large ensemble of formate-containing configurations, Zr 3 O n (OH) m (OCHO) l , making the evaluation of the catalytic activity very challenging. We found that the sites on the various catalyst configurations exhibit markedly different activities for formate hydrogenation, despite their similar free energy and composition. To understand these differences in reactivity, we examined the structural and electronic nature of the low free-energy catalyst configurations and identified that the energy of the lowest unoccupied orbital of the reacting formate, modified by its binding with the catalytic site, is a descriptor for the reaction energy of the formate hydrogenation step. From there, we screened an ensemble of catalyst structures using this descriptor to predict highly active metastable catalyst configurations and computed the reaction pathways and transition states for formate hydrogenation. From this investigation, we distinguished reactive from nonreactive sites and formate species on the ZrO 2 /Cu inverse catalyst based on structural and electronic features. We showed that rare metastable configurations control the activity. Additionally, an efficient method for examining the reactivity of a large number of coexisting catalyst structures was developed.

catalysts↗
DOE OSTI2025

Diffusion Model-Guided Inverse Design of Bimetallic Catalysts for Ammonia Decomposition

In the past decade, artificial intelligence and deep learning have played increasingly prominent roles in materials design and discovery. Among these, generative AI models, known for their ability to create unique and complex structures, have emerged as state-of-the-art tools for materials screening due to their high efficiency and low computational cost. In catalysis, one of the major challenges is identifying promising material candidates within an immense chemical space. This challenge can be addressed using generative approaches, such as diffusion-based inverse design models. In this study, we present a machine learning-guided workflow that employed a diffusion model for the inverse design of bimetallic alloy catalysts for low-carbon ammonia decomposition, a key reaction for ammonia emission control and sustainable hydrogen production. Catalyst candidates were evaluated using nitrogen adsorption energy as the key descriptor, inspired by multiscale modeling. The proposed workflow identified low-cost, environmentally friendly catalysts with excellent catalytic performance, which have been validated theoretically and experimentally. Our framework decoupled the generative and property-prediction components, enhancing both flexibility and accuracy in the catalytic material design process.

Adsorption↗
DOE OSTI2025

Structure of the Ecuadorian Upper Plate From a Joint Seismic‐Gravity Inversion

The Ecuadorian portion of the South American subduction zone presents an interesting case study in the structure and complex evolution of an upper plate. There are outstanding questions about its tectonic history, composition, and magmatic processes. While previous studies have employed ambient noise tomography to image the Ecuadorian upper plate, surface wave inversions alone often lack sensitivity at relevant shallow depths. This limitation can be overcome with an independent, complementary data set, such as gravity. We have jointly inverted Rayleigh wave phase velocities and Bouguer gravity anomalies to provide a more detailed seismic velocity model of the Ecuadorian upper plate. Our joint inversion has yielded several key improvements from previous models. First, we observe much shallower slow velocities beneath major basins (the Manabí, Progreso, and Gulf of Guayaquil), better aligning with expected basin structure. Second, we identify a high-velocity block beneath the entire forearc, corresponding to the Piñon Terrane, with velocities suggesting the presence of ultramafic material. Third, we highlight a new narrow swath of slow velocities beneath the Ecuadorian Andes, which closely follows the active volcanoes along the Eastern Cordillera. The extent of these slow velocities coincides with the termination of active arc volcanism and the predicted location of the subducted Carnegie Ridge. The predicted compositions for the mid to lower crust in the region preclude a purely compositional explanation for these velocities, suggesting that some level of partial melt is necessary.

Birkey, Andrew [Univ. of Delaware, Newark, DE (Uni↗
DOE OSTI2025

Customizable wave tailoring nonlinear materials enabled by bilevel inverse design

Abstract Passive wave transformation via nonlinearity is ubiquitous in settings from acoustics to optics and electromagnetics. It is well known that different nonlinearities yield different effects on propagating signals, which raises the question of “what precise nonlinearity is the best for a given wave tailoring application?” In this work, considering a one-dimensional spring-mass chain connected by polynomial springs (a variant of the Fermi-Pasta-Ulam-Tsingou system), we introduce a bilevel inverse design method which couples the shape optimization of structures for tailored constitutive responses with reduced-order nonlinear dynamical inverse design. We apply it to two qualitatively distinct problems—minimization of peak transmitted kinetic energy from impact, and pulse shape transformation—demonstrating our method’s breadth of applicability. For the impact problem, we obtain two fundamental insights. First, small differences in nonlinearity can drastically change the dynamic response of the system, from severely under- to outperforming a comparative linear system. Second, the oft-used strategy of impact mitigation via “energy locking” bistability can be significantly outperformed by our optimal nonlinearity. We validate this case with impact experiments and find excellent agreement. This study establishes a framework for broader passive nonlinear mechanical wave tailoring material design, with applications to computing, signal processing, shock mitigation, and autonomous materials.

Science & Technology - Other Topics↗
DOE OSTI2026

Stable and tunable MeV $$\gamma$$-ray generation via dual-laser inverse Thomson scattering from a laser-plasma accelerator

Abstract Inverse Thomson scattering from laser-plasma accelerators offers a pathway to compact, tunable MeV $$\gamma$$ -ray sources for reduced-dose radiography and enhanced performance in nuclear resonance fluorescence (NRF)-based isotope identification. However, photon yield and spectral quality are often limited by constraints on interaction geometry and scatter-laser tunability. Here we demonstrate a MeV $$\gamma$$ -ray source based on a dual-laser inverse Thomson scattering configuration driven by a 100-TW laser-plasma accelerator. Electron beams tunable from 122 to 204 MeV with $$<5$$ mrad divergence and $$<1$$ mrad pointing stability generate $$\gamma$$ rays with peak energies from 276 keV to 1.2 MeV and yields up to $$2\times 10^{7}$$ photons per shot. By independently controlling the interaction position and the scatter-pulse duration, we experimentally match the scatter pulse to the walk-off-limited interaction length. Extending the scatter pulse to 200 fs increases photon production by approximately $$15\%$$ while maintaining operation in the linear Thomson regime, thereby preserving narrow spectral bandwidth and controlled radiation divergence. Radiographic characterization demonstrates MeV-level penetration and $$\approx 0.1$$ mm spatial resolution, while stable operation is sustained over multi-hour timescales across multiple days. These results show that interaction-length optimization provides a scalable strategy for improving photon yield, spectral control, and operational stability in compact laser-plasma-accelerator-driven $$\gamma$$ -ray sources.

Tsai, Hai-En↗
DOE OSTI2025

Dynamics of inverse metal oxides on metal catalysts using spectro-kinetics: reversible Brønsted acid site formation and irreversible reduction

Brønsted acid sites (BASs) in inverse catalysts are vital for the selective hydrogenolysis of polyols, specifically cleaving secondary C–O bonds. These BASs form dynamically in situ in an H 2 environment. While H 2 enables rapid BAS generation on short timescales, it reduces the catalyst at prolonged exposures. The active center for BAS generation, the kinetics of BAS formation, its reverse decomposition, and the irreversible oxide reduction have lacked direct experimental evidence. Here, aided by advanced spectro-kinetic studies, we identify trimeric W 3 O x sites on Pt as the active centers for BAS generation, whereas isolated WO x species on SiO 2 act merely as spectator species, demonstrated using an inverse WO x /Pt catalyst as a representative system. A detailed kinetic profile capturing the dynamics of W 3 O x sites on Pt is also established. The rate constant for BAS formation is two orders of magnitude higher than for its decomposition, which is one order of magnitude faster than the irreversible site reduction. Co-fed H 2 O suppresses the site reduction by ∼50%. Furthermore, the H 2 partial pressure plays an important role. While lower gas-phase H 2 partial pressure does not influence the reversible BAS formation, it can significantly (∼3×) suppress catalyst reduction. Finally, these findings offer critical insights into optimizing reaction conditions through periodic H 2 pulsing, enhancing catalyst stability and performance in hydrogenolysis reactions.

Sourav, Sagar [Indian Inst. of Technology (IIT), M↗
DOE OSTI2025

The influence of protein electrostatics on potential inversion in flavoproteins

Biology uses relatively few electron-transfer cofactors, tuning their potentials, electronic couplings, and reorganization energies to carry out the required chemistry. It is remarkable that the potential ordering of two-electron transfer active flavins can be normal (first oxidation at low potential and second oxidation at high potential) or inverted, and the gap between the potentials can be as large as one volt. Analysis based on structural bioinformatics and electrostatics indicates that the ordering of the flavin redox potential is influenced by protein electrostatics. In all 36 flavoproteins examined, the introduction of a negative charge near the flavin in silico increases the extent of potential inversion (by lowering the electrochemical potential of the second electron-transfer step); the introduction of a positive charge near the flavin favors normally ordered potentials. We also find that the addition of positive charges increases the electrochemical potential for the naturally occurring one-electron transition in flavodoxins (between deprotonated hydroquinone and neutral semiquinone) and also increases the second one-electron transition in bifurcating flavins (between anionic semiquinone and fully oxidized flavin). Finally, we find that proximity of a proton acceptor, notably conserved arginine, supports proton-coupled electron transfer because it may act as a proton acceptor, promoting potential inversion. This key arginine residue may enable two-electron transfer chemistry by promoting the proton-coupled electron transfer process over the pure electron transfer process, suggesting how a protein's flavin environment may influence one- or two-electron chemistry in flavoproteins.

Singh, Niven [Duke Univ., Durham, NC (United State↗
DOE OSTI2026

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗
DOE OSTI2025

Scattering-based structural inversion of soft materials via Kolmogorov–Arnold networks

Small-angle scattering techniques are indispensable tools for probing the structure of soft materials. However, traditional analytical models often face limitations in structural inversion for complex systems, primarily due to the absence of closed-form expressions of scattering functions. To address these challenges, we present a machine learning framework based on the Kolmogorov–Arnold Network (KAN) for directly extracting real-space structural information from scattering spectra in reciprocal space. This model-independent, data-driven approach provides a versatile solution for analyzing intricate configurations in soft matter. By applying the KAN to lyotropic lamellar phases and colloidal suspensions—two representative soft matter systems—we demonstrate its ability to accurately and efficiently resolve structural collectivity and complexity. Here, our findings highlight the transformative potential of machine learning in enhancing the quantitative analysis of soft materials, paving the way for robust structural inversion across diverse systems.

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

Statistical modelling and Bayesian inversion for a Compton imaging system: application to radioactive source localization

Abstract This paper presents a statistical forward model for a Compton imaging system, called Compton imager. This system, under development at the University of Illinois Urbana Champaign, is a variant of Compton cameras with a single type of sensors which can simultaneously act as scatterers and absorbers. This imager is convenient for imaging situations requiring a wide field of view. The proposed statistical forward model is then used to solve the inverse problem of estimating the location and energy of point-like sources from observed data. This inverse problem is formulated and solved in a Bayesian framework by using a Metropolis within Gibbs algorithm for the estimation of the location, and an expectation-maximization algorithm for the estimation of the energy. This approach leads to more accurate estimation when compared with the deterministic standard back-projection approach, with the additional benefit of uncertainty quantification in the low photon imaging setting.

Tarpau, Cécilia (ORCID:0000000286539490)↗
DOE OSTI2024

One-shot omnidirectional pressure integration through matrix inversion

In this work, we present a method to perform 2D and 3D omnidirectional pressure integration from velocity measurements with a single-iteration matrix inversion approach. This work builds upon our previous work, where the rotating parallel ray approach was extended to the limit of infinite rays by taking continuous projection integrals of the ray paths and recasting the problem as an iterative matrix inversion problem. This iterative matrix equation is now 'fast-forwarded' to the 'infinity' iteration, leading to a different matrix equation that can be solved in a single step, thereby presenting the same computational complexity as the Poisson equation. We observe computational speedups of ~10 6 when compared to brute-force omnidirectional integration methods, enabling the treatment of grids of ~10 9 points and potentially even larger in a desktop setup at the time of publication. Further examination of the boundary conditions of our one-shot method shows that omnidirectional pressure integration implements a boundary condition where the boundary points are treated as interior points to the extent that information is available. Finally, we show how the method can be extended from the regular grids typical of particle image velocimetry to the unstructured meshes characteristic of particle tracking velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Divertor heat load estimates on NSTX and DIII-D using new and open-source 2D inversion analysis code

A thermography inversion algorithm has been developed in the open-source Python-based computer code, HYPERION, to calculate the heat flux incident on plasma-facing components (PFCs) in axisymmetric tokamaks. The chosen mesh size at the surface significantly affects the calculated transient heat flux results. The calculated transient heat flux will exceed the real value when the mesh size tends to zero but will underestimate the real value when the mesh size is large. A criterion for determining the appropriate mesh size for the transient heat flux calculation will be discussed. The numerical scheme for HYPERION uses a 2D fully implicit finite-difference approach, allowing temperature-dependent thermal properties of PFC materials. The inversion algorithm is benchmarked against established heat flux calculation codes, TACO and THEODOR, based on thermography data from NSTX and DIII-D respectively. The primary benefits of HYPERION compared to TACO and THEODOR are that it is open-source and it allows for the optimization of mesh thickness along the substrate. The algorithm also accounts for the thermal properties of thin surface layers that characteristically form on PFCs due to plasma-material interactions. The agreement between HYPERION and THEODOR is excellent, as the percent difference between the codes is ~5% on average in the case of the DIII-D data for moderate to high heat flux. Verification tests with TACO show slightly higher average percent differences of 8% and 12%. In using HYPERION to study filaments in heat flux, the initial results indicate that small ELMs filaments significantly broaden the divertor heat flux, and decrease divertor peak flux. Compared to the inter-ELM, the small ELM filaments decrease the divertor peak surface temperature. With intermittent divertor filaments, the divertor heat flux width is comparable with that found in L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗