Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable 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 127 records · Page 7

Source-to-Source Automatic Differentiation of OpenMP Parallel Loops

This article presents our work toward correct and efficient automatic differentiation of OpenMP parallel worksharing loops in forward and reverse mode. Automatic differentiation is a method to obtain gradients of numerical programs, which are crucial in optimization, uncertainty quantification, and machine learning. The computational cost to compute gradients is a common bottleneck in practice. For applications that are parallelized for multicore CPUs or GPUs using OpenMP, one also wishes to compute the gradients in parallel. Here, we propose a framework to reason about the correctness of the generated derivative code, from which we justify our OpenMP extension to the differentiation model. We implement this model in the automatic differentiation tool Tapenade and present test cases that are differentiated following our extended differentiation procedure. Performance of the generated derivative programs in forward and reverse mode is better than sequential, although our reverse mode often scales worse than the input programs.

97 MATHEMATICS AND COMPUTING↗

Dynamic Programming and Error Estimates for Stochastic Control Problems with Maximum Cost

This work is concerned with stochastic optimal control for a running maximum cost. A direct approach based on dynamic programming techniques is studied leading to the characterization of the value function as the unique viscosity solution of a second order Hamilton–Jacobi–Bellman (HJB) equation with an oblique derivative boundary condition. A general numerical scheme is proposed and a convergence result is provided. Error estimates are obtained for the semi-Lagrangian scheme. These results can apply to the case of lookback options in finance. Moreover, optimal control problems with maximum cost arise in the characterization of the reachable sets for a system of controlled stochastic differential equations. Some numerical simulations on examples of reachable analysis are included to illustrate our approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗

Certification for Small Wind Turbine Installers: What's the Hang Up?; Preprint

Several programs have been implemented to support the advancement of a professional, mature small wind industry and to ensure that this industry moves forward in a sustainable direction. The development of a standard for small wind turbine systems and the creation of the Small Wind Certification Council support small wind technology that is reliable and safe. Consumers and incentive programs willultimately rely on certification to differentiate among systems sold in the U.S. market. Certification of small wind installers is yet another component deemed necessary for this industry to expand. The National Renewable Energy Laboratory, under the guidance and funding support of the U.S. Department of Energy, supported the development of small wind system installer certification provided viathe North American Board of Certified Energy Practitioners. However, the small wind community is not supportive of the installer certification. There are currently only nine certified installers in the U.S. pool. This paper provides an overview of the installer certification program and why more small wind turbine installers are not pursuing this certification.

17 WIND ENERGY↗

Integrated Data Collection Analysis (IDCA) Program - SSST Testing Methods

The Integrated Data Collection Analysis (IDCA) program is conducting a proficiency study for Small- Scale Safety and Thermal (SSST) testing of homemade explosives (HMEs). Described here are the methods used for impact, friction, electrostatic discharge, and differential scanning calorimetry analysis during the IDCA program. These methods changed throughout the Proficiency Test and the reasons for these changes are documented in this report. The most significant modifications in standard testing methods are: 1) including one specified sandpaper in impact testing among all the participants, 2) diversifying liquid test methods for selected participants, and 3) including sealed sample holders for thermal testing by at least one participant. This effort, funded by the Department of Homeland Security (DHS), is putting the issues of safe handling of these materials in perspective with standard military explosives. The study is adding SSST testing results for a broad suite of different HMEs to the literature. Ultimately the study will suggest new guidelines and methods and possibly establish the SSST testing accuracies needed to develop safe handling practices for HMEs. Each participating testing laboratory uses identical test materials and preparation methods wherever possible. The testing performers involved are Lawrence Livermore National Laboratory (LLNL), Los Alamos National Laboratory (LANL), Indian Head Division, Naval Surface Warfare Center, (NSWC IHD), Sandia National Laboratories (SNL), and Air Force Research Laboratory (AFRL/RXQL). These tests are conducted as a proficiency study in order to establish some consistency in test protocols, procedures, and experiments and to compare results when these testing variables cannot be made consistent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Repository performance confirmation.

