Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 811 records · Page 45

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↗

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↗

Gas Hydrate Characterization in the GoM using Marine EM Methods

In spite of the importance of gas hydrate as a low-carbon fuel, a possible contributor to rapid climate change, and a significant natural hazard, our current understanding about the amount and distribution of submarine gas hydrate is somewhat poor; estimates of total volume vary by at least an order of magnitude, and commercially useful concentrations of hydrate have remained an elusive target. This is largely because conventional geophysical tools have intrinsic limitations in their ability to quantitatively image hydrate. It has long been known from well logs that gas hydrate is resistive compared to the host sediments, and electrical and electromagnetic methods have been proposed and occasionally used to image hydrates. This project seeks to expand our capabilities to use electromagnetic methods to explore for gas hydrate in the marine environment. An important basic science aspect of our work was to quantify the resistivity of pure gas hydrate as a function of temperature at seafloor pressures. We designed, constructed, and tested a highpressure cell in which hydrate could be synthesized and then subjected to electrical conductivity measurements. Impedance spectroscopy at frequencies between 20 Hz and 2 MHz was used to separate the effect of the blocking electrodes from the intrinsic conductivity of the hydrate. We obtained very reproducible results that showed that pure methane hydrate was several times more resistive than the water ice that seeded the synthesis, 20,000 {Ohm}m at 0{degrees} C, and that the activation energy is 30.6 kJ/mol over the temperature range of -15 to 15{degrees} C. Adding silica sand to the hydrate, however, showed that the addition of the extra phase caused the conductivity of the assemblage to increase in a counterintuitive way. The fact that the increased conductivity collapsed after a percolation threshold was reached, and that the addition of glass beads does not produce a similar increase in conductivity, together suggest that while the surface of the gas hydrate grains are not intrinsically conductive, the presence of sand does increase their conductivity. In the field component of this project, we carried out an 18day cruise on the R.V. Roger Revelle in the Gulf of Mexico from 7th-26th October 2008 to collect controlled-source electromagnetic (CSEM) data over four hydrate prospects; blocks AC 818, WR 313, GC 955, and MC 118. During these surveys we deployed 30 ocean bottom electromagnetic (OBEM) recorders a total of 94 times at four survey areas and towed the Scripps Undersea Electromagnetic Source Instrument (SUESI) a total of 103 hours. SUESI transmission was 200 A on a 50 m dipole antenna at heights of 70-100 m above the seafloor. We also towed a neutrally buoyant 3-axis electric field recorder behind the SUESI antenna at a constant offset of 300 m. The use of a towed receiver that is "flown" above the seafloor allowed us to operate in areas where seafloor infrastructure such as wellheads, pipelines, and installed scientific equipment existed. We reduced the data to apparent resistivity psuedosections. The most compelling results come from the hydrate observatory at MC 118, where a localized resistivity anomaly is clearly identified under the southeast crater in an otherwise uniform 1 {Ohm}m background. The data from MC 118 also show that authigenic carbonate does not necessarily express itself as a confounding resistor, as was feared at the start of this project. While the results from the other prospects are much more complicated, the data are well correlated with known geology, and line to line agreement is good. Although these data are not amenable to 1D inversion as was initially hoped, we expect to use a newly developed 2D CSEM inversion code to continue to get useful information from this rich data set.

03 NATURAL GAS↗

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↗

Test of a prototype neutron spectrometer based on diamond detectors in a fast reactor

A prototype of neutron spectrometer based on diamond detectors has been developed. This prototype consists of a {sup 6}Li neutron converter sandwiched between two CVD diamond crystals. The radiation hardness of the diamond crystals makes it suitable for applications in low power research reactors, while a low sensitivity to gamma rays and low leakage current of the detector permit to reach good energy resolution. A fast coincidence between two crystals is used to reject background. The detector was read out using two different electronic chains connected to it by a few meters of cable. The first chain was based on conventional charge-sensitive amplifiers, the other used a custom fast charge amplifier developed for this purpose. The prototype has been tested at various neutron sources and showed its practicability. In particular, the detector was calibrated in a TRIGA thermal reactor (LENA laboratory, University of Pavia) with neutron fluxes of 10{sup 8} n/cm{sup 2}s and at the 3 MeV D-D monochromatic neutron source named FNG (ENEA, Rome) with neutron fluxes of 10{sup 6} n/cm{sup 2}s. The neutron spectrum measurement was performed at the TAPIRO fast research reactor (ENEA, Casaccia) with fluxes of 10{sup 9} n/cm{sup 2}s. The obtained spectra were compared to Monte Carlo simulations, modeling detector response with MCNP and Geant4. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

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↗

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↗

