Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local Method”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Global Aerodynamic Modeling Using Automated Local Model Networks in Real Time

A novel method is presented for automated real-time global aerodynamic modeling using local model networks, known as Smoothed Partitioning with Localized Trees in Real Time (SPLITR), as part of NASA’s Learn-to-Fly technology development initiative. The global nonlinear aerodynamics are partitioned into several local regions known as cells, with the dimension, location, and timing of each partition automatically selected based on a residual characterization procedure, under the constraints of real-time operation. Regression trees represent the successive partitioning of the global flight envelope and describe the evolution of the cell structure. Recursive equation-error least-squares parameter estimation in the time domain is used to estimate a model that represents the local aerodynamics in each region, so that it can be updated independently with non-contiguous data in the range of each cell over time. A weighted superposition of these piecewise local models across the flight envelope forms a global nonlinear model that also accurately captures the local aerodynamics. The SPLITR approach is demonstrated using both simulation and flight data, and the results are analyzed in terms of model predictive capabilities as well as interpretability. The results show that SPLITR can be used to automatically partition complex nonlinear aerodynamic behavior, produce an accurate model, and provide valuable physical insight into the local and global aerodynamics.

Weinstein, Rose↗
DOE OSTI2023

Data-Driven Smoothers for Extreme-Scale Computing

Patch-based relaxation refers to a family of methods for solving linear systems which partitions the matrix into smaller pieces often corresponding to groups of adjacent degrees of freedom residing within patches of the computational domain. The two most common families of patch-based methods are block-Jacobi and Schwarz methods, where the former typically corresponds to non-overlapping domains and the later implies some overlap. We focus on cases where each patch consists of the degrees of freedom on a finite element method mesh cell. Patch methods often capture complex local physics much more effectively than simpler point-smoothers such as Jacobi; however, forming, inverting, and applying each patch can be prohibitively expensive in terms of both storage and computation time. To this end, we propose several approaches for performing analysis on these patches and constructing a reduced representation. The compression techniques rely on either matrix norm comparisons or unsupervised learning via a clustering approach. We illustrate how it is frequently possible to retain/factor less than 5% of all patches and still develop a method that converges only a little slower than when all patches are stored/factored.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Resolution of the 1D regularized Burgers equation using a spatial wavelet approximation

The Burgers equation with a small viscosity term, initial and periodic boundary conditions is resolved using a spatial approximation constructed from an orthonormal basis of wavelets. The algorithm is directly derived from the notions of multiresolution analysis and tree algorithms. Before the numerical algorithm is described these notions are first recalled. The method uses extensively the localization properties of the wavelets in the physical and Fourier spaces. Moreover, the authors take advantage of the fact that the involved linear operators have constant coefficients. Finally, the algorithm can be considered as a time marching version of the tree algorithm. The most important point is that an adaptive version of the algorithm exists: it allows one to reduce in a significant way the number of degrees of freedom required for a good computation of the solution. Numerical results and description of the different elements of the algorithm are provided in combination with different mathematical comments on the method and some comparison with more classical numerical algorithms.

Liandrat, J.↗
NASA NTRS2006

Ground-Based Localization of Mars Rovers

The document discusses a procedure for localizing the Mars rovers in site frame, a locally defined reference frame on the Martian surface. MER onboard position within a site frame is estimated onboard and is based on wheel odometry. Odometry estimation of rover position is only reliable over relatively short distances assuming no wheel slip, sinkage, etc. As the rover traverses, its onboard estimate of position in the current site frame accumulates errors and will need to be corrected on occasions via relocalization on the ground (mission operations). The procedure provides a systematic process for ground operators to localize the rover. The method focuses on analysis of acquired images used to declare a site frame and images acquired post-drive. Target selection is performed using two main steps. In the first step, the user identifies features of interest from the images used to declare the current site. Each of the selected target s position in site frame is recorded. In the second step, post-traverse measurements of the selected features positions are recorded again, this time in rover frame, using images acquired post-traverse. In the third step, we transform the post-traverse target s positions to local level frame. In the fourth step, we compute the delta differences in the pre- and post-traverse target s position. In the fifth step, we analyze the delta differences with techniques that compute their statistics to determine the rover s position in the site frame.

Trebi-Ollennu, Ashitey↗
DOE OSTI2025

Weakly supervised anomaly detection for resonant new physics in the dijet final state using proton-proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

