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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Rigorous Two-Step Activation for Fusion Systems with PyNE - Paper 127

In fusion systems, high energy neutrons activate system components. After shutdown, the biological dose rate from photons emitted from these components must be estimated as a function of position and time for maintenance planning. High resolution results can be obtained using the mesh-based Rigorous Two-Step (R2S) method. A neutron transport simulation provides mesh-based neutron fluxes. Using these fluxes and an appropriate irradiation scenario, a nuclear inventory analysis code is used to determine the photon source density within each mesh volume element at any decay time after shutdown. The photon source density mesh is used as a source for photon transport. Photon fluxes are tallied at positions of interest and flux-to-dose rate conversion factors are used to estimate the biological dose rate. Previous work by the authors involved the creation of R2S-ACT, an open source R2S workflow, using the Direct Accelerated Geometry Monte Carlo (DAGMC) version of MCNP5 (DAG-MCNP5) for both transport steps and the ALARA nuclear inventory analysis code. Various use cases may require different combinations of transport and activation codes, and nearly all of the capabilities of R2S-ACT have broader utility. For these reasons, the individual components of R2S-ACT were reimplemented within PyNE. PyNE is a collaborative, trans-institutional, open source project consisting of Python, Cython, and C++ code pertinent to nuclear simulations and analysis. This work involved additions to the PyNE mcnp, mesh, material, alara, dagmc, and sampling modules. A lightweight r2s module was then added to connect the code from these modules. The r2s and sampling modules, together with DAG-MCNP5 and ALARA form a complete R2S workflow for obtaining shutdown dose rate results. Mesh operations, geometry discretization, physics code input/output parsing, and mesh-based source sampling capabilities have been added to PyNE generically, eliminating the need to duplicate this effort across multiple applications. The PyNE r2s module uses these capabilities to couple MCNP and ALARA to form an R2S workflow, which has been validated with the FNG shutdown dose rate benchmark. Future work will involve benchmarking the workflow using tetrahedral mesh and also non-analog sampling options. added so that they can be used interchangeably. This will allow for comparisons between results obtained by different physics codes with minimal human effort. Work will also be done to create and implement an importance sampling method for biasing neutron transport towards regions of phase space of high importance to photon dose rate tallies. This will improve the efficiency of obtaining shutdown dose rate results in large, highly-attenuating geometries. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Code for the manuscript "Mori-Zwanzig Modal Decomposition"

We would like to create an open source repository in LANL's github on code written in Julia, in which we implement and extend the data-driven Mori-Zwanzig method for extracting large-scale spatio-temporal structures from data, which we call MZMD. This method is an extension of Dynamic Mode Decomposition (DMD) in which Mori-Zwanzig memory kernels are included into the associated companion matrix. In the code we would like to release, we apply MZMD to a flow over a cylinder with Reynolds number 100 rather than the much larger data set used in the associated manuscript. DMD is used extensively in the fluid dynamics community mainly for extracting large scale spatio-temporal structures (patters) from flow data. This is useful for understanding the key mechanisms that generate certain complex dynamical process relevant in engineering design. In MZMD, we improve upon DMD by adding the Mori-Zwanzig memory kernels, and show this improvement is especially important in strongly nonlinear regions of the flow.

Woodward, Michael↗
DOE OSTI2018

MCNP Version 6.2 Release Notes

Monte Carlo N-Particle or MCNP ® is a general-purpose Monte Carlo radiation-transport code designed to track many particle types over broad ranges of energies. This MCNP Version 6.2 follows the MCNP6.1.1 beta version and has been released in order to provide the radiation transport community with the latest feature developments and bug fixes for MCNP. Since the last release of MCNP major work has been conducted to improve the code base, add features, and provide tools to facilitate ease of use of MCNP version 6.2 as well as the analysis of results. These release notes serve as a general guide for the new/improved physics, source, data, tallies, unstructured mesh, code enhancements and tools. For more detailed information on each of the topics, please refer to the appropriate references or the user manual which can be found at http://mcnp.lanl.gov. This release of MCNP version 6.2 contains 39 new features in addition to 172 bug fixes and code enhancements. There are still some 33 known issues the user should familiarize themselves with (see Appendix).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

