SEARCH · Search NASA

Results for “linear optimization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗

DOE OSTI2026

A Discrete Hankel Transform Approach to Nuclear Data Processing for Fusion Applications

This study introduces advancements to the numerical solutions employed in the processing of nuclear data for fusion applications. It leverages the convolution theorem and Fourier transform techniques to enhance computational efficiency and broaden applicability. Building upon a previously reported discrete Hankel transform approach for Doppler broadening, this work refines the solution of convolution integrals central to these applications. The methodology provides a general and unified framework for evaluating any convolution operation, regardless of whether the underlying problem involves temperature effects in nuclear reactions. The applicability to the nuclear data processing for fusion is demonstrated by deriving the convolution integrals for some of the fusion-related quantities. As before, the convolution operation utilizes a Gaussian-based kernel; however, the discrete Hankel transform of order $𝛼$ = $\frac{1}{2}$ is now applied to the forward Fourier transform of the nonkernel argument, rather than the inverse Fourier transform. This modification eliminates the need for the integration of the nonkernel, cross section–based function, which is a step that posed challenges for certain pointwise cross-section representations. It also removes the requirement for cross-section linearization. Optimized for graphics processing unit architectures, the approach significantly improves computational performance. These advancements are currently under evaluation as the foundation for the next-generation thermonuclear data file processing codes being developed at Lawrence Livermore National Laboratory.

Nuclear science and engineering↗

DOE OSTI2021

Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with 𝑚 𝑢⁢𝑑 ≈ 𝑚$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of 𝑎 ≈0.086 fm, lattice volume of 𝐿 = 48⁢𝑎 ≈ 4.1 fm and pion mass 𝑚 𝜋 ≈ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Foreground-immune CMB lensing reconstruction with polarization

Extragalactic foregrounds are known to generate significant biases in temperature based CMB lensing reconstruction. Several techniques, which include "source hardening'" and "shear-only estimators'' have been proposed to mitigate contamination and have been shown to be very effective at reducing foreground-induced biases. Here we extend both techniques to polarization, which will be an essential component of CMB lensing reconstruction for future experiments and investigate the "large lens'' limit analytically to gain insight on the origin and scaling of foreground biases, as well as the sensitivity to their profiles. Furthermore, using simulations of polarized point sources, we estimate the expected bias to both Simons Observatory and CMB-S4 like (polarization-based) lensing reconstruction, finding that biases to the former are minuscule while those to the latter are potentially non-negligible at small scales ($L$ ~ $1000-2000$). In particular, we show that for a CMB-S4 like experiment, an optimal linear combination of point-source hardened estimators can reduce the (point-source induced) bias to the CMB lensing power spectrum by up to two orders of magnitude, at a ~ 4% noise cost relative to the global minimum variance estimator.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Resilient Operation of Power Distribution Systems Using MPC-Based Critical Service Restoration

Power distribution systems are more prone to disruptions and cause most power system outages. We propose a service restoration technique to recover the system service (electricity delivery) following an extreme event-triggered substation outage. The proposed technique considers the problem of controlling distributed energy resources (DERs) of a distribution system with the objective of achieving maximum load pick up while satisfying network flow and voltage constraints. The problem is formulated as a model predictive control (MPC), where a linearized optimal power flow (OPF) model is employed to describe the network. The formulation is augmented with a ramping (up) reserve product for the DERs to ensure an upward monotonic load restoration as time evolves. We perform simulations considering the IEEE 13-bus test feeder integrated with wind, solar, microturbine, and energy storage devices. We demonstrate the efficacy of the devised technique in restoring the system loads monotonically, without shedding previously restored loads. We also show the benefit of co-optimization of power and reserve products for DERs on service restoration. In addition, the capability of the technique in regulating nodal voltages and reducing renewable power curtailment is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

How state transitions balance photosynthetic electron transport in plants – a quantitative study

