Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “C CODES”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Large language model evaluation for high–performance computing software development

We apply AI-assisted large language model (LLM) capabilities of GPT-3 targeting high-performance computing (HPC) kernels for (i) code generation, and (ii) auto-parallelization of serial code in C ++, Fortran, Python and Julia. Our scope includes the following fundamental numerical kernels: AXPY, GEMV, GEMM, SpMV, Jacobi Stencil, and CG, and language/programming models: (1) C++ (e.g., OpenMP [including offload], OpenACC, Kokkos, SyCL, CUDA, and HIP), (2) Fortran (e.g., OpenMP [including offload] and OpenACC), (3) Python (e.g., numpy, Numba, cuPy, and pyCUDA), and (4) Julia (e.g., Threads, CUDA.jl, AMDGPU.jl, and KernelAbstractions.jl). Kernel implementations are generated using GitHub Copilot capabilities powered by the GPT-based OpenAI Codex available in Visual Studio Code given simple + + prompt variants. To quantify and compare the generated results, we propose a proficiency metric around the initial 10 suggestions given for each prompt. For auto-parallelization, we use ChatGPT interactively giving simple prompts as in a dialogue with another human including simple “prompt engineering” follow ups. Results suggest that correct outputs for C++ correlate with the adoption and maturity of programming models. For example, OpenMP and CUDA score really high, whereas HIP is still lacking. We found that prompts from either a targeted language such as Fortran or the more general-purpose Python can benefit from adding language keywords, while Julia prompts perform acceptably well for its Threads and CUDA.jl programming models. Finally, we expect to provide an initial quantifiable point of reference for code generation in each programming model using a state-of-the-art LLM. Overall, understanding the convergence of LLMs, AI, and HPC is crucial due to its rapidly evolving nature and how it is redefining human-computer interactions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

A real-time expert system for self-repairing flight control

An integrated environment for specifying, prototyping, and implementing a self-repairing flight-control (SRFC) strategy is described. At an interactive workstation, the user can select paradigms such as rule-based expert systems, state-transition diagrams, and signal-flow graphs and hierarchically nest them, assign timing and priority attributes, establish blackboard-type communication, and specify concurrent execution on single or multiple processors. High-fidelity nonlinear simulations of aircraft and SRFC systems can be performed off-line, with the possibility of changing SRFC rules, inference strategies, and other heuristics to correct for control deficiencies. Finally, the off-line-generated SRFC can be transformed into highly optimized application-specific real-time C-language code. An application of this environment to the design of aircraft fault detection, isolation, and accommodation algorithms is presented in detail.

Gaither, S. A.↗
DOE OSTI2022

Uncertainty analysis for VERA problem 2 using the cell-code Condor v2.8.05

Condor is a cell-level neutronic calculation code that applies multi-group collision probabilities with heterogeneous response coupling method within generic geometry configurations. Under the Condor's code continuous development, the incorporation of up-to-date methodologies and state-of-the-art practices in reactor analysis represents a driving force. In this work, the capabilities of Condor v2.8.05 to develop an uncertainty analysis for realistic PWR-kind fuel assemblies are studied. The Total Monte Carlo approach is applied to quantify the impact of fabrication tolerances in the code's results for reactivity and power distributions, by means of randomly sampled input values using the VERA problem 2 as basis. The VERA problem 2 proposes a series of Westinghouse 2D 17 x 17-type fuel lattices, to be calculated reflected at beginning-of-life without Xe. The configurations correspond to a modern PWR. Selected neutronic parameters from Condor runs are thus analyzed in terms of the observed spread as well as the obtained distributions for the randomly perturbed cases, showing the capability of the code to handle the required input data, as well as its ability to provide valuable insights regarding uncertainty quantification.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Reactivity initiated accident uncertainty quantification for fuel assembly with subchannel code

UNIST CORE lab has developed a multiphysics coupling framework (MPCORE) consisting of Neutronics, Thermal Hydraulics and Fuel Performance modules. It can accommodate one-dimensional as well as sub-channel code Thermal Hydraulics (TH) module. Generally, running a transient requires more computational power due to the convergence of modules with each other. The difference between one-dimensional and sub-channel TH modules is studied in this research for a Reactivity Initiated Accident (RIA). Both TH modules are compared for a RIA uncertainty propagation in a single VERA fuel assembly with 2.11% enrichment. MPCORE is capable of analyzing the transient at any burnup point but for current work, only fresh fuel has been considered. The results have been obtained using dynamic gap heat conductance in FRAPTRAN. Peak centerline temperature, fuel enthalpy and DNBR are compared for both the approaches. The results indicate that the use of sub-channel code lead to greater safety margin for critical parameters. Computation time comparison is also presented for both the cases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Peregrine Software Development: Report on the Code Conversion From Python to C++

