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,495 records · Page 83

MODELING THE AMBIENT CONDITION EFFECTS OF AN AIR-COOLED NATURAL CIRCULATION SYSTEM

The Reactor Cavity Cooling System (RCCS) is a passive safety concept under consideration for the overall safety strategy of advanced reactors such as the High Temperature Gas-Cooled Reactor (HTGR). One such variant, air-cooled RCCS, uses natural convection to drive the flow of air from outside the reactor building to remove decay heat during normal operation and accident scenarios. The Natural convection Shutdown heat removal Test Facility (NSTF) at Argonne National Laboratory (“Argonne”) is a half-scale model of the primary features of one conceptual air-cooled RCCS design. The facility was constructed to carry out highly instrumented experiments to study the performance of the RCCS concept for reactor decay heat removal that relies on natural convection cooling. Parallel modeling and simulation efforts were performed to support the design, operation, and analysis of the natural convection system. Throughout the testing program, strong influences of ambient conditions were observed in the experimental data when baseline tests were repeated under the same test procedures. Thus, significant analysis efforts were devoted to gaining a better understanding of these influences and the subsequent response of the NSTF to ambient conditions. It was determined that air humidity had negligible impacts on NSTF system performance and therefore did not warrant consideration in the models. However, temperature differences between the building exterior and interior air, along with the outside wind speed, were shown to be dominant factors. Combining the stack and wind effects together, an empirical model was developed based on theoretical considerations and using experimental data to correlate zero-power system flow rates with ambient meteorological conditions. Some coefficients in the model were obtained based on best fitting the experimental data. The predictive capability of the empirical model was demonstrated by applying it to the new set of experimental data. The empirical model was also implemented in the computational models of the NSTF using both RELAP5-3D and STARCCM+ codes. Accounting for the effects of ambient conditions, simulations from both codes predicted the natural circulation flow rates very well.

Hu, Rui↗

Unbalancing Symbiotic Nitrogen Fixation: Can We Make Effectiveness More Effective?

One of the most critical uses of carbon fuels is in generating nitrogen fertilizer. Atmospheric dinitrogen is too stable to be used directly by plants, so plants need other chemical forms of nitrogen. Nitrogen fertilizer production uses ~4% of the world’s natural gas. Making N fertilizers leverages methane’s energy content by enabling additional energy to be captured into food and fiber through increased photosynthesis. Manufacturing fertilizer is essential—without it there would only be enough combined nitrogen to feed ~3 billion of the world’s ~7 billion people. How to alter this addiction is not obvious, but the best chance to secure substantial and sustainable amounts of future N might be through symbiotic nitrogen fixation (SNF). SNF describes mutualistic interactions where nitrogen-fixing bacteria provide fixed nitrogen to their plant hosts, freeing them from the need for nitrogen fertilizers. SNF already occurs on a large scale; "biological nitrogen” contributes about 5 Tg N/yr to American agriculture, a quarter of the nitrogen needed. Improving SNF and replacing inorganic fertilizers are both important goals in establishing sustainable and more energy-efficient farming. SNF has been used in agriculture for millennia, largely by using legumes like beans or alfalfa in mixed cropping systems. Legume SNF occurs in root nodules, organs that develop after a growing root is infected by bacteria called rhizobia. In the plant cells of a legume nodule, specialized forms of rhizobia fix nitrogen and can produce enough ammonia to supply the growing plant. The exchange of photosynthetically derived plant carbon compounds for nitrogen reduced by the bacteria is the engine that drives the symbiosis. Important details about the organization of development and metabolism in SNF remain to be determined, such as how plant and bacterial metabolism are coordinated as nitrogen is being fixed and what governs the overall level of nitrogen fixation. Most bacterial mutations that affect SNF completely block fixation (e.g. are Fix – ), providing a limited window into the feedback interactions that must be part of the process. We have investigated an unusual mutant of Sinorhizobium meliloti, a symbiotic partner of alfalfa, that fixes dinitrogen at a normal rate (i.e., it is Fix + ) but is not effective in supporting plant growth (i.e., it is ineffective or Eff – ). We have shown that the mutant has a specific mutation in the glnD gene, which codes for a protein that regulates the bacterial response to nitrogen stress, among other key pathways. A Fix + Eff – phenotype contains an unavoidable puzzle—how is it possible for the bacteria to fix nitrogen at a normal rate and without benefiting the nitrogen starved plant host? After eliminating other possibilities, we proposed that the bacteria synthesize a nitrogen-containing compound that the plant can’t use so that acquiring this compound does not relieve plant nitrogen stress. Extracts of nodules made by glnD mutant strains have high concentrations of pyruvate canaline oxime (PCO), a derivative of the plant defensive compound canaline, a toxic analog of ornithine, an amino acid. Many legumes produce canaline and the related arginine mimic canavanine. For reasons we do not yet understand, the mutation in glnD appears to trigger significant overproduction of canaline in a defense response that leads to substantial synthesis of PCO. We expected the bacteria to be able to use PCO, but we have not been able to grow S. meliloti on PCO under several conditions. PCO thus appears to be a metabolic dead end for both the plant and bacteria. We also examined bacterial catabolic pathways that degrade arginine and canavanine and showed that mutation of these did not alter the symbiosis in an obvious way. In related experiments to examine the metabolic development of root nodules, we carried out an ambitious experiment to simultaneously measure metabolites and protein changes during nodule maturation. For these experiments we used a related symbiotic interaction between S. medicae and Medicago truncatula. M. truncatula has a simpler genome than alfalfa so we could identify plant proteins based on predictions from the DNA sequence. Nodules formed on Medicago grow linearly, with the least mature tissues at the tip. Because development of nitrogen fixation is accompanied by the production of the pigmented plant oxygen carrier protein leghemoglobin, we were able to locate and manually separate the tip, transitional, and mature nitrogen-fixing regions from each other. Advances in proteomic and metabolomic techniques allowed us to analyze these tissues from a pool of as few as 10 nodules. We observed correlated transitions of various enzymes in the plant and bacterial proteomes from those characteristic of free-living metabolism to the microaerobic metabolism of the nitrogen-fixing tissues. The metabolome could not be separated into plant and bacterial domains but was consistent with this overall transition. A subsequent paper examined proteolysis in nodule bacteria. Protein turnover in mature regions of the nodule is significant. We were able to show complexes between various proteins and important proteases, using precipitation capture techniques to isolate the proteases and sensitive proteomic analysis to reveal the associated proteins.

