SEARCH · Engineering Papers

Results for “constrained 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 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Analysis and Integration of the Hydraulic Fracturing Test Site-2 (HFTS-2) Comprehensive Dataset

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Permian (Delaware) Basin. The unique aspect of this program was the acquisition of a unique, comprehensive, diagnostic dataset. Additionally, shorter parent wells drilled three years before the child wells offered clear distinction between the stages influenced by parent-child effects and the stages without any effects. The goal of this study was to analyze and integrate this comprehensive diagnostic dataset to understand the areal and vertical extent of hydraulic fractures (HF). The paper also provides insights on the effects of parent wells’ depletion on child well HF geometry based on various monitoring methods and subsurface models. Areal and vertical coverage for all HFTS-2 wells during stimulation and depletion was estimated based on analysis and interpretation of diagnostics and advanced modeling results. HFTS-2 diagnostics included microseismic (MS), pre- and post-stimulation logs and cores, bottomhole gauges, and fiber optic (FO) data. The diagnostics results (MS, FO, image logs) were integrated and used to calibrate subsurface models. Additional field tests were designed and implemented for depletion monitoring. The tailored program for monitoring depletion included vertical and slant well pressures, interference testing, and a vertical strain depletion trial. Areal Coverage: Conventional MS (and FO MS) were used to compute HF dimensions, which were compared with diagnostics (FO strain, gauge, image logs) observations and calibrated subsurface models. A post-production interference test did not show offset well communication. Vertical Coverage: Vertical coverage during stimulation was monitored using a vertical monitoring well. The stronger mechanical strain signals showed good correlation with MS event intensities, geomechanical properties, and gauge inferences. Vertical depletion was estimated based on vertical/slant well gauges and strain depletion tests. Parent-Child Effects: Diagnostics and calibrated subsurface models show asymmetry in child well fracture geometries for stages that overlap parent wells. Child well image logs serve as a good indicator for parent Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/2-21URTC/D021S031R004/2477423/urtec-2021-5241-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5241 well HF tracking. Child well MS events had an eastward bias, in line with pre-stimulation image logs, and was confirmed by parent well frac hits. Novel/Additive Information: The dataset presents a unique, over-constrained problem space to compare independent techniques to arrive at HF metrics (i.e., stimulation height and/or half-length), unlike a single source dataset, in which calibration is done using available data to guide predictions. Here, the asymmetry in HF geometry seen in the stages influenced by parent-child effects offers unique insights into well spacing and landing, which are key capital decisions the unconventional resources industry is seeking to optimize.

58 GEOSCIENCES↗

DOE OSTI2023

Understanding the Structure and Dynamics of Conjugated Polymers by Advancing Deuteration Chemistry and Neutron Scattering (Final Report)

The overarching goal of the proposed work was to set up a partnership between the University of Southern Mississippi (USM) and Oak Ridge National Laboratory (ORNL) to develop novel approaches to measure the backbone rigidity of conjugated polymers (CPs) and understand the critical role of sidechains on the backbone conformation and the materials macroscopic property. The backbone rigidity greatly influences the electronic properties of CPs, which ultimately determines the functionality and performance of these materials. Improvements in the electronic properties of CPs would allow for enhanced charge transport in semiconductor devices, improved photovoltaic performance, recycling of waste heat in thermoelectrics, and discovery of new phenomena that will enable the next generation of energy technologies. Although significant progress has been made to optimize the optical and electronic properties of CPs, largely through Edisonian methodologies, it remains a challenge to experimentally characterize conjugated backbone conformation (chain rigidity, torsion, planarity, and short-range order) and relate these to the fundamental optical and electronic properties (electronic coupling, charge transport, etc.). This has left fundamental gaps in our knowledge of the most basic structure/property relationships within these systems, precluded the study of fundamental physical phenomena, and constrained the design and realization of new electronic and device functionalities. Thus, the major goal of this work is to use novel deuteration methodologies via systematic synthetic approaches, and neutron scattering techniques to comprehensively characterize the structural and dynamic properties of CPs in contrast-matching solvents. Our work would, for the first time, elucidate the relationship between backbone rigidity and macroscopic properties. They will also allow a rational formulation of design principles for next-generation CPs that are resilient to disorder through precise control of the delocalized electrons along the polymer backbone. Overall, this project will advance our understanding of the structure, dynamics, and fundamental physics of these materials, which is crucial for enabling the prediction, design, control, and manipulation of current and emerging material electronic properties.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Optimal Control of Differentially Private EV Charging: A Scalable Learning Approach Under Uncertainty

