SEARCH · Engineering Papers

Results for “reach code”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Added value for integrated marine energy data systems

Launched in 2019, the Portal and Repository for Information on Marine Renewable Energy (PRIMRE) provides centralized access, standardization, community building, and integration of United States (U.S.) databases, tools and codes, and other resources that cover a range of marine energy information. The PRIMRE universe contains a series of Knowledge Hubs that represent data and information on testing of marine energy devices (MHKDR); environmental effects (Tethys); engineering and technical papers (Tethys Engineering); descriptions of marine energy companies and technologies (Marine Energy Projects Database); codes and models (Marine Energy Software); and guidance on testing and measurements (Telesto). Content is added to PRIMRE by applying a set of Guidelines and Best Practices (PRIMRE Guidelines). An aggregate search across the PRIMRE site enables users to find data and information from all the PRIMRE Knowledge Hubs simultaneously, using a single entry-point (PRIMRE Search). In addition to providing access to a range of data and information on marine energy development, testing, and effects, PRIMRE allows for the development of value added products and processes that will help move the marine energy industry forward. The PRIMRE team has recently launched two key initiatives in the U.S.—Signature Projects and Lessons Learned. Outputs and outcomes from these two initiatives will be highlighted in this paper. The Signature Projects initiative is intended to bring focus to a selection of ongoing and completed marine energy projects funded by the U.S. Department of Energy’s Water Power Technologies Office (WPTO), and to inform the marine energy community of what investigations have been undertaken, what tools are available, and where gaps in information persist. Each Signature Project tags papers, reports, and data from large marine energy research projects, providing easy access and attention to all the output and associated products from each project. The Lessons Learned initiative is intended to ensure that hard-won achievements are recognized and available for those who come later, that missteps and unfortunate outcomes can be prevented in future, and that efficiencies and effective shortcuts can be publicized and used as the marine energy industry moves forward. This initiative builds off the Knowledge Hubs and reaches out to members of the marine energy community, particularly technology developers and researchers, to integrate experience in the development, deployment, assessment, success, and challenges while creating this new industry and field of study.

Copping, Andrea E.↗

DOE OSTI2021

Added value for integrated marine energy data systems

Copping, Andrea E.↗

DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗

DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2017

Hybrid-PIC modeling of laser-plasma interactions and hot electron generation in gold hohlraum walls

The walls of the hohlraum used in experiments at the national ignition facility are heated by laser beams with intensities ∼1015 W/cm2, a wavelength of ∼1/3 μm, and pulse lengths on the order of a ns, with collisional absorption believed to be the primary heating mechanism. X-rays generated by the hot ablated plasma at the gold walls are then used to implode a target in the hohlraum interior. In addition to the collisional absorption of laser energy at the walls, non-linear laser-plasma interactions (LPI), such as stimulated Raman scattering and two plasmon decay, are believed to generate a population of supra-thermal electrons which, if present in the hohlraum, can have a deleterious effect on target implosion. We describe results of hohlraum modeling using a hybrid particle-in-cell code. To enable this work, new particle-based algorithms for a multiple-ion magneto-hydrodynamic (MHD) treatment, and a particle-based ray-tracing model were developed. The use of such hybrid methods relaxes the requirement to resolve the laser wavelength, and allows for relatively large-scale hohlraum simulations with a reasonable number of cells. But the non-linear effects which are believed to be the cause of hot electron generation can only be captured by fully kinetic simulations with good resolution of the laser wavelength. For this reason, we employ a two-tiered approach to hohlraum modeling. Large-scale simulations of the collisional absorption process can be conducted using the fast quasi-neutral MHD algorithm with fluid particle species. From these simulations, we can observe the time evolution of the hohlraum walls and characterize the density and temperature profiles. From these results, we can transition to smaller-scale highly resolved simulations using traditional kinetic particle-in-cell methods, from which we can fully model all of the non-linear laser-plasma interactions, as well as assess the details of the electron distribution function. We find that vacuum hohlraums should be stable to both two plasmon decay and stimulated Raman scattering instabilities for intensities ≤1015 W/cm2. In gas-filled hohlraums, shocks may be induced in the blowoff gold plasma, which leads to more complex density and temperatures profiles. The resulting effect on LPI stability depends strongly on the details of the profile, and it is possible for the gas-filled hohlraum to become unstable to two plasmon decay at 1015 W/cm2 if the quarter-critical surface reaches temperatures exceeding 1 keV.

