Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 469 records · Page 26

System Analysis of an Internal Combustion Engine (ICE) - Solid Oxide Fuel Cell (SOFC) Hybrid Cycle

The variability of renewable energy sources poses challenges for reliable grid operation. Conventional thermal power sources, though reliable, often lack operational flexibility. Hybrid energy systems that integrate Solid Oxide Fuel cells (SOFC) with Internal Combustion Engines (ICE) offer a promising solution by achieving relatively higher efficiency and grid-following capability. This study investigates performance of a 100-kW pressurized SOFC-ICE hybrid cycle. In this configuration, unutilized SOFC fuel is used to drive the engine, with a turbocharger providing air supply and an external reformer generating syngas. System components were numerically modeled using MATLAB/SIMULINK for the SOFC and reformer, and EBSILON® for the ICE and balance of plant. Parametric studies varied fuel utilization (70-90%), reformer temperature (600 – 1000K), anode off-gas recirculation (0 – 70%), and current density (0.2 – 0.55 A/cm2). Results show that the SOFC and ICE operate as thermally independent topping and bottoming cycles, achieving peak efficiency of 62% under optimized conditions.

hybrid↗

Predicting Open Quantum Dynamics with Data-Informed Quantum-Classical Dynamics

We introduce a data-informed quantum-classical dynamics (DIQCD) approach for predicting the evolution of an open quantum system. The equation of motion in DIQCD is a Lindblad equation with a flexible, time-dependent Hamiltonian that can be optimized to fit sparse and noisy data from local observations of an extensive open quantum system. We demonstrate the accuracy and efficiency of DIQCD for both experimental and simulated quantum devices. We show that DIQCD can predict entanglement dynamics of ultracold molecules (calcium fluoride) in optical tweezer arrays. DIQCD also successfully predicts carrier mobility in organic semiconductors (rubrene) with accuracy comparable to nearly exact numerical methods.

Lindblad equation↗

Evaluation of the Energy, Hygrothermal, and Thermal Capacity Performance of Cross-Laminated Timber

Cross-laminated timber (CLT) construction is gaining momentum in the US because it offers multiple advantages over traditional construction methods. Benefits that have received the most attention focus on constructability, the environment, and protection (e.g., blast resistance), although CLT construction is likely to offer other benefits, as well. Still, these have not been studied at length because such evaluations are costly, requiring long-term assessments in an actual building and specialized technical knowledge. Among the possible benefits, CLT construction likely provides a higher-performing building envelope. Using CLT panels to enclose a building means fewer joints in the opaque envelope than what is required in traditional stick-framed construction. Fewer joints mean fewer locations where the air- and water-resistive barrier (WRB) could be compromised; thus, a CLT building enclosure may require less maintenance and have a longer lifespan than a traditionally built structure because of fewer air and water leaks. In addition, CLT’s thermal mass moderates indoor temperatures, allowing the heating, ventilation, and air conditioning (HVAC) system to operate more efficiently during peak hours, reducing operational energy consumption throughout the lifetime of the CLT building (Salonvaara et al., 2022). Furthermore, more stable indoor temperatures can increase occupant comfort. The CLT’s thermal mass can also reduce energy costs by adjusting to utility time-of-use pricing without affecting occupant comfort. The ability of CLT buildings to bridge periods without HVAC operation prepares them for future grid interaction and provides a certain level of resilience against power outages. Researchers have attempted to quantify these benefits; however, their work is based on simplified simulations with numerous assumptions. To correctly understand the benefits, an actual building must be monitored. Therefore, information needs to be gathered on indoor and outdoor temperatures, HVAC energy consumption, thermostat setpoints, temperatures, and thermal transport in CLT components to comprehend how these parameters are affected by the CLT’s thermal mass. These data are needed to reduce the number of assumptions and calibrate simulation models to optimize HVAC controls to minimize overall energy consumption, reduce energy use and higher fees during peak demand, and maintain occupant comfort. Additionally, the calibrated simulation model allows the optimization exercise to be repeated in various US climates. Potential benefits can be tailored to buildings in various locations, and decisions can be made on where CLT construction could be most advantageous. Furthermore, monitoring and simulation results are needed to evaluate the durability of the CLT structures in different climates. This project’s researchers gathered information to help understand and quantify the benefits of CLT buildings concerning operational energy, moderated indoor temperatures, and comfort; the dynamic operation to provide grid services; and resilience in times of power outage. Through the corroboration of simulation models with real-world measurements, this study paves the way for extrapolating findings to other climatic zones and building typologies, thereby broadening the understanding of CLT’s multifaceted benefits and reinforcing its position as a material of choice in sustainable construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Multi-parametric analysis for mixed integer linear programming: An application to transmission upgrade and congestion management

