Time-Parallel Multigrid Preconditioning for KKT Systems in Constrained Optimization
presenting at SIAM OP23 about our work on time parallel preconditioners.
Vuchkov, Radoslav Georgiev↗
SEARCH · Engineering Papers
Results for “timing system”
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.
Data-Driven Learning of Non-Autonomous Systems
We present a numerical framework for recovering unknown non-autonomous dynamical systems with time-dependent inputs. To circumvent the difficulty presented by the non-autonomous nature of the system, our method transforms the solution state into piecewise integration of the system over a discrete set of time instances. The time-dependent inputs are then locally parameterized by using a proper model, for example, polynomial regression, in the pieces determined by the time instances. This transforms the original system into a piecewise parametric system that is locally time invariant. We then design a deep neural network structure to learn the local models. Once the network model is constructed, it can be iteratively used over time to conduct global system prediction. We provide theoretical analysis of our algorithm and present a number of numerical examples to demonstrate the effectiveness of the method.
97 MATHEMATICS AND COMPUTING↗
Guerbet upgrading of ethanol to n -butanol using Ru( iii ) catalysts under air
A series of in situ prepared Ru(III) complexes supported by easily accessible N-donor organic pincer ligands were used as catalysts in the Guerbet upgrading reaction of ethanol under aerobic conditions. Tridentate bis(benzimidazole) ligand systems containing amino-dimethyl (L1) and pyridine (L2) backbones were found to make more efficient catalyst systems as compared to the bidentate bis(benzimidazole) ligand systems containing phenyl (L3) and ethene (L4) backbones. Potassium t-butoxide was found to be the most compatible base for this catalyst system. Reaction with 0.1 mol% of the catalyst and 10 mol% of potassium t-butoxide yielded 27% of n-butanol at 71% selectivity (150 °C, 24 h). Increase in the catalyst or base loading mostly resulted in increased reactivity but selectivity towards the key product n-butanol was found to decrease. On the other hand, reducing the reaction period to 12 h resulted in slightly decreased reactivity but the reaction provides n-butanol with high selectivity (76%). Contrarily, increased reaction period resulted in enhanced conversion of ethanol to higher alcohols. Under moderate and aerobic reaction conditions, the catalytic system was found to efficiently upgrade ethanol to higher alcohols. Furthermore, the decrease in the catalytic activity of the system over time was speculated to be due to the gradual deactivation of the base upon reaction with water (by-product of the Guerbet reaction). Poor solubility of the catalytic system in aqueous solutions makes it unsuitable for direct Guerbet reaction of fermentation broth.
09 BIOMASS FUELS↗
Optimal renormalization of multiscale systems
While model order reduction is a promising approach in dealing with multi-scale time-dependent systems that are too large or too expensive to simulate for long times, the resulting reduced order models can suffer from instabilities. We have recently developed a time-dependent renormalization approach to stabilize such reduced models. In the current work, we extend this framework by introducing a parameter that controls the time-decay of the memory of such models and optimally selecting this parameter based on limited fully resolved simulations. First, we demonstrate our framework on the inviscid Burgers equation whose solution develops a finite-time singularity. Our renormalized reduced order models are stable and accurate for long times while using for their calibration only data from a full order simulation before the occurrence of the singularity. Furthermore, we apply this framework to the 3D Euler equations of incompressible fluid flow, where the problem of finite-time singularity formation is still open and where brute force simulation is only feasible for short times. Our approach allows us to obtain for the first time a perturbatively renormalizable model which is stable for long times and includes all the complex effects present in the 3D Euler dynamics. We find that, in each application, the renormalization coefficients display algebraic decay with increasing resolution, and that the parameter which controls the time-decay of the memory is problem-dependent.
Price, Jacob↗
Large Scale Simulations of Quantum Systems on HPC with Analytics for HEP Algorithms
The objective of this work is to combine HEP simulation experience with HPC quantum systems simulations to provide an environment for executing large-scale experiments with advanced analysis capabilities for studying multi-qubit and multi-level reduced-noise, long coherence time quantum systems such as 3D Superconducting RF (SRF) cavities integrated with transmon qubits.
97 MATHEMATICS AND COMPUTING↗
Resilient Operations of Networked Microgrids (RONM)
Objective: Improve the resiliency of power systems with optimization-based methods that leverage advanced microgrid technologies to reduce system recovery times after extreme event induced outages. Outcome: First-of-kind, high-fidelity physics-based optimization method for modeling networked microgrids which includes key engineering constraints associated with system recovery after extreme events.
24 POWER TRANSMISSION AND DISTRIBUTION↗
PNNL DSSE Docker Image
This is the docker image for Pacific Northwest National Laboratory's DSSE (distribution system state estimator) used for the demo of OEDI-SI platform. To support the operation of modern distribution systems, operators require real-time visibility into system states. Due to a lack of measurements and unbalanced operation, the state estimation in distribution systems is challenging. This submission is related to an OEDI-SI use case which demonstrates the application of an extended Kalman filter-based state estimator on an IEEE-123 bus system. The state estimator uses the measurements to generate voltage estimates for the system.
Array↗
Measuring and Simulating T1 and T2 for Qubits
Quantum computers perform computations exploiting quantum mechanics to a possible advantage, allowing us to prepare and manipulate states that do not have a classical equivalent. In particular, phenomena like superposition and entanglement may enable quantum computers to outperform their classical counterparts in certain applications. Implementing these useful quantum algorithms is contingent upon building accurate quantum hardware that is not affected by noise. Environmental noise decreases coherence time of qubits, meaning that qubits do not stay in a desired state long enough to carry out a complex computation. To that end, harnessing the full power of quantum computers necessitate characterization of noise sources and how they impact a given quantum system. Often times, T1and T2 are used to quantify noise. In this project, we provide an approach as to how T1 and T2 values are calculated and simulated for quantum systems. In addition, we compare simulated values of T1 and T2 with those of real quantum computer’s measurements. IBMQ Experience, an open source software allowing users to simulate and use real quantum hardware, is used. QuTip, a Python-based toolbox offering quantum simulation tools for open quantum systems, is also used.
Youssef, Rahaf↗
Digital Grid Twin–Direct Communication Scheme Test Bed for Assessing Relay-to-Relay Radio Antenna and Optical Fiber Performance and Misoperations
This study introduces a novel “Digital Grid Twin–Direct Communication Scheme” test bed. This advanced platform evaluates point-to-point communication between transmitter and receiver relays with optical fiber and radio omnidirectional antenna systems, implemented at the Advanced Protection lab in the Grid Research Innovation and Development Center at Oak Ridge National Laboratory. The increased diversity of energy sources has led to more protective relay misoperations. In North America, microgrid protection schemes now use point-to-point communication along distribution lines between relays to implement advanced logic in nonradial grids that include both high- and low-inertia generators. This trend challenges utilities to minimize misoperations while ensuring rapid fault clearance and accurate selectivity coordination between primary and backup relays. This study assesses relay-to-relay communication schemes by introducing an advanced testing platform based on a digital grid twin protection test bed using a synchronized time source system. The platform evaluates the communication system using radio antennas or optical fiber links by integrating protective relays that operate breakers within the digital twin and record relay events and communication signals. In the experiments, transmitter and receiver relays were configured with inverse time overcurrent and breaker trip detection logic to assess the total time of the communication protection schemes based on the sum of the relay protection element operating time, radio latency, propagation delay, baud rate delay, and relay processing time. These delays were derived from recorded relay events and communication signals from the interface of a real-time simulator set as a digital grid twin. The test bed successfully simulated various electrical faults along a distribution line while ensuring effective and reliable point-to-point communication between transmitter and receiver relays. The radio antenna communication system exhibited latency because of the radio. This latency depends on the baud rate setting and type of radio application; in general, the higher the baud rate, the lower the radio latency. The measured radio latency (for Mirrored Bits with an encryption card at 9,600 bps) was about 9–10 ms. Additionally, calculated propagation delay per mile for radio antennas and optical fiber was 5.36 µs/mi and 8.04 µs/mi, respectively. Optical fiber communication did not demonstrate radio latency. Instead, the protection element operating time depends mainly on the protection logic function set in the relay, and the relay processing time depends on the processing rate of the relay in samples per power system cycle.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Nonlinear Elimination Applied to Radiation Diffusion
We apply a nonlinearly preconditioned, quasi-Newton framework to accelerate the numerical solution of the thermal radiative transfer (TRT) equations. This framework was inspired by the unpublished method that has existed for years in Teton, Lawrence Livermore National Laboratory’s deterministic TRT code. In this paper, we cast this iteration scheme within a formal nonlinear preconditioning framework and compare its performance against other iteration schemes in the framework. With proper choices of iteration controls for the various levels of the solver, we can recover the standard linearized one-step method, a full nonlinear Newton scheme, as well as the method in Teton. In brief, the nonlinear preconditioning TRT scheme formally eliminates the material temperature equation from the nonlinear system in a nonlinear analog of a Schur complement. This nonlinear elimination step involves solving a decoupled nonlinear equation for each spatial degree of freedom and is therefore inexpensive. By applying a quasi-Newton iteration scheme on the new system, we obtain a three-level iteration scheme that is at least as efficient as commonly used TRT schemes. The new method allows full convergence to the nonlinear backward Euler time-discretized system, increasing accuracy and robustness, while using a similar number of linear iterations as the more common linearized one-step methods Eq. (4).
77 NANOSCIENCE AND NANOTECHNOLOGY↗
GPS Spoofing Mitigation and Timing Risk Analysis in Networked Phasor Measurement Units via Stochastic Reachability
To address phasor measurement unit (PMU) vulnerability to spoofing, we propose the use of a set-valued state estimation technique known as stochastic reachability (SR)-based distributed Kalman filter (DKF) that computes secure global positioning system (GPS) timing across a network of receivers. Utilizing SR, we estimate not only GPS time but also its stochastic reachable set, which is parameterized by probabilistic zonotope (p-Zonotope). While requiring known measurement error bounds in only non-spoofed conditions, we designed a two-tiered approach. We first performed measurement-level spoofing mitigation via deviation of a measurement innovation from its expected p-Zonotope. We then performed state-level timing risk analysis via a determination of the intersection probability of the estimated p-Zonotope with an unsafe set that violates IEEE C37.118.1a-2014 standards. Finally, we validated our SR-DKF algorithm by subjecting it to a simulated receiver network to coordinate signal-level spoofing. We demonstrate improved timing accuracy and successful spoofing mitigation via the use of our SR-DKF algorithm. We also validated the robustness of the estimated timing risk as the number of receivers were varied.
47 OTHER INSTRUMENTATION↗
Summary of GNU Octave Lady Godiva Simulator
A pseudo real-time interactive reactor simulator for the Lady Godiva reactor has been developed using the GNU Octave software package. The simulator tracks the time-dependent power and fuel temperature behavior of the Lady Godiva reactor during sub-prompt critical transients. Users interact with the simulator by controlling reactivity input. This input occurs through users repositioning a control rod as desired during simulator operation. Various reactor parameters are displayed graphically in near “real time,” facilitating user assessment and feedback input. The reactor model consists of a lumped-fuel core, a single control rod, and a simple heat transfer mechanism that allows heat to be rejected from the core. The goal of the simulator is to demonstrate in real-time the system response to reactivity perturbations, the effect of reactivity feedback, and the effect of temperature-dependent equilibrium conditions. The dynamic behavior of the reactor is primarily driven by the relative characteristic time-constants between the delayed neutrons and the overall heat transfer parameter. The simulator will allow users to gain experience performing standard reactor operations, such as establishing steady-state at different fuel temperatures and making power level changes
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Real-World Cyber Security Demonstration for Networked Electric Drives
In this article, we present the design and implementation of a cyber-physical security testbed for networked electric drive systems, aimed at conducting real-world security demonstrations. To our knowledge, this is one of the first security testbeds for networked electric drives, seamlessly integrating the domains of power electronics and computer science, and cybersecurity. By doing so, the testbed offers a comprehensive platform to explore and understand the intricate and often complex interactions between cyber and physical systems. The core of our testbed consists of four electric machine drives, meticulously configured to emulate small-scale but realistic information technology (IT) and operational technology (OT) networks. This setup both provides a controlled environment for simulating a wide array of cyber-attacks, and mirrors potential real-world attack scenarios with a high degree of fidelity. The testbed serves as an invaluable resource for the study of cyber-physical security, offering a practical and dynamic platform for testing and validating cybersecurity measures in the context of networked electric drive systems. As a concrete example of the testbed's capabilities, we have developed and implemented a Python-based script designed to execute step-stone attacks over a wireless local area network (WLAN). This script leverages a sequence of target IP addresses, simulating a real-world attack vector that could be exploited by adversaries. To counteract such threats, we demonstrate the efficacy of our developed cyber-attack detection algorithms, which are integral to our testbed's security framework. Furthermore, the testbed incorporates a real-time visualization system using InfluxDB and Grafana, providing a dynamic and interactive representation of networked electric drives and their associated security monitoring mechanisms. This visualization component not only enhances the testbed's usability but also offers insightful, real-time data for researchers and practitioners, thereby facilitating a deeper understanding of cyber-physical security dynamics in networked electric drive systems.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Experimental and Modeling Review of the Plutonium-Zirconium (Pu-Zr) System: Lost in Translation and Over Time?
This article presents a review of the phase equilibria of the Pu-Zr system and discusses the contradictory reports on the experimental phase diagram and misinterpretations that have led to confusion over time. In addition, a review of the few Pu-Zr CALPHAD (CALculation of PHAse Diagrams) assessments is presented with emphasis on the heat of formation of the bcc phase, highlighting a disagreement between CALPHAD and ab initio calculations. Based on the information gathered in this study, a critical re-assessment of the Pu-Zr system involving the ζ (Pu 28 Zr), θ (Pu 4 Zr), kinetically hindered κ (PuZr 2 ), and once observed ι phases is called for.
36 MATERIALS SCIENCE↗
Microcontroller-based aerosol jet printer control software
This software is used with a microcontroller to control an aerosol jet printing system. It was written for a 32-bit ARM microcontroller, providing greater processing speed and complexity relative to more conventional low-cost printer controllers. The microcontroller software was written using a real-time operating system (FreeRTOS), which simplifies parallel operation of multiple functions and further customization. It supports simultaneous motion planning, pulse generation for stepper motors, reading encoder feedback, communicating with multiple mass flow controllers, controlling solid state relays, handling a software-based emergency stop, and sending data to the client PC. The program is written in C and communicates with the client PC via USB connection. Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-0436 O
Secor, EthanBenjamin↗
The Fermilab SSR1 and HB650 Synoptic and EPICS remote cryogenic control system
The 325MHz Single Spoke Resonator 1 (SSR1) Cryomodule (CM) has a 5.3-meter length, and it has 288-liter liquid helium @ 4.5K mode. The High Beta 650 MHz (HB650) CM contains six elliptical cavity cryomodules, has a 9.9-meter overall length and it has 530-liter liquid helium. The SSR1 and HB650 CM had been tested at Fermilab Cryomodule Test Facility (CMTF). Their cryogenic control system includes the Siemens Process Control System S7-400, Automation Direct DL205 PLC, and remote synoptic HMI/EPICS Phoebus GUI. This paper presents a method which has been successfully used by Fermilab PIP2 IT cryogenic remote on-line, real-time control systems.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
The Impact of Inverter Dead-Time in Single-Phase Wireless Power Transfer Systems
In this paper, the effect of the dead-time on a single-phase wireless power transfer system is studied in detail. In practice, the dead-time is always placed between the complementary switching pulses of the inverter phase-leg. At higher operating frequencies, the dead-time issues in the resonant inverter become critical, especially as the power level increases. The detailed analysis of the dead-time on a wireless power transfer system is discussed for different operating conditions of the inverter duty-cycle and power factor. The switching characteristics of the wireless power transfer system inverter are analyzed, and the notch phenomenon that appears at the output of the inverter is also discussed. A notch equation based on the observations is derived to predict the notch occurrence during the system operation. Furthermore, the mathematical expressions are presented for different notch conditions. Subsequently, the effect of the notches on the sensitivity and the power transfer of the series-series compensated wireless power transfer system is analyzed. Finally, the approach is verified experimentally on an 8 kW wireless power transfer system prototype, and the results are compared with the theoretical analysis.
42 ENGINEERING↗
A portable and reusable sensor system based on graphene for real-time and sensitive detection of lead ions in water
Long-term exposure to Pb 2+ can cause irreversible damage to the nervous, cardiovascular, and reproductive systems. Therefore, developing a fast and sensitive detection system capable of monitoring minuscule concentrations of Pb 2+ is essential. In this study, we demonstrated a fully portable sensor system enabling rapid, sensitive, and real-time monitoring of Pb 2+ . The sensor system adopted the remote-gate field-effect transistor (RGFET) detection scheme and was easy to operate, even for non-experts. The sensor system comprised two printed circuit boards (PCBs): a sensor PCB with a remote-gate electrode and an analyzer PCB with a metal-oxide-semiconductor field-effect transistor (MOSFET) transducer and peripheral electronics to manage sensor signals. To achieve a high sensitivity for Pb 2+ , we utilized graphene ink drop-casted on the sensor PCB as a sensing membrane. The graphene film was easy to deposit and remove, enabling the sensor PCB to be reused multiple times. The sensor system was further linked to a smartphone application that instantly monitors the sensor response, allowing for rapid point-of-use detection. The sensor exhibited a high sensitivity of 21.7% when the limit of detection (LOD) value of 1 nM (∼0.2 ppb) was detected, and the typical detection time for each sample was approximately 60 seconds. This portable sensor system advances sensing technologies and could potentially supplement expensive, laborious conventional sensing equipment.
54 ENVIRONMENTAL SCIENCES↗