Thoma, C.↗

NASA NTRS2011

Dust Processing in Supernova Remnants: Spitzer MIPS SED and IRS Observations

We present Spitzer MIPS SED and IRS observations of 14 Galactic Supernova Remnants previously identified in the GLIMPSE survey. We find evidence for SNR/molecular cloud interaction through detection of [OI] emission, ionic lines, and emission from molecular hydrogen. Through black-body fitting of the MIPS SEDs we find the large grains to be warm, 29-66 K. The dust emission is modeled using the DUSTEM code and a three component dust model composed of populations of big grains, very small grains, and polycyclic aromatic hydrocarbons. We find the dust to be moderately heated, typically by 30-100 times the interstellar radiation field. The source of the radiation is likely hydrogen recombination, where the excitation of hydrogen occurred in the shock front. The ratio of very small grains to big grains is found for most of the molecular interacting SNRs to be higher than that found in the plane of the Milky Way, typically by a factor of 2--3. We suggest that dust shattering is responsible for the relative over-abundance of small grains, in agreement with prediction from dust destruction models. However, two of the SNRs are best fit with a very low abundance of carbon grains to silicate grains and with a very high radiation field. A likely reason for the low abundance of small carbon grains is sputtering. We find evidence for silicate emission at 20 $\mu$m in their SEDs, indicating that they are young SNRs based on the strong radiation field necessary to reproduce the observed SEDs.

Hewitt, John W.↗

DOE OSTI2017

Advanced Simulation Aided Metallic Fuel Design

Nuclear engineering is reaching a domain space in which computational power and methods are starting to surpass the models that are used to simulate nuclear fuel rods. As these models are updated, simulations utilizing different compositions, operating conditions, fabrication techniques, and geometries can provide core designers with an important first step towards fuel qualification. Due to the cost associated with fuel irradiations, it can be expected that time and money spent building better fuel models may provide an attractive return during initial studies of advanced fuels, especially with regards to fuels beyond the extensively studied uranium oxide system. In general, the primary U.S. experience with irradiated fuel rods beyond traditional oxide fuel is with binary (U-Zr) and ternary (U-Pu-Zr) metal fuels. The bulk of the experimental data for metal fuels is obtained from rods irradiated in the EBR-II sodium fast reactor (SFR), thus many of the correlations developed for metal fuels comes from the post-irradiation (PIE) examination of EBR-II rods. One of the unique characteristics that set U-Pu-Zr fuels apart from traditional oxide fuels is the presence of different phases that develop in the fuel during irradiation. Each of these different phases has unique thermal, mechanical, and diffusive behavior. The transitions between these phases are not static as the initially homogenous zirconium diffuses towards the center and edge of the fuel rod, resulting in a low-Zr, high-U ring. Along with the differing material properties that occur in each phase, the increase in uranium corresponds to an increase in localized fission and power production. The direct coupling between the thermal and diffusion solves highlights one of the many examples of tightly coupled phenomenon that are present in nuclear fuel, and motivation for the development of codes that simultaneously solve each piece of physics in a fully-coupled manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2016

A Numerical Approach to Estimate the Ballistic Coefficient of Space Debris from TLE Orbital Data