Upgrading the capacity of existing transmission lines is essential for meeting the growing energy demands, facilitating the integration of renewable energy, and ensuring the security of the transmission system. This study focuses on the selection of lines whose capacities and by how much should be expanded from the perspective of the Independent System Operators (ISOs) to minimize the total system cost. We employ advanced multi-parametric programming and an enhanced branch-and-bound algorithm to address complex mixed-integer linear programming (MILP) problems, considering multi-period time constraints and physical limitations of generators and transmission lines. To characterize the various decisions in transmission expansion, we model the increased capacity of existing lines as parameters within a specified range. This study first relaxes the binary variables to continuous variables and applies the Lagrange method and Karush-Kuhn-Tucker (KKT) conditions to obtain optimal solutions and identify critical regions associated with active and inactive constraints. Moreover, we extend the traditional branch-and-bound (B&B) method by determining the problem’s upper and lower bounds at each node of the B&B decision tree, helping to manage computational challenges in large-scale MILP problems. Here, we compare the difference between the upper and lower bounds to obtain an approximate optimal solution within the decision-makers’ tolerable error range. In addition, the first derivative of the objective function on the parameters of each line is used to inform the selection of lines for easing congestion and maximizing social welfare. Finally, the capacity upgrades are selected by weighing the reductions in system costs against the expense of upgrading line capacities. The findings are supported by numerical simulations and provide transmission-line planners with decision-making guidance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Studies of the outer-off-midplane lower hybrid wave launch scenario for plasma start-up on the TST-2 spherical tokamak

Abstract Establishment of an efficient central solenoid (CS) free tokamak plasma start-up method may lead to an economical fusion reactor. CS-free start-up using lower hybrid (LH) waves has been studied on the TST-2 spherical tokamak. Plasma current of about a quarter of CS-driven discharges has been obtained fully non-inductively using the outer-midplane and top LH launchers. Recently, an outer-off-midplane LH launcher was developed to achieve higher plasma current by optimizing for core absorption and minimal fast electron losses. Using the (outer-)off-midplane launcher, fully non-inductive plasma current start-up up to about 8 kA was achieved. Coupled ray-tracing and Fokker–Planck simulation was performed on equilibria reconstructed with an extended MHD model. It was found that the experimentally observed plasma current was in reasonable agreement with the numerical simulation. The simulation predicted appreciable orbit losses for the off-midplane launcher driven discharge at the present parameters, which was consistent with the experimentally observed x-ray radiation characteristics. The simulation showed that the current density was saturated for the present off-midplane launcher discharges and higher density and higher LH power was necessary to achieve higher plasma current.

Physics↗

Quantum Davidson algorithm for excited states

Abstract Excited state properties play a pivotal role in various chemical and physical phenomena, such as charge separation and light emission. However, the primary focus of most existing quantum algorithms has been the ground state, as seen in quantum phase estimation and the variational quantum eigensolver (VQE). Although VQE-type methods have been extended to explore excited states, these methods grapple with optimization challenges. In contrast, the quantum Krylov subspace (QKS) method has been introduced to address both ground and excited states, positioning itself as a cost-effective alternative to quantum phase estimation. However, conventional QKS methodologies depend on a pre-generated subspace through real or imaginary-time evolutions. This subspace is inherently expansive and can be plagued with issues like slow convergence or numerical instabilities, often leading to relatively deep circuits. Our research presents an economic QKS algorithm, which we term the quantum Davidson (QDavidson) algorithm. This innovation hinges on the iterative expansion of the Krylov subspace and the incorporation of a pre-conditioner within the Davidson framework. By using the residues of eigenstates to expand the Krylov subspace, we manage to formulate a compact subspace that aligns closely with the exact solutions. This iterative subspace expansion paves the way for a more rapid convergence in comparison to other QKS techniques, such as the quantum Lanczos. Using quantum simulators, we employ the novel QDavidson algorithm to delve into the excited state properties of various systems, spanning from the Heisenberg spin model to real molecules. Compared to the existing QKS methods, the QDavidson algorithm not only converges swiftly but also demands a significantly shallower circuit. This efficiency establishes the QDavidson method as a pragmatic tool for elucidating both ground and excited state properties on quantum computing platforms.

