Search NASA⌕ Search

SEARCH · Search NASA

Results for “code modernization”

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

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗
DOE OSTI2017

Dependency graph for code analysis on emerging architectures

Direct acyclic dependency (DAG) graph is becoming the standard for modern multi-physics codes.The ideal DAG is the true block-scheme of a multi-physics code. Therefore, it is the convenient object for insitu analysis of the cost of computations and algorithmic bottlenecks related to statistical frequent data motion and dymanical machine state.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Improved Fracture Models for Relocation Modeling

The BISON fuel performance code is being developed to provide a modern tool that has the flexibility to analyze a wide variety of fuel forms and to model conditions and phenomena that could not be represented in legacy tools. There are a number of motivations for this, including providing support for development of advanced fuel with improved accident tolerance for existing light water (LWR) reactors, improved understanding of mechanisms in fuel designs in current use in a wider variety of conditions, and facilitating the development of fuel for advanced reactor designs. To accomplish these goals, it is clear that BISON must rely on models of fuel behavior that are based on fundamental physical behavior, rather than on empirical correlations that represent that behavior in a simplified fashion. BISON still does employ many empirical models that were originally developed for other fuel performance codes, but efforts are underway to replace these with models that are more physically based. Radial relocation in LWR fuel is an example of a phenomenon that is currently represented by an empirical model, but which is ripe for replacement by physically based models. During normal operation, ceramic LWR fuel experiences significant fracturing that is driven by spatially nonuniform volumetric expansion. This occurs due to the significant thermal gradients within the fuel that occur in fresh and irradiated fuel, as well as nonuniform swelling due to fission products that occurs over longer-term irradiation. Fracture and fragmentation of fuel allows the outer radius of the fuel pellet to expand due to the loss of mechanical constraint and outward radial migration of fragments. This radial expansion has a significant effect on the fuel system response because it decreases the size of the gap between the fuel and cladding. This decreases the thermal resistance across that gap, which leads to decreased fuel centerline temperatures. This report documents recent work to improve the ability of BISON to model radial relocation using the extended finite element method (XFEM). These enhancements include the ability to include cohesive zone models, improved code architecture for propagating cracks based on fracture integrals, and improved handling of fuel and cladding interfaces with XFEM. These capabilities are demonstrated on a simulation of radial relocation including a surface interaction model that enforces residual opening.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1994

Electromagnetic on-aircraft antenna radiation in the presence of composite plates

The UTD-based NEWAIR3 code is modified such that it can model modern aircraft by composite plates. One good model of conductor-backed composites is the impedance boundary condition where the composites are replaced by surfaces with complex impedances. This impedance-plate model is then used to model the composite plates in the NEWAIR3 code. In most applications, the aircraft distorts the desired radiation pattern of the antenna. However, test examples conducted in this report have shown that the undesired scattered fields are minimized if the right impedance values are chosen for the surface impedance plates.

Kan, S. H-T.↗
NASA NTRS2007

General Mission Analysis Tool (GMAT) Architectural Specification. Draft

Early in 2002, Goddard Space Flight Center (GSFC) began to identify requirements for the flight dynamics software needed to fly upcoming missions that use formations of spacecraft to collect data. These requirements ranged from low level modeling features to large scale interoperability requirements. In 2003 we began work on a system designed to meet these requirement; this system is GMAT. The General Mission Analysis Tool (GMAT) is a general purpose flight dynamics modeling tool built on open source principles. The GMAT code is written in C++, and uses modern C++ constructs extensively. GMAT can be run through either a fully functional Graphical User Interface (GUI) or as a command line program with minimal user feedback. The system is built and runs on Microsoft Windows, Linux, and Macintosh OS X platforms. The GMAT GUI is written using wxWidgets, a cross platform library of components that streamlines the development and extension of the user interface Flight dynamics modeling is performed in GMAT by building components that represent the players in the analysis problem that is being modeled. These components interact through the sequential execution of instructions, embodied in the GMAT Mission Sequence. A typical Mission Sequence will model the trajectories of a set of spacecraft evolving over time, calculating relevant parameters during this propagation, and maneuvering individual spacecraft to maintain a set of mission constraints as established by the mission analyst. All of the elements used in GMAT for mission analysis can be viewed in the GMAT GUI or through a custom scripting language. Analysis problems modeled in GMAT are saved as script files, and these files can be read into GMAT. When a script is read into the GMAT GUI, the corresponding user interface elements are constructed in the GMAT GUI. The GMAT system was developed from the ground up to run in a platform agnostic environment. The source code compiles on numerous different platforms, and is regularly exercised running on Windows, Linux and Macintosh computers by the development and analysis teams working on the project. The system can be run using either a graphical user interface, written using the open source wxWidgets framework, or from a text console. The GMAT source code was written using open source tools. GSFC has released the code using the NASA open source license.

