Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

USA Crop Yield Estimation with MODIS NDVI: Are Remotely Sensed Models Better Than Simple Trend Analyses?

Crop yield forecasting is performed monthly during the growing season by the United States Department of Agriculture’s National Agricultural Statistics Service. The underpinnings are long-established probability surveys reliant on farmers’ feedback in parallel with biophysical measurements. Over the last decade though, satellite imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) has been used to corroborate the survey information. This is facilitated through the Global Inventory Modeling and Mapping Studies/Global Agricultural Monitoring system, which provides open access to pertinent real-time normalized difference vegetation index (NDVI) data. Hence, two relatively straightforward MODIS-based modeling methods are employed operationally. The first model constitutes mid-season timing based on the maximum peak NDVI value, while the second is reflective of late-season timing by integrating accumulated NDVI over a threshold value. Corn model results nationally show the peak NDVI method provides a R^(2) of 0.88 and a coefficient of variation (CV) of 3.5%. The accumulated method, using an optimally derived 0.58 NDVI threshold, improves the performance to 0.93 and 2.7%, respectively. Both these models outperform simple trend analysis, which is 0.48 and 7.4%, correspondingly. For soybeans the R^(2) results of the peak NDVI model are 0.62, and 0.73 for the accumulated using a 0.56 threshold. CVs are 6.8% and 5.7%, respectively. Spring wheat’s R2performance with the accumulated NDVI model is 0.60 but just 0.40 with peak NDVI. The soybean and spring wheat models perform similarly to trend analysis. Winter wheat and upland cotton show poor model performance, regardless of method. Ultimately, corn yield forecasting derived from MODIS imagery is robust, and there are circumstances when forecasts for soybeans and spring wheat have merit too.

crop yield↗
DOE OSTI2025

Applying Linear Optics from Closed Orbits Modulation for finding beam-based alignment of harmonic sextupoles

A fast and accurate beam-based alignment (BBA) method for harmonic sextupoles has been devel oped at NSLS-II using Linear Optics from Closed Orbit Modulation (LOCOM). The approach excites the beam with simultaneous sine-wave signals at two fast correctors, chosen with an appropriate phase advance to span the full betatron phase space. This strategy suppresses systematic errors from hysteresis, while additional errors from orbit drift and power-supply calibration are minimized by the short measurement time (a few minutes) and reliance solely on beam-based current-to-field conver sion of the sextupoles, with hysteresis explicitly included. Simulations indicate that Linear Optics from Closed Orbits (LOCO) combined with 0.5 mm local orbit bumps can resolve relative sextupole field offsets (∆k₂) with precision better than 10% of k₂, reflecting to beam-based alignment (BBA) accuracy finer than 50 µm. Moreover, employing machine-learning-optimized local orbit bumps en hances sextupole-induced quadrupole signals in a deterministic manner and maintains orbit stability under large sextupole strength variations (±40%). The proposed method is experimentally validated through proof-of-principle measurements at NSLS-II, demonstrating its potential as a fast, precise, and robust tool for harmonic sextupole alignment.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

CONCURRENT, CONDENSED STEIN VARIATIONAL GRADIENT DESCENT FOR UNCERTAINTY QUANTIFICATION OF NEURAL NETWORKS

In this work, we propose a Stein variational gradient descent (SVGD) method to concurrently sparsify, train, and provide uncertainty quantification (UQ) of a complexly parameterized model, such as a neural network (NN). It employs a graph reconciliation and condensation process to reduce complexity and increase similarity in the Stein ensemble of parameterizations. Therefore, the proposed concurrent, condensed SVGD (ccSVGD) method can provide UQ on parameters, not just outputs. Furthermore, the parameter reduction speeds up the convergence of the Stein gradient descent as it reduces the combinatorial complexity by aligning and differentiating the sensitivity to parameters. These properties are demonstrated with an illustrative example and an application to a mechanical response representation problem in solid mechanics.

42 ENGINEERING↗
DOE OSTI2024

Convergence of variational Monte Carlo simulation and scale-invariant pre-training

