Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coding”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

photon-plasma: A modern high-order particle-in-cell code

We present the photon-plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach, the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Thresholds of surface codes on the general lattice structures suffering biased error and loss

A family of surface codes with general lattice structures is proposed. We can control the error tolerances against bit and phase errors asymmetrically by changing the underlying lattice geometries. The surface codes on various lattices are found to be efficient in the sense that their threshold values universally approach the quantum Gilbert-Varshamov bound. We find that the error tolerance of the surface codes depends on the connectivity of the underlying lattices; the error chains on a lattice of lower connectivity are easier to correct. On the other hand, the loss tolerance of the surface codes exhibits an opposite behavior; the logical information on a lattice of higher connectivity has more robustness against qubit loss. As a result, we come upon a fundamental trade-off between error and loss tolerances in the family of surface codes with different lattice geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Tokamak plasma high field side response to an n = 3 magnetic perturbation: a comparison of 3D equilibrium solutions from seven different codes

In comparing equilibrium solutions for a DIII-D shot that is amenable to analysis by both stellarator and tokamak three-dimensional (3D) equilibrium codes, a significant disagreement has been seen between solutions of the VMEC stellarator equilibrium code and solutions of tokamak perturbative 3D equilibrium codes. In this research, the source of that disagreement has been investigated, and that investigation has led to new insights into the domain of validity of the different equilibrium calculations, and to a finding that the manner in which localized screening currents at low order rational surfaces are handled can affect global properties of the equilibrium solution. The perturbative treatment has been found to break down at surprisingly small perturbation amplitudes due to overlap of the calculated perturbed flux surfaces, and that treatment is not valid in the pedestal region of the DIII-D shot studied. The perturbative treatment is valid, however, further into the interior of the plasma, and flux surface overlap does not account for the disagreement investigated here. Calculated equilibrium solutions for simple model cases and comparison of the 3D equilibrium solutions with those of other codes indicate that the disagreement arises from a difference in handling of localized currents at low order rational surfaces, with such currents being absent in VMEC and present in the perturbative codes. The significant differences in the global equilibrium solutions associated with the presence or absence of very localized screening currents at rational surfaces suggests that it may be possible to extract information about localized currents from appropriate measurements of global equilibrium plasma properties. Therefore, that would require improved diagnostic capability on the high field side of the tokamak plasma, a region difficult to access with diagnostics.

3-DIMENSIONAL MAGNETOHYDRODYNAMIC EQUILIBRIA↗
DOE OSTI2013

YNOGK: A NEW PUBLIC CODE FOR CALCULATING NULL GEODESICS IN THE KERR SPACETIME

Following the work of Dexter and Agol, we present a new public code for the fast calculation of null geodesics in the Kerr spacetime. Using Weierstrass's and Jacobi's elliptic functions, we express all coordinates and affine parameters as analytical and numerical functions of a parameter p, which is an integral value along the geodesic. This is the main difference between our code and previous similar ones. The advantage of this treatment is that the information about the turning points does not need to be specified in advance by the user, and many applications such as imaging, the calculation of line profiles, and the observer-emitter problem, become root-finding problems. All elliptic integrations are computed by Carlson's elliptic integral method as in Dexter and Agol, which guarantees the fast computational speed of our code. The formulae to compute the constants of motion given by Cunningham and Bardeen have been extended, which allow one to readily handle situations in which the emitter or the observer has an arbitrary distance from, and motion state with respect to, the central compact object. The validation of the code has been extensively tested through applications to toy problems from the literature. The source FORTRAN code is freely available for download on our Web site http://www1.ynao.ac.cn/{approx}yangxl/yxl.html.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

RAY-RAMSES: a code for ray tracing on the fly in N-body simulations

We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zel'dovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Logical error rates for the surface code under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions

With fault-tolerant quantum computing (FTQC) on the horizon, it is critical to understand sources of logical errors in plausible hardware implementations of quantum error-correcting codes. Detailed error modeling of computational instructions on particular FTQC architectures will enable the better prediction of error propagation in FT-encoded quantum circuits while revealing where greater attention is needed in hardware design. In this work, we consider logical error rates for the surface code implemented on a hypothetical grid-based trapped-ion quantum charge-coupled device architecture. Specifically, we construct logical channels for the idling surface code and examine its diamond error under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions. We include the coherent dephasing noise that is known to accumulate during physical qubit idling and transport in these systems, determining idling and transport durations using the time-resolved output of an open-source trapped-ion surface code compiler. To estimate expectation values of logical Pauli observables following hardware circuits containing non-Clifford sources of noise, we utilize a Monte Carlo technique to sample from an underlying quasiprobability distribution of Clifford circuits that we independently simulate in a phase-sensitive fashion. We verify error suppression up to code distance 𝑑 = 11 at coherent dephasing rates near and below those of current-generation trapped-ion quantum computers and find that logical error rates align with those of analogous fully stochastic simulations in this regime. Exploring higher dephasing rates at 𝑑 = 3−5, we find evidence for growing coherent rotations about all three logical Pauli axes, increased diagonal logical error process matrix elements relative to those of stochastic simulations, and a reduced dephasing rate threshold. Overall, our work paves a way toward realistic hardware emulation of small fault-tolerant quantum processes, e.g., members of an FTQC instruction set.