Repository performance confirmation links the technical bases of repository science and societal acceptance. This paper explores the myriad aspects of what has been labeled performance confirmation in U.S. programs, which involves monitoring as a collection of distinct activities combining technical and social significance in radioactive waste management. This paper is divided into four parts: (1) A distinction is drawn between performance confirmation monitoring and other testing and monitoring objectives; (2) A case study illustrates confirmation activities integrated within a long-term testing and monitoring strategy for Yucca Mountain; (3) A case study reviews compliance monitoring developed and implemented for the Waste Isolation Pilot Plant; and (4) An approach for developing, evaluating and implementing the next generation of performance confirmation monitoring is presented. International interest in repository monitoring is exhibited by the European Commission Seventh Framework Programme 'Monitoring Developments for Safe Repository Operation and Staged Closure' (MoDeRn) Project. The MoDeRn partners are considering the role of monitoring in a phased approach to the geological disposal of radioactive waste. As repository plans advance in different countries, the need to consider monitoring strategies within a controlled framework has become more apparent. The MoDeRn project pulls together technical and societal experts to assimilate a common understanding of a process that could be followed to develop a monitoring program. A fundamental consideration is the differentiation of confirmation monitoring from the many other testing and monitoring activities. Recently, the license application for Yucca Mountain provided a case study including a technical process for meeting regulatory requirements to confirm repository performance as well as considerations related to the preservation of retrievability. The performance confirmation plan developed as part of the Yucca Mountain license application identified a broad suite of monitoring activities. A revision of the plan was expected to winnow the number of activities down to a manageable size. As a result, an objective process for the next stage of performance confirmation planning was developed as an integral part of an overarching long-term testing and monitoring strategy. The Waste Isolation Pilot Plant compliance monitoring program at once reflects its importance to stakeholders while demonstrating adequate understanding of relevant monitoring parameters. The compliance criteria were stated by regulation and are currently monitored as part of the regulatory rule for disposal. At the outset, the screening practice and parameter selection were not predicated on a direct or indirect correlation to system performance metrics, as was the case for Yucca Mountain. Later on, correlation to performance was established, and the Waste Isolation Pilot Plant continues to monitor ten parameters originally identified in the compliance certification documentation. The monitoring program has proven to be effective for the technical intentions and societal or public assurance. The experience with performance confirmation in the license application process for Yucca Mountain helped identify an objective, quantitative methodology for this purpose. Revision of the existing plan would be based on findings of the total system performance assessment. Identification and prioritization of confirmation activities would then derive from performance metrics associated with performance assessment. Given the understanding of repository performance confirmation, as reviewed in this paper, it is evident that the performance confirmation program for the Yucca Mountain project could be readily re-engaged if licensing activities resumed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characterization and thermal behavior of PrMO{sub 3} (M = Co or Ni) ceramic materials obtained from gelatin

Graphical abstract: The micrograph in figure shows sample calcined at temperature 900 °C. The sample exhibits morphology with considerable porosity and the formation of agglomerated nanometric particles. Gelatin provides the system with a large amount of organic matter, which is then removed during calcinations, favoring the appearance of pores in the material. Highlights: ► Oxides with PrNiO{sub 3} and PrCoO{sub 3} were prepared by new method synthesis. ► The gelatin, through its carboxylate groups and amine, is an efficient director. ► The obtained materials have magnetic properties and application in catalysis. ► The decomposition kinetic study of bonding groups of gelatin with metallic ions that takes part in the synthesis of PrMO{sub 3}. -- Abstract: Metal oxides with perovskite-type structure have attracted considerable interest in recent years due to their magnetic and electrical properties, as well as their catalytic activity. In this study, oxides with PrNiO{sub 3} and PrCoO{sub 3} composition were prepared by using gelatin powder as a precursor agent for its use as a catalyst. The powders obtained were calcined at 700 °C and 900 °C and characterized using the X-ray diffraction, thermal analysis (thermogravimetry and differential thermal analysis), infrared spectroscopy, temperature programed reduction and scanning electron microscopy techniques. Thermogravimetric data using the non-isothermal kinetic models of Flynn and Wall and “Model-free Kinetics” were used to determine the activation energy to study the decomposition kinetics of the ligand groups with system's metallic ions that takes part in the synthesis of PrMO{sub 3} (M = Ni or Co).

36 MATERIALS SCIENCE↗

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗

Metric Learning to Accelerate Convergence of Operator Splitting Methods

