Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elemental profiling”

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

Steady state thermo-mechanics and material property definition framework for analyzing DCLL blanket in the fusion nuclear science facility

In this work, a thermo-mechanics model that relies on creating the material property definition framework (MPDF) and multiphysics coupling of the heat transfer and the solid mechanics modules is developed to determine the structural integrity of the recently designed dual cooled lead lithium (DCLL) inboard blanket (IB) for the Fusion Nuclear Science Facility under steady state loads. The MPDF is called to supply fusion relevant neutron irradiation and temperature induced changes in material properties during multiphysics finite element runs, and PbLi temperature profiles are used to approximate Magnetohydrodynamics effect and the nuclear volumetric heating on the PbLi. Neutron irradiation and temperature induced reduction of the yield and ultimate strengths of F82H steel at the first wall (FW) are quantified for one year. A blanket in an assembly with gaps between blanket sectors and another blanket in an assembly with no gaps between blanket sectors, both exposed to radiation damage that lasted for one year are analyzed. Analysis using the elastic ITER structural design criteria for in-vessel components (ITER SDC-IC) design rules and a linear isotropic-hardening-type elastoplastic material model are used where most appropriate. The IB blanket with gaps between blanket sectors will withstand the steady state combined thermal and coolant loads for one year operational period but will fail if no gaps are allowed between blanket sectors. It is recommended that a gap of about 7.62 mm should be provided between IB blanket sectors during assembly which would close up during service, stop neutron streaming, reduce stresses and reduce bending of the FW into the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing the properties and mechanisms of failure waves in soda-lime glass

Soda-lime glass (SLG) and other silica glasses exhibit the failure wave phenomenon under shock compression. The mechanism responsible for this peculiar behavior of glasses is still unresolved. In this study, a series of plate impact experiments was performed at three different impact stresses of 6.4, 8.3, and 10.8 GPa to better understand the mechanisms underlying the failure wave phenomenon. Specifically, spall experiments were conducted to probe the speed and existence of failure waves at different stresses in SLG. A layered glass target was used to probe the possibility of a “renucleation” of the wave at the SLG–SLG interface. When it existed, the failure wave was inferred to propagate at a speed of 1.3 km/s. However, it was observed that the failure wave phenomenon ceases to exist for impact stresses higher than 10 GPa. In experiments with a 6.4 GPa impact stress, the peak free surface velocity was significantly less than what is predicted by stress-Hugoniot calculations. This velocity deficit and other important features of the measured free surface velocity profiles were simulated using finite element analysis by incorporating an abrupt densification of SLG at a critical stress in the equation of state. Furthermore, this densification feature is similar to what would be expected of a phase transition. Although unable to unambiguously reveal the mechanism causing the failure wave phenomenon, the results of the present work clearly indicate that the failure wave causes a secondary compression and densification in SLG.

36 MATERIALS SCIENCE↗

Various facets of magnetic charge correlation: Micromagnetic and distorted-wave Born approximation simulations study

The emergent concept of the magnetic charge quasiparticle provides a new realm to study the evolution of magnetic properties in two-dimensional artificially frustrated magnets. Here we report on the exploration of magnetic phases due to various magnetic charge correlation using the complementary numerical techniques of micromagnetic and distorted-wave Born approximation simulations in artificial permalloy honeycomb lattice. The honeycomb element length varies between 10 nm and 100 nm, while the width and thickness are kept within the single domain limit. In addition to the charge ordered loop state, we observe disordered charge arrangement, characterized by the random distribution of ±Q charges, in single domain size honeycomb lattice. As the length of the honeycomb element increases, low multiplicity magnetic charges tend to form contiguous bands in thinner lattice. Thin honeycomb lattice with 100 nm element length exhibits a perfect spin ice pattern, which remains unaffected by the modest increase in the width of element size. We simulate scattering profiles under the pretext of distorted-wave Born approximation formalism for the micromagnetic phases. The results are expected to provide useful guidance in the experimental investigation of magnetic phases in an artificial honeycomb magnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Digital Twin Approach to Study Additive Manufacturing Processing Using Embedded Optical Fiber Sensors and Numerical Modeling

