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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Outer Radiation Belt Dynamics During the October 2012 Storm Revisited: Rapid Inward Radial Transport From a Dynamic Outer Boundary

Earth's outer radiation belt electron flux is highly variable and can be enhanced by over an order of magnitude over timescales less than one day, as observed during the October 2012 storm. Previous studies of this storm (e.g., Reeves et al., 2013, https://doi.org/10.1126/science.1237743) have invoked local acceleration to explain this. However, here, we argue that the observations can instead be explained by fast inward radial transport. One method often invoked to distinguish between these two acceleration processes is the existence of local peaks in electron phase space density (PSD) as a function of L* at fixed first, M, and second, K, adiabatic invariants. However, this method relies on the assumption that the evolution of the PSD as a function of L* occurs over timescales slower than the satellite orbital period. Here, high spatiotemporal resolution data from the Global Positioning System (GPS) spacecraft constellation is used to show that enhancements in the PSD occur during the October 2012 storm over short timescales not resolvable by the Van Allen Probes. In addition, Geostationary Operational Environmental Satellite spacecraft data also indicate that these enhancements are consistent with relativistic electron injections. A radial diffusion model is shown to reproduce the PSD dynamics observed by the Van Allen Probes, once rapid variations at the simulation outer boundary are included, consistent with GPS data. This verifies that apparently “locally growing” peaks in PSD along high apogee satellite orbits can be produced by fast inward radial transport without requiring the action of any local acceleration processes.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Investigations on the multiple-sector hard-copper X-band accelerating structures

The development of advanced, high gradient accelerating structures is one of the leading activity of the particle accelerator community. In the technological research of new construction methods for these devices, high-power testing is a critical step for the verification of their viability. Recent experiments showed that accelerating cavities made out of hard copper, fabricated without high-temperature processes, can achieve better performance as compared with soft copper ones. Recently, we have built cavities using Tungsten Inert Gas welding and the high-power experiments confirmed that this joining process is a robust and low-cost alternative to brazing or diffusion bonding. This is a good solution for high-gradient operation, with a gradient of about 150 MV/m in X-band, at a breakdown rate of 10 -3 /pulse/meter using a shaped RF pulse with a 150 ns flat part. We continue the design, construction and high power tests of three-cells standing-wave X-band accelerating cavities fabricated with a split-open geometry, made of hard copper and vacuum-sealed by welding. Our aim is the fabrication of accelerating structures made out of hard copper alloys by using innovative cost-effective technologies. Moreover, our method of multiple-sector structures opens the way to new technological approaches and design methods, for example providing the possibility of placing dampers of the parasitic higher-order modes. Such structures with dampers allow for acceleration of multi-bunch beams, both in standing and in travelling-wave accelerating modality. Here, this paper describes the design of two cavities made of multiple sectors that were fabricated at the Italian Institute for Nuclear Physics in Frascati, Italy (INFN-LNF), assembled by means of clamping and then joined with the Tungsten Inert Gas welding in order to preserve the hardness of the metal. These are three-cell standing-wave X-band accelerating cavities, to be operated in the π-mode. We report the low-power RF tests of the first one, made of two halves, and the RF design of the second one, made of four quadrants. Both structures were built for high-gradient tests at SLAC National Accelerator Laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics

Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

Liu, Wen‐Hao↗
DOE OSTI2024

Democratizing uncertainty quantification

Uncertainty Quantification (UQ) is vital to safety-critical model-based analyses, but the widespread adoption of sophisticated UQ methods is limited by technical complexity. In this paper, we introduce UM-Bridge (the UQ and Modeling Bridge), a high-level abstraction and software protocol that facilitates universal interoperability of UQ software with simulation codes. It breaks down the technical complexity of advanced UQ applications and enables separation of concerns between experts. UM-Bridge democratizes UQ by allowing effective interdisciplinary collaboration, accelerating the development of advanced UQ methods, and making it easy to perform UQ analyses from prototype to High Performance Computing (HPC) scale. In addition, we present a library of ready-to-run UQ benchmark problems, all easily accessible through UM-Bridge. These benchmarks support UQ methodology research, enabling reproducible performance comparisons. We demonstrate UM-Bridge with several scientific applications, harnessing HPC resources even using UQ codes not designed with HPC support.

Benchmarks↗
DOE OSTI2024

A fast, matrix-based method to perform omnidirectional pressure integration

Abstract Experimentally-measured pressure fields play an important role in understanding many fluid dynamics problems. Unfortunately, pressure fields are difficult to measure directly with non-invasive, spatially resolved diagnostics, and calculations of pressure from velocity have proven sensitive to error in the data. Omnidirectional line integration methods are usually more accurate and robust to these effects as compared to implicit Poisson equations, but have seen slower uptake due to the higher computational and memory costs, particularly in 3D domains. This paper demonstrates how omnidirectional line integration approaches can be converted to a matrix inversion problem. This novel formulation uses an iterative approach so that the boundary conditions are updated each step, preserving the convergence behavior of omnidirectional schemes while also keeping the computational efficiency of Poisson solvers. This method is implemented in Matlab and also as a GPU-accelerated code in CUDA-C++. The behavior of the new method is demonstrated on 2D and 3D synthetic and experimental data. Three-dimensional grid sizes of up to 125 million grid points are tractable with this method, opening exciting opportunities to perform volumetric pressure field estimation from 3D PIV measurements.