We provide theoretical convergence bounds for the variational Monte Carlo (VMC) method as applied to optimize neural network wave functions for the electronic structure problem. Here, we study both the energy minimization phase and the supervised pre-training phase that is commonly used prior to energy minimization. For the energy minimization phase, the standard algorithm is scale-invariant by design, and we provide a proof of convergence for this algorithm without modifications. The pre-training stage typically does not feature such scale-invariance. We propose using a scale-invariant loss for the pretraining phase and demonstrate empirically that it leads to faster pre-training.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

Determination of efficiencies, loss mechanisms, and performance degradation factors in chopper controlled dc vehical motors. Section 2: The time dependent finite element modeling of the electromagnetic field in electrical machines: Methods and applications

The time dependent solution of the magnetic field is introduced as a method for accounting for the variation, in time, of the machine parameters in predicting and analyzing the performance of the electrical machines. The method of time dependent finite element was used in combination with an also time dependent construction of a grid for the air gap region. The Maxwell stress tensor was used to calculate the airgap torque from the magnetic vector potential distribution. Incremental inductances were defined and calculated as functions of time, depending on eddy currents and saturation. The currents in all the machine circuits were calculated in the time domain based on these inductances, which were continuously updated. The method was applied to a chopper controlled DC series motor used for electric vehicle drive, and to a salient pole sychronous motor with damper bars. Simulation results were compared to experimentally obtained ones.

Hamilton, H. B.↗
NASA NTRS1984

A method for interface shape control during Bridgman type crystal growth of HgCdTe alloys

A method is described for interface shape control of Hg(1-x)Cd(x)Te alloys during Bridgman-type crystal growth. The authors have used this method for several years to obtain the flat or convex interface shapes required for simultaneously reducing radial compositional variations and crystal defects. The method turns to advantage the thick-walled, fused-silica ampules required for growing these alloys. A proper combination of furnace geometry (two isothermal zones separated by a thin barrier), upper- and lower-zone temperatures, and growth rate is required. Examples of results are reported.

Szofran, F. R.↗
NASA NTRS1987

Equation-of-state spinning fluids in the Einstein-Cartan theory

The relativistic fluid equations may be completed in two physically distinct methods. One method assumes the mass, rho, (or particle number) is conserved, while the other method assumes an equation of state of the form P = P(rho). A variational principle for the mass conservation method both with and without an intrinsic spin for the fluid was constructed earlier (Ray and Smalley, 1982 and 1983). A variational principle for the fluid described by an equation of state both with and without spin is formulated. In all cases the variational principle is set in the Einstein-Cartan metric-torsion U4 geometry. The results for general relativity follow as a special case.

Ray, John R.↗
DOE OSTI2023

BUQEYE guide to projection-based emulators in nuclear physics

The BUQEYE collaboration (Bayesian Uncertainty Quantification: Errors in Your effective field theory) presents a pedagogical introduction to projection-based, reduced-order emulators for applications in low-energy nuclear physics. The term emulator refers here to a fast surrogate model capable of reliably approximating high-fidelity models. As the general tools employed by these emulators are not yet well-known in the nuclear physics community, we discuss variational and Galerkin projection methods, emphasize the benefits of offline-online decompositions, and explore how these concepts lead to emulators for bound and scattering systems that enable fast and accurate calculations using many different model parameter sets. We also point to future extensions and applications of these emulators for nuclear physics, guided by the mature field of model (order) reduction. All examples discussed here and more are available as interactive, open-source Python code so that practitioners can readily adapt projection-based emulators for their own work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

Variable Sampling Mapping

