Search NASA⌕ Search

SEARCH · Search NASA

Results for “C 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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Wasatch: An architecture-proof multiphysics development environment using a Domain Specific Language and graph theory

To address the coding and software challenges of modern hybrid architectures, we propose an approach to multiphysics code development for high-performance computing. This approach is based on using a Domain Specific Language (DSL) in tandem with a directed acyclic graph (DAG) representation of the problem to be solved that allows runtime algorithm generation. When coupled with a large-scale parallel framework, the result is a portable development framework capable of executing on hybrid platforms and handling the challenges of multiphysics applications. In addition, we share our experience developing a code in such an environment – an effort that spans an interdisciplinary team of engineers and computer scientists.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Magnetic plasma expulsion

Magnetic confinement of plasmas is widely studied for restricting a plasma's extent to a given region. Presented is one scheme for what is referred to as magnetic plasma expulsion, which would prevent the plasma from entering a given region. Magnetic plasma expulsion techniques may allow instrumentation access to the interiors of plasmas. Magnetic plasma expulsion techniques may also make alternative approaches for magnetic plasma confinement possible. Particle-in-cell simulation with the Warp code is used to study magnetic plasma expulsion. Some conditions for achieving magnetic plasma expulsion are reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Theoretical quantification of shock-timing sensitivities for direct-drive inertial confinement fusion implosions on OMEGA

Using temporally shaped laser pulses, multiple shocks can be launched in direct-drive inertial confinement fusion implosion experiments to set the shell on a desired isentrope or adiabat. The velocity of the first shock and the times at which subsequent shocks catch up to it are measured through the VISAR diagnostic on OMEGA. Simulations reproduce these velocity and shock-merger time measurements when using laser pulses designed for setting mid-adiabat (α ~ 3) implosions, but agreement degrades for lower-adiabat (α ~ 1) designs. Several possibilities for this difference are studied: errors in placing the target at the center of irradiation (target offset), variations in energy between the different incident beams (power imbalance), and errors in modeling the laser energy coupled into the capsule. Simulation results indicate that shock timing is most sensitive to details of the density and temperature profiles in the coronal plasma, which influences the laser energy coupled into the target, and only marginally sensitive to target offset and beam power imbalance. A new technique under development to infer coronal profiles using x-ray self-emission imaging can be applied to the pulse shapes used in shock-timing experiments. In conclusion, this will help identify improved physics models to implement in codes and consequently enhance shock-timing predictive capability for low-adiabat pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Toroidal injection angle dependence of EC assisted plasma initiation at DIII-D

An experimental scan of the electron cyclotron waves (EC) toroidal injection angle in plasma breakdown is performed at the DIII-D tokamak. The second harmonic, extraordinary mode EC is used for the study. The dependence of n e and T e on the EC injection angle cannot be conclusively drawn from this study due to the large error bars in the n e and T e measurements. On the other hand, high T e data points are observed in some discharges which can be explained by nonlinear heating. The D α emission measurement shows a clear relation between the breakdown time and the injection angle. An experimental investigation of the cause of the dependence of breakdown delay on the EC injection angle suggests that when the injection angle is oblique, the EC heating after the reflection at the inboard wall may become ineffective and cause the breakdown delay even when the EC heating before and directly upon injection remains effective. A preliminary run of the heat and transport balance code DYON indicates that the obtained dataset is suitable for a quantitative validation of EC absorption models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Towards Resiliency Evaluation of Vector Programs

The systems resilience research community has developed methods to manually insert additional source-program level assertions to trap errors, and also devised tools to conduct fault injection studies for scalar program codes. In this work, we contribute the first vector oriented LLVM-level fault injector VULFI to help study the effects of faults in vector architectures that are of growing importance, especially for vectorizing loops. Using VULFI, we conduct a resiliency study of nine real-world vector benchmarks using Intel’s AVX and SSE extensions as the target vector instruction sets, and offer the first reported understanding of how faults affect vector instruction sets. We take this work further toward automating the insertion of resilience assertions during compilation. This is based on our observation that during intermediate (e.g., LLVM-level) code generation to handle full and partial vectorization, modern compilers exploit (and explicate in their code-documentation) critical invariants. These invariants are turned into error-checking code. We confirm the efficacy of these automatically inserted low-overhead error detectors for vectorized for-loops.