42 ENGINEERING↗
DOE OSTI2021

Neutron generation using pyroelectric crystals

According to one embodiment, a method for producing a directed neutron beam includes producing a voltage of negative polarity of at least −100 keV on a surface of a deuterated or tritiated target in response to a temperature change of a pyroelectric crystal of less than about 40° C., the pyroelectric crystal having the deuterated or tritiated target coupled thereto, 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 at least one of a voltage of the pyroelectric crystal, and a high gradient insulator (HGI) surrounding the pyroelectric crystal. The directionality of the neutron beam is controlled by changing the accelerating voltage of the system. Other methods are presented as well.

Tang, Vincent↗
DOE OSTI2023

Accelerated fission rate irradiation design, pre-irradiation characterization, and adaptation of conventional PIE methods for U-10Mo and U-17Mo

Metallic U alloys have high U density and thermal conductivity and thus have been explored since the beginning of nuclear power research. Alloys of U with modest amounts of Mo, such as U-10 wt % Mo (U-10Mo), are of particular interest because the γ-U crystal structure in this alloying addition shows prolonged stability in reactor service. Historically, radiation data on U-10Mo fuels were collected in Na fast reactors or lower temperature research reactor conditions, but little is known about irradiation behavior, particularly swelling and creep, at irradiation temperatures between 250 and 500°C. This work discusses the methodology and pre-irradiation characterization results from a U-Mo irradiation campaign performed in the High Flux Isotope Reactor at Oak Ridge National Laboratory. U-10Mo and U-17Mo samples irradiations are being completed at temperatures ranging from 250 to 500°C to three targeted fission densities between 2 × 10 20 and 1.5 × 10 21 fissions per cubic centimeter. Swelling measurement of the specimen sizes studied here required development and assessment of new methods for volume determination before and after irradiation. Laser profilometry and X-ray computation tomography (XCT) were used to provide preirradiation characterization of samples to determine the error and applicability of each to determine swelling following irradiation. These outcomes are contextualized through use of BISON simulations performed to assess the predicted expansion of U-Mo fuels subjected to the irradiation conditions of this work. Use of existing BISON fuel performance models predicted a maximum of 7% swelling under the irradiation conditions of this study. Pre-irradiation characterization revealed the as-cast U-Mo fuel samples were uniformly large-grained fully cubic U crystals with small U-C/N bearing precipitates and pores distributed throughout. Samples were found to contain a bulk porosity between .4 and 3% because of the casting process. Local porosity in areas far from large, interconnected pores was found by Slice-and-View to be under .2%. Nanometer-sized precipitates rich in C and N were identified in all samples, likely because of impurities during the fabrication process. Dendritic bands were also observed throughout the samples. These bands were characterized by variable Mo content that deviated from the overall Mo content by 2–3 wt %. No other microstructural features were correlated to these bands. Mechanical properties were found to be slightly strengthened compared to literature reports of bulk U-Mo fuels due to the nano-scale precipitates throughout the sample.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Machine-learning-accelerated multimodal characterization and multiobjective design optimization of natural porous materials

Natural porous materials such as nanoporous clays are used as green and low-cost adsorbents and catalysts. The key factors determining their performance in these applications are the pore morphology and surface activity, which are typically represented by properties such as specific surface area, pore volume, micropore content and pH. The latter may be modified and tuned to specific applications through material processing and/or chemical treatment. Characterization of the material, raw or processed, is typically performed experimentally, which can become costly especially in the context of tuning of the properties towards specific application requirements and needing numerous experiments. In this work, we present an application of tree-based machine learning methods trained on experimental datasets to accelerate the characterization of natural porous materials. The resulting models allow reliable prediction of the outcomes of experimental characterization of processed materials (R2 from 0.78 to 0.99) as well as identification of key factors contributing to those properties through feature importance analysis. Furthermore, the high throughput of the models enables exploration of processing parameter–property correlations and multiobjective optimization of prototype materials towards specific applications. We have applied these methodologies to pinpoint and rationalize optimal processing conditions for clays exploitable in acid catalysis. One of such identified materials was synthesized and tested revealing appreciable acid character improvement with respect to the pristine material. Specifically, it achieved 79% removal of chlorophyll-a in acid catalyzed degradation.

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

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗
DOE OSTI2022

Differentiable Preisach Modeling for Characterization and Optimization of Particle Accelerator Systems with Hysteresis