In plants, the process of state transition regulates the allocation of sunlight energy between Photosystem II (PSII) and PSI. However, the implications of state transitions for harmonizing electron transport rates between photosystems, and a full quantitative picture of this process, remain underexplored. We integrated quantitative biology (biochemical and biophysical approaches) with in vivo spectroscopy on wild-type Arabidopsis and protein phosphorylation mutants. This combination facilitated monitoring of Chl redistribution and its functional implications for light harvesting and electron transport. Our findings demonstrate the reallocation of 12% of highly phosphorylated ‘extra’ light-harvesting complex II under state 2 from stacked to unstacked thylakoids. This reduces the number of Chls per PSII from 216 to 182, while increasing the number in PSI from 187 to 223. Such Chl redistribution compensates for differences in photosystem stoichiometry and photochemical quantum efficiencies, thereby precisely synchronizing electron transport rates in both photosystems. Mutant analyses corroborate that this regulatory mechanism involves reversible phosphorylation. We inferred that state transitions optimize linear electron transport, leaving no additional capacity for cyclic electron transport. Furthermore, the results suggest that the controversies about long-range migration of LHCII from stacked to unstacked thylakoid domains arise from differences in phosphorylation levels.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

GET-Solar (Generation forecasting and parameter Estimation Tool for Solar) [SWR-20-115]

GET-Solar (Generation forecasting and parameter Estimation Tool for Solar) can be used to estimate behind the meter (BTM) solar generation from limited datasets. Using distributed sets of diverse solar sites, but with close spatial proximity, and/or limited historical PV generation information can be used to generate of solar profiles. This is useful for applications where partial solar generation data is available, and on-site irradiance data is not available. GET-Solar calculates PV technical data parameters and calculates output by modeling shading and cloud cover impacts. Clear sky irradiance profiles are used to calculate base generation profiles on which the local shading and cloud cover layers are added. In order to arrive at clear sky generation data, PV mounting information like panel tilt and and angle is needed. For cases where PV mounting information is missing, panel tilt and azimuth angle can be estimated using a non-linear optimization algorithm. GET-Solar can also be used for forecasting applications where historical PV data is available to estimate local shading at various solar positions. This overall framework can be used for estimation of solar parameters, and assisting in analyzing and modeling diverse sets of DER where solar generation is also present in utility meter data. GET-Solar has been developed in Python and uses the Pyomo optimization language.

Abraham, Sherin Ann↗

DOE OSTI2023

TwoFold: Highly accurate structure and affinity prediction for protein-ligand complexes from sequences

We describe our development of ab initio protein-ligand binding pose prediction models based on transformers and binding affinity prediction models based on the neural tangent kernel (NTK). Folding both protein and ligand, the TwoFold models achieve efficient and quality predictions matching state-of-the-art implementations while additionally reconstructing protein structures. In conclusion, solving NTK models points to a new use case for highly optimized linear solver benchmarking codes on HPC.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Are coupled renewable-battery power plants more valuable than independently sited installations?

Coupled renewable-battery powerplants differ from the traditional concept of independent siting of electricity resources within transmission networks. Prior research on the value proposition and cost savings from coupling did not consider the geographic constraint of co-location. This paper fills the gap by assessing how pricing volatility differences between nodes within electricity markets impact the system value of coupled renewable-battery projects as compared to independent VRE and battery installations. We use wholesale power market prices from 2012–2019 across the seven main U.S. independent system operators (ISOs) with a linear optimization program to compare the electricity market value of coupled projects to the value of the same underlying sub-components, deployed separately. We find that additional value from adding a 4-hour battery sized to 50% of renewable-plant nameplate capacity is $\$$10/MWh across ISOs on average. The highest boost occurs in California ($\$$15/MWh), where the value of adding storage to solar rises over time in tandem with increased solar penetration in the region. If renewables and batteries are deployed independently, we estimate that $\$$12.5/MWh of additional value could be achieved because of more flexibility on battery siting and operation. The $\$$12.5/MWh coupling penalty is reduced to $\$$1.6/MWh when considering alternative approaches to integrating battery storage. This result implies that renewable-battery power plants will play an increasing role in electricity systems if they can be built for $\$$2–$\$$13/MWh less than independent projects of comparable size. However, the wide regional variation in coupling penalties, along with the importance of conditions captured in our sensitivity cases, suggests the tradeoff between coupling penalties and savings will vary by situation. Therefore, roles exist for independent and coupled projects from a system optimization perspective.