Internet of Things (IoT)-enabled electric vehicles (IoEVs) enable intelligent charging coordination that accounts for grid congestion. However, increased data exchange raises privacy concerns, as charging patterns can reveal sensitive driver behavior to grid operators. Here, we propose a differentially private (DP) EV charging framework that enables coordinated control while protecting driver data with theoretical privacy guarantees. Nevertheless, integrating DP inevitably introduces uncertainty into the control strategy for EVs, which can lead to infeasible solutions. To tackle this challenge, we develop a feasible and scalable control algorithm based on constrained reinforcement learning (CRL) and convex hulls. While our framework is designed to handle the uncertainty introduced by DP, it is general and also applicable to other sources of uncertainty in EV charging, such as the stochastic nature of driver behavior and renewable variability. This ensures feasible and privacy-preserving coordination of EV charging at scale. Our method constructs convex hulls within the action space to guarantee feasibility under stochastic constraints and incorporates constraint reduction techniques to improve scalability. Case studies based on IEEE benchmark systems demonstrate that the proposed approach effectively balances feasibility under uncertainty, scalability, and privacy in large-scale EV charging control.

Engineering - Power transmission and distribution↗

DOE OSTI2021

Optimization-based algorithms for nonlinear mechanics and frictional contact

An optimization-based strategy for solving nonlinear mechanics problems is proposed. In contrast to typical nonlinear equation solver algorithms that aim to find zeros in the residual force function, we minimize an energy (or energy-like) function to encourage solutions which are locally stable equilibria. These smooth and potentially non-convex objective functions are minimized using a preconditioned conjugate-gradient trust-region algorithm. Contact is formulated as an inequality constrained minimization problem, and is solved with an augmented Lagrangian algorithm. Friction is included in the approach via a regularized quasi-potential energy, and other dissipative behavior is included through the use of variational constitutive updates. Finally, to accelerate convergence rates for the Lagrange multipliers, we propose a novel multiplier update algorithm utilizing the Fischer-Burmeister function, and demonstrate super-linear solver convergence for some applications.

42 ENGINEERING↗

DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Characterization and Optimization of Skipper CCDs for the SOAR Integral Field Spectrograph

We present results from the characterization and optimization of Skipper charge-coupled devices (CCDs) for use in a focal plane prototype for the Southern Astrophysical Research Integral Field Spectrograph (SIFS). We tested eight Skipper CCDs and selected six for SIFS based on performance results. The Skipper CCDs are 6k × 1k, 15 μm pixels, thick, fully depleted, p-channel devices that have been thinned to ~250 μm, backside processed, and treated with an anti-reflective coating. We demonstrate a single-sample readout noise of <4.3 e - rms pixel -1 in all amplifiers. We optimize the readout sequence timing to achieve a readout noise of 0.5 e - rms pixel -1 after 74 non-destructive measurements, which can be accomplished in a region covering 5% of the detector area in a readout time of <4 minutes. We demonstrate single-photon-counting in all 24 amplifiers (four amplifiers per detector) with a readnoise of σ N ~ 0.18 e - rms pixel -1 after N samp = 400 samples, and we constrain the degree of nonlinearity to be ≲1% at low signal levels (0 e - to 50 e - ). Clock-induced charge (CIC) remains an important issue when the Skipper CCD is configured to provide a large full-well capacity. We achieve a CIC rate of <1.45 × 10 -3 e - pixel -1 frame -1 for a full-well capacity of ~900 e - , which increases to a CIC rate of ~3 e - pixel -1 frame -1 for full-well capacities ~40,000-65,000 e - . We also perform conventional CCD characterization measurements such as charge transfer inefficiency (3.44 × 10 -7 on average), dark current (~2 × 10 -4 e - pixel -1 s -1 ), photon transfer curves, cosmetic defects (<0.45% “bad” pixels), and charge diffusion (point-spread function < 7.5 μm) to verify that these properties are consistent with expectations from conventional p-channel CCDs used for astronomy. Furthermore, we provide the first measurements of the brighter-fatter effect and absolute quantum efficiency (≳80% between 450 and 980 nm; ≳90% between 600 and 900 nm) using Skipper CCDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Search for the nonresonant production of Higgs boson pairs via gluon fusion and vector-boson fusion in the $b\bar{b}$τ + τ – final state in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for the nonresonant production of Higgs boson pairs in the HH → $b\overline{b}$$\tau^{+}\tau^{-}$ channel is performed using 140 fb -1 of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. The analysis strategy is optimized to probe anomalous values of the Higgs boson self-coupling modifier $\kappa_{\tau}$ and of the quartic HHVV (v = W,Z) coupling modifier $\kappa$ 2⁢V . No significant excess above the expected background from Standard Model processes is observed. An observed (expected) upper limit μ HH < 5.9⁢(3.3) is set at 95% confidence-level on the Higgs boson pair production cross section normalized to its Standard Model prediction. The coupling modifiers are constrained to an observed (expected) 95% confidence interval of -3.1 <$\kappa_{\tau}$ < 9.0 (-2.5 < $\kappa_{\tau}$ < 9.3) and -0.5 <$\kappa$ 2⁢V < 2.7 (-0.2 <$\kappa$ 2⁢V < 2.4), assuming all other Higgs boson couplings are fixed to the Standard Model prediction. The results are also interpreted in the context of effective field theories via constraints on anomalous Higgs boson couplings and Higgs boson pair production cross sections assuming different kinematic benchmark scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Fast Adaptive Neural Control of Resonant Extraction at Fermilab