Hughes, Steven P.↗
DOE OSTI2020

High-performance computing in water resources hydrodynamics

In this work, we present a vision of future water resources hydrodynamics codes that can fully utilize the strengths of modern high-performance computing. The advances to computing power, formerly driven by the improvement of central processing unit processors, now focus on parallel computing and, in particular, the use of graphics processing units (GPUs). However, this shift to a parallel framework requires refactoring the code to make efficient use of the data as well as changing even the nature of the algorithm that solves the system of equations. These concepts along with other features such as the precision for the computations, dry regions management, and input/output data are analyzed in this paper. A 2D multi-GPU flood code applied to a large-scale test case is used to corroborate our statements and ascertain the new challenges for the next-generation parallel water resources codes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

InSAR Scientific Computing Environment

This computing environment is the next generation of geodetic image processing technology for repeat-pass Interferometric Synthetic Aperture (InSAR) sensors, identified by the community as a needed capability to provide flexibility and extensibility in reducing measurements from radar satellites and aircraft to new geophysical products. This software allows users of interferometric radar data the flexibility to process from Level 0 to Level 4 products using a variety of algorithms and for a range of available sensors. There are many radar satellites in orbit today delivering to the science community data of unprecedented quantity and quality, making possible large-scale studies in climate research, natural hazards, and the Earth's ecosystem. The proposed DESDynI mission, now under consideration by NASA for launch later in this decade, would provide time series and multiimage measurements that permit 4D models of Earth surface processes so that, for example, climate-induced changes over time would become apparent and quantifiable. This advanced data processing technology, applied to a global data set such as from the proposed DESDynI mission, enables a new class of analyses at time and spatial scales unavailable using current approaches. This software implements an accurate, extensible, and modular processing system designed to realize the full potential of InSAR data from future missions such as the proposed DESDynI, existing radar satellite data, as well as data from the NASA UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar), and other airborne platforms. The processing approach has been re-thought in order to enable multi-scene analysis by adding new algorithms and data interfaces, to permit user-reconfigurable operation and extensibility, and to capitalize on codes already developed by NASA and the science community. The framework incorporates modern programming methods based on recent research, including object-oriented scripts controlling legacy and new codes, abstraction and generalization of the data model for efficient manipulation of objects among modules, and well-designed module interfaces suitable for command- line execution or GUI-programming. The framework is designed to allow users contributions to promote maximum utility and sophistication of the code, creating an open-source community that could extend the framework into the indefinite future.

Rosen, Paul A.↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Agentic AI vs ML-Based Autotuning: A Comparative Study for Loop Reordering Optimization

High Performance Computing (HPC) applications rely heavily on code optimizations to achieve good performance on modern CPU and GPU architectures. Traditional Machine Learning auto-tuning approaches have demonstrated success in exploring high-dimensional spaces, but they often require expensive compile-run evaluations and lack adaptability for large HPC applications. The recent advances in Large Language Models (LLMs) and Agentic AI systems raise intriguing questions about the potential of these approaches to address specific optimization methodologies. This work aims to answer an essential question for the HPC community: “How Agentic AI Systems Compare to Traditional ML Autotuning Techniques?” To address this question, we present a comparative analysis between a traditional ML-based optimization approach and an Agentic AI system, evaluating their respective capabilities and limitations for loop-level optimization. In addition, we introduced a new Agentic AI system named LoopGen-AI using three different Large Language Models: GPT-4.1, Claude 4.0, and Gemini 2.5. A key finding is that LoopGen-AI achieves competitive per-formance with only a few program runs, the reasoning logs from the agents revealed that their decisions rely heavily on the combination of semantic understanding of the target kernel with dynamic feedback from the environment, highlighting a promising new dimension in performance tuning. In contrast, ML-based autotuners focus on statistical exploration, and require orders of magnitude more runs to reach peak performance. Additionally, our analysis shows that prompt engineering, particularly using Persona + Context Manager patterns, significantly impacts the effectiveness of Agentic AI. Our results indicate that while Agentic AI systems are not yet a complete replacement for ML-based autotuners, it can effectively complement traditional methods.

Rosas, Miguel Romero↗
DOE OSTI2021

Nomenclatural issues concerning cultured yeasts and other fungi: why it is important to avoid unneeded name changes