This work package seeks to convert the Peregrine software tool from its original Python implementation to a production version based on the C++ language. Peregrine is a powerful research platform with a multitude of advanced data analytics and data visualization functionalities. Developed by scientists to explore multimodal and multidimensional data related to the production of components using powder bed additive manufacturing processes, the tool implements state-of-the-art algorithms to assist machine users in making build or part quality determinations. Given that Peregrine is data-intensive, the goal of this conversion is to enhance the tool’s flexibility and interactivity and reduce the number of code dependencies to facilitate its deployment as part of the ongoing technology transfer campaign. This brief document provides an overview of Peregrine’s functionalities and capabilities, along with a detailed description of the core functionalities that have been implemented to date in the new C++ version. This document serves as a development update at the end of the first year of the ongoing conversion and will be regularly updated as progress continues.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2007

Driver Code for Adaptive Optics

A special-purpose computer code for a deformable-mirror adaptive-optics control system transmits pixel-registered control from (1) a personal computer running software that generates the control data to (2) a circuit board with 128 digital-to-analog converters (DACs) that generate voltages to drive the deformable-mirror actuators. This program reads control-voltage codes from a text file, then sends them, via the computer s parallel port, to a circuit board with four AD5535 (or equivalent) chips. Whereas a similar prior computer program was capable of transmitting data to only one chip at a time, this program can send data to four chips simultaneously. This program is in the form of C-language code that can be compiled and linked into an adaptive-optics software system. The program as supplied includes source code for integration into the adaptive-optics software, documentation, and a component that provides a demonstration of loading DAC codes from a text file. On a standard Windows desktop computer, the software can update 128 channels in 10 ms. On Real-Time Linux with a digital I/O card, the software can update 1024 channels (8 boards in parallel) every 8 ms.

Rao, Shanti↗
DOE OSTI2021

Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle

Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Characterizing Supernovae Neutrinos in MARLEY

This summer I worked on characterizing the uncertainties in neutrino-argon interactions due to the optical model potential. Functionally, this was a great deal of new experience coding in C++. I wrote a script using MARLEY(Gardiner 2019) to compute both the branching ratios and exclusive cross sections for a whole host of different decay processes from the unstable Potassium nucleus created by a v+Ar charged current interaction. These decays were computed through MARLEY, governed by a global parametrization of the Optical Model potential(Koning-Delaroche 2003). Helpfully, Pruitt+ 2023 provided public supplemental material outlining 400 separate "universes" for the optical model potential parameters. These "universes" represent uncertainties on the optical potential itself, and I wrote code to incorporate them into the MARLEY framework and quantify their impact on the predictions of neutrino-argon cross sections. This stands as the first computation of nuclear de-excitation systematic errors on cross sections in this energy band for v+Ar interactions, a critical step to understanding the science we can do with DUNE in the tens-of-MeV band.

Abu El-Haj, Luca↗
DOE OSTI2021

Review of Codes and Standards for Energy Storage Systems

This article identifies several examples of industry efforts and successes in removing gaps in energy storage (ES) Codes & Standards (C&S) by updating or creating and publishing new standards. A particular challenge discussed in this article is that while modern battery technologies including lithium ion (Li-ion) increase technical and economic viability of grid energy storage, newer battery technologies also present new or unknown risks to managing the safety of energy storage systems (ESS). There has been progress in filling gaps in published ES C&S that recognize and address the expanding range of technologies and their unique characteristics. However, there remain significant need and opportunity for researchers to contribute to the underlying knowledge base that informs development of technical references and standards, and ultimately the application of published standards for the effective and safe design and use of modern ESS.

Energy Storage, Codes and Standards↗
NASA NTRS2008

Guidance and Control Software Project Data - Volume 2: Development Documents

The Guidance and Control Software (GCS) project was the last in a series of software reliability studies conducted at Langley Research Center between 1977 and 1994. The technical results of the GCS project were recorded after the experiment was completed. Some of the support documentation produced as part of the experiment, however, is serving an unexpected role far beyond its original project context. Some of the software used as part of the GCS project was developed to conform to the RTCA/DO-178B software standard, "Software Considerations in Airborne Systems and Equipment Certification," used in the civil aviation industry. That standard requires extensive documentation throughout the software development life cycle, including plans, software requirements, design and source code, verification cases and results, and configuration management and quality control data. The project documentation that includes this information is open for public scrutiny without the legal or safety implications associated with comparable data from an avionics manufacturer. This public availability has afforded an opportunity to use the GCS project documents for DO-178B training. This report provides a brief overview of the GCS project, describes the 4-volume set of documents and the role they are playing in training, and includes the development documents from the GCS project. Volume 2 contains three appendices: A. Guidance and Control Software Development Specification; B. Design Description for the Pluto Implementation of the Guidance and Control Software; and C. Source Code for the Pluto Implementation of the Guidance and Control Software

