Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational 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.

618 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Variational Quantum Circuits to Prepare Low Energy Symmetry States

We explore how to build quantum circuits that compute the lowest energy state corresponding to a given Hamiltonian within a symmetry subspace by explicitly encoding it into the circuit. We create an explicit unitary and a variationally trained unitary that maps any vector output by ansatz A(α → ) from a defined subspace to a vector in the symmetry space. The parameters are trained varitionally to minimize the energy, thus keeping the output within the labelled symmetry value. The method was tested for a spin XXZ Hamiltonian using rotation and reflection symmetry and H 2 Hamiltonian within S z = 0 subspace using S 2 symmetry. We have found the variationally trained unitary gives good results with very low depth circuits and can thus be used to prepare symmetry states within near term quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

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

Robust Spectral Anomaly Detection in EELS Spectral Images via 3D Convolutional Variational Autoencoders

Abstract A 3D Convolutional Variational Autoencoder (3D‐CVAE) is introduced for automated anomaly detection in electron energy‐loss spectroscopy spectrum imaging (EELS‐SI) data. This approach leverages the full 3D structure of EELS‐SI data to detect subtle spectral anomalies while preserving both spatial and spectral correlations across the datacube. By employing cross‐entropy loss and training on bulk spectra, the model learns to reconstruct bulk features characteristic of the defect‐free material. In exploring methods for anomaly detection, both the 3D‐CVAE approach and principal component analysis (PCA) are evaluated, testing their performance using FeL‐edge ΔEpeak shifts designed to simulate material defects. These results show that 3D‐CVAE achieves superior anomaly detection and maintains consistent performance across various shift magnitudes. The method demonstrates clear bimodal separation between bulk and anomalous spectra, enabling reliable classification. Further analysis verifies that lower‐dimensional representations are robust to anomalies in the data. While performance advantages over PCA diminish with decreasing anomaly concentration, our method maintains high reconstruction quality even in challenging, noise‐dominated spectral regions. This approach provides a robust framework for unsupervised automated detection of spectral anomalies in EELS‐SI data, particularly valuable for analyzing complex material systems.

Chemistry↗
NASA NTRS1965

Computation of Optimal Interplanetary Low-Thrust Trajectories With Bounded Thrust Magnitude By Means of the Generalized Newton-Raphson Method

The generalized Newton-Raphson method, an iterative procedure for solving nonlinear operator equations, has been extended in application to variational problems with bounded control variables. A minimum fuel interplanetary low thrust orbital transfer problem is worked out in detail to demonstrate the practical aspects of the algorithm as well as its computational effectiveness. The control variables are the thrust magnitude, limited from zero to some prescribed maximum value, and the thrust steering angle.

Computation↗
DOE OSTI2024

Uncovering obscured phonon dynamics from powder inelastic neutron scattering using machine learning

The study of phonon dynamics is pivotal for understanding material properties, yet it faces challenges due to the irreversible information loss inherent in powder inelastic neutron scattering spectra and the limitations of traditional analysis methods. In this study, we present a machine learning framework designed to reveal obscured phonon dynamics from powder spectra. Using a variational autoencoder, we obtain a disentangled latent representation of spectra and successfully extract force constants for reconstructing phonon dispersions. Notably, our model demonstrates effective applicability to experimental data even when trained exclusively on physics-based simulations. The fine-tuning with experimental spectra further mitigates issues arising from domain shift. Analysis of latent space underscores the model’s versatility and generalizability, affirming its suitability for complex system applications. Furthermore, our framework’s two-stage design is promising for developing a universal pre-trained feature extractor. This approach has the potential to revolutionize neutron measurements of phonon dynamics, offering researchers a potent tool to decipher intricate spectra and gain valuable insights into the intrinsic physics of materials.

domain adaptation↗
NASA NTRS2026

Characterizing the subglacial environment of lower Thwaites Glacier using radar modeling