We present the development of a machine learning (ML) based regulation system for third-order resonant beam extraction in the Mu2e experiment at Fermilab. Classical and ML-based controllers have been optimized using semi-analytic simulations and evaluated in terms of regulation performance and training efficiency. We compare several controller architectures and discuss the integration of neural control into an adaptive framework. We also present progress on surrogate models that predict the controller response given a spill intensity and controller action history. To enable real-time deployment, we report progress on implementing low-latency, edge-based inference suitable for hardware-constrained environments. Our results demonstrate the feasibility and advantages of ML-based control in managing complex, time-varying physical systems, with broader implications for accelerator operations and other domains requiring fast, adaptive regulation.

Berlioz, Jose Rene [Fermilab]↗

DOE OSTI2024

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^{2}$ over 5 yr to constrain the cosmic expansion history through precise measurements of baryon acoustic oscillations (BAO). The scientific program for DESI was evaluated during a 5 month survey validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar Milky Way Survey (MWS), Bright Galaxy Survey (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the 5 yr program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a One-Percent Survey conducted at the conclusion of SV covering 140 deg$^{2}$ using the final target selection algorithms with exposures of a depth typical of the main survey. The SV indicates that DESI will be able to complete the full 14,000 deg$^{2}$ program with spectroscopically confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval z < 1.1, 0.39% over the redshift interval 1.1 < z < 1.9, and 0.46% over the redshift interval 1.9 < z < 3.5.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Using a surrogate-assisted Bayesian framework to calibrate the runoff-generation scheme in the Energy Exascale Earth System Model (E3SM) v1

Abstract. Runoff is a critical component of the terrestrial water cycle, and Earth system models (ESMs) are essential tools to study its spatiotemporal variability. Runoff schemes in ESMs typically include many parameters so that model calibration is necessary to improve the accuracy of simulated runoff. However, runoff calibration at a global scale is challenging because of the high computational cost and the lack of reliable observational datasets. In this study, we calibrated 11 runoff relevant parameters in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) using a surrogate-assisted Bayesian framework. First, the polynomial chaos expansion machinery with Bayesian compressed sensing is used to construct computationally inexpensive surrogate models for ELM-simulated runoff at 0.5∘ × 0.5∘ for 1991–2010. The error metric between the ELM simulations and the benchmark data is selected to construct the surrogates, which facilitates efficient calibration and avoids the more conventional, but challenging, construction of high-dimensional surrogates for the ELM simulated runoff. Second, the Sobol' index sensitivity analysis is performed using the surrogate models to identify the most sensitive parameters, and our results show that, in most regions, ELM-simulated runoff is strongly sensitive to 3 of the 11 uncertain parameters. Third, a Bayesian method is used to infer the optimal values of the most sensitive parameters using an observation-based global runoff dataset as the benchmark. Our results show that model performance is significantly improved with the inferred parameter values. Although the parametric uncertainty of simulated runoff is reduced after the parameter inference, it remains comparable to the multimodel ensemble uncertainty represented by the global hydrological models in ISMIP2a. Additionally, the annual global runoff trend during the simulation period is not well constrained by the inferred parameter values, suggesting the importance of including parametric uncertainty in future runoff projections.