14 SOLAR ENERGY↗

DOE OSTI2021

Keep it short: Exploring the impacts of configuration choices on the recent economics of solar-plus-battery and wind-plus-battery hybrid energy plants

Coupled renewable-battery powerplants differ from the traditional concept of independent siting of electricity resources within transmission networks. Prior research on the value proposition and cost savings from coupling did not consider the geographic constraint of co-location. This paper fills the gap by assessing how pricing volatility differences between nodes within electricity markets impact the system value of coupled renewable-battery projects as compared to independent VRE and battery installations. We use wholesale power market prices from 2012–2019 across the seven main U.S. independent system operators (ISOs) with a linear optimization program to compare the electricity market value of coupled projects to the value of the same underlying sub-components, deployed separately. We find that additional value from adding a 4-hour battery sized to 50% of renewable-plant nameplate capacity is $\$10$/MWh across ISOs on average. The highest boost occurs in California ($\$15$/MWh), where the value of adding storage to solar rises over time in tandem with increased solar penetration in the region. If renewables and batteries are deployed independently, we estimate that $\$12.5$/MWh of additional value could be achieved because of more flexibility on battery siting and operation. The $\$12.5$/MWh coupling penalty is reduced to $\$1.6$/MWh when considering alternative approaches to integrating battery storage. This result implies that renewable-battery power plants will play an increasing role in electricity systems if they can be built for $\$2$–$\$13$/MWh less than independent projects of comparable size. However, the wide regional variation in coupling penalties, along with the importance of conditions captured in our sensitivity cases, suggests the tradeoff between coupling penalties and savings will vary by situation. Therefore, roles exist for independent and coupled projects from a system optimization perspective.

14 SOLAR ENERGY↗

DOE OSTI2023

A General Method for Estimating Zonal Transmission Interface Limits from Nodal Network Data: Preprint

Capacity expansion models for the electric power system often employ zonal (rather than nodal) resolution, necessitating estimates of aggregate power transfer limits across the interfaces between model zones. Interface limits between planning areas are sometimes published, but they are not generalizable to arbitrary zone shapes. There is thus a need for a reproducible method for estimating interface transfer limits (ITLs) between user-defined zones directly from nodal transmission system data. Here, we present a simple method for estimating ITLs using a DC power flow approximation via the power transfer distribution factor (PTDF) matrix. Linear optimization is performed to identify the distribution of power flows that maximizes the total flow on interface-crossing lines, subject to individual line ratings, limits on bus injection/withdrawal, and the relationships among flows, injections, and withdrawals imposed by the PTDF matrix. We demonstrate the application of the method on a 134-zone ~65000-bus system, and we explore the influence of flow direction, contingency level, and zone size on the estimated ITLs. There is significant heterogeneity in the ratio of the ITL to the sum of interface-crossing line ratings, which highlights the importance of accounting for the physical constraints on power flows imposed by Kirchhoff's laws when estimating zonal ITLs.

17 WIND ENERGY↗

DOE OSTI2021

Development of an in situ dual-channel thermal desorption gas chromatography instrument for consistent quantification of volatile, intermediate-volatility and semivolatile organic compounds

