Search NASA⌕ Search

SEARCH · Search NASA

Results for “test case”

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,765 records · Page 98

ISR-1 Seminar: Nuclear Cloud Lofting

This report goes into detail about motivation, background, lofting methodologies (empirical, parcel methods, Navier-Stokes), DELFIC Model, validation, DIORAMA integration, and has the following conclusions: A parcel methodology was applied to develop the DIORAMA cloud lofting module. It is based on the DELFIC model and treats the cloud as a homogenous unit. It solves the set of 8 ODEs for cloud properties. It outputs the time history of the cloud height, radius, and other parameters. The cloud lofting module was tuned with 54 test cases. An iterative brute force search was carried out to find the best fit tuning parameters. The best fit parameters yielded average cloud height errors of 12.9% and 10.3% for the bottom and top, respectively. This all allows for more accurate modeling of the propagation of delayed gamma rays in DIORAMA.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Verification Test of the SURF and SURFplus Models in xRage: Part III Affect of Mesh Alignment

The previous studies used an underdriven detonation wave in 1-dimension (steady ZND reaction zone profile followed by a scale-invariant rarefaction wave) for PBX 9502 as a verification test of the implementation of the SURF and SURFplus models in the xRage code. Since the SURF rate is a function of the lead shock pressure, the question arises as to the effect on accuracy of variations in the detected shock pressure due to the alignment of the shock front with the mesh. To study the effect of mesh alignment we simulate a cylindrically diverging detonation wave using a planar 2-D mesh. The leading issue is the magnitude of azimuthal asymmetries in the numerical solution. The 2-D test case does not have an exact analytic solution. To quantify the accuracy, the 2-D solution along rays through the origin are compared to a highly resolved 1-D simulation in cylindrical geometry.

97 MATHEMATICS AND COMPUTING↗

Warthog: Progress on Coupling BISON and PROTEUS

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program from the Office of Nuclear Energy at the Department of Energy (DOE) provides a robust toolkit for modeling and simulation of current and future advanced nuclear reactor designs. This toolkit provides these technologies organized across product lines, with two divisions targeted at fuels and end-to-end reactor modeling, and a third for integration, coupling, and high-level workflow management. The Fuels Product Line (FPL) and the Reactor Product Line (RPL) provide advanced computational technologies that serve each respective field effectively. There is currently a lack of integration between the product lines, impeding future improvements of simulation solution fidelity. In order to mix and match tools across the product lines, a new application called Warthog was produced. Warthog is built on the Multi-physics Object-Oriented Simulation Environment (MOOSE) framework developed at Idaho National Laboratory (INL). This report details the continuing efforts to provide the Integration Product Line (IPL) with interoperability using the Warthog code. Currently, this application strives to couple the BISON fuel performance application from the FPL using the PROTEUS Core Neutronics application from the RPL. Warthog leverages as much prior work from the NEAMS program as possible, enabling interoperability between the independently developed MOOSE and SHARP frameworks, and the libMesh and MOAB mesh data formats. Previous work performed on Warthog allowed it to couple a pin cell between the two codes. However, as the temperature changed due to the BISON calculation, the cross sections were not recalculated, leading to errors as the temperature got further away from the initial conditions. XSProc from the SCALE code suite was used to calculate the cross sections as needed. The remainder of this report discusses the changes to Warthog to allow for the implementation of XSProc as an external code. It also discusses the changes made to Warthog to allow it to fit more cleanly into the MultiApp syntax of the MOOSE framework. The capabilities, design, and limitations of Warthog will be described, in addition to some of the test cases that were used to demonstrate the code. Future plans for Warthog will be discussed, including continuation of the modifications to the input and coupling to other SHARP codes such as Nek5000.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Turbulence statistics in a spectral element code: a toolbox for High-Fidelity Simulations

