Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Analytic partial derivatives for estimating low-thrust parameters.

Analytic partial derivatives for estimating orbital low-thrust parameters via differential correction are developed and compared with two different numerical methods. The formulation is independent of the particular thrust model used and is applicable to all physically possible elliptic orbits. The starting point for the development is the set of variational equations of the elliptic orbital elements in the form due to Lagrange. The first time derivatives of the elements are transformed to derivatives with respect to the space variable, true anomaly, and integrated to first order in closed form in a straightforward general perturbations approach, with one exception: particular attention is given to the mean anomaly as influenced by thrust perturbations in the semimajor axis so that the complete first-order effect is included. The partials of the elements are then taken with respect to any given thrust parameter. Two comparisons are made with numerical methods for computing these thrust partials: numerical quotients and numerical integration of the variational equations for thrust.

Cunningham, G. W.↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
NASA NTRS1981

Investigation of control, display, and crew-loading requirements for helicopter instrument approach

A ground simulation experiment was conducted on a flight simulator for advanced aircraft to investigate the influence and interaction of flight-control system, flight-director display, and crew-loading situation on helicopter flying qualities during terminal-area operations in instrument conditions. Six levels of control complexity were implemented on a representative helicopter model. The six levels of augmentation were examined with display variations consisting of raw elevation and azimuth data only and of raw data plus one-, two-, and three-cue flight directors. Crew-loading situations simulated for the control-display combinations were dual-pilot operation and single-pilot operation. Four pilots performed a total of 150 evaluations of combinations of these parameters for a representative microwave landing system approach task. Pilot rating results indicated the existence of a control display trade-off for ratings of satisfactory, whereas ratings of adequate-but-unsatisfactory depended primarily on the control system; the control system required for ratings of adequate-but-unsatisfactory was clearly more complex for the single-pilot situation than that for the dual-pilot situation.

Lebacqz, J. V.↗
DOE OSTI2025

Application of Lagrangian techniques for calculating the on-axis rotational transform

The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6$ X $6$ matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Finding simplicity: unsupervised discovery of features, patterns, and order parameters via shift-invariant variational autoencoders *

Abstract Recent advances in scanning tunneling and transmission electron microscopies (STM and STEM) have allowed routine generation of large volumes of imaging data containing information on the structure and functionality of materials. The experimental data sets contain signatures of long-range phenomena such as physical order parameter fields, polarization, and strain gradients in STEM, or standing electronic waves and carrier-mediated exchange interactions in STM, all superimposed onto scanning system distortions and gradual changes of contrast due to drift and/or mis-tilt effects. Correspondingly, while the human eye can readily identify certain patterns in the images such as lattice periodicities, repeating structural elements, or microstructures, their automatic extraction and classification are highly non-trivial and universal pathways to accomplish such analyses are absent. We pose that the most distinctive elements of the patterns observed in STM and (S)TEM images are similarity and (almost-) periodicity, behaviors stemming directly from the parsimony of elementary atomic structures, superimposed on the gradual changes reflective of order parameter distributions. However, the discovery of these elements via global Fourier methods is non-trivial due to variability and lack of ideal discrete translation symmetry. To address this problem, we explore the shift-invariant variational autoencoders (shift-VAEs) that allow disentangling characteristic repeating features in the images, their variations, and shifts that inevitably occur when randomly sampling the image space. Shift-VAEs balance the uncertainty in the position of the object of interest with the uncertainty in shape reconstruction. This approach is illustrated for model 1D data, and further extended to synthetic and experimental STM and STEM 2D data. We further introduce an approach for training shift-VAEs that allows finding the latent variables that comport to known physical behavior. In this specific case, the condition is that the latent variable maps should be smooth on the length scale of the atomic lattice (as expected for physical order parameters), but other conditions can be imposed. The opportunities and limitations of the shift VAE analysis for pattern discovery are elucidated.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Robust Thermal Control of Propulsion Lines for Space Missions