Recent developments in machine learning have led to promising advances in accelerating the solution of constrained optimization problems. Increasing demand for real-time decision-making capabilities in applications such as artificial intelligence and optimal control has led to a variety of proposed strategies for learning to produce fast solutions to optimization problems. For example, recent works have shown that it is possible to accelerate the convergence of optimization algorithms by learning to select their parameters, such as gradient descent stepsizes. This work proposes a new approach, in which the underlying metric spaces of proximal operator splitting algorithms are learned to maximize convergence rate. While prior works in optimization theory have derived optimal metrics in simple cases, no such result exists for many practical problem forms including general Quadratic Programming (QP). This paper shows how differentiable optimization can enable the end-to-end learning of proximal metrics, enhancing the convergence of proximal algorithms for QP problems beyond what is possible based on known theory. Additionally, the results illustrate a strong connection between the learned proximal metrics and active constraints at the optima, leading to an interpretation in which the predicted proximal metrics can be viewed as a form of active set prediction.

King, Ethan [BATTELLE (PACIFIC NW LAB)]↗

ASC CSSE Level 2 Milestone #6362: Resilient Asynchronous Many Task Programming Model

This report is an outcome of the ASC CSSE Level 2 Milestone 6362: Analysis of Re- silient Asynchronous Many-Task (AMT) Programming Model. It comprises a summary and in-depth analysis of resilience schemes adapted to the AMT programming model. Herein, performance trade-offs of a resilient-AMT prograrnming model are assessed through two ap- proaches: (1) an analytical model realized by discrete event simulations and (2) empirical evaluation of benchmark programs representing regular and irregular workloads of explicit partial differential equation solvers. As part of this effort, an AMT execution simulator and a prototype resilient-AMT programming framework have been developed. The former permits us to hypothesize the performance behavior of a resilient-AMT model, and has undergone a verification and validation (V&V) process. The latter allows empirical evaluation of the perfor- mance of resilience schemes under emulated program failures and enabled the aforementioned V&V process. The outcome indicates that (1) resilience techniques implemented within an AMT framework allow efficient and scalable recovery under frequent failures, that (2) the abstraction of task and data instances in the AMT programming model enables readily us- able Application Program Interfaces (APIs) for resilience, and that (3) this abstraction enables predicting the performance of resilient-AMT applications with a simple simulation infrastruc- ture. This outcome will provide guidance for the design of the AMT programming model and runtime systems, user-level resilience support, and application development for ASC's next generation platforms (NGPs).

97 MATHEMATICS AND COMPUTING↗

THE YOUNG SOLAR ANALOGS PROJECT. I. SPECTROSCOPIC AND PHOTOMETRIC METHODS AND MULTI-YEAR TIMESCALE SPECTROSCOPIC RESULTS

This is the first in a series of papers presenting methods and results from the Young Solar Analogs Project, which began in 2007. This project monitors both spectroscopically and photometrically a set of 31 young (300–1500 Myr) solar-type stars with the goal of gaining insight into the space environment of the Earth during the period when life first appeared. From our spectroscopic observations we derive the Mount Wilson S chromospheric activity index (S{sub MW}), and describe the method we use to transform our instrumental indices to S{sub MW} without the need for a color term. We introduce three photospheric indices based on strong absorption features in the blue-violet spectrum—the G-band, the Ca i resonance line, and the Hydrogen-γ line—with the expectation that these indices might prove to be useful in detecting variations in the surface temperatures of active solar-type stars. We also describe our photometric program, and in particular our “Superstar technique” for differential photometry which, instead of relying on a handful of comparison stars, uses the photon flux in the entire star field in the CCD image to derive the program star magnitude. This enables photometric errors on the order of 0.005–0.007 magnitude. We present time series plots of our spectroscopic data for all four indices, and carry out extensive statistical tests on those time series demonstrating the reality of variations on timescales of years in all four indices. We also statistically test for and discover correlations and anti-correlations between the four indices. We discuss the physical basis of those correlations. As it turns out, the “photospheric” indices appear to be most strongly affected by emission in the Paschen continuum. We thus anticipate that these indices may prove to be useful proxies for monitoring emission in the ultraviolet Balmer continuum. Future papers in this series will discuss variability of the program stars on medium (days–months) and short (minutes to hours) timescales.

79 ASTRONOMY AND ASTROPHYSICS↗

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗

Assembly and Delivery of Rabbit Capsules for Irradiation of Silicon Carbide Cladding Tube Specimens in the High Flux Isotope Reactor