09 BIOMASS FUELS↗

Design and implementation of wireless dose logger network for radiological emergency decision support system

A decision support system (DSS) is implemented in Radiological Safety Division, Indira Gandhi Centre for Atomic Research for providing guidance for emergency decision making in case of an inadvertent nuclear accident. Real time gamma dose rate measurement around the stack is used for estimating the radioactive release rate (source term) by using inverse calculation. Wireless gamma dose logging network is designed, implemented, and installed around the Madras Atomic Power Station reactor stack to continuously acquire the environmental gamma dose rate and the details are presented in the paper. The network uses XBee–Pro wireless modules and PSoC controller for wireless interfacing, and the data are logged at the base station. A LabView based program is developed to receive the data, display it on the Google Map, plot the data over the time scale, and register the data in a file to share with DSS software. The DSS at the base station evaluates the real time source term to assess radiation impact.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Scalable line and plane relaxation in a parallel structured multigrid solver

The efficient solution of sparse, linear systems that arise through the discretization of partial differential equations remains a key challenge for a range of high performance scientific simulations. One approach for reducing data movement and improving performance is by exposing and exploiting structure in a problem through the use of robust structured multilevel solvers. By choosing coarsening that preserves the structure of the problem, these methods maintain efficient structured computation and communication throughout the multigrid hierarchy. However, when coarsening is not permitted to be dependent on the operator, anisotropy must be addressed by the smoother — producing error compatible for coarse-grid correction with structured coarsening. Here, the components required in a scalable parallel structured solver are described with a focus on memory and communication efficiency of robust smoothers. While the implementation of communication and memory reduction techniques in smoothers integrated in a complete 3D solver present a significant engineering challenge, a novel approach is proposed that addresses these challenges systematically through a change to the solver’s execution model. Enabled by user-level threading paired with a set of data and communication abstractions, this approach permits seamless aggregation of communication in plane smoothers — directly reusing code for a 2D distributed multilevel cycle. Results show an effective reduction in communication costs for coarse-grid problems, and result in a speedup of 8.7x in smoothing routines shown in Fig. 12 using this approach. This produces a significant improvement to strong scalability while maintaining favorable weak scaling behavior. Finally, a parallel scaling study using a series of refined meshes is included that demonstrates the effectiveness of this approach in an application of interest.

