Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEAR CHEMISTRY”

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.

1,702 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗
DOE OSTI2026

Liquified SO 2 induced solid/cathode electrolyte interphase for lithium ion batteries

Formation of robust solid/cathode electrolyte interphases (S/CEI) is vital for long-term stability and high-performance operation of lithium-ion batteries (LIBs), particularly under high voltage regimes. However, engineering electrochemically stable S/CEIs that effectively suppress interfacial side reactions remains a key challenge. Herein, we introduce a liquefied sulfur dioxide (SO 2 )– ionic liquid complex as a fluorine-free multifunctional electrolyte additive for the first time that significantly improves the formation of sulfate/sulfite-rich S/CEI layers at both graphite and NMC811 interfaces. The unique SO 2 -N coordination with a 1,2,4-triazolide-based ionic liquid enables homogeneous SO 2 dissolution, resulting in controlled SO 2 decomposition during the initial electrochemical cycle. This decomposition yields sulfur-rich interphase species that stabilize the electrolyte-electrode interface, reduce impedance growth, and lessen electrolyte decomposition. Electrochemical tests show significantly improved cycle life, reduced polarization, and increased Coulombic efficiency for both anodes and cathodes. XPS confirms the presence of SO 2 -derived surface species that contribute to interfacial stability. In conclusion, this approach highlights a new direction for interphase engineering using liquefied gas additives and opens pathways for sulfur-based S/CEI chemistry in advanced battery systems.

Graphite↗
DOE OSTI2026

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗
DOE OSTI2024

Charmonium Transport in Heavy-Ion Collisions at the LHC

We provide an update on our semi-classical transport approach for quarkonium production in high-energy heavy-ion collisions, focusing on 𝐽/𝜓 and 𝜓(2𝑆) mesons in 5.02 TeV Pb-Pb collisions at the Large Hadron Collider (LHC) at both forward and mid-rapidity. In particular, we employ the most recent charm-production cross sections reported in pp collisions, which are pivotal for the magnitude of the regeneration contribution, and their modifications due to cold-nuclear-matter (CNM) effects. Multi-differential observables are calculated in terms of nuclear modification factors as a function of centrality, transverse momentum, and rapidity, including the contributions from feeddown from bottom hadron decays. For our predictions for 𝜓(2𝑆) production, the mechanism of sequential regeneration relative to the more strongly bound 𝐽/𝜓 meson plays an important role in interpreting recent ALICE data.

Astronomy & Astrophysics↗
DOE OSTI2026

NOx mediated formation of ε – UO 3 during thermal decomposition of uranyl nitrate hexahydrate under static air conditions

The thermal decomposition of uranyl nitrate hexahydrate in air at temperatures of 650 – 800 °C is expected to yield α – U 3 O 8 through intermediate UO 3 phases. Here, n this study, Raman spectroscopy complimented by powder X-ray diffraction revealed the unexpected formation of ε – UO 3 as an accompanying phase during thermal decomposition of UNH at 700 °C under static air conditions. Experiments demonstrated that ε – UO 3 does not form during the initial ramp up stage nor during the high temperature heating period; but instead forms during the cooling stage between 250 – 400 °C. Additional experiments revealed that UNH initially decomposes through an amorphous UO 3 intermediate prior to α – U 3 O 8 formation. Attempts to bypass the U 3 O 8 precursor in the formation of ε – UO 3 were unsuccessful, supporting the necessity of U 3 O 8 (α – U 3 O 8 in this study) in the formation of ε – UO 3 . The observed phase evolution is proposed to result from the NOx species generated during UNH thermal decomposition, which create localized oxidizing conditions within the furnace under static air conditions. Furthermore, a predominantly phase pure ε – UO 3 was synthesized under static air conditions through the addition of an extra plateau at 250 °C during the cooling step. These findings demonstrate the importance of gas-phase chemistry and cooling conditions in uranium oxide phase evolution and provide additional insight into ε – UO 3 formation pathways.

Epsilon uranium trioxide↗
DOE OSTI2026

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗
DOE OSTI2026