Low Earth Orbit (LEO) is full of space debris, which consist of spent rocket stages, old satellites and fragments from explosions and collisions. As of 2009, more than 21,000 orbital debris larger than 10 cm are known to exist], and while it is hard to track anything smaller than that, the estimated population of particles between 1 and 10 cm in diameter is approximately 500,000, whereas small as 1 cm exceeds 100 million. These objects orbit Earth with huge kinetic energies speeds usually exceed 7 kms. The shape of their orbit varies from almost circular to highly elliptical and covers all LEO, a region in space between 160 and 2,000 km above sea level. Unfortunately, LEO is also the place where most of our active satellites are situated, as well as, International Space Station (ISS) and Hubble Space Telescope, whose orbits are around 400 and 550 km above sea level, respectively.This poses a real threat as debris can collide with satellites and deal substantial damage or even destroy them.Collisions between two or more debris create clouds of smaller debris, which are harder to track and increase overall object density and collision probability. At some point, the debris density couldthen reach a critical value, which would start a chain reaction and the number of space debris would grow exponentially. This phenomenon was first described by Kessler in 1978 and he concluded that it would lead to creation of debris belt, which would vastly complicate satellite operations in LEO. The debris density is already relatively high, as seen from several necessary debris avoidance maneuvers done by Shuttle, before it was discontinued, and ISS. But not all satellites have a propulsion system to avoid collision, hence different methods need to be applied. One of the proposed collision avoidance concepts is called LightForce and it suggests using photon pressure to induce small orbital corrections to deflect debris from colliding. This method is very efficient as seen from theoretical simulations, even few continuous mode 10 kW ground-based lasers, focused by 1.5 m telescopes with adaptive optics, were enough to prevent significant amount of the debris collisions. Simulations were done by propagating all space objects in LEO by 1 year into the future and checking whether the probability of collision was high. For those space objects different ground-based lasers were used to divert them, afterwards collision probabilities were reevaluated. However, the actual accuracy of the LightForce software, which has been developed at NASA AmesResearch Center, depends on the veracity of the input parameters, one of which is the objects ballistic coefficient. It is a measure of bodys ability to overcome air resistance, which has a significant impact on the debris in LEO, and thus it is responsible for the shape of the trajectory of the debris. Having the exact values of the ballistic coefficient would make significantly better collision predictions, unfortunately, we do not know what are the values for most of the objects.In this research, we were working with part of LightForce code, which estimates the ballistic coefficient from ephemerides. Previously used method gave highly inaccurate values, when compared to known objects, and it needed to be changed. The goal of this work was to try out a different method of estimating the ballistic coefficient and to check whether or not it gives noticeable improvements.

Collisions↗

DOE OSTI2019

Pushered single shell implosions for mix and radiation trapping studies using high-Z layers on National Ignition Facility

Pushered Single Shells (PSSs) are another approach to Inertial Confinement Fusion implosions that employ high-Z materials in the innermost capsule layer (pusher) as a means to enhance radiation trapping and lower core ignition requirements. However, adding high-Z materials can also increase losses due to mix, provide extra tamping, and make the capsule emission opaque to x-ray diagnostics. The first PSS implosions performed on the National Ignition Facility use plastic ablators with a germanium (Ge) dopant as a high-Z surrogate in the pusher to isolate the effects of high-Z mix and radiation trapping without changing tamping. Using a 2-shock laser pulse, the PSS implosions are designed and symmetrized to reach 3.7 keV core temperatures. A low concentration (2.8%) Ge dopant is added to the innermost layer, and the resulting effects on mix and x-ray opacity are observed. The method of separated reactants is used to infer information about mixing between the deuterated plastic pusher and the capsule fill gas (25% tritium) from the resulting nuclear DT reactions. Radiation transport is studied via capsule emission x-ray spectroscopy and imaging. Both nuclear and x-ray data corroborate the hypothesis that the addition of Ge strongly affects the mix region through radiation losses but has a minimal effect on the core and the warm, unmixed regions. Simulations using diffusive and turbulent mix models agree qualitatively with data, but quantitative agreement may require hybrid mix models that can model the transitional regime between turbulence and diffusion. Simulations matching the observables show increased core radiation trapping when Ge is added.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗

DOE OSTI2017

Spanning the Core Design Space with Binary Perturbations - A Mathematical Investigation

