Search NASA⌕ Search

SEARCH · Search NASA

Results for “Physical Algorithm”

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 739 records · Page 41

Algorithms for the Computation of Debris Risks

Determining the risks from space debris involve a number of statistical calculations. These calculations inevitably involve assumptions about geometry - including the physical geometry of orbits and the geometry of non-spherical satellites. A number of tools have been developed in NASA's Orbital Debris Program Office to handle these calculations; many of which have never been published before. These include algorithms that are used in NASA's Orbital Debris Engineering Model ORDEM 3.0, as well as other tools useful for computing orbital collision rates and ground casualty risks. This paper will present an introduction to these algorithms and the assumptions upon which they are based.

Matney, Mark↗

A 3-Channel Algorithm for Retrieving Spatially and Temporally Continuous Cloud Properties Across Different Geostationary Satellite Imagers

Cloud property retrieval algorithms for passive satellite imagers are generally designed to take advantage of all the useful spectral information available for a particular satellite. This strategy optimizes accuracy and reduces misidentification and retrieval biases, particularly for modern satellites with many spectral channels. However, the application of dissimilar algorithms tailored for different satellite sensor scan present a problem within the climate data record (CDR). Algorithm inconsistencies can introduce artificial trends in the CDR that are tied to instrument changes rather than physical changes, especially when older satellites with limited spectral information are included. The NASA CERES (Clouds and the Earth’s Radiant Energy System) data record provides global cloud property retrievals across 23 years and more than 25 satellites. With the goal of producing a spatially and temporally continuous record of cloud properties, the CERES cloud working group has developed algorithms that use only 3 channels that are common to most geostationary satellite imagers: 0.65, 3.9, and 10.8 μm.

Sarah Bedka↗

Direct structural retrieval from gas-phase ultrafast diffraction data using a genetic algorithm

Ultrafast scattering techniques such as ultrafast electron diffraction and ultrafast x-ray diffraction have been utilized to elucidate the structural dynamics, reaction intermediates, and final products in molecular reactions following photoexcitation. The time-dependent structures are typically not directly retrieved from the experimental data, but they rely on comparison with calculations. The genetic algorithm (GA), a global optimization strategy, can be used to retrieve the molecular structures directly from diffraction patterns without any theoretical input. However, the robustness of the GA with respect to real experimental conditions such as a limited momentum transfer range, noise, and artifacts has not been studied in detail. In this work, we characterize the performance of the GA with simulated data that mimic realistic experimental conditions. We have developed and implemented a variant of the GA specific to diffraction measurements which performs better in the presence of imperfect data compared to the standard implementation of the GA. We demonstrate this method with both synthetic data and experimental ultrafast electron diffraction data on the UV-induced photodissociation of trifluoroiodomethane (C⁢F 3⁡ I) molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

A physical split window technique for the retrieval of precipitable water from satellite measurements

The split-window precipitable water algorithm developed by Jedlovec (1987) is described, and its advantages are examined. In particular, it is noted that the algorithm can be applied to multispectral imaging data, available half-hourly from VAS; it uses only limited spatial averaging and can be applied to instruments that lack sounding channels. Here, previous work is reviewed, and the procedure is extended to the Advanced Very High Resolution Radiometer (AVHRR) and GOES-NEXT using simulated data. Root mean square errors for the simulated and actual retrievals are below +/- 5.8 mm. Of the instruments examined, AVHRR appears to have the greatest potential, with the rms errors below +/- 2.3 mm.

Guillory, Anthony R.↗

Link!: Potential Field Guidance Algorithm for In-Flight Linking of Multi-Rotor Aircraft

Link! is a multi-center NASA e ort to study the feasibility of multi-aircraft aerial docking systems. In these systems, a group of vehicles physically link to each other during flight to form a larger ensemble vehicle with increased aerodynamic performance and mission utility. This paper presents a potential field guidance algorithm for a group of multi-rotor vehicles to link to each other during flight. The linking is done in pairs. Each vehicle first selects a mate. Then the potential field is constructed with three rules: move towards the mate, avoid collisions with non-mates, and stay close to the rest of the group. Once a pair links, they are then considered to be a single vehicle. After each pair is linked, the process repeats until there is only one vehicle left. The paper contains simulation results for a system of 16 vehicles.

Cooper, John R.↗

