Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear programming”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Analysis of SNL/MSU/DOE Fatigue Database Trends for Wind Turbine Blade Materials 2010-2015.

Wind turbine blades are designed to several major structural conditions, including tip deflection, strength and b uckling during severe loading, as well as very high numbers of fatigue cycles and various service environments. The MSU Database Program has, since 1989, addressed the broad range of properties needed for current and potential blade materials through stati c and fatigue testing and test development in cooperation with Sandia National Laboratories and wind industry and supplier partners. This report is the latest in a series, giving test results and analysis for the period 2010 - 2015. Program data are compiled in a public database [1] and other reports and publications given in the cited references. The report begins with an executive summary and introductory material including background discussion of previous related studies. Section 3 describes experimental methods including processing, test methods, instrumentation and test development. Section 4 provides static tension, compression and shear stress - strain properties in three directions using coupons sectioned from a thick infused unidirectional glass/epoxy laminate. The nonlinear, shear dominated static properties were characterized with loading - u nloading - reloading (LUR) tests in tension and compression to increasing load levels, for +-45O laminates. Section 5 explores the origins of tensile fatigue sensitivity in glass fiber dominated laminates with variations in fabric architecture including speci ally prepared fabrics and aligned strand laminates. Several types of resins are considered, with variations in resin toughness and bonding to fibers, as well as cure cycle variations for an epoxy. Conclusions are drawn as to the limits of tensile fatigue r esistance and the effects of resin type and fabric architecture, including the behavior of a commercial aligned glass strand product. Interactions between cyclic fatigue response and creep are addressed for off - axis (+-45O) glass/epoxy laminates in Sectio n 6. The nonlinear fatigue and creep stress - strain and cumulative strain response are characterized in tension and compression as a function of stress level, cycles and cumulative time, using square and sinewave loading over a broad range of frequency. The results are analyzed in terms of the cycles and cumulative time under load. A cumulative strain failure criterion is established, and used to construct shear and tension constant life diagrams (CLD's) with data for nine R - values. The effects of a more duc tile urethne resin are also explored. A previous study of thick adhesives testing is extended to mixed mode fracture mechanics testing in Section 7. Mechanisms of static and fatigue crack extension near the laminate adherend interface are reported in deta il. Data are presented for mixed mode adhesive fracture, compared to mixed mode fracture in ply delamination. Fatigue crack growth exponents are also developed for a mixed mode cracked lap shear coupon. The data for fatigue trends and relative failure stra ins and exponents are compared for various blade component materials in Section 8. The effects of temperature and seawater saturation are considered for selected materials of interest for wind and hydrokinetic turbine blades in Section 9. Section 10 gives detailed conclusions for each section. A cknowledgements The research presented in this report was carried out under Sandia National Laboratories purchase orders 1325028 an d 1543945 between 2010 and 2015, with support from the DOE Wind and Water Technologies Office . In addition to the authors listed, significant contributions were made by Patrick Flaherty, Pancastya Agastra, Michael Schuster, and Michael Voth. Industry m aterials suppliers include Vectorply, Saertex, OCV, AGY, Bayer, Ashland, 3M and Nextel. Industry suppliers with significant contributions to the study were Hexion, PPG, Reichhold, Gurit and NEPTCO. Intentionally Left Blank

17 WIND ENERGY↗
DOE OSTI2012

Potential for luminosity improvement for low-energy RHIC operation

At the Brookhaven National Laboratory, a physics program, motivated by the search of the QCD phase transition critical point, requires operation of the Relativistic Heavy Ion Collider (RHIC) with heavy ions at very low beam energies corresponding to 2.5-20 GeV/n. Several physics runs were already successfully performed at these low energies. However, the luminosity is very low at lowest energies of interest (< 10 GeV/n) limited by the intra-beam scattering and space-charge, as well as by machine nonlinearities. At these low energies, electron cooling could be very effective in counteracting luminosity degradation due to the IBS, while it is less effective against other limitations. Overall potential luminosity improvement for low-energy RHIC operation from cooling is summarized for various energies, taking into account all these limitations as well as beam lifetime measured during the low-energy RHIC runs. We also explore a possibility of further luminosity improvement under the space-charge limitation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Advances in Multi-Physics Modeling for Parallel Finite Element Code Suite ACE3P

SLAC developed ACE3P code suite is a comprehensive set of parallel finite-element codes for multi-physics modeling of particle accelerators. Running on massively parallel computer platforms for high fidelity and high accuracy simulation, ACE3P enables rapid virtual prototyping of accelerator and RF component design, optimization and analysis. Furthermore, recent advances of ACE3P have been focused on the implementation of advanced numerical algorithms, enhancement of multi-physics capabilities, and improvement of code performance on state-of-the-art high performance computing (HPC) platforms for large-scale RF modeling in accelerator applications and beyond. In this paper, we will present a nonlinear electromagnetic (EM) eigensolver for damping calculation of resonant modes for cavities coupled with waveguides, the thermal solver for accurate evaluation of heat loads in superconducting RF (SRF) cavity power coupler, the elastic solver for investigation of SRF cavity cryomodule (CM) bowing due to thermal effect, the frequency domain mechanical solver for studying SRF cavity mechanical oscillation modes, the particle tracking code for dark current simulation in SRF CM using a hybrid MPI+OpenMP parallel programing model, the time domain solver for high fidelity modeling of EM wave propagation for medical device in human body.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
DOE OSTI2021

Solving the Dynamics-Aware Economic Dispatch Problem with the Koopman Operator

The dynamics-aware economic dispatch (DED) problem embeds low-level generator dynamics and operational constraints to enable near real-time scheduling of generation units in a power network. DED produces a more dynamic supervisory control policy than traditional economic dispatch (T-ED) that reduces overall generation costs. However, in contrast to T-ED, DED is a nonlinear, non-convex optimization problem that is computationally prohibitive to solve. We introduce a machine learning-based operator-theoretic approach for solving the DED problem efficiently. Specifically, we develop a novel discrete-time Koopman Operator (KO) formulation that embeds domain information into the structure of the KO to learn high-fidelity approximations of the generator dynamics. Using the KO approximation, the DED problem can be reformulated as a computationally tractable linear program (abbreviated DED-KO). We demonstrate the high solution quality and computational-time savings of the DED-KO model over the original DED formulation on a 9-bus test system.

King, Ethan↗
DOE OSTI2019

Economic Optimal Control for Minimizing Fuel Consumption of Heavy-Duty Trucks in a Highway Environment

This paper provides a comparative assessment of three economic optimal control strategies, aimed at minimizing the fuel consumption of heavy-duty trucks in a highway environment, under a representative lead vehicle model informed by traffic data. These strategies fuse a global, offline dynamic programming (DP) optimization with online model predictive control (MPC). We then show how two of the three strategies can be adapted to accommodate the presence of traffic and optimally navigate signalized intersections using infrastructure-to-vehicular (I2V) communication. The MPC optimization, which is local in nature, makes refinements to a coarsely (but globally, subject to grid resolution) optimized target velocity profile from the DP optimization. The three candidate economic MPC formulations that are evaluated include: a nonlinear time-based formulation that directly penalizes predicted fuel consumption, a nonlinear time-based formulation that penalizes braking effort as a surrogate for fuel consumption, and a linear distance-based convex formulation that maintains a tradeoff between energy expenditure and tracking of the coarsely optimized velocity profile obtained from DP. Using a medium-fidelity Simulink model, based on a Volvo truck's longitudinal and engine dynamics, we analyze the optimization's performance on four highway routes under various traffic scenarios. Results demonstrate 3.7-8.3% fuel economy improvement on highway routes without traffic and 6.5-10% on the same routes with traffic included. Furthermore, we present a detailed analysis of energy usage by "type" (aerodynamic losses, braking losses, and comparison of brake-specific fuel consumption), under each candidate control strategy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

One- and Two-Qubit Gates: Rabi Technique and Single Unipolar Pulses

Methods for controlling states of interacting superconducting flux qubits using power-efficient devices of fast single-quantum logic (Josephson nonlinearity cavities) are studied. One- and two-qubit quantum logical operations performed within the conventional control technique using Rabi pulses and using picosecond single unipolar magnetic field pulses are comparatively analyzed. It is shown that all main operations can be implemented with an accuracy of better than 97% due to optimization of the shape and parameters of unipolar control pulses (associated with, e.g., propagation of fluxons in transmission lines). The efficiency of the developed technique for programming a two-qubit quantum processor implementing the simplest Deutsch–Jozsa algorithm is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Seismic Retrofit for Electric Power Systems

Our paper develops a two-stage stochastic program and solution procedure to optimize the selection of seismic retrofit strategies to increase the resilience of electric power systems against earthquake hazards. The model explicitly considers the range of earthquake events that are possible and, for each, an approximation of the distribution of damage experienced. Furthermore, this is important because electric power systems are spatially distributed and so their performance is driven by the distribution of component damage. We also test this solution procedure against the nonlinear integer solver in LINGO 13 and apply the formulation and solution strategy to the Eastern Interconnection, where seismic hazard stems from the New Madrid seismic zone.

58 GEOSCIENCES↗
DOE OSTI2011

CESR Conversion Damping Ring Studies of Electron Cloud Instabilities (CESR-TA)

In the International Linear Collider, two linear accelerators will accelerate bunches of positrons and electrons to over a hundred billion electron volts and collide them in a central detector. In order to obtain useful collision rates, the bunches, each containing twenty billion particles, must be compressed to a cross section of a few nanometers by a few hundred nanometers. In order to prepare these ultra high density bunches, damping rings (DRs) are employed before the linear accelerators. The DRs take the high emittance bunches that are provided by the electron and positron sources and, through the process of radiation damping, squeeze them into ultra low emittance beams that are ready for the main linear accelerators. In the damping rings, a number of effects can prevent the successful preparation of the beams. In the electron ring, an effect known as the fast ion instability can lead to beam growth and, in the positron ring, the build-up of an electron cloud (EC), which interacts with the circulating bunches, can produce the same effect. EC build-up and the subsequent interaction of the cloud with the positron beam in the DR have been identified as major risks for the successful construction of a linear collider. The CESRTA research program at the Cornell Electron Storage Ring (CESR) was developed in order to study the build-up of the EC, the details of its impact on ultra low emittance beams, as well as methods to mitigate the impact of the cloud. In the DR, the EC forms when synchrotron photons radiated from the circulating beam strike the walls of the vacuum chamber, resulting in the emission of photoelectrons. These low energy electrons can be accelerated across the vacuum chamber by the electric field of the beam, and strike the walls, causing the emission of secondary electrons. The secondary electrons are subsequently accelerated into the walls yet again via the same mechanism. The result is that the EC can rapidly begin to fill the vacuum chamber. In an electron DR, the EC build-up is limited by the Coulomb repulsion. But in a positron ring, the electrons are pulled into the potential well of the beam. The resulting interaction of the circulating bunches of positrons with the EC that ultimately limits DR performance. Typically we store long trains of closely space bunches in the damping ring. The interaction of stored beam and the EC that is generated by the long train is manifested by different mechanisms: (1) The cloud focuses the beam, which causes a tune shift that increases along the bunch train as the cloud density increases. (2) The cloud electrons couple the motion of bunches along a train. Transverse motion of a leading bunch is transferred to the cloud, and subsequently to a trailing bunch which can result in a multi-bunch instability. (3) The cloud couples the positrons in the head of the bunch to those in the tail of the same bunch, which can excite a 'head-tail' instability. (4) The nonlinear fields of the EC can lead to emittance growth before the onset of instabilities. The CESRTA collaboration, which includes researchers from Cornell University's Laboratory for Elementary-Particle Physics as well as more than 50 senior staff members from over a dozen accelerator laboratories and universities around the world, has operated CESR as a damping ring for the past three years to study these EC effects. A range of specialized instrumentation has been deployed to study the local build-up of the cloud in the vacuum chambers as well as the complicated dynamics exhibited when the beam and the EC interact. The program has significantly advanced our understanding of these issues and has helped identify the most promising methods to mitigate the impact of the EC on the DR beams. It has pointed the way towards a DR design that can meet the stringent specifications of the ILC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Automated property optimization via ab initio O(N) elongation method: Application to (hyper-)polarizability in DNA

An automated property optimization method was developed based on the ab initio O(N) elongation (ELG) method and applied to the optimization of nonlinear optical (NLO) properties in DNA as a first test. The ELG method mimics a polymerization reaction on a computer, and the reaction terminal of a starting cluster is attacked by monomers sequentially to elongate the electronic structure of the system by solving in each step a limited space including the terminal (localized molecular orbitals at the terminal) and monomer. The ELG-finite field (ELG-FF) method for calculating (hyper-)polarizabilities was used as the engine program of the optimization method, and it was found to show linear scaling efficiency while maintaining high computational accuracy for a random sequenced DNA model. Furthermore, the self-consistent field convergence was significantly improved by using the ELG-FF method compared with a conventional method, and it can lead to more feasible NLO property values in the FF treatment. The automated optimization method successfully chose an appropriate base pair from four base pairs (A, T, G, and C) for each elongation step according to an evaluation function. From test optimizations for the first order hyper-polarizability (β) in DNA, a substantial difference was observed depending on optimization conditions between “choose-maximum” (choose a base pair giving the maximum β for each step) and “choose-minimum” (choose a base pair giving the minimum β). In contrast, there was an ambiguous difference between these conditions for optimizing the second order hyper-polarizability (γ) because of the small absolute value of γ and the limitation of numerical differential calculations in the FF method. It can be concluded that the ab initio level property optimization method introduced here can be an effective step towards an advanced computer aided material design method as long as the numerical limitation of the FF method is taken into account.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

NEW NEUTRON-CAPTURE MEASUREMENTS IN 23 OPEN CLUSTERS. I. THE r -PROCESS

Neutron-capture elements, those with Z > 35, are the least well understood in terms of nucleosynthesis and formation environments. The rapid neutron-capture, or r -process, elements are formed in the environments and/or remnants of massive stars, while the slow neutron-capture, or s -process, elements are primarily formed in low-mass AGB stars. These elements can provide much information about Galactic star formation and enrichment, but observational data are limited. We have assembled a sample of 68 stars in 23 open clusters that we use to probe abundance trends for six neutron-capture elements (Eu, Gd, Dy, Mo, Pr, and Nd) with cluster age and location in the disk of the Galaxy. In order to keep our analysis as homogeneous as possible, we use an automated synthesis fitting program, which also enables us to measure multiple (3–10) lines for each element. We find that the pure r -process elements (Eu, Gd, and Dy) have positive trends with increasing cluster age, while the mixed r - and s -process elements (Mo, Pr, and Nd) have insignificant trends consistent with zero. Pr, Nd, Eu, Gd, and Dy have similar, slight (although mostly statistically significant) gradients of ∼0.04 dex kpc{sup −1}. The mixed elements also appear to have nonlinear relationships with R {sub GC}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Final Seismic Shake Table Test Plan

The Spent Fuel Waste Disposition (SFWD) program is planning to conduct a full-scale seismic shake table test on the dry storage systems of spent nuclear fuel (SNF) to close the gap related to seismic loads on fuel assemblies in dry storage systems. This test will allow for quantifying the strains and accelerations on surrogate fuel assembly hardware and cladding during earthquakes of different magnitudes and frequency content. Full-scale testing is needed because a dry storage system is a complex and highly nonlinear system making it hard to predict (model) the responses to seismic excitations. The non-linearity arises from the multiple spatial gaps in the system – between fuel rods and the basket, between the basket and dry storage canister, between the dry storage canister and the storage cask (overpack), and ventilation gaps. The non-linearities pose significant limitations on the value of tests with scaled systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Nonlinear Multiobjective MPC-Based Optimal Operation of a High Consistency Refining System in Papermaking

As one of the most important unit in the papermaking industry, the high consistency (HC) refining system is confronted with challenges such as improving pulp quality, energy saving, and emissions reduction in its operation processes. Here in this correspondence, an optimal operation of HC refining system is presented using nonlinear multiobjective model predictive control strategies that aim at set-point tracking objective of pulp quality, economic objective, and specific energy (SE) consumption objective, respectively. First, a set of input and output data at different times are employed to construct the subprocess model of the state process model for the HC refining system, and then the Wiener-type model can be obtained through combining the mechanism model of Canadian Standard Freeness and the state process model that determines their structures based on Akaike information criterion. Second, the multiobjective optimization strategy that optimizes both the set-point tracking objective of pulp quality and SE consumption is proposed simultaneously, which uses NSGA-II approach to obtain the Pareto optimal set. Furthermore, targeting at the set-point tracking objective of pulp quality, economic objective, and SE consumption objective, the sequential quadratic programming method is utilized to produce the optimal predictive controllers. In conclusion, the simulation results demonstrate that the proposed methods can make the HC refining system provide a better performance of set-point tracking of pulp quality when these predictive controllers are employed. In addition, while the optimal predictive controllers orienting with comprehensive economic objective and SE consumption objective, it has been shown that they have significantly reduced the energy consumption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

AN ESTIMATE OF SPENT NUCLEAR FUEL MECHANICAL LOADS IN THE GENERAL 30 CM PACKAGE DROP SCENARIO

The US Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is performing research to determine the mechanical loading conditions applied to spent nuclear fuel (SNF) during normal conditions of transport to inform mechanical tests of SNF and close an important knowledge gap related to the practical disposition of SNF in the US. Researchers at Pacific Northwest National Laboratory (PNNL) have completed an extensive finite element study to characterize and estimate the potential mechanical loads on SNF during a hypothetical 30 cm drop of an SNF transportation package. This modeling study is validated with test data collected by the SFWST program during a physical test campaign that included one-third scale package drop tests and full-scale single fuel assembly drop tests. The test campaign was led by Sandia National Laboratories (SNL) and included international collaboration with Equipos Nucleares S.A, S.M.E (ENSA) and Bundesanstalt für Materialforschung und -prüfung (BAM). The key contribution of the modeling study is to go beyond the limitations of the limited number of physical tests to estimate the impact response to variations in impact angle, initial gap conditions, fuel assembly design, burnup and other parameters that affect the mechanical loads. The methodology of this study included a classic parametric study to calculate the impact response of highly detailed fuel assemblies over many combinations of parameters. Models of a 17x17 pressurized water reactor fuel assembly and a generic 10x10 boiling water reactor fuel assembly were both used in this study to cover the major fuel assembly types in the US inventory. Over 2,000 impact responses were calculated. The results of the parametric study were evaluated using traditional methods and basic statistics. The results were also used to construct a damage model using multiple nonlinear regression techniques to predict the mechanical loads over the full range of all input parameters. The damage model was found to work very well for all impact angle cases where the cask came to rest on its side. It was concluded that end drop cases where the cask remained vertical (instead of tipping over onto its side) were not sufficiently characterized by the current set of parametric study cases to include in the damage model, but it was not a priority to fully investigate that range because the highest mechanical loads were observed in the broader range of side impact cases. This modeling work provides sufficient insight into the mechanical loads on SNF during a hypothetical 30 cm package drop that, when considered along with the physical test data collected by the SNL-led team, the SFSWT program can consider the knowledge gap closed.

Klymyshyn, Nicholas A.↗
DOE OSTI2022

Efficient Region of Attraction Characterization for Control and Stabilization of Load Tap Changer Dynamics

In this article, we study the monitoring and control of long-term voltage stability considering load tap changer (LTC) dynamics. We show that under generic conditions, the LTC dynamics admit a unique stable equilibrium. For the stable equilibrium, we characterize an explicit inner approximation of its largest region of attraction (ROA). Compared to existing results, the computational complexity of the ROA characterization is drastically reduced. We propose a quadratically constrained linear program formulation for the ROA characterization problem. In addition, we formulate a second-order cone program for online voltage stability monitoring and control exploiting the proposed ROA characterization. Finally, we demonstrate the efficacy of the proposed formulations on the ROA characterization and stability monitoring and control using a standard IEEE test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A Performance Portable, Fully Implicit Landau Collision Operator with Batched Linear Solvers

Modern accelerators use hierarchical parallel programming models that enable massive multithreading within a processing element (PE), with multiple PEs per device driven by traditional processes. Batching is a technique for exposing PE-level parallelism in algorithms that have traditionally run on MPI processes or multiple threads within a single process. Opportunities for batching arise in, for example, kinetic discretizations of magnetized plasmas where collisions are advanced in velocity space at each spatial point independently. This paper builds on previous work on a high-performance, fully nonlinear, Landau collision operator by batching the linear solver, as well as batching the spatial point problems and adding new support for multiple grids for multiscale, multispecies problems. An anisotropic relaxation verification test that agrees well with previously published results and analytical models is presented. The performance results from NVIDIA A100 and AMD MI250X nodes are presented with hardware utilization analysis for each architecture. Finally, the entire implicit Landau operator time advance is implemented in Kokkos for performance portability, running entirely on the device and is available in the PETSc numerical library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Spatial and Spin Symmetry Breaking in Semidefinite-Programming-Based Hartree–Fock Theory

The Hartree–Fock problem was recently recast as a semidefinite optimization over the space of rank-constrained two-body reduced-density matrices (RDMs) [ Phys. Rev. A 2014, 89, 010502(R)]. This formulation of the problem transfers the nonconvexity of the Hartree–Fock energy functional to the rank constraint on the two-body RDM. We consider an equivalent optimization over the space of positive semidefinite one-electron RDMs (1-RDMs) that retains the nonconvexity of the Hartree–Fock energy expression. The optimized 1-RDM satisfies ensemble N-representability conditions, and ensemble spin-state conditions may be imposed as well. The spin-state conditions place additional linear and nonlinear constraints on the 1-RDM. We apply this RDM-based approach to several molecular systems and explore its spatial (point group) and spin (Ŝ 2 and Ŝ 3 ) symmetry breaking properties. When imposing Ŝ 2 and Ŝ 3 symmetry but relaxing point group symmetry, the procedure often locates spatial-symmetry-broken solutions that are difficult to identify using standard algorithms. For example, the RDM-based approach yields a smooth, spatial-symmetry-broken potential energy curve for the well-known Be–H 2 insertion pathway. We also demonstrate numerically that, upon relaxation of Ŝ 2 and Ŝ 3 symmetry constraints, the RDM-based approach is equivalent to real-valued generalized Hartree–Fock theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A Mixed integer linear programming‐based distributed energy management for networked microgrids considering network operational objectives and constraints

Abstract Mixed integer linear programming (MILP)–based distributed energy management for networked microgrids embedded modern distribution systems is proposed. Considering the diverse ownership of microgrids, distributed energy resources (DERs) that interface directly with utilities and responsive loads, an alternating direction method of multipliers–based distributed framework was formulated for the scheduling of networked microgrids embedded modern distribution systems by adjusting nodal price signals iteratively. In addition, to make the formulated optimization problems resolvable through more accessible and popular MILP solvers, different linearisation techniques were employed to transform the nonlinear terms into linear or mixed integer linear formats. The proposed MILP‐based distributed method preserves all participants' autonomy (e.g., microgrids, DERs that interface directly with utilities and responsive loads), while incentivising them to actively participate in the distribution system operation with price signals. The proposed method is validated with results of numerical simulation using a modern distribution system consisting of multiple networked microgrids, DERs that interface directly with utilities, as well as responsive loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