A reactor core design engineer, even one with many years of experience or one using advanced fuel pattern optimization software, will likely fine tune the pattern before finalizing the design. Fine tuning can further minimize local enthalpy peaks, optimize reactivity to extend burnup, and add additional margin to limits on reactivity coefficients or heat flux. For reactors that have reached an equilibrium fuel cycle, the number of feeds each cycle is typically the similar to previous cycles as is their placement in the core. It is then the engineer's responsibly to efficiently select the placement of burned assemblies to ensure optimal fuel utilization, satisfy safety limits, and meet cycle energy requirements. One of the engineer's most commonly used tools in this fine tuning investigation is to swap the positions of two assemblies, most likely two previously irradiated ('burned') assemblies since their reactivity is lower than fresh fuel. Such an operation is referred to in this work as a binary perturbation or a 'swap'. If quarter core symmetry is maintained, the perturbation translates to the other three quadrants. Likewise if eighth core symmetry is maintained, the perturbation translates to other seven octants. The question is then, how much of the core design space can binary perturbations possibly span? This work seeks to quantitatively define the space spanned by m number of binary perturbations as a function of variable size through combinatorial principles and computer simulation. The results of this study are useful in understanding the in-core fuel management problem faced by both professionals who employ hand manipulations and those who attempt stochastic solutions of the problem.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Sandia Cloud Forensics and Incident Response Platform v 1.0

As systems and devices become virtualized and deployed in the cloud, the hypervisor becomes an increasingly appropriate place to collect performance data, system state, system landscape, function calls, transaction traces, and other characteristics. We developed a method by which an introspection application may be coupled with a hypervisor to “reach into” the VM with minimal intrusiveness to collect data critical to the reconstruction of events, files, and operations. Such a capability is required to take advantage of the hypervisor as an instrumentation platform and to integrate that data with more traditional collection mechanisms. The concept of a VM serviced by a lightweight hypervisor is a relatively new paradigm for forensic practitioners. Traditional forensic techniques, based on assumptions that the filesystem was directly interacting with the hardware through an abstraction, afforded the forensic practitioner the assumption that there was nothing controlling the application below the filesystem. This is not the case when using virtualized technologies. Hypervisors can covertly monitor, introspect and interact with the guest in a transparent fashion. The problems of storage and collection of actionable data are exhausting. The current challenge is most hypervisors do not expose a useful application programming interface (API) at a sufficient level to do transparent, fine-grained and customizable introspection. Scalable VM instrumentation and introspection at an in-depth level requires fast handling of events, as well as direct access to VM state. Furthermore, deep introspection benefits greatly from the ability to gather data from the hardware during the VM’s exit to the hypervisor. All of this requires identical access to the system as the hypervisor itself; improper use of this ability could easily cause system instability. It is for this reason we believe that the hypervisor developers have been hesitant to grant this much control through their APIs. However, our approach leverages other means to collect and monitor the guests in a targeted fashion. Virtual machine introspection (VMI) is a technique used to monitor the runtime state of a system-level virtual machine. The runtime state can be may include processor registers, memory, disk, network, and any other hardware-level events A review of research literature and current VMI technologies exposed several limitations and tradeoffs in VMI approaches, including: the use of in-guest agents; kernel to user space transitions (dramatically slowing down processing); VMI tool pre-configuration requirements; hypervisor version lock-in or source code patching; reliance on operating system symbols; limited processor features due to hypervisor (even if the hardware could do more). To address these constraints, a VMI tool was envisioned to provide the cloud forensic capabilities while having as few of these limitations as possible. The Kernel-based Virtual Machine Introspection (KVMi) tool was developed for hypervisors on Intel’s x86-64 architecture. To meet performance, scoping and use-case demands, the follow criteria were applied to KVMi to support live, forensic data collection: • Shall not require in-guest agents. • Shall work with any recent version of modern hypervisors. • Shall introspect upon any guest (VM) running modern versions of Windows, Linux and OS X. • Shall have the ability to find and track operating system artifacts • Shall not require operating system symbol files. • Shall be able to fully handle VM-exits, bypassing execution of the hypervisor if necessary, to facilitate new features the hypervisor may not support. • Shall be compilable/loadable on a running system with standard build tools.