One of the major challenges for metal-powder-based additive manufacturing is measuring and mitigating residual strain induced during the manufacturing processes. This article reports distributed fiber optic sensors embedded in Inconel alloy components as experimental means to validate numerical models of additive manufacturing process. Electroplating was used to deposit a metal protective jacket onto standard telecom single-mode fibers for strain measurements, Fiber sensors were embedded in an Inconel alloy substrate using the laser engineered net shaping (LENS) process. Here by using a Rayleigh-scattering optical frequency domain reflectometer (OFDR), temperature changes, and residual strain in the metal substrate were monitored with 5 mm spatial resolution during the LENS process. Using finite element analysis, temperature and strain profiles induced by the LENS deposition processes were also numerically studied. Discrepancies between the simulated temperature and strain profiles and those measured directly were less than 10%. Results presented in this article demonstrates a digital twin approach to fuse modeling results with distributed fiber sensor measurement data to study additive manufacturing process toward design and fabrication process optimization.

36 MATERIALS SCIENCE↗

THE DESIGN OF LAYER TIME OPTIMIZATION IN LARGE SCALE ADDITIVE MANUFACTURING WITH FIBER REINFORCED POLYMER COMPOSITES

In this study, we have developed a method to optimize a layer deposition time (a.k.a. layer time) for large-scale additive manufacturing (AM) via physics-based simulations. A long layer time leads to an over-cooled surface on which a new layer is deposited, and therefore, it may result in a weak bonding or debonding between layers, cracking, or warping. A short layer time leads to a high temperature of the structure due to insufficient cooling, and therefore, the structure may not be stiff enough and may collapse during manufacturing. Therefore, it is important to estimate the optimal layer time in additive manufacturing for a high-quality product. The temperature of a top layer right before deposition is recommended to be slightly higher than the glass temperature of the material. A temperature cooling was approximated to an exponential function of time, and the optimized layer time was obtained based on a target temperature while maintaining a minimal printing time. The material used is carbon fiber-reinforced polycarbonate (CF/PC), and the large-scale deposition system used is LSAM TM from Thermwood Corporation. Three different layer time cases were used for experiments, and a series of thermal images were obtained via an infra-red (IR) camera during the entire AM processes. AM process simulations were performed using a finite element method and the temperature profiles from the simulation were in good agreement with those from experiments. The layer time optimization was performed based on the temperature profiles from the simulations. A layer temperature with the optimal layer time was confirmed as the target temperature through simulation.

Jo, Eonyeon↗

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗

High-resolution Compton spectroscopy using x-ray microcalorimeters

X-ray Compton spectroscopy is one of the few direct probes of the electron momentum distribution of bulk materials in ambient and operando environments. Herein, we report high-resolution inelastic x-ray scattering experiments with high momentum and energy transfer performed at a storage-ring-based high-energy x-ray light source facility using an x-ray transition-edge sensor (TES) microcalorimeter detector. The performance was compared with a silicon drift detector (SDD), an energy-resolving semiconductor detector, and Compton profiles were measured for lithium and cobalt oxide powders relevant to lithium-ion battery research. Spectroscopic analysis of the measured Compton profiles demonstrates the high-sensitivity to the low- Z elements and oxidation states. The line shape analysis of the measured Compton profiles in comparison with computed Hartree–Fock profiles is usually limited by the resolution of the semiconductor detector. We have characterized an x-ray TES microcalorimeter detector for high-resolution Compton scattering experiments using a bending magnet source at the Advanced Photon Source with a double crystal monochromator, providing monochromatic photon energies near 27.5 keV. The momentum resolution below 0.16 atomic units (a.u.) was measured, yielding an improvement of more than a factor of 7 over a state-of-the-art SDD for the same scattering geometry. Furthermore, the lineshapes of narrow valence and broad core electron profiles of sealed lithium metal were clearly resolved using an x-ray TES compared to smeared and broadened lineshapes observed when using the SDD. High-resolution Compton scattering using the energy-resolving area detector shown here presents new opportunities for spatial imaging of electron momentum distributions for a wide class of materials with applications ranging from electrochemistry to condensed matter physics.

47 OTHER INSTRUMENTATION↗

Effects of midplane density gradient in the propagation of high-harmonic fast waves considering high temperature anisotropy in the scrape-off layer of NSTX-U