Possible dynamical explanations for Paltridge's principle of maximum entropy production

Throughout the history of non-equilibrium thermodynamics a number of theories have been proposed in which complex, far from equilibrium flow systems are hypothesised to reach a steady state that maximises some quantity. Perhaps the most celebrated is Paltridge's principle of maximum entropy production for the horizontal heat flux in Earth's atmosphere, for which there is some empirical support. There have been a number of attempts to derive such a principle from maximum entropy considerations. However, we currently lack a more mechanistic explanation of how any particular system might self-organise into a state that maximises some quantity. This is in contrast to equilibrium thermodynamics, in which models such as the Ising model have been a great help in understanding the relationship between the predictions of MaxEnt and the dynamics of physical systems. In this paper we show that, unlike in the equilibrium case, Paltridge-type maximisation in non-equilibrium systems cannot be achieved by a simple dynamical feedback mechanism. Nevertheless, we propose several possible mechanisms by which maximisation could occur. Showing that these occur in any real system is a task for future work. The possibilities presented here may not be the only ones. We hope that by presenting them we can provoke further discussion about the possible dynamical mechanisms behind extremum principles for non-equilibrium systems, and their relationship to predictions obtained through MaxEnt.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigations of Dual-Purpose Canister Direct Disposal Feasibility (FY14)

Results reported here continue to support the FY13 conclusion that direct disposal of DPCs is technically feasible, at least for some DPCs, and for some disposal concepts (geologic host media). Much of the work performed has reached a point where site-specific information would be needed for further resolution. Several activities in FY14 have focused on clay/shale media because of potential complications resulting from low thermal conductivity, limited temperature tolerance, and the need to construct hundreds of kilometers of emplacement drifts that will remain stable for at least 50 years. Technologies for rapid excavation and liner installation have significantly advanced in the past 20 years. Tunnel boring machines are the clear choice for large-scale excavation. The first TBM excavations, including some constructed in clay or shale media, are now approaching 50 years of service. Open-type TBMs are a good choice but the repository host formation would need to have sufficient compressive strength for the excavation face to be self-supporting. One way to improve the strength-stress relationship is to reduce the repository depth in soft formations (e.g., 300 m depth). The fastest construction appears to be possible using TBMs with a single-pass liner made of pre-fabricated concrete segments. Major projects have been constructed with prefabricated segmented liner systems, and with cast-in-place concrete liners. Cost comparisons show that differences in project management and financing may be larger cost factors than the choice of liner systems. Costs for large-scale excavation and construction in clay/shale media vary widely but can probably be limited to $10,000 per linear meter, which is similar to previous estimates for repository construction. Concepts for disposal of DPC-based waste packages in clay/shale media are associated with thermal management challenges because of the relatively low thermal conductivity and limited temperature tolerance. Peak temperature limits of 100°C or lower for clay-rich materials have been selected by some international programs, but a limit above 100°C could help to shorten the duration of surface decay storage and repository ventilation. The effects of locally higher peak temperatures on repository performance need to be evaluated (in addition to the effects at lower temperatures). This report describes a modeling approach that couples the TOUGH2 and FLAC3D codes to represent thermally driven THM processes, as a demonstration of the types of models needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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.↗

A Prototype {sup 212}Pb Medical Dose Calibrator for Alpha Radioimmunotherapy

