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 703 records · Page 39

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

Summary of completed project

The primary goal of the work in this proposal has been the determination of the requirement for acceleration of particles, especially electrons, during solar flares. This is one of the important and yet unsolved problems in solar flares and other astronomical objects. This is done by detailed theoretical investigation of interaction of accelerated particles with ambient plasma including interaction with particles, waves, and large scale magnetic fields. Comparison with observed impulsive phase radiation such as X-rays, gamma-rays, microwaves, and type III radiation is used for constraining the models and for determination of the characteristics of the accelerated particles. Steady state and time-dependent numerical codes have been developed based on the Fokker-Planck method for study of the particle acceleration and transport during a flare. Using these, the correlation between characteristics of radiation at various wavelengths have been investigated and some constraints have been set on the model parameters. During the last three year period of this grant, dissertation research of three graduate students have been partially supported by this grant. Each of these theses has dealt with various observed properties of flare emission with particular emphasis on the interrelationship between radiation in at least two energy bands: X-rays and gamma-rays, X-rays and microwaves, and X-rays and type III radiation. This has allowed considerable tightening of constraints on models, and inclusion of wave particle interaction has allowed direct calculations of acceleration processes.

Petrosian, Vahe↗

Comparison of three explicit multigrid methods for the Euler and Navier-Stokes equations

Three explicit multigrid methods, Ni's method, Jameson's finite-volume method, and a finite-difference method based on Brandt's work, are described and compared for two model problems. All three methods use an explicit multistage Runge-Kutta scheme on the fine grid, and this scheme is also described. Convergence histories for inviscid flow over a bump in a channel for the fine-grid scheme alone show that convergence rate is proportional to Courant number and that implicit residual smoothing can significantly accelerate the scheme. Ni's method was slightly slower than the implicitly-smoothed scheme alone. Brandt's and Jameson's methods are shown to be equivalent in form but differ in their node versus cell-centered implementations. They are about 8.5 times faster than Ni's method in terms of CPU time. Results for an oblique shock/boundary layer interaction problem verify the accuracy of the finite-difference code. All methods slowed considerably on the stretched viscous grid but Brandt's method was still 2.1 times faster than Ni's method.

Chima, Rodrick V.↗

Comparison of three explicit multigrid methods for the Euler and Navier-Stokes equations

Three explicit multigrid methods, Ni's method, Jameson's finite-volume method, and a finite-difference method based on Brandt's work, are described and compared for two model problems. All three methods use an explicit multistage Runge-Kutta scheme on the fine grid, and this scheme is also described. Convergence histories for inviscid flow over a bump in a channel for the fine-grid scheme alone show that convergence rate is proportional to Courant number and that implicit residual smoothing can significantly accelerate the scheme. Ni's method was slightly slower than the implicitly-smoothed scheme alone. Brandt's and Jameson's methods are shown to be equivalent in form but differ in their node versus cell-centered implementations. They are about 8.5 times faster than Ni's method in terms of CPU time. Results for an oblique shock/boundary layer interaction problem verify the accuracy of the finite-difference code. All methods slowed considerably on the stretched viscous grid but Brandt's method was still 2.1 times faster than Ni's method.

Chima, Rodrick V.↗

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗

Nonparametric identification of a class of nonlinear close-coupled dynamic systems

A nonparametric identification technique for the identification of close coupled dynamic systems with arbitrary memoryless nonlinearities is presented. The method utilizes noisy recorded data (acceleration, velocity and displacement) to identify the restoring forces in the system. The masses in the system are assumed to be known (or fairly well estimated from the design drawings). The restoring forces are expanded in a series of orthogonal polnomials and the coefficients of these polynomial expansions are obtained by using least square fit method. A particularly simple and computationally efficient method is proposed for dealing with separable restoring forces. The identified results are found to be relatively insensitive to measurement noise. An analysis of the effects of measurement noise on the quality of the estimates is given. The computations are shown to be relatively quick (when compared say to the Wiener identification method) and the core storage required relatively small, making the method suitable for onboard identification of large space structures.

Udwadia, F. E.↗

Rolling Optimization of Transmission Network Recovery and Load Restoration Considering Hybrid Wind-Storage System and Cold Load Pickup

