Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHARP”

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

Preserving the Josephson Coupling of Twisted Cuprate Junctions via Tailored Silicon Nitride Circuits Boards

Controlled fabrication of twisted van der Waals heterostructures is essential to unlock the full potential of moiré materials. However, achieving reproducibility remains a major challenge, particularly for air-sensitive materials such as Bi 2 Sr 2 CaCu 2 O 8 + δ (BSCCO), where it is crucial to preserve the intrinsic and delicate superconducting properties of the interface throughout the entire fabrication process. Here, a dry, inert and cryogenic assembly method is presented that combines silicon nitride nanomembranes (NMBs) with pre-patterned electrodes and the cryogenic stacking technique (CST) to fabricate high-quality twisted BSCCO Josephson junctions (JJs). This protocol prevents thermal and chemical degradation during both interface formation and electrical contact integration. It is also found that asymmetric membrane designs, such as a double cantilever, effectively suppress vibration-induced disorder due to wire bonding, resulting in sharp and hysteretic current–voltage characteristics. The junctions exhibit a twist-angle-dependent Josephson coupling with magnitudes comparable to the highest-performing devices reported to date, but achieved through a straightforward and versatile contact method, offering a scalable and adaptable platform for future applications. These findings highlight the importance of both interface and contact engineering in addressing reproducibility in superconducting van der Waals heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Corners and islands in the S-matrix bootstrap of the open superstring

We bootstrap the Veneziano superstring amplitude in 10 dimensions from the bottom-up. Starting with the most general maximally supersymmetric Yang-Mills EFT, we input information about the lowest-lying massive states, which we assume contribute via tree-level exchanges to the 4-point amplitude. We show the following: (1) if there is only a single state at the lowest mass, it must be a scalar. (2) Assuming a string-inspired gap between the mass of this scalar and any other massive states, the allowed region of Wilson coefficients has a new sharp corner where the Veneziano amplitude is located. (3) Upon fixing the next massive state to be a vector, the EFT bounds have a one-parameter family of corners; these would correspond to models with linear Regge trajectories of varying slopes, one of which is the open superstring. (4) When the ratio between the massive scalar coupling and the tr F 4 coefficient is fixed to its string value, the spin and mass of the second massive state is determined by the bootstrap and the Veneziano amplitude is isolated on a small island in parameter space. Finally, we compare with other recent bootstraps approaches, both the pion model and imposing Regge-inspired maximal spin constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Splitting regions and shrinking islands from higher point constraints

We study constraints from higher-point amplitudes on 2 → 2 scattering in the context of effective field theory (EFT) using the perturbative numerical S-matrix bootstrap. Specifically, we investigate the class of weakly coupled EFTs with amplitudes that obey the hidden zero and split conditions that are known to hold both for Tr(Φ 3 ) theory and for certain string tree amplitudes, including at 4-point the beta function. Requiring the splitting condition for the 5-point amplitude not only fixes nearly all its contact terms, but it also imposes non-linear constraints among the 4-point EFT Wilson coefficients. When included in the bootstrap, the resulting allowed region consistent with positivity is no longer convex but is restricted to a smaller non-convex region — which has a sharp corner near the string beta function! Assuming the absence of an infinite spin tower at the mass gap, the allowed region bifurcates into a trivial region (with states only above a chosen cutoff) and an island that continues to shrink around the string as more constraints are included in the bootstrap. The numerics indicate that in the absence of single-mass infinite spin towers the string beta function is the unique 4-point amplitude compatible with hidden zero and the 5-point splitting constraints. The analysis provides a prototype example for how features of higher-point amplitudes constrain the bootstrap of 4-point amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Spectral Bounds on Hyperbolic 3-Manifolds: Associativity and the Trace Formula

We constrain the low-energy spectra of Laplace operators on closed hyperbolic manifolds and orbifolds in three dimensions, including the standard Laplace--Beltrami operator on functions and the Laplacian on powers of the cotangent bundle. Our approach employs linear programming techniques to derive rigorous bounds by leveraging two types of spectral identities. The first type, inspired by the conformal bootstrap, arises from the consistency of the spectral decomposition of the product of Laplace eigensections, and involves the Laplacian spectra as well as integrals of triple products of eigensections. We formulate these conditions in the language of representation theory of PSL 2 (C) and use them to prove upper bounds on the first and second Laplacian eigenvalues. The second type of spectral identities follows from the Selberg trace formula. We use them to find upper bounds on the spectral gap of the Laplace--Beltrami operator on hyperbolic 3-orbifolds, as well as on the systole length of hyperbolic 3-manifolds, as a function of the volume. Further, we prove that the spectral gap λ 1 of the Laplace--Beltrami operator on all closed hyperbolic 3-manifolds satisfies λ 1 < 47.32. Along the way, we use the trace formula to estimate the low-energy spectra of a large set of example orbifolds and compare them with our general bounds, finding that the bounds are nearly sharp in several cases.