97 MATHEMATICS AND COMPUTING↗

Collective neutrino oscillations in three flavors on qubit and qutrit processors

Collective neutrino flavor oscillations are of primary importance in understanding the dynamic evolution of core-collapse supernovae and subsequent terrestrial detection, but also among the most challenging aspects of numerical simulations. This situation is complicated by the quantum many-body nature of the problem due to neutrino-neutrino interactions, which demands a quantum treatment. An additional complication is the presence of three flavors, which often is approximated by the electron flavor and a heavy lepton flavor. In this work, we provide both qubit and qutrit encodings for all three flavors, and develop optimized quantum circuits for the time evolution and analyze the Trotter error. We conclude our study with a hardware experiment of a system of two neutrinos with superconducting hardware: the IBM Torino device for qubits and Advanced Quantum Testbed device at the Lawrence Berkeley National Laboratory for qutrits. We find that error mitigation greatly helps in obtaining a signal consistent with simulations. Finally, while hardware results are comparable at this stage, we expect the qutrit setup to be more convenient for large-scale simulations since it does not suffer from probability leakage into nonphysical qubit space, unlike the qubit setup.

Neutrino oscillations↗

Single-Crystal Diffuse Neutron Scattering Study of the Dipole-Octupole Quantum Spin-Ice Candidate Ce 2⁢ Zr 2 ⁢O 7 : No Apparent Octupolar Correlations Above 𝑇 = 0.05 K

The insulating magnetic pyrochlore Ce 2 ⁢Zr 2 ⁢O 7 has gained attention as a quantum spin-ice candidate with dipole-octupole character that arises from the crystal-electric-field ground-state doublet for the Ce 3+ Kramers ion. This dipole-octupole character permits both spin-ice phases based on magnetic dipoles and those based on more-exotic octupoles. This work reports low-temperature neutron diffraction measurements on single-crystal Ce 2 ⁢Zr 2⁢ O 7 with 𝑄 coverage both at low 𝑄, where the magnetic form factor for dipoles is near maximal, and at high 𝑄, covering the region where the magnetic form factor for Ce 3+ octupoles is near maximal. This study was motivated by recent powder neutron diffraction studies of other Ce-based dipole-octupole pyrochlores, Ce 2 ⁢Sn 2 ⁢O 7 and Ce 2 ⁢Hf 2 ⁢O 7 , which each showed temperature-dependent diffuse diffraction at high 𝑄, interpreted as arising from octupolar correlations. Our measurements use an optimized single-crystal diffuse scattering instrument that allows us to screen against strong Bragg scattering from Ce 2 ⁢Zr 2 ⁢O 7 . The temperature-difference neutron diffraction reveals a low-𝑄 peak consistent with dipolar spin-ice correlations reported in previous work, and an alternation between positive and negative net intensity at higher 𝑄. These features are consistent with our numerical-linked-cluster calculations using pseudospin interaction parameters previously reported for Ce 2 ⁢Zr 2⁢ O 7 , Ce 2 ⁢Sn 2 ⁢O 7 , and Ce 2 ⁢Hf 2 ⁢O 7 . Importantly, neither the measured data nor any of the NLC calculations show evidence for increased scattering at high 𝑄 resulting from octupolar correlations. We conclude that at the lowest attainable temperature for our measurements (𝑇 = 0.05 K), scattering from octupolar correlations in Ce 2 ⁢Zr 2 ⁢O 7 is not present in the neutron diffraction signal on the level of our observation threshold of around 0.1% of the low-𝑄 dipole scattering. We compare these results to those obtained earlier on powder Ce 2 ⁢Sn 2 ⁢O 7 and Ce 2⁢ Hf 2⁢ O 7 , and to low-energy inelastic neutron scattering from single-crystal Ce 2 ⁢Zr 2 ⁢O 7 .

36 MATERIALS SCIENCE↗

Memory-Aware External Facelist Calculation: A Data-Parallel Atomic Hash Counting Approach

