Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Comparative Analysis of HEATNETS for Geothermal Network Performance

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification and validation process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study focuses on a validation of HEATNETS using results from TRNSYS. HEATNETS is not intended to replace other modeling tools, but this work demonstrates, via a comparison with an industry standard code, that HEATNETS can be a unique, high-level and rapid modeling tool for estimating the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

GANpiler

Focal Area(s): Rather than augment or replace physical models with machine learning, we instead preserve the existing models and augment the underlying code with faster surrogate models created by generative adversarial networks (GAN). To leverage the performance optimization, we also propose a runtime system and user interfaces that allow prediction and tracking of accumulated error as well as dynamic, per-process decision making as to which model (if any) to use for each iteration. Science Challenge: This proposal sits at the nexus of two hard problems. First, climate models based on machine learning will be, by their nature, difficult to trust once their predictions begin diverging from the consensus. Second, compilers and hardware have been making only incremental performance gains for decades. GPGPUs have provided a welcome performance boost for codes that can take advantage of them, but there is no similar technology on the horizon to provide the next performance leap.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗
DOE OSTI2023

Implementation and verification of a model linearized multi-species collision operator in the COGENT code

We implement the multi-species linearized model collision operator, based on the approach proposed by Kolesnikov, in the full-f continuum gyrokinetic code COGENT. We describe the operator and use it in several COGENT simulations of highly collisional plasmas. We analyze simulation results with the Braginskii fluid model to illustrate that COGENT recovers both friction and thermal forces. Finally, we simulate the neoclassical transport of heavy impurities with COGENT and illustrate agreement with previously published theoretical results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Prototype Thick-Target Bremsstrahlung Model with Angularly-Dependent Emission in the MCNP6 ® Code

This document summarizes the current thick-target bremsstrahlung (TTB) model in MCNP and provides test results for an alternative implementation to improve the accuracy with reduced cost compared to full electron transport. It has been observed that the current TTB model produces inaccurate results in problems where the medium is thick with respect to electrons, but the photon distribution in the problem has a strong directionality. An example of such a simulation is detectors surrounding a metal target irradiated with a radiographic beam of high energy photons. The primary cause of this discrepancy is the current TTB method emits all bremsstrahlung photons in the same direction as the primary electron produced from each (γ, e ± ) interaction, leading to artificially forward peaked photon distributions for intermediate to high-energy incident photons. To improve the TTB model, we have implemented an angularly-dependent TTB model in a developer version of the MCNP6 ® code; for developers, this was done on the branch prototype/angular_ttb in the mcnp6 repo on bitbucket. The angularly-dependent TTB model accounts for the energy and scattering of electrons as they slow down in the current material, but does not sample the computationally expensive energy straggling, secondary electron events, and tracking electrons; this approach is significantly less computationally expensive than full electron transport and can be comparable to the original TTB method for problems with sufficiently complex materials and geometry. To evaluate the method, we have modeled a simple problem of a beam of 5 MeV photons incident on a sphere of plutonium surrounded by detectors at different deflection angles. For this problem, the angularly-dependent TTB produces a photon flux within 8.0% for a 90 degree deflection angle and 0.8% along the beam axis, as compared to the electron transport solution. This is an improvement compared to a 43% and 81% discrepancy with the original TTB method, respectively. The rest of this work includes the following: the first section details the current TTB treatment in MCNP, which has not been well documented elsewhere. Then, the modified TTB algorithm is detailed and the approximations compared to the condensed history algorithm are compared. Results are given comparing the two TTB methods to the condensed history transport algorithm. The appendix includes details for code developers on relevant electron transport implementation details and potential code improvements for future work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Multiphysics Simulation of KRUSTY Warm Critical Experiments Using MOOSE Tools

A series of critical experiments were performed during the Kilopower Reactor Using Stirling TechnologY (KRUSTY) project, producing valuable and unique data for the validation of corresponding microreactor models. In this paper, the development of a high-fidelity multiphysics model of the KRUSTY microreactor within the MOOSE ecosystem is presented, with a focus on the microreactor configuration during the warm critical experiments. The model was established by coupling a deterministic neutronics model using the Griffin code, with a BISON thermomechanical model. The multiphysics model successfully simulated the KRUSTY 15 ¢ and 30 ¢ warm critical reactivity insertion experiments in satisfactory agreement with the experimental measurement data. In conclusion, the modeling results demonstrated the performance of MOOSE-based multiphysics tools in microreactor simulations, and spotlighted the useful applications of these tools in supporting microreactor design, qualification, and licensing activities.

KRUSTY↗
DOE OSTI2025

Amorphization and siliconization of silicon carbide as a first wall material

The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

TMS XTRRA experiments on NIF