A common solution to deal with the stochasticity introduced by fast-ramping wind power integration is to equip wind farms (WFs) with energy storage systems (ESSs) to formulate hybrid WF-ESSs. In addition to leveling off wind power fluctuations during normal operations, a hybrid WF-ESS can be a flexible power source to accumulate the power system restoration. In this paper, we propose a rolling optimization model for transmission network recovery and load restoration considering the contributions of WF-ESSs. The proposed model is formulated as a mixed integer linear programming problem that simultaneously optimizes the amount and location of restorable load blocks as well as the restoration lines. The cold load pickup features of interrupted loads considering the outage duration are modeled in detail. A chance-constrained method is employed to deal with the uncertainty of wind power, and a rolling horizon-based framework is adopted to reduce the influence of forecast error. Case studies are conducted on both New England 39-bus system and part of a provincial power system in China. The results show that the load restoration process can be significantly accelerated by employing the proposed method and contributions of hybrid WF-ESSs to power system restoration are validated.

chance-constrained optimization↗

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗

Real World Applications of AI/ML in Optimizing Airspace Operations

As National Airspace System (NAS) is going through the Digital Transformation journey, data science and analytics methods can significantly contribute to improving the traditional physics-based decision-making tools. The adoption of AI/ML methods will not only help accelerate Federal Aviation Administration (FAA)’s vision of Info-centric NAS but also contribute to the overall objective of sustainable aviation. AI/ML can improve the ground and airspace operations by enhancing the accuracy of current decision-making tools used by the airlines and the FAA to manage traffic on the ground and in the air. Huge amount of data that gets collected during a flight. AI/ML methods can extract information from this data and provide valuable insights to make better operational decisions. NASA has partnered with the FAA and commercial airlines such as American and Southwest Airlines on this effort and has successfully demonstrated the benefits of using ML in real world environment by reducing delays and optimizing ground operations at the US airports. In 2022 itself, NASA demonstrated real-world benefits (over 24K lbs. of fuel savings, over 76.6K lbs. CO2 emission savings, and several hours of delay savings) by deploying ML based prediction models to optimize ground operations at Dallas/Fort Worth International and Dallas Love Field Airports in Texas. These tools are being deployed on the cloud for broader deployment, adaptability, and scalability. NASA is developed a reference implementation of the cloud-based platform to significantly lower the bar to development and distribution of these digital services for aviation. In this talk, I will share information about the Digital Information Platform project, the novel AI/ML based approaches used for optimizing ground operations and the opportunities to partner with NASA on these demonstrations.

air traffic management↗

Highly accelerated life testing (HALT): A review from a statistical perspective

Despite its use in one form or another for at least four decades, HALT and related techniques [e.g., highly accelerated-stress screening (HASS) and stress audits (HASA)] are not well understood within the statistical community and remain controversial. This largely reflects a conflict in motivation between engineers, testing under harsh conditions to discover and eliminate failure modes, and statisticians, taking a more cautious approach to develop quantitative estimates of parameters such as mean time between failures (MTBF). Here, this review article will clarify HALT concepts and methods and explain where it fits within the universe of methods that involve the application of accelerating factors to compress the time required to evaluate or enhance product reliability. A major distinction is between methods such as HALT, a high-stress test-analyze-fix-test iterative process directed at improving reliability by discovering and fixing weak points in a design, and quantitative accelerated life testing (QALT), whose goal is the estimation of product life for a fixed design. We discuss methods such as physics of failure that offer some hope of bridging the gap between the qualitative nature of HALT, and purely quantitative statistical methods. We present a variety of engineering applications of HALT including metal fatigue, piping and pressure vessels, structural damage, radiation damage, and rotating machinery. We also discuss potential synergies between HALT and QALT, such as rapid identification, through HALT, of failure modes requiring quantitative analysis. For further study, extensive references to the applicable literature are provided as well as an appendix that describes related methods.

97 MATHEMATICS AND COMPUTING↗

A multi-level solution algorithm for steady-state Markov chains