Abstract. Aerosols are a source of great uncertainty in radiative forcing predictions and have poorly understood health impacts. Most aerosol mass is formed in the atmosphere from reactive gas-phase organic precursors, forming secondary organic aerosol (SOA). Semivolatile organic compounds (SVOCs) (effective saturation concentration, C*, of 10−1–103 µg m−3) comprise a large fraction of organic aerosol, while intermediate-volatility organic compounds (IVOCs) (C* of 103–106 µg m−3) and volatile organic compounds (VOCs) (C* ≥ 106 µg m−3) are gas-phase precursors to SOA and ozone. The Comprehensive Thermal Desorption Aerosol Gas Chromatograph (cTAG) is the first single instrument simultaneously quantitative for a broad range of compound-specific VOCs, IVOCs and SVOCs. cTAG is a two-channel instrument which measures concentrations of C5–C16 alkane-equivalent-volatility VOCs and IVOCs on one channel and C14–C32 SVOCs on the other coupled to a single high-resolution time-of-flight mass spectrometer, achieving consistent quantification across 15 orders of magnitude of vapor pressure. cTAG obtains concentrations hourly and gas–particle partitioning for SVOCs every other hour, enabling observation of the evolution of these species through oxidation and partitioning into the particle phase. Online derivatization for the SVOC channel enables detection of more polar and oxidized species. In this work we present design details and data evaluating key parameters of instrument performance such as I/VOC collector design optimization, linearity and reproducibility of calibration curves obtained using a custom liquid evaporation system for I/VOCs and the effect of an ozone removal filter on instrument performance. Example timelines of precursors with secondary products are shown, and analysis of a subset of compounds detectable by cTAG demonstrates some of the analytical possibilities with this instrument.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Development of observation-based global multilayer soil moisture products for 1970 to 2016