ICARUS at the Short-Baseline Neutrino Program: First Results

The ICARUS collaboration employed the 760-ton T600 detector in a successful three-year physics run at the underground LNGS laboratory. In 2021, ICARUS started its new operation at Fermilab, collecting a substantial amount of neutrino events from the Booster Neutrino Beam (BNB) and the neutrinos at the Main Injector (NuMI) beam off-axis. These were used to test the ICARUS event selection, reconstruction, and analysis algorithms. ICARUS successfully completed its commissioning phase in June 2022, moving then to data taking for neutrino oscillation physics, aiming at first to either confirm or refute the claim by the Neutrino-4 short-baseline reactor experiment. ICARUS will also perform measurements of neutrino cross sections in LAr with the NuMI beam and several Beyond Standard Model studies. After the first year of operations, ICARUS will search for evidence of sterile neutrinos jointly with the Short-Baseline Near Detector, within the Short-Baseline Neutrino program. In this work, preliminary results from the ICARUS data with the BNB and NuMI beams are presented, both in terms of the performance of all ICARUS subsystems and the capability to select and reconstruct neutrino events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Physics of Acoustic Radiation from Jet Engine Inlets

Numerical simulations of acoustic radiation from a jet engine inlet are performed using advanced computational aeroacoustics (CAA) algorithms and high-quality numerical boundary treatments. As a model of modern commercial jet engine inlets, the inlet geometry of the NASA Source Diagnostic Test (SDT) is used. Fan noise consists of tones and broadband sound. This investigation considers the radiation of tones associated with upstream propagating duct modes. The primary objective is to identify the dominant physical processes that determine the directivity of the radiated sound. Two such processes have been identified. They are acoustic diffraction and refraction. Diffraction is the natural tendency for an acoustic wave to follow a curved solid surface as it propagates. Refraction is the turning of the direction of propagation of sound waves by mean flow gradients. Parametric studies on the changes in the directivity of radiated sound due to variations in forward flight Mach number and duct mode frequency, azimuthal mode number, and radial mode number are carried out. It is found there is a significant difference in directivity for the radiation of the same duct mode from an engine inlet when operating in static condition and in forward flight. It will be shown that the large change in directivity is the result of the combined effects of diffraction and refraction.

Tam, Christopher K. W.↗

Likelihood-Based Particle Identification in the Short-Baseline Near Detector

Accurate particle identification is crucial in any high-energy physics experiment, allowing scientists to understand the unique interactions and mechanisms at play in a detector. In this project, I develop and study a new particle identification (PID) algorithm for the Short-Baseline Near Detector, a likelihood-based approach, different from out current $\chi^2$ method. A likelihood estimation offers a more physically motivated strategy for PID. The distribution random energy losses of charged particles traveling through a medium are described by the Vavilov probability density function. By using this model, we can account for random energy losses and construct likelihood functions specific to each particle type, potentially enabling a more accurate method for PID.