In spacecraft that have propulsion lines that are located externally with open bus architecture, the lines are typically insulated by Multi Layer Insulation (MLI) blankets to protect them thermally from the cold space environment. In addition to heat loss through the insulation, mechanical supports used to attach the lines to the spacecraft structure also create heat leaks from the lines. These lines typically have very low thermal conduction in the axial direction, so the heat balance in the lines tends to be very local without much heat spreading. The typical allowable temperature range for hydrazine-based lines is +15/+50°C. This tight temperature range has to be maintained for every location on these lines. For typical spacecraft, these lines can be several meters long. Temperature control is typically achieved by closed loop monitoring of temperatures along the lines and the corresponding powering of the heaters in a bang-bang approach to maintain the temperatures within the dead band of the control loop. The temperatures of propulsion lines are a function of several parameters with heat loss characteristics of the MLI being the key one. Unfortunately, this same key characteristic (MLI effective emittance) has a large variation along its length due to its dependence on workmanship, which in turn leads to large uncertainties in the propulsion lines’ local temperatures. Because of the poor conduction along the axial direction, heat balance along the length varies dramatically from one location to the next, even few inches apart, depending on the combination of the controlling parameters. This paper describes various robust design and implementation approaches that have been investigated to greatly reduce the randomness associated with predicting the temperature of these propulsion lines.

Bhandari, Pradeep↗
NASA NTRS2002

Better Weather Prediction and Climate Diagnostics Using Rainfall Measurements from Space

Progress in understanding of the role of water in global weather and climate is currently limited by our knowledge of the spatial and temporal variability of primary hydrological fields such as precipitation and evaporation. The Tropical Rainfall Measuring Mission (TRMM) has recently demonstrated that use of microwave-based rainfall observations from space in data assimilation can provide better climate data sets and improve short-range weather forecasting. At NASA, we have been exploring non-traditional approaches to assimilating TRMM Microwave Imager (TMI) and Special Sensor Microwavehager (SSM/I) surface rain rate and latent heating profile information in global systems. In this talk we show that assimilating microwave rain rates using a continuous variational assimilation scheme based on moisture tendency corrections improves quantitative precipitation estimates (QPE) and related clouds, radiation energy fluxes, and large-scale circulations in the Goddard Earth Observing System (GEOS) reanalyses. Short-range forecasts initialized with these improved analyses also yield better QPE scores and storm track predictions for Hurricanes Bonnie and Floyd. We present a status report on current efforts to assimilate convective and stratiform latent heating profile information within the general variational framework of model parameter estimation to seek further improvements. Within the next 5 years, there will be a gradual increase in microwave rain products available from operational and research satellites, culminating to a target constellation of 9 satellites to provide global rain measurements every 3 hours with the proposed Global Precipitation Measurement (GPM) mission in 2007/2008. Based on what has been learned from TRMM, there is a high degree of confidence that these observations can play a'major role in improving weather forecasts and producing better global datasets for understanding the Earth's water and energy cycle. The key to success is to adopt an integrated approach to retrieval, validation, modeling, and data assimilation in a coordinated end-to-end observation-application program.

Hou, Arthur↗
DOE OSTI2020

A non-orthogonal variational quantum eigensolver

Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. This strategy doubles the size of the state preparation circuits but not their depth, while adding a small number of additional two-qubit gates relative to standard variational quantum eigensolver. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H4 and the π-system of hexatriene (C6H8).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

OR-AGENT framework – Architecting electrified heavy-duty drayage applications