AREVA Med, an AREVA group subsidiary, is developing innovative cancer-fighting therapies involving the use of {sup 212}Pb for alpha radioimmunotherapy. Canberra Industries, the nuclear measurement subsidiary of AREVA, has been working with AREVA Med to develop a prototype measurement system to assay syringes containing a {sup 212}Pb solution following production by an elution system. The relative fraction of emitted radiation from the source associated directly with the {sup 212}Pb remains dynamic for approximately 6 hours after the parent is chemically purified. A significant challenge for this measurement application is that the short half-life of the parent nuclide requires assay prior to reaching equilibrium with progeny nuclides. A gross counting detector was developed to minimize system costs and meet the large dynamic range of source activities. Prior to equilibrium, a gross counting system must include the period since the {sup 212}Pb was pure to calculate the count rate attributable to the parent rather than the progeny. The dynamic state is determined by solving the set of differential equations, or Bateman equations, describing the source decay behavior while also applying the component measurement efficiencies for each nuclide. The efficiencies were initially estimated using mathematical modeling (MCNP) but were then benchmarked with source measurements. The goal accuracy of the system was required to be within 5%. Independent measurements of the source using a high resolution spectroscopic detector have shown good agreement with the prototype system results. The prototype design was driven by cost, compactness and simplicity. The detector development costs were minimized by using existing electronics and firmware with a Geiger-Mueller tube derived from Canberra's EcoGamma environmental monitoring product. The acquisition electronics, communications and interface were controlled using Python with the EcoGamma software development kit on a Raspberry Pi Linux computer mounted inside a standard project box. The results of initial calibration measurements are presented. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Inductively Driven, 3D Liner Compression of a Magnetized Plasma to Megabar Energy Densities

To take advantage of the smaller scale, higher density regime of fusion an efficient method for achieving the compressional heating required to reach fusion gain conditions must be found. What is proposed is a more flexible metallic liner compression scheme that minimizes the kinetic energy required to reach fusion. It is believed that it is possible to accomplish this at sub-megajoule energies. This however will require operation at very small scale. To have a realistic hope of inexpensive, repetitive operation, it is essential to have the liner kinetic energy under a megajoule which allows for the survivability of the vacuum and power systems. At small scale the implosion speed must be reasonably fast to maintain the magnetized plasma (FRC) equilibrium during compression. For limited liner kinetic energy, it becomes clear that the thinnest liner imploded to the smallest radius consistent with the requirements for FRC equilibrium lifetime is desired. The proposed work is directed toward accomplishing this goal. Typically an axial (Z) current is employed for liner compression. There are however several advantages to using a θ-pinch coil. With the θ-pinch the liner currents are inductively driven which greatly simplifies the apparatus and vacuum system, and avoids difficulties with the post implosion vacuum integrity. With fractional flux leakage, the foil liner automatically provides for the seed axial compression field. To achieve it with optimal switching techniques, and at an accelerated pace however will require additional funding. This extra expense is well justified as the compression technique that will be enabled by this funding is unique in the ability to implode individual segments of the liner at different times. This is highly advantageous as the liner can be imploded in a manner that maximizes the energy transfer to the FRC. Production of shaped liner implosions for additional axial compression can thus be readily accomplished with the modified power modules. The additional energy and switching capability proposed will thus provide for optimal utilization of the liner energy. The following tasks were outlined for the three year effort: (1) Design and assemble the foil liner compression test structure and chamber including the compression bank and test foils [Year 1]. (2) Perform foil liner compression experiments and obtain performance data over a range on liner dimensions and bank parameters [Year 2]. (3) Carry out compression experiments of the FRC plasma to Megagauss fields and measure key fusion parameters [Year 3]. (4) Develop numerical codes and analyze experimental results, and determine the physics and scaling for future work [Year 1-3]. The principle task of the project was to design and assemble the foil liner FRC formation chamber, the full compression test structure and chamber including the compression bank. This task was completed successfully. The second task was to test foils in the test facility constructed in year one and characterize the performance obtained from liner compression. These experimental measurements were then compared with analytical predictions, and numerical code results. The liner testing was completed and compared with both the analytical results as well as the code work performed with the 3D structural dynamics package of ANSYS Metaphysics®. This code is capable of modeling the dynamic behavior of materials well into the non-linear regime (e.g. a bullet hit plate glass). The liner dynamic behavior was found to be remarkably close to that predicted by the 3D structural dynamics results. Incorporating a code that can also include the magnetics and plasma physics has also made significant progress at the UW. The remaining test bed construction and assembly task is was completed, and the FRC formation and merging experiments were carried out as planned. The liner compression of the FRC to Megagauss fields was not performed due to not obtaining a sufficiently long lived FRC during the final year of the grant. Modifications planned to correct this deficiency included a larger FRC source as well as a much larger liner driver energy storage system. Due to discontinuation of the grant neither of these improvements were carried out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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↗