Neutron irradiation of silicon carbide (SiC)-based fuel cladding under a high radial heat flux presents a critical challenge for SiC cladding concepts in light water reactors (LWRs). Fission heating in the fuel provides a high heat flux through the cladding, which, combined with the degraded thermal conductivity of SiC under irradiation, results in a large temperature gradient through the thickness of the cladding. The strong temperature dependence of swelling in SiC creates a complex stress profile in SiCbased cladding tubes as a result of differential swelling. The Nuclear Science User Facilities (NSUF) Program within the US Department of Energy Office of Nuclear Energy is supporting research efforts to improve the scientific understanding of the effects of irradiation on SiC cladding tubes. Ultimately, the results of this project will provide experimental validation of multi-physics models for SiC-based fuel cladding during LWR operation. The first objective of this project is to irradiate tube specimens using a previously developed design that allows for irradiation testing of miniature SiC tube specimens subjected to a high radial heat flux. The previous “rabbit” capsule design uses the gamma heating in the core of the High Flux Isotope Reactor (HFIR) to drive a high heat flux through the cladding tube specimens. A compressible aluminum foil allows for a constant thermal contact conductance between the cladding tubes and the rabbit housing despite swelling of the SiC tubes. To allow separation of the effects of irradiation from those due to differential swelling under a high heat flux, a new design was developed under the NSUF program. This design allows for irradiation of similar SiC cladding tube specimens without a high radial heat flux. This report briefly describes the irradiation experiment design concepts, summarizes the irradiation test matrix, and reports on the successful delivery of six rabbit capsules to the HFIR. Rabbits of both low and high heat flux configurations have been assembled, welded, evaluated, and delivered to the HFIR along with a complete quality assurance fabrication package. These rabbits contain a wide variety of specimens including monolith tubes, SiC fiber SiC matrix (SiC/SiC) composites, duplex specimens (inner composite, outer monolith), and specimens with a variety of metallic or ceramic coatings on the outer surface. The rabbits are targeted for insertion during HFIR cycle 475, which is scheduled for September 2017.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hes6 is required for actin cytoskeletal organization in differentiating C2C12 myoblasts

Hes6 is a member of the hairy-enhancer-of-split family of transcription factors that regulate proliferating cell fate in development and is known to be expressed in developing muscle. Here we investigate its function in myogenesis in vitro. We show that Hes6 is a direct transcriptional target of the myogenic transcription factors MyoD and Myf5, indicating that it is integral to the myogenic transcriptional program. The localization of Hes6 protein changes during differentiation, becoming predominantly nuclear. Knockdown of Hes6 mRNA levels by siRNA has no effect on cell cycle exit or induction of myosin heavy chain expression in differentiating C2C12 myoblasts, but F-actin filament formation is disrupted and both cell motility and myoblast fusion are reduced. The knockdown phenotype is rescued by expression of Hes6 cDNA resistant to siRNA. These results define a novel role for Hes6 in actin cytoskeletal dynamics in post mitotic myoblasts.

60 APPLIED LIFE SCIENCES↗

pyomo.dae: a modeling and automatic discretization framework for optimization with differential and algebraic equations

We describe pyomo.dae, an open source Python-based modeling framework that enables high-level abstract specification of optimization problems with differential and algebraic equations. The pyomo.dae framework is integrated with the Pyomo open source algebraic modeling language, and is available at http://www.pyomo.org. One key feature of pyomo.dae is that it does not restrict users to standard, predefined forms of differential equations, providing a high degree of modeling flexibility and the ability to express constraints that cannot be easily specified in other modeling frameworks. Other key features of pyomo.dae are the ability to specify optimization problems with high-order differential equations and partial differential equations, defined on restricted domain types, and the ability to automatically transform high-level abstract models into finite-dimensional algebraic problems that can be solved with off-the-shelf solvers. Moreover, pyomo.dae users can leverage existing capabilities of Pyomo to embed differential equation models within stochastic and integer programming models and mathematical programs with equilibrium constraint formulations. Collectively, these features enable the exploration of new modeling concepts, discretization schemes, and the benchmarking of state-of-the-art optimization solvers.

97 MATHEMATICS AND COMPUTING↗

Single-Crystal Sapphire Optical Fiber Sensor Instrumentation