In the present document we describe a toolbox for the spectral-element code Nek5000, aimed at computing turbulence statistics. The toolbox is presented for a small test case, namely a square duct with L x = 2h, L y = 2h and L z = 4h, where x, y and z are the horizontal, vertical and streamwise directions, respectively. The number of elements in the xy-plane is 16 X 16 = 256, and the number of elements in z is 4, leading to a total of 1,204 spectral elements. A polynomial order of N = 5 is chosen, and the mesh is generated using the Nek5000 tool genbox. The toolbox presented here allows to compute mean-velocity components, the Reynolds-stress tensor as well as turbulent kinetic energy (TKE) and Reynolds-stress budgets. Note that the present toolbox allows to compute turbulence statistics in turbulent flows with one homogeneous direction (where the statistics are based on time-averaging as well as averaging in the homogeneous direction), as well as in fully three-dimensional flows (with no periodic directions, where only time-averaging is considered).

97 MATHEMATICS AND COMPUTING↗

Device Scale Modeling of Solvent Absorption using MFIX-TFM

Recent climate change is largely attributed to greenhouse gases (e.g., carbon dioxide, methane) and fossil fuels account for a large majority of global CO 2 emissions. That said, fossil fuels will continue to play a significant role in the generation of power for the foreseeable future. The extent to which CO 2 is emitted needs to be reduced, however, carbon capture and sequestration are also necessary actions to tackle climate change. Different approaches exist for CO 2 capture including both post-combustion and pre-combustion technologies, oxy-fuel combustion and/or chemical looping combustion. The focus of this effort is on post-combustion solvent-absorption technology. To apply CO 2 technologies at commercial scale, the availability and maturity and the potential for scalability of that technology need to be considered. Solvent absorption is a proven technology but not at the scale needed by typical power plant. The scale up and down and design of laboratory and commercial packed bed reactors depends heavily on the specific knowledge of two-phase pressure drop, liquid holdup, the wetting efficiency and mass transfer efficiency as a function of operating conditions. Simple scaling rules often fail to provide proper design. Conventional reactor design modeling approaches will generally characterize complex non-ideal flow and mixing patterns using simplified and/or mechanistic flow assumptions. While there are varying levels of complexity used within these approaches, none of these models resolve the local velocity fields. Consequently, they are unable to account for important design factors such as flow maldistribution and channeling from a fundamental perspective. Ideally design would be aided by development of predictive models based on truer representation of the physical and chemical processes that occur at different scales. Computational fluid dynamic (CFD) models are based on multidimensional flow equations with first principle foundations. CFD models can include a more accurate physical description of flow processes and be modified to include more complex behavior. Wetting performance and spatial liquid distribution inside the absorber are recognized as weak areas of knowledge requiring further investigation. CFD tools offer a possible method to investigating such topics and gaining a better understanding of their influence on reactor performance. This report focuses first on describing a hydrodynamic model for countercurrent gas-liquid flow through a packed column and then on the chemistry, heat and mass transfer specific to CO 2 absorption using monoethanolamine (MEA). The indicated model is implemented in MFIX, a CFD open source software package. The user defined functions needed to build this model are described in detail along with the keywords for the corresponding input file. A test case is outlined along with a few results. The example serves to briefly illustrate the developed CFD tool and its potential capability to investigate solvent absorption.

09 BIOMASS FUELS↗

Additional Model Datasets and Results to Accelerate the Verification and Validation of RELAP-7

The RELAP-7 code verification and validation activities are ongoing under the code assessment plan proposed in the previous document (INL-EXT-16-40015). Among the list of V&V test problems in the ‘RELAP-7 code V&V RTM (Requirements Traceability Matrix)’, the RELAP-7 7-equation model has been tested with additional demonstration problems and the results of these tests are reported in this document. In this report, we describe the testing process, the test cases that were conducted, and the results of the evaluation.

97 MATHEMATICS AND COMPUTING↗

Modification of the SAS4A Safety Analysis Code for Integration with the ADAPT Discrete Dynamic Event Tree Framework.

