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

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Analytical expression of a finite, long, conical canted-cosine-theta coil for particle collider interaction regions

Magnets in the accelerator interaction region (IR) present significant challenges because of high field requirements and limited available space. Conical-shaped magnets offer advantages in these environments by allowing closer placement to the interaction point while maintaining clearance from synchrotron radiation. Interestingly, numerical studies have shown that conical canted-cosine-theta (CCT) designs produce a constant field distribution along the axial direction in the IR quadrupoles for the Electron-Ion Collider (EIC) at Brookhaven National Laboratory. However, the field harmonics generated by conical CCT windings are not yet fully understood. This paper presents an analytical approach to describe the magnetic field produced by a conical surface current and proposes a method for designing conical CCT magnets for accelerator applications. First, we begin with a surface current sheet having a general cosine-theta distribution in spherical coordinates and solve the vector potential using the Green’s function. The magnetic fields generated by the conical current sheet are expressed using associated Legendre polynomials. These results are then related to circular field harmonics and integral field harmonics for designing a coil that produces a pure multipole field. Next, a single layer of the conical CCT winding path is produced based on the cosine-theta current distribution. Finally, the magnetic field quality of dipole and quadrupole conical CCT coils with multiple layers is verified using the Biot-Savart law.

Yang, Ye↗
DOE OSTI2026

Navigating the Path to Autonomy: Real-World Lessons from an Air-Free Self-Driving Laboratory

While autonomous experimentation has promise to accelerate discovery in physcial sciences, the real-world integration of predictive models and experimentation is non-trivial. Here we describe the genesis of a self-driving laboratory (SDL) for air-sensitive chemistry at Argonne National Laboratory and demonstrate the experimental design considerations needed for high-throughput experiments before predictive models can lead to scientific discovery. Our SDL was designed to explore battery electrolyte stability. Our final SDL utilized plate readers in a glovebox with a nitrogen atmosphere to perform kinetic assays and screen hundreds of battery-relevant solvents. However, the roadmap to autonomy and airfree-friendly experimentation required the complex evaluation of several spectroscopic and chromatographic methods. The greatest experimental challenges were (a) developing long-term sampling methods that remained air-free; (b) accelerating kinetics to advance reactivity projections; and (c) ensuring labware compatibility with nonaqueous solvents used in battery chemistry. Our experiences highlight the practical gap between closed-loop aspirations and the realities of chemical discovery, offering lessons on the challenges of transferring every day laboratory workflows to autonomy. These results suggest a more realistic blueprint for autonomy in chemistry—one that balances thoughtful and realistic experimental formulation.

Robertson, Lily A.↗
DOE OSTI2015

Simplified Predictive Models for CO 2 Sequestration Performance Assessment: Research Topical Report on Task #4 - Reduced-Order Method (ROM) Based Models

Reduced-order models provide a means for greatly accelerating the detailed simulations that will be required to manage CO 2 storage operations. In this work, we investigate the use of one such method, POD-TPWL, which has previously been shown to be effective in oil reservoir simulation problems. This method combines trajectory piecewise linearization (TPWL), in which the solution to a new (test) problem is represented through a linearization around the solution to a previously-simulated (training) problem, with proper orthogonal decomposition (POD), which enables solution states to be expressed in terms of a relatively small number of parameters. We describe the application of POD-TPWL for CO 2 -water systems simulated using a compositional procedure. Stanford’s Automatic Differentiation-based General Purpose Research Simulator (AD-GPRS) performs the full-order training simulations and provides the output (derivative matrices and system states) required by the POD-TPWL method. A new POD-TPWL capability introduced in this work is the use of horizontal injection wells that operate under rate (rather than bottom-hole pressure) control. Simulation results are presented for CO 2 injection into a synthetic aquifer and into a simplified model of the Mount Simon formation. Test cases involve the use of time-varying well controls that differ from those used in training runs. Results of reasonable accuracy are consistently achieved for relevant well quantities. Runtime speedups of around a factor of 370 relative to full- order AD-GPRS simulations are achieved, though the preprocessing needed for POD-TPWL model construction corresponds to the computational requirements for about 2.3 full-order simulation runs. A preliminary treatment for POD-TPWL modeling in which test cases differ from training runs in terms of geological parameters (rather than well controls) is also presented. Results in this case involve only small differences between training and test runs, though they do demonstrate that the approach is able to capture basic solution trends. The impact of some of the detailed numerical treatments within the POD-TPWL formulation is considered in an Appendix.