An anomaly detection search for narrow-width resonances beyond the Standard Model that decay into a pair of jets is presented. The search is based on 139 fb −1 of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded during 2015–2018 with the ATLAS detector at the Large Hadron Collider. The analysis is optimized without a particular signal model and aims to be sensitive to a broad range of new physics. It uses two different machine learning strategies to estimate the background in different signal regions. In each region, a weakly supervised classifier is trained to distinguish this background model from data. The analysis focuses on events with high transverse momentum jets reconstructed as large-radius jets. The mass and substructure of these jets are used as inputs to the classifiers. After a classifier-based selection, the distribution of the invariant mass of the two jets is used to search for potential local excesses. The model-independent results of both the anomaly detection methods show no signs of significant local excesses. In addition to model-independent results, a representative set of signal models is injected into the data, and the sensitivity of the methods to these scenarios is reported.

Aad, G. [Aix-Marseille Université] (ORCID:00000002↗
NASA NTRSDate not supplied

Methods for System-Level Multidisciplinary Uncertainty Analysis of Low-Boom Flight Vehicles

Current research supporting NASA’s Commercial Supersonic Technology project is focused on the efficient prediction of uncertainty in sonic boom loudness generated by low-boom aircraft concepts. This paper focuses on research incorporating aircraft trim and aerostructural analysis into a multidisciplinary system-level uncertainty analysis. This enables the modeling of a steady-state representation of a point in the uncertainty space, simulating the vehicle as it would be flown. This approach also enables multiple uncertain parameters defining the configuration of the vehicle to be reduced to three: Mach number, altitude, and aircraft weight. To demonstrate this methodology, a case study exploring a conceptual low-boom supersonic aircraft is performed. Two different approaches are used to model the interactions between nearfield pressure signature analysis and sonic boom propagation, and their performance is evaluated in terms of accuracy and computational expense. One method uses a set of local surrogate models to generate a large number of nearfield signatures and perform Monte Carlo analysis. This method is found to produce, at a lower expense, uncertainty metrics that are comparable to the second method, in which uncertainty metrics are computed based on loudness metric values obtained directly from simulated nearfield signatures.

UQ↗
NASA NTRSDate not supplied

Method for System-Level Multidisciplinary Uncertainty Analysis of Low-Boom Flight Vehicles

Current research supporting NASA’s Commercial Supersonic Technology project is focused on the efficient prediction of uncertainty in sonic boom loudness generated by low-boom aircraft concepts. This paper focuses on research incorporating aircraft trim and aerostructural analysis into a multidisciplinary system-level uncertainty analysis. This enables the modeling of a steady-state representation of a point in the uncertainty space, simulating the vehicle as it would be flown. This approach also enables multiple uncertain parameters defining the configuration of the vehicle to be reduced to three: Mach number, altitude, and aircraft weight. To demonstrate this methodology, a case study exploring a conceptual low-boom supersonic aircraft is performed. Two different approaches are used to model the interactions between nearfield pressure signature analysis and sonic boom propagation, and their performance is evaluated in terms of accuracy and computational expense. One method uses a set of local surrogate models to generate a large number of nearfield signatures and perform Monte Carlo analysis. This method is found to produce, at a lower expense, uncertainty metrics that are comparable to the second method, in which uncertainty metrics are computed based on loudness metric values obtained directly from simulated nearfield signatures.

Supersonics↗
NASA NTRS2012

Implementation and Testing of the First Prompt Search for Electromagnetic Counterparts to Gravitational Wave Transients

A transient astrophysical event observed in both gravitational wave (OW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient OW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize OW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 172009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify OW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency OW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with approx.50% or better probability with a few pointings of wide-field telescopes. Key words. gravitational waves - methods: observational

Camp, Jordan↗
NASA NTRS1992

Constraints on baryonic dark matter in the Galactic halo and Local Group

A four-color method and deep CCD data are used to search for very faint metal-poor stars in the direction of the south Galactic pole. The results make it possible to limit the contribution of ordinary old, metal-poor stars to the dynamical halo of the Galaxy or to the Local Group. The ratio of the mass of the halo to its ordinary starlight must be more than about 2000, unless the halo is very small. For the Local Group, this ratio is greater than about 400. If this local dark matter is baryonic, the process of compact-object formation must produce very few 'impurities' in the form of stars similar to those found in globular clusters. The expected number of unbound stars with MV not greater than 6 within 100 pc of the sun is less than 1 based on the present 90-percent upper limit to the Local Group starlight.

Richstone, Douglas↗
DOE OSTI2020

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Topological Data Analysis for Particulate Gels

Soft gels, formed via the self-assembly of particulate materials, exhibit intricate multiscale structures that provide them with flexibility and resilience when subjected to external stresses. Here, this work combines particle simulations and topological data analysis (TDA) to characterize the complex multiscale structure of soft gels. Our TDA analysis focuses on the use of the Euler characteristic, which is an interpretable and computationally scalable topological descriptor that is combined with filtration operations to obtain information on the geometric (local) and topological (global) structure of soft gels. We reduce the topological information obtained with TDA using principal component analysis (PCA) and show that this provides an informative low-dimensional representation of the gel structure. We use the proposed computational framework to investigate the influence of gel preparation (e.g., quench rate, volume fraction) on soft gel structure and to explore dynamic deformations that emerge under oscillatory shear in various response regimes (linear, nonlinear, and flow). Our analysis provides evidence of the existence of hierarchical structures in soft gels, which are not easily identifiable otherwise. Moreover, our analysis reveals direct correlations between topological changes of the gel structure under deformation and mechanical phenomena distinctive of gel materials, such as stiffening and yielding. In summary, we show that TDA facilitates the mathematical representation, quantification, and analysis of soft gel structures, extending traditional network analysis methods to capture both local and global organization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Instability of non-zonal baroclinic flows - Multiple-scale analysis

The linear instability of a non-zonal flow can be reduced to an eigenvalue-eigenfunction problem, governed by a nonseparable partial differential equation (Niehaus, 1980). Approximate solutions, found by the method of multiple scales, are derived here and compared with earlier results found using a spectral method. The amplitude maxima are correctly located. The zonal variations of local wavenumber and of amplitude are qualitatively correct, but not sufficiently extreme. Because the method is oversensitive to local conditions, and less sensitive to global constraints, this comparison provides theoretical limits to the possibility of parameterizing transient eddies in terms of the local time mean state of the atmosphere. The method can be extended easily to flows with more realistic vertical structure.

Niehaus, M. C. W.↗
DOE OSTI2020

Application of Self-Interaction Corrected Density Functional Theory to Early, Middle, and Late Transition States

Density functional theory (DFT)-based methods often significantly underpredict chemical reaction barriers compared with experiments because of the tendency of DFT to over-stabilize transition states with stretched bonds due to the impact of unphysical electron self-interaction. But many reactions have early or late transition states where the transition state geometry closely resembles the reactants or products, respectively. The role of self-interaction in those cases is not known. Here we compare the performance of DFT with and without self-interaction correction (SIC) for describing the hydrogenation of CO and CO 2 catalyzed by a Lewis acid-base pair incorporated onto an aromatic cluster, using CCSD(T) results for reference. The three elementary steps in these reactions consist of an early, a middle, and a late transition. Our results show that the Perdew-Zunger SIC (PZ-SIC), implemented in the Fermi-L Lowdin orbital SIC (FLO-SIC) approach, qualitatively improves the description of the forward and reverse reaction barriers relative to uncorrected DFT for the middle transition, but not the early or late transitions. By contrast, the local scaling SIC (LSIC) method, also implemented in the FLO-SIC framework, significantly improves the calculated barriers over DFT and PZ-SIC in all but one case. Lastly, the results also show how the FLO-SIC approach can provide insight into the bonding in aromatic systems.

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

Heterogeneity in MacMullin Number of Li-Ion Battery Electrodes Studied by Means of an Aperture Probe

Heterogeneity of MacMullin number within battery electrodes is a key metric affecting cell performance. To characterize this heterogeneity, an aperture probe was developed. This probe, coupled with a newly developed transmission-line model, allows for measurements of tortuosity, represented by the MacMullin number, on millimeter length scales. Local MacMullin number values of seven electrodes were measured, and the ionic resistance profiles of these electrodes are given through contour maps of the MacMullin number. The method is validated by comparing the average MacMullin number to the value obtained through other measurement methods. The results show significant local MacMullin number variation in such electrodes on a millimeter length scale. Furthermore, this method will allow battery manufacturers and researchers to better quantify sources of heterogeneity and improve electrode quality.

25 ENERGY STORAGE↗
NASA NTRS1982

Marching grid generation using parabolic partial differential equations

The feasibility of using parabolic partial differential equations for grid generation is examined. Source terms in the form of a linear interpolation between the current grid and the outer boundary are used in generating a grid for a two-dimensional airfoil flow. Grid generation equations are derived by assuming that grid spacings are locally nonuniform on the computational domain. The grid spacing control method is described in detail, and the local orthogonality of the grid lines is discussed. An O-mesh and an H-mesh generated by the described method are shown.

Nakamura, S.↗
DOE OSTI2021

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