It is difficult to assess the consequences of a transient in a sodium-cooled fast reactor (SFR) using traditional probabilistic risk assessment (PRA) methods, as numerous safety-related sys- tems have passive characteristics. Often there is significant dependence on the value of con- tinuous stochastic parameters rather than binary success/failure determinations. One form of dynamic PRA uses a system simulator to represent the progression of a transient, tracking events through time in a discrete dynamic event tree (DDET). In order to function in a DDET environment, a simulator must have characteristics that make it amenable to changing physical parameters midway through the analysis. The SAS4A SFR system analysis code did not have these characteristics as received. This report describes the code modifications made to allow dynamic operation as well as the linking to a Sandia DDET driver code. A test case is briefly described to demonstrate the utility of the changes.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Block-Parallel Data Analysis with DIY2

DIY2 is a programming model and runtime for block-parallel analytics on distributed-memory machines. Its main abstraction is block-structured data parallelism: data are decomposed into blocks; blocks are assigned to processing elements (processes or threads); computation is described as iterations over these blocks, and communication between blocks is defined by reusable patterns. By expressing computation in this general form, the DIY2 runtime is free to optimize the movement of blocks between slow and fast memories (disk and flash vs. DRAM) and to concurrently execute blocks residing in memory with multiple threads. This enables the same program to execute in-core, out-of-core, serial, parallel, single-threaded, multithreaded, or combinations thereof. This paper describes the implementation of the main features of the DIY2 programming model and optimizations to improve performance. DIY2 is evaluated on benchmark test cases to establish baseline performance for several common patterns and on larger complete analysis codes running on large-scale HPC machines.

97 MATHEMATICS AND COMPUTING↗

Reference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate Boiling

The objective of this effort is to establish a strategy and process for generation of suitable computational mesh for computational fluid dynamics simulations of departure from nucleate boiling in a 5 by 5 fuel rod assembly held in place by PWR mixing vane spacer grids. This mesh generation process will support ongoing efforts to develop, demonstrate and validate advanced multi-phase computational fluid dynamics methods that enable more robust identification of dryout conditions and DNB occurrence.Building upon prior efforts and experience, multiple computational meshes were developed using the native mesh generation capabilities of the commercial CFD code STAR-CCM+. These meshes were used to simulate two test cases from the Westinghouse 5 by 5 rod bundle facility. The sensitivity of predicted quantities of interest to the mesh resolution was then established using two evaluation methods, the Grid Convergence Index method and the Least Squares method. This evaluation suggests that the Least Squares method can reliably establish the uncertainty associated with local parameters such as vector velocity components at a point in the domain or surface averaged quantities such as outlet velocity magnitude. However, neither method is suitable for characterization of uncertainty in global extrema such as peak fuel surface temperature, primarily because such parameters are not necessarily associated with a fixed point in space. This shortcoming is significant because the current generation algorithm for identification of DNB event conditions relies on identification of such global extrema. Ongoing efforts to identify DNB based on local surface conditions will address this challenge.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

PORFLOW 6.42.9 Testing and Verification Document

This report is a continuation of the series of PORFLOW QA documents including (Aleman, 2007; Whiteside, 2007, 2008, 2010a,c,b, 2016b,a, 2017a,b). In this report we have updated the test cases where necessary to account for new features in PORFLOW and to ensure those features perform as described in the documentation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Investigation of Coupled Processes in Argillite Rock: FY18 Progress

The report presents the current status of the TOUGH-FLAC and TOUGH-RBSN simulators for modeling of coupled THM processes in argillite, including fracturing, and verification of the TOUGH-FLAC model against major in situ heater experiments at the Mont Terri Underground Research Laboratory (URL) in Switzerland and at the URL in Bure, France. The heater experiments modeled are the Mont Terri FE (Full-scale Emplacement) Experiment, conducted as part of the Mont Terri Project, and the TED and ALC experiments conducted in Callovo-Oxfordian claystone (COx) at the Meuse/Haute-Marne (MHM) URL in Bure, France. The modeling of the TED and ALC heater experiments is part of a modeling task (Task E) of the international DECOVALEX-2019 project. DECOVALEX, which stands for DEvelopment of COupled Models and their VALidation against EXperiments, is an international collaborative activity, in which DOE and LBNL gain access to unique laboratory and field data defined as modeling test cases. DECOVALEX tasks are studied collectively among several international groups to better understand the processes and to improve numerical models that could eventually be applied in the performance assessment for nuclear waste disposal in clay host rocks and bentonite backfill. The report also includes a description of LBNL’s activities for modeling gas migration in bentonite related to Task A of the international DECOVALEX project, and presents new FY2018 work related to monitoring and characterization of faults using a unique borehole tool.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Predicting Ice Sheet and Climate Evolution at Extreme Scales (PISCEES) Final Technical Report for effort at the University of Texas at Austin Award DE-SC0008083 March 2016 - August 2017