Vanderwaal, Sophia [U. Alabama, Huntsville] (ORCID↗

Real-time plasma monitoring framework for advanced plasma control and ML-research in DIII-D

Real-time and adaptive plasma control is crucial for robust tokamak operation, requiring sensitivity and tolerance measurements of the plasma state. This paper presents the implementation of an integrated real-time plasma monitoring framework on the DIII-D tokamak to support advanced control approaches, including machine-learning (ML) methods. The system is built on the SHIELD framework, a high-performance modular architecture that provides a unified pipeline for integrating diverse diagnostics. The framework leverages high-bandwidth digitizers, fast numerical processing, and deterministic, low-latency interconnects to stream high-fidelity data from diagnostics such as electron cyclotron emission (ECE), beam emission spectroscopy (BES), CO interferometers, and a visible tangential divertor camera (TangTV). The system’s validity is demonstrated through direct comparisons of real-time and offline data. Furthermore, we present two key applications of the developed plasma monitoring system with ML-based plasma control strategies, including real-time divertor detachment and active Alfvén Eigenmode control. As a result, this work presents a robust and scalable approach for integrating high-frequency, multidimensional diagnostics into advanced control algorithms for future fusion devices.

AI/ML↗

Predictor-Corrector Entry Guidance for Reusable Launch Vehicles

An online entry guidance algorithm has been developed using a predictor-corrector approach. The algorithm is designed for the Reusable Launch Vehicle (RLV) and is demonstrated by using, the X-33 model. The objective of the design is to handle widely dispersed entry conditions and deliver the vehicle at the Terminal Area Energy Management (TAEM) interface box within an acceptable tolerance and without violating any of the vehicle physical constraints. Combination of several control variables is used in testing the performance and computational requirement of the algorithm. The control variables are the bank angle, angle-of-attack and the time for roll reversal. The bank angle and angle-of-attack profiles are the nominal profiles plus the perturbations in each direction. The initial guess of the bank profile is a 45 degrees bank angle with reversal at 360 seconds from liftoff. A six-element state vector is propagated to the TAEM interface box through the integration of the equations of motion (EOM). Altitude, heading and range errors are computed between the desired and the achieved state at the TAEM interface. These errors are used to correct the initial guess of the control variables. This process is repeated until the errors meet an acceptable level at the TAEM interface. Several numerical optimization methods are used to evaluate the convergent property of the predictor-predictor methodology. Successful results are demonstrated using the X-33 model.

Youssef, Hussein↗

Domain Knowledge Guided Bayesian Optimization For Autonomous Alignment Of Complex Scientific Instruments

Bayesian Optimization (BO) is a powerful tool for optimizing complex non-linear systems. However, its performance degrades in high-dimensional problems with tightly coupled parameters and highly asymmetric objective landscapes, where rewards are sparse. In such needle-in-a-haystack scenarios, even advanced methods like trust-region BO (TurBO) often lead to unsatisfactory results. We propose a domain knowledge guided Bayesian Optimization approach, which leverages physical insight to fundamentally simplify the search problem by transforming coordinates to decouple input features and align the active subspaces with the primary search axes. We demonstrate this approach's efficacy on a challenging 12-dimensional, 6-crystal Split-and-Delay optical system, where conventional approaches, including standard BO, TuRBO and multi-objective BO, consistently led to unsatisfactory results. When combined with an reverse annealing exploration strategy, this approach reliably converges to the global optimum. The coordinate transformation itself is the key to this success, significantly accelerating the search by aligning input co-ordinate axes with the problem's active subspaces. As increasingly complex scientific instruments, from large telescopes to new spectrometers at X-ray Free Electron Lasers are deployed, the demand for robust high-dimensional optimization grows. Our results demonstrate a generalizable paradigm: leveraging physical insight to transform high-dimensional, coupled optimization problems into simpler representations can enable rapid and robust automated tuning for consistent high performance while still retaining current optimization algorithms.

FOS: Computer and information sciences↗

From Latent Dynamics to Meaningful Representations

While representation learning has been central to the rise of machine learning and artificial intelligence, a key problem remains in making the learnt representations meaningful. For this the typical approach is to regularize the learned representation through prior probability distributions. However such priors are usually unavailable or are ad hoc. To deal with this, recent efforts have shifted towards leveraging the insights from physical principles to guide the learning process. In this spirit, we propose a purely dynamics-constrained representation learning framework. Instead of relying on predefined probabilities, we restrict the latent representation to follow overdamped Langevin dynamics with a learnable transition density — a prior driven by statistical mechanics. We show this is a more natural constraint for representation learning in stochastic dynamical systems, with the crucial ability to uniquely identify the ground truth representation. We validate our framework for different systems including a real-world fluorescent DNA movie dataset. Here, we show that our algorithm can uniquely identify orthogonal, isometric and meaningful latent representations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An adaptive and stability-promoting layerwise training approach for sparse deep neural network architecture

This work presents a two-stage adaptive framework for progressively developing deep neural network (DNN) architectures that generalize well for a given training data set. In the first stage, a layerwise training approach is adopted where a new layer is added each time and trained independently by freezing parameters in the previous layers. We impose desirable structures on the DNN by employing manifold regularization, sparsity regularization, and physics-informed terms. We introduce a ε – δ – stability-promoting concept as a desirable property for a learning algorithm and show that employing manifold regularization yields a ε – δ stability-promoting algorithm. Further, we also derive the necessary conditions for the trainability of a newly added layer and investigate the training saturation problem. In the second stage of the algorithm (post-processing), a sequence of shallow networks is employed to extract information from the residual produced in the first stage, thereby improving the prediction accuracy. Numerical investigations on prototype regression and classification problems demonstrate that the proposed approach can outperform fully connected DNNs of the same size. Moreover, by equipping the physics-informed neural network (PINN) with the proposed adaptive architecture strategy to solve partial differential equations, we numerically show that adaptive PINNs not only are superior to standard PINNs but also produce interpretable hidden layers with provable stability. As a result, we also apply our architecture design strategy to solve inverse problems governed by elliptic partial differential equations.

42 ENGINEERING↗

Using Markov Models of Fault Growth Physics and Environmental Stresses to Optimize Control Actions

A generalized Markov chain representation of fault dynamics is presented for the case that available modeling of fault growth physics and future environmental stresses can be represented by two independent stochastic process models. A contrived but representatively challenging example will be presented and analyzed, in which uncertainty in the modeling of fault growth physics is represented by a uniformly distributed dice throwing process, and a discrete random walk is used to represent uncertain modeling of future exogenous loading demands to be placed on the system. A finite horizon dynamic programming algorithm is used to solve for an optimal control policy over a finite time window for the case that stochastic models representing physics of failure and future environmental stresses are known, and the states of both stochastic processes are observable by implemented control routines. The fundamental limitations of optimization performed in the presence of uncertain modeling information are examined by comparing the outcomes obtained from simulations of an optimizing control policy with the outcomes that would be achievable if all modeling uncertainties were removed from the system.

Bole, Brian↗

Smart Pixel Sensors for the HL-LHC

Large-scale particle physics experiments produce tens of terabytes of data every second. Innovative methods to manage the data rate at the HL-LHC, which expects to operate at 10x the luminosity of what the LHC was initially designed for, are needed. AI-on the chip provides a way to intelligently filter out low momentum clusters in the pixel detector. This will open up an opportunity to use the pixel detector for the first time in the CMS Level-1 trigger, and lead to increased sensitivity to new physics measurements and searches. We have taped out our first chip, which incorporates a $p_T$ filtering algorithm on an ASIC chip. Our initial $p_T$ filtering algorithm considers clusters that are tracked by CMS. We will report on ongoing studies seeking to enhance the performance of our filter by utilizing unsupervised learning on untracked clusters, thus increasing background rejection.

43 PARTICLE ACCELERATORS↗

Mapping Forest Disturbance Intensity in North and South Carolina Using Annual Landsat Observations and Field Inventory Data

Disturbance and regrowth are vital processes in determining the roles of forest ecosystem in the carbon and biogeochemical cycles. Using time series observations, the vegetation change tracker (VCT) algorithm was designed to map the location, timing, and spectral magnitudes of forest disturbance events. While these spectral disturbance magnitudes are indicative of physical changes in tree cover or biomass, their quantitative relationships have yet to be established. This study focuses on estimating disturbance intensity as measured by percent basal area removal using spectral indices from the VCT algorithm over North and South Carolina. Repeat measurements on Forest Service Forest Inventory Analysis (FIA) ground plots, which provide changes in basal area between multiple dates at precise locations, are used for training and validation of the model. The overall R2 between predicted disturbance intensity and reference data is 0.66, and cross-validation prediction uncertainty is 14% in North Carolina. Possible causes of this uncertainty could be site heterogeneity and the temporal offset between ground measurements and satellite observations. Results show the area of stand clearing disturbances remains relatively stable around 1143 km2 yr−1 in North and South Carolina throughout the period of observations (1985–2015). The average amount of forest area affected by partial disturbance is much higher at 3287 km2 yr−1. The area of partial disturbances has strong inter-annual variability with a high value of 6000 km2 in 2007 and a low value of 1919 km2 in 2013.

Xin Tao↗

Attribution of heterogeneous stress distributions in low-grain polycrystals under conditions leading to damage

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain’s propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain’s mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure’s propensity to produce elevated stress states. As a result, these statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

36 MATERIALS SCIENCE↗

Pairing a Global Optimization Algorithm with EXAFS to Characterize Lanthanide Structure in Solution

Ensemble-average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. Here, to accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO 3 ) 3 dissolved in acetonitrile with a terpyridine ligand (terpyNO 2 ), and Nd(NO 3 ) 3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO 3 ) 3 (terpyNO 2 ) and Nd(NO 3 )(acetonitrile) 3 , with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