An Educational Guide for 2D Stellar Structure Calculations of Rapidly Rotating Stars using the ESTER code

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Measurement and analysis of GPU-accelerated applications with HPCToolkit

To address the challenge of performance analysis on the US DOE’s forthcoming exascale supercomputers, Rice University has been extending its HPCToolkit performance tools to support measurement and analysis of GPU-accelerated applications. To help developers understand the performance of accelerated applications as a whole, HPCToolkit’s measurement and analysis tools attribute metrics to calling contexts that span both CPUs and GPUs. To measure GPU-accelerated applications efficiently, HPCToolkit employs a novel wait-free data structure to coordinate monitoring and attribution of GPU performance. To help developers understand the performance of complex GPU code generated from high-level programming models, HPCToolkit constructs sophisticated approximations of call path profiles for GPU computations. To support fine-grained analysis and tuning, HPCToolkit uses PC sampling and instrumentation to measure and attribute GPU performance metrics to source lines, loops, and inlined code. To supplement fine-grained measurements, HPCToolkit can measure GPU kernel executions using hardware performance counters. To provide a view of how an execution evolves over time, HPCToolkit can collect, analyze, and visualize call path traces within and across nodes. Finally, on NVIDIA GPUs, HPCToolkit can derive and attribute a collection of useful performance metrics based on measurements using GPU PC samples. Here, we illustrate HPCToolkit’s new capabilities for analyzing GPU-accelerated applications with several codes developed as part of the Exascale Computing Project.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Lightweight runtime checking of C programs with RTC

The C Programming Language is known for being an efficient language that can be compiled on almost any architecture and operating system. However the absence of dynamic safety checks and a relatively weak type system allows programmer oversights that are hard to spot. In this paper, we present RTC, a runtime monitoring tool that instruments unsafe code and monitors the program execution. RTC is built on top of the ROSE compiler infrastructure. RTC finds memory bugs and arithmetic overflows and underflows, and run-time type violations. Most of the instrumentations are directly added to the source file and only require a minimal runtime system. As a result, the instrumented code remains portable. In tests against known error detection benchmarks, RTC found 98% of all memory related bugs and had zero false positives. Finally, in performance tests conducted with well known algorithms, such as binary search and MD5, we determined that our tool has an average run-time overhead rate of 9.7× and memory overhead rate of 3.5×.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Final Report for ALCC Allocation: Predictive Simulation of Complex Flow in Wind Farms

This report documents work performed using ALCC computing resources granted under a proposal submitted in February 2016, with the resource allocation period spanning the period July 2016 through June 2017. The award allocation was 10.7 million processor-hours at the National Energy Research Scientific Computing Center. The simulations performed were in support of two projects: the Atmosphere to Electrons (A2e) project, supported by the DOE EERE office; and the Exascale Computing Project (ECP), supported by the DOE Office of Science. The project team for both efforts consists of staff scientists and postdocs from Sandia National Laboratories and the National Renewable Energy Laboratory. At the heart of these projects is the open-source computational-fluid-dynamics (CFD) code, Nalu. Nalu solves the low-Mach-number Navier-Stokes equations using an unstructured- grid discretization. Nalu leverages the open-source Trilinos solver library and the Sierra Toolkit (STK) for parallelization and I/O. This report documents baseline computational performance of the Nalu code on problems of direct relevance to the wind plant physics application - namely, Large Eddy Simulation (LES) of an atmospheric boundary layer (ABL) flow and wall-modeled LES of a flow past a static wind turbine rotor blade. Parallel performance of Nalu and its constituent solver routines residing in the Trilinos library has been assessed previously under various campaigns. However, both Nalu and Trilinos have been, and remain, in active development and resources have not been available previously to rigorously track code performance over time. With the initiation of the ECP, it is important to establish and document baseline code performance on the problems of interest. This will allow the project team to identify and target any deficiencies in performance, as well as highlight any performance bottlenecks as we exercise the code on a greater variety of platforms and at larger scales. The current study is rather modest in scale, examining performance on problem sizes of O(100 million) elements and core counts up to 8k cores. This will be expanded as more computational resources become available to the projects.