58 GEOSCIENCES↗

DOE OSTI2026

Applications of Multibody for Everybody (M4E) in Marine Energy

Multibody for Everybody (M4E) is an open-source symbolic dynamics modeling framework designed to automate the derivation of equations of motion and simulation of constrained multibody systems using the joint coordinate formulation.

16 TIDAL AND WAVE POWER↗

DOE OSTI2026

Quantum Routing and Entanglement Dynamics Through Bottlenecks

To implement arbitrary quantum circuits in architectures with restricted interactions, one may effectively simulate all-to-all connectivity by routing quantum information. We consider the entanglement dynamics and routing between two regions only connected through an intermediate “bottleneck” region with few qubits. In such systems, where the entanglement rate is restricted by a vertex boundary rather than an edge boundary of the underlying interaction graph, existing results such as the small incremental entangling theorem give only a trivial constant lower bound on the routing time (the minimum time to perform an arbitrary permutation). We significantly improve the lower bound on the routing time in systems with a vertex bottleneck. Specifically, for any system with two regions 𝐿,𝑅 with 𝑁 𝐿 ,𝑁 𝑅 qubits, respectively, coupled only through an intermediate region 𝐶 with 𝑁 𝐶 qubits, for any 𝛿 > 0 we show a lower bound of Ω⁢(𝑁$^{1−𝛿}_{𝑅}$/√𝑁 𝐿⁢ 𝑁 𝐶 ) on the Hamiltonian quantum routing time when using piecewise time-independent Hamiltonians, or time-dependent Hamiltonians subject to a smoothness condition. We also prove an upper bound on the average amount of bipartite entanglement between 𝐿 and 𝐶,𝑅 that can be generated in time 𝑡 by such architecture-respecting Hamiltonians in systems constrained by vertex bottlenecks, improving the scaling in the system size from 𝑂⁡(𝑁 𝐿⁢ 𝑡) to 𝑂⁡(√𝑁 𝐿⁢ 𝑡). As a special case, when applied to the star graph (i.e., one vertex connected to 𝑁 leaves), we obtain an Ω⁡(√𝑁 1−𝛿 ) lower bound on the routing time and on the time to prepare 𝑁/2 Bell pairs between the vertices. We also show that, in systems of free particles, we can route optimally on the star graph in time Θ⁡(√𝑁) using Hamiltonian quantum routing, obtaining a speedup over gate-based routing, which takes time Θ⁡(𝑁).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Data-Driven Compression of Electron-Phonon Interactions

First-principles calculations of electron interactions in materials have seen rapid progress in recent years, with electron-phonon ( e − ph ) interactions being a prime example. However, these techniques use large matrices encoding the interactions on dense momentum grids, which reduces computational efficiency and obscures interpretability. For e − ph interactions, existing interpolation techniques leverage locality in real space, but the high dimensionality of the data remains a bottleneck to balance cost and accuracy. Here we show an efficient way to compress e − ph interactions based on singular value decomposition (SVD), a widely used matrix and image compression technique. Leveraging (un)constrained SVD methods, we accurately predict material properties related to e − ph interactions—including charge mobility, spin relaxation times, band renormalization, and superconducting critical temperature—while using only a small fraction (1%–2%) of the interaction data. These findings unveil the hidden low-dimensional nature of e − ph interactions. Furthermore, they accelerate state-of-the-art first-principles e − ph calculations by about 2 orders of magnitude without sacrificing accuracy. Our Pareto-optimal parametrization of e − ph interactions can be readily generalized to electron-electron and electron-defect interactions, as well as to other couplings, advancing quantitative studies of condensed matter. Published by the American Physical Society 2024

Physics↗

DOE OSTI2025