This report summarizes technical progress on the program “Single-Crystal Sapphire Optical Fiber Sensor Instrumentation,” funded by the National Energy Technology Laboratory of the U.S. Department of Energy, and performed by the Center for Photonics Technology of the Bradley Department of Electrical and Computer Engineering at Virginia Tech. This project was completed in three phases, each with a separate focus. Phase I of the program, from October 1999 to April 2002, was devoted to development of sensing schema for use in high temperature, harsh environments. Different sensing designs were proposed and tested in the laboratory. Phase II of the program, from April 2002 to April 2009, focused on bringing the sensor technologies, which had already been successfully demonstrated in the laboratory, to a level where the sensors could be deployed in harsh industrial environments and eventually become commercially viable through a series of field tests. Also, a new sensing scheme was developed and tested with numerous advantages over all previous ones in Phase II. Phase III of the program, September 2009 to December 2013, focused on development of the new sensing scheme for field testing in conjunction with materials engineering of the improved sensor packaging lifetimes. In Phase I, three different sensing principles were studied: sapphire air-gap extrinsic Fabry-Perot sensors; intensity-based polarimetric sensors; and broadband polarimetric sensors. Black body radiation tests and corrosion tests were also performed in this phase. The outcome of the first phase of this program was the selection of broadband polarimetric differential interferometry (BPDI) for further prototype instrumentation development. This approach is based on the measurement of the optical path difference (OPD) between two orthogonally polarized light beams in a single-crystal sapphire disk. At the beginning of Phase II, in June 2004, the BPDI sensor was tested at the Wabash River coal gasifier facility in Terre Haute, Indiana. Due to business conditions at industrial partner and several logistical problems, this field test was not successful. An alternative high-temperature sensing system using sapphire wafer-based extrinsic Fabry-Perot interferometry was then developed as a significant improvement over the BPDI solution. From June 2006 to June 2008, three consecutive field tests were performed with the new sapphire wafer sensors at the TECO coal gasifier in Tampa, Florida. One of the sensors survived in the industrial coal gasifier for 7 months, over which time the existing thermocouples were replaced twice. The outcome of these TECO field tests suggests that the sapphire wafer sensor has very good potential to be commercialized. However packaging and sensor protection issues need additional development. During Phase III, several major improvements in the design and fabrication process of the sensor have been achieved through experiments and theoretical analysis. Studies on the property of the key components in the sensor head, including the sapphire fiber and sapphire wafer, were also conducted, for a better understanding of the sensor behavior. A final design based on all knowledge and experience has been developed, free of any issues encountered during the entire research. Sensors with this design performed well as expected in lab long-term tests, and were deployed in the sensing probe of the final coal-gasifier field test. Sensor packaging and protection was improved through materials engineering through testing of packaging designs in two blank probe packaging tests at Eastman Chemical in Kingsport, TN. Performance analysis of the blank probe packaging resulted in improve package designs culminating in a 3rd generation probe packaging utilized for the full field test of the sapphire optical sensor and materials designed sensor packaging.

47 OTHER INSTRUMENTATION↗

State elimination for mixed‐integer optimal control of partial differential equations by semigroup theory

Abstract Mixed‐integer optimal control problems governed by partial differential equations (MIPDECOs) are powerful modeling tools but also challenging in terms of theory and computation. We propose a highly efficient state elimination approach for MIPDECOs that are governed by partial differential equations that have the structure of an abstract ordinary differential equation in function space. This allows us to avoid repeated calculations of the states for all time steps, and our approach is applied only once before starting the optimization. The presentation of theoretical results is complemented by numerical experiments.

97 MATHEMATICS AND COMPUTING↗

Applying Quantum Computing to Simulate Power System Dynamics

Power system dynamics are generally modeled by high dimensional nonlinear differential-algebraic equations due to a large number of generators, loads, and transmission lines. Thus, its computational complexity grows exponentially with the system size. This paper demonstrates the potential use of quantum computing algorithms to model the power system dynamics. Leveraging a symbolic programming framework, we equivalently convert the power system dynamics’ differential algebraic equations (DAEs) into ordinary differential equations (ODEs), where the data of the state vector can be encoded into quantum computers via amplitude encoding. The system's nonlinearity is captured by Taylor polynomial expansion, the quantum state tensor, and Hamiltonian simulation, whereas state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can simulate the dynamics of the power system with high accuracy, whereas its complexity is polynomial in the logarithm of the system dimension. Our work also illustrates the use of scientific machine learning tools for implementing scientific computing concepts, e.g., Taylor expansion, DAEs/ODEs transform, and quantum computing solver, in the field of power engineering.

Tran, Huynh↗