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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

High-Fidelity Simulation Aerodynamics Dataset of NACA 0012 and 0021 Airfoils

The repository contains time series data of pressure, viscous and moment forces for the NACA 0012 and 0021 airfoils at angles of attack (AOA), 5, 17, 30, 45, 60 and 90 degrees. The data was generated during the high-fidelity simulations of the airfoils using the open-source CFD code, Nalu-Wind (https://github.com/Exawind/nalu-wind). The research objective was to study the effects of mesh resolution and turbulence model on the three-dimensional deep stall aerodynamics of the airfoils, which has been published in the Journal of Turbulence (https://doi.org/10.1080/14685248.2023.2225141). Several meshes of the O-grid type were generated using the modified Pointwise Glyph script of Carrigan (https://github.com/pointwise/AirfoilMesh) by varying the wall normal and spanwise resolutions. Note that, the geometry was extruded by four chord lengths in the spanwise direction. In the paper, we showed that Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS/LES turbulence model with a minimum of 24 cells per chord length in the spanwise direction is necessary to correctly predict the loads in the deep stall regime. Under the NACA 0012 airfoil directory, the data for several combinations of wall-normal and spanwise resolutions are presented in each of the AOA subdirectories. The corresponding experimental data for the lift and drag polars can be found in the "exp_data" subdirectory. The directory for the NACA 0021 airfoil consists of CFD force data and airfoil polars for Reynolds numbers 2.7e5 and 2.0e6.

17 WIND ENERGY↗
DOE OSTI2019

Improvements to PVWatts for Fixed and One-Axis Tracking Systems

This work presents improvements to the widely used NREL PVWatts photovoltaic system energy model to improve modeling accuracy for typical fixed and one axis system designs. These improvements allow users to more credibly and quickly evaluate competing system designs in early stage feasibility, and are shown to improve PVWatts' system performance prediction capabilities without major impact to ease of use.

14 SOLAR ENERGY↗
DOE OSTI2018

Multiphase Large Eddy Simulation of a Pulsed Sieve-Plate Extraction Column - 18125

Presented here is research resulting from a framework of single-phase and multiphase modelling trials. These are intended to accurately simulate the internal flow phenomena arising in pulsed sieve-plate extraction columns representative of those found in industry. Three-dimensional unsteady flow calculations using large eddy simulation coupled with dynamic sub-grid scale modelling have been used to capture the transient effects of sinusoidal flow conditions using the open source CFD code OpenFOAM{sup R}. Multiphase interface interactions are modelled using the volume of fluid method with appropriate heavy and light phase fluidic properties capturing surface tension, density and viscous effects. The system studied consists of a 150 mm diameter column with two perforated plates with a fractional free area of 25 %, operational under 1 Hz pulsing frequency at a 50 mm amplitude and plate spacing of 300 mm. The materials simulated include water (single-phase) and 3 M nitric acid and dodecane / tributyl phosphate mixture at 30 vol% (multiphase). Results are compared with appropriate correlations from literature. The dispersive mixing efficiency of the system was also evaluated via a mixing index value which allows quantification of the ratio of rotational and irrotational flow components. This study has found that the predictive methods employed are effective at capturing the different scales of turbulence present within pulsed sieve-plate extraction columns, and its influence on their operation. As such, the work highlights the ability of large eddy simulation-based modelling to predict the complex flows in such column geometries. Furthermore, the techniques employed provide valuable insight into the hydrodynamic flow characteristics and mixing mechanisms arising within pulsed sieve-plate extraction columns. (authors)

42 ENGINEERING↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

Hakobyan, Hayk [Flatiron Institute, New York, NY (↗
DOE OSTI2015

Enhancements to the MCNP6 background source

The particle transport code MCNP has been used to produce a background radiation data file on a worldwide grid that can easily be sampled as a source in the code. Location-dependent cosmic showers were modeled by Monte Carlo methods to produce the resulting neutron and photon background flux at 2054 locations around Earth. An improved galactic-cosmic-ray feature was used to model the source term as well as data from multiple sources to model the transport environment through atmosphere, soil, and seawater. A new elevation scaling feature was also added to the code to increase the accuracy of the cosmic neutron background for user locations with off-grid elevations. Furthermore, benchmarking has shown the neutron integral flux values to be within experimental error.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Proof of Concept Coded Aperture Miniature Mass Spectrometer Using a Cycloidal Sector Mass Analyzer, a Carbon Nanotube (CNT) Field Emission Electron Ionization Source, and an Array Detector

Despite many potential applications, miniature mass spectrometers have had limited adoption in the field due to the tradeoff between throughput and resolution that limits their performance relative to laboratory instruments. Recently, a solution to this tradeoff has been demonstrated by using spatially coded apertures in magnetic sector mass spectrometers, enabling throughput and signal-to-background improvements of greater than an order of magnitude with no loss of resolution. This paper describes a proof of concept demonstration of a cycloidal coded aperture miniature mass spectrometer (C-CAMMS) demonstrating use of spatially coded apertures in a cycloidal sector mass analyzer for the first time. C-CAMMS also incorporates a miniature carbon nanotube (CNT) field emission electron ionization source and a capacitive transimpedance amplifier (CTIA) ion array detector. Results confirm the cycloidal mass analyzer’s compatibility with aperture coding. A >10× increase in throughput was achieved without loss of resolution compared with a single slit instrument. Several areas where additional improvement can be realized are identified. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Automatic translation of MPI source into a latency-tolerant, data-driven form

Hiding communication behind useful computation is an important performance programming technique but remains an inscrutable programming exercise even for the expert. We present Bamboo, a code transformation framework that can realize communication overlap in applications written in MPI without the need to intrusively modify the source code. Bamboo reformulates MPI source into the form of a task dependency graph that expresses a partial ordering among tasks, enabling the program to execute in a data-driven fashion under the control of an external runtime system. Experimental results demonstrate that Bamboo significantly reduces communication delays while requiring only modest amounts of programmer annotation for a variety of applications and platforms, including those employing co-processors and accelerators. Moreover, Bamboo's performance meets or exceeds that of labor-intensive hand coding. The translator is more than a means of hiding communication costs automatically; it demonstrates the utility of semantic level optimization against a wellknown library.

Communication-computation↗
DOE OSTI2014

SCALE Enhancements for Detailed Cask Dose Rate Analysis - Paper 37

The Nuclear Fuels Storage and Transportation Planning Project and the Used Fuel Disposition Campaign, both under the US Department of Energy (DOE) Office of Nuclear Energy, require the ability to predict k{sub eff}, heat generation and temperature distributions, and external dose rates from any combination of commercial spent nuclear fuel assemblies loaded into any type of licensed cask. Work in the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory has begun to develop the UNF-ST[and]DARDS system (Used Nuclear Fuel Storage, Transportation and Disposal Analysis Resource and Data System) that uses many of the nuclear analysis codes from the SCALE package, a comprehensive database of used nuclear fuel inventories from utilities, and a library of model template files to produce complete input files for the various types of UNF computational analyses. UNF-ST[and]DARDS uses the template files and input parameters that are applicable to a safety analysis to assemble and execute a complete input file, collect the output files and relevant information for subsequent calculations, and illustrate calculation results as a function of decay time and other relevant analysis parameters. This paper will focus on the shielding aspects and the ability to predict external dose rates from as-loaded casks. Realistic dose rate values based on actual canister fuel contents may be used in developing appropriate plans, procedures, and methods to keep personnel radiation exposures as low as reasonably achievable during cask repackaging and transportation operations, as well as in performing total system analyses and transportation risk evaluations. The MAVRIC shielding sequence of the SCALE package was created for the US Nuclear Regulatory Commission for analyzing dose rates outside of storage, transfer, or transportation casks loaded with spent nuclear fuel. MAVRIC contains automated variance reduction capabilities based on the CADIS and FW-CADIS methods that are effective for optimizing deep penetration Monte Carlo shielding problems for achieving low relative uncertainties in short computing times. MAVRIC uses the Denovo discrete ordinates code to compute importance maps and biased source terms for the Monaco fixed-source Monte Carlo code, which can use multi-group or continuous-energy cross section data. ORIGEN (Oak Ridge Isotope Generation and Depletion) is the SCALE computer code for calculating time-dependent isotopic concentrations during irradiation and decay for more than 2,200 nuclides that can be produced by nuclear fuel irradiation and activation. A series of SCALE enhancements was implemented to enable automatic decay heat and radiation source term generation and shielding safety analysis of as-loaded casks. The enhancements to SCALE made under this project will be a part of the SCALE 6.2 release expected in fall 2014. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Property of neoclassical GAMs induced by pellet generated plasma perturbations in the gyrokinetic code XGC

A pellet source model has been implemented into the gyrokinetic code XGC and applied to neoclassical geodesic acoustic mode (GAM) study as the first step to investigate potentially rapid kinetic spread of pellet-born particles following their injection into plasma. Here in this study, GAM oscillations of the radial electric field, as well as of their effect on the radial particle and energy fluxes, are studied and the observed frequencies agree with theoretical expectations. It is verified that (i) GAM oscillations are driven both within the pellet-fuelled region and for locations at a larger radius than this region and not for locations radially inside the pellet-fuelled region, and that, (ii) even though the poloidal spread along the magnetic field lines is fast, the time-averaged radial plasma transport under GAM oscillations is kept to roughly the same level as the neoclassical transport in the absence of plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Parser for Sabin-to-Mahoney Transition Model of Quasispecies Replication

This code is a data parse for preparing output from the Qspp agent-based stochastic simulation model for plotting in Excel. This code is specific to a set of simulations that were run for the purpose of preparing data for a publication. It is necessary to make this code open-source in order to publish the model code (Qspp), which has already been released. There is a necessity of assuring that results from using Qspp for a publication

Ecale Zhou, Carol↗
DOE OSTI2016

Parser for Sabin-to-Mahoney Transition Model of Quasispecies Replication

This code is a data parse for preparing output from the Qspp agent-based stochastic simulation model for plotting in Excel. This code is specific to a set of simulations that were run for the purpose of preparing data for a publication. It is necessary to make this code open-source in order to publish the model code (Qspp), which has already been released. There is a necessity of assuring that results from using Qspp for a publication

Ecale Zhou, Carol↗
DOE OSTI2015

Development of the computer code for multi-calculations of source term of spent fuel - 15351

Many PWR spent fuels have been transferred between the neighboring nuclear power plants at Kori, Hanbit, Hanwul site, and from nuclear power plant to the KAERI(Korea Atomic Energy Research Institute) for the research purpose in Korea. In case of the PHWR type spent fuels, these have been transferred from the spent fuel pool to an on-site dry storage facility. Total radioactivity of spent fuels has to be calculated before transport in advance by the regulation requirement. In case of PWR, KN-12 transport cask can transfer 12 spent fuel assemblies at one time and in case of PHWR, Hi-STAR 63 transport cask can transfer 120 spent fuel assemblies at one time. Therefore, total radioactivity of all transported spent fuel assemblies that are transported generally are calculated by using the ORIGEN-ARP to meet this regulatory requirement. Basically, ORIGEN-ARP calculates the one fuel assembly at one time. SCALE code including this code provides the batch mode to perform the multi-calculations of the source term. In this study, in order to use the batch mode of ORIGEN-ARP, ARP-auto program was developed to generate input files and batch file of ORIGEN-ARP automatically and to make the summary file extracting the result from output files of ORIGEN-ARP. This code enhances the efficiency at the field. The working time to complete the calculation of source terms of 120 spent fuel assemblies take only 20 minutes. In order to confirm the credibility of this code, total radioactivity of 12 spent fuel assemblies of 17x17 Westinghouse type were calculated and compared to the previous result. The initial enrichments of these spent fuels are 3.46 ∼ 3.50 w/o and burnups are 29 ∼ 43 GWD/MTU. As the results, the relative errors are about 0.9 ∼ 3.0%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Data and Scripts associated with “Lambda-PFLOTRAN: Workflow for Incorporating Organic Matter Chemistry Informed by Ultra High Resolution Mass Spectrometry into Biogeochemical Modeling.”

This data package is associated with the publication “Lambda-PFLOTRAN: Workflow for Incorporating Organic Matter Chemistry Informed by Ultra High Resolution Mass Spectrometry into Biogeochemical Modeling” submitted to Geoscientific Model Development (Muller et al., 2024). In this manuscript, organic matter chemistry and thermodynamics are directly connected to reactive transport simulators through the newly developed Lambda-PFLOTRAN (Parallel Reactive Flow and Transport model) workflow tool that succinctly incorporates organic matter chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate aerobic respiration of the organic matter and the resulting biogeochemistry. Lambda-PFLOTRAN is a python-based workflow, executed through a Jupyter Notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow, and transport code PFLOTRAN. This data package contains Jupyter Notebook based workflows for two test cases for running biogeochemical simulations of organic matter oxidation identified by FTICR-MS. It contains four primary folders (workflow, data, src, and analysis), a file-level metadata file (Muller_2024_Lambda_PFLOTRAN_Manuscript_Data_Package_flmd.csv) that lists all the files contained in this data package with a short description of each, and a data dictionary (Muller_2024_Lambda_PFLOTRAN_Manuscript_Data_Package_dd.csv) file that describes the tabular column headers. The ‘workflow’ folder contains the Jupyter Notebook based workflows for running the lambda analysis, PFLOTRAN simulation, sensitivity analysis and parameter estimation. The ‘data’ folder contains the FTICR-MS data, initial conditions, and incubation data for test cases 1 and 2 in folders titled ‘WHONDRS’ and ‘Colloids’, respectively. The data folder also has a ‘Database’ folder containing a reaction network for bulk organic matter (assumed to be CH2O) and a general database for PFLOTRAN (hanford_rxn_network). The CH2O reaction network defines bulk organic matter oxidation. Biogeochemical simulations are completed for both the lambda binned organic matter and bulk organic matter reaction networks. The ‘hanford_rxn_network’ database includes information required for PFLTORAN simulations including ion size, molar mass, and charge of the aqueous species, gases, and minerals phases. The ‘src’ folder contains python source codes for performing lambda analysis, PFLOTRAN simulation, sensitivity analysis and parameter estimation. The ‘analysis’ folder contains outputs from the test cases 1 and 2 including lambda analysis, PFLOTRAN runs and the calibration results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD

We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD.The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters.Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances:(i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community;(ii) the AD tool used, OpenAD, is an open-source tool;(iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations.A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large.As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet.Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that surface temperature sensitivities are higher broadly across the floating ice shelves, while basal temperature sensitivities are highest at the grounding lines of floating ice shelves and outlet glaciers. A uniformly perturbed region of East Antarctica reveals that, among the four control variables tested here, total ice sheet volume is the most sensitive to variations in austral summer precipitation as formulated in SICOPOLIS.Comparison between adjoint- and finite-difference-derived sensitivities shows good agreement, lending confidence that the AD tool is producing correct adjoint code.The new modeling infrastructure is freely available at http://www.sicopolis.net (last access: 2 April 2020) under the development trunk.

58 GEOSCIENCES↗