Stout, WilliamM.S.↗

DOE OSTI2024

Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM

The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.

58 GEOSCIENCES↗

DOE OSTI2021

AGR-5/6/7 Irradiation As Run Predictions Using PARFUME

The PARticle FUel ModEl (PARFUME), a fuel performance modeling code used for high-temperature gas-cooled reactors, was used to model the Advanced Gas Reactor (AGR)-5/6/7 irradiation test using as-run physics and thermal data. The AGR-5/6/7 irradiation test consists of the combined fifth, sixth, and seventh planned irradiations of the AGR Fuel Development and Qualification Program. The AGR-5/6/7 test train is a multi-capsule, instrumented experiment that is designed for irradiation in the 133.4-mm diameter northeast flux trap position of the Advanced Test Reactor (ATR) at Idaho National Laboratory. Each capsule contains compacts filled with uranium oxycarbide unaltered fuel particles. This report documents the calculations performed to predict the failure probability of tristructural isotropic (TRISO)-coated fuel particles during the AGR-5/6/7 experiment. In addition, this report documents the calculated fission product release fraction from the fuel. The calculations include modeling of the AGR 5/6/7 irradiation that occurred from February 2018 to July 2020 over nine ATR cycles, including six normal cycles and three power axial locator mechanism cycles, for a total of approximately 376 effective full power days (EFPD). The irradiation conditions and material properties of the AGR-5/6/7 test predicted zero fuel particle failures in Capsules 1, 3, and 4. Fuel particle failures were predicted in two of the compacts in Capsule 2 and one particle failure is predicted in each one of the compacts in Capsule 5. All compacts that exhibited fuel particle failures predicted by PARFUME were caused by localized stress concentrations in the silicon carbide (SiC) layer caused by cracking in the inner pyrolytic (IPyC) layer. In addition, shrinkage of the buffer and IPyC layer during irradiation resulted in formation of a buffer-IPyC gap. Compacts with a lower irradiation temperature and fluence experienced the smallest buffer-IPyC gap formation. Conversely, higher irradiated temperature compacts with a high fluence experienced the largest buffer-IPyC gap formation. Compact 3-6-3 experienced the largest buffer IPyC gap formation of just under 21.7 µm. The release fraction of fission products silver (Ag), cesium (Cs), and strontium (Sr) vary depending on capsule location and irradiation temperature. The maximum release fraction of Ag occurs in Capsule 3, reaching up to 59.5% for the TRISO fuel particles (compact 3-6-3). The release fraction of the other two fission products, Cs and Sr, are much smaller. A maximum Cs release fraction of 1.1% occurred in compact 3-4-3 and 4.4% for Sr in compact 3-6-3.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2018

Aerial Lidar and Geomorphic Change Detection Software- Tools for Monitoring Canyon Sediments at LANL - 18668

