Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIPS”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2026

High-Purity Quantum Emission from an Au 24 (S-CH 2 Ph- t Bu) 20 Nanocluster at Room Temperature

Atomically precise gold nanoclusters have garnered significant attention for their diverse applications, ranging from biological labeling to optoelectronics. Their potential in optical quantum computing, which calls for ideal single-photon sources, has recently become a key area of interest. In the current work, we use photon antibunching experiments to explore the single-photon emission efficiency of atomically precise Au 24 nanoclusters protected by 4-tert-butylbenzyl mercaptan ligands (Au 24 (TBBM) 20 ). This cluster exhibits quantum emission with good photostability and without any observable blinking or spectral drift at room temperature under an inert gas atmosphere, with antibunching dips (g 2 (0)) as low as 0.07 in the solid state or, equivalently, a single-photon purity of 93% under time-gated conditions. Transient absorption and time-gated antibunching studies reveal that the short emission lifetime of this cluster and its high photoluminescence quantum yield in the solid state play critical roles in enhancing the emitted single-photon purity. This research advances the understanding of single-emitter behavior in atomically precise gold nanoclusters, contributing to the development of stable quantum emitters that are essential for quantum computing and cryptography.

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

Excitation Laser Energy Dependence of the Gap-Mode TERS Spectra of WS 2 and MoS 2 on Silver

In this work we present a systematic study of the dependence of the gap mode tip-enhanced Raman scattering (TERS) response of the mono- and bi-layer WS 2 and MoS 2 on silver as a function of the excitation laser energy in a broad spectral range from 473nm to 830 nm. For this purpose, we collected consecutive TERS maps of the same area in the sample containing mono-and bi-layer regions with the same TERS probe with 6 different excitation lasers. To decrease the number of collected TERS maps, we used for the first time to the best of our knowledge, concurrent excitation and collection with two lasers simultaneously. We found that the E 2g /A 1g peak intensity ratio for the bilayer WS 2 @Ag and the ratio of the A’/A 1g peak intensity of the out-of-plane mode for the mono- and the bilayer, change in a significantly non-monotonous way as the excitation laser energy is swept from 1.58 eV to 2.62 eV. The former ratio increases at energies corresponding to A and B excitons (~2.0 eV and 2.4 eV correspondingly) in bilayer WS 2 . The absolute intensity of the A’ peak in the monolayer, and correspondingly the A’/A 1g ratio, is surprisingly high at lower excitation energies, but dips dramatically at the energy corresponding to the A exciton, being restored partially in between A and B excitons, but still showing the descending trend as the excitation laser energy increases. A somewhat similar picture was observed in mono- and bi-layers of MoS 2 @Ag, though the existing set of excitation lasers did not match the excitonic profile of this material as nicely as for the case of WS 2 . We attribute the observed behavior to the presence of intermediate (Fano resonance) or strong (Rabi splitting) coupling between the excitons in transition metal dichalcogenides (TMDs) and the plasmons in the tip-substrate nanocavity. This is akin to the so-called Fano (Rabi) transparency experimentally observed in far field scattering from transition metal dichalcogenides between two plasmonic metals. The possibility of the formation of intermediate/strong coupling between the excitonic resonances in TMDs and the nanocavity re-evaluates the role of various resonances in the gap-mode TERS and should become an important factor to be considered by TERS practitioners during planning the experiments. Lastly, based on observed phenomena and its explanation, we propose the “ideal” substrate for efficient TERS and tip enhanced photoluminescence (TEPL) measurements.

Raman spectroscopy↗
DOE OSTI2022

MPCVD-Grown Nanodiamond Microelectrodes with Oxygen Plasma Activation for Neurochemical Applications