Quantum benchmarking↗
DOE OSTI2016

Towards performance-portability of the Albany Finite Element analysis code using the Kokkos library of Trilinos

Performance portability on heterogeneous high-performance computing (HPC) systems is a major challenge faced today by code developers: parallel code needs to execute correctly as well as with high performance on machines with different architectures, operating systems, and software libraries. The Finite Element Method (FEM) is a popular and flexible method for discretizing partial differential equations arising in a wide variety of scientific, engineering, and industry applications that require HPC. This paper presents some preliminary results pertaining to our development of a performance portable implementation of the FEM-based Albany code. Performance portability is achieved using the Kokkos library of Trilinos. We present performance results for two different physics simulations modules in Albany: the Aeras global atmosphere dynamical code and the FELIX land-ice solver. As a result, numerical experiments show that our single code implementation gives reasonable performance across two multi-core/many-core architectures: NVIDIA GPUs and multi-core CPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Reactor Fuel Isotopics and Code Validation for Nuclear Applications

Experimentally measured isotopic concentrations of well characterized spent nuclear fuel (SNF) samples have been collected and analyzed by previous researchers. These sets of experimental data have been used extensively to validate the accuracy of depletion code predictions for given sets of burnups, initial enrichments, and varying power histories for different reactor types. The purpose of this report is to present the diversity of data in a concise manner and summarize the current accuracy of depletion modeling. All calculations performed for this report were done using the Oak Ridge Isotope GENeration (ORIGEN) code, an internationally used irradiation and decay code solver within the SCALE comprehensive modeling and simulation code. The diversity of data given in this report includes key actinides, stable fission products, and radioactive fission products. In general, when using the current ENDF/B-VII.0 nuclear data libraries in SCALE, the major actinides are predicted to within 5% of the measured values. Large improvements were seen for several of the curium isotopes when using improved cross section data found in evaluated nuclear data file ENDF/B-VII.0 as compared to ENDF/B-V-based results. The impact of the flux spectrum on the plutonium isotope concentrations as a function of burnup was also shown. The general accuracy noted for the actinide samples for reactor types with burnups greater than 5,000 MWd/MTU was not observed for the low-burnup Hanford B samples. More work is needed in understanding these large discrepancies. The stable neodymium and samarium isotopes were predicted to within a few percent of the measured values. Large improvements were seen in prediction for a few of the samarium isotopes when using the ENDF/B-VII.0 libraries compared to results obtained with ENDF/B-V libraries. Very accurate predictions were obtained for 133 Cs and 153 Eu. However, the predicted values for the stable ruthenium and rhodium isotopes varied as much as 10% of the measured values, and 109 Ag was consistently over-predicted by as much as 170%. In general, there is a larger uncertainty for modeling radioactive fission products when compared to either the actinides or the stable fission products in SNF. The relative C/E ratios ranged from a few percent for 137 Cs up to 60% and 100% for 106 Ru and 125 Sb, respectively. Limited or no radioactive fission products data exist in the current data sets for reactor types other than PWRs and BWRs. More work is needed in obtaining a greater diversity of radioactive fission product data. While performing this survey, issues leading to inconsistencies in nuclear fission yield data were discovered that specifically impacted the fission product noble gases. Emphasis was given to this legacy data, and corrective actions were taken as described in this report. After the fission yield data were corrected, the stable xenon and krypton fission products were predicted to within 5% of their measurements. However, preliminary results not explicitly given in this report indicate that the relative C/E ratio for the radioactive isotope 85 Kr varied as much as 10%. Due to the complex migration and the difficulty in measuring noble gases in the fuel, a more thorough investigation is needed to understand how accurately depletion codes can calculate these gas concentrations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Verification and Validation of the PLTEMP/ANL Code for Thermal-Hydraulic Analysis of Experimental and Test Reactors