The Los Alamos National Laboratory Environmental Remediation program has conducted extensive geomorphic surveys of sediment deposits in the canyon bottoms within and surrounding the facility. The goal of the geomorphic surveys was to provide basic information needed to characterize chemical inventories of the sediments downstream of historical release sites. Sediment units were mapped based on depositional setting and age. Representative characterization samples were collected for each unit and submitted for chemical analysis. Chemical inventories were estimated by calculating the total volume of each depositional unit and multiplying by average chemical concentration. This method of characterization provides a spatial distribution of chemical data based on fluvial processes. In general, for the canyon bottom sediments, the higher concentrations of industrial chemicals were associated with fine-grained floodplain deposits closest to the source. Since initial characterization, many of the reaches have undergone changes in response to post-wildfire floods and severe rain events. These events impacted the channel and floodplain deposits. Much of the elevated concentrations are confined to the over-bank deposits distal from the primary channel. In order to ensure that the sediment units containing elevated chemical concentrations are stable, repeat aerial light detection and ranging (LiDAR) surveys are conducted and changes are assessed using the Geomorphic Change Detection (GCD) plug-in for ArcGIS. These tools provide high confidence uncertainty estimates and extensive areal coverage. This is a valuable method for long-term monitoring of surface sediments with low-level but widespread contaminant distributions. All of the monitored canyon bottom reaches pass risk assessment and present technically infeasible sediment removal options. Following the 2014 monsoon season, aerial LiDAR data has been collected annually and the resulting digital elevation models (DEMs) processed using the GCD method. Changes are identified using GCD and are verified by field walk downs. If erosion is identified in areas which are vulnerable to mobilizing contaminated sediment, enhanced controls can be installed to stabilize banks or other areas with high erosion potential. Annual sediment monitoring samples are collected where there is evidence for flow and deposition within the past year. Since 2014, no severe floods have been catalogued using the LiDAR surveys and GCD software. The recent rainfall/runoff events have not caused major changes to the canyon bottom channels and floodplains. Geomorphic processes identified in the monitoring results include channel aggradation and incision that over the course of a monsoon season results in minor changes to the system. Other active processes contributing to observed changes are minor slides, flows, slumps, and falls of unconsolidated sediment on steep bedrock or soil surfaces. A comprehensive baseline has been established by geomorphic mapping and characterization sampling. Periodic monitoring consists of LiDAR surveys and GCD analysis which offers an effective and defensible methodology for documenting stability of surface sediments. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2017

AGR-5/6/7 Irradiation Test Predictions using PARFUME

PARFUME, (PARticle FUel ModEl) a fuel performance modeling code used for high temperature gas-cooled reactors (HTGRs), was used to model the Advanced Gas Reactor (AGR)-5/6/7 irradiation test using predicted physics and thermal hydraulics data. The AGR-5/6/7 test consists of the combined fifth, sixth, and seventh planned irradiations of the AGR Fuel Development and Qualification Program. The AGR-5/6/7 test train is a multi-capsule, instrumented experiment that is designed for irradiation in the 133.4-mm diameter north east flux trap (NEFT) position of Advanced Test Reactor (ATR). Each capsule contains compacts filled with uranium oxycarbide (UCO) unaltered fuel particles. This report documents the calculations performed to predict the failure probability of tristructural isotropic (TRISO)-coated fuel particles during the AGR-5/6/7 experiment. In addition, this report documents the calculated source term from the driver fuel. The calculations include modeling of the AGR-5/6/7 irradiation that is scheduled to occur from October 2017 to April 2021 over a total of 13 ATR cycles, including nine normal cycles and four Power Axial Locator Mechanism (PALM) cycle for a total between 500 – 550 effective full power days (EFPD). The irradiation conditions and material properties of the AGR-5/6/7 test predicted zero fuel particle failures in Capsules 1, 2, and 4. Fuel particle failures were predicted in Capsule 3 due to internal particle pressure. These failures were predicted in the highest temperature compacts. Capsule 5 fuel particle failures were due to inner pyrolytic carbon (IPyC) cracking causing localized stresses concentrations in the SiC layer. This capsule predicted the highest particle failures due to the lower irradiation temperature. In addition, shrinkage of the buffer and IPyC layer during irradiation resulted in formation of a buffer-IPyC gap. The two capsules at the two ends of the test train, Capsules 1 and 5 experienced the smallest buffer-IPyC gap formation due to the lower irradiation fluences and temperatures. Capsule 3 experienced the largest buffer-IPyC gap formation of just under 24 µm. The release fraction of fission products Ag, Cs, and Sr silver (Ag), cesium (Cs), and strontium (Sr) vary depending on capsule location and irradiation temperature. The maximum release fraction of Ag occurs in Capsule 3, reaching up to 84.8% for the TRISO fuel particles. The release fraction of the other two fission products, Cs and Sr are much smaller and, in most cases, less than 1%. The notable exception is again in Capsule 3, where the release fraction for Cs and Sr reach up to 9.7% and 19.1%, respectively.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