The unambiguous application of fungal names is important to communicate scientific findings. Names are critical for (clinical) diagnostics, legal compliance, and regulatory controls, such as biosafety, food security, quarantine regulations, and industrial applications. Consequently, the stability of the taxonomic system and the traceability of nomenclatural changes is crucial for a broad range of users and taxonomists. The unambiguous application of names is assured by the preservation of nomenclatural history and the physical organisms representing a name. Fungi are extremely diverse in terms of ecology, lifestyle, and methods of study. Predominantly unicellular fungi known as yeasts are usually investigated as living cultures. Methods to characterize yeasts include physiological (growth) tests and experiments to induce a sexual morph; both methods require viable cultures. Thus, the preservation and availability of viable reference cultures are important, and cultures representing reference material are cited in species descriptions. Historical surveys revealed drawbacks and inconsistencies between past practices and modern requirements as stated in the International Code of Nomenclature for Algae, Fungi, and Plants (ICNafp). Improper typification of yeasts is a common problem, resulting in a large number invalid yeast species names. With this opinion letter, we address the problem that culturable microorganisms, notably some fungi and algae, require specific provisions under the ICNafp. We use yeasts as a prominent example of fungi known from cultures. But viable type material is important not only for yeasts, but also for other cultivable Fungi that are characterized by particular morphological structures (a specific type of spores), growth properties, and secondary metabolites. We summarize potential proposals which, in our opinion, will improve the stability of fungal names, in particular by protecting those names for which the reference material can be traced back to the original isolate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

SW4 final report for iCOE

The seismic simulation code SW4 was adapted to run on modern supercomputers. SW4 solves the elastic wave equation in second order formulation. It uses a fourth order accurate finite difference method with explicit time stepping. The approximation conserves energy perfectly, and is stable under a CFL restriction on the time step.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Laying the Foundation for a Predictive Theory of Electroweak Nuclear Phenomena

Reliable model-independent predictions of nuclear processes such as beta-decay and neutrinonucleus scattering are important for many fundamental science and mission-related investigations; For example, accurate descriptions of as beta-decay are required for modeling the anti-neutrino spectrum of nuclear reactor, whether for trying to detect new physics such as sterile neutrinos, or remotely identifying the purpose of that reactor. For these activities, we need accurate calculations, as well as careful and thorough analyses of theoretical error and uncertainty in the underlying models and the many-body methods and codes used in solving that model. Unfortunately, modern theoretical calculations of these processes are inadequate for this task. At the heart of this inadequacy is the interaction between nucleons and other particles known as nuclear currents. Historically, descriptions of nuclear currents took into account only the simplest of particle interactions. Recently, the development of an effective field theory known as chiral EFT has been used to derive large corrections to these particle interactions. However, the complex nature of these corrections inhibited their use in ab initio methods of predicting nuclear processes. We developed a novel method to address these complexities and applied it in a series of benchmark and demonstrative calculations using state-of-the-art quantum few-nucleon solvers. Specifically, we studied the role of two-nucleon currents in electroweak processes using light nuclei (i.e., nuclei containing fewer than five nucleons) as a workbench to compute these processes and quantify their uncertainties using Bayesian extensions of the EFT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

CORRIGENDUM to Figures S3-S4 of "Quasi-stationary Planetary Scale Waves in Modern Climate Models"

In addition to fixing the color-coding legend, the new figures plot Fourier phase in a more easily recognizable format: polar coordinates rather than Cartesian coordinates. The same phase data is displayed—indicating the position of the first ridge to the east of the prime meridian—but the new plots avoid unphysical jumps, e.g. from +180° to -180° as longitude increases across the International Dateline. A slight disadvantage of the new format is that each hemisphere requires a separate plot. The wavenumber-1 Northern Hemisphere phase plots give, in effect, a view from above the North Pole, with Eurasia at the top and North America at the bottom. Wavenumber-1 Southern Hemisphere plots give an analogous view from the South Pole. Plots for higher wavenumbers, however, wrap around at fractions of 360° (180°, 120°, ...) and cannot be interpreted as a geographical view.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

L'Arlesienne de ROOT

Over many years, ROOT users have repeatedly stumbled over—and loudly rediscovered—the infamous 1 GB limit on individual I/O operations, a constraint that somehow survived long past the era when anyone thought a gigabyte was “a lot.” As experiments embraced ever-larger objects and collections, this limit became an increasingly unavoidable rite of passage. This contribution recounts the sustained, multi-year quest by ROOT I/O developers to finally retire this relic, navigating a maze of legacy APIs, memory-management assumptions, and integer boundaries that seemed determined to preserve the status quo. We describe how internal interfaces were carefully modernized to introduce fully 64-bit–capable code paths without breaking the mountains of existing user code that would definitely have noticed. With the limit now lifted, ROOT can finally handle multi-gigabyte objects in a single read or write operation, even when splitting them into an RNTuple is not an option (we’re looking at you, large RooWorkspaces and giant histograms), liberating users from yet another “fun” debugging adventure and clearing the way for the massive analyses of the HL-LHC and beyond.