Nanodiamonds (NDs) are a carbon nanomaterial that has a diamond core with heteroatoms and defects at the surface. The large surface area, defect sites, and functional groups on NDs make them a promising material for electrochemical sensing. Previously, we dip-coated ND onto carbon-fiber microelectrodes (CFMEs) and found increases in sensitivity, but the coating was sparse. Here, for this work, we directly grew thin films of ND on niobium wires using microwave plasma chemical vapor deposition (MP-CVD) to provide full surface coverage. ND microelectrodes show a reliable performance in neurotransmitter detection with good antifouling properties. To improve sensitivity, we oxygen plasma etched ND films to activate the surface and intentionally add defects and oxygen surface functional groups. For fast-scan cyclic voltammetry detection of dopamine, oxygen plasma-etching increases the sensitivity from 21 nA/μM to 90 nA/μM after treatment. Fouling was tested by repeated injections of serotonin or tyramine, and both ND and plasma oxidized nanodiamond (NDO) microelectrodes maintain their currents better compared to CFMEs and therefore are more antifouling. Furthermore, a biofouling test in brain slices shows that ND microelectrodes barely have any current drop, while the more hydrophilic NDO microelectrodes decrease more, but still not as much as CFMEs. Overall, grown ND microelectrodes are promising in neurotransmitter detection with excellent fouling resistance, whereas oxygen plasma etching slightly lowers the fouling resistance but dramatically increases sensitivity.

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

Remote Hydroamination of Disubstituted Alkenes by a Combination of Isomerization and Regioselective N–H Addition

Remote hydrofunctionalizations of alkenes incorporate functional groups distal to existing carbon–carbon double bonds. While remote carbonylations are well-known, remote hydrofunctionalizations are most common for addition of relatively nonpolar B–H, Si–H, and C–H bonds with alkenes. We report a system for the remote hydroamination of disubstituted alkenes to functionalize an alkyl chain selectively at the subterminal, unactivated, methylene position. Critical to the high regioselectivity and reaction rates are the electronic properties of the substituent on the amine and the development of the ligand DIP-Ad-SEGPHOS by evaluating the steric and electronic effects of ligand modules on reactivity and selectivity. The remote hydroamination is compatible with a broad scope of alkenes and aminopyridines and enables the regioconvergent synthesis of amines from an isomeric mixture of alkenes. In conclusion, the products can be derivatized by nucleophilic aromatic substitution on the amino substituent with a variety of nucleophiles.

Hydrocarbons↗
DOE OSTI2020

The Role of Premagmatic Rifting in Shaping a Volcanic Continental Margin: An Example From the Eastern North American Margin

Both magmatic and tectonic processes contribute to the formation of volcanic continental margins. Such margins are thought to undergo extension across a narrow zone of lithospheric thinning (~100 km). New observations based on existing and reprocessed data from the Eastern North American Margin contradict this hypothesis. With ~64,000 km of 2-D seismic data tied to 40 wells combined with published refraction, deep reflection, receiver function, and onshore drilling efforts, we quantified along-strike variations in the distribution of rift structures, magmatism, crustal thickness, and early post-rift sedimentation under the shelf of Baltimore Canyon Trough (BCT), Long Island Platform, and Georges Bank Basin (GBB). Results indicate that BCT is narrow (80–120 km) with a sharp basement hinge and few rift basins. The seaward dipping reflectors (SDR) there extend ~50 km seaward of the hinge line. In contrast, the GBB is wide (~200 km), has many syn-rift structures, and the SDR there extend ~200 km seaward of the hinge line. Early post-rift depocenters at the GBB coincide with thinner crust suggesting “uniform” thinning of the entire lithosphere. Models for the formation of volcanic margins do not explain the wide structure of the GBB. We argue that crustal thinning of the BCT was closely associated with late syn-rift magmatism, whereas the broad thinning of the GBB segment predated magmatism. Correlation of these variations to crustal terranes of different compositions suggests that the inherited rheology determined the premagmatic response of the lithosphere to extension.

58 GEOSCIENCES↗
DOE OSTI2022

Shallow Aseismic Slip in the Delaware Basin Determined by Sentinel-1 InSAR