The study of x-ray induced rapid heating and the resulting material responses - (thermomechanical shock (TMS) and thermo-structural response (TSR)) - is important to our understanding of material behavior in extreme environments. A series of high conversion efficiency (XRCE) x-ray sources have been developed for the NIF laser facility ranging in fluence and spectral content. In conjunction with source development, we have developed a new XTRRA (X-ray Transport and Radiation Response Analysis) test cassette to hold a set of samples, each at equal distance from target chamber center (TCC) where the x-ray source is placed. Design and performance details of the x-ray sources is discussed elsewhere. Here, we describe the development and use of this TMS XTRRA test cassette. The XTRRA test cassette is a snout mounted on a NIF Diagnostic Instrument Manipulator (DIM) designed to position six samples equidistant from the x-ray source to assure uniform irradiation of all samples. Photon Displacement interferometry (PDI) is used to measure the displacement of the rear surface of each sample resulting from x-ray exposure on the front of the sample. This will help determine the x-ray generated impulse on the sample. The samples are retrieved from the snout after the shot for inspection and post-shot analysis including the measurement of sample ablated mass and engineering tests examining changes in the material structure. The data from experiments are compared to predictions from LLNL’s ASC codes and used to inform models implemented in these codes. Thus, XRCE x-ray sources with the XTRRA snout at the National Ignition Facility (NIF) laser at the Lawrence Livermore National Laboratory [1,2] provide the ability to expose a broad range of material samples directly to high x-ray fluence, measure shock transits, and recover the sample for further analysis post-shot.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermo-Mechanical Shock, X-Ray Transport and Radiation Response Analysis Experiments on NIF

The study of x-ray induced rapid heating and the resulting material responses - (thermomechanical shock (TMS) and thermo-structural response (TSR)) - is important to our understanding of material behavior in extreme environments. A series of high conversion efficiency (XRCE) x-ray sources have been developed for the NIF laser facility ranging in fluence and spectral content. In conjunction with source development, we have developed a new XTRRA (X-ray Transport and Radiation Response Analysis) test cassette to hold a set of samples, each at equal distance from target chamber center (TCC) where the x-ray source is placed. Design and performance details of the x-ray sources is discussed elsewhere. Here, we describe the development and use of this TMS XTRRA test cassette. The XTRRA test cassette is a snout mounted on a NIF Diagnostic Instrument Manipulator (DIM) designed to position six samples equidistant from the x-ray source to assure uniform irradiation of all samples. Photon Displacement interferometry (PDI) is used to measure the displacement of the rear surface of each sample resulting from x-ray exposure on the front of the sample. This will help determine the x-ray generated impulse on the sample. The samples are retrieved from the snout after the shot for inspection and post-shot analysis including the measurement of sample ablated mass and engineering tests examining changes in the material structure. The data from experiments are compared to predictions from LLNL’s ASC codes and used to inform models implemented in these codes. Thus, XRCE x-ray sources with the XTRRA snout at the National Ignition Facility (NIF) laser at the Lawrence Livermore National Laboratory provide the ability to expose a broad range of material samples directly to high x-ray fluence, measure shock transits, and recover the sample for further analysis post-shot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development of a Turbulent Liquid Spray Atomization Model for Diesel Engine Simulations (Final Technical Report)

This project addresses the systematic lack of predictive capabilities by spray models within engine CFD codes. We develop a new modeling approach to predict the breakup of diesel sprays based on recent literature showing that liquid turbulence plays a fundamental role in spray atomization. A new body of quantitative validation data is also developed as a critical element of the project, leveraging the joint capabilities of Georgia Tech’s high-pressure continuous-flow spray chamber and Argonne National Lab’s near-nozzle x-ray diagnostics at the Advanced Photon Source. This project contributes spatially-resolved measurements of drop size distribution within well-characterized diesel injectors, Spray A and D, from the Engine Combustion Network (ECN) to the engine combustion community for the first time. Utilizing this new body of measurements, we validate and demonstrate a new spray model for diesel sprays, termed the KH-Faeth model, that predicts global and local spray characteristic more accurately than the widely adopted and employed KH model. Predicted drop size distributions are seen to predict measured drops sizes both quantitatively and predictively, with accurate response in droplet size distributions over a wide range of ambient density, injection pressure, and injector nozzle size (Spray A and D) without model tuning. The KH-Faeth model can reduce error in the predicted centerline droplet size profile by up to 80% for ECN Spray D simulations when compared to use of the widely employed KH model.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

TRISO Fuel Performance Modeling

Overview of TRISO fuel performance modeling TRISO fuel performance code PARFUME TRISO fuel performance modeling with BISON PARFUME and BISON comparison

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Standards, dissemination, and best practices in systems biology

In this study, the reproducibility of scientific research is crucial to the success of the scientific method. Here, we review the current best practices when publishing mechanistic models in systems biology. We recommend, where possible, to use software engineering strategies such as testing, verification, validation, documentation, versioning, iterative development, and continuous integration. In addition, adhering to the Findable, Accessible, Interoperable, and Reusable modeling principles allows other scientists to collaborate and build off of each other’s work. Existing standards such as Systems Biology Markup Language, CellML, or Simulation Experiment Description Markup Language can greatly improve the likelihood that a published model is reproducible, especially if such models are deposited in well-established model repositories. Where models are published in executable programming languages, the source code and their data should be published as open-source in public code repositories together with any documentation and testing code. For complex models, we recommend container-based solutions where any software dependencies and the run-time context can be easily replicated.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Modeling heat pipe startup and noncondensable gases in Sockeye