The widespread adoption of zero-emission vehicles in heavy-duty (HD) commercial freight transportation faces considerable technoeconomic challenges. For heavy-duty trucks, ensuring high uptime, cost parity with diesel, and safety standards is especially critical as these vehicles operate over long distances with heavy loads, where any downtime or off-nominal behaviors significantly impacts logistics, productivity, and the total cost of ownership. Unlike traditional diesel refueling, BEV charging infrastructure must be co-optimized with vehicle deployment, operational demands, and grid capacity to ensure cost-effective and reliable freight operations. However, the lack of a standardized ownership and service model has led to a fragmented approach—where commercial vehicle operators may invest in, own, and maintain both vehicle/batteries and charging/energy infrastructure. This disconnect may exclude energy service providers from the equation, forcing fleet operators to explore ‘behind-the-fence’ energy solutions that increase capital investment, operational downtime, overhead costs, and, in some cases, net carbon emissions. To address these issues, this study introduces OR-AGENT (Optimal Regional Architecture Generation for Efficient National Transport), a comprehensive modeling framework that integrates powertrain architectures, charging infrastructures, and energy backbone systems into a cohesive strategy. In this paper, OR-AGENT is applied to develop an interconnected systems architecture for energy efficiency and resiliency enhancement of heavy-duty drayage vehicles at the Port of Savannah, GA. This framework showcases an interconnected systems approach to electrifying heavy-duty drayage vehicles at the Port of Savannah, GA. The study assessed BEVs with 400–1200 kWh battery capacities, accounting for seasonal variations in weather and freight routing. A diverse charging mix (150 kW–1250 kW) was evaluated alongside grid capacity constraints, cost, and carbon intensity analysis, leading to the development of a strategic microgrid/Distributed Energy Resources (DER) deployment architecture to ensure a reliable and sustainable transition. However, the findings also highlight the need for alternative zero-emission solutions for remaining trips, such as larger batteries, electrified roadways, hydrogen powertrains, or net-zero emission fuels. In conclusion, the findings are incorporated into a Total Cost of Ownership (TCO) model to identify optimal architectures for an interconnected electrified ecosystem.

Commercial vehicles↗
DOE OSTI2023

A state-of-the-art review of experimental and computational studies of granular materials: Properties, advances, challenges, and future directions

Modeling of heterogeneous materials and media is a problem of fundamental importance to a wide class of phenomena and systems, ranging from condensed matter physics, soft materials, and composite media to porous media, biological systems, geosystems, ceramic engineering, pharmaceutical science and even in space discoveries. Among the most important materials are granular systems, which have received intense interest from the engineering, physics, and mathematics communities. In this review paper, the recent developments and new advances in experimental, and computational methods on a variety of scales and physics that extend understanding to a wide range of materials and phenomena are reviewed. Experimental advances include computed neutron and nanometer-scale tomography, magnetic resonance imaging, refractive index matching, digital image correlation, acoustic emission analysis, and the most recent 4D techniques. Furthermore, a tremendous shift has occurred from the continuum scale to micro-scale and developing multiscale approaches. As such, various computational methods, including, constitutive modeling, discrete modeling, and multiscale approaches, have been developed. In conclusion, aside from all these evolutions, more complicated modeling called coupled, or multiphysics, systems representing a simultaneous presence of heat, fluid, chemical variation, and mechanical effect are also explored.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probing the Mechanism of Selective Phosphate Adsorption from Wastewater Using Aqueous and Synchrotron X-ray Characterization

Ion exchange shows promise for recovering phosphate from wastewater as value-added products, but requires high phosphate selectivity to compete with conventional treatment. Hybrid anion exchange (HAIX) resins, which contain non-selective basic functional groups and selective iron oxide nanoparticles (FeOnp), can effectively remove phosphate from wastewater. However, knowledge gaps remain regarding the mechanisms of phosphate selectivity and influence of competing ions, hindering needed efforts to model adsorption dynamics and design scalable adsorption processes for varying wastewaters. To address these gaps, we integrated aqueous-phase adsorption analysis with solid-phase, synchrotron-based X-ray characterization; this integration facilitated elucidation of the distribution and speciation of iron, phosphate, and competing anions on HAIX resins. We compared a quaternary ammonium-functionalized HAIX resin (SBA) to a tertiary amine version (WBA) to determine the role of functional groups. X-ray radiography revealed differences in FeOnp speciation (goethite vs. ferrihydrite) and distribution (peripheral vs. homogeneous) between the resins, resulting in varied phosphate affinity and intraparticle diffusion resistance. Using micro-X-ray fluorescence (μ-XRF) and micro-X-ray absorption near-edge structure (μ-XANES) spectroscopy, we identified differences in where and how phosphate binds across resin types and wastewaters. Across wastewater compositions, FeOnp sites in WBA contribute more to phosphate adsorption than in SBA, possibly due to variations in Fe distribution and speciation. Phosphate adsorption densities calculated from quantitative μ-XRF maps matched those from aqueous analysis, demonstrating the effectiveness of this integrated approach. Altogether, results demonstrate the use of synchrotron-based X-ray characterization for investigating adsorption mechanisms and advance HAIX as a phosphate recovery technology from complex wastewaters.