The Delaware Basin, Texas is currently a hot-spot of induced seismicity and ground deformation due to fluid extraction and injection associated with horizontal drilling techniques; however, the driving mechanism behind the seismicity and deformation remains under debate. Here, using vertical and east-west horizontal surface deformation measurements derived from Sentinel-1 interferometric synthetic aperture radar (InSAR), we show that the subsurface responds differently to oil and gas activity in the northern and southeastern portions of the basin. In the north, where there is little seismicity, deformation patterns display long-wavelengths and equidimensional patterns. In contrast, the southeast region hosts most of the seismicity and displays spatial deformation patterns with narrow linear features that strike parallel to the maximum principal horizontal stress and to trends in seismicity, suggesting movement along normal faults. We model a linear deformation feature using edge dislocations and show that the InSAR observations can be reproduced by slip on normal faults contained within the Delaware Mountain Group (DMG), the formation that hosts local wastewater injection and the majority of earthquakes. Our model consists of three parallel, high-angle normal faults, with two dipping toward one another in a graben structure. Slip magnitudes reach up to 25 cm and are spatially correlated with injection wells. Measured seismicity can only explain ~2% of the fault motion predicted by our fault model, suggesting that slip leading to the deformation is predominantly aseismic. We conclude that seismic and aseismic fault motion in the southeastern Delaware Basin is likely driven by wastewater injection near critically-stressed normal faults within the DMG.

58 GEOSCIENCES↗
DOE OSTI2022

Bedrock Vadose Zone Storage Dynamics Under Extreme Drought: Consequences for Plant Water Availability, Recharge, and Runoff

Bedrock vadose zone water storage (i.e., rock moisture) dynamics are rarely observed but potentially key to understanding drought responses. Exploiting a borehole network at a Mediterranean blue oak savanna site—Rancho Venada—we document how water storage capacity in deeply weathered bedrock profiles regulates woody plant water availability and groundwater recharge. The site is in the Northern California Coast Range within steeply dipping turbidites. In a wet year (water year 2019; 647 mm of precipitation), rock moisture was quickly replenished to a characteristic storage capacity, recharging groundwater that emerged at springs to generate streamflow. In the subsequent rainless summer growing season, rock moisture was depleted by about 93 mm. In two drought years that followed (212 and 121 mm of precipitation) the total amount of rock moisture gained each winter was about 54 and 20 mm, respectively, and declines were documented exceeding these amounts, resulting in progressively lower rock moisture content. Oaks, which are rooted into bedrock, demonstrated signs of water stress in drought, including reduced transpiration rates and extremely low water potentials. In the 2020–2021 drought, precipitation did not exceed storage capacity, resulting in variable belowground water storage, increased plant water stress, and no recharge or runoff. Rock moisture deficits (rather than soil moisture deficits) explain these responses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Influence of Fault Architecture on Induced Earthquake Sequence Evolution Revealed by High-Resolution Focal Mechanism Solutions

The increasing seismicity and improved seismic observation network in recent years provide an opportunity to explore factors that influence the triggering processes, spatiotemporal evolution, and maximum magnitude of induced sequences. We map the fault architecture and stress state of four induced sequences in Oklahoma to determine their influence on the seismicity. We systematically relocate the earthquakes and compute hundreds of focal mechanisms of small to medium events (1.0 < M < 5.1) using various techniques, including machine learning, for the Guthrie, Woodward, Cushing, and Fairview sequences. The detailed fault geometry and spatiotemporal evolution of seismicity and stress states reveal different dominant driving forces for each sequence. In Cushing and Fairview (largest event ≥M5.0), the main fault structures are near-vertical narrow strike-slip faults, with most of the small earthquake fault planes optimally oriented. The two sequences exhibit discontinuous temporal migration but strong earthquake self-driven rupture growth. In Guthrie and Woodward (largest event <M5.0), the two sequences show more complex diffuse fault structures with varying dipping angles along depth. The inverted focal mechanisms show a mix of strike-slip faulting and normal faulting in both sequences, and the normal faulting events are less optimally oriented than strike-slip events. The two sequences are dominated by continuous diffusive migration in time driven by pore pressure propagation. The above results suggest that fault architecture and stress state influence sequence evolution, major driving forces, and possibly maximum magnitude.

58 GEOSCIENCES↗
DOE OSTI2023

