Search NASA⌕ Search

SEARCH · Search NASA

Results for “SCALE Code System”

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 1,387 records · Page 77

A Novel Auto-Scaling MPPT Algorithm based on Perturb and Observe Method for Photovoltaic Modules under Partial Shading Conditions

Maximum Power Point Tracking (MPPT) controller is required in a solar photovoltaic (PV) system to deliver the maximum power to load from PV module. This paper proposes a novel stepped MPPT method to realize a simple MPPT controller, which can track the real maximum power point (RMPP) even under partial shading conditions. The proposed algorithm is started by scanning the characteristic curve of the PV modules to detect the global maximum power point and then the algorithm will be switched to the conventional P&O algorithm to track the true maximum power point. The obtained simulation results, using Power electronic simulation software (PSIM), are compared with those found using the P&O method to confirm the performance of our proposed MPPT method even under non-uniform solar irradiation.

14 SOLAR ENERGY↗

Information criteria for quantifying loss of reversibility in parallelized KMC

Parallel Kinetic Monte Carlo (KMC) is a potent tool to simulate stochastic particle systems efficiently. However, despite literature on quantifying domain decomposition errors of the particle system for this class of algorithms in the short and in the long time regime, no study yet explores and quantifies the loss of time-reversibility in Parallel KMC. Inspired by concepts from non-equilibrium statistical mechanics, we propose the entropy production per unit time, or entropy production rate, given in terms of an observable and a corresponding estimator, as a metric that quantifies the loss of reversibility. Typically, this is a quantity that cannot be computed explicitly for Parallel KMC, which is why we develop a posteriori estimators that have good scaling properties with respect to the size of the system. Through these estimators, we can connect the different parameters of the scheme, such as the communication time step of the parallelization, the choice of the domain decomposition, and the computational schedule, with its performance in controlling the loss of reversibility. From this point of view, the entropy production rate can be seen both as an information criterion to compare the reversibility of different parallel schemes and as a tool to diagnose reversibility issues with a particular scheme. As a demonstration, we use Sandia Lab's SPPARKS software to compare different parallelization schemes and different domain (lattice) decompositions.

97 MATHEMATICS AND COMPUTING↗

Resiliency of Degraded Built Infrastructure

Infrastructure resiliency depends on the ability of infrastructure systems to withstand, adapt, and recover from chronic and extreme stresses. In this white paper, we address the resiliency of infrastructure assets and discuss improving infrastructure stability through development of our understanding of cement and concrete degradation. The resiliency of infrastructure during extreme events relies on the condition, adaptability, and recoverability of built infrastructure (roads, bridges, dams), which serves as the backbone of existing infrastructure systems. Much of the built infrastructure in the US has consistently been rated D+ by the American Society of Civil Engineers (ASCE). Aged infrastructure introduces risk to the system, since unreliable infrastructure increases the likelihood of failures under chronic and extreme stress and are particularly concerning when extreme events occur. To understand and account for this added risk from poor infrastructure quality, more research is needed on (i) how the changing environment alters the aging of new and existing built infrastructure and (ii) how degradation causes unique failure mechanisms. The aging of built infrastructure is based on degradation of the structural materials, such as concrete and steel supports, which causes failure. Current work in cement/concrete degradation is based on (i) the development of high strength and degradation resistance concrete mixtures, (ii) methods of assessing the age and reliability of existing structures, and (3) modeling of structural stability and the microstructural evolution of concrete/cement from degradation mechanisms (sulfide attack, carbonation, decalcification). Sandia National Laboratories (SNL) has made several investments in studying the durability and degradation of cement based materials, including using SNL-developed codes and methodologies (peridynamics, PFLOTRAN) to focus on chemo-mechanical fracture of cement for energy applications. Additionally, a recent collaboration with the University of Colorado Boulder has included fracture of concrete gravity dams, scaling the existing work to applications in full sized infrastructure problems. Ultimately, SNL has the experience in degradation of cementitious materials to extend the current research portfolio and answer concerns about the resilience of aging built infrastructure.