The document compiles in a single volume several verification and validation works done for the PLTEMP/ANL code during the years of its development and improvement. Some works that are available in the open literature are simply referenced at the outset, and are not included in the document. PLTEMP has been used in conversion safety analysis reports of several US and foreign research reactors that have been licensed and converted. A list of such reactors is given. Each chapter of the document deals with the verification or validation of a specific model. The model verification is usually done by comparing the code with hand calculation, Microsoft spreadsheet calculation, or Mathematica calculation. The model validation is done by comparing the code with experimental data or a more validated code like the RELAP5 code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

High-fidelity plasma codes for burn physics

Accurate predictions of equation of state (EOS), ionic and electronic transport properties are of critical importance for high-energy-density plasma science. Transport coefficients inform radiation-hydrodynamic codes and impact diagnostic interpretation, which in turn impacts our understanding of the development of instabilities, the overall energy balance of burning plasmas, and the efficacy of self-heating from charged-particle stopping. Important processes include thermal and electrical conduction, electron-ion coupling, inter-diffusion, ion viscosity, and charged particle stopping. However, uncertainties in these coefficients are not well established. Fundamental plasma science codes, also called high-fidelity plasma codes, are a relatively recent computational tool that augments both experimental data and theoretical foundations of transport coefficients. This paper addresses the current status of HFPC codes and their future development, and the potential impact they play in improving the predictive capability of the multi-physics hydrodynamic codes used in HED design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Code Verification Capabilities and Assessments in Support of ASC V&V Level 2 Milestone #6035

This document provides a summary of the code verification activities supporting the FY17 Level 2 V&V milestone entitled “Deliver a Capability for V&V Assessments of Code Implementations of Physics Models and Numerical Algorithms in Support of Future Predictive Capability Framework Pegposts.” The physics validation activities supporting this milestone are documented separately. The objectives of this portion of the milestone are: 1) Develop software tools to support code verification analysis; 2) Document standard definitions of code verification test problems; and 3) Perform code verification assessments (focusing on error behavior of algorithms). This report and a set of additional standalone documents serve as the compilation of results demonstrating accomplishment of these objectives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The Alabama Energy Code Compliance Project

The U.S. Department of Energy (DOE) Building Energy Codes Program launched a series of research studies in 2014 across eight U.S. states (shown in Figure 1) to investigate typical energy efficiency practices in residential buildings and identify opportunities to reduce household utility bills. A primary goal of the study was to help document baseline practices for new single-family construction, target areas for improvement relative to the state energy code, and quantify the related savings potential. This information is intended to assist states in measuring compliance with their energy codes and to identify areas of focus for future education and training initiatives. Following the original baseline studies, project teams in each state conducted a series of education, training, and outreach activities aimed at improving compliance with the state residential energy code. A second baseline study was then implemented in each state with the objective of determining whether the investment in education, training, and outreach programs had produced a significant and measurable change in statewide residential energy use. Through the program, DOE planned to establish a sufficient data set to represent baseline construction trends across states and detect significant changes in energy use resulting from education and training activities. This report outlines the application of the study in the state of Alabama, which was led by the Institute for Market Transformation (IMT). Data gathered have been submitted to Pacific Northwest National Laboratory (PNNL) for analysis. Findings available to date are published as PNNL technical reports on the DOE Building Energy Codes Program website.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Building Energy Codes Program: National Benefits Assessment (1992-2040)

Building energy codes and standards set minimum requirements for energy-efficient design and construction for new and renovated buildings, assuring reductions in energy use and greenhouse gas emissions over the life of buildings. DOE, through the Building Energy Codes Program, supports the improvement of energy efficiency in buildings. BECP periodically assesses the impacts of its activities by estimating historical and projected energy savings, consumer savings, and avoided emissions. The Pacific Northwest National Laboratory conducted the codes benefits assessment in support of the BECP. Underlying the assessment is a series of calculations that estimate and compare energy savings under two scenarios: "with BECP" and "without BECP." The analysis covers the years 1992-2040 and includes comparing the nominal energy savings (assuming 100 percent adoption and compliance) attributable to different code versions, determining the applicable floor space (both residential and commercial) subject to the code, and estimating the final energy savings by adjusting nominal energy savings for the applicable floor space according to the estimated actual adoption and compliance levels. The resulting estimates of energy consumption for each scenario are compared, and the difference equals the impact of BECP activities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Impacts of Model Building Energy Codes – Interim Update

The U.S. Department of Energy (DOE) Building Energy Codes Program (BECP) periodically evaluates state and national impacts associated with energy codes in residential and commercial buildings. Pacific Northwest National Laboratory (PNNL), funded by DOE, assessed the prospective impacts of national model building energy codes from 2010 through 2040. The model codes of interest in this report are the International Energy Conservation Code (IECC) for residential buildings and ASHRAE Standard 90.1 for commercial buildings. This report describes the methodology used for the assessment and presents the impacts in terms of energy savings, consumer cost savings, and reduced CO 2 emissions at the state and national levels.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Impacts of Model Building Energy Codes