Mishra, Srikanta↗
DOE OSTI2016

SU-C-BRC-07: Parametrized GPU Accelerated Electron Monte Carlo Second Check

Purpose: I am presenting a parameterized 3D GPU accelerated electron Monte Carlo second check program. Method: I wrote the 3D grid dose calculation algorithm in CUDA and utilized an NVIDIA GeForce GTX 780 Ti to run all of the calculations. The electron path beyond the distal end of the cone is governed by four parameters: the amplitude of scattering (AMP), the mean and width of a Gaussian energy distribution (E and α), and the percentage of photons. In my code, I adjusted all parameters until the calculated PDD and profile fit the measured 10×10 open beam data within 1%/1mm. I then wrote a user interface for reading the DICOM treatment plan and images in Python. In order to verify the algorithm, I calculated 3D dose distributions on a variety of phantoms and geometries, and compared them with the Eclipse eMC calculations. I also calculated several patient specific dose distributions, including a nose and an ear. Finally, I compared my algorithm’s computation times to Eclipse’s. Results: The calculated MU for all of the investigated geometries agree with the TPS within the TG-114 action level of 5%. The MU for the nose was < 0.5 % different while the MU for the ear at 105 SSD was ∼2 %. Calculation times for a 12MeV 10×10 open beam ranged from 1 second for a 2.5 mm grid resolution with ∼15 million particles to 33 seconds on a 1 mm grid with ∼460 million particles. Eclipse calculation runtimes distributed over 10 FAS workers were 9 seconds to 15 minutes respectively. Conclusion: The GPU accelerated second check allows quick MU verification while accounting for patient specific geometry and heterogeneity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Green's function fast multipole method for continuum mechanics (SM-FMM)

Solid Mechanics Fast Multipole Method based on elastic Green's function and accelerated using FFTs The code calculates the mechanical fields (stress and strain) for a heterogeneous elasto-plastic material unit cell under quasi-static conditions (no dynamic effects). The fields are calculated using a Green's function method, where the strain field is given by a discrete convolution of Green's operator with an auxiliary stress field. The convolution is calculated using the fast multipole method.

Zecevic, Miroslav↗
DOE OSTI2017

Neutron interrogation systems using pyroelectric crystals and methods of preparation thereof

According to one embodiment, an apparatus includes a pyroelectric crystal, a deuterated or tritiated target, an ion source, and a common support coupled to the pyroelectric crystal, the deuterated or tritiated target, and the ion source. In another embodiment, a method includes producing a voltage of negative polarity on a surface of a deuterated or tritiated target in response to a temperature change of a pyroelectric crystal, pulsing a deuterium ion source to produce a deuterium ion beam, accelerating the deuterium ion beam to the deuterated or tritiated target to produce a neutron beam, and directing the ion beam onto the deuterated or tritiated target to make neutrons using a voltage of the pyroelectric crystal and/or an HGI surrounding the pyroelectric crystal. The directionality of the neutron beam is controlled by changing the accelerating voltage of the system. Other apparatuses and methods are presented as well.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Radial acceleration of ions during adiabatic expansion of a multicomponent cylindrical plasma

The methods of modern group analysis allow an analytic solution of the Cauchy problem to be constructed for the system of kinetic equations for a fully ionised electron – ion plasma, describing the acceleration of ions during the adiabatic expansion of a cylindrical plasma. Time and spatial dependences of the distribution functions of particles are obtained and their integral characteristics, such as density, average velocity, temperature, and energy spectrum, are found. The formation of the energy spectrum of accelerated ions, asymptotically repeating the spatial distribution of their density, and the cooling of electrons in the process of ion acceleration are analytically described. Particular attention is paid to the investigation of the influence of the heavy ionic component on the dynamics of the light component. The features of ion acceleration in the case of a two-temperature electron distribution function that describes the presence of hot and cold electron components are studied, which corresponds to the typical conditions of the experiment on plasma heating by intense laser radiation. (paper)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis Background & Noise in Stretched Wire Alignment Technique Measurements