High-Harmonic Fast Wave (HHFW) heating experiments in NSTX have shown that up to 60% of the injected power can be lost in the Scrape-Off Layer (SOL) when the density is above the fast wave cutoff density in front of the antenna, for which the fast wave is able to propagate into the plasma. This work models HHFW propagation in the SOL plasmas of NSTX-U using a 2D divertor SOL profile derived from the pressure assumption and the finite element temperature solution, which accounts for the high anisotropy of heat conduction in a 2D axisymmetric geometry. In this work, the two-dimensional axisymmetric SOL temperature profile is first evaluated by solving the steady-state non-linear heat conduction equation, in which thermal conductivity depends on temperature, using a finite element approach in the Petra-M workbench. A 2D density profile is then obtained from a prescribed density profile at the outer midplane assuming pressure is uniform along a flux tube in the SOL. This approach results in density and temperature profiles in which the strong asymmetric nature of heat conduction between the parallel and perpendicular background magnetic field is successfully captured. Furthermore, this work focuses on investigating the effect of the SOL plasma density profile on parasitic HHFW propagation in the SOL. The simulation results show that the radial gradient of the density profile affects the wavefield propagation in the SOL. As the density profile broadens, the wavefield intensity and its poloidal extent are reduced in the SOL, and the core coupling increases. Similarly, as shown by collisional power deposition, a proxy for power absorption, the fraction of power deposited in the SOL decreases with a broadening profile.

HHFW↗

Probing lanmodulin's mechanisms of rare-earth selectivity for protein-based bioseparations

Our BES Separation Science program project, DE-SC0021007, supported our efforts to begin to understand the mechanisms underlying selectivity of a novel class of lanthanide-binding proteins discovered by our laboratory, called lanmodulin (LanM), and to leverage these proteins for recovery and separations of trivalent rare earth elements (REEs) as well as of trivalent actinides. Overall, our work provides important insights into how higher-order (e.g., secondary, tertiary, and quaternary) protein structure modulates selectivity profiles of proteins that bind f-elements highly selectively. These results are important for advancing the concept of protein-based separations of REEs and, perhaps, of other critical minerals.

Lanmodulin, rare earth elements, protein-based met↗

Evaluating the diffusion of Kr in UO 2 and ADOPT TM using time-of-flight elastic recoil detection analysis (ToF-erda)

A combination of 300 keV 84 Kr ion implantation and Time-of-Flight Elastic Recoil Detection Analysis is utilized to investigate the diffusion of Kr in UO 2 and ADOPT TM fuels. Composition depth-profiles on the nanometer scale were obtained, both for as-implanted samples and after annealing at 800°C for 1 hour. Observed drifts in the 84 Kr profiles could be associated with short-range diffusion mechanisms. The approach employed here provides the possibility to make direct comparisons with atomistic scale modelling data, and can be of service as a separate effect test in line with the Accelerated Fuel Qualification initiative.

ADOPT UO2↗

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗

System, apparatus, and method for determining elemental composition using 4D STEM

The present disclosure relates to transmission electron microscopy for evaluation of biological matter. According to an embodiment, the present disclosure further relates to an apparatus for determining the structure and/or elemental composition of a sample using 4D STEM, comprising a direct bombardment detector operating with global shutter readout, processing circuitry configured to acquire images of bright-field disks using either a contiguous array or non-contiguous array of detector pixel elements, correct distortions in the images, align each image of the images based on a centroid of the bright-field disk, calculate a radial profile of the images, normalize the radial profiles by a scaling factor, calculate the rotationally-averaged edge profile of the bright-field disk, and determine elemental composition within the specimen based on the characteristics of the edge profile of the bright-field disk corresponding to each specimen location.

Bammes, Benjamin↗

Targeted synthesis of predicted metastable compounds using modulated elemental reactants