For this work, a one-dimensional gas mixture flow model was developed and implemented in the heat pipe code Sockeye to model the effects of noncondensable gases. Additionally, a startup model based on the dusty gas model was implemented to model the transition from rarefied gas dynamics to continuum flow, which occurs during the frozen startup of high-temperature heat pipes. Multiple startup and noncondensable gas models were tested against experimental data for sodium heat pipes, showing excellent agreement. Additionally, the newly developed gas mixture model for modeling noncondensable gas is further tested with a theoretical case study with arbitrary heating configurations. Finally, several recommendations and conclusions are made from the studies in this work to guide future heat pipe modeling efforts.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Modeling Heat Pipe Startup And Noncondensable Gases In Sockeye

A one-dimensional gas mixture flow model was developed and implemented in the heat pipe code Sockeye to model the effects of non-condensable gases. Additionally a startup model based on the dusty gas model was implemented to model the transition from rarefied gas dynamics to continuum flow, which occurs during the frozen startup of high-temperature heat pipes. Multiple startup models and non-condensable gas models were tested against experimental data for sodium heat pipes, showing excellent agreement. Additionally, the newly developed gas mixture model for modeling non-condensable gas is further tested with a theoretical case study with arbitrary heating configurations. Finally, several recommendations and conclusions are made from the studies in this work to guide future heat pipe modeling efforts.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Status of the OECD/NEA Watts Bar unit 1 benchmark and calculations of local reactor core power data of the Zero Power physics tests

The paper aims to provide an update on the status of the Organization for Economic Cooperation and Development (OECD) / the Nuclear Energy Agency (NEA) Tennessee Valley Authority (TVA) Watts Bar 1 (WB1) Multi-Physics Multi-Cycle depletion benchmark and the results associated with Exercise 1. The benchmark relies on the set of benchmark progression problems developed by the Department of Energy (DOE) Consortium for the Advanced Simulation of Light Water Reactors (CASL) for the Virtual Environment for Reactor (VERA), which are based on real plant design and operational data. The OECD/NEA TVA WB1 benchmark is designed for validation of both traditional and novel high-fidelity multi-physics codes to analyze Pressurized Water Reactors (PWR) depletion cycles. The activities are conducted under the Expert Group on Reactor Systems Multi-Physics (EGMUP) at NEA/OECD. In this work, we analyze Exercise 1 of the benchmark: stand-alone Three-Dimensional (3D) neutronics at Start-up Zero Power Physics Test (ZPPT) at Hot Zero Power Conditions (HZP) and the resulting power maps. The code used to model the exercise is the continuous energy Monte Carlo code Serpent 2.1.31 along with the nuclear data library ENDF/B-VII.1. Serpent results are compared with the high-fidelity deterministic code MPACT, which is part of VERA. The paper presents three selected results from power calculations required output. The results compared are the normalized axially integrated radial core power maps, the normalized axial averaged core power shapes, and the normalized core hottest and coldest assemblies radial power maps. The Root Mean Square Deviation (RMSD) between Serpent and MPACT is 0.51 % for the axially integrated radial core power map, 1.41 % for the axial averaged core power shape, 2.42 % for the hottest assembly axial power shape, 0.88 % for the coldest assembly axial integrated hottest assembly radial pin power map, and 0.17 % for the axially integrated coldest assembly radial pin power map. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

RMACXX: An Efficient High-Level C++ Interface over MPI-3 RMA

Parallel scientific applications can benefit from de- coupling communication and synchronization. One-sided pro- gramming abstractions, which separate communication from syn- chronization, have in fact served as a motivation for partitioned global address space (PGAS) models. However, the use of PGAS models in application codes in a manner that fully exploits the benefit of these programming models requires significant development effort. Meanwhile, a vast majority of scientific codes already use the Message Passing Interface (MPI) and need convenient features to support application-specific one-sided communication scenarios. MPI Remote Memory Access (RMA) can be employed for this purpose. MPI is a low-level API, however, and developing applications with MPI RMA requires programmers to be well versed in its nuances. We present RMACXX, a compact set of C++ bindings to MPI-3 RMA, to ease the use of MPI RMA. Unlike other PGAS models, which may have interoperability issues with MPI, RMACXX is written on top of MPI and uses the same runtime as MPI. The basic functionality of RMACXX adds only a relatively small number of extra instructions (about 20) to the critical communication path. Moreover, RMACXX provides an intuitive API for building a wide variety of scientific applications while enjoying performance matching handwritten MPI-3 RMA codes.

MPI-3 RMA, one-sided communication, PGAS, C++, exp↗