Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modelling Capabilities”

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,531 records · Page 85

NEAMS Burnup Extension Accomplishments and Remaining Modeling Gaps

The economic viability of light-water reactors (LWRs) in the United States is declining in heavily subsidized markets, and as a result, the nuclear industry is looking for opportunities to enhance the economic competitiveness of nuclear power. This is not a foreign concept to the nuclear industry: in the mid-2000s, the nuclear industry set out to achieve zero fuel failures by 2010. The goal in this effort was to drive down the cost of reactor shut down by replacing a pin or bundle in response to fuel rod failure. 2010 brought about the initiative to deliver the nuclear promise to reduce operating cost by 30% to improve nuclear energy’s economic competitiveness before 2020. The emergence of accident-tolerant fuel also offers the nuclear industry an opportunity to build on these past successes and deliver affordable, clean energy. Accident-tolerant fuel has been shown to provide superior performance compared to traditional Zircaloy/UO2 fuel concepts, offering the unique ability to remove operational limitations that inhibit the economic viability of nuclear power. This has led the industry to begin building a technical case to extend the peak rod average burnup beyond 62 GWd/tU to extend pressurized water reactor cycle lengths to 24 months and to develop more efficient boiling water reactor core designs. The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program mission is to develop advanced modeling and simulation tools and capabilities to accelerate the deployment of advanced nuclear energy technologies. The primary safety concern inhibiting the nuclear industry from extending burnup is related to high-burnup fuel fragmentation, relocation, and dispersal. Therefore, the NEAMS program developed a targeted 5-year plan to support the industry’s efforts to extend burnup. This milestone report summarizes the 5-year plan that was enacted in FY20, followed by a discussion of the ongoing activates required to fulfill the 5-year plan, as well as the approach to address the current modeling gaps. Additionally, an LWR stakeholder meeting was held to communicate work performed in the NEAMS program over the past three years, to assess the LWR community’s perspective on the impact of the program, and to identify remaining significant gaps in the NEAMS suite of capabilities. This engagement will be documented by the Electric Power Research Institute and used by NEAMS to redirect current LWR scope as needed and to develop the next phase for LWR research and development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A mathematical model of asynchronous data flow in parallel computers *

Abstract We present a simplified model of data flow on processors in a high-performance computing framework involving computations necessitating inter-processor communications. From this ordinary differential model, we take its asymptotic limit, resulting in a model which treats the computer as a continuum of processors and data flow as an Eulerian fluid governed by a conservation law. We derive a Hamilton–Jacobi equation associated with this conservation law for which the existence and uniqueness of solutions can be proven. We then present the results of numerical experiments for both discrete and continuum models; these show a qualitative agreement between the two and the effect of variations in the computing environment’s processing capabilities on the progress of the modelled computation.

Barnard, Richard C.↗

Enhancement of phase change material hysteresis model: A case study of modeling building envelope in EnergyPlus

Nowadays, buildings are expected to offer demand side services to the power grid to enhance the electrical load flexibility, which leads to the concepts of grid-interactive efficient buildings (GEBs). Phase change material (PCM)-based thermal energy storage has seen increasing attention in recent years for peak load shifting of grid-interactive efficient buildings (GEBs). Numerical models are critical tools for design and evaluation of PCM-integrated systems. Most industrial-grade PCMs are reported to melt/freeze over a temperature range instead of at a unique temperature. Such thermal hysteresis effect significantly affects the reliability of simulation results because not only the heat transfer process depends on melting and freezing temperatures, the PCM thermal properties change significantly during the phase change process as well. This study is aimed to develop a model for the PCMs used in the building envelope with the capability to accurately simulate hysteretic behaviors. Further, this model is based on a two-phase assumption and is implemented in a whole building energy performance simulation program (i.e., EnergyPlus). A comparison between numerical results and experimental data shows that during a complete phase transition, the two-phase model could achieve a good agreement with the experimental data. During a partial phase transition, the two-phase model could lead to significant improvements compared to other alternative PCM models, including the existing PCM model in EnergyPlus. Last, whole building simulations were performed to study this model's performance regarding heating/cooling loads and zone mean air temperature of a given building. The results show that the difference in hourly heating/cooling loads introduced by the models was less than 1% in design conditions, while significant changes were observed in both hourly heating/cooling loads and zone mean air temperature when the PCM envelope underwent partial phase transition processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

RELAP5-3D Capabilities for High-Temperature Gas-Cooled Reactor Analysis