97 MATHEMATICS AND COMPUTING↗

Sci-Sat AM: Radiation Dosimetry and Practical Therapy Solutions - 03: Energy dependence of a clinical probe-format calorimeter and its pertinence to absolute photon and electron beam dosimetry

Purpose: To evaluate the intrinsic and absorbed-dose energy dependence of a small-scale graphite calorimeter probe (GPC) developed for use as a routine clinical dosimeter. The influence of charge deposition on the response of the GPC was also assessed by performing absolute dosimetry in clinical linac-based electron beams. Methods: Intrinsic energy dependence was determined by performing constant-temperature calorimetry dose measurements in a water-equivalent solid phantom, under otherwise reference conditions, in five high-energy photon (63.5 < %dd(10){sub X} < 76.3), and five electron (2.3 cm < R{sub 50} < 8.3 cm) beams. Reference dosimetry was performed for all beams in question using an Exradin A19 ion chamber with a calibration traceable to national standards. The absorbed-dose component of the overall energy dependence was calculated using the EGSnrc egs-chamber user code. Results: A total of 72 measurements were performed with the GPC, resulting in a standard error on the mean absorbed dose of better than 0.3 % for all ten beams. For both the photon and electron beams, no statistically-significant energy dependence was observed experimentally. Peak-to-peak, variations in the relative response of the GPC across all beam qualities of a given radiation type were on the order of 1 %. No effects, either transient or permanent, were attributable to the charge deposited by the electron beams. Conclusions: The GPC’s apparent energy-independence, combined with its well-established linearity and dose rate independence, make it a potentially useful dosimetry system capable measuring photon and electron doses in absolute terms at the clinical level.

61 RADIATION PROTECTION AND DOSIMETRY↗

A Fault-oblivious Extreme-scale Execution Environment

Exascale computing systems will provide a thousand-fold increase in parallelism and a proportional increase in failure rate relative to today's machines. Systems software for exascale machines must provide the infrastructure to support existing applications while simultaneously enabling efficient execution of new programming models that naturally express dynamic, adaptive, irregular computation; coupled simulations; and massive data analysis in a highly unreliable hardware environment with billions of threads of execution. We propose a new approach to the data and work distribution model provided by system software based on the unifying formalism of an abstract file system. The proposed hierarchical data model provides simple, familiar visibility and access to data structures through the file system hierarchy, while providing fault tolerance through selective redundancy. The hierarchical task model features work queues whose form and organization are represented as file system objects. Data and work are both first class entities. By exposing the relationships between data and work to the runtime system, information is available to optimize execution time and provide fault tolerance. The data distribution scheme provides replication (where desirable and possible) for fault tolerance and efficiency, and it is hierarchical to make it possible to take advantage of locality. The user, tools, and applications, including legacy applications, can interface with the data, work queues, and one another through the abstract file model. This runtime environment will provide multiple interfaces to support traditional Message Passing Interface applications, languages developed under DARPA's High Productivity Computing Systems program, as well as other, experimental programming models. We will validate our runtime system with pilot codes on existing platforms and will use simulation to validate for exascale-class platforms. In this final report, we summarize research results from the work done at the Ohio State University towards the larger goals of the project listed above.

97 MATHEMATICS AND COMPUTING↗

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗

Preliminary Logistics Modeling of the Integrated Wash Aid Treatment Emergency Reuse System - 19495

We employed Analysis of Mobility Platform (AMP) and GoldSim software to model the logistics of deploying materials to implement the Integrated Wash Aid Treatment Emergency Reuse System (IWATERS) in a hypothetical cesium-contamination scenario covering a four-block strip of downtown Chicago, IL. In this deployment of IWATERS, buildings were washed down via fire-hosing and contaminated wash solution was collected in the city's waste/stormwater collection system. The contaminated wash water was treated to remove Cs via ad hoc sand/clay filtration beds, and the treated wash solution was recycled during continuing decontamination operations. A detailed task timeline for the IWATERS deployment and information (or assumptions) on the availability and location of required resources for operations provided input parameters, as well as a framework, to enable the AMP and GoldSim simulations. From the decision to deploy IWATERS to the end of operations, including the completion of processing all contaminated wash water, took 49 days with decontamination teams active for 35 days and treatment bed teams active for 37 days. The modeling results showed the availability of vermiculite clay (based on data supplied by its commercial vendors/distributors) as the limiting factor for deploying treatment bed teams and, consequentially, the recycling rate of wash water and the timeline for decontamination operations. These simulations represent the first attempt at understanding the material and personnel requirements for large-scale IWATERS operation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Context-aware learning of hierarchies of low-fidelity models for multi-fidelity uncertainty quantification