Unstructured volumetric meshes serve as fundamental data representations in various scientific simulations and analyses. They play a crucial role in representing complex computational domains and are essential for important numerical techniques, such as finite element analysis. Whenever such a mesh is read from a file, streamed in-situ, or generated by algorithms, scientific visualization libraries rely on calculating the external surface of a geometry, named “external facelist”, to produce a polygonal mesh for rendering. Consequently, external facelist calculation has become one of the most widely used algorithms in the scientific visualization domain, necessitating optimal performance. In this paper, we explore relevant work on external facelist calculation algorithms in two common visualization libraries, VTK and Viskores, assess their performance and memory constraints, and introduce a novel memory-aware external facelist calculation algorithm employing an atomic hash counting approach. This algorithm fully leverages Viskores' data-parallel primitive operations, facilitating its execution across diverse many-core architectures. Our algorithm features the lowest memory footprint on the GPU and the second-lowest on the CPU among all evaluated methods, and it also delivers the fastest performance on both CPU and GPU. It has been made available under an open-source license in the VTK and Viskores visualization systems.

Tsalikis, Spiros [Kitware] (ORCID:0000000151137195↗

Using "AI Poincare" to analyze non-linear integrable optics

This study dives into the applicability of using automated discovery of conserved quantities in dynamical systems relevant to accelerator physics. Specifically, we explore the performance of AI Poincaré in analyzing numerical trajectory data obtained using the McMillan system of non-linear integrable optics. A comprehensive evaluation of the algorithm's performance is conducted through diverse methodologies. These include the analysis of the estimated number of conserved quantities embedded in a dataset and the deviation of interpolated points on the inferred manifold with respect to points in actually in the dataset. the investigation identifies an optimal range of perturbation distances where the underlying manifold extraction algorithm inside AI Poincaré exhibits optimal performance. Additionally, an improved neural network architecture is proposed based on the observed results. Finally, we apply the algorithm to preliminary experimental data from the Integrable Optics Test Accelerator at Fermilab to successfully infer the number of conserved quantities even in the presence of fast decoherence of the measured signal.

Osmanov, Lazare [Free U. Tbilisi]↗

Optimal Operation of Residential High Performance Water Heater for Reduction of Electricity Cost and Peak Demand Through Field Validation

Water heating accounts for about 18% of a typical US home’s energy use. Modern water heaters have enabled control options through APIs, offering customers the opportunity to reduce their energy cost and peak demand by dynamically adjusting settings. A water heater’s capacity to store energy using its storage tank makes it an asset for peak demand reduction and energy cost savings. For this reason, a mixed-integer linear programming model is proposed to minimize the energy cost of a high-performance water heater while also reducing the peak demand of the residential household under a time-of-use utility rate by dynamically changing the water heater’s running mode. Specifically, a multi-objective optimization model is formulated to determine the mode settings of the water heater considering hot water use, time-of-use rate, and peak demand limit of the residential household. The mode settings are associated with different dead bands of water temperature for triggering on/off action of the heat pump and heating element. A 66-gal hybrid electric high performance water heater was used for numerical simulation and practical experiments. The simulation results were well aligned with measurements of practical experiments, validating the soundness of the thermodynamic model. In addition, reductions of energy cost, enabling affordability, and reducing peak demand are demonstrated. The research team also developed a software framework with dashboards to automatically and continuously monitor and manage devices.

Liu, Guodong [ORNL] (ORCID:0000000213498608)↗

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗

Methods and Tools To Assess Robustness of Nuclear Plant Outages

Refueling outages are one of the most challenging phases in a nuclear power plant (NPP) operating cycle. Refueling outages are extremely costly for an NPP due to the large amount of required resources and because of lost revenue due to plant being off the grid. Outage durations have steadily decreased across the industry over that last few decades primarily due to improved planning and coordination, but there are still many plants that struggle to meet the performance metrics accomplished by other utilities. Schedule resilience is one of the issues. NPP outages require scheduling thousands of activities within 30 days on average. Despite detailed planning, once the outage starts, numerous emergent issues typically appear along with schedule delays requiring continuous replanning and adjustment. When schedule disruption occurs during an outage, plant staff make urgent efforts to recover but are often not able to maintain the planned outage duration. These outage delays can cost a utility several million dollars per day. Tools that could help outage schedulers create a more resilient schedule and allow them to optimally reschedule emergent work could significantly reduce outage delays. One key aspect of creating a resilient schedule is to have accurate estimates for activity duration. Another important outage scheduling capability is the ability to schedule emergent work with minimal disruption. This paper focuses on developing tools and methods to support NPPs with outage schedule optimization and it describes the initial development of tools to support outage management that leverage computational and machine learning methods.

97 - MATHEMATICS AND COMPUTING↗

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗

A Model for Optimally Allocating Curbside Space Among Competing Uses

The emergence of various new forms of urban mobility services in recent years is leading to new pressures on curbside space. Municipalities, the entities typically responsible for managing the curbside, are in many instances handling these growing pressures by reallocating portions of the curbside away from traditional uses (such as metered and residential parking) in favor of uses such as ridehailing, scooter and bike-share corrals. As yet, however, such actions are being undertaken on an ad-hoc basis, due to the rapidly growing complexity of the curbside and the lack of standard analytical approaches. This lack of analytical capability is due to the traditional focus of transportation network modeling being focused predominantly on the interaction of supply and demand on links and nodes, with limited focus on link edges (the curbside). In this paper we address this research need by proposing a framework for modeling inter-modal competition for curbside space, inspired by the classical Bid-Rent Model of urban land use, intended to support curb managers to move towards maximizing the aspects of economic welfare that relate to curb access. In the bi-level model, choices made by the curbside manager impact travelers’ mode choices, and vice versa. We then present a simple numerical case study to demonstrate the properties of the proposed model, showing its tractability, flexibility, and intuitive sensitivity to systematic variation in inputs. The framework demonstrates the type of adaptive and evolving approach needed to maximize benefits from increasingly dynamic curb management strategies. The paper concludes with a brief discussion of future research needs to advance this line of inquiry.

33 ADVANCED PROPULSION SYSTEMS↗

Active operator learning with predictive uncertainty quantification for partial differential equations

With the increased prevalence of neural operators being used to provide rapid solutions to partial differential equations (PDEs), understanding the accuracy of model predictions and the associated error levels is necessary for deploying reliable surrogate models in scientific applications. Existing uncertainty quantification (UQ) frameworks employ ensembles or Bayesian methods, which can incur substantial computational costs during both training and inference. Here, we propose a lightweight predictive UQ method tailored for Deep operator networks (DeepONets) that also generalizes to other operator networks. Numerical experiments on linear and nonlinear PDEs demonstrate that the framework’s uncertainty estimates are unbiased and provide accurate out-of-distribution uncertainty predictions with a sufficiently large training dataset. Our framework provides fast inference and uncertainty estimates that can efficiently drive outer-loop analyses that would be prohibitively expensive with conventional solvers. We demonstrate how predictive uncertainties can be used in the context of Bayesian optimization and active learning problems to yield improvements in accuracy and data-efficiency for outer-loop optimization procedures. In the active learning setup, we extend the framework to Fourier Neural Operators (FNO) and describe a generalized method for other operator networks. To enable real-time deployment, we introduce an inference strategy based on precomputed trunk outputs and a sparse placement matrix, reducing evaluation time by more than a factor of five. Our method provides a practical route to uncertainty-aware operator learning in time-sensitive settings.

97 MATHEMATICS AND COMPUTING↗

Data-driven reduced-order models for port-Hamiltonian systems with operator inference

Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.

97 MATHEMATICS AND COMPUTING↗

On the role of Battery Energy Storage Systems in the day-ahead Contingency-Constrained Unit Commitment problem under renewable penetration

The integration of variable Renewable Energy Sources (vRES) to alleviate greenhouse gas emissions has introduced significant challenges for power systems operations. These challenges include high levels of uncertainty due to the intermittence associated with vRES and therefore impose the need to devise a reliable and cost-effective day-ahead unit commitment and power and reserves scheduling for real-time operations. Also, this increasing penetration of vRES requires higher ramping capabilities from units originally designed for other purposes (e.g., base-load generation), which might be exacerbated during contingency states. Hence, in this work, we propose a methodology to address the day-ahead Contingency-Constrained Unit Commitment (CCUC) problem that leverages the participation of Battery Energy Storage Systems (BESSs) to address load-following and post-contingency management, therefore alleviating the ramping burden on conventional thermal generators. To do so, we formulate a three-level optimization problem that represents the decision-making process of obtaining the least-cost commitment, generation and reserves scheduling, while restricting the Conditional Value-at-Risk (CVaR) of the system imbalance at real-time operations to user-defined tolerance levels. In addition, we devise a computationally efficient solution approach for the proposed problem based on the Column-and Constraint Generation (CCG) algorithmic framework. Two numerical experiments are conducted to empirically illustrate the benefits of the proposed methodology. Key results indicate a reduction in real-time ramping needs and a better usage of the system resources, with a reduction in the overall system commitment levels and reserve scheduling costs when compared to a benchmark case in which storage is not available.

Moreira, Alexandre↗