Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source 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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Radiation Source Localization Using Surrogate Models Constructed from 3-D Monte Carlo Transport Physics Simulations

Recent research has focused on the development of surrogate models for radiation source localization in a simulated urban domain. We employ the Monte Carlo N-Particle (MCNP) code to provide high- delity simulations of radiation transport within an urban domain. The model is constructed to employ a source location (x, y, z) as input and return the estimated count rate for a set of speci ed detector locations. Because MCNP simulations are computationally expensive, we develop e cient and accurate surrogate models of the detector responses. We construct surrogate models using Gaussian processes (GP) and neural networks (NN) that we train and verify using the MCNP simulations. The trained surrogate models provide an e cient framework for Bayesian inference and experimental design. We employ Delayed Rejection Adaptive Metropolis (DRAM), a Markov Chain Monte Carlo (MCMC) algorithm, to infer the location and intensity of an unknown source. The DRAM results yield a posterior probability distribution for the source's location conditioned on the observed detector count rates. The posterior distribution exhibits regions of high and low probability within the simulated environment identifying potential source locations. In this manner, we can quantify the source location to within at least one of these regions of high probability in the considered cases. Employing these methods, we are able to reduce the space of potential source locations by at least 60%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

EAM-X Potential Generator

SAND2023-06608O The EAM-X Potential Generator code establishes parameters that describe interactions between atoms in pure metals and alloys. The open-source software allows for the understanding of the relationship between general properties of atoms and macroscopic properties of metals. The code reads in values for the parameters of the model and outputs a table that describes the interactions. This code will be published by Acta Materialia along with two papers describing the method. The output of this code is meant for use by the LAMMPS simulation code created by Sandia. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Chandross, Michael↗
DOE OSTI2015

A new line-of-sight approach to the non-linear Cosmic Microwave Background

We develop the transport operator formalism, a new line-of-sight integration framework to calculate the anisotropies of the Cosmic Microwave Background (CMB) at the linear and non-linear level. This formalism utilises a transformation operator that removes all inhomogeneous propagation effects acting on the photon distribution function, thus achieving a split between perturbative collisional effects at recombination and non-perturbative line-of-sight effects at later times. The former can be computed in the framework of standard cosmological perturbation theory with a second-order Boltzmann code such as SONG, while the latter can be treated within a separate perturbative scheme allowing the use of non-linear Newtonian potentials. We thus provide a consistent framework to compute all physical effects contained in the Boltzmann equation and to combine the standard remapping approach with Boltzmann codes at any order in perturbation theory, without assuming that all sources are localised at recombination.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Impurity mixing and radiation asymmetry in massive gas injection simulations of DIII-D

Simulations of neon massive gas injection into DIII-D are performed with the 3D MHD code NIMROD. The poloidal and toroidal distribution of the impurity source is varied. This report will focus on the effects of the source variation on impurity mixing and radiated power asymmetry. Even toroidally symmetric impurity injection is found to produce asymmetric radiated power due to asymmetric convective heat flux produced by the 1/1 mode. When the gas source is toroidally localized, the phase relationship between the mode and the source location is important, affecting both radiation peaking and impurity mixing. Under certain circumstances, a single, localized gas jet could produce better radiation symmetry during the disruption thermal quench than evenly distributed impurities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Using Griffin's Transmutation Solver to Calculate Radiation Damage

Displacement Radiation Damage originates from all nuclides, not just those that are naturally occurring. Currently only damage from naturally-occurring nuclides, or sometimes damage from one transmutation product is considered. It is proposed that the transmutation solvers implemented in many codes be used to calculate this radiation damage. This would explicitly treat damage from all sources without additionally burdening the user. A proof-of-concept implementation was created in Griffin. The implementation showed that only minor code modifications are necessary to add this feature. When compared against an analytical benchmark the results Griffin could calculate were very accurate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications

Here, the analysis of many fusion applications such as liquid-metal blankets requires application of computational fluid dynamics (CFD) methods for electrically conductive liquids in geometrically complex regions and in the presence of a strong magnetic field. A current state of the art general purpose CFD code allows modeling of the flow in complex geometric regions, with simultaneous conjugated heat transfer analysis in liquid and surrounding solid parts. Together with a magnetohydrodynamics (MHD) capability, the general purpose CFD code will be a valuable tool for the design and optimization of fusion devices. This paper describes an introduction of MHD capability into the general purpose CFD code CFX, part of the ANSYS Workbench. The code was adapted for MHD problems using a magnetic induction approach. CFX allows introduction of user-defined variables using transport or Poisson equations. For MHD adaptation of the code three additional transport equations were introduced for the components of the magnetic field, in addition to the Poisson equation for electric potential. The Lorentz force is included in the momentum transport equation as a source term. Fusion applications usually involve very strong magnetic fields, with values of the Hartmann number of up to tens of thousands. In this situation a system of MHD equations become very rigid with very large source terms and very strong variable gradients. To increase system robustness, special measures were introduced during the iterative convergence process, such as linearization using source coefficient for momentum equations. The MHD implementation in general purpose CFD code was tested against benchmarks, specifically selected for liquid-metal blanket applications. Results of numerical simulations using present implementation closely match analytical solutions for a Hartmann number of up to 1500 for a 2-D laminar flow in the duct of square cross section, with conducting and nonconducting walls. Results for a 3-D test case are also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

SuperNu Version 4.x

We seek to release SuperNu, Version 4.x, as a continuation of development for the open source SuperNu software. The SuperNu, Version 3.x Monte Carlo radiative transfer code for astrophysical transients was released with GPLv3 copyright, asserted by LANL in 2015. For the next release we have features planned for development , including: opacity implementation (including non-local thermodynamic equilibrium effects), generalized source implementation (e.g. for emulating shock heating in Type II supernovae), 3T (electron, ion, radiation) internal energy update, special relativity corrections through O(v^2/c^2), light polarization (e.g. for comparison to spectroplarimetry observations of supernovae and kilonovae), and infrastructure features (checkpoint and restart of simulations, and tools including setup and Slurm batch scripts and simulation post-processing/analysis scripts). These features are intended to improve the fidelity and/or better understand uncertainty in supernova and kilonova light curve calculations.

Wollaeger, Ryan↗
DOE OSTI2015

COG11.1 description, new features, and development activities - 14521

The latest COG code version, COG11.1, is now available. New features include: hybrid criticality detector variance reduction (CritDetVR); delayed fission gammas (DFG); nuclear resonance fluorescence (NRF); radiation sources (RadSrc); an updated COG dictionary (COGLEX); and additional nuclear data libraries. Ongoing software quality assurance and code development activities are also described. (authors)

42 ENGINEERING↗
DOE OSTI2022

Multiphysics analysis system for heat pipe cooled micro-reactors employing PRAGMA-OpenFOAM-ANLHTP

A multiphysics analysis system for neutronics/thermo-mechanical/heat-pipe analysis of heat pipe cooled micro-reactors was developed using the PRAGMA code as the neutronics engine. PRAGMA, which used to be a GPU-based continuous-energy MC code for power reactor applications, now has an extended geometry package to handle geometries with unstructured meshes generated by the ANSYS Design-Modeler and Meshing. The NVIDIA ray tracing engine OptiX was exploited for efficient neutron transport on unstructured geometry. On the multiphysics side, the open-source CFD tool OpenFOAM and one-dimensional heat pipe analysis code ANLHTP were adopted. The manager-worker system based on the MPI dynamic process management (DPM) model enables efficient coupling of codes employing different parallelization schemes. With all the features, the multiphysics analysis of the one-sixth symmetrical MegaPower 2D core was performed and it demonstrated the benefit of the tight integration of the three-way coupling system and one-to-one geometry coupling strategy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

GAMSOR: Gamma Source Preparation and DIF3D Flux Solution