Understanding the spatial heterogeneity beneath Thwaites Glacier, West Antarctica, is vital to projecting its impact on future sea levels. Radar-echo sounding (RES) is commonly used to infer subglacial conditions, but these data can be challenging to interpret. We assess basal heterogeneity across Thwaites Glacier by comparing RES returns to a radar backscattering simulator for over 400 km of RES data. The modeled variations in bed returned power exhibited a strong correlation with actual RES data in 40% of our simulated flight segments, which we consider evidence for a relatively homogeneous glacier bed. Other sites (40%) demonstrated improved fit quality when hydrology or substrate transitions were introduced in the bed material model. The remaining simulated segments (20%) were diagnosed as having more complex basal heterogeneity. The spatial distribution of complex heterogeneity appears to coincide with asymmetric patterns in the RES specularity content, which has been interpreted in previous studies as a signature for channelized hydrology. Conversely, the homogeneous substrate locations coincide with areas of fast-moving ice in western Thwaites. Our simulation method can isolate power variations induced by material heterogeneity vs topography, which is an important limitation of existing RES analysis methods.

Antarctic glaciology↗
DOE OSTI2025

Environmental drivers of spatial variation in tropical forest canopy height: Insights from NASA’s GEDI spaceborne LiDAR

Forest canopy height is a fundamental ecosystem property—influencing patterns of forest carbon storage and forest ecosystem responses to climate variability and change. Previous studies have analyzed environmental drivers influencing spatial variation in canopy height at landscape-to-regional scales; however, far less is known about the environmental determinants underlying regional and global scale variation in forest canopy height. Using the canopy height metrics products from Global Ecosystem Dynamics Investigation (GEDI), a space-borne Light Detection and Ranging (LiDAR) instrument specifically designed to characterize forest structure, we analyze the environmental correlates of spatial variation of global tropical forest canopy height. Our study demonstrates that climate, topography, and soil properties account for 75% of the variation in tropical forest canopy height. Elevation, dry season length, and solar radiation are the most important drivers in determining canopy height both locally and regionally. These results emphasize the vulnerability of tropical forest structure to ongoing changes in the earth’s climate and provide a valuable empirical baseline for tropical forest management.

Science & Technology - Other Topics↗
NASA NTRS1993

Use of Global Positioning System Measurements to Determine Geocentric Coordinates and Variations in Earth Orientation

Geocentric tracking station coordinates and short-period Earth-orientation variations can be measured with Global Positioning System (GPS) measurements. Unless calibrated, geocentric coordinate errors and changes in Earth orientation can lead to significant deep-space tracking errors. Ground-based GPS estimates of daily and subdaily changes in Earth orientation presently show centimeter-level precision. Comparison between GPS-estimated Earth-rotation variations, which are the differences between Universal Time 1 and Universal Coordinated Time (UT1-UTC), and those calculated from ocean tide models suggests that observed subdaily variations in Earth rotation are dominated by oceanic tidal effects. Preliminary GPS estimates for the geocenter location (from a 3-week experiment) agree with independent satellite laser-ranging estimates to better than 10 cm. Covariance analysis predicts that temporal resolution of GPS estimates for Earth orientation and geocenter improves significantly when data collected from low Earth-orbiting satellites as well as from ground sites are combined. The low Earth GPS tracking data enhance the accuracy and resolution for measuring high-frequency global geodynamical signals over time scales of less than 1 day.

R P Malla↗
DOE OSTI2026

SQUID Readout of a High-$Q$ Superconducting $LC$ Resonator

We demonstrate the readout of a superconducting $LC$ resonator, with $Q \approx 2 \times 10^{6}$ at approximately 250 kHz, using a dc SQUID followed by a SQUID series-array amplifier readout chain. We find that $Q$ depends on the SQUID flux-bias point, increasing near the shallow-slope point of the SQUID modulation curve and decreasing near the steep-slope point, consistent with the previously observed SQUID damping effects. From this variation, we infer the SQUID effective input impedance. We further infer the effective temperature of the resonator circuit from the resonance peak in the noise spectrum and show that it also depends on the SQUID flux-bias point, suggesting a possible contribution from SQUID back-action noise. This system provides a working prototype for future experiments based on lumped-element resonators with SQUID readout, in particular low-mass axion searches within the DMRadio program.

FOS: Physical sciences↗
NASA NTRSDate not supplied

Impact of Technology and Mission Variations on Aircraft Designed for Advanced Air Mobiity