36 MATERIALS SCIENCE↗

Isotopic and Fuel Lattice Parameter Trends in Extended Enrichment and Higher Burnup LWR Fuel Vol I: PWR fuel

Commercial light water reactor (LWR) operators and fuel vendors in the United States are pursuing changes to nuclear fuel that include extended enrichment (EE) and accident-tolerant fuel (ATF) designs. The term EE (8% > 235 U > 5%) is used in this report to refer to a subset of high assay low-enriched uranium (HALEU) that is considered usable in commercial US LWRs in the near term. ATF features are designed to improve fuel system performance under accident conditions. One goal of EE is to improve fuel cycle economy by enabling fuel to be depleted to higher burnup than the typical current maximum pin burnup limits (62 gigawatt-days per metric ton of uranium [GWd/MTU]). Adoption of EE, ATF, and high burnup (HBU) fuels in the US commercial fleet requires a clear understanding of the effects on core physics parameters and used fuel isotopic content, as well as confidence in the accuracy of computer code predictions over an expanded range of materials, enrichment, and burnup. A thorough understanding of the applicability and adequacy of benchmark data (e.g., criticality, decay heat, isotopic content) for computer code validation is necessary to ensure that appropriate safety margins are maintained. To prepare for and support these potential changes, the effects of EE, ATF, and HBU are being assessed for selected representative LWR fuel designs. The project is divided into phases: this report summarizes the findings of Phase 1, which focuses on the lattice physics parameter and used fuel isotopic changes for a conventional Westinghouse 17×17 pressurized water reactor (PWR) design. The primary investigation tool is the SCALE Polaris code using the SCALE 56-group Evaluated Nuclear Data File (ENDF)/B-VII.1 cross sections. The goal of the current work is to (1) identify and explain important effects of EE and HBU (reactivity, lattice physics, and isotopic effects) assuming that PWR fuel design and usage remain similar to those for current enrichment fuel, (2) provide limited code-to-code comparisons with higher order cross section libraries and/or codes, and (3) identify any apparent anomalous trends in the results for further investigation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Preconditioning for the Navier-Stokes equations with finite-rate chemistry

The extension of Van Leer's preconditioning procedure to generalized finite-rate chemistry is discussed. Application to viscous flow is begun with the proper preconditioning matrix for the one-dimensional Navier-Stokes equations. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from nearly stagnant flow to hypersonic. Specific benefits are realized at the low and transonic flow speeds typical of complete propulsion-system simulations. The extended preconditioning matrix necessarily accounts for both thermal and chemical nonequilibrium. Numerical analysis reveals the possible theoretical improvements from using a preconditioner for all Mach number regimes. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number areas. Van Leer, Lee, and Roe recently developed an optimal, analytic preconditioning technique to reduce eigenvalue stiffness over the full Mach-number range. By multiplying the flux-balance residual with the preconditioning matrix, the acoustic wave speeds are scaled so that all waves propagate at the same rate, an essential property to eliminate inherent eigenvalue stiffness. This session discusses a synthesis of the thermochemical nonequilibrium flux-splitting developed by Grossman and Cinnella and the characteristic wave preconditioning of Van Leer into a powerful tool for implicitly solving two and three-dimensional flows with generalized finite-rate chemistry. For finite-rate chemistry, the state vector of unknowns is variable in length. Therefore, the preconditioning matrix extended to generalized finite-rate chemistry must accommodate a flexible system of moving waves. Fortunately, no new kind of wave appears in the system. The only existing waves are entropy and vorticity waves, which move with the fluid, and acoustic waves, which propagate in Mach number dependent directions. The nonequilibrium vibrational energies and species densities in the unknown state vector act strictly as convective waves. The essential concept for extending the preconditioning to generalized chemistry models is determining the differential variables which symmetrize the flux Jacobians. The extension is then straight-forward. This algorithm research effort will be released in a future version of the production level computational code coined the General Aerodynamic Simulation Program (GASP), developed by Walters, Slack, and McGrory.

Godfrey, Andrew G.↗