Multi-fidelity Monte Carlo methods leverage low-fidelity and surrogate models for variance reduction to make tractable uncertainty quantification even when numerically simulating the physical systems of interest with high-fidelity models is computationally expensive. This work proposes a context-aware multi-fidelity Monte Carlo method that optimally balances the costs of training low-fidelity models with the costs of Monte Carlo sampling. It generalizes the previously developed context-aware bi-fidelity Monte Carlo method to hierarchies of multiple models and to more general types of low-fidelity models. When training low-fidelity models, the proposed approach takes into account the context in which the learned low-fidelity models will be used, namely for variance reduction in Monte Carlo estimation, which allows it to find optimal trade-offs between training and sampling to minimize upper bounds of the mean-squared errors of the estimators for given computational budgets. This is in stark contrast to traditional surrogate modeling and model reduction techniques that construct low-fidelity models with the primary goal of approximating well the high-fidelity model outputs and typically ignore the context in which the learned models will be used in upstream tasks. Further, the proposed context-aware multi-fidelity Monte Carlo method applies to hierarchies of a wide range of types of low-fidelity models such as sparse-grid and deep-network models. Numerical experiments with the gyrokinetic simulation code Gene show speedups of up to two orders of magnitude compared to standard estimators when quantifying uncertainties in small-scale fluctuations in confined plasma in fusion reactors. This corresponds to a runtime reduction from 72 days to four hours on one node of the Lonestar6 supercomputer at the Texas Advanced Computing Center.

42 ENGINEERING↗

A fault-tolerant neutral-atom architecture for universal quantum computation

Quantum error correction (QEC) is essential for the realization of large-scale quantum computers. However, owing to the complexity of operating on the encoded ‘logical’ qubits, understanding the physical principles for building fault-tolerant quantum devices and combining them into efficient architectures is an outstanding scientific challenge. Here we use reconfigurable arrays of up to 448 neutral atoms to implement the key elements of a universal, fault-tolerant quantum processing architecture and experimentally explore their underlying working mechanisms. We first use surface codes to study how repeated QEC suppresses errors, demonstrating 2.14(13)x below-threshold performance in a four-round characterization circuit by leveraging atom loss detection and machine learning decoding. We then investigate logical entanglement using transversal gates and lattice surgery and extend it to universal logic through transversal teleportation with three-dimensional [[15,1,3]] codes, enabling arbitrary-angle synthesis with polylogarithmic overhead. Finally, we develop mid-circuit qubit reuse16, increasing experimental cycle rates by two orders of magnitude and enabling deep-circuit protocols with dozens of logical qubits and hundreds of logical teleportations with [[7,1,3]] and high-rate [[16,6,4]] codes while maintaining constant internal entropy. Our experiments show key principles for efficient architecture design, involving the interplay between quantum logic and entropy removal, judiciously using physical entanglement in logic gates and magic state generation, and leveraging teleportations for universality and physical qubit reset. These results establish foundations for scalable, universal error-corrected processing and its practical implementation in neutral atom systems.

atomic and molecular physics↗

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders

The SCALE code package offers multiple nuclear data libraries and sensitivity and uncertainty (S/U) methods supporting and derived from Monte Carlo (MC) transport. The CSAS and TSUNAMI-3D sequences use KENO MC, utilizing either continuous-energy (CE) cross sections or multigroup (MG) cross section libraries. TSUNAMI-3D has two CE calculational methods: the iterated fission probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications (e.g., fissionable and polyethylene reflectors). The HEU-MET-FAST-084 (HMF-084) International Criticality Safety Benchmark Evaluation Project evaluation consists of 27 cylindrical highly enriched uranium metal cores with 14 unique reflector materials of 0.5 and 1 in. thicknesses. Included in this list of reflector materials are natural uranium and polyethylene. This work utilized SCALE 6.2.4 models of the HMF-084 evaluation, with additional non-physical configurations to test both the MG bias and CLUTCH functionality across a variety of reflector material thicknesses. The evaluation’s use of concentric cylinders allowed for examination of several MG self-shielding methods: infinite homogenous, cylindrical, and spherical. The results indicate that the use of polyethylene reflectors with CLUTCH is not fundamentally impossible but sensitive to geometry. The poor performance of CLUTCH with fissionable reflectors was reaffirmed. The 2 in. and greater polyethylene-reflected calculations demonstrate the necessity of using the 302-group library for fast systems. The nickel MG bias was substantial, as discussed in a companion paper, as were cobalt and iron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupled Experimental/Computational Investigation of the Dynamics of Interacting Magnetized Plasmas