Mass loss from the Greenland and Antarctic ice sheets is accelerating. Although ice sheet models have improved in recent years, much work is needed to make these models robust and efficient on continental scales and to quantify uncertainties in their projected outputs. The University of Texas at Austin's role in the PISCEES project is to provide guidance on how best to apply uncertainty quantification strategies to evaluating uncertainties in simulations of Greenland and Antarctic contributions to sea level rise in response to environmental forcing. Here we report on several activities that involved project members at UTA including Section 2.1 where we provide a broader overview of challenges to use UQ methods on ice ow modeling. Within Section 2.2 we describe the progress and challenges that remain for implementing important aspects of our "end-to-end" UQ work flow. In sum, we had a break-through in realizing the potential for implementing Active Subspace strategies making use of existing experiments exploring the sensitivity of the Greenland ice sheet flow to uncertain basal traction parameters. Section 2.3 presents progress that is being made to explore the effect of uncertainties in bed topography on stability of Thwaites basin, West Antarctica, using the BISICLES ice flow model. This effort is joint with Dan Martin at LBNL and a graduate student being co-mentored by Jackson. The results provide a compelling test case for the BISICLES ice flow model to address an important science question concerning the role of uncertain boundary conditions on ice sheet contributions to sea level rise.

54 ENVIRONMENTAL SCIENCES↗

Research, Development, and Testing of a Master Controller with Applications to the Bronzeville Community Microgrid System

The objective of this project was to research, develop and test a Microgrid Master Controller (MMC) with applications to the Bronzeville Community Microgrid (BCM), capable of monitoring and controlling a microgrid cluster, as well as to demonstrate its performance in an environment that will closely resemble actual operation of the microgrid. During the project, a hierarchical control strategy of the MMC consisting of primary, secondary, and tertiary layers of control was developed, to operate the MMC in conjunction with the IIT microgrid controller for the coordinated operation of the microgrid cluster, including the control of DERs and facilitation of power flow between the microgrids. A detailed laboratory test plan, including methodology, test setup, test scenarios, and techniques for data gathering and analysis, was developed to verify the MMC functionalities and to validate the performance of the control algorithms under various microgrid operating conditions and potential system contingencies. The MMC underwent a thorough lab testing with more than 320 test cases assigned to 14 performance categories. The final project review and technology feasibility demonstration was held with Department of Energy (DOE) team, in the ComEd’s Grid of the Future lab. During this review, technical performance of the MMC was demonstrated in all key scenarios where the proper performance of the MMC was witnessed under a wide array of scenarios, including planned islanding, unplanned islanding, re-synchronization, load restoration, and basic dispatch in grid-connected and islanded mode. The demonstration was based on use cases focused on main microgrid capabilities and features, which further and the microgrid performance. The following cases were shown to DOE: Planned Islanding; Re-synchronization; Load Restoration; Basic Dispatch (grid-connected mode); Basic Dispatch (islanded mode); Unplanned Islanding. The outcome of this project is a viable microgrid control technology that served as a basis for further enhancement and commercialization of the utility-grade microgrid controller that will be employed for control of the BCM.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Ice-Sheet / Ocean Interaction Model for Greenland Fjords Using High-Order Discontinuous Galerkin Methods (Final Report)