RELAP5-3D is the premier systems thermal-hydraulic code in the world. The code has been used significantly for high-temperature gas-cooled reactor (HTGR) analyses over the last 10 years, but there has not been a comprehensive summary of RELAP5-3D capabilities and methods for HTGR modeling. In this paper, we provide a guide for how novice users can build models for prismatic HTGRs. This guide includes best practices for developing heat structures and setting up radiation and conduction heat transfer. We also provide a summary of relevant validation studies. These studies are divided into three categories: separate effects tests, integral effects tests, and nuclear reactor studies. We present a comprehensive review of RELAP5-3D HTGR applications for reactor modeling over the last 10 years. These applications include work for the Modular High-Temperature Gas-Cooled Reactor benchmark, design activities for a prismatic HTGR as part of the U.S. Department of Energy’s Advanced Demonstration and Test Reactor Options Study, and studies performed to design gas-cooled microreactors. Following the comprehensive literature review, we present a gap analysis highlighting capabilities of RELAP5-3D that are available and missing for a complete analysis of HTGR thermal hydraulics. We also discuss the limitations of existing code capabilities and potential improvements that could improve analyses. Current research efforts have yet to fully validate RELAP5-3D for a prismatic HTGR analysis; however, some of this is attributable to confounding factors in the available experimental data as opposed to code capabilities. In conclusion, this paper presents the most comprehensive review to date on RELAP5-3D capabilities, limitations, and applications for HTGR analysis.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Real-Time Anomaly Detection for Searches Beyond the Standard Model in the ProtoDUNE Horizontal Drift Detector

This paper summarizes work conducted throughout a SULI internship at Fermi National Accelerator Laboratory focused on building an unsupervised machine learning model for real-time anomaly detection in ProtoDUNE Horizontal Drift. Using simulated data, we trained an autoencoder model on a pure cosmic dataset, and evaluated it on both cosmic and neutrino events—making the model an anomaly detector. The goal was to make a model which matches or exceeds the current ADC Simple Window trigger algorithm so that our model can perform at the same rate but provide sensitivity to potential beyond-the-Standard-Model (BSM) signatures. In the end, we were able to construct a model which slightly exceeds the capabilities of the ADC Simple Window while remaining completely unsupervised, achieving 31.9 ± 0.2% (26.6 ± 0.2%) ν efficiency at 5 Hz (2 Hz), a 3.6 (3.2) percentage point increase. Additionally, 17.5 ± 0.3% (18.3 ± 0.3%) of the events that passed the autoencoder at 5 Hz (2 Hz) were missed by the current trigger algorithm. Future work will investigate alternative normalization methods, including quantile transformation, and evaluate the model on ProtoDUNE-HD detector-glitch data if that data becomes available.

Wilson, C. [Cincinnati U., RWC]↗

Modelling the Value of Networked, Small-Hydro Generators for Storage and Dispatch Capability on the Grid

Natel Energy, Inc. (Natel), Alameda, California, has developed a new type of fish-friendly, low-head hydropower turbine, the hydroEngine (R), that is designed to provide reliable, dispatchable power and ancillary services at significantly reduced civil works costs. For the U.S. Department of Energy's (DOE) Small Business Voucher for Hydropower, Third Round, Natel and NREL will model and evaluate the potential value streams of a portfolio of networked small hydropower generators. The project is designed to advance the understanding of the value of small hydro generation to the grid of the future. It will help support dialog with potential customers and stakeholders regarding the value of the Natel's technologies and inform future R&D efforts in this area.

13 HYDRO ENERGY↗

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗

Real-Time Anomaly Detection for Beyond Standard Model Searches in ProtoDUNE Horizontal Drift

This paper summarizes work conducted throughout a SULI internship at Fermi National Accelerator Laboratory focused on building an unsupervised machine learning model for real-time anomaly detection in ProtoDUNE Horizontal Drift. Using simulated data, we trained an autoencoder model on a pure cosmic dataset, and evaluated it on both cosmic and neutrino events---making the model an anomaly detector. The goal was to make a model which matches or exceeds the current ADC Simple Window trigger algorithm so that our model can perform at the same rate but provide sensitivity to potential beyond-the-Standard-Model (BSM) signatures. In the end, we were able to construct a model which slightly exceeds the capabilities of the ADC Simple Window while remaining completely unsupervised, achieving $31.9 \pm 0.2$\% ($26.6 \pm 0.2$\%) $\nu$ efficiency at 5 Hz (2 Hz), a 3.6 (3.2) percentage point increase. Additionally, $17.5 \pm 0.3$\% ($18.3 \pm 0.3$\%) of the events that passed the autoencoder at 5 Hz (2 Hz) were missed by the current trigger algorithm. Future work will investigate alternative normalization methods, including quantile transformation, and evaluate the model on ProtoDUNE-HD detector-glitch data if that data becomes available.

Wilson, Cameron C. [Cincinnati U., RWC]↗

Real-Time Anomaly Detection for Beyond Standard Model Searches in ProtoDUNE Horizontal Drift