Canal, Philippe G. [Fermilab] (ORCID:0000000277487↗
DOE OSTI2026

L'Arlesienne de ROOT

Over many years, ROOT users have repeatedly stumbled over—and loudly rediscovered—the infamous 1 GB limit on individual I/O operations, a constraint that somehow survived long past the era when anyone thought a gigabyte was “a lot.” As experiments embraced ever-larger objects and collections, this limit became an increasingly unavoidable rite of passage. This contribution recounts the sustained, multi-year quest by ROOT I/O developers to finally retire this relic, navigating a maze of legacy APIs, memory-management assumptions, and integer boundaries that seemed determined to preserve the status quo. We describe how internal interfaces were carefully modernized to introduce fully 64-bit–capable code paths without breaking the mountains of existing user code that would definitely have noticed. With the limit now lifted, ROOT can finally handle multi-gigabyte objects in a single read or write operation, even when splitting them into an RNTuple is not an option (we’re looking at you, large RooWorkspaces and giant histograms), liberating users from yet another “fun” debugging adventure and clearing the way for the massive analyses of the HL-LHC and beyond.

Canal, Philippe G. [Fermilab] (ORCID:0000000277487↗
NASA NTRS1987

Fifty Year Canon of Solar Eclipses: 1986 - 2035

A complete catalog is presented, listing the general characteristics of every solar eclipse from 1901 through 2100. To complement this catalog, a detailed set of cylindrical projection world maps shows the umbral paths of every solar eclipse over the 200 year interval. Focusing in on the next 50 years, accurate geodetic path coordinates and local circumstances for the 71 central eclipses from 1987 through 2035 are tabulated. Finally, the geodetic paths of the umbral and penumbral shadows of all 109 solar eclipses in this period are plotted on orthographic projection maps of the Earth. Appendices are included which discuss eclipse geometry, eclipse frequency and occurrence, modern eclipse prediction and time determination. Finally, code for a simple Fortran program is given to predict the occurrence and characteristics of solar eclipses.

Espenak, Fred↗
DOE OSTI2017

Development and Assessment of CTF for Pin-resolved BWR Modeling

CTF is the modernized and improved version of the subchannel code, COBRA-TF. It has been adopted by the Consortium for Advanced Simulation for Light Water Reactors (CASL) for subchannel analysis applications and thermal hydraulic feedback calculations in the Virtual Environment for Reactor Applications Core Simulator (VERA-CS). CTF is now jointly developed by Oak Ridge National Laboratory and North Carolina State University. Until now, CTF has been used for pressurized water reactor modeling and simulation in CASL, but in the future it will be extended to boiling water reactor designs. This required development activities to integrate the code into the VERA-CS workflow and to make it more ecient for full-core, pin resolved simulations. Additionally, there is a significant emphasis on producing high quality tools that follow a regimented software quality assurance plan in CASL. Part of this plan involves performing validation and verification assessments on the code that are easily repeatable and tied to specific code versions. This work has resulted in the CTF validation and verification matrix being expanded to include several two-phase flow experiments, including the General Electric 3 3 facility and the BWR Full-Size Fine Mesh Bundle Tests (BFBT). Comparisons with both experimental databases is reasonable, but the BFBT analysis reveals a tendency of CTF to overpredict void, especially in the slug flow regime. The execution of these tests is fully automated, analysis is documented in the CTF Validation and Verification manual, and the tests have become part of CASL continuous regression testing system. This paper will summarize these recent developments and some of the two-phase assessments that have been performed on CTF.

Salko, Robert K↗
DOE OSTI2018

ANALYSIS OF COMMUNICATION PERFORMANCE DEGRADATION OF THE RADIATION TRANSPORT CODE PIDOTS ON HIGH-UTILIZATION, MULTI-USER HPC SYSTEMS

The PIDOTS radiation transport code implements a spatially decomposed Integral Transport Matrix Method (ITMM) solver, which is intended to fully utilize the capabilities of modern, massively parallel high-performance computing (HPC) systems. While the code shows promising results, there was also an unexpected loss of parallel efficiency on the test systems as the number of participating processors grew and the parallelization grain size got finer. We seek to identify and characterize the communication-latency based effects that contribute to this slowdown. Through the creation of a high-level, parameterized representation of the workload of the code, the communication methodology was identified as the likely reason. From there, low-level, InfiniBand communication performance data measured on the Falcon HPC at Idaho National Laboratory (INL) was used to produce a discrete-event model of the communication scheme. The data, as well as the model, show that infrequent, large spikes in send/receive latency drive an unexpectedly substantial increase in projected runtime of a code. Furthermore, these effects are more pronounced under the high-utilization, scattered-workload environment of Falcon. Based on results from the model, it can be shown that a high-communication volume, small message-size regime can easily lead to slowdowns of ~5-10x, when compared to ideal results, as the processor count increases.

97 MATHEMATICS AND COMPUTING↗