This report summarizes the development and validation of the Non-hydrostatic Unified Model of the Ocean (NUMO). The overarching science goal of the project was to advance Earth System science by proposing a more realistic representation of physical processes in a Greenland fjord. To achieve this objective, we proposed to build a new, separate, high-resolution ocean model to realistically represent the complex three-dimensional physics of the fjord, its bathymetry, coastlines using continuous/discontinuous Galerkin (CG/DG) methods. The report describes the technical details of model components, i.e., choice of equations set, time integration scheme, boundary conditions, spatial discretization methods, mesh generation. We show a suite of validation results for test cases ranging from the Kovasznay flow, through the lock-exchange, simplified two-dimensional fjord, to a three-dimensional density current on a slope. Results are compared with the literature. The report is concluded with a discussion of efforts taken to perform a simulation of Sermilik Fjord with realistic bathymetry data.

58 GEOSCIENCES↗

Final Report: Modeling coupled ice sheet-ocean interactions in the Model for Prediction Across Scales (MPAS) and in DOE Earth System Models

Ice sheet-ocean interactions are involved in a potential climate tipping point in the climate system. Weakening of dynamical ocean barriers that currently hold warmer waters off the Antarctic continental shelf could lead to warming under Antarctica’s largest ice shelves, causing as much as a hundredfold increase in sub-ice-shelf melting. Increased ice-sheet melting could trigger dynamical instabilities in the ice sheet itself if sufficient buttressing from ice shelves is lost. This project developed and tested several modeling capabilities needed to investigate these phenomena. In collaboration with a number of international groups, notably the climate modeling group at Los Alamos National Laboratory (LANL), the following capabilities were added to the Energy Exascale Earth System Model (E3SM) and its ocean component (the Model for Prediction Across Scales Ocean, MPAS-Ocean): 1) a new vertical coordinate appropriate for ice-shelf cavities; 2) an initialization procedure appropriate for damping waves from imbalances in ocean pressure and that of the overlying ice shelf; 3) A suite of idealized test cases related to ice shelf-ocean interactions; 4) A set of diagnostics for tracking freshwater fluxes and heat budgets from ice-shelf melting; 5) infrastructure for computing ice sheet-ocean fluxes in the E3SM coupler; and 6) A set of analysis and visualization tools for E3SM ocean and sea-ice validation with a particular focus on Antarctic and Southern Ocean climate. Additionally, the project included participation and leadership in a number of international model intercomparison projects (MIPs). Most notably, the PI was a co-chair of the Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP), designing the experiment, coordinating participation in the MIP and performing the bulk of the analysis. Although the project was cut short when the PI accepted a staff scientist position at LANL, the project played a role in seven publications and has segued smoothly into ongoing Department of Energy projects in the “ecosystem” surrounding E3SM.

54 ENVIRONMENTAL SCIENCES↗

A comparison of anisotropic and isotropic grain growth simulation results of Uranium dioxide and Copper

Grain boundaries migrate to reduce the total grain boundary energy in a material. However, grain boundary energy is anisotropic, such that some grain boundaries are less energy favorable in the system and tend to disappear. Idaho National Laboratory’s (INL) mesoscale fuel performance code MARMOT currently assumes isotropic (non-directionally dependent) grain boundary energy in all phase field sim- ulations of grain boundary migration. As a preliminary step to incorporating fully anisotropic grain boundary energies into MARMOT, a comparison of the results of grain growth simulations for UO 2 nuclear fuels utilizing the old isotropic MARMOT code and new code incorporating anisotropic grain boundary energies was performed. Comparisons of Copper simulations were also performed. The re- sults analyzed were from simple circle and thumb-shaped grain test cases, and polycrystalline systems using artificially generated initial conditions. The metrics used to compare the grains from isotropic and anisotropic cases include grain area and average grain boundary energy. Anisotropic behavior was observed in the systems of circle and thumb-shaped grains. Polycrystalline systems do not behave as expected, suggesting missing terms in the equation implemented in MARMOT which governs the mi- crostructural evolution of the grain boundaries.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computational Assessment of Hydrodynamic Loads on Rockeries for River Bank Protection