NASA is conducting investigations in Advanced Air Mobility (AAM) aircraft and operations, including the development of Urban Air Mobility (UAM) aircraft designs that can be used to focus and guide research activities in support of AAM. This report is an investigation of the impact of technology and mission variations on several of the NASA AAM concept aircraft: quadrotor, quiet single main rotor, side-by-side, and tiltrotor configurations, with turboshaft and electric propulsion variants for each. First, the mission and aircraft models of the baseline designs were reassessed and updated, including rotor geometry optimization, update of the rotor performance models, and disk loading optimization. For these eight designs, technology and mission excursions were performed. Relative to the calibration cases that can be considered examples of good design practice, the impact of the weight technology factors is significant. For the electric aircraft, there is a very large impact of battery specific energy (Wh/kg), and correspondingly a very large impact of mission range. The vision of Advanced Air Mobility is driven by missions that will enable new transportation capabilities. Hence it is appropriate to compare Concept Vehicles of different lift and propulsive architectures, all designed to accomplish the same UAM mission. It is also useful however to consider specific missions that can take advantage of the strengths of individual aircraft configurations. So alternate designs were also developed for the concept vehicles: for turboshaft aircraft, longer unrefueled range, including faster cruise speed for the tiltrotor; for electric aircraft, shorter range and more realistic battery weight.

Technology↗
DOE OSTI2026

Fast Active-Set Thresholding Method for Nonnegative Least Squares

Nonnegative Least Squares (NNLS) is a fundamental constrained optimization problem encountered in many applications such as image deblurring, signal processing, nonnegative matrix factorization, magnetic microscopy, and hyperspectral imaging. Active-set based methods are a common class of algorithms for solving NNLS which identify the optimal variable set of the NNLS solution. They do so by iteratively solving a series of unconstrained least squares problems, identifying which variables violate the nonnegativity constraints, and then swapping variables in/out of consideration until the optimal set of variables is found. Several variations improving upon this method exist in the literature. In this work, we propose an active-set swap heuristic which further improves upon existing active-set based methods for NNLS. Our optimizations are based upon adding multiple variables to the passive set within a threshold of the smallest gradient value and removing variables within a similar threshold of the closest boundary constraint. We leverage these optimizations to yield a Fast Active-Set Thresholding NNLS (FAST-NNLS) algorithm which significantly outperforms the existing state-of-the-art NNLS algorithms for a wide range of problems. Rigorous convergence guarantees are proven for the proposed method. We demonstrate the effectiveness of our proposed method on multiple synthetic datasets and two realworld text analysis applications. In doing so, we present the most comprehensive NNLS solver comparison in the literature to date.

Cobb, Benjamin [Georgia Institute of Technology]↗
DOE OSTI2025

Probing the role of local tunnel variations in early-stage lithiation of α-MnO₂ nanowires via in situ TEM

Understanding lithium-ion transport in tunnel-structured manganese oxides is essential for designing high-performance lithium-ion battery electrode materials. Here, we elucidate the early-stage lithiation mechanism of potassium-stabilized α-MnO 2 nanowires using in situ transmission electron microscopy (TEM) coupled with electron energy-loss spectroscopy (EELS), high-resolution TEM (HRTEM), and geometric phase analysis (GPA). Real-time TEM imaging reveals clear volume expansion at the reaction front, while EELS analysis uncovers lithium-ion diffusion far beyond this region, where no visible expansion is observed, indicating fast, defect-assisted transport. GPA and HRTEM analyses show that localized tensile and compressive strain fields, originating from pre-existing local tunnel structural variations, persist after lithiation. The tensile-strained regions enable lithium-ion insertion with minimal lattice distortion, offering additional free volume that facilitates rapid lithium-ion accommodation ahead of the structural transformation. Our results demonstrate a local tunnel variation-mediated fast diffusion pathway that precedes bulk reaction, underscoring the critical role of local strain in enabling early-stage lithium transport. Given the structural versatility of MnO 2 and its ability to accommodate diverse atomic arrangements beyond the well-known tunnel phases (β-, γ-, δ-, λ-, R-phases), our findings highlight the importance of understanding and engineering local structural environments. This work provides fundamental insights into the interplay between defects, strain, and ion dynamics, and presents defect engineering as a promising approach to enhance both rate performance and structural stability in manganese-based cathodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Kalman Filtering and RTS Smoothing for Arc-Jet Sample Edge Tracking