Abstract. Soil moisture (SM) datasets are critical to understanding the global water, energy, and biogeochemical cycles and benefit extensive societal applications. However, individual sources of SM data (e.g., in situ and satellite observations, reanalysis, offline land surface model simulations, Earth system model – ESM – simulations) have source-specific limitations and biases related to the spatiotemporal continuity, resolutions, and modeling and retrieval assumptions. Here, we developed seven global, gap-free, long-term (1970–2016), multilayer (0–10, 10–30, 30–50, and 50–100 cm) SM products at monthly 0.5∘ resolution (available at https://doi.org/10.6084/m9.figshare.13661312.v1; Wang and Mao, 2021) by synthesizing a wide range of SM datasets using three statistical methods (unweighted averaging, optimal linear combination, and emergent constraint). The merged products outperformed their source datasets when evaluated with in situ observations (mean bias from −0.044 to 0.033 m3 m−3, root mean square errors from 0.076 to 0.104 m3 m−3, Pearson correlations from 0.35 to 0.67) and multiple gridded datasets that did not enter merging because of insufficient spatial, temporal, or soil layer coverage. Three of the new SM products, which were produced by applying any of the three merging methods to the source datasets excluding the ESMs, had lower bias and root mean square errors and higher correlations than the ESM-dependent merged products. The ESM-independent products also showed a better ability to capture historical large-scale drought events than the ESM-dependent products. The merged products generally showed reasonable temporal homogeneity and physically plausible global sensitivities to observed meteorological factors, except that the ESM-dependent products underestimated the low-frequency temporal variability in SM and overestimated the high-frequency variability for the 50–100 cm depth. Based on these evaluation results, the three ESM-independent products were finally recommended for future applications because of their better performances than the ESM-dependent ones. Despite uncertainties in the raw SM datasets and fusion methods, these hybrid products create added value over existing SM datasets because of the performance improvement and harmonized spatial, temporal, and vertical coverages, and they provide a new foundation for scientific investigation and resource management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Robust historical evapotranspiration trends across climate regimes

Evapotranspiration (ET) links the hydrological, energy and carbon cycles on the land surface. Quantifying ET and its spatio-temporal changes is also key to understanding climate extremes such as droughts, heatwaves and flooding. Regional ET estimates require reliable observation-based gridded ET datasets, and while many have been developed using physically based, empirically based and hybrid techniques, their efficacy, and particularly the efficacy of their uncertainty estimates, is difficult to verify. In this work, we extend the methodology used in Hobeichi et al. (2018) to derive two new versions of the Derived Optimal Linear Combination Evapotranspiration (DOLCE) product, with observationally constrained spatio-temporally varying uncertainty estimates, higher spatial resolution, more constituent products and extended temporal coverage (1980–2018). After demonstrating the efficacy of these uncertainty estimates with out-of-sample testing, we derive novel ET climatology clusters for the land surface, based on the magnitude and variability of ET at each location on land. The new clusters include three wet and three dry regimes and provide an approximation of Köppen–Geiger climate classes. The verified uncertainty estimates and extended time period then allow us to examine the robustness of historical trends spatially and in each of these six ET climatology clusters. We find that despite robust decreasing ET trends in some regions these do not correlate with behavioural ET clusters. Each cluster, and the majority of the Earth's surface, shows clear robust increases in ET over the recent historical period. The new datasets DOLCE V2.1 and DOLCE V3 can be used for benchmarking global ET estimates and for examining ET trends respectively.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Resilient Operation of Power Distribution Systems Using MPC-Based Critical Service Restoration: Preprint

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2021

Resilient Operation of Power Distribution Systems Using MPC-Based Critical Service Restoration

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2020

Identification and Quantification of Beta-Emitters in Waste Packages by Segmented Gamma-Scanning - 20114

A commonly applied method for the non-destructive characterization of the radioactive inventory of waste packages (e.g. 200-l drums) is segmented gamma-scanning. The gamma-emitting inventory is identified via their characteristic lines in the measured spectra. The corresponding net peak areas allow for quantification applying appropriate models. The identification and quantification of beta-emitters present in the waste package require further evaluation. While some of them emit gamma-rays, too, which can be used for identification and quantification, pure beta-emitters contribute to the gamma-spectra only by their Bremsstrahlung-spectra. {sup 32}P, {sup 90}Sr and {sup 90}Y are typical representatives of these isotopes. A method is developed for identification and quantification of these pure beta-emitters in radioactive waste packages. Origin of the procedure is the evaluation of the measured gamma-spectrum. After the identification of all lines present in the spectrum (characteristic gamma-lines, single- and double-escape-lines, X-ray-lines, annihilation peak etc.) this information is used together with calibration data of the detection system and further a-priori information like matrix composition, cross-section data from databases etc. to simulate the spectra without the Bremsstrahlung-components by an appropriate method. This may be Monte Carlo simulations or, like in our case, a computer code taking into account mathematical descriptions of the physical effects like photo effect, (multiple) Compton scattering, backscattering etc. and the attenuation effects of the waste matrix. In any case, the resulting spectrum must simulate the real spectrum without its Bremsstrahlung-component as good as possible. Subtracting the simulated from the measured spectrum then results in the Bremsstrahlung-components caused by the beta-emitters. Prior to evaluation, a set of Bremsstrahlung-spectra is calculated and stored for different beta-emitters in different matrices using the Bethe-Heitler equation. This data set can be extended for new conditions, if necessary, and is allocated in a data base for future application. Applying a non-linear optimization using the generalized reduced gradient method for comparing the calculated with the derived Bremsstrahlung-spectra results in the identification of the beta-emitter and a factor, which is proportional to its activity. The proportionality is determined by previous calibration measurements. The application of the method is demonstrated by measurements on a 200-l mock-up drum and a real waste package. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2022

Approximating a linear multiplicative objective in watershed management optimization

Implementing management practices in a cost-efficient manner is critical for regional efforts to reduce the amount of pollutants entering the Chesapeake Bay. We study the problem of selecting a subset of practices that minimizes pollutant load—subject to budgetary and environmental constraints—as simulated in a widely used regulatory watershed model. Mimicking the computation of pollutant load in the regulatory model, we formulate this problem as a continuous optimization model with a linear multiplicative objective function and linear constraints. To lay the groundwork for incorporating additional stakeholder requirements in the future, especially those that would require integer variables, we present and study a continuous linear optimization model that approximates the nonlinear model. The linear model, which requires an exponential number of variables, arises naturally as an alternative model for the same underlying physical process. We examine the theoretical behavior of these optimization models and investigate restrictions of the linear model to handle its large number of variables. Through extensive computational tests on real and randomly generated instances, we demonstrate that the linear model and its restrictions provide optimal solutions close to those of the nonlinear model in practice, despite poor approximation properties in the worst case. We conclude that the linear model—together with our approach to handling its large number of variables—provides a viable framework from which to extend the optimization model to better meet the needs of the Chesapeake Bay watershed management stakeholders.

54 ENVIRONMENTAL SCIENCES↗