Nuclear reactors that rely upon the fission reaction have two modes of thermal energy deposition in the reactor system: neutron absorption and gamma absorption. The gamma rays are typically generated by neutron absorption reactions or during the fission process which means the primary driver of energy production is of course the neutron interaction. In conventional reactor physics methods, the gamma heating component is ignored such that the gamma absorption is forced to occur at the gamma emission site. For experimental reactor systems like EBR-II and FFTF, the placement of structural pins and assemblies internal to the core leads to problems with power heating predictions because there is no fission power source internal to the assembly to dictate a spatial distribution of the power. As part of the EBR-II support work in the 1980s, the GAMSOR code was developed to assist analysts in calculating the gamma heating. The GAMSOR code is a modified version of DIF3D and actually functions within a sequence of DIF3D calculations. The gamma flux in a conventional fission reactor system does not perturb the neutron flux and thus the gamma flux calculation can be cast as a fixed source problem given a solution to the steady state neutron flux equation. This leads to a sequence of DIF3D calculations, called the GAMSOR sequence, which involves solving the neutron flux, then the gamma flux, then combining the results to do a summary edit. In this manuscript, we go over the GAMSOR code and detail how it is put together and functions. We also discuss how to setup the GAMSOR sequence and input for each DIF3D calculation in the GAMSOR sequence. With the GAMSOR capability, users can take any valid steady state DIF3D calculation and compute the power distribution due to neutron and gamma heating. The MC2-3 code is the preferable companion code to use for generating neutron and gamma cross section data, but the GAMSOR code can accept cross section data from other sources. To further this aspect, an additional utility code was created which demonstrates how to merge the neutron and gamma cross section data together to carry out a simultaneous solve of the two systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. II. General Relativistic Module

Black hole (BH) environments often host plasmas that are fully collisionless or contain intrinsically collisionless regions, including relativistic jets and coronae, where particle energization is ubiquitous. Capturing the physics of these systems requires numerical methods capable of modeling relativistic, magnetized, collisionless plasmas in strong gravitational fields. In this work, we introduce the general relativistic module for Entity—the first open-source, coordinate-agnostic, performance-portable, particle-in-cell code. The code enables fast axisymmetric simulations of collisionless plasmas around BHs on any modern high-performance computing architecture (both GPUs and CPUs).