The Variable Continuous Bimaterial Interface in the San Jacinto Fault Zone Revealed by Dense Seismic Array Analysis of Fault Zone Head Waves

Key factors controlling earthquake ruptures include fault geometry, continuity, and seismic velocity structure around the fault. We present a novel tool that better informs deep bimaterial fault geometry embedded in distributed damage and seismicity, associated velocity contrasts across the fault, and their correlations with surface complexities. The method employs fault zone head and direct body waves and is applied to recordings from five spatiotemporally different seismic arrays along the complex San Jacinto fault zone (SJFZ) in southern California. We detect and distinguish these signals based on instantaneous phase coherence and relative energy in a cascading manner from one scale array to another. The analysis reveals a >70-km long continuous bimaterial interface within the SJFZ with several deep northeast dipping fault segments. The northern SJFZ, for instance, locates ~7 km northeast of its surface expression at 18-km depth. P-wave velocity contrasts range from near 0% to >15%, consistent with other bimaterial faults, and differ by a few % depending on fault-array azimuth, implying directional-dependent velocity contrasts. S-wave head waves and velocity contrasts are also imaged for the first time at the southern SJFZ, averaging to 2.9% in agreement with tomography results. The imaged geometry and continuity suggest the SJFZ initiated along remnant tectonic structures and translates to a rupture potential of M > 7.2, i.e., the sizes of its largest paleo-earthquakes. The P and S contrasts, and their ratios, have important implications for earthquake rupture speed, mode, directivity, and frictional heating along the SJFZ and other major faults globally.

58 GEOSCIENCES↗
DOE OSTI2023

Enhanced Microseismicity During Production Pumping Cessation at the San Emidio Geothermal Field (Nevada, USA) in December 2016

Abstract Tectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

High Values of the Arctic Amplification in the Early Decades of the 21st Century: Causes of Discrepancy by CMIP6 Models Between Observation and Simulation

Arctic Amplification (AA) in the first decade of the 21st century has reached values between 4 and 5, with a subsequent decrease to current values of about 3.6, while the value was from 2 to 3 during the twentieth century. The ensemble mean of the CMIP6 models has difficulty in reproducing the recently observed high values of the AA. In this report, we identify the main reason for this difficulty to be the CMIP6 models overestimate of the mean global temperature trend since about 1990. The largest values of the AA are observed in winter and spring. A sharp AA peak in 1987 spring was caused by a peak in the Arctic temperature trend occurring at the same time as a dip in the trend of mean global temperature. The winter AA has increased almost monotonically since 1990. Dividing the AA between the Arctic land and ocean areas shows that the ocean area makes a larger contribution to the AA. Our future projection of the AA suggests an increasing AA for about the next decade, followed by a slow decrease to about 3.5 in the 2050s.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.

Environmental Sciences & Ecology↗
DOE OSTI2026

Sequential Fracture Activation and Stress Evolution During EGS Stimulation at Utah FORGE Revealed by Waveform Cross‐Correlation

Mapping fracture networks in Enhanced Geothermal Systems (EGS) is essential for optimizing reservoir performance, yet complex fracture evolution during stimulation remains difficult to resolve. This study examines the evolution of microseismicity and fracture networks during stage 3 of the 2022 EGS stimulation at the Utah Frontier Observatory for Research in Geothermal Energy site. We map the fracture network represented by 20 clusters of seismic events identified by waveform similarities with cross-correlation. We characterize their geometric properties such as strike, dip, length, and width, and analyze the time evolution of activated fractures. The results reveal a systematic fracture evolution: early activation of pre-existing natural fractures, complex network development during peak injection, and continued activation of less favorably oriented fractures post-injection. Magnitude calibration using the Principal Component Analysis of cross-correlated waveforms improves relative amplitude measurements, refining estimations of the Gutenberg-Richter b-values with spatial variations in b-values suggesting stress re-distribution across the stimulated area. Analysis of the stress state of selected fractures further shows that fractures requiring higher excess pore pressure primarily activate at the end of injection and post-injection, highlighting stress transfer due to pore pressure as a dominant triggering mechanism. These findings provide insights into fracture propagation, stress evolution, and seismic hazard assessment in EGS reservoirs.