Crack Growth Rate and Large Plate Demonstration of Chloride-induced Stress Corrosion Cracking in Spent Nuclear Fuel Storage Canisters

This report describes the body of work involving testing and analyses performed in FY18 and FY19 at the Savannah River National Laboratory and at Korea University to determine chloride-induced stress corrosion crack (CISCC) growth rates in weldments of austenitic stainless steel plate material. The primary work involved: several separate crack growth rate test campaigns with new testing setups using laboratory-scale corrosion specimens; analyses to estimate the stress intensity factors for flaws postulated to be present in the residual stress fields of axial and circumferential weldments in a prototypic dry storage canister; and establishment of a “large plate” (approximately 51×46 cm with thickness 16 mm and weight about 30 kg) CISCC test that uses a plate with a circumferential weldment harvested from the Sandia National Laboratories’ full-size mockup dry storage canister. The tests for crack growth rate were conducted on ASTM E1861 bolt-load compact tension specimens in a setup designed to allow initially dried salt deposits to deliquesce and infuse the brine to the crack front under conditions relevant to the canister storage environments (e.g., temperature and humidity). The crack growth rate test results matched previous experimental data in the open literature. These literature data had been compiled by SNL and EPRI to develop a temperature-dependent (K-insensitive) crack growth rate model to support flaw disposition in the ASME Boiler and Pressure Vessel Section XI Code Case N-860, “Examination Requirements and Acceptance Standards for Spent Nuclear Fuel Storage and Transportation Containment Systems.” The knowledge of the salt, temperature, humidity conditions for CISCC including the threshold stress intensity factor, $K_{ISCC}$, were used to setup the large plate test in which a set of partthrough- wall and through-wall electrical discharge machined (EDM) defects were positioned in the weld residual stress field of the welded plate material from the mockup canister. A salt solution was applied to this plate specimen with an air-brush spray, then dried, and then exposed to a constant 75% RH condition at an ambient room temperature of 22 °C. The salt loading and environmental condition is expected to be aggressive to cause CISCC. The visual results at the surface of the plate show corrosion product staining (brown rust appearance) but with no CISCC emanating from the machined notches after approximately 5 months exposure to date. The plate test will continue in FY20 and phased-array ultrasonic testing and other NDE techniques will be applied to identify and characterize CISCC that may emanate from the machined defects. The notches and any CISCC from them will be sectioned at the conclusion of the test to verify NDE characterization, and to develop effective through-wall crack growth rates from cracks that would likely be branched.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Systematic and objective evaluation of Earth system models: PCMDI Metrics Package (PMP) version 3

Systematic, routine, and comprehensive evaluation of Earth system models (ESMs) facilitates benchmarking improvement across model generations and identifying the strengths and weaknesses of different model configurations. By gauging the consistency between models and observations, this endeavor is becoming increasingly necessary to objectively synthesize the thousands of simulations contributed to the Coupled Model Intercomparison Project (CMIP) to date. The Program for Climate Model Diagnosis and Intercomparison (PCMDI) Metrics Package (PMP) is an open-source Python software package that provides quick-look objective comparisons of ESMs with one another and with observations. The comparisons include metrics of large- to global-scale climatologies, tropical inter-annual and intra-seasonal variability modes such as the El Niño–Southern Oscillation (ENSO) and Madden–Julian Oscillation (MJO), extratropical modes of variability, regional monsoons, cloud radiative feedbacks, and high-frequency characteristics of simulated precipitation, including its extremes. The PMP comparison results are produced using all model simulations contributed to CMIP6 and earlier CMIP phases. An important objective of the PMP is to document the performance of ESMs participating in the recent phases of CMIP, together with providing version-controlled information for all datasets, software packages, and analysis codes being used in the evaluation process. Among other purposes, this also enables modeling groups to assess performance changes during the ESM development cycle in the context of the error distribution of the multi-model ensemble. Quantitative model evaluation provided by the PMP can assist modelers in their development priorities. In this paper, we provide an overview of the PMP, including its latest capabilities, and discuss its future direction.