Galishnikova, Alisa [Flatiron Institute, New York,↗
DOE OSTI2014

Directional gamma sensing from covariance processing of inter-detector Compton crosstalk energy asymmetries

Energy asymmetry of inter-detector crosstalk from Compton scattering can be exploited to infer the direction to a gamma source. A covariance approach extracts the correlated crosstalk from data streams to estimate matched signals from Compton gammas split over two detectors. On a covariance map the signal appears as an asymmetric cross diagonal band with axes intercepts at the full photo-peak energy of the original gamma. The asymmetry of the crosstalk band can be processed to determine the direction to the radiation source. The technique does not require detector shadowing, masking, or coded apertures, thus sensitivity is not sacrificed to obtain the directional information. An angular precision of better than 1° of arc is possible, and processing of data streams can be done in real time with very modest computing hardware.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Fast and accurate CMB computations in non-flat FLRW universes

We present a new method for calculating CMB anisotropies in a non-flat Friedmann universe, relying on a very stable algorithm for the calculation of hyperspherical Bessel functions, that can be pushed to arbitrary precision levels. We also introduce a new approximation scheme which gradually takes over in the flat space limit and leads to significant reductions of the computation time. Our method is implemented in the Boltzmann code class. It can be used to benchmark the accuracy of the camb code in curved space, which is found to match expectations. For default precision settings, corresponding to 0.1% for scalar temperature spectra and 0.2% for scalar polarisation spectra, our code is two to three times faster, depending on curvature. We also simplify the temperature and polarisation source terms significantly, so the different contributions to the C{sub ℓ} 's are easy to identify inside the code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Using Python For 3D Visualizations of MCNP output

MCNP is ubiquitously used within industries, universities and national laboratories as a tool for accurately modeling radiation. With so much focus on the details of transport, attention is seldom given to the visualizations of these meticulously constructed outputs. Python is a programing language designed to promote code readability and has steadily grown in popularity since the 1990's. Python is also a heavily crowd-sourced programing language which generally means that the problems that you are trying to solve with computer codes are already solved; The user is often left with the now-simplified task of cobbling together a good and useful compilation of these tools and adapting them to solve the problem at hand. I have taken advantage of this situation and constructed some visualization tools in python that can be used to duplicate some of the 2D radiation fields that MCNP produces and also produce a meaningful 3D representation of a radiation field. (author)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Two-Scale 13 C Metabolic Flux Analysis for Metabolic Engineering

Accelerating the Design-Build-Test-Learn (DBTL) cycle in synthetic biology is critical to achieving rapid and facile bioengineering of organisms for the production of, e.g., biofuels and other chemicals. The Learn phase involves using data obtained from the Test phase to inform the next Design phase. As part of the Learn phase, mathematical models of metabolic fluxes give a mechanistic level of comprehension to cellular metabolism, isolating the principle drivers of metabolic behavior from the peripheral ones, and directing future experimental designs and engineering methodologies. Furthermore, the measurement of intracellular metabolic fluxes is specifically noteworthy as providing a rapid and easy-to-understand picture of how carbon and energy flow throughout the cell. Here, we present a detailed guide to performing metabolic flux analysis in the Learn phase of the DBTL cycle, where we show how one can take the isotope labeling data from a 13 C labeling experiment and immediately turn it into a determination of cellular fluxes that points in the direction of genetic engineering strategies that will advance the metabolic engineering process.For our modeling purposes we use the Joint BioEnergy Institute (JBEI) Quantitative Metabolic Modeling (jQMM) library, which provides an open-source, python-based framework for modeling internal metabolic fluxes and making actionable predictions on how to modify cellular metabolism for specific bioengineering goals. It presents a complete toolbox for performing different types of flux analysis such as Flux Balance Analysis, 13 C Metabolic Flux Analysis, and it introduces the capability to use 13 C labeling experimental data to constrain comprehensive genome-scale models through a technique called two-scale 13 C Metabolic Flux Analysis (2S- 13 C MFA) [1]. In addition to several other capabilities, the jQMM is also able to predict the effects of knockouts using the MoMA and ROOM methodologies. The use of the jQMM library is illustrated through a step-by-step demonstration, which is also contained in a digital Jupyter Notebook format that enhances reproducibility and provides the capability to be adopted to the user's specific needs. As an open-source software project, users can modify and extend the code base and make improvements at will, providing a base for future modeling efforts.

13C Metabolic flux analysis↗
DOE OSTI2018

Genome-Scale 13 C Fluxomics Modeling for Metabolic Engineering of Saccharomyces cerevisiae .

Synthetic biology is a rapidly developing field that pursues the application of engineering principles and development approaches to biological engineering. Synthetic biology is poised to change the way biology is practiced, and has important practical applications: for example, building genetically engineered organisms to produce biofuels, medicines, and other chemicals. Traditionally, synthetic biology has focused on manipulating a few genes (e.g., in a single pathway or genetic circuit), but its combination with systems biology holds the promise of creating new cellular architectures and constructing complex biological systems from the ground up. Enabling this merge of synthetic and systems biology will require greater predictive capability for modeling the behavior of cellular systems, and more comprehensive data sets for building and calibrating these models. The so-called "-omics" data sets can now be generated via high throughput techniques in the form of genomic, proteomic, transcriptomic, and metabolomic information on the engineered biological system. Of particular interest with respect to the engineering of microbes capable of producing biofuels and other chemicals economically and at scale are metabolomic datasets, and their insights into intracellular metabolic fluxes. Metabolic fluxes provide a rapid and easy to understand picture of how carbon and energy flow throughout the cell. Furthermore, we present a detailed guide to performing metabolic flux analysis and modeling using the open source JBEI Quantitative Metabolic Modeling (jQMM) library. This library allows the user to transform metabolomics data in the form of isotope labeling data from a 13 C labeling experiment into a determination of cellular fluxes that can be used to develop genetic engineering strategies for metabolic engineering.The jQMM library presents a complete toolbox for performing a range of different tasks of interest in metabolic engineering. Various different types of flux analysis and modeling can be performed such as flux balance analysis, 13 C metabolic flux analysis, and two-scale 13 C metabolic flux analysis (2S- 13 C MFA). 2S- 13 C MFA is a novel method that determines genome-scale fluxes without the need of every single carbon transition in the metabolic network. In addition to several other capabilities, the jQMM library can make model based predictions for how various genetic engineering strategies can be incorporated toward bioengineering goals: it can predict the effects of reaction knockouts on metabolism using both the MoMA and ROOM methodologies. In this chapter, we will illustrate the use of the jQMM library through a step-by-step demonstration of flux determination and knockout prediction in a complex eukaryotic model organism: Saccharomyces cerevisiae ( S. cerevisiae ). Included with this chapter is a digital Jupyter Notebook file that provides a computable appendix showing a self-contained example of jQMM usage, which can be changed to fit the user's specific needs. As an open source software project, users can modify and extend the code base to make improvements at will, allowing them to share their development work and contribute back to the jQMM modeling community.

-Omics data↗
DOE OSTI2015

Validation of Monte Carlo simulation of neutron production in a spallation experiment

A renewed interest in experimental research on Accelerator-Driven Systems (ADS) has been initiated by the global attempt to produce energy from thorium as a safe(r), clean(er) and (more) proliferation-resistant alternative to the uranium-fuelled thermal nuclear reactors. The ADS research has been actively pursued at the Joint Institute for Nuclear Research (JINR), Dubna, since decades. Most recently, the emission of fast neutrons was experimentally investigated at the massive (m = 512 kg) natural uranium spallation target QUINTA. The target has been irradiated with the relativistic deuteron beams of energy from 0.5 AGeV up to 4 AGeV at the JINR Nuclotron accelerator in numerous experiments since 2011. Neutron production inside the target was studied through the gamma-ray spectrometry measurement of natural uranium activation detectors. Experimental reaction rates for (n,γ), (n,f) and (n,2n) reactions in uranium have provided valuable information about the neutron distribution over a wide range of energies up to some GeV. The experimental data were compared to the predictions of Monte Carlo simulations using the MCNPX 2.7.0 code. In conclusion, the results are presented and potential sources of partial disagreement are discussed later in this work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Benchmark of the KGMf with a coupled Boltzmann equation solver

The Kinetic Global Model framework (KGMf) is an open-source general-purpose global model (spatially averaged) simulation code developed to explore the reaction kinetics and pathways in plasma discharge systems. It contains species continuity and electron energy balance equations, with a time-dependent evaluated electron energy distribution function (EEDF) for electron impact reactions. The EEDF is utilized to determine the rate coefficients for electron impact reactions, which can have profound impact on the temporal evolutions of plasma parameters. Previously, the EEDF was commonly assumed as Maxwellian or an analytical function of the effective electron temperature. In this work, the KGMf is coupled with a Boltzmann equation (BE) solver to self-consistently compute the EEDF. The EEDF evolution frequency is determined based on relative changes of the reduced electric field. The KGMf is benchmarked with the ZDPlasKin code based on high-pressure low-temperature argon plasma discharge cases. Additionally, the temporal evolutions of reduced electric field, electron temperature, EEDF, reaction rates, and species densities, are obtained and compared under different discharge conditions, showing good agreement between the KGMf and the ZDPlasKin simulations. The application of the KGMf for predicting breakdown times in high power microwave discharges is also presented, which shows qualitative agreement with particle-in-cell simulations. The KGMf can be further applied for more complicated plasma discharge systems, where the reaction kinetics are more intricate (e.g., plasma-assisted combustion systems).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