A new iterative algorithm, the multi-level algorithm, for the numerical solution of steady state Markov chains is presented. The method utilizes a set of recursively coarsened representations of the original system to achieve accelerated convergence. It is motivated by multigrid methods, which are widely used for fast solution of partial differential equations. Initial results of numerical experiments are reported, showing significant reductions in computation time, often an order of magnitude or more, relative to the Gauss-Seidel and optimal SOR algorithms for a variety of test problems. The multi-level method is compared and contrasted with the iterative aggregation-disaggregation algorithm of Takahashi.

Horton, Graham↗

Towards Agentic AI on Particle Accelerators

As particle accelerators grow in complexity, traditional control methods face increasing challenges in achieving optimal performance. This paper envisions a paradigm shift: a decentralized multi-agent framework for accelerator control, powered by Large Language Models (LLMs) and distributed among autonomous agents. We present a proposition of a self-improving decentralized system where intelligent agents handle high-level tasks and communication and each agent is specialized control individual accelerator components. This approach raises some questions: What are the future applications of AI in particle accelerators? How can we implement an autonomous complex system such as a particle accelerator where agents gradually improve through experience and human feedback? What are the implications of integrating a human-in-the-loop component for labeling operational data and providing expert guidance? We show two examples, where we demonstrate viability of such architecture.

43 PARTICLE ACCELERATORS↗

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. Here, in this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sensitivity of liquid bridges subject to axial residual acceleration

It has become evident that the float zone crystal growth method and isothermal liquid bridges may be very sensitive to the residual acceleration environment of a spacelab. Using a slender-body approximation, the problem of determining the axisymmetric response of the shape of the free surface of a cylindrical liquid column bounded by two solid regions is modeled by a one-dimensional system of nonlinear equations. It is found that the sensitivity of the zone shape depends on the static Bond number, BO, aspect ratio, and viscosity, as well as the amplitude and frequency of the disturbance. The general trend is an increase in tolerable residual gravity with increasing frequency. At the eigenfrequencies of the zone, however, there are dramatic deviations from this trend. At these frequencies the tolerable residual gravity level can be two orders of magnitude lower at this frequency.

Zhang, Yiqiang↗

Industrialized and Robotic Construction Advances in Terrestrial Construction and Opportunities in Space Construction

Advanced Building Construction (ABC) methods such as industrialized prefabrication and robotically controlled additive deposition are becoming promising technologies to support rapid construction of buildings at volume and scale. However, apart from the structural shell of the building, optimal integration of other systems – mechanical, electrical, plumbing, thermal, power generation, storage and distribution, is key to enabling functional buildings. Leveraging ABC methods to integrate such systems in the buildings requires the building design to be optimized for constructability constraints imposed by the industrialized and robotically controlled production setups entailing such ABC methods. The Industrialized Construction Innovation (ICI) team at the National Renewable Energy Laboratory (NREL) has been focusing on research and development of such ABC technologies and processes that can lead to production of buildings which are resilient to extreme climatic conditions at reduced cost and time. There is a potential synergy between such advancements in terrestrial construction, which can be leveraged for the benefit of space (orbital and extra-terrestrial) construction. Similarly, advances in space construction can lend itself to accelerated adoption of robotic construction methods in terrestrial construction. This presentation will give an overview of selective research efforts of the ICI team at NREL, along with a vision for potential collaborations with agencies like NASA.

3D printing↗

Radionuclide measurements by accelerator mass spectrometry at Arizona

Over the past years, Tandem Accelerator Mass Spectrometry (TAMS) has become established as an important method for radionuclide analysis. In the Arizona system the accelerator is operated at a thermal voltage of 1.8MV for C-14 analysis, and 1.6 to 2MV for Be-10. Samples are inserted into a cesium sputter ion source in solid form. Negative ions sputtered from the target are accelerated to about 25kV, and the injection magnet selects ions of a particular mass. Ions of the 3+ charge state, having an energy of about 9MeV are selected by an electrostatic deflector, surviving ions pass through two magnets, where only ions of the desired mass-energy product are selected. The final detector is a combination ionization chamber to measure energy loss (and hence, Z), and a silicon surface-barrier detector which measures residual energy. After counting the trace iosotope for a fixed time, the injected ions are switched to the major isotope used for normalization. These ions are deflected into a Faraday cup after the first high-energy magnet. Repeated measurements of the isotope ratio of both sample and standards results in a measurement of the concentration of the radionuclide. Recent improvements in sample preparation for C-14 make preparation of high-beam current graphite targets directly from CO2 feasible. Except for some measurements of standards and backgrounds for Be-10 measurements to date have been on C-14. Although most results have been in archaeology and quaternary geology, studies have been expanded to include cosmogenic C-14 in meteorites. The data obtained so far tend to confirm the antiquity of Antarctic meteorites from the Allan Hills site. Data on three samples of Yamato meteorites gave terrestrial ages of between about 3 and 22 thousand years.