54 ENVIRONMENTAL SCIENCES↗

A Digital Control Algorithm for Magnetic Suspension Systems

An ongoing program exists to investigate and develop magnetic suspension technologies and modelling techniques at NASA Langley Research Center. Presently, there is a laboratory-scale large air-gap suspension system capable of five degree-of-freedom (DOF) control that is operational and a six DOF system that is under development. Those systems levitate a cylindrical element containing a permanent magnet core above a planar array of electromagnets, which are used for levitation and control purposes. In order to evaluate various control approaches with those systems, the Generic Real-Time State-Space Controller (GRTSSC) software package was developed. That control software package allows the user to implement multiple control methods and allows for varied input/output commands. The development of the control algorithm is presented. The desired functionality of the software is discussed, including the ability to inject noise on sensor inputs and/or actuator outputs. Various limitations, common issues, and trade-offs are discussed including data format precision; the drawbacks of using either Direct Memory Access (DMA), interrupts, or program control techniques for data acquisition; and platform dependent concerns related to the portability of the software, such as memory addressing formats. Efforts to minimize overall controller loop-rate and a comparison of achievable controller sample rates are discussed. The implementation of a modular code structure is presented. The format for the controller input data file and the noise information file is presented. Controller input vector information is available for post-processing by mathematical analysis software such as MATLAB1.

Britton, Thomas C.↗

Grid Sensitivity Study for Slat Noise Simulations

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations in conjunction with a Ffowcs Williams-Hawkings acoustics solver. Many previous simulations have been performed for the configuration, and the case was introduced as a new category for the Second AIAA workshop on Benchmark problems for Airframe Noise Configurations (BANC-II). However, the cost of the simulations has restricted the study of grid resolution effects to a baseline grid and coarser meshes. In the present study, two different approaches are being used to investigate the effect of finer resolution of near-field unsteady structures. First, a standard grid refinement by a factor of two is used, and the calculations are performed by using the same CFL3D solver employed in the majority of the previous simulations. Second, the OVERFLOW code is applied to the baseline grid, but with a 5th-order upwind spatial discretization as compared with the second-order discretization used in the CFL3D simulations. In general, the fine grid CFL3D simulation and OVERFLOW calculation are in very good agreement and exhibit the lowest levels of both surface pressure fluctuations and radiated noise. Although the smaller scales resolved by these simulations increase the velocity fluctuation levels, they appear to mitigate the influence of the larger scales on the surface pressure. These new simulations are used to investigate the influence of the grid on unsteady high-lift simulations and to gain a better understanding of the physics responsible for the noise generation and radiation.

Lockard, David P.↗

Sensitivity Study of Multiscale and Phenomenological Elasto-Viscoplastic Grade 91 Material Models for Component-Scale Response

Many advanced nuclear reactor concepts currently being developed are targeting higher operating temperatures relative to the current fleet of light water nuclear reactors, for efficiency gains and other operational considerations. The design of high temperature structural components with reliable long-term operational performance will depend on material models that accurately capture the inelastic deformation mechanisms active in these environments. In this work, we perform a detailed parameter sensitivity analysis of two unified elasto-viscoplastic Grade 91 material models capable of capturing long term high temperature creep deformation. The first model is a phenomelogical material model from the Nuclear Engineering Material Library (NEML) developed at Argonne National Lab. The NEML model parameters and their uncertainty were fit to a range of Grade 91 experimental data using Bayesian Markov Chain Monte Carlo analysis. The second model is a LAROMance data-driven surrogate material model developed at Los Alamos National Lab. The LAROMance model is fit to a large database of responses produced by a mechanistic crystal plasticity based polycrystal model. Parameters for the LAROMance surrogate material model reflect the pedigree of the Grade 91 microstructure. Both material models have been integrated into the Grizzly code, based on the open-source MOOSE multiphysics simulation framework, to simulate both the progression of aging mechanisms and the effects of that aging on nuclear power plant structures. Grizzly is used analyze a three-dimensional Grade 91 piping system to compare the long-term inelastic response predicted by these two fundamentally different models and assess the sensitivity of the material model input parameters on this quantity of interest.