MO-PIS-Exhibit Hall-01: Imaging: CT Dose Optimization Technologies I

Partners in Solutions is an exciting new program in which AAPM partners with our vendors to present practical “hands-on” information about the equipment and software systems that we use in our clinics. The imaging topic this year is CT scanner dose optimization capabilities. Note that the sessions are being held in a special purpose room built on the Exhibit Hall Floor, to encourage further interaction with the vendors. Dose Optimization Capabilities of GE Computed Tomography Scanners Presentation Time: 11:15 – 11:45 AM GE Healthcare is dedicated to the delivery of high quality clinical images through the development of technologies, which optimize the application of ionizing radiation. In computed tomography, dose management solutions fall into four categories: employs projection data and statistical modeling to decrease noise in the reconstructed image - creating an opportunity for mA reduction in the acquisition of diagnostic images. Veo represents true Model Based Iterative Reconstruction (MBiR). Using high-level algorithms in tandem with advanced computing power, Veo enables lower pixel noise standard deviation and improved spatial resolution within a single image. Advanced Adaptive Image Filters allow for maintenance of spatial resolution while reducing image noise. Examples of adaptive image space filters include Neuro 3-D filters and Cardiac Noise Reduction Filters. AutomA adjusts mA along the z-axis and is the CT equivalent of auto exposure control in conventional x-ray systems. Dynamic Z-axis Tracking offers an additional opportunity for dose reduction in helical acquisitions while SmartTrack Z-axis Tracking serves to ensure beam, collimator and detector alignment during tube rotation. SmartmA provides angular mA modulation. ECG Helical Modulation reduces mA during the systolic phase of the heart cycle. SmartBeam optimization uses bowtie beam-shaping hardware and software to filter off-axis x-rays - minimizing dose and reducing x-ray scatter. The DICOM Radiation Dose Structured Report (RDSR) generates a dose report at the conclusion of every examination. Dose Check preemptively notifies CT operators when scan parameters exceed user-defined dose thresholds. DoseWatch is an information technology application providing vendor-agnostic dose tracking and analysis for CT (and all other diagnostic x-ray modalities) SnapShot Pulse improves coronary CTA dose management. VolumeShuttle uses two acquisitions to increase coverage, decrease dose, and conserve on contrast administration. Color-Coding for Kids applies the Broselow-Luten Pediatric System to facilitate pediatric emergency care and reduce medical errors. FeatherLight achieves dose optimization through pediatric procedure-based protocols. Adventure Series scanners provide a child-friendly imaging environment promoting patient cooperation with resultant reduction in retakes and patient motion. Philips CT Dose Optimization Tools and Advanced Reconstruction Presentation Time: 11:45 ‘ 12:15 PM The first part of the talk will cover “Dose Reduction and Dose Optimization Technologies” present in Philips CT Scanners. The main Technologies to be presented include: DoseRight and tube current modulation (DoseRight, Z-DOM, 3D-DOM, DoseRight Cardiac) Special acquisition modes Beam filtration and beam shapers Eclipse collimator and ClearRay collimator NanoPanel detector DoseRight will cover automatic tube current selection that automatically adjusts the dose for the individual patient. The presentation will explore the modulation techniques currently employed in Philips CT scanners and will include the algorithmic concepts as well as illustrative examples. Modulation and current selection technologies to be covered include the Automatic Current Selection component of DoseRight, ZDOM longitudinal dose modulation, 3D-DOM (combination of longitudinal and rotational dose modulation), Cardiac Dose right (an ECG based dose modulation scheme), and the DoseRight Index (DRI) IQ index. The special acquisition modes covers acquisition techniques such as prospective gating that is designed to reduce exposure to the patient through the Cardiac Step and Shoot scan mode. This mode can substitute the much higher dose retrospective scan modes for certain types of cardiac imaging. The beam filtration and beam shaper portion will discuss the variety of filtration and beam shaping configurations available on Philips scanners. This topic includes the x-ray beam characteristics, tube filtration as well as dose compensator characteristics. The Eclipse collimator, ClearRay collimator and the NanoPanel detector portion will discuss additional technologies specific to wide coverage CT that address some of the unique challenges encountered and techniques employed to optimize image quality and optimize dose utilization. The Eclipse collimator reduces extraneous exposure by actively blocking the radiation tails at either end of helical scans that do not contribute to the image generation. The ClearRay collimator and the NanoPanel detector optimize the quality of the signal that reaches the detectors by addressing the increased scattered radiation present in wide coverage and the NanoPanel detector adds superior electronic noise characteristics valuable when imaging at a low dose level. The second part of the talk will present “Advanced Reconstruction Technologies” currently available on Philips CT Scanners. The talk will cover filtered back projection (FBP), iDose4 and Iterative Model Reconstruction (IMR). Each reconstruction method will include a discussion of the algorithm as well as similarities and differences between the algorithms. Examples illustrating the merits of each algorithm will be presented, and techniques and metrics to characterize the performance of each type of algorithm will be presented. The Filtered Back projection portion will discuss and provide a brief summary of relevant standard image reconstruction techniques in common use, and discuss the common tradeoffs when using the FBP algorithm. The iDose4 portion will present the algorithms used for iDose4 as well the different levels. The meaning of different levels of iDose4 available will be presented and quantified. Guidelines for selection iDose4 parameters based on the imaging need will be explained. The different image quality goals available with iDose4 and specifically how iDose4 enables noise reduction, spatial resolution improvement or both will be explained. The approaches to leveraging the benefits of iDose4 such as improved spatial resolution, decreased noise, and artifact prevention will be described and quantified; and measurements and metrics behind the improvements will be presented. The image quality benefits in specific imaging situations as well as how to best combine the technology with other dose reduction strategies to ensure the best image quality at a given dose level will be presented. Insight into the IMR algorithm as well as contrast to the iDose4 techniques and performance characteristics will be discussed. Metrics and techniques for characterizing this class of algorithm and IQ performance will be presented. The image quality benefits and the dose reduction capabilities of IMR will be explored. Illustrative examples of the noise reduction, spatial resolution improvement, and low contrast detectability improvements of the reconstruction method will be presented: clinical cases and phantom measurements demonstrating the benefits of IMR in the areas of low dose imaging, spatial resolution and low contrast resolution are discussed and the technical details behind the measurements will be presented compared to both iDose4 and traditional filtered back projection (FBP)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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↗