The Stretched-Wire Alignment Technique (SWAT) is one method of magnet alignment for linear induction accelerators. The applications of SWAT have been implemented for aligning solenoid magnets on the Scorpius linear induction accelerator which will be sited at the Nevada National Security Site and the Flash X-Ray (FXR) linear induction accelerator at Lawrence Livermore National Laboratory’s Contained Firing Facility. This article describes both systematic (repeatable) and random sources of background and noise as well as practical ways to eliminate or reduce them to acceptable levels. Systematic sources include reflections from wire ends, rapid sag due to ohmic heating of the wire, magnetic materials, and shot rate. Random sources include air currents, vibration of nearby equipment, mechanical stability of test equipment, and the instruments used to measure the wire motion. Mitigations include curve fitting and adaptive noise signal cancellation, and mechanical damping. Finite Element Analysis (FEA) was used to identify and resolve a repeatable wire vibration frequency interfering with the signal resolution. Two stretched wire alignment technique set ups from Sandia National Labs and Lawrence Livermore National Lab have shown background noise sources and ways of mitigating them by either analysis methods or change of mechanical configuration. Conclusions that were drawn included the severe sensitivity of the deflection to even small external interferences of the SWAT wire such that it requires attention to detail in mechanical set up and analysis.

Linear Inductive Accelerator↗
DOE OSTI2024

Mechanical Performance of Additively Manufactured ODS 316L and 316H Stainless Steels

This work package within the Advanced Materials and Manufacturing Technologies (AMMT) program focused on the mechanical and microstructural characterization of oxide dispersion strengthened (ODS) stainless steels produced using advanced manufacturing techniques. The research aimed to identify an accelerated development path for ODS alloys by integrating additive manufacturing (AM) technologies with recent advancements in ODS materials and traditional manufacturing methods. For FY 24, the research specifically targeted exploring an accelerated development path for ODS alloys by combining AM technologies with these advancements. Novel AM and post-build processing routes were developed for ODS austenitic alloys, including Fe-Cr-Ni alloys such as 316L and 316H. Electron microscopy and mechanical characterizations were conducted to evaluate the impact of process variables on microstructure and properties, with the goal of optimizing these properties economically. Traditionally, producing ODS alloys involves multi-day, high-energy mechanical milling of alloy powder with yttria powder, followed by milled-powder consolidation through extrusion or other methods, and additional thermomechanical processing (TMP) for property control. To overcome the challenges associated with this complex and costly approach, we propose exploring alternative, cost-effective processing routes that emphasize AM and traditional TMP methods. The new processing routes for ODS alloys have achieved significantly higher strengths—several times greater than those of wrought stainless steels—while maintaining substantial ductility and fracture toughness. This report outlines the novel and economical AM-based processing routes for ODS austenitic alloys, combined with post-build TMPs, and discusses the mechanical and microstructural characteristics of the developed materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting System Response at Unmeasured Locations

Traditional techniques to derive dynamic specification for components have a great deal of uncertainty. One of the major sources of uncertainty is that the number of response measurements in the operational system environment is insufficient to determine the component motion. This inadequacy is due to logistical limitations for data recording in field testing and space limitations for accelerometers, strain gages and associated wiring. Available measurements are often some distance from the component and therefore do not represent component motion. Typical straight-line envelopes of these unrepresentative measurements guarantee an increase in the uncertainty. In this paper multiple methods are attempted to expand a sparse set of field test measurements on a system to responses of interest that cannot be measured in the field due to the limitations. Proof of concept is demonstrated on the Modal Analysis Test Vehicle (MATV). The responses of interest, known as “truth responses”, are measured in a system vibration environment along with an optimized sparse set of 30 field responses. Methods to expand the field responses to the truth responses are demonstrated by comparing the acceleration spectral density of the expanded response to the measured response. Two methods utilize a validated finite element model of the MATV. One is developed from purely experiment based frequency response functions of a laboratory pre-test. Finally, these approaches are designed to drastically reduce the uncertainty of the component in-service motion as a basis for developing specifications that are guaranteed to be conservative with a known (instead of unknown) conservatism.