Three metastable compounds predicted to be kinetically stable using an “island” approach were successfully synthesized from designed modulated elemental reactants. Fe0.8V0.2Se2 was synthesized by depositing ultrathin elemental layers in a V|Fe|Se sequence to control the local composition. An alloyed rock salt structured Pb3Mn2Se5 constituent layer, which does not exist as a bulk compound, was synthesized in the heterostructure (Pb3Mn2Se5)0.6VSe2 by depositing a precursor with a V|Se|Pb|Se|Mn|Se|Pb|Se|Mn|Se|Pb|Se sequence of elemental layers that mimicked the compositional profile of the targeted heterostructure. The heterostructure (PbSe)1+δ(FeSe2)2 was prepared by depositing a precursor with a repeating layering sequence of Fe|Pb|Fe|Se, where each sequence contains the number of atoms required to form a single unit cell. In all three systems, the local compositions in the layer sequence kinetically favored the nucleation and growth of the targeted products during the deposition. The diffusion lengths to form the targeted compounds were short, and the diffusion was limited by postdeposition low temperature annealing to favor the growth of the targeted compounds and avoid the decomposition into a mixture of thermodynamically stable compounds.

Lemon, Mellie (ORCID:0000000328705881)↗

All-ceramic channel waveguides fabricated via 3D printing

All-ceramic channel waveguides (CWGs) in Yb:YAG transparent ceramics have been fabricated for the first time, to the best of our knowledge, via direct ink write (DIW) and their laser performance has been demonstrated. Single filaments of Yb:YAG nanoparticle-loaded ink were extruded into undoped YAG; the Yb:YAG filaments formed the CWGs, surrounded by undoped cladding. Elemental mapping confirmed the Yb doping profile and waveguide integrity. Optical characterization showed low cladding scatter losses (<1.3%/cm at 1.3 µm), and laser testing with a 940 nm Ti:sapphire pump demonstrated efficient lasing at 1030 nm. The best-performing waveguide, with an elliptical cross-section (100 µm × 60 µm and a length of 1.4 cm), achieved a slope efficiency of 61% and a roundtrip loss of 12.4%. These results identify DIW as a promising approach for fabricating high-performance channel waveguides in transparent ceramics.

36 MATERIALS SCIENCE↗

Section-level genome sequencing and comparative genomics of Aspergillus sections Cavernicolus and Usti.

The genus Aspergillus is diverse, including species of industrial importance, human pathogens, plant pests, and model organisms. Aspergillus includes species from sections Usti and Cavernicolus, which until recently were joined in section Usti, but have now been proposed to be non-monophyletic and were split by section Nidulantes, Aenei and Raperi. To learn more about these sections, we have sequenced the genomes of 13 Aspergillus species from section Cavernicolus (A. cavernicola, A. californicus, and A. egyptiacus), section Usti (A. carlsbadensis, A. germanicus, A. granulosus, A. heterothallicus, A. insuetus, A. keveii, A. lucknowensis, A. pseudodeflectus and A. pseudoustus), and section Nidulantes (A. quadrilineatus, previously A. tetrazonus). We compared these genomes with 16 additional species from Aspergillus to explore their genetic diversity, based on their genome content, repeat-induced point mutations (RIPs), transposable elements, carbohydrate-active enzyme (CAZyme) profile, growth on plant polysaccharides, and secondary metabolite gene clusters (SMGCs). All analyses support the split of section Usti and provide additional insights: Analyses of genes found only in single species show that these constitute genes which appear to be involved in adaptation to new carbon sources, regulation to fit new niches, and bioactive compounds for competitive advantages, suggesting that these support species differentiation in Aspergillus species. Sections Usti and Cavernicolus have mainly unique SMGCs. Section Usti contains very large and information-rich genomes, an expansion partially driven by CAZymes, as section Usti contains the most CAZyme-rich species seen in genus Aspergillus. Section Usti is clearly an underutilized source of plant biomass degraders and shows great potential as industrial enzyme producers. Citation: Nybo JL, Vesth TC, Theobald S, Frisvad JC, Larsen TO, Kjaerboelling I, Rothschild-Mancinelli K, Lyhne EK, Barry K, Clum A, Yoshinaga Y, Ledsgaard L, Daum C, Lipzen A, Kuo A, Riley R, Mondo S, LaButti K, Haridas S, Pangalinan J, Salamov AA, Simmons BA, Magnuson JK, Chen J, Drula E, Henrissat B, Wiebenga A, Lubbers RJM, Müller A, dos Santos Gomes AC, Mäkelä MR, Stajich JE, Grigoriev IV, Mortensen UH, de Vries RP, Baker SE, Andersen MR (2025). Section-level genome sequencing and comparative genomics of Aspergillus sections Cavernicolus and Usti. Studies in Mycology 111: 101-114. doi: 10.3114/sim.2025.111.03.