The interaction, or interpenetration, of magnetized plasmas of different density and/or pressure occurs in a wide variety of natural and man-made systems. Such systems include extragalactic jets propagating into the intergalactic medium, solar coronal mass ejections into background solar wind, compact toroid (CT) fueling of magnetic fusion plasmas, and jets of capsule shell impurities into DT fusion fuel, which can lead to enhanced impurity mix in inertial fusion implosions. These plasmas may take the form of jets, with open, helical magnetic structure, or plasma “bubbles” with closed magnetic fields (B-fields), such as spheromaks or CT’s. Such structures, both open and closed B-field cases, can transport heat, particles and magnetic flux or magnetic helicity into background plasma regions. For example, the origin of extragalactic magnetic fields may be due, at least in part, to transport by astrophysical jets. The goal of this proposed work was to elucidate the detailed plasma and magnetic field dynamics of high-density plasma jets (open B-field) and bubbles (closed B-field) propagating into lower density background magnetized plasma through controlled laboratory experiments and closely coupled nonlinear MHD modeling. These experiments were conducted in the HelCat (Helicon-Cathode) linear plasma device at the University of New Mexico (UNM). Plasma jets and bubbles were launched via an existing compact coaxial plasma gun, mounted on the HelCat device. This gun produced plasmas tens of cm in scale and lasting tens of microseconds, thereby allowing detailed multipoint, space- and time-resolved measurements to be made routinely. The experiments were directly modeled using the extended magnetohydrodynamic (XMHD) PERSEUS code, developed at Cornell University [23,24]. Both experimental and numerical modeling work are ongoing. The main results to date are reported here. Additional supplemental funding for one year (8/1/2019 – 7/31/2020) supported numerical investigation of photoionization processes important in many low temperature plasmas, including the HelCat device. Initial results of this modeling work is also reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates

Progressive damage and failure analysis (PDFA) tools are needed to predict the nonlinear response of advanced fiber-reinforced composite structures. Predictive tools should incorporate the underlying physics of the damage and failure mechanisms observed in the composite, and should utilize as few input parameters as possible. The purpose of the Enhanced Schapery Theory (EST) was to create a PDFA tool that operates in conjunction with a commercially available finite element (FE) code (Abaqus). The tool captures the physics of the damage and failure mechanisms that result in the nonlinear behavior of the material, and the failure methodology employed yields numerical results that are relatively insensitive to changes in the FE mesh. The EST code is written in Fortran and compiled into a static library that is linked to Abaqus. A Fortran Abaqus UMAT material subroutine is used to facilitate the communication between Abaqus and EST. A clear distinction between damage and failure is imposed. Damage mechanisms result in pre-peak nonlinearity in the stress strain curve. Four internal state variables (ISVs) are utilized to control the damage and failure degradation. All damage is said to result from matrix microdamage, and a single ISV marks the micro-damage evolution as it is used to degrade the transverse and shear moduli of the lamina using a set of experimentally obtainable matrix microdamage functions. Three separate failure ISVs are used to incorporate failure due to fiber breakage, mode I matrix cracking, and mode II matrix cracking. Failure initiation is determined using a failure criterion, and the evolution of these ISVs is controlled by a set of traction-separation laws. The traction separation laws are postulated such that the area under the curves is equal to the fracture toughness of the material associated with the corresponding failure mechanism. A characteristic finite element length is used to transform the traction-separation laws into stress-strain laws. The ISV evolution equations are derived in a thermodynamically consistent manner by invoking the stationary principle on the total work of the system with respect to each ISV. A novel feature is the inclusion of both pre-peak damage and appropriately scaled, post-peak strain softening failure. Also, the characteristic elements used in the failure degradation scheme are calculated using the element nodal coordinates, rather than simply the square root of the area of the element.

Pineda, Evan J.↗

ASCR/HEP Exascale Requirements Review Report

This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.

97 MATHEMATICS AND COMPUTING↗

Drop Analysis of Department of Energy Standard Canister with Fort Saint Vrain SNF