Probing the Higgs boson CP properties in vector-boson fusion production in the H → τ + τ − channel with the ATLAS detector

The CP properties of the Higgs boson are studied in the vector-boson fusion production mode. The analysis exploits the decay mode of the Higgs boson into two τ-leptons using 140 fb −1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Results are obtained using the Optimal Observable method. CP-violating interactions between the Higgs boson and electroweak gauge bosons are considered in the effective field theory framework, with the interaction strength described in the HISZ basis by $\tilde{d}$, and in the Warsaw basis by c H$\tilde{W}$ , c H$\tilde{B}$ , and c H$\tilde{W}$B . No deviations relative to the Standard Model are observed, and limits are obtained on the strength parameters. The $\tilde{d}$ parameter is constrained to the interval [−0.012, 0.044] at the 95% confidence level while c H$\tilde{W}$ is constrained to [−0.24, 0.83], when considering both linear and quadratic effects of physics beyond the Standard Model.

Hadron-Hadron Scattering↗

DOE OSTI2025

Optimizing lossy state preparation for quantum sensing using Hamiltonian engineering

One of the most prominent platforms for demonstrating quantum sensing below the standard quantum limit is the spinor Bose–Einstein condensate. While a quantum advantage using several tens of thousands of atoms has been demonstrated in this platform, it faces an important challenge: atom loss. Atom loss is a Markovian error process modeled by Lindblad jump operators, and a no-go theorem, which we also show here, states that the loss of atoms in all spin components reduces the quantum advantage to a constant factor. Here, we show that this no-go theorem can be circumvented if we constrain atom losses to a single spin component. Moreover, we show that in this case, the maximum quantum Fisher information with N atoms scales as N 3/2 , establishing that a scalable quantum advantage can be achieved despite atom loss. Although Lindblad jump operators are generally non-Hermitian and non-invertible, we use their Moore–Penrose inverse to develop a framework for constructing several states with this scaling of Fisher information in the presence of losses. We use Hamiltonian engineering with realistic Hamiltonians to develop experimental protocols for preparing these states. Finally, we discuss possible experimental techniques to constrain the losses to a single spin mode.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Predicting multiphase flow and tracer transport for an underground chemical explosive test

Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Data-centric machine learning in quantum information science

Abstract We propose a series of data-centric heuristics for improving the performance of machine learning systems when applied to problems in quantum information science. In particular, we consider how systematic engineering of training sets can significantly enhance the accuracy of pre-trained neural networks used for quantum state reconstruction without altering the underlying architecture. We find that it is not always optimal to engineer training sets to exactly match the expected distribution of a target scenario, and instead, performance can be further improved by biasing the training set to be slightly more mixed than the target. This is due to the heterogeneity in the number of free variables required to describe states of different purity, and as a result, overall accuracy of the network improves when training sets of a fixed size focus on states with the least constrained free variables. For further clarity, we also include a ‘toy model’ demonstration of how spurious correlations can inadvertently enter synthetic data sets used for training, how the performance of systems trained with these correlations can degrade dramatically, and how the inclusion of even relatively few counterexamples can effectively remedy such problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Data-Driven Unit Commitment Refinement - a Scalable Approach for Complex Modern Power Grids

Integration of renewable generation, which is often intermittent and decentralized, substantially increases the stochasticity and complexity of power grid operations. Future power systems planning will require significant computational capability to evaluate balance between demand and supply under varying conditions, both temporally and spatially. The standard approach for generation unit commitment is to use mixed-integer linear programming to find the optimal generation schedule considering ramping and generator constraints. In the future grid this poses computational scalability challenges because generation and demand are not known with certainty due to stochasticity in weather and complexity of the grid. To address this challenge, we present a data-driven unit commitment approach that can efficiently include stochastic weather impacts and contingency considerations to improve unit commitment. Our approach uses graph-based data analytics techniques on solutions to the security constrained (and possibly stochastic) economic dispatch problem to identify potential improvements to a given unit commitment. Recent breakthroughs in fully-parallel stochastic economic dispatch software allow this approach to be scalably deployed. Simulations on synthetic South Carolina and Texas grids show this method can improve grid reliability with security constraints over a set of contingencies, while also meaningfully lowering total generation cost.

Holt, Timothy↗