59 BASIC BIOLOGICAL SCIENCES↗

Long range frequency chirping of Alfvén eigenmodes

In this work, a theoretical framework has been developed for an NBI scenario to model the hard non-linear evolution of global Alfvén eigenmodes (GAEs) where the adiabatic motion of phase-space structures (holes and clumps), associated with frequency chirping, occurs in generalized phase-space of slowing down energetic particles. The radial profile of the GAE is expanded using finite elements, which allows an update of the mode structure as the mode frequency chirps. Constants of motion are introduced to track the dynamics of energetic particles during frequency chirping by implementing proper action-angle variables and canonical transformations which reduce the dynamics essentially to 1D. Consequently, we specify whether the particles are drifting inward/outward as the frequency deviates from the initial MHD eigenfrequency. Using the principle of least action, we have derived the non-linear equation describing the evolution of the radial profile by varying the total Lagrangian of the system with respect to the weights of the finite elements. For the choice of parameters in this work, it is shown that the peak of the radial profile is shifted and also broadens due to frequency chirping. The time rate of frequency change is also calculated using the energy balance and we show that the adiabatic condition remains valid once it is satisfied. This model clearly illustrates the theoretical treatment to study the long-range adiabatic frequency sweeping events observed for Alfvén gap modes in real experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TCF Base Technology-Specific: Real-time soil core scanner for carbon and other elements (Final Report)

The aim of this project was to commercialize an instrument that will accurately and precisely measure carbon in soil, providing a breakthrough in scalable soil core analysis and enabling a robust measurement of the total amount of carbon per acre. The instrument will measure the whole soil core, rather than just a few sub-samples, and provide data for carbon depth profiling, as well as other key soil elements such as Al, Si, O, Fe. This instrument will provide faster turnaround times and more precise results than current methods (such as loss on ignition), facilitating the analysis of many soil cores per acre. Ultimately, this will enable farmers to improve their carbon-in-soil content, thus contributing to an improved soil health and product yield. To achieve this goal, we divided the project into three main tasks: 1) commercialization research, 2) technical development, and 3) instrument development, demonstration, and benchmarking.

47 OTHER INSTRUMENTATION↗

As-Built Simulation of the High Flux Isotope Reactor

The Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) is an 85 MWt flux trap-type research reactor that supports key research missions, including isotope production, materials irradiation, and neutron scattering. The core consists of an inner and an outer fuel element containing 171 and 369 involute-shaped plates, respectively. The thin fuel plates consist of a U 3 O 8 -Al dispersion fuel (highly enriched), an aluminum-based filler, and aluminum cladding. The fuel meat thickness is varied across the width of the involute plate to reduce thermal flux peaks at the radial edges of the fuel elements. Some deviation from the designed fuel meat shaping is allowed during manufacturing. A homogeneity scan of each fuel plate checks for potential anomalies in the fuel distribution by scanning the surface of the plate and comparing the attenuation of the beam to calibration standards. While typical HFIR simulations use homogenized fuel regions, explicit models of the plates were developed under the Low-Enriched Uranium Conversion Program. These explicit models typically include one inner and one outer fuel plate with nominal fuel distributions, and then the plates are duplicated to fill the space of the corresponding fuel element. Therefore, data extracted from these simulations are limited to azimuthally averaged quantities. To determine the reactivity and physics impacts of an as-built outer fuel element and generate azimuthally dependent data in the element, 369 unique fuel plate models were generated and positioned. This model generates the three-dimensional (i.e., radial–axial–azimuthal) plate power profile, where the azimuthal profile is impacted by features within the adjacent control element region and beryllium reflector. For an as-built model of the outer fuel element, plate-specific homogeneity data, 235 U loading, enrichment, and channel thickness measurements were translated into the model, yielding a much more varied azimuthal power profile encompassed by uncertainty factors in analyses. These models were run with the ORNL-TN and Shift Monte Carlo tools, and they contained upwards of 500,000 cells and 100,000 unique tallies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