The performance of an optical system (for example, a telescope) is limited by the misalignments and manufacturing imperfections of the optical elements in the system. The impact of these misalignments and imperfections can be quantified by the phase variations imparted on light traveling through the system. Phase retrieval is a methodology for determining these variations. Phase retrieval uses images taken with the optical system and using a light source of known shape and characteristics. Unlike interferometric methods, which require an optical reference for comparison, and unlike Shack-Hartmann wavefront sensors that require special optical hardware at the optical system's exit pupil, phase retrieval is an in situ, image-based method for determining the phase variations of light at the system s exit pupil. Phase retrieval can be used both as an optical metrology tool (during fabrication of optical surfaces and assembly of optical systems) and as a sensor used in active, closed-loop control of an optical system, to optimize performance. One class of phase-retrieval algorithms is the iterative transform algorithm (ITA). ITAs estimate the phase variations by iteratively enforcing known constraints in the exit pupil and at the detector, determined from modeled or measured data. The Variable Sampling Mapping (VSM) technique is a new method for enforcing these constraints in ITAs. VSM is an open framework for addressing a wide range of issues that have previously been considered detrimental to high-accuracy phase retrieval, including undersampled images, broadband illumination, images taken at or near best focus, chromatic aberrations, jitter or vibration of the optical system or detector, and dead or noisy detector pixels. The VSM is a model-to-data mapping procedure. In VSM, fully sampled electric fields at multiple wavelengths are modeled inside the phase-retrieval algorithm, and then these fields are mapped to intensities on the light detector, using the properties of the detector and optical system, for comparison with measured data. Ultimately, this model-to-data mapping procedure enables a more robust and accurate way of incorporating the exit-pupil and image detector constraints, which are fundamental to the general class of ITA phase retrieval algorithms.

Smith, Jeffrey, S.↗
DOE OSTI2020

A Facile Approach to High Precision Detection of Cell-to-Cell Variation for Li-ion Batteries

Abstract Over the past decade, it has been repeatedly demonstrated that homogeneity in electrochemical performance of lithium-ion cells plays a major role in determining the life and safety of lithium-ion battery modules or packs. Generally, the homogeneity of a battery pack is evaluated by characterizing the cells individually in terms of capacity, mass, impedance. Particularly, high quality electrochemical data heavily relies on the availability of high precision current source to minimize the discrepancy induced by the channel-to-channel variation. Here, a facile and precise measurement method is reported for screening cell-to-cell variations, in which voltage is the only indicator parameter independent of high precision current source. In detail, by connecting the cells in series (CiS), the measurement error of electrochemical data caused by stability and discrepancy of current sources among different charge/discharge equipment can be effectively avoided. The findings of this work showed that the cell-to-cell variations can be simply and sensitively detected with CiS configuration. For example, the relative standard deviation, which is the evaluation criterion of battery homogeneity, was 2.14% based on CiS while it was 0.43% based on individual measurements. The simple and precise CiS measurement is promising for evaluation of cell quality or module integration quality. In addition, this work can also provide a solid foundation for the development of detection algorithms for battery management systems to rapidly monitor battery homogeneity.

25 ENERGY STORAGE↗
DOE OSTI2022

Predicting variations of the least principal stress with depth: Application to unconventional oil and gas reservoirs using a log-based viscoelastic stress relaxation model

Knowledge of layer-to-layer variations of the least principal stress, S hmin , with depth is essential for optimization of multi-stage hydraulic fracturing in unconventional reservoirs. Utilizing a geomechanical model based on viscoelastic stress relaxation in relatively clay rich rocks, we present a new method for predicting continuous S hmin variations with depth. The method utilizes geophysical log data and S hmin measurements from routine diagnostic fracture injection tests (DFITs) at several depths for calibration. We consider a case study in the Wolfcamp formation in the Midland Basin, where both geophysical logs and values of S hmin from DFITs are available. We compute a continuous stress profile as a function of the well logs that fits all of the DFITs well. We utilized several machine learning technologies, such as bootstrap aggregation (or bagging), to improve the generalization of the model and demonstrate that the excellent fit between predicted and observed stress values is not the result of over-fitting the calibration points. The model is then validated by accurately predicting hold-out stress measurements from four wells within the study area and, without recalibration, accurately predicting stress as a function of depth in an offset pad about 6 miles away.

58 GEOSCIENCES↗
NASA NTRS2019

The Role of Thermodynamic Phase Shifts in Cloud Optical Depth Variations with Temperature

We present a novel method that identifies the contributions of thermodynamic phase shifts and processes governing supercooled liquid and ice clouds to cloud optical depth variations with temperature using MODIS observations. Our findings suggest that thermodynamic phase shifts outweigh the net influence of processes governing supercooled liquid and ice clouds in causing increases in mid-latitudinal cold cloud optical depth with temperature. Cloud regime analysis suggests that dynamical conditions appear to have less influence on the contribution of thermodynamic phase shifts to cloud optical depth variations with temperature. Thermodynamic phase shifts also contribute more to increases in cloud optical depth during colder seasons due to the enhanced optical thickness contrast between liquid and ice clouds. The results of this study highlight the importance of thermodynamic phase shifts in explaining cold cloud optical depth increases with temperature in the current climate and may elucidate their role in the cloud optical depth feedback.