Reconstruction and Dissolution of Copper Catalysts during Electrocatalytic Nitrate Reduction

Copper-based electrocatalysts are widely explored for electrochemical nitrate remediation, yet their stability under operating conditions remains poorly understood. While nanoscale and subnanoscale Cu motifs are known to restructure during the nitrate reduction reaction (NO 3 RR), how these transformations translate into irreversible material loss remains unclear. Here, we quantify Cu dissolution during NO 3 RR as a function of catalyst architecture and electrolyte chemistry using two model systems: single-atom Cu (Cu 1 ) and Cu nanoparticles (Cu NP ). Time-resolved leaching measurements, integrated with in situ X-ray absorption spectroscopy (XAS), reveal measurable Cu loss for both catalysts during NO 3 RR. Dissolution is concentrated at the initiation of electrolysis, coinciding with rapid restructuring, and depends strongly on morphology, with Cu NP consistently exhibiting greater Cu loss than Cu 1 . In situ XAS reveals that Cu 1 forms transient metallic clusters under reduction that largely redisperse upon returning to open-circuit voltage, whereas Cu NP undergoes reduction of an oxidized surface layer accompanied by sustained Cu loss during electrolysis. Notably, the presence of nitrate significantly intensifies restructuring and Cu loss, highlighting the critical role of electrolyte composition. Furthermore, these findings establish a direct link between electrochemical restructuring and Cu dissolution, unveiling electrolyte-dependent interactions as key determinants of catalyst durability in electrochemical nitrate conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Patchy nanoparticles by atomic stencilling

Stencilling, in which patterns are created by painting over masks, has ubiquitous applications in art, architecture and manufacturing. Modern, top-down microfabrication methods have succeeded in reducing mask sizes to under 10 nm, enabling ever smaller microdevices as today’s fastest computer chips. Meanwhile, bottom-up masking using chemical bonds or physical interactions has remained largely unexplored, despite its advantages of low cost, solution-processability, scalability and high compatibility with complex, curved and three-dimensional (3D) surfaces. Here we report atomic stencilling to make patchy nanoparticles (NPs), using surface-adsorbed iodide submonolayers to create the mask and ligand-mediated grafted polymers onto unmasked regions as ‘paint’. We use this approach to synthesize more than 20 different types of NP coated with polymer patches in high yield. Polymer scaling theory and molecular dynamics (MD) simulation show that stencilling, along with the interplay of enthalpic and entropic effects of polymers, generates patchy particle morphologies not reported previously. These polymer-patched NPs self-assemble into extended crystals owing to highly uniform patches, including different non-closely packed superlattices. We propose that atomic stencilling opens new avenues in patterning NPs and other substrates at the nanometre length scale, leading to precise control of their chemistry, reactivity and interactions for a wide range of applications, such as targeted delivery, catalysis, microelectronics, integrated metamaterials and tissue engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2026

High-Throughput Uniformity and Defect Monitoring in Low-Temperature Electrolysis Porous Transport Layers Using X-Ray Radiography

Effective quality control (QC) for manufacturing proton exchange membrane water electrolysis (PEMWE) components is critical to enabling widespread adoption of the technology for hydrogen generation. This study investigates X-ray radiography as a novel, high-throughput, potentially in-line QC technique for detecting defects and assessing material property distributions in titanium-based porous transport layers (PTLs) which constitute a crucial component of low temperature PEMWE stacks. We obtain radiographs of a set of fifteen PTLs and model their absorbance of the broadband radiation as a second-order polynomial to account for the non-monoenergetic radiation source used in this study. The resulting model serves as a basis for predicting the areal density and porosity distributions of the PTLs. We find radiography successful in detecting multiple instances of defects, including holes/depressions, cracks, and excess material on the surface or in the pores of the material, demonstrating its potential as a robust in-line QC tool for PTL manufacturing.

08 HYDROGEN↗
DOE OSTI2025

Mesoscale Magnetostructural Phase Separation in Fe‐deficient Fe 5 GeTe 2