fault injection↗
DOE OSTI2016

VCTR: A Versatile Coupled Test Reactor Concept

The United States has currently no domestic fast neutron test and irradiation capabilities necessary to pursue a strong research and development program related to fast neutron reactors. This lack prevents the testing of fast reactor fuels in a prototypic environment, the testing of new safety features and of new coolants. A new irradiation test reactor is a costly investment that will operate for several decades and the initial specifications are crucial to ensure the relevance of such a facility in the long run. This new test reactor should be highly reconfigurable and allow operation as a fast test reactor or a thermal test reactor or a cou-pled fast-thermal test reactor. This report presents the on-going effort at Idaho National Laboratory (INL) to design a versatile coupled test reactor (VCTR—pronounced Victor) that could fulfill the R&D needs requiring high fast neutron fluxes but also high thermal neutron fluxes in configurations that could not be met in current thermal neutron test reactors such as ATR. The objective of this report is to summarize the efforts undertaken at INL from May 2015 to March 2016 aimed at understanding the main options for designing and building a versatile test reactor, with a primary mission to provide thermal and fast neutrons for irradiation of fuels and materials in a manner complementary to the ATR and HFIR reactors. Other missions can also be considered, including the validation of computer codes and the development of specific reactor technologies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Assessing the ASME Section III, Division 5, Class A Primary Load Design Rules Against Creep Notch Effects

This report assesses the ASME Section III, Division 5, Subsection HB, Subpart B rules covering the design and construction of high temperature Class A nuclear reactor components for their robustness against creep notch effects. The creep notch effect combines the effect of multiaxial stresses on material creep deformation and damage. This report considers both effects, including the potential for creep mechanism shifts affecting the extrapolation of design data from high stress, short-time experimental data to low stress, long time operating component conditions. The report summarizes the history of and literature on multiaxial creep and surveys the current ASME design rules dealing with multiaxial effects. The report then describes the results of two dedicated numerical studies, one assessing the robustness of the ASME rules against uncertainty in extrapolating from uniaxial creep test data to multiaxial component conditions and a second study examining the potential effects of a mechanism shift from high stress, dislocation-mediated creep to diffusion-dominated creep at lower stresses. The final conclusion of the report is that the ASME rules adequately guard against multiaxial creep failure, though there are several aspects of the Code that could be optimized to provide less over-conservative design predictions or to provide a more consistent design margin as a function of temperature and stress. In a few areas, particularly the potential for mechanism shifts outside the currently-available experimental data, the Code rules should be revaluated as additional experimental data and new modeling and simulation results become available

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Comparison of RESRAD-ONSITE and PRG/DCC Calculators and Application to Navy Environmental Media

The purpose of this project is to develop an improved Navy-wide approach for estimating risk and calculating cleanup goals for radiologically impacted environmental media (e.g., soils). This comparative analysis focused on the new features added (or improved) to the PRG/DCC Calculators since the 2015 comparison study (Yu et al. 2015a). RESRAD-ONSITE code Version 7.2 and the PRG/DCC Calculators available on the USEPA website between January and June of 2024 were used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

User’s Manual for RESRAD-RDD&IND Code Version 2: Vol. 2—User’s Guide for RESRAD-RDD&IND Code