Hayhurst, Kelly J.↗
DOE OSTI2023

Nonlinear magnetohydrodynamic modeling of current-drive-induced sawtooth-like crashes in the W7-X stellarator

Sawtooth-like core electron temperature crashes have been observed in W7-X experiments with electron cyclotron current drive. Here we present nonlinear single-fluid magnetohydrodynamic simulations of this phenomenon using the newly developed stellarator modeling capability of the M3D-C 1 code. The near-axis current drive gives rise to two ι = 1 resonances in the equilibrium rotational transform profile so that two consecutive (1, 1) internal kink modes are seen in the simulations. A small-amplitude crash at the inner resonance occurs first, which may correspond to the sawtooth precursors observed in the experiments. A bigger crash at the outer resonance then flattens the core temperature profile, which shows semi-quantitative agreements with experimental measurements on certain metrics such as the crash amplitude and the inversion radius of the temperature change. These results illustrate a likely mechanism of the current-drive-induced sawtooth-like crashes in W7-X and, to some extent, validate the stellarator modeling capability of M3D-C 1 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Coupled neutronic-thermal-hydraulic simulations of the European SFR core

Within the European SFR (Sodium-cooled Fast Reactor) - Safety Measures Assessment and Research Tools (ESFR-SMART) project, steady-state coupled simulation of the ESFR core has been performed using several core analysis packages, with the objective of quantifying the coupling effect. Focus is on the fuel Doppler effect and coolant expansion effect. Standalone neutronics calculations in TRACE/PARCS (PSI), DYN3D (HZDR) and WIMS (Jacobs) showed superb agreement with the reference Serpent power distribution (root mean square - rms discrepancies of 1.3%, 1.5% and 0.7% respectively). Results for COUNTHER (CIEMAT) were also in reasonable agreement, but with a somewhat higher discrepancy of 3.7%. Temperature distributions from thermal-hydraulic calculations were also compared and are found to be in good agreement. The effect of Doppler and coolant density feedback on core power distribution was predicted by TRACE/PARCS, DYN3D and WIMS to be between 0.4% and 0.8% rms difference in assembly powers. Reactivity coefficients for perturbations in the inlet temperature, flow rate and core power were shown to be negative for these three codes, with values of roughly -0.5 pcm/C. degree, -0.3 pcm/C. degree and -3.5 pcm/C. degree respectively. A preliminary investigation of differential thermal expansion effects indicates that this may have a significant effect on core power distribution of a few %, greater than anticipated a priori and may warrant inclusion in coupled core analysis to ensure the accurate calculation of power distributions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Coupling of nTRACER to COBRA-TF for full core high-fidelity analysis of VVERs

The evolution of computing clusters allows the use of high-resolution multi-physics solvers for reactor analysis. Despite the continuous expansion of VVER technology, only one other high-fidelity sub-pin multi-physics core solver is currently developed for safety analysis. To that end the Laboratory of Reactor Physics and Thermal-Hydraulics (LRT) of Paul Scherrer Institut is developing such a tool with the coupling of the neutronic code nTRACER and the sub-channel code COBRA-TF. This work follows the initial steps of the coupling and focuses to the extension of the core solver to full core VVER geometries. The X2 VVER-1000 benchmark is modeled with the novel core solver. The results are compared to the ones of a standalone nTRACER calculation where the feedbacks are provided by a simplified 1D thermal-hydraulic solver. Despite relatively good agreement in power distribution, the limitations of the simplified solver become apparent especially when comparing global temperature profiles. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Nonlinear magnetohydrodynamic modeling of ideal ballooning modes in high- β Wendelstein 7-X plasmas

We present nonlinear magnetohydrodynamic (MHD) simulations of high- β Wendelstein 7-X plasmas using the stellarator extension of the M3D- C 1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D- C 1 can play an instrumental role.

Zhou, Yao [Shanghai Jiao Tong University (China)] ↗
DOE OSTI2022

ANTARES: a new IRSN computational chain for safety analysis