Future improvements in particle accelerator performance are predicated on increasingly accurate online modeling of accelerators. Hysteresis effects in magnetic, mechanical, and material components of accelerators are often neglected in online accelerator models used to inform control algorithms, even though reproducibility errors from systems exhibiting hysteresis are not negligible in high precision accelerators. Here, we combine the classical Preisach model of hysteresis with machine learning techniques to efficiently create nonparametric, high-fidelity models of arbitrary systems exhibiting hysteresis. We experimentally demonstrate how these methods can be used in situ, where a hysteresis model of an accelerator magnet is combined with a Bayesian statistical model of the beam response, allowing characterization of magnetic hysteresis solely from beam-based measurements. Finally, we explore how using these joint hysteresis-Bayesian statistical models allows us to overcome optimization performance limitations that arise when hysteresis effects are ignored.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Matrix-Free High-Performance Saddle-Point Solvers for High-Order Problems in \(\boldsymbol{H}(\operatorname{\textbf{div}})\)

Here, this work describes the development of matrix-free GPU-accelerated solvers for high-order finite element problems in H(div). The solvers are applicable to grad-div and Darcy problems in saddle-point formulation, and have applications in radiation diffusion and porous media flow problems, among others. Using the interpolation–histopolation basis, efficient matrix-free preconditioners can be constructed for the (1, 1)-block and Schur complement of the block system. With these approximations, block-preconditioned MINRES converges in a number of iterations that is independent of the mesh size and polynomial degree. The approximate Schur complement takes the form of an M-matrix graph Laplacian and therefore can be well-preconditioned by highly scalable algebraic multigrid methods. High-performance GPU-accelerated algorithms for all components of the solution algorithm are developed, discussed, and benchmarked. Numerical results are presented on a number of challenging test cases, including the “crooked pipe” grad-div problem, the SPE10 reservoir modeling benchmark problem, and a nonlinear radiation diffusion test case.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Mechanism of Peptide Agonist Binding in CXCR4 Chemokine Receptor

Chemokine receptors are key G-protein-coupled receptors (GPCRs) that control cell migration in immune system responses, development of cardiovascular and central nervous systems, and numerous diseases. In particular, the CXCR4 chemokine receptor promotes metastasis, tumor growth and angiogenesis in cancers. CXCR4 is also used as one of the two co-receptors for T-tropic HIV-1 entry into host cells. Therefore, CXCR4 serves as an important therapeutic target for treating cancers and HIV infection. Apart from the CXCL12 endogenous peptide agonist, previous studies suggested that the first 17 amino acids of CXCL12 are sufficient to activate CXCR4. Two 17-residue peptides with positions 1–4 mutated to RSVM and ASLW functioned as super and partial agonists of CXCR4, respectively. However, the mechanism of peptide agonist binding in CXCR4 remains unclear. Here, we have investigated this mechanism through all-atom simulations using a novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method. The Pep-GaMD simulations have allowed us to explore representative binding conformations of each peptide and identify critical low-energy states of CXCR4 activated by the super versus partial peptide agonists. Our simulations have provided important mechanistic insights into peptide agonist binding in CXCR4, which are expected to facilitate rational design of new peptide modulators of CXCR4 and other chemokine receptors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Harmonic analysis of discrete tracers of large-scale structure

It is commonplace in cosmology to analyze fields projected onto the celestial sphere, and in particular density fields that are defined by a set of points e.g. galaxies. When performing an harmonic-space analysis of such data (e.g. an angular power spectrum) using a pixelized map one has to deal with aliasing of small-scale power and pixel window functions. We compare and contrast the approaches to this problem taken in the cosmic microwave background and large-scale structure communities, and advocate for a direct approach that avoids pixelization. We describe a method for performing a pseudo-spectrum analysis of a galaxy data set and show that it can be implemented efficiently using well-known algorithms for special functions that are suited to acceleration by graphics processing units (GPUs). The method returns the same spectra as the more traditional map-based approach if in the latter the number of pixels is taken to be sufficiently large and the mask is well sampled. The method is readily generalizable to cross-spectra and higher-order functions. It also provides a convenient route for distributing the information in a galaxy catalog directly in harmonic space, as a complement to releasing the configuration-space positions and weights, and a route to spectral apodization. Finally, we make public a code enabling the application of our method to existing and upcoming datasets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Accelerated Aging of Fast Pyrolysis Bio-Oil Using Carbonyl Titration: Laboratory Analytical Procedure (LAP)

This laboratory analytical procedure covers the accelerated aging of fast pyrolysis bio-oils. Bio-oils undergo reactions that result in physical and chemical changes over time. These changes typically result in an increase in molecular weight, decrease in some functional groups such as carbonyls, and an increase in viscosity; additionally, the aging process often leads to phase separation. Studies have shown that accelerated aging of bio-oils using this method closely mimics room temperature aging for long periods of time (over 3 years). This procedure has been developed for the accelerated aging of fast pyrolysis bio-oils only. Fast pyrolysis bio-oils are more reactive than other bio-oils (e.g., catalytic fast pyrolysis) and undergo aging much more rapidly. Therefore, this procedure should only be used for fast pyrolysis bio-oil samples.

09 BIOMASS FUELS↗
DOE OSTI2022

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. Here, these methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