Nutrient recovery↗
DOE OSTI2025

Updates on MURAVES Project at Mt. Vesuvius

The MUon RAdiography of VESuvius (MURAVES) project aims to employ muography imaging techniques to investigate the internal structure of the summit of Mount Vesuvius, an active volcano located near Naples, Italy. This paper reports recent advancements in data analysis and simulation tools that significantly improve the quality and reliability of the experiment’s results. A new track selection method, referred to as the Golden Selection, has been developed to identify high-quality muon tracks by applying an improved χ 2 -based criterion. This method enhances the signal-to-background ratio and improves the resolution of the resulting muographic images. Moreover, the simulation framework has been upgraded through the integration of the MULDER (MUon simuLation for DEnsity Reconstruction) library, which consolidates the functionalities of previously used libraries into a single, unified platform. MULDER enables efficient and accurate modeling of muon flux variations induced by topographical features. A good agreement is observed between the simulated and measured muon flux maps, validating the effectiveness of the new analysis and simulation approaches.

Cosmic rays↗
DOE OSTI2022

Inferring safety critical events from vehicle kinematics in naturalistic driving environment: Application of deep learning Algorithms

Advances in sensing technology has enabled the collection of countless terabytes of second-by-second kinematics data. Such data provides opportunities for real-time monitoring of driving behavior and identification of safety critical events (SCEs) including crashes and near crashes. The concept of volatility is relevant in this context, which identifies instability and erratic variations in driving behavior prior to involvement in SCEs. This study utilized vehicle kinematics from a large-scale naturalistic driving data to develop a deep learning approach based on 1D convolutional neural networks (CNN) for inferring SCEs. The data are unique in the sense that such accurate pre-crash data at high fidelity are not available in traditional crash repositories. This study contributes to the literature by providing a first attempt at predicting responses to SCEs by developing deep learning-based CNN architectures using novel driving volatility based kinematic thresholds for a sample of 9553 events. The key contribution lies in developing a volatility-based CNN input layout that is acceptable to CNN schemes and represents the motion kinematics such as speed, acceleration and volatility measures. Several 1D-CNN architectures were developed using layers, numbers of convolutions, layer patterns, and kernels. Shallow and deep architectures were tested, revealing higher accuracy of shallow architectures in detecting SCEs. The optimal number of epochs were identified using an early stopping method while the CNN performance was improved by increasing the number of epochs. The ensemble CNN had the highest predictive accuracy of 95.6% for detection of crashes and near crashes, which was 2.5% higher than the optimal CNN using 20% hold out test data. The ensemble CNN also outperformed classical machine learning models and model performance reported in past studies on detection of SCEs. Finally, these results have implications for identification of safety hotspots and providing real-time alerts and warnings in connected and highly automated vehicle environment including society of automotive engineers levels 3–5.

42 ENGINEERING↗
DOE OSTI2025

Improved heavy-ion PID using scintillation light detector with neural network analysis: a Monte Carlo simulation study