Two-dimensional van der Waals ferromagnet Fe 5-x GeTe 2 (F5GT) is promising for spintronic applications due to its high Curie temperature, layered structure, and ability to host complex magnetic textures. However, the origin of its sample-dependent magnetic anisotropy remains unclear, hindering control of its magnetic behavior. Here, we use spatially resolved cryogenic scanning transmission electron microscopy (STEM) to correlatively map magnetism, lattice structure, and chemistry across atomic-to-micron scales. We reveal that only mesoscale, not nanoscale, inclusions of a Fe-deficient secondary phase significantly modify magnetic behavior, establishing a previously unrecognized critical length scale. This phase separation, induced by quenching, leads to in-plane magnetic anisotropy, while slow cooling confines separation to a few nanometers and preserves out-of-plane anisotropy. These findings reconcile prior inconsistencies and establish a predictive framework for tuning magnetism in F5GT through thermal processing, with broader implications for controlling anisotropy in other two-dimensional magnetic materials.

2D ferromagnets↗
DOE OSTI2026

Optical Vibrational Spectroscopic Investigation of Natural and Synthetic Analogs of Uranyl Oxyhydroxyhydrate Minerals

Uranyl oxy-hydroxy-hydrate minerals are common alteration products of uraninite (UO2+x) which is chemically and structurally analogous to uranium dioxide nuclear fuel. Therefore, structural and spectroscopic investigations of these alteration minerals and their analogous anthropogenic counterparts can provide insight into the environmental behavior of nuclear fuel cycle materials. Previously, we compiled available vibrational spectroscopic data for uranyl minerals in the Compendium of Uranium Raman and Infrared Experimental Spectra (CURIES) and found that only 37% of known uranyl oxy-hydroxy-hydrate minerals had spectra readily available in the literature and existing databases for inclusion therein. Furthermore, no available infrared spectra for this mineral group were included in CURIES. To expand our understanding of the spectroscopic features of uranyl oxy-hydroxy-hydrates and the structural origins thereof, we collected, and now include in CURIES, Raman and infrared spectra for an additional 12 uranyl hydroxide phases. To better understand the impact of structural and compositional variations of these phases on their spectroscopic features, we compare Raman spectra of different anion sheet topological groups and of analogous phases hosting different counter cations. We identify spectroscopic variations related to differences in equatorial bonding and structural changes as a result of cation substitution. We also prepare a uranyl hydroxide phase via hydrolysis of uranyl fluoride (UO2F2) as an analog of hydrolysis reactions that occur in nuclear fuel cycle materials; and we find that the alteration product of UO2F2, despite chemical and structural similarities to uranyl oxy-hydroxy-hydrate minerals, is readily distinguishable from related mineral phases using Raman spectroscopy. In this work, we provide new insights into the structural origins of spectroscopic features in uranyl oxy-hydroxy-hydrate minerals, improve the average Raman spectrum for this group of minerals, and thereby improve capabilities for identifying these mineral species and related anthropogenic phases using Raman spectroscopy.

Barth, Brodie [ORNL] (ORCID:0000000256142601)↗
DOE OSTI2026

The Future Polarized Target Program at Jefferson Lab

Polarized targets have played a crucial role in Jefferson Lab's exploration of nuclear structure over the past four decades. The three original experimental halls have seen 19 separate installations of polarized solid or gas targets for use in the particle physics scattering experiments, and this trend will continue in the next decade. Five polarized target systems are in preparation for use at JLab in the coming years. Hall B will see the use of two solid polarized targets, one longitudinally polarized to the beam, the other transversely, as well as a novel 3He gas polarized target. In Hall C, new experiments will augment the tensor polarization in dynamically polarized solids. Plans are under development to bring a polarized solid target to Jefferson Lab's photon beam hall, Hall D, for the first time. To support these efforts, the JLab polarized target group is building a test laboratory to develop dynamic nuclear polarization techniques, as well as an apparatus to irradiate target material using electrons from JLab's injector test facility. In this talk, we will explore the development progress and plans for each of these efforts.