DOE manages over 300 types of SNF, most of which are located at the INL site. The Idaho Cleanup Project and INL are collaborating on the Road-Ready Capability Demonstration Project, which will develop and demonstrate the designs, technology, processes, and regulatory framework for packaging DOE-managed SNF at the INL site for ?road-ready dry storage?. ?Road-ready dry storage? is a SNF management concept where SNF is packaged into dry and sealed canisters, which are then placed in on-site storage in anticipation of later transportation. The forward-looking goal of the Demonstration is establishing the foundation for a large-scale road-ready dry storage program at the INL site. The Demonstration will first package Fort Saint Vrain SNF currently stored at INL into several DOE Standard Canisters. These Standard Canisters will then be loaded into another commercial transportation or storage containment system (e.g., storage cask or transportation package). The Standard Canister is a class of standardized canisters designed for containing the large variety of DOE-managed SNF during interim storage, transportation and/or disposal at a geological repository. One critical aspect of road-ready dry storage is the ability to license the DOE Standard Canisters and its associated transportation package to 10 CFR 71. Depending on the SNF and transportation strategy, the Standard Canisters may have to maintain structural integrity under normal conditions of transport and hypothetical accident scenarios (i.e., drop events). The Standard Canisters have been tested and analyzed under various SNF loading configurations and accident drop events in support of the Idaho Spent Fuel Facility and other DOE programs. However, no analysis has been completed to support the recent Demonstration. This analysis will consider the Ø0.5 m × 5.1 m Standard Canister under drop scenario(s) considered in previous INL tests and analyses, including the 9 m drop at 80 degree off vertical. However, this analysis will consider the more recent Fort Saint Vrain loading configurations proposed for the Demonstration. This analysis will performed using strain-based acceptance criteria established by the American Society of Mechanical Engineers Boiler and Pressure Vessel Code, Section III, Division 3. It will be compared to previous analyses and form the foundation of further formal calculations that will be used to support licensing efforts of the road-ready dry storage system at INL.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Center of Excellence Collaboration Projects: Second Wave

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. In particular the NEAMS program is supporting the development of novel thermal-hydraulic codes. The Center of Excellence for Thermal Fluids Application in Nuclear Energy, launched in 2018, has as its key goals to serve as a front door to industry. The Center of Excellence for Thermal Fluids Applications in Nuclear Energy has recently launched a program to start collaborative efforts between the laboratories and industry with the objective of stimulating cooperation and increasing adoption of T/H tools developed under NEAMS by the industry-at-large. In particular, two industry partners agreed to participate in a second wave of short-term collaborations aimed at demonstrating the value of NEAMS tools to their designs. These are Framatome and TerraPower who proposed projects related to LWR (Light Water Reactor) and MSFR (Molten Salt Fast Reactor) designs respectively. In this report we present the results of these two collaborations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerated Adaptive MGS Phase Retrieval

The Modified Gerchberg-Saxton (MGS) algorithm is an image-based wavefront-sensing method that can turn any science instrument focal plane into a wavefront sensor. MGS characterizes optical systems by estimating the wavefront errors in the exit pupil using only intensity images of a star or other point source of light. This innovative implementation of MGS significantly accelerates the MGS phase retrieval algorithm by using stream-processing hardware on conventional graphics cards. Stream processing is a relatively new, yet powerful, paradigm to allow parallel processing of certain applications that apply single instructions to multiple data (SIMD). These stream processors are designed specifically to support large-scale parallel computing on a single graphics chip. Computationally intensive algorithms, such as the Fast Fourier Transform (FFT), are particularly well suited for this computing environment. This high-speed version of MGS exploits commercially available hardware to accomplish the same objective in a fraction of the original time. The exploit involves performing matrix calculations in nVidia graphic cards. The graphical processor unit (GPU) is hardware that is specialized for computationally intensive, highly parallel computation. From the software perspective, a parallel programming model is used, called CUDA, to transparently scale multicore parallelism in hardware. This technology gives computationally intensive applications access to the processing power of the nVidia GPUs through a C/C++ programming interface. The AAMGS (Accelerated Adaptive MGS) software takes advantage of these advanced technologies, to accelerate the optical phase error characterization. With a single PC that contains four nVidia GTX-280 graphic cards, the new implementation can process four images simultaneously to produce a JWST (James Webb Space Telescope) wavefront measurement 60 times faster than the previous code.

Lam, Raymond K.↗