Colorado experienced an extreme flooding event in September of 2013 that ravaged canyon roadways southeast of Estes Park, CO. Central Federal Lands was asked to rebuild SH 34 and CO 43 after the floods. To minimize river encroachments in narrow canyons, roadway fill embankments must be over steepened. Using a rockery to stabilize the roadway embankment rather than traditional vertical, concrete structures was the solution. However, no analytical techniques were available for evaluating the hydraulic stability of such structures. Argonne National Laboratory was employed to estimate the hydrodynamic forces by using advanced, computational fluid dynamics (CFD). An investigation was conducted on a section of a model rockery comprised of large boulders exposed to a variety of hydraulic flow conditions. A repeatable geometric section of a typical rockery structure was modeled and used to create boundaries for the CFD model. The test case matrix covers computer simulations of various rockery and channel geometries. The geometrical parameters taken into account are: rockery height, outside face batter angle of the rockery, rockery base width, river bed width, and the hydraulic parameters are: flow depth, flow velocity, and hydraulic angle of attack. Two failure mechanisms of rockery structures, sliding and overturning, are considered. The existing design equations for both modes are modified by addition of hydrodynamic force terms. The new terms were calculated for a parametric set of geometry and condition variation using the CFD model. Moreover, formulas for dimensionless hydraulic force coefficients are proposed that can be used in the design of new rockeries or in evaluation of existing dry-stack rockeries in river environments.

42 ENGINEERING↗

Demonstration of Fully Coupled Calculations in Helical Steam Generator: Toward Predictions of Fluid Elastic Instability

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. Flow-induced vibration (FIV) is widespread problem in energy systems because they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occurs. Given the importance of reliable components in the nuclear industry, flow-induced vibration has long been a major concern in safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from flow-induced vibration and related failures. Advanced reactors, such as integral Pressurized Water Reactors (PWRs) considered for Small Modular Reactors (SMR), often rely on innovative component designs to meet cost and safety targets. One component that is the subject of advanced designs is the steam generator, some designs of which forego the usual shell-and-tube architecture in order to fit within the primary vessel. In addition to being more cost- and space-efficient, such steam generators need to be more reliable, since failure of the primary vessel represents a potential loss of coolant and a safety concern. A significant amount of data exists on flow-induced vibration in shell-and-tube heat exchangers, and heuristic methods are available to predict their occurrence based on a set of given assumptions. In contrast, advanced designs have far less data available. Advanced modeling and simulation based on coupled structural and fluid simulations have the potential to predict flow-induced vibration in a variety of designs, reducing the need for expensive experimental programs, especially at the design stage. Over the past seven years, the Reactor Product Line has developed the integrated multi-physics code suite SHARP. The goal of developing such a tool is to perform multi-physics neutronics, thermal/fluid, and structural mechanics modeling of the components inside the full reactor core or portions of it with a user-specified fidelity. In particular SHARP contains high-fidelity single-physics codes Diablo for structural mechanics and Nek5000 for fluid mechanics calculations. Both codes are state-of-the-art, highly scalable tools that have been extensively validated. These tools form a strong basis on which to build a flow-induced vibration modeling capability. In the last two years a significant effort was devoted to 1) develop and validate a fully coupled capability in SHARP for Fluid-Structure calculations (FSI) and 2) validate the use of SHARP for one-way coupled calculations in helical steam generators on legacy data available as well as novel experiments. In this report we discuss an extension of previous one-way coupled calculations performed with Nek5000 and Diablo aimed at simulating available FIV experiments in helical steam generators in the turbulent buffeting regime. It was demonstrated that in the buffeting regime one-way coupling was judged sufficient because the pressure loads do not cause substantial displacements. However, higher speeds led to a reduction in accuracy and the need for a fully coupled capability for velocity or approaching the critical velocity and fluid elastic instability. In this report we demonstrate such fully coupled capability on the Argonne test case discussed in previous reports.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