17 WIND ENERGY↗
DOE OSTI2023

Code-verification techniques for the method-of-moments implementation of the magnetic-field integral equation

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementation of the magnetic-field integral equation, are frequently used to solve Maxwell's equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches for cases with and without coding errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Agile Software Components and Libraries for Large-Scale Reservoir Simulations

The objective of this project was for Sandia National Laboratories (SNL) and ExxonMobil Research Company (URC) to collaborate on the development and integration of numerical algorithms and software components for the solution of partial differential equations (PDEs) that address the high fidelity simulation of reservoir flows. SNL’s role was to consult and advise URC on the best practices and use of SNL opensource software components as part of URC's work to incorporate these components into their PDE-based URC reservoir simulation code. A desired outcome was for URC to leverage SNL’s open-source software components as within their own codes, rather than to have to develop from scratch those same capabilities themselves, and for SNL to harden and mature those software components based on usage and feedback from URC.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH

We present here RICH, a state-of-the-art two-dimensional hydrodynamic code based on Godunov’s method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time-advancement schemeS as well as a novel parallelization scheme based on Voronoi tessellation. Using our code, we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time-advancement scheme is more accurate and robust than is AREPO when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Lloyd iterations) and its effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. We perform an analytic analysis for finite difference schemes that reveals that a Lagrangian simulation is better than a Eulerian simulation in the case of a highly supersonic flow. Moreover, we show that Voronoi-based moving mesh schemes suffer from an error, which is resolution independent, due to inconsistencies between the flux calculation and the change in the area of a cell. Our code is publicly available as open source and designed in an object-oriented, user-friendly way that facilitates incorporation of new algorithms and physical processes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Status of PyGriffin Development for Integration into the NEAMS Workbench

The integration of Griffin into Workbench, by way of the PyGriffin code package, was initiated by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY-2022. PyGriffin was developed as a Python wrapper to Griffin to streamline the complex workflow involving mesh, cross section, and Griffin input generation, code execution, and results postprocessing, by improving user experience with high fidelity neutronics analysis through the Workbench GUI interface. PyGriffin can be used as a standalone application or through the PyARC code suite. FY-2023 saw the continued development of PyGriffin to expand its capabilities including advanced post-processing of simulation results and initial development of the Monte Carlo (MC) cross section generation workflow with the Shift MC code. PyGriffin was also approved for open-source software status in FY-2023, initiating plans to move the code repository into an open location. In addition to PyGriffin development, there were many improvements and developments made in PyARC leading to several software releases (latest version is v2.3.0), but the focus of this report is PyGriffin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Pseudoatom molecular dynamics plasma microfields

Spectral lines are powerful diagnostic tools for both laboratory and astrophysical plasmas, as their shape is sensitive to the plasma environment. The low-frequency component of the electric microfield is an important input for semi-analytic line broadening codes. Here, in this paper, we detail a new method of calculating plasma microfields using configuration-resolved pseudoatom molecular dynamics. This approach accounts for both quantum atomic structure and N-body effects, similar to density functional theory molecular dynamics, but with less computational cost. We present pseudoatom microfields at plasma conditions relevant for recent high energy density laboratory astrophysics experiments conducted at the Sandia Z-Machine, National Ignition Facility, and Linac Coherent Light Source. Compared to established microfield codes we find moderate deviations at solid density conditions and strong agreement at lower plasma densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

CARDINAL

Cardinal is a code that automates the workflow for iteratively solving neutral particle transport via the OpenMC code and computational fluid dynamics via the NekRS code within the Multiphysics Object Oriented Simulation Environment (MOOSE), a finite element framework for solving general partial differential equations. Cardinal establishes a mapping between the Monte Carlo geometry representation used in OpenMC and the spectral element mesh used in NekRS to other physics applications built upon the MOOSE framework, allowing high-resolution particle transport and fluid dynamics physics feedback to the MOOSE application ecosystem. Heat generation rates as solved by OpenMC are used as heat sources in MOOSE, while the temperature and density fields from NekRS are used as conjugate heat transfer boundary conditions and/or source fields in MOOSE. OpenMC, Nek5000, and MOOSE are open source, community-developed simulation codes.