This paper summarizes work conducted throughout a SULI internship at Fermi National Accelerator Laboratory focused on building an unsupervised machine learning model for real-time anomaly detection in ProtoDUNE Horizontal Drift. Using simulated data, we trained an autoencoder model on a pure cosmic dataset, and evaluated it on both cosmic and neutrino events---making the model an anomaly detector. The goal was to make a model which matches or exceeds the current ADC Simple Window trigger algorithm so that our model can perform at the same rate but provide sensitivity to potential beyond-the-Standard-Model (BSM) signatures. In the end, we were able to construct a model which slightly exceeds the capabilities of the ADC Simple Window while remaining completely unsupervised, achieving $31.9 \pm 0.2$\% ($26.6 \pm 0.2$\%) $\nu$ efficiency at 5 Hz (2 Hz), a 3.6 (3.2) percentage point increase. Additionally, $17.5 \pm 0.3$\% ($18.3 \pm 0.3$\%) of the events that passed the autoencoder at 5 Hz (2 Hz) were missed by the current trigger algorithm. Future work will investigate alternative normalization methods, including quantile transformation, and evaluate the model on ProtoDUNE-HD detector-glitch data if that data becomes available.

Wilson, Cameron C. [Cincinnati U., RWC]↗

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale Aerosol Modeling Across Space and Composition (Final Report)

The objective of this project was to harness the capabilities of our unique particle-resolved aerosol model PartMC as an integrating tool between models and measurements across scales. This particle-resolved model explicitly resolves aerosol aging processes that must be approximated in large-scale models, making it ideal for quantifying structural uncertainty and rigorously deriving simplified parameterizations (“coarse graining”) to be used in large-scale models. It also directly connects with laboratory and field measurements at both particle and bulk scales. The significant accomplishments of this project are as follows: 1. Advancing the state of the art in particle-resolved modeling; 2. Connecting modeling with ASR/ARM measurements; 3. Bridging from the process level to aerosol impacts on the global scale; 4. Community engagement through leading and participating in ASR program activities and publishing a review paper on aerosol mixing state that synthesizes experimental and modeling perspectives.

58 GEOSCIENCES↗

Future Projections of Lifecycle Cost and Greenhouse Gas Emissions of Light-Duty Vehicles

Vehicles with electrified powertrains carry the promise of significant reductions in greenhouse gas (GHG) emissions from a lifecycle analysis (LCA) standpoint compared to conventional internal combustion engine (CICE) vehicles. However, trade-offs exist between different types of electrified powertrains in terms of cost, consumer acceptance, and GHG reduction efficacy for different operating conditions. The open-source tool CarGHG was developed with an aim to enable the exploration of a plethora of parametric study scenarios, including the cost of electrification technologies, different driving patterns and charging habits, and the cost and carbon intensity of electricity and fuel blends. This paper introduces the framework of CarGHG, then showcases total cost of ownership (TCO) and LCA GHG results for select models of light-duty vehicles. Another capability of CarGHG, which is the ability to estimate the performance of “virtual” vehicle models (perceived vehicle design specifications not yet on the market), is utilized to explore future scenarios of electrification and low-carbon fuel blends for Small Sports Utility Vehicles (SUVs), a popular light-duty vehicle segment in North America. With opportunities, but also uncertainties, in future scenarios, it is likely wise to continue pursuing multiple ways towards the reduction of LCA GHG.

Hamza, Karim↗

Modeling Dense Particle Flow in Multistage and Obstructed Flow Receivers Using High Fidelity Simulations

Particles are a leading contender for next-generation, concentrating solar power technologies, and the design of the particle receiver is critical to minimize the levelized cost of electricity. Falling particle receivers (FPRs) are a viable receiver concept, but many new designs feature complex particle obstructions that include dense discrete phase flows. This creates additional challenges for modeling as particle-to-particle interactions (i.e., collisions) and particle drag become more complex. To improve upon existing modeling strategies, a CFD-DEM simulation capability was created by coupling two independent codes: Sierra/Fuego and LAMMPS. A suitable receiver model was then defined using a traditional continuum-based model for the air and a granular model for the particle curtain. A sensitivity study was executed using this model to determine the relevance of different granular model inputs on important quantities of interest in obstructed flow FPRs: the particle velocity and curtain opacity. The study showed that the granular model inputs had little effect on the particle velocity magnitude and curtain opacity after an obstruction.

Mills, Brantley↗

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhanced mechanical property evaluation using innovative data analytics capability

This report focuses on efforts to improve Multiphysics Object-Oriented Simulation Environ-ment (MOOSE) for mechanical property evaluation using data analytics. These efforts involveimprovements to the stochastic tools module (STM) for stochastic simulations of MOOSE mul-tiphysics model and the development of inverse optimization capabilities. The report gives anoverview of the STM and describes recent updates to its core capabilities and theory on itsreduced-order model (ROM) schemes. Examples are also provided showing the impact of theseupdates and exhibits the usefulness of ROMs. An overview of the gradient based inverse opti-mziation algorithms are given along with examples of their application to source identification.Inverse optimization will provide a new methodology in STM for fitting model parameters toexperimental data.

97 MATHEMATICS AND COMPUTING↗

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

FY20 Verification of BISON Using Analytic and Manufactured Solutions

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