42 ENGINEERING↗

Electric Propulsion Interactions Code (EPIC): Recent Enhancements and Goals for Future Capabilities

The Electric Propulsion Interactions Code (EPIC) is the leading interactive computer tool for assessing the effects of electric thruster plumes on spacecraft subsystems. EPIC, developed by SAIC under the sponsorship of the Space Environments and Effects (SEE) Program at the NASA Marshall Space Flight Center, has three primary modules. One is PlumeTool, which calculates plumes of electrostatic thrusters and Hall-effect thrusters by modeling the primary ion beam as well as elastic scattering and charge-exchange of beam ions with thruster-generated neutrals. ObjectToolkit is a 3-D object definition and spacecraft surface modeling tool developed for use with several SEE Program codes. The main EPIC interface integrates the thruster plume into the 3-D geometry of the spacecraft and calculates interactions and effects of the plume with the spacecraft. Effects modeled include erosion of surfaces due to sputtering, re-deposition of sputtered materials, surface heating, torque on the spacecraft, and changes in surface properties due to erosion and deposition. In support of Prometheus I (JIMO), a number of new capabilities and enhancements were made to existing EPIC models. Enhancements to EPIC include adding the ability to scale and view individual plume components, to import a neutral plume associated with a thruster (to model a grid erosion plume, for example), and to calculate the plume from new initial beam conditions. Unfortunately, changes in program direction have left a number of desired enhancements undone. Variable gridding over a surface and resputtering of deposited materials, including multiple bounces and sticking coefficients, would significantly enhance the erosion/deposition model. Other modifications such as improving the heating model and the PlumeTool neutral plume model, enabling time dependent surface interactions, and including EM1 and optical effects would enable EPIC to better serve the aerospace engineer and electric propulsion systems integrator. We review EPIC S overall capabilities and recent modifications, and discuss directions for future enhancements.

Gardner, Barbara M.↗

Libpanda: A High Performance Library for Vehicle Data Collection

Cyber-Physical Systems (CPS) generally involve time-critical components due to physical dynamics, therefore necessitating high-performance subsystems. This is also true in data collection scenarios to infer physical phenomena. This paper covers Libpanda as an example of a component that has been designed to address performance issues in CPS implementations. Libpanda is a C++ library that interfaces software with a Comma.ai Panda device. Pandas are used for installation in modern vehicles to read the vehicle CAN bus, providing rich sensor data and limited vehicle control through message injection. The motivation to design lib-panda stems from the lack of performance in Python-based code that runs on inexpensive hardware like a Raspberry Pi. In such situations, Python code would result in utilizing 92% CPU while also dropping around 40% of the CAN packet due to bottlenecks. Without using different tools, inconsistent data collection means a loss of time-based vehicle state interpretation. Libpanda addresses these issues through implementation in a different language and implementation of different design paradigms involving asynchronous calls and multithreading. The Panda also features a GPS module that allows multiple instances to synchronize clocks for large-scale data collection scenarios. Libpanda has been designed with time-synchronization in mind to aid in the measurement of inter-vehicle dynamics. The performance improvements of libpanda have resulted in it becoming an important component in automotive dynamics research that requires a higher technical performance in large-scale experiments.

Bunting, Matthew↗

A Safe Response to Renewable Energy Hazards