Since some years, the IRSN has developed a new project to simulate accurately the interaction between neutronic and thermal hydraulic phenomena. This project called ANTARES (Advanced Neutronics and Thermal-hydraulic for the Analysis of Reactor Safety) couples the thermal hydraulic code CATHARE-3 with the 3D neutronic nodal code PARCS. The purpose of this project is to take the advantage of the complementary of these 2 codes, by the knowledge of the 3D core power distribution in CATHARE-3 (computed by PARCS) and improving the thermal hydraulic boundaries conditions in PARCS. In this way, the IRSN diversifies and enlarges its available tools to perform safety analysis with improved accuracy. The current technique usually adopted in France for the safety demonstrations for Reactivity Insertion Accident (RIA) is the so-called 'conservative' approach, which consists of reducing all the feedback (Doppler and moderator effects) and in modifying some physical quantities in such a way to increase a power peak. For this reason, these facilities ('penalties') have been implemented in ANTARES. In this paper we will give two examples of accidental transients that can be simulated with ANTARES: a Main Steam Line Break (MSLB) and an inadvertent boron dilution event. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development of a high-fidelity multi-physics coupling between MCNP6.2 and CTF4.0 for VVER applications

Ensuring system safety in the design, licensing, and operation phases is a priority in the nuclear industry. Performing extensive, full-scale reactor safety experiments is often prohibitive due to the large associated costs. Computational simulations offer an alternative safety analysis method, typically with significant cost reductions. Recent high-level developments in technology and increased availability of computational resources have allowed the development of high-fidelity, high-resolution multi-physics coupled codes. Such developments may be used to generate reference models for deterministic core calculations. Under this framework, the high-fidelity continuous energy Monte Carlo-based neutron transport code MCNP6.2 was coupled externally with the state-of-the-art thermal-hydraulics subchannel code, CTF4.0 for VVER (Water-Water Energetic Reactor) applications. A VVER-1000 fuel assembly model was used to demonstrate the capability of the coupled code. The converged coupled solution is compared to initial results, consisting of the first MCNP evaluation and first CTF evaluation after initialization. The VVER-1000 type assembly results are compared to other evaluations of the same assembly model. The findings indicate good agreement with expectations and reference cases, where available. The initial results of the coupled MCNP6.2/CTF4.0 calculations at the assembly level for steady-state calculations are presented in this study, which may support future work toward high-fidelity coupled full-core modeling capabilities. The multi-physics model may be further improved, enhanced, and expanded for both cycle depletion and transient applications. Such a multi-physics system will also be applicable to the VVER-1200 and other triangular lattice designs and support their deployment and operation safely and economically. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Quantification of the importance of the modeling fidelity in the Best Estimate Plus Uncertainty analysis

In this work, the multi-physics coupling of CTF-PARCS for both steady-state and transient simulations is evaluated focusing on the quantification of the importance of the modeling fidelity when applied to the Best Estimate Plus Uncertainty safety analysis. For the steady-state, a Picard iteration approach is adopted in the coupled code while for the transient an explicit approach has been developed as a first step. The CTF-PARCS coupling is applied to the REA at HFP BOC scenario proposed in the exercise I-2c of LWR-UAM benchmark Phase III. Two different modeling choices are investigated for the nominal calculations. The first consists of including the 3D burnup information in CTF, and the second of modeling the spacers grids in both PARCS and CTF. It is observed that due to the degradation of the fuel thermal conductivity, the 3D burnup inclusion in CTF leads to significant higher centerline temperatures than the base case, with a difference up to 300 K. The spacer grids modeling shows the expected power and temperature decrease in the axial positions of the spacer grids. For the uncertainty quantification study, the multi-physics modeling including the burnup information in CTF is used. The results show that the inserted reactivity has a relative standard deviation of 13.6% mainly attributed to the nuclear data. The maximum linear power has a relative standard deviation of 16.2% due to the nuclear data as well, while the maximum centerline temperature has a relative standard deviation of 3.5% primarily caused by the fuel thermal conductivity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Predicting deviations in ex-core detector response in Krsko NPP during an earthquake with Monte-Carlo neutron transport methods

To study the response of ex-core detectors during an earthquake, an approach using Monte Carlo code MCNP was developed. As a first step to predict the response of the ex-core detectors a detailed MCNP model of the Krsko reactor core was created. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krsko nuclear power plant. Due to the outside the core position of the neutron detectors, the hybrid code ADVANTG is used to generate weight windows to speed up neutron transport outside the reactor core. To be able to use ADVANTG, fixed neutron source had to be reconstructed from the criticality core calculation. It was found that the power range ex-core detector movement has a negligible effect on the value of the thermal neutron flux at the active part of the detector. However, the 5 cm radial movement of the intermediate-range detector leads to 7 % - 8 % change in thermal neutron flux within the active intermediate-range detector region. The analysis continued with the evaluation of the effect of core barrel movement on the ex-core detector response. It was determined that the 2 mm core barrel radial oscillation can lead to ∼ 4% change in thermal neutron flux within the active detector region. It was evaluated that reactor core movement inside the baffle can lead up to ∼ 1.9% and ∼ 1.3% deviation in power-range and intermediate-range detector signal respectively. Analysis showed that the mechanical movement of ex-core neutron detectors could not explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