42 ENGINEERING↗
DOE OSTI2023

10 Be at iThemba LABS using a silicon nitride membrane stack as absorber for isobar suppression

10 Be is an important isotope for accelerator mass spectrometry (AMS) because of the demand for cosmogenic radionuclide dating methods in the earth science and paleo-sciences community. At the iThemba Laboratory for Accelerator Based Science (iThemba LABS) we implemented full suppression of the interfering isobar 10 B using a silicon nitride foil-stack, utilizing the 2+ charge state for high efficiency. We demonstrate the performance of this newly established AMS system using standards and test samples. Here we further present the results of an inter-comparison between iThemba LABS and the Center for Accelerator Mass Spectrometry/Lawrence Livermore National Laboratory, on AMS samples prepared from purified quartz at the University of Vermont. The results for 10 Be from the laboratories are in close agreement, fully consistent with cross-calibration between them. AMS results for 26 Al are in similarly good agreement, demonstrating the performance and accuracy of iThemba LABS for the most commonly measured in situ produced cosmogenic nuclides.

10-Be↗
DOE OSTI2025

Reconstruction of beam parameters and betatron radiation spectra measured with a Compton spectrometer

The photon flux resulting from high-energy electron beam interactions with high-field systems, such as those found in the upcoming FACET-II experiments at the SLAC National Accelerator Laboratory, yields deep insight into the electron beam’s underlying dynamics during the interaction. However, extracting this information is an intricate process. To demonstrate how to approach this challenge using modern methods, this paper utilizes simulated data that models plasma wakefield acceleration-derived betatron radiation in experiments to determine reliable methods of reconstructing key beam and beam-plasma interaction properties. For betatron radiation measurements, translating the observed 200⁢ keV to 30⁢ MeV photon double-differential energy-angle spectra obtained from an advanced Compton spectrometer requires testing multiple methods to optimize the pipeline from its response to incident electron beam information. The paper compares maximum likelihood estimation and machine learning to refine the translation of photon spectra into precise electron beam metrics, such as spot size, energy, and emittance, enhancing the understanding of beam behavior within these dense, high-field environments. We also introduce machine learning and the expected maximization algorithm to reconstruct the primary photon spectrum, employing a multilayer neural network for regression analysis of the energy and angle spectra. With appropriate modifications, the advanced methods reproduce relevant incident beam parameters with high accuracy, even for beam sizes in the <10 μ⁢m range. This capacity is critical to understanding intense beam propagation and its optimization in plasma.

Beam code development & simulation techniques↗
DOE OSTI2025

Hardware acceleration for HPS algorithms in two and three dimensions

We provide a flexible, open-source framework for hardware acceleration, namely massively-parallel execution on general-purpose graphics processing units (GPUs), applied to the hierarchical Poincaré–Steklov (HPS) family of algorithms for building fast direct solvers for linear elliptic partial differential equations. To take full advantage of the power of hardware acceleration, we propose two variants of HPS algorithms to improve performance on two- and three-dimensional problems. In the two-dimensional setting, we introduce a novel recomputation strategy that minimizes costly data transfers to and from the GPU; in three dimensions, we modify and extend the adaptive discretization technique of Geldermans and Gillman [1] to greatly reduce peak memory usage. We provide an open-source implementation of these methods written in JAX, a high-level accelerated linear algebra package, which allows for the first integration of a high-order fast direct solver with automatic differentiation tools. We conclude with extensive numerical examples showing our methods are fast and accurate on two- and three-dimensional problems.

Fast direct solvers↗
DOE OSTI2021

Galactic Potential and Dark Matter Density from Angular Stellar Accelerations