cloud optical depth↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Uranium Oxide Synthetic Pathway Discernment through Unsupervised Morphological Analysis

We present a novel unsupervised machine learning method for quantitative representation of scanning electron micrographs and its applications and performance for nuclear forensic analysis of uranium ore concentrates. The method uses a vector quantizing variational autoencoder followed by a histogram operation to encode a micrograph into a single dimensional representation, called the latent vector. The method requires no extant labeling of the data and can be applied over large datasets of micrographs with minimal human interaction. The representations generated are broadly descriptive of each micrograph and the microstructure of the material imaged. In the case of uranium ore concentrate analysis, the representations were amenable to processing reagent and ore concentrate species classification with accuracy of 81:8%, which is competitive with state-of-the-art supervised networks. The representations were also used to classify previously unseen processing routes, were able to classify imaging parameters such as magnification (to 76:0% accuracy), were able to classify fine grained process parameters such as calcining temperature (to 74:4% accuracy), and their informatic properties indicate that they are generally descriptive of the image represented. This method can be applied across microstructure analysis fields to perform quantitative analysis without the need for labor intensive and possibly biased human analysis.

Scanning Electron Microscopy, Vector Quantizing Va↗
NASA NTRS1957

Rocket thrust variation with foamed liquid propellants

An analysis is presented on a method for varying rocket thrust by varying the bulk density of the propellants. This density variation was accomplished by uniformly dispersing an inert, insoluble gas in the liquid propellants. Only qualitative agreement with theory was obtained from preliminary experiments with a 1000-pound-thrust ammonia - nitric acid rocket engine; the required experimental gas-flow rates were two to six times greater than those predicted by theory. It was demonstrated, however, that this method of rocket-thrust variation is feasible.

Morrell, G↗
DOE OSTI2024

Contrasting roles of rice root iron plaque in retention and plant uptake of silicon, phosphorus, arsenic, and selenium in diverse paddy soils

Abstract Background and aims Iron (Fe) plaque on rice roots is a mixture of Fe oxide and oxyhydroxide minerals thought to protect rice from high levels of arsenic (As) in flooded paddy soils. Silicon (Si), phosphorus (P), and selenium (Se) also exist as oxyanions in rice paddies, but the impacts of Fe plaque on uptake of these nutrients are unknown. Methods We used natural variation in paddy soil chemistry to test how Si, P, As, and Se move from porewater to plaque to plant via multiple techniques. In a pot study, we monitored Fe plaque deposition and porewater chemistry in 5 different soils over time and measured plaque/plant chemistry and Fe plaque mineralogy at harvest. We normalized oxyanion concentrations by Fe to determine the preferential retention on plaque or plant uptake. Results Low phosphorus availability increased root Fe-oxidizing activity, while Fe, Si, P, As, and Se concentrations in plaque were strongly correlated with porewater. Plaque did not appreciably retain Si and Se, and the oxyanions did not compete for adsorption sites on the Fe plaque. Conclusion Root Fe plaque seems to protect rice from As uptake, does not interfere with Si and Se uptake, and roots adapt to maintain P nutrition even with retention of porewater P on plaque.

Agriculture↗
DOE OSTI2025

Theory and numerics of subspace approximation of eigenvalue problems

Large-scale eigenvalue problems arise in various fields of science and engineering and demand computationally efficient solutions. In this study, we investigate the subspace approximation for parametric linear eigenvalue problems, aiming to mitigate the computational burden associated with high-fidelity systems. Furthermore, we provide general error estimates under non-simple eigenvalue conditions, establishing some theoretical foundations for understanding the convergence behavior of subspace approximations. Numerical examples, including problems with one-dimensional to three-dimensional spatial domain and one-dimensional to two-dimensional parameter domain, are presented to demonstrate the efficacy of reduced basis method in handling parametric variations in boundary conditions and coefficient fields to achieve significant computational savings while maintaining high accuracy, making them promising tools for practical applications in large-scale eigenvalue computations.

Eigenvalue problems↗