Asirifi, Richard [Texas A & M Univ., College Stati↗
DOE OSTI2026

From Points to Planes: A Workflow for Converting Three‐Dimensional Point Cloud Data Into Discrete Fracture Network Flow and Transport Models

We present the Point cLoud Algorithm for NEtwork Extraction of Discrete Fracture Networks (PLANE-DFN), a point cloud–based algorithm for automatic fracture network extraction designed to support discrete fracture network (DFN) modeling workflows. PLANE-DFN segments three-dimensional fracture planes from raw point cloud data using RANdom SAmple Consensus coupled with statistical outlier removal and density-based clustering to isolate individual fracture features. Each candidate plane is constrained against site-specific structural constraints based on strike and dip. After segmentation, each fracture is converted into a 2-D convex polygon suitable for meshing and simulation. The PLANE-DFN algorithm is validated by comparing geometric and flow and transport data against data from dfnWorks simulations with ensembles of plane-fit networks. We find that the flow and transport in plane-fit networks are comparable to dfnWorks-generated networks when realistic network geometry is maintained. The PLANE-DFN algorithm provides an automated and streamlined workflow to transform point clouds of data into DFN network geometry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effect of tectonic processes on biosphere–geosphere feedbacks across a convergent margin

The subsurface is among Earth’s largest biomes, but the extent to which microbial communities vary across tectonic plate boundaries or interact with subduction-scale geological processes remains unknown. Here we compare bacterial community composition with deep-subsurface geochemistry from 21 hot springs across the Costa Rican convergent margin. We find that cation and anion compositions of the springs reflect the dip angle and position of the underlying tectonic structure and also correlate with the bacterial community. Co-occurring microbial cliques related to cultured chemolithoautotrophs that use the reverse tricarboxylic acid cycle (rTCA) as well as abundances of metagenomic rTCA genes correlate with concentrations of slab-volatilized carbon. This, combined with carbon isotope evidence, suggests that fixation of slab-derived CO2 into biomass may support a chemolithoautotrophy-based subsurface ecosystem. We calculate that this forearc subsurface biosphere could sequester 1.4 × 10 9 to 1.4 × 10 10 mol of carbon per year, which would decrease estimates of the total carbon delivered to the mantle by 2 to 22%. Furthermore, based on the observed correlations, we suggest that distribution and composition of the subsurface bacterial community are probably affected by deep tectonic processes across the Costa Rican convergent margin and that, by sequestering carbon volatilized during subduction, these chemolithoautotrophic communities could in turn impact the geosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Introduction of Chalcogenide Glasses to Additive Manufacturing: Nanoparticle Ink Formulation, Inkjet Printing, and Phase Change Devices Fabrication

Abstract Chalcogenide glasses are one of the most versatile materials that have been widely researched because of their flexible optical, chemical, electronic, and phase change properties. Their application is usually in the form of thin films, which work as active layers in sensors and memory devices. In this work, we investigate the formulation of nanoparticle ink of Ge–Se chalcogenide glasses and its potential applications. The process steps reported in this work describe nanoparticle ink formulation from chalcogenide glasses, its application via inkjet printing and dip-coating methods and sintering to manufacture phase change devices. We report data regarding nanoparticle production by ball milling and ultrasonication along with the essential characteristics of the formed inks, like contact angle and viscosity. The printed chalcogenide glass films were characterized by Raman spectroscopy, X-ray diffraction, energy dispersive spectroscopy and atomic force microscopy. The printed films exhibited similar compositional, structural, electronic and optical properties as the thermally evaporated thin films. The crystallization processes of the printed films are discussed compared to those obtained by vacuum thermal deposition. We demonstrate the formation of printed thin films using nanoparticle inks, low-temperature sintering and proof for the first time, their application in electronic and photonic temperature sensors utilizing their phase change property. This work adds chalcogenide glasses to the list of inkjet printable materials, thus offering an easy way to form arbitrary device structures for optical and electronic applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