Version 2.0 of the RESRAD-RDD&IND computer code is designed to support the implementation of protective action guides (PAGs) after a nuclear emergency incident including a radiological dispersal device (RDD) and/or an improvised nuclear device (IND) incident (EPA 2017). Eight different group types, addressing various decisions, are available for selection. The RESRAD-RDD&IND code calculates radiological doses, stay times, etc., for the selected group that the user wishes to focus on. (That is, the results for all the groups are not calculated simultaneously, and the input for those other groups do not matter, although some parameter values are shared between groups.) Version 2.0 has a user-friendly interface so that the RESRAD-RDD&IND code can be used with minimal training. For example, the user can select the major characteristics of the problem-event type, source term, and decision type from the left side of the interface and then calculate the results with the default assumptions for the exposure scenarios. More in-depth analysis would include specifying site-specific exposure scenario characteristics in the right side of the interface. The procedures for data entry and results viewing are self-explanatory. This is because common window maneuvering features and text instructions were incorporated in the interface design. General and context-specific help are available to aid users entering parameter values, as well. The RESRAD-RDD&IND computer code gives the user the option to select either an RDD or IND incident for analysis. For an RDD event analysis, 11 radionuclides (Am-241, Cf-252, Cm-244, Co-60, Cs-137, Ir-192, Po-210, Pu-238, Pu-239, Ra-226, and Sr-90) are included. These 11 radionuclides are the radionuclides most likely used for an RDD. More than 90 radionuclides can be selected for an IND event analysis. Initial default concentrations are provided for 44 radionuclides for a uranium-fueled IND event. These 44 radionuclides are those that would contribute significantly to the radiation dose associated with a uranium-fueled bomb detonation. The radionuclides generated from ingrowth of these 44 initial radionuclides are also automatically included in the analysis. Pu-239, Cs-134m, Ru-105, and Rb-89 and their progeny can be selected for analysis if they are detected and their concentrations are determined. This user’s guide, which is Volume 2 of the User’s Manual for RESRAD-RDD&IND Code Version 2, provides instructions to users on how to install the RESRAD-RDD&IND code, navigate the interface, and use the various features, including those discussed above, to set up an analysis and view/print the results in text outputs. Volume 1 of the User’s Manual for RESRAD-RDD&IND Code Version 2 (Yu et al. 2026), which contains descriptions of the methodology and theoretical basis for dose modeling and the mathematical equations implemented in the code, can be accessed and viewed through the Help menu in the code or can be downloaded from the RESRAD website (https://resrad.evs.anl.gov).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Dataset_for_Conserved_macromolecular_architecture_of_Poplar_secondary_cell_walls_revealed_by_ssNMR_and_atomistic_modeling

This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).

09 BIOMASS FUELS↗
DOE OSTI2024

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

WOMBAT: A Scalable and High-performance Astrophysical Magnetohydrodynamics Code

We present a new code for astrophysical magnetohydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged from a collaboration between Cray, Inc. and the University of Minnesota. This design utilizes MPI-RMA optimized for thread scaling, which allows the code to run extremely efficiently at very high thread counts ideal for the latest generation of multi-core and many-core architectures. Such performance characteristics are needed in the era of “exascale” computing. We describe and demonstrate our high-performance design in detail with the intent that it may be used as a model for other, future astrophysical codes intended for applications demanding exceptional performance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Hydrodynamic and Radiographic Toolbox (HART)

With a multi-lab and university team, we propose to develop new methods for a Hydrodynamic and Radiographic Toolbox (HART) that will enable a fuller and more extensive use of experimental radiographic data towards better characterizing and reducing uncertainties in predictive modeling of weapons performance. We will achieve this by leveraging recent developments in areas of computational imaging, statistical and machine learning, and reduced order modeling of hydrodynamics. In terms of software practices, by partnering with XCP (RISTRA Project) we will conform to recent XCP standards that are inline with modern software practices and standards and ensure compatibility and ease of inter-operability with existing codes. Three new activities under this proposal include (a) the development and use of deep learning-based surrogates to accelerate reconstruction and variational inference of density fields from radiographs of hydrotests, (b) a model-data fusion strategy that couples deep learning-based density reconstructions with fast hydrodynamics simulators to better constrain the reconstruction, and (c) a method for treating asymmetries using techniques adopted from limited view tomography. All three new activities will be based on improved treatment of scatter, noise, beam spot movement, detector blur, and flat fielding in the forward model, and a use of sophisticated priors to aid in the re construction. The improvements to the forward model and improved algorithmic design of the reconstruction when complete will be contained in the iterative reconstruction code SHIVA—a code project that we have recently initiated. The many ways in which machine learning can be used in the reconstruction work will be contained in a code HERMES that has been initiated with DTRA support. For example, the significant levels of acceleration that will likely be achieved by the use of machine learning techniques will permit us (and are required) to quantify uncertainties in density retrievals. Next, the two-way coupling between density reconstruction and model-based simulation of the hydrodynamics will be contained in code EREBUS, and will permit a fuller realization of the potential of the data to constrain the hydrodynamic model and better address issues related to asymmetries in the problem. Finally, we anticipate that the better consistency with physics achieved in our reconstructions will allow them to be used by X-Division more so than today.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Thermal Aging Effects on the Yield and Tensile Strength of 9Cr-1Mo-V (Grade 91)