The Department of Energy (DOE) Building Energy Codes Program (BECP) periodically evaluates state and national impacts associated with energy codes in residential and commercial buildings. Pacific Northwest National Laboratory (PNNL), funded by DOE, assessed the prospective impacts of national model building energy codes from 2010 through 2040. The model codes of interest in this report are the International Energy Conservation Code (IECC) for residential buildings and ASHRAE Standard 90.1 for commercial buildings. This report describes the methodology used for the assessment and presents the impacts in terms of energy savings, consumer cost savings, and reduced CO 2 emissions at the state and national levels.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Code of Record: DOE Standard Canister (DOESC)

The United States Department of Energy (DOE) Spent Nuclear Fuel (SNF) Packaging Demonstration seeks to develop and demonstrate the designs, technology, processes, and regulatory framework for packaging DOE SNF for road-ready dry storage (RRDS); and establish the processes that will be used in a future production facility. The Packaging Demonstration will utilize the DOE Standard Canister (DOESC) for packaging select DOE-managed SNF types for interim storage, transportation and disposal as part of RRDS. Placing DOE-managed SNF into RRDS is part of the strategic framework for SNF on the Idaho National Laboratory (INL) site. To comply with DOE, INL and Nuclear Regulatory Commission (NRC) requirements, this Code of Record establishes the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC), Section III, Division 3 as the Code that will govern DOESC and internal support structure constructioni, with certain clarifications. This Code of Record establishes a rationale for proceeding without construction certification of the DOESC (i.e., “N-stamping) and the extent to which a Registered Professional Engineer is required for the DOE Spent Fuel Packaging Demonstration. Given the (i) standard industry practice to pursue independent licensure of commercial storage casks and transportation packages by the NRC in lieu of ASME certification and (ii) guidance provided by 10 CFR 830, DOESC construction activities need not be certified (i.e., “stamped”) to the ASME BPVC. However, additional quality assurance requirements will apply as outlined in this Code of Record.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

TEA: A CODE CALCULATING THERMOCHEMICAL EQUILIBRIUM ABUNDANCES

We present an open-source Thermochemical Equilibrium Abundances (TEA) code that calculates the abundances of gaseous molecular species. The code is based on the methodology of White et al. and Eriksson. It applies Gibbs free-energy minimization using an iterative, Lagrangian optimization scheme. Given elemental abundances, TEA calculates molecular abundances for a particular temperature and pressure or a list of temperature–pressure pairs. We tested the code against the method of Burrows and Sharp, the free thermochemical equilibrium code Chemical Equilibrium with Applications (CEA), and the example given by Burrows and Sharp. Using their thermodynamic data, TEA reproduces their final abundances, but with higher precision. We also applied the TEA abundance calculations to models of several hot-Jupiter exoplanets, producing expected results. TEA is written in Python in a modular format. There is a start guide, a user manual, and a code document in addition to this theory paper. TEA is available under a reproducible-research, open-source license via https://github.com/dzesmin/TEA.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Long-range correlation properties of coding and noncoding DNA sequences: GenBank analysis

An open question in computational molecular biology is whether long-range correlations are present in both coding and noncoding DNA or only in the latter. To answer this question, we consider all 33301 coding and all 29453 noncoding eukaryotic sequences--each of length larger than 512 base pairs (bp)--in the present release of the GenBank to dtermine whether there is any statistically significant distinction in their long-range correlation properties. Standard fast Fourier transform (FFT) analysis indicates that coding sequences have practically no correlations in the range from 10 bp to 100 bp (spectral exponent beta=0.00 +/- 0.04, where the uncertainty is two standard deviations). In contrast, for noncoding sequences, the average value of the spectral exponent beta is positive (0.16 +/- 0.05) which unambiguously shows the presence of long-range correlations. We also separately analyze the 874 coding and the 1157 noncoding sequences that have more than 4096 bp and find a larger region of power-law behavior. We calculate the probability that these two data sets (coding and noncoding) were drawn from the same distribution and we find that it is less than 10(-10). We obtain independent confirmation of these findings using the method of detrended fluctuation analysis (DFA), which is designed to treat sequences with statistical heterogeneity, such as DNA's known mosaic structure ("patchiness") arising from the nonstationarity of nucleotide concentration. The near-perfect agreement between the two independent analysis methods, FFT and DFA, increases the confidence in the reliability of our conclusion.

Non-NASA Center↗