SU-E-I-49: Simulation Study for Removing Scatter Radiation in Cesium-Iodine Based Flat Panel Detector System

Purpose: This study aims to identify the feasibility of a novel cesium-iodine (CsI)-based flat-panel detector (FPD) for removing scatter radiation in diagnostic radiology. Methods: The indirect FPD comprises three layers: a substrate, scintillation, and thin-film-transistor (TFT) layer. The TFT layer has a matrix structure with pixels. There are ineffective dimensions on the TFT layer, such as the voltage and data lines; therefore, we devised a new FPD system having net-like lead in the substrate layer, matching the ineffective area, to block the scatter radiation so that only primary X-rays could reach the effective dimension.To evaluate the performance of this new FPD system, we conducted a Monte Carlo simulation using MCNPX 2.6.0 software. Scatter fractions (SFs) were acquired using no grid, a parallel grid (8:1 grid ratio), and the new system, and the performances were compared.Two systems having different thicknesses of lead in the substrate layer—10 and 20μm—were simulated. Additionally, we examined the effects of different pixel sizes (153×153 and 163×163μm) on the image quality, while keeping the effective area of pixels constant (143×143μm). Results: In case of 10μm lead, the SFs of the new system (∼11%) were lower than those of the other system (∼27% with no grid, ∼16% with parallel grid) at 40kV. However, as the tube voltage increased, the SF of new system (∼19%) was higher than that of parallel grid (∼18%) at 120kV. In the case of 20μm lead, the SFs of the new system were lower than those of the other systems at all ranges of the tube voltage (40–120kV). Conclusion: The novel CsI-based FPD system for removing scatter radiation is feasible for improving the image contrast but must be optimized with respect to the lead thickness, considering the system’s purposes and the ranges of the tube voltage in diagnostic radiology. This study was supported by a grant(K1422651) from Institute of Health Science, Korea University.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