Jull, A. J. T.↗

Three Component Velocity and Acceleration Measurement Using FLEET

The femtosecond laser electronic excitation and tagging (FLEET) method has been used to measure three components of velocity and acceleration for the first time. A jet of pure N2 issuing into atmospheric pressure air was probed by the FLEET system. The femtosecond laser was focused down to a point to create a small measurement volume in the flow. The long-lived lifetime of this fluorescence was used to measure the location of the tagged particles at different times. Simultaneous images of the flow were taken from two orthogonal views using a mirror assembly and a single intensified CCD camera, allowing two components of velocity to be measured in each view. These different velocity components were combined to determine three orthogonal velocity components. The differences between subsequent velocity components could be used to measure the acceleration. Velocity accuracy and precision were roughly estimated to be +/-4 m/s and +/-10 m/s respectively. These errors were small compared to the approx. 100 m/s velocity of the subsonic jet studied.

Danehy, Paul M.↗

Variable Depth Bragg Peak Method for Single Event Effects Testing

Traditionally, accelerator SEE testing is accomplished by removing the tops of packages so that the IC chips are accessible to heavy ions. However, ICs in some advanced packages cannot be de-lidded so a different approach is used that involves grinding and/or chemically etching away part of the package and the chip from the back side. The parts are then tested from the back side with ions having sufficient range to reach the sensitive volume. More recently, the entire silicon substrate in an SOI/SRAM was removed, making it possible to use low-energy ions with shorter ranges. Where removal of part of the package is not possible, facilities at Michigan State, NASA Space Radiation Laboratory, GANIL (France) and GSI (Germany) offer high-energy heavy ions with long ranges so that the ions can reach the devices' sensitive volumes without much change in the LET. Unfortunately, a run will typically involve only one ion species having a single energy and LET due to the long time it takes to tune a new energy. The Variable Depth Bragg Peak (VDBP) method is similar to the above method in that it involves the use of high-energy heavy ions that are able to pass through the packaging material and reach the device, obviating the need to remove the package. However, the method provides a broad range of LETs from a single ion by inserting degraders in the beam that modify the ion energy and, therefore, the LET. The crux of the method involves establishing a fiduciary point for degrader thickness, i.e., where the Bragg peak is located precisely at the sensitive volume in the device, for which the measured SEU cross-section and the ion LET are both also maxima and can be calculated using a Monte-Carlo program, TRIM. Once the fiduciary point has been established, calibrated high density polyethylene (HDPE) degraders are inserted into or removed from the beam to vary the ion LET at the device in a known manner. After each change of degrader thickness, the SEU cross-section is measured and the corresponding LET calculated from the change in degrader thickness. That information is used to generate a plot of cross-section as a function of ion LET. The advantages of this approach are that the part does not have to be de-lidded and a broad range of LETs is available from a single heavy ion without having to go to non-normal angles of incidence to change the "effective" LET. As we will show, it is possible to obtain an entire curve of cross-section versus LET using just two or three ions. Fig. 1 shows curves of cross-section vs LET for a Freescale 4 Mbit SOI/SRAM measured at the 88" Cyclotron at Berkeley and at NSRL. The open symbols are the data obtained from Berkeley for top-side and back-side irradiation. The solid data points are for the data obtained at NSRL using a device for which the package was intact. The data are for Iron and Gold and cover a range of LETs from 4 MeV.cm2/mg to 84 MeV.cm2/mg. The agreement between the data obtained from Berkeley and from NSRL is excellent, demonstrating that the VDBP method is capable of providing accurate values of cross-section versus LET, at least for the 4 Mbit SRAM. Details of the technique will be included in the final presentation.

Buchner, S.↗