The long-term exposure of 9Cr-1Mo-V (Grade 91) steel to elevated temperature can have a significant effect on reducing its yield and tensile strength. The yield and tensile strength changes, in turn, have potential implications to the assurance of the integrity of components constructed with this material over their design or intended lifetime. The ASME Boiler & Pressure Vessel Code (BPVC), Section III, Division 5 (III-5, high temperature reactors) provides tabulated reduction factors for Grade 91 yield and tensile strength as a function of exposure temperature and time up to 300,000 h. ASME BPVC III-5 s intent to extend these factors to an exposure duration of 500,000 h, the lack of available historic information to support the existing factors, and the recent development of a physics-based prediction model for ASME BPVC application are prime motivation for this study. This work describes results of the conventional time-temperature Hollomon–Jaffe parameter, strength reduction ratio prediction method using an updated, extensive Grade 91 unaged and related aged material strength database. The method, previously used by Oak Ridge National Laboratory in its evaluation of Grade 91 and likely used in development of the existing BPVC III-5 reduction factors, provides strength reduction ratio predictions useful for general component fitness-for-service assessments and for computing BPVC III-5 reduction factors as defined. Specific reduction factors to 500,000 h at 650 °C applicable to ASME BPVC III-5 are computed from the strength reduction ratios.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Exploration of the Super H-mode regime on DIII-D and potential advantages for burning plasma devices

A new high pedestal regime (“Super H-mode”) has been predicted and accessed on DIII-D. Super H-mode was first achieved on DIII-D using a quiescent H-mode edge, enabling a smooth trajectory through pedestal parameter space. By exploiting Super H-mode, it has been possible to access high pedestal pressures at high normalized densities. While elimination of Edge localized modes (ELMs) is beneficial for Super H-mode, it may not be a requirement, as recent experiments have maintained high pedestals with ELMs triggered by lithium granule injection. Simulations using TGLF for core transport and the EPED model for the pedestal find that ITER can benefit from the improved performance associated with Super H-mode, with increased values of fusion power and gain possible. Similar studies demonstrate that the Super H-mode pedestal can be advantageous for a steady-state power plant, by providing a path to increasing the bootstrap current while simultaneously reducing the demands on the core physics performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Results of the KRUSTY Warm Critical Experiments

The Kilowatt Reactor Using Stirling TechnologY (KRUSTY) was a prototypic nuclear-powered test of a 5-kW(thermal) Kilopower space reactor. This paper presents results from the KRUSTY warm critical experiments, which were completed prior to the final system test. The first set of criticals comprised cold or zero-power criticals; i.e., the core was not heated by fission power. These were followed by three warm criticals, where fission power heated the core to 200°C, 300°C, and 450°C, respectively. These criticals provided the data, confidence, and regulatory framework that were needed to proceed with the KRUSTY nuclear system test. The criticals also provided valuable data for the benchmarking of codes applicable to all nuclear systems. Finally, a comparison of KRUSTY results to pretest predictions is provided, and overall, the models matched the experimental results very closely.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Ab initio structure determination of nanocrystals of organic pharmaceutical compounds by electron diffraction at room temperature using a Timepix quantum area direct electron detector

A specialized quantum area detector for electron diffraction studies makes it possible to solve the structure of small organic compound nanocrystals in non-cryo conditions by direct methods. Until recently, structure determination by transmission electron microscopy of beam-sensitive three-dimensional nanocrystals required electron diffraction tomography data collection at liquid-nitrogen temperature, in order to reduce radiation damage. Here it is shown that the novel Timepix detector combines a high dynamic range with a very high signal-to-noise ratio and single-electron sensitivity, enabling ab initio phasing of beam-sensitive organic compounds. Low-dose electron diffraction data (∼0.013 e{sup −} Å{sup −2} s{sup −1}) were collected at room temperature with the rotation method. It was ascertained that the data were of sufficient quality for structure solution using direct methods using software developed for X-ray crystallography (XDS, SHELX) and for electron crystallography (ADT3D/PETS, SIR2014)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