Accurate and temporally consistent measurements of test-article recession are required to characterize material response during arc-jet ablation experiments. However, image-based boundary measurements are often affected by segmentation noise, brightness variations, and frame-to-frame variability. This work extends arcjetCV [1] with filtering methods for tracking the evolving boundaries of hemispherical and ISO-Q test articles. Two boundary-tracking approaches based on Kalman filtering [2] were implemented. The first applies a point-wise Kalman filter followed by a Rauch–Tung–Striebel smoother [3] to individual boundary-point locations. The second applies the same filtering and smoothing framework to a reduced set of geometry-dependent shape parameters. The point-wise method reduces local frame-to-frame fluctuations while preserving spatial details along the detected boundary. The shape-parameter method provides a compact and geometrically constrained estimate of the sample contour. Figure 1 illustrates the two approaches for a hemispherical test article. Both methods improve temporal consistency and support more robust estimation of surface recession from arc-jet imagery. The two filtering approaches provide complementary representations of boundary evolution and improve the reliability of image-based recession measurements during arc-jet experiments.

recession measurement↗
NASA NTRSDate not supplied

Effect of Composition Variation on Isothermal Rapid Curing Resins for Aerospace Applications

To meet the projected demand for single-aisle composite aircraft in 2040, the production rate of these aircraft will need to increase by four to six times what is currently achievable. A composite manufacturing process with the potential to enable the targeted rates is isothermal resin infusion where infusion and cure occur at the same temperature. However, presently there is not an accepted aerospace-grade resin that can support this isothermal infusion process and maintain acceptable resin properties. Therefore, resin development initiatives have designed anionic catalyzed epoxy systems capable of isothermal infusion and rapid cure at a temperature below 100 ˚C in one hour. High-performance aerospace properties are achieved after undergoing a freestanding post-cure. The resins facilitate faster, energy efficient production cycles by reducing the time and temperature required for cure and eliminating temperature ramps typically required of materials used in the aerospace industry. This work discusses several of the developed resin systems and the impact of varying composition on rheological, thermal, and mechanical properties. Carbon fiber composites were fabricated for one composition by resin transfer molding to demonstrate the potential of these isothermal, rapid cure systems for meeting industry needs.

resin transfer molding↗
DOE OSTI2024

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗
DOE OSTI2026

Seasonal variation of the surface cross-shelf exchange in the northern South China Sea: a Lagrangian perspective

Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.

Hao, Dongliang↗
NASA NTRS2026

Ground-Based Microwave Measurements of ClO and O3 Over Scott Base, Antarctica

We present measurements from a new microwave instrument that has been deployed at Scott Base, Antarctica, primarily to measure ClO. The new instrument allows for the extension of the time series of lower stratospheric ClO measurements begun in 1996 at this site, while some of the technical improvements offer the possibility of making additional measurements. We show measurements of ClO and O 3 in both the upper and lower stratosphere near both the fall and spring equinox and document the diurnal variations of these measurements throughout the season. The ground-based measurements are compared with Aura Microwave Limb Sounder (MLS) measurements, which are available twice per day, and with a NASA Goddard 2D model (GSFC2D) simulation. We show that, in agreement with the GSFC2D model diurnal variations, upper stratospheric ClO mixing ratios gradually increase throughout the daylight hours, while in the lower stratospheric austral spring, when heterogeneous chemistry associated with polar stratospheric clouds occurs, daily maximum ClO mixing ratios of near 2 ppbv occur near noon. We also show measurements of stratospheric O 3 , primarily in the upper stratosphere, which, while exhibiting a low bias of ∼15%–20% relative to Aura MLS, show temporal variations in good agreement with MLS.

stratospheric ozone↗
NASA NTRS1982

Dynamic environments for space shuttle payloads

Payload bay dynamic data from the first two space shuttle flights are summarized and evaluated. Development of dynamic environment design and test criteria for shuttle payloads from measured flight data is discussed. Factors that must be considered are flight to flight variations, spatial variations, temporal variations, measurement bias errors and the degree of confidence desired that a predicted environment will not be exceeded in flight. Summary and conclusion reports will be published after STS-4 and at appropriate intervals thereafter. The nature of these future reports and their impact on the user community is discussed.

Dennis L Kern↗