The International Association of Fire Fighters (IAFF) and Underwriters Laboratories, LLC (UL) in conjunction with UL Solutions initiated a joint project in 2022 under an agreement with the United States Department of Energy-Office of Energy Efficiency and Renewable Energy (DOE-EERE). This project focused on two separate and important initiatives related to energy efficiency in residential buildings. Initiative 1: Fire Performance on Energy Efficient Exterior Walls. Initiative 2: Firefighting Tactics in Residential Properties with Building Energy Storage Systems (BESS). The project’s first initiative addresses concerns surrounding new technologies with enhanced, energy-efficient exterior walls installed on residential properties. The concerns of fire rapidly traveling vertically up the exterior of these walls were examined. This addressed a growing concern from the first responder community that many times, the fires on the exterior of residential buildings have already evolved into an attic fire by the time of arrival – making it problematic to address the fire scenario. The test plan for Initiative 1 incorporated a modified version of an American Society for Testing and Materials (ASTM) test method, ASTM E2707, Standard Test Method for Determining Fire Penetration of Exterior Wall Assemblies Using a Direct Flame Impingement Exposure, as the foundation of the research. The test method involved a wall structure intended to represent retrofit construction to evaluate how fire would spread vertically or laterally. The second aspect of the UL-IAFF Project focuses on the fire service response to Residential Battery Energy Storage System (RBESS) incidents. These simulation tests were constructed in the large-scale fire test facility at UL Solutions’ Northbrook, IL campus. A baseline test was conducted that involved a test structure with no batteries—shelving units populated with standardized commodities, representing a typical garage with cellulosic and plastic contents. Three additional tests have been conducted to generate data with the contribution of energy storage system (ESS) batteries to compare fire and explosion hazards against the baseline test. Through this work, fire service tactical considerations can be explored. From the data, the team can determine 1) the visual indicators of a residential fire that has involved an RBESS (or, potentially, other large batteries) and 2) the impact of fire service-initiated ventilation of the structure on the fire conditions and explosion risks.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Conservation equations and physical models for hypersonic air flows over the aeroassist flight experiment vehicle

The code development and application program for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA), with emphasis directed toward support of the Aeroassist Flight Experiment (AFE) in the near term and Aeroassisted Space Transfer Vehicle (ASTV) design in the long term is reviewed. LAURA is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3-D, viscous, hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over an AFE wind tunnel model. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium over the full scale AFE configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗

ParFlow Sand Tank: A tool for groundwater exploration

The ParFlow Sand Tank model is an open source application designed to allow users to interactively simulate and visualize groundwater movement through the subsurface. The app is designed for both research and education; teaching hydrogeology concepts and making it easy explore and run sophisticated groundwater simulations. Our goal is to support increased accessibility and usability of research grade hydrology tools for research and teaching. The Sand Tank application simulates groundwater and surface water fluxes as well as contaminant transport in real time using the integrated physical hydrology model ParFlow (Kollet & Maxwell, 2006; Maxwell & Miller, 2005; Osei-Kuffuor et al., 2014) and the particle tracking code EcoSlim (Maxwell et al., 2019). ParFlow is a numerical hydrology model that simulates spatially distributed groundwater and surface water flow. It is a well established research tool with more than 90 publications documenting its development use to advance our understanding of groundwater dynamics and groundwater surface water interactions from the hillslope to the continental scale e.g. (Condon et al., 2020; Condon & Maxwell, 2019; Maxwell & Condon, 2016). It is designed for efficient parallel computation and has been run on many platforms spanning from laptops to supercomputers. However, one of the challenges of ParFlow is that it requires significant training and hydrologic expertise to develop simulations. The Sand Tank application makes this model accessible to anyone for education and exploration. Our application uses ParFlow for its simulation backend and ParaView for the data loading and processing. The communication infrastructure relies on the ParaViewWeb framework. We use model templates deployed in Docker images to setup the Sand Tank framework. Users can build the application locally or interact with it through our web deployment. When interacting with a template users can interactively change model parameters like subsurface processes or pump/inject water into the subsurface and watch the system respond to their changes in real time as the simulation runs. Additionally, our template setup will allow more advanced users to build custom templates of increasing complexity for both research and educational purposes.

54 ENVIRONMENTAL SCIENCES↗

Modeling and Simulation of Xe-100-type Pebble Bed Gas-Cooled Reactor with SCALE