NOVAK, APRILJEAN RODGERS↗
DOE OSTI2015

Beam emittance measurements and simulations of injector line for radio frequency quadrupole

A 400 keV deuteron (D{sup +}) radio frequency quadrupole (RFQ) accelerator has been designed, built, and commissioned at the Bhabha Atomic Research Centre, India. A beam line has been developed for injecting deuterons into the 400 keV RFQ. This comprises of a RF plasma ion source and a low energy beam transport (LEBT) system, consisting of two solenoid magnets and two steerer magnets. The ion source is characterized in terms of transverse beam emittance. A slit-wire scanner based emittance measurement setup has been developed for the transverse emittance measurements of H{sup +} and D{sup +} beams. The measured emittance values are found to be well within the acceptance value for the RFQ. These measured emittance parameters are used to optimize the solenoid fields in LEBT to match the beam from the ion source to RFQ. TRACEWIN simulation code is used for the beam transport simulations. The simulations show 99% transmission of D{sup +} beam through the RFQ, while 95% transmission has been measured experimentally.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Analytical formula for pellet fuel source density in toroidal plasma configurations based on an areal deposition model

An analytical formula for the fuel source density η F (Ψ) over a magnetic flux surface Ψ following pellet injection is derived. During pellet motion the ablated and ionized fuel will be initially deposited nearby the pellet, prior to subsequent homogenization within the magnetic flux surfaces. Here in our research, it is assumed that the deposition area projected normal to the magnetic field is on the order of the cross-sectional area of the ionized ablation column whose axis is aligned with the magnetic field. Previous pellet deposition codes [1] employed a point source deposition approximation, which overpredicts the surface-averaged density, particularly when the obliquity between the pellet trajectory and the flux surface normal approaches 90 degrees. A key finding of this study is that as the obliquity approaches 90 degrees, grazing trajectory, the areal model predicts a well behaved η F and removes the unphysical divergence stemming from the point model. To illustrate the difference exhibited between areal and point models for near grazing trajectories, we considered vertical pellet injection into a CFETR-like tokamak with a negative triangularity plasma cross section using a pellet ablation/deposition module.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A ray-tracing algorithm for ab-initio calculation of thermal load in undulator-based synchrotron beamlines

The OASYS suite and its powerful integration features are used to implement a ray-tracing algorithm to accurately calculate the thermal load in any component of an undulator-based synchrotron beamline. This is achieved by sampling and converting the SRW source of a given energy into a Shadow source and using the latter code to ray trace the full beamline. The accuracy of the algorithm is proved by reconstructing the full undulator radiation distribution through an aperture and comparing the result with direct calculaton of the total power using SRW. The algorithm is particularly suited to analyze cases with complex beamline layouts and optical elements, such as crystals, multilayers, and compound refractive lenses. Examples of its use to calculate the power load on elements of two of the feature beamlines at the Advanced Photon Source Upgrade Project and a comparison of the results with analytical calculations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Spatial distribution of the plasma parameters in the RF negative ion source prototype for fusion

A numerical model, based on the fluid plasma theory, has been used for description of the spatial distribution of the plasma parameters (electron density and temperature, plasma potential as well as densities of the three types of positive hydrogen ions) in the IPP prototype RF negative hydrogen ion source. The model covers the driver and the expansion plasma region of the source with their actual size and accounts for the presence of the magnetic filter field with its actual value and location as well as for the bias potential applied to the plasma grid. The obtained results show that without a magnetic filter the two 2D geometries considered, respectively, with an axial symmetry and a planar one, represent accurately the complex 3D structure of the source. The 2D model with a planar symmetry (where the E×B and diamagnetic drifts could be involved in the description) has been used for analysis of the influence, via the charged-particle and electron-energy fluxes, of the magnetic filter and of the bias potential on the spatial structure of the plasma parameters in the source. Benchmarking of results from the code to experimental data shows that the model reproduces the general trend in the axial behavior of the plasma parameters in the source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