Bonifacio, James [University of Mississippi, MS (U↗
DOE OSTI2026

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
DOE OSTI2024

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems—like fixation of cells coupled with fluorescence microscopy—provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. In this study, these snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Laser-induced selective local patterning of vanadium oxide phases

The same elements can form different compounds with widely different physical properties. Synthesis of a single-phase material is commonly achieved by controlling experimental conditions. Synthesizing materials that incorporate multiple specific spatially distributed chemical phases is often challenging, especially if different phases must be organized into well-defined spatial patterns. Here, we present an efficient solid reaction laser annealing (SRLA) approach to directly write regions of different local chemical compositions. We demonstrate the practical utility of our approach by locally writing microscale patterns of distinct chemical phases in vanadium oxide thin films. Specifically, we achieved the controlled local recrystallization of a uniform V 2 O 3 matrix into VO 2 , V 3 O 5 , and V 4 O 7 regions exhibiting sharp 1st- and 2nd-order metal–insulator phase transitions over a wide range of critical temperatures, i.e., a characteristic feature of select vanadium oxides that is extremely sensitive to even minute structural or compositional imperfections. We utilized the local chemical phase writing to pattern spiking oscillators with distinct electrical behavior directly in the thin film sample without employing elaborate lithography fabrication. Our laser tuning local chemical composition opens a pathway to synthesize a wide range of artificially micropatterned composite materials, with precision and control unattainable in conventional material synthesis methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Causality-respecting adaptive refinement for PINNs: enabling precise interface evolution in phase field modeling

Physics-informed neural networks (PINNs) have emerged as a powerful tool for solving physical systems described by partial differential equations (PDEs). However, their accuracy in dynamical systems, particularly those involving sharp moving boundaries with complex initial morphologies, remains a challenge. Here, this study introduces an approach combining residual-based adaptive refinement (RBAR) with causality-informed training to enhance the performance of PINNs in solving spatio-temporal PDEs. Our method employs a three-step iterative process: initial causality-based training, RBAR-guided domain refinement, and subsequent causality training on the refined mesh. Applied to the Allen-Cahn equation, a widely-used model in phase field simulations, our approach demonstrates significant improvements in solution accuracy and computational efficiency over traditional PINNs. Notably, we observe an ‘overshoot and relocate’ phenomenon in dynamic cases with complex morphologies, showcasing the method’s adaptive error correction capabilities. This synergistic interaction between RBAR and causality training enables accurate capture of interface evolution, even in challenging scenarios where traditional PINNs fail. Our framework not only resolves the limitations of uniform refinement strategies but also provides a generalizable methodology for solving a broad range of spatio-temporal PDEs. The enhanced performance of the RBAR–causality combined framework demonstrates its strong potential for advancing PINN-based modeling of physical systems characterized by complex, evolving interfaces.

Allen-Cahn equations↗
DOE OSTI2026

E-PINNs: Epistemic Physics-Informed Neural Networks

Physics-informed neural networks (PINNs) have demonstrated promise as a framework for solving forward and inverse problems involving partial differential equations. Despite recent progress in the field, it remains challenging to quantify uncertainty in these networks. While techniques such as Bayesian PINNs (B-PINNs) provide a principled approach to capturing epistemic uncertainty through Bayesian inference, they can be computationally expensive for large-scale applications. In this work, we propose Epistemic Physics-Informed Neural Networks (E-PINNs), a framework that uses a small network, the epinet, to efficiently quantify epistemic uncertainty in PINNs. The proposed approach works as an add-on to existing, pre-trained PINNs with a small computational overhead. We demonstrate the applicability of the proposed framework in various test cases and compare the results with B-PINNs using Hamiltonian Monte Carlo (HMC) posterior estimation and dropout-equipped PINNs (Dropout-PINNs). In our experiments, E-PINNs achieve calibrated coverage with competitive sharpness at substantially lower cost. We demonstrate that when B-PINNs produce narrower bands, they under-cover in our tests. E-PINNs also show better calibration than Dropout-PINNs in these examples, indicating a favorable accuracy-efficiency trade-off.

AI for Science↗
DOE OSTI2026

Quantitative characterization of gradient microstructures: A study on friction stir spot processing of pure cobalt

Heterogeneous microstructures in polycrystalline metals can enhance the strength and ductility, outperforming homogeneous structures of similar composition. This study investigates deformed cobalt via friction stir spot processing with varying dwell times to uncover the effects of plastic deformation and heat generation on the formation of morphological, phase, and grain boundary character gradients. A new approach to quantify the morphological gradients in materials, which describes grain morphology in terms of density followed by parametric regression, enables direct quantification of processing depth and gradient sharpness. Results show that longer processing times increase the steepness of morphological gradients and reduce the deformation depth for friction stir spot processing with low plunge depths and high tool rotational speeds. The amount of retained FCC is increased in the shorter processing conditions, primarily due to refined grain size, increased defect content, and reduced heat generation. Crystallographic texture analysis of the HCP phase indicated a dominant B-fiber described by (0001) ∥ shear plane normal in the extreme processing conditions and the formation of a P-fiber, shear direction ∥ ⟨11$\bar2$0⟩ for intermediate dwell times. The texture of the FCC phase for low processing times was a C texture {100}⟨011⟩ where longer processing times were dominated by a [001] fiber texture with a main {110}⟨100⟩ orientation and emergence of a slight [111] fiber in the longest processing condition. The approaches outlined in this work give insight into quantifying gradients and improve the understanding of highly deformed cobalt.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High compressive energy absorption and shape recovery behavior of additively manufactured textile-inspired cylindrical braided metamaterials

Mechanical metamaterials (MMs) are engineered structures with unique mechanical properties that arise from their unique spatial arrangement or lattice-like structure. The most commonly designed MMs such as honeycomb and re-entrant auxetics are prone to failure at the sharp corners and weak joints due to the increased stress concentration under deformation. To mitigate this challenge, braided MM structures involving intertwining threads of nylon—forming curved unit cells—have been studied. These textile-inspired cylindrical braided metamaterials (CBMMs) with contrasting unit cells, namely diamond and regular CBMMs, were fabricated by 3D printing. The layer-by-layer deposited structure built by fused filament fabrication delivered an assembly of overlapped threads that are fused at the contact point. To understand deformation behavior of these MMs, finite element models were developed for various load scenarios including quasi-static compression, cyclic and creep loads at room temperature. Stress distribution, deformation mechanisms, and failure modes were analyzed and validated by experiments to analyze the geometries and associated performance. The diamond CBMMs showed stress softening at 30 % compressive strain, withstanding a load of ∼440 N, whereas the regular CBMMs at 50 % strain experienced ∼250 N. The diamond CBMMs delivered higher creep resistance under sustained load and better energy absorption under cyclic loading than the regular CBMMs. The latter, however, exhibited 94 % shape recovery in contrast to 88 % recovery in former prototype during their first cyclic load. In conclusion, this study helps design mechanical lightweight devices that endure significant sustained load and exhibit enhanced energy absorption and shape recovery characteristics in cyclic loading.

Creep↗
DOE OSTI2024

A Modeling framework for flocculated cohesive sediment transport in the current bottom boundary layer

Cohesive sediment transport, where its settling velocity is controlled by the flocculation process, is a crucial component in determining biochemical cycles, fate of pollutants, and morphodynamics in many aquatic ecosystems. In this study, a modeling framework is presented to investigate how flocculation influences cohesive sediment transport in the current bottom boundary layer in dilute conditions, consistent with the calibration range of the flocculation model. From a local analysis of floc dynamics in homogenous turbulence, we identify that the floc size distribution is mainly controlled by floc cohesion and yield strength. The uncertainty in fractal dimension plays a minor role for the floc size but it influences the resulting floc density and settling velocity. The transport analysis in the current boundary layer shows that the flocculation process alters the vertical distribution of the settling velocity and hence the sediment concentration with a strong dependence on cohesion, floc yield strength, and floc structure. When the flocs are more susceptible to breaking, a well-mixed concentration profile is obtained. In contrast, for flocs with higher cohesion or yield strength, higher concentration with a sharp gradient is observed close to the bed. Overall, the settling velocity exhibits a low vertical variability within 20 % of the depth-averaged value except near the bed. Further, this suggests that using a depth-averaged settling velocity yields acceptable predictions of the sediment concentration profiles, especially for flocs with lower cohesion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Unveiling the drivers contributing to global wheat yield shocks through quantile regression

Sudden reductions in crop yield (i.e., yield shocks) severely disrupt the food supply, intensify food insecurity, depress farmers' welfare, and worsen a country's economic conditions. Here, we study the spatiotemporal patterns of wheat yield shocks, quantified by the lower quantiles of yield fluctuations, in 86 countries over 30 years. Furthermore, we assess the relationships between shocks and their key ecological and socioeconomic drivers using quantile regression based on statistical (linear quantile mixed model) and machine learning (quantile random forest) models. Using a panel dataset that captures spatiotemporal patterns of yield shocks and possible drivers in 86 countries, we find that the severity of yield shocks has been increasing globally since 1997. Moreover, our cross-validation exercise shows that quantile random forest outperforms the linear quantile regression model. Despite this performance difference, both models consistently reveal that the severity of shocks is associated with higher weather stress, nitrogen fertilizer application rate, and gross domestic product (GDP) per capita (a typical indicator for economic and technological advancement in a country). While the unexpected negative association between more severe wheat yield shocks and higher fertilizer application rate and GDP per capita does not imply a direct causal effect, they indicate that the advancement in wheat production has been primarily on achieving higher yields and less on lowering the possibility and magnitude of sharp yield reductions. Hence, in the context of growing extreme weather stress, there is a critical need to enhance the technology and management practices that mitigate yield shocks to improve the resilience of the world food systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

High-current-density electrosynthesis of formate from captured CO 2 solution by MOF-derived bismuth nanosheets

Greenhouse gas emissions present a significant challenge to humanity, and utilizing renewable electricity to convert emitted CO 2 into value-added products offers a promising solution; however, traditional CO 2 capture and regeneration processes remain energy-intensive, restricting the overall system efficiency and decarbonization efficacy. In this study, an advanced direct reduction of captured CO 2 with large current densities for formate electrosynthesis was demonstrated without the need for CO 2 regeneration or compression. The bismuth nanosheet (DRM-BiNS) was synthesized by direct reduction of a Bi-based MOF, representing a new class of catalytic materials with a large surface area and interconnected pores, suitable for the direct reduction of captured CO 2 . By seamlessly combining experimentation and simulation, insights into the structure-parameter-performance relation were acquired in a flow cell setting, including critical membrane-electrode distance, cell orientation, and pumping flow rate. Important flow-cell components, such as catholyte volume, electrode substrate, membrane choice, and ionomer type, were also carefully examined to enhance the cell performance. In sharp contrast to prior studies limited to current densities below 20 mA/cm 2 in bicarbonate-based captured CO 2 solutions, this work demonstrates a remarkable current density of 300 mA/cm 2 with an FE to formate comparable to the case with gas-fed CO 2 reduction. Moreover, the process sustained an FE above 50% at a high current density of 500 mA/cm 2 . The DRM-BiNS catalyst exhibited outstanding selectivity, activity, and stability, significantly outperforming oxide-derived bismuth nanosheets (OD-BiNS) in captured CO 2 reduction. Furthermore, these findings offer critical insights into the development of sustainable and scalable CO 2 utilization technologies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Simulated moving bed-inspired method for continuous adsorptive denitrogenation of model fuel

Efficient, continuous routes for removing nitrogen-containing compounds from hydrothermal liquefaction-derived synthetic aviation fuel are needed to enable direct blending with conventional jet fuels. Here, we report a simulated moving bed-inspired process for adsorptive denitrogenation of a model fuel. Unlike conventional simulated moving bed systems, which are designed for sharp separations between similar solutes, this approach was run deliberately outside the classical separation region so that both pyridine and indole were removed together from the hydrocarbon stream. Alcohol solvents were used to regenerate the silica adsorbent, maintaining performance over extended operation and avoiding the downtime and energy demand associated with calcination. Under these conditions, the system demonstrates removal of more than 98% of nitrogen while cutting solvent use by 28% compared to batch operation. Classical modeling tools predicted column concentration profiles even in this nontraditional regime, suggesting a straightforward path to scaling. Together, these results motivate solvent-efficient, continuous denitrogenation strategies that could be integrated with biorefinery processes.

Adsorption↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