The US Department of Energy (DOE) announced the Advanced Reactor Demonstration Program (ARDP) to accelerate the deployment of advanced reactor concepts. Awardees of ARDP funds are expected to demonstrate the operation of an advanced reactor within 7 years of receiving the award. X-Energy’s advanced reactor concept, the Xe-100, was selected as one of two advanced reactor concepts to receive funding to demonstrate the operation of its high-temperature gas-cooled pebble-bed reactor before the end of this decade. As a result of this push to bring advanced reactors to maturation and commercialization, transition and deployment scenario studies are being performed under the Systems Analysis and Integration (SA&I) campaign within the DOE Office of Nuclear Energy (DOE-NE) to evaluate the transition of the current US commercial fleet of light-water reactors (LWRs) to a future fleet of advanced reactors consisting of a mix of ARDP type reactor concepts and advanced LWRs. To accurately evaluate the front- and back-end resource requirements, it is important to perform reactor physics calculations to determine the discharge burnup and isotopic content, fuel residence time, as well as other parameters. For this purpose, a summer project funded by the SA&I campaign allowed for the setup of SCALE models for full-core Xe-100 type high-temperature gas-cooled pebble-bed reactor and a Xe-100 type slice using publicly available information. The core-averaged equilibrium compositions and zone-wise equilibrium compositions for the slice and 3D models, respectively, were obtained following an iterative depletion method developed by Bostelmann et al. using SCALE’s reactor physics sequence TRITON. The slice model was used with TRITON to generate burnup-dependent cross section libraries at different temperatures which can be used with SCALE’s ORIGAMI code to rapidly determine fuel inventory and therefore to perform quick sensitivity studies on parameters such as the pebble location in the core. The SCALE/TRITON transport and depletion calculation for the Xe-100 type slice model indicates that the isotopic concentrations are in good agreement at 1,300 effective full power days (EFPD) for 235 U. An analysis of 236 U results match 239 Pu results would seem to indicate a typographical error in Mulder and Boyes wherein the reported results of 236 U and 239 Pu are reversed. In addition to SCALE/TRITON calculations, a new capability within SCALE/ORIGAMI for the simulation of pebble-bed reactors was used to study the burnup sensitivity with respect to the pebble pathway through the core. The SCALE/ORIGAMI results show that pebbles that travel closer to the reflector for the entire depletion history have a higher burnup than pebbles that travel through the middle of the core because of the higher thermal to fast flux ratio near the reflector. Consequently, a pebble’s burnup is strongly affected by the pebble’s pathway for each pass. Additional phenomena such as temperature distributions in the core and different travel times of the pebbles in the individual radial zones further affect the burnup distribution. The sensitivity of the discharge vector to the pebble pathways taken during each pass can be evaluated in the future using SCALE/ORIGAMI now that the SCALE inputs have been established.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CMU DeepLens: deep learning for automatic image-based galaxy–galaxy strong lens finding

Galaxy-scale strong gravitational lensing can not only provide a valuable probe of the dark matter distribution of massive galaxies, but also provide valuable cosmological constraints, either by studying the population of strong lenses or by measuring time delays in lensed quasars. Due to the rarity of galaxy-scale strongly lensed systems, fast and reliable automated lens finding methods will be essential in the era of large surveys such as Large Synoptic Survey Telescope, Euclid and Wide-Field Infrared Survey Telescope. To tackle this challenge, we introduce CMU DeepLens, a new fully automated galaxy-galaxy lens finding method based on deep learning. This supervised machine learning approach does not require any tuning after the training step which only requires realistic image simulations of strongly lensed systems. We train and validate our model on a set of 20 000 LSST-like mock observations including a range of lensed systems of various sizes and signal-to-noise ratios (S/N). We find on our simulated data set that for a rejection rate of non-lenses of 99 per cent, a completeness of 90 per cent can be achieved for lenses with Einstein radii larger than 1.4 arcsec and S/N larger than 20 on individual g-band LSST exposures. Finally, we emphasize the importance of realistically complex simulations for training such machine learning methods by demonstrating that the performance of models of significantly different complexities cannot be distinguished on simpler simulations. We make our code publicly available at https://github.com/McWilliamsCenter/CMUDeepLens.

gravitational lensing: strong↗