Here, we present an approach to measure the Milky Way (MW) potential using the angular accelerations of stars in aggregate as measured by astrometric surveys like Gaia. Accelerations directly probe the gradient of the MW potential, as opposed to indirect methods using, e.g., stellar velocities. We show that end-of-mission Gaia stellar acceleration data may be used to measure the potential of the MW disk at approximately 3σ significance and, if recent measurements of the solar acceleration are included, the local dark matter density at ~2σ significance. Since the significance of detection scales steeply as t 5/2 for observing time t, future surveys that include angular accelerations in the astrometric solutions may be combined with Gaia to precisely measure the local dark matter density and shape of the density profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Excitation of {sup 111m}Cd, {sup 113m}In, and {sup 115m}In isomeric states by photons of energy up to 8 MeV

Experimental results on the excitation of isomeric nuclear states in the reactions {sup 111}Cd(γ, γ′){sup 111m}Cd, {sup 113}In(γ, γ′){sup 113m}In, and {sup 115}In(γ, γ′){sup 115m}In studied by the activation method at photon energies of up to 8 MeV are presented. The respective experiments were performed at the LUE-8 linear electron accelerator of the Institute for Nuclear Research (INR, Moscow). A method for calculating cross sections per equivalent photon was developed with the aid of a simulation on the basis of the GEANT-4 code. The data obtained in this way are discussed from the point of view of the possibility of studying collective nuclear excitations of various multipolarity near the threshold.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

SU-F-T-547: Off-Isocenter Winston-Lutz Test for Stereotactic Radiosurgery/stereotactic Body Radiotherapy

Purpose: To perform a quantitative study to verify that the mechanical field center coincides with the radiation field center when both are off from the isocenter during the single-isocenter technique in linear accelerator-based SRS/SBRT procedure to treat multiple lesions. Methods: We developed an innovative method to measure this accuracy, called the off-isocenter Winston-Lutz test, and here we provide a practical clinical guideline to implement this technique. We used ImagePro V.6 to analyze images of a Winston-Lutz phantom obtained using a Varian 21EX linear accelerator with an electronic portal imaging device, set up as for single-isocenter SRS/SBRT for multiple lesions. We investigated asymmetry field centers that were 3 cm and 5 cm away from the isocenter, as well as performing the standard Winston-Lutz test. We used a special beam configuration to acquire images while avoiding collision, and we investigated both jaw and multileaf collimation. Results: For the jaw collimator setting, at 3 cm off-isocenter, the mechanical field deviated from the radiation field by about 2.5 mm; at 5 cm, the deviation was above 3 mm, up to 4.27 mm. For the multileaf collimator setting, at 3 cm off-isocenter, the deviation was below 1 mm; at 5 cm, the deviation was above 1 mm, up to 1.72 mm, which is 72% higher than the tolerance threshold. Conclusion: These results indicated that the further the asymmetry field center is from the machine isocenter, the larger the deviation of the mechanical field from the radiation field, and the distance between the center of the asymmetry field and the isocenter should not exceed 3 cm in of our clinic. We recommend that every clinic that uses linear accelerator, multileaf collimator-based SRS/SBRT perform the off-isocenter Winston-Lutz test in addition to the standard Winston-Lutz test and use their own deviation data to design the treatment plan.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

HyKKT: a hybrid direct-iterative method for solving KKT linear systems

Here, we propose a solution strategy for the large indefinite linear systems arising in interior methods for nonlinear optimization. The method is suitable for implementation on hardware accelerators such as graphical processing units (GPUs). The current gold standard for sparse indefinite systems is the LBLT factorization where L is a lower triangular matrix and B is 1×1 or 2×2 block diagonal. However, this requires pivoting, which substantially increases communication cost and degrades performance on GPUs. Our approach solves a large indefinite system by solving multiple smaller positive definite systems, using an iterative solver on the Schur complement and an inner direct solve (via Cholesky factorization) within each iteration. Cholesky is stable without pivoting, thereby reducing communication and allowing reuse of the symbolic factorization. We demonstrate the practicality of our approach on large optimal power flow problems and show that it can efficiently utilize GPUs and outperform LBL T factorization of the full system.

97 MATHEMATICS AND COMPUTING↗