Maxwell, James [Thomas Jefferson National Accelera↗
DOE OSTI2026

EvoDiffMol: evolutionary diffusion framework for 3D molecular design with optimized properties

Designing molecules with specific target properties remains a fundamental challenge in computational chemistry. While existing approaches show promise, most rely on simplified representations like SMILES strings or 2D graphs that lack essential three-dimensional geometric information. We present EvoDiffMol, a computational framework that integrates evolutionary algorithms with three-dimensional diffusion models for property-driven molecular generation. The method operates through adaptive evolutionary optimization, where population-based selection guides the generation process toward desired property landscapes. EvoDiffMol supports both unconstrained molecular design and scaffold-constrained generation that preserves fixed substructures while optimizing complementary regions. Comprehensive evaluation demonstrates exceptional performance, achieving the highest drug-likeness score (0.94) among all compared state-of-the-art methods while maintaining excellent validity, uniqueness, and novelty. Beyond single property optimization, the framework demonstrates flexible multi-property optimization capabilities, simultaneously controlling multiple molecular descriptors including synthetic accessibility, lipophilicity, topological polar surface area, and clinically relevant ADMET properties such as cardiotoxicity (hERG) and intestinal permeability (Caco-2). This adaptability spans from simple descriptors to practical pharmaceutical endpoints without requiring complete model retraining. The framework achieves precise control over target property values, generating molecules with properties closely matching specified targets for both single and multiple descriptors. Scaffold-constrained experiments preserve fixed molecular cores while maintaining effective property optimization. The three-dimensional representation offers advantages in maintaining structural validity during iterative optimization, with potential for geometry-aware applications in materials science and drug discovery.

3D molecular generation↗
DOE OSTI2026

Dielectric Spectroscopy Cable NDE for Unequal Aging Over Different Lengths

This Pacific Northwest National Laboratory milestone report assesses the effect of unequal aging over different lengths of cable on dielectric spectroscopy (DS) bulk impedance measurements from the cable end. DS measurements can indicate the remaining useful life of cables; however, most benchmark tests are based on accelerated aging of the entire cable. If only a portion of the cable is exposed to a harsh environment (as is frequently the case), the DS measurement will indicate a significantly less-aged cable than if the entire cable were exposed to a uniform environmental stress. The DS measurement is primarily related to the intrinsic insulation material relative permittivity—typically between 1.5 and 3.0. Insulation aging increases permittivity and, correspondingly, the cable’s overall capacitance measured from the cable end. To judge the health of the overall cable system, operators are primarily interested in the condition of the most severely aged section since that is where a failure is most likely to occur. This can be expressed as the dimensionless relative permittivity of the severely aged section divided by the permittivity of a pristine cable. This relative permittivity is an intrinsic material parameter independent of cable length and can serve as a quantitative indication of the insulation condition. If the relative percentages of cable exposed to a harsh and benign environment are known, the DS measurement of the entire cable length, including both pristine and aged sections, can be compensated to predict the relative permittivity (relative to pristine permittivity) of the most severly aged segment. A test was performed on three 50 ft (16 m) cables with 10%, 50%, and 90% of the cable inside a thermal aging oven. Predictably, the aging effect on the DS response was greatest in the 90% sample, followed by the 50% sample, and then the 10% sample. The more interesting result is the ratio of the aged insulation permittivity to the pristine insulation permittivity. Based on prior work, an end-of-life cable has a value of approximately 1.35. If the relative percentages of pristine and severely aged cable lengths are known, and the total cable-length capacitance of the initial pristine cable and the combined pristine plus aged cable are known, the relative aged permittivity (ratio of aged section permittivity / pristine section permittivity) can be estimated. The relative cable lengths exposed to a harsh environment may be known or estimated based on plant layout drawings or cable reflectometry tests that can locate the oven entry and exit points. The compensation was verified using a lumped-parameter model to predict the aged-segment permittivity e2/e1, where e2 is the aged-segment permittivity and e1 is the pristine-segment permittivity. For all three percentages of aged cable lengths, the ratio was quite similar, as would be expected since all cable segments were exposed to the same aging environment. This ratio of aged permittivity to pristine permittivity can be compared to other damage indicating tests, including elongation at break, to assess cable condition.

ARENA Test Bed↗
DOE OSTI2026

Evaluation of 3D pixel silicon sensors for the CMS Phase-2 Inner Tracker

The high-luminosity upgrade of the CERN LHC requires the replacement of the CMS tracking detector to cope with the increased radiation fluence while maintaining its excellent performance. An extensive R&D program, aiming at using 3D pixel silicon sensors in the innermost barrel layer of the detector, has been carried out by CMS in collaboration with the FBK (Trento, Italy) and CNM (Barcelona, Spain) foundries. The sensors will feature a pixel cell size of 25 × 100 µm 2 , with a centrally located electrode connected to the readout chip. The sensors are read out by the RD53A and CROCv1 chips, developed in 65 nm CMOS technology by the RD53 Collaboration, a joint effort between the ATLAS and CMS groups. This paper reports the results achieved in beam test experiments before and after irradiation, up to a fluence of approximately 2 . 6 × 1 0 16 n eq /cm 2 . Measurements of assemblies irradiated to a fluence of 1 × 10 16 n˙eq/cm 2 show a hit detection efficiency higher than 96% at normal incidence, with fewer than 2% of channels masked, across a bias voltage range greater than 50 V . Even after irradiation to a higher fluence of 1.6 × 10 16 n˙eq/cm 2 , similar performance is maintained over a bias voltage range of 30 V , remaining well within CMS requirements.

3D pixel↗
DOE OSTI2026

Emergence of low-energy spin waves in superconducting electron-doped cuprates

In order to fully utilize the technological potential of unconventional superconductors, an enhanced understanding of the superconducting mechanism is necessary. In the best performing superconductors, the cuprates, superconductivity is intimately linked with magnetism, although the details of this coupling remain elusive. Here, we address this gap by studying the electron-doped cuprate Nd 1.85 Ce 0.15 CuO 4−δ that has an antiferromagnetic ground state when synthesized and only becomes superconducting after a reductive annealing process. Using neutron spectroscopy, we show that the as-grown crystal exhibits a large spin pseudogap in the magnetic fluctuation spectrum. Annealing removes defects introduced by the commonly employed synthesis method and significantly reduces the spin pseudogap. While the spin pseudogap in the annealed sample likely arises from superconductivity, in the as-grown sample it results from the absence of long-wavelength spin waves. These results reveal a direct connection between defects, magnetism, and superconductivity, offering new insight into the mechanisms underlying high-temperature superconductivity and guiding the design of improved superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation↗
DOE OSTI2026

Mapping the Microhardness of Al Matrix in U-7Mo/Al Dispersion Fuels at Medium and High Burn-up by in situ Nanoindentation

The mechanical degradation of the Al matrix in U-7Mo dispersion fuels was detected after irradiation in previous studies, but more details of the degradation mechanism are not clear yet to clarify the root cause. In this study, microhardness mapping in fresh and irradiated U-7Mo dispersion fuels coated with ZrN, with burn-up of 3.35 × 10 21 and 6.28 × 10 21 f/cm 3 , are measured to reveal the hardening of the Al matrix as a function of distance to fuel particles. Micro cubes were pulled by plasma focused ion beam scanning electron microscopy and tested by in situ nanoindenter with Berkovitch tip. Size effects are measured in medium and high burn-up samples. The higher burn-up specimen exhibits a more pronounced size effect when indent load is lower than 20 mN. Size effect correction models are also calculated by plots fitting, which predicts the true microhardness without size effect. Results shows that the hardness of the Al matrix near the ZrN coating has the highest value, then decreases moving away from the coating and become stable when the distance reaches ∼10 µm. The hardness starts increasing again getting closer to the next fuel particle, which makes the hardness distribution between two fuel particles to be “U” shape.

36 - MATERIALS SCIENCE↗