The photon collection efficiency of gaseous scintillator detectors varies according to the position of the impinging charged particles in the medium that generates scintillation light. Thus, when impinging particles are distributed over a large area, the intrinsic photon-number resolution of the system is affected by a large variation. This work presents and discusses a method for adjusting the total number of detected photons to account for variation in the photon collection efficiency as a function of the position of the light source within the scintillating medium. The method was developed and validated by processing data from systematic simulation studies based on GEANT4 that model the response of the Energy Loss Optical Scintillation System (ELOSS) detector. The position of the charged particle is calculated using a deep neural network algorithm. This is accomplished by analyzing the distribution of scintillation light recorded by the array of photosensors. The estimated particle position is then used to calculate the correction factor and adjust the amount of captured light to account for variations in the photon collection efficiency. The neural network algorithm provides excellent tracking capabilities, achieving sub-millimeter position resolution and an angular resolution of 12 mrad, approaching the performance of traditional tracking detectors (e.g., drift chambers). The present method can be generalized to any optical scintillation system where the photon collection efficiency depends on the position of the impinging particle.

Heavy-ion detectors↗
DOE OSTI2024

Two-dimensional coherent spectrum of high-spin models via a quantum computing approach

Here in this work we present and benchmark a quantum computing approach to calculate the two-dimensional coherent spectrum (2DCS) of high-spin models. Our approach is based on simulating their real-time dynamics in the presence of several magnetic field pulses, which are spaced in time. We utilize the adaptive variational quantum dynamics simulation algorithm for the study due to its compact circuits, which enables simulations over sufficiently long times to achieve the required resolution in frequency space. Specifically, we consider an antiferromagnetic quantum spin model that incorporates Dzyaloshinskii-Moriya interactions and single-ion anisotropy. The obtained 2DCS spectra exhibit distinct peaks at multiples of the magnon frequency, arising from transitions between different eigenstates of the unperturbed Hamiltonian. By comparing the one-dimensional coherent spectrum with 2DCS, we demonstrate that 2DCS provides a higher resolution of the energy spectrum. We further investigate how the quantum resources scale with the magnitude of the spin using two different binary encodings of the high-spin operators: the standard binary encoding and the Gray code. At low magnetic fields both encodings require comparable quantum resources, but at larger field strengths the Gray code is advantageous. Numerical simulations for spin models with increasing number of sites indicate a polynomial system-size scaling for quantum resources. Lastly, we compare the numerical 2DCS with experimental results on a rare-earth orthoferrite system. The observed strength of the magnonic high-harmonic generation signals in the 2DCS of the quantum high-spin model aligns well with the experimental data, showing significant improvement over the corresponding mean-field results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Geographical Insights into Suicide Mortality Through Spatial Machine Learning

Suicide mortality is a leading cause of death in the United States, with an upward trend that emphasizes its significance as a public health issue. Previous research has employed global models like ordinary least squares (OLS) regression and local models such as geographically weighted regression (GWR). While local models are useful for analyzing spatial variations in suicide mortality, they share limitations with traditional global models, particularly about their inability to handle multi-collinearity and non-linear relationships. Machine learning approaches, like random forests (RF), can address some of these limitations but often fail to account for spatial variability. This gap highlights the need for spatial ML models specifically designed to tackle suicide mortality. This research seeks to fill this void by using a geographically weighted random forest model (GWRF) to examine the associations between county-level suicide mortality in the U.S. from 2010 to 2020 and various social and environmental determinants of health. A key aspect of our methodology is disciplined feature selection, which reduces the pool of explanatory variables by about 90%. This refinement enhances the explanatory power of both global (R2 improved from 0.59 to 0.67) and local (R2 improved from 0.64 to 0.67) RF models while reducing their run times. An analysis of the importance scores for these selected features reveals that the drivers of suicide mortality vary by context. Thus, to effectively address regional disparities and inform targeted public health interventions, a holistic approach that incorporates multiple county-level characteristics is essential.

Lebakula, Viswadeep [ORNL] (ORCID:0000000152935914↗