Search NASA⌕ Search

SEARCH · Search NASA

Results for “Acceleration methods”

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,351 records · Page 75

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

An Investigation of Energy Transmission Due to Flexural Wave Propagation in Lightweight, Built-Up Structures

A technique to measure flexural structure-borne noise intensity is investigated. Two accelerometers serve as transducers in this cross-spectral technique. The structure-borne sound power is obtained by two different techniques and compared. In the first method, a contour integral of intensity is performed from the values provided by the two-accelerometer intensity technique. In the second method, input power is calculated directly from the output of force and acceleration transducers. A plate and two beams were the subjects of the sound power comparisons. Excitation for the structures was either band-limited white noise or a deterministic signal similar to a swept sine. The two-accelerometer method was found to be sharply limited by near field and transducer spacing limitations. In addition, for the lightweight structures investigated, it was found that the probe inertia can have a significant influence on the power input to the structure. In addition to the experimental investigation of structure-borne sound energy, an extensive study of the point harmonically forced, point-damped beam boundary value problem was performed to gain insight into measurements of this nature. The intensity formulations were also incorporated into the finite element method. Intensity mappings were obtained analytically via finite element modeling of simple structures.

Mickol, John Douglas↗

On acceleration of MacCormack's scheme

An acceleration of MacCormack's scheme due to Desideri and Tannehill is analyzed. It is found that for hyperbolic problems one cannot improve upon the efficiency of MacCormack's method. For parabolic problems the time step can be chosen arbitrarily large without loss of stability by an appropriate choice of the acceleration parameters. When applied to the heat equation this method is equivalent to both the Dufort-Frankel scheme and to MSOR.

Gottlieb, D.↗

Evaluation of ASCE 4-16 and AISC 43-18 (Draft) for use in the Risk-Informed Performance-Based Seismic Design of Nuclear Power Plant Structures, Systems, and Components

This report describes the assessment of ASCE Standards 4-16 and 43-18 (Draft) for use in the Risk-Informed Performance-Based (RIPB) seismic design of structures, systems, and components (SSCs) at nuclear power plants. This work was performed for the U.S. Nuclear Regulatory Commission (NRC) Office of Regulatory Research (RES), to support potential endorsement of these industry standards for the design of nuclear power plants based on the RIPB approach. Currently, the NRC endorses a deterministic design approach for demonstrating the design adequacy of SSCs based on the Standard Review Plan (NUREG-0800) and NRC Regulatory Guides. In the RIPB approach, the design criteria are developed to achieve a target performance goal, which is defined by the annual frequency of occurrence of the design basis earthquake (i.e., Seismic Design Category) and the acceptable level of structural performance (i.e., Limit State) for the SSCs. ASCE 4-16 provides methods for performing a seismic analysis of structures to obtain the seismic response of these structures (e.g., building displacements, accelerations, in-structure response spectra) which are used in the design of the SSCs. It also provides methods for performing seismic analysis of SSCs to determine the seismic demands (e.g. member forces and displacements) needed to design individual SSCs. ASCE 43-18 (Draft) provides the criteria for the seismic design of SSCs using the seismic demands developed in ASCE 4-16. The use of ASCE Standard 43-18 (Draft), along with ASCE 4-16, provides the criteria for the seismic design of the SSCs. ASCE 43-18 (DRAFT) relies on other consensus codes and standards such as ACI 349 for reinforced concrete, AISC/N690 for steel structures, ASME Section III for pressure-retaining mechanical components and Containments, and IEEE-344 for Class 1E equipment. The goal of this technical review was to assess the adequacy of the provisions in these standards for use by the NRC in developing regulatory guidance for design of SSCs in nuclear power plants, based on the RIPB approach. The research reported herein describes the basis for acceptance of the new standards and identifies areas where additional staff guidance is needed for the seismic design of SSCs at nuclear power plants. This technical review has determined that ASCE 4-16 and ASCE 43-18 (Draft) provide an appropriate framework for the seismic design of SSCs at nuclear power plants using a Risk-Informed Performance-Based approach. However, some of the criteria therein warrant exceptions, qualifications, and/or clarifications.

42 ENGINEERING↗

Rapid Alloy Surface Engineering through Closed-Vessel Reagent Pyrolysis

For rapid surface engineering of Cr-containing alloys by low-temperature nitrocarburization, we introduce a process based on pyrolysis of solid reagents, e.g., urea, performed in an evacuated closed vessel. Upon heating to temperatures high enough for rapid diffusion of interstitial solute, but low enough to avoid second-phase precipitation, the reagent is pyrolyzed to a gas atmosphere containing molecules that (i) activate the alloy surface by stripping away the passivating Cr 2 O 3 -rich surface film (diffusion barrier) and (ii) rapidly infuse carbon and nitrogen into the alloy. We demonstrate quantitatively that this method can generate a subsurface zone with concentrated carbon and nitrogen comparable to what can be accomplished by established (e.g., gas-phase- or plasma-based) methods, but with significantly reduced processing time. As another important difference to established gas-phase processing, the interaction of gas molecules with the alloy surface can have auto-catalytic effects by altering the gas composition in a way that accelerates solute infusion by providing a high activity of HNCO. The new method lends itself to rapid experimentation with a minimum of laboratory equipment.

36 MATERIALS SCIENCE↗

Hardware-Accelerated Ray Tracing of CAD-Based Geometry for Monte Carlo Radiation Transport

Monte Carlo radiation transport (MCRT) methods have been used to simulate radiation environments for many decades by tracking individual particles through a model to accumulate statistical information. MCRT geometry is historically formed using the constructive solid geometry (CSG). Recently, significant work has been performed to support simulations using computer-aided design (CAD)-based tessellated surfaces to support highly complex geometries. Ray tracing acceleration data structures from the rendering and visualization community are applied to accelerate particle tracking in CAD-based models. Despite these efforts, CSG representations provide the superior performance in surface intersection operations during particle flight. Concurrently, pseudo Monte Carlo methods have become prevalent in rendering applications to support more realistic models for scattering media, motivating innovations that are advantageous for MCRT simulations. Finally, the authors’ work extends these innovations by employing Intel’s Embree ray tracing kernel within a geometry toolkit for Monte Carlo to improve the simulation performance using CAD-based models by factors of 1.5 to 2.

42 ENGINEERING↗

Early-time γ -ray constraints on cosmic-ray acceleration in the core-collapse SN 2023ixf with the Fermi Large Area Telescope

Context. While supernova remnants (SNRs) have been considered the most relevant Galactic cosmic ray (CR) accelerators for decades, core-collapse supernovae (CCSNe) could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV γ -rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. Aims. The multi-wavelength observation of SN 2023ixf, a Type II supernova (SN) in the nearby galaxy M 101 (at a distance of 6.85 Mpc), opens the possibility to constrain CR acceleration within a few days after the collapse of the red super-giant stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. Methods. We performed a maximum-likelihood analysis of γ -ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy, non-thermal emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. Results. We do not find significant γ -ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Conclusions. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. Consequently, a more sophisticated, inhomogeneous characterisation of the shock and the progenitor’s environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times.

79 ASTRONOMY AND ASTROPHYSICS↗

An Accelerated Creep Testing (ACT) Program for Advanced Creep Resistant Alloys for High Temperature Fossil Energy (FE) Applications (Final Report)

The research objective of this study is to develop an accelerated creep testing (ACT) program for advanced creep resistant alloys for high temperature FE applications. This research produced a new ACT for metallic materials based on Time-Temperature-Stress Superposition (TTSSP). In TTSSP, there exists a fundamental relationship between creep resistance, time, temperature, and stress that can be manipulated to reduce test durations. In our study, the stepped isostress method was employed where step increases of load during creep testing accelerated the time-to-rupture. Well-posed constitutive models including sin-hyperbolic and Wilshire-Cano-Stewart where calibrated using the SSM data and employed to predict the conventional creep response of the material. In our best experiments, <100-hour tests on lnconel 718 at 750 deg C produced over a 64x acceleration of data when compared to conventional data from Japan's National Institute of Materials Science. The ACTs enable the rapid screening of candidate materials, where a small test matrix of short duration tests can be employed to determine the creep resistance (minimum-creep-strain-rate, stress-rupture, and creep deformation) of a material across decades.

20 FOSSIL-FUELED POWER PLANTS↗

Laser ion-shotgun acceleration

Relativistic transparency induced by an ultra-intense laser ($5 \times 10^{22}$ W/cm$^2$) offers an attractive method for generating extreme electric (50 TV/m) and magnetic (0.5 MT) fields by enabling laser pulse propagation through a very dense and otherwise opaque material. A slowly evolving structure that effectively accelerates ions is generated by employing a target ($\sim 10$ $\mu$m thick) with a pre-filled channel that guides the laser pulse. As the pulse exits the channel driving a strong electron current, strong accelerating longitudinal and focusing transverse electric fields emerge at the rear side of the target. Three-dimensional kinetic simulations predict a low-diverging dense mono-energetic proton bunch with energies of 200-300 MeV. The charge in an opening angle of $10^{\circ}$ is tens of nC. This scheme potentially allows high repetition rate operation and relaxes laser contrast requirements.

47 OTHER INSTRUMENTATION↗

A Positive-Definite, WENO-Limited, High-Order Finite Volume Solver for 2-D Transport on the Cubed Sphere Using an ADER Time Discretization

Modern computer architectures reward added computation if it reduces algorithmic dependence, reduces data movement, increases accuracy/robustness, and improves memory accesses. The driving motive for this study is to develop a numerical algorithm that respects these constraints while improving accuracy and robustness. This study introduces the ADER-DT (Arbitrary DERivatives in time and space-differential transform) time discretization to positive-definite, weighted essentially nonoscillatory (WENO)-limited, finite volume transport on the cubed sphere in lieu of semidiscrete integrators. The cost of the ADER-DT algorithm is significantly improved from previous implementations without affecting accuracy. A new function-based WENO implementation is also detailed for use with the ADER-DT time discretization. While ADER-DT costs about 1.5 times more than a fourth-order, five-stage strong stability preserving Runge-Kutta (SSPRK4) method, it is far more computationally dense (which is advantageous on accelerators such as graphics processing units), and it has a larger effective maximum stable time step. ADER-DT errors converge more quickly with grid refinement than SSPRK4, giving 6.5 times less error in the L ∞ norm than SSPRK4 at the highest refinement level for smooth data. For nonsmooth data, ADER-DT resolves C 0 discontinuities more sharply. For a complex flow field, ADER exhibits less phase error than SSPRK4. In conclusion, improving both accuracy and robustness as well as better respecting modern computational efficiency requirements, we believe the method presented herein is competitive for efficiently transporting tracers over the sphere for applications targeting modern computing architectures.

97 MATHEMATICS AND COMPUTING↗

Multi-Pass Sequential Mini-Batch Stochastic Gradient Descent Algorithms for Noise Covariance Estimation in Adaptive Kalman Filtering

Estimation of unknown noise covariances in a Kalman filter is a problem of significant practical interest in a wide array of applications. Although this problem has a long history, reliable algorithms for their estimation were scant, and necessary and sufficient conditions for identifiability of the covariances were in dispute until recently. Necessary and sufficient conditions for covariance estimation and a batch estimation algorithm were presented in our previous study. This paper presents stochastic gradient descent algorithms for noise covariance estimation in adaptive Kalman filters that are an order of magnitude faster than the batch method for similar or better root mean square error. More significantly, these algorithms are applicable to non-stationary systems where the noise covariances can occasionally jump up or down by an unknown magnitude. The computational efficiency of the new algorithms stems from adaptive thresholds for convergence, recursive fading memory estimation of the sample cross-correlations of the innovations, and accelerated stochastic gradient descent algorithms. The comparative evaluation of the proposed methods on a number of test cases demonstrates their computational efficiency and accuracy.

Adaptive Kalman filtering↗

Adaptive machine learning for time-varying systems: low dimensional latent space tuning

Machine learning (ML) tools such as encoder-decoder convolutional neural networks (CNN) can represent incredibly complex nonlinear functions which map between combinations of images and scalars. For example, CNNs can be used to map combinations of accelerator parameters and images which are 2D projections of the 6D phase space distributions of charged particle beams as they are transported between various particle accelerator locations. Despite their strengths, applying ML to time-varying systems, or systems with shifting distributions, is an open problem, especially for large systems for which collecting new data for re-training is impractical or interrupts operations. Particle accelerators are one example of large time-varying systems for which collecting detailed training data requires lengthy dedicated beam measurements which may no longer be available during regular operations. We present a novel method of adaptive ML for time-varying systems. Our approach is to map very high (N ≈ 100k) dimensional inputs (a combination of scalar parameters and images) into the low dimensional (N ≈ 2) latent space at the output of the encoder section of an encoder-decoder CNN. We then actively tune the low dimensional latent space-based representation of complex system dynamics by the addition of an adaptively tuned feedback vector directly before the decoder sections builds back up to our image-based high-dimensional phase space density representations. This method allows us to learn correlations within and to quickly tune the characteristics of incredibly large parameter space systems and to track their evolution in real time based on feedback without massive new data sets for re-training. We demonstrate that our method can accurately predict and track the phase space of charged particle beams at various locations in a particle accelerator by adaptively adjusting in real-time while the unknown input beam distribution of the accelerator is changing in shape, charge, and offset and while the RF system of the accelerator itself is also changing in an unpredictable way. For FACET-II we demonstrate that such an approach has the potential to use transverse deflecting cavity and energy spread spectrum beam measurements to accurately predict 2D projections of the 6D phase space of the electron beam at the plasma wakefield acceleration interaction point where such diagnostics are unavailable.

47 OTHER INSTRUMENTATION↗

Particle-in-cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Abstract Relativistic magnetic reconnection is a nonideal plasma process that is a source of nonthermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a nonstandard finite-difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that, for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than that of Yee, resulting in a ∼40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss’s law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4–3 times over the standard Yee step still yields accurate results, but it only translates to modest performance improvements over CKC, due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

79 ASTRONOMY AND ASTROPHYSICS↗

Application of Particle Accelerators to Mitigate Energy and Climate Change Problems Facing America

Reliable CO 2 -free baseload power is needed to address ever-increasing demands for electricity while minimizing adverse climate change. Diversified power supply provided by solar, wind, geothermal, and nuclear reactors can displace the use of fossil fuels. The United States (U.S.) is taking a new look at nuclear power as a source of electrical energy and the nuclear industry is proposing new approaches which may minimize capital costs. However, nuclear power comes with a variety of technical problems. Chief among those is managing the used fuel from nuclear reactors. No long-term, practicable solution to this problem is available. The lack of progress on a comprehensive waste management strategy restricts growth in the nuclear power industry and minimizes the role that nuclear energy may serve as part of a zero-carbon future (Bahr, 2021). According to the Nuclear Waste Policy Act of 1982, as amended, the U.S. Government has possession of the used reactor fuel and incurs large annual storage fees paid to utilities to store and safeguard the accumulated used fuel. In effect, short term on-site storage of nuclear waste is the current waste management plan. The amount of used reactor fuel in the U.S. is approximately 80,000 metric tons. While no geologic repository exists within the U.S., the potential site at Yucca Mountain, NV would accommodate 70,000 metric tons, which fails to meet current and future needs. Given the technical and political challenges associated with establishing a single geological storage site, it is necessary for the U.S. to implement technologies that improve the suitability of geological storage by reducing the volume and radiotoxicity of stored material. No extant technology meets this need. Accelerator driven waste burners have been proposed as a scalable method to process used nuclear fuel, however considerable technical challenges remain which impede commercialization. Innovations in design require investigation of novel materials, development of accurate models and simulations, and a comprehensive assessment of safety and performance. This proposal elaborates how LANL is uniquely positioned to make key contributions to this area of research.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Virtual Engineering of Low-Temperature Conversion

In this work, we present the development of an overarching software framework and supporting multiphysics models to simulate the end-to-end process of biomass conversion. This virtual engineering (VE) software is designed with the goal of accelerating research and development and reducing risk for market-relevant biomass conversion processes. We currently support multiple models, computing paradigms, and fidelities representing the steps of feedstock pretreatment, enzymatic hydrolysis, and bioconversion. Although this VE approach was developed to support a biomass workflow, we have designed each component in a way that allows us to easily support new domains, unit models, and feedstocks. We begin by presenting the user-facing aspects of the VE software and highlight how simulated elements are defined and linked by VE functions. We then present an overview of the high-fidelity computational fluid dynamics (CFD) models developed to support our target domain before segueing into our efforts to develop accurate and fast surrogate models, capturing the salient outcomes from the CFD simulations in a significantly less computationally demanding manner. We then present how VE calculations interface with a commercial techno-economic analysis software, Aspen Plus. We conclude by presenting VE case studies that leverage these methods and discuss how our methods can be extended to support a wide variety of accelerated biofuel commercialization pathways in the future.

biofuel↗

Brochure for the DOE Office of Science Workshop on Envisioning Frontiers in AI and Computing for Biological Research

In February of 2025 a joint ASCR/BER workshop was held to identify key transformational research directions for understanding biology using artificial intelligence (AI), digital twins and high-performance (HPC) computational methods to facilitate scientific discovery and innovation in support of the Department of Energy mission. AI technologies offer exciting new groundbreaking methods to analyze large volumes of complex biological data, thereby greatly accelerating the ability to understand, predict, and design biological processes for beneficial purposes. In the laboratory, the bridging of AI-enabled automated experimental technologies, HPC and digital twins will provide potent tools for researchers to explore the fundamental nature of biology and harness its inherent metabolic potential for a variety of beneficial purposes. The focus of this workshop was on how high-performance computational methods can impact this objective by exploring digital twins, foundational models, and data-driven approaches with applications to advance automated laboratory experiments, modeling of complex living systems and engineering new functions into plants and microbial systems relevant to DOE mission. Workshop attendees with expertise in plant science, microbiology, mathematics, computer science, and AI assessed the current state of the science, trends, and AI challenges at the interface of plant and microbial systems biology and computational science to identify opportunities for high-impact research. This collaborative effort capitalized on ASCR's advancements in applied mathematics, computer science, and Exascale systems, and BER's expertise in basic genomics-enabled research on DOE relevant plant and microbial systems. The workshop culminated in four key priority research directions to guide future research and development within DOE Office of Science programs.

59 BASIC BIOLOGICAL SCIENCES↗

Independent Confirmation of the Pioneer 10 Anomalous Acceleration

I perform an independent analysis of radio Doppler tracking data from the Pioneer 10 spacecraft for the time period 1987-1994. All of the tracking data were taken from public archive sources, and the analysis tools were developed independently by myself. I confirm that an apparent anomalous acceleration is acting on the Pioneer 10 spacecraft, which is not accounted for by present physical models of spacecraft navigation. My best fit value for the acceleration, including corrections for systematic biases and uncertainties, is (8.60 plus or minus 1.34) x 10(exp -8) centimeters per second, directed towards the Sun. This value compares favorably to previous results. I examine the robustness of my result to various perturbations of the analysis method, and find agreement to within plus or minus 5%. The anomalous acceleration is reasonably constant with time, with a characteristic variation time scale of greater than 70 yr. Such a variation timescale is still too short to rule out on-board thermal radiation effects, based on this particular Pioneer 10 data set.

Markwardt, Craig B.↗

Time-temperature-superposition analysis of diverse datasets by the minimum-arclength method: long-term prediction with uncertainty margins

In a recent publication, we carried out an extensive analysis of an unsupervised method of determining optimum shift factors in time-temperature-superposition of accelerated-aging data that involves minimizing the vertical arclength to obtain the master curve. For synthetic Arrhenius data with a variety of noise distributions, the work showed that, in conjunction with bootstrap-resampling, the method can produce reliable estimates of the mean activation energy along with uncertainty quantification. Here, we apply the above method to six different datasets taken from the published literature and demonstrate accurate prediction of mean activation energy from the data as-is without the need for any pre-processing or fitting. We also compare uncertainty margins computed by second-order bootstrap with that by linear regression theory and show that the former appears to provide consistent margins in the presence of common noise types in real data, including intra-isotherm measurement-errors, sample-to-sample variations, and intrinsic deviation from perfect Arrhenius behavior.

36 MATERIALS SCIENCE↗