Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Real numbers”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Adaptive Sampling for In Situ Cloud Probe (Final Report)

Clouds play a leading role in the Earth's global energy and solar radiation balance and hydrological cycle. Improving cloud models requires detailed information on the cloud microphysical properties, such as droplet size distribution and number density, liquid water content and cloud composition (droplets, ice particles), which can only be provided by aerial in situ measurements. However, for many atmospheric measurement instruments, the lack of flexibility in selecting the operational mode during operation can lead to uncertainties in sampling and measurement characteristics under continuously varying atmospheric conditions. This SBIR project is developing an advanced, compact optical imaging technology for in situ characterization of cloud hydrometeors. The development involves a deep modification of the existing Mesa Photonics’ Cloud Droplet Measurement System (CDMS) in order to implement real-time automatic adaptive sampling based on the acquired in situ data and environmental parameters. The new system, CDMS-2, implements two measurement modes: side-scatter imaging for smaller hydrometeors and direct bright-field-illumination imaging for larger hydrometeors in a significantly larger sample volume. The system measures the droplet size distribution (DSD) and number density with an added capability of discriminating between liquid water and ice hydrometeors (based on polarization-resolved side-scatter imaging). The instrument will implement automatic switching or alternating between the regular side-scatter imaging mode and sparse/large hydrometeor mode (based on the acquired data). Other adaptive sampling capabilities include variable sample volume and dynamic range (based on the measured DSD). The preferred deployment platforms are uncrewed aircraft systems (UAS) and tethered balloon/kite systems (TBS). The Phase I project achieved (or exceeded) the goals listed in the Work Plan. A CDMS-2 laboratory prototype implementing the polarization-resolved side-scatter imaging mode and direct bright-field-illumination imaging mode was designed and built. Additional capabilities included the variable illumination pulse energy and sample volume. The smallest detectable droplet diameter was improved to 3–4 μm (from the nominal 10 μm value specified for the original CDMS). Discrimination between water droplets and ice particles was experimentally demonstrated. The Phase I prototype was extensively tested and calibrated in the laboratory and also tested in the Pi Cloud Chamber at Michigan Technological University (MTU). The two intensive experimental campaigns at MTU provided unique opportunities of testing the CDMS-2 laboratory prototype under realistic warm and mixed-phase cloud conditions (stable for long periods of time), testing different sampling modes and intercomparing the CDMS-2 prototype to other co-located cloud characterization instruments. The Phase I project successfully demonstrated the feasibility of the proposed technology and identified the engineering challenges of designing a field deployable prototype instrument in Phase II. The Phase I study provides a solid basis for development, characterization and field-testing of the proposed advanced cloud probe with adaptive sampling in Phase II followed by commercialization of the technology in Phase III.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Quantification of Lanthanides in a Molten Salt Reactor Surrogate Off-Gas Stream Using Laser-Induced Breakdown Spectroscopy

To enable the deployment of molten salt reactor technology, the development of off-gas treatment systems and advanced monitoring tools capable of operating with high temperatures and radiation fields while delivering near real-time information is necessary. Here, we aim to fulfill this requirement and proposes laser-induced breakdown spectroscopy (LIBS) for monitoring molten salt aerosol streams. A sheath gas measuring method was developed to protect optical elements from aerosol particles and to ensure a relatively constant aerosol stream for measurement. An aqueous system was studied to demonstrate the utility of LIBS for monitoring possible fission products in an aerosol stream: Gd, Nd, and Sm up to 2000 parts per million (ppm). A calibration model was built using partial least squares (PLS) regression with five, six, and nine latent variables for Gd, Nd, and Sm, respectively. This calibration model successfully estimated the concentrations of three test samples, which were validated with inductively charged plasma optical emission spectroscopy measurements at a 99.9% confidence interval. To enhance these models, a genetic algorithm was used to filter the spectra before entering the PLS model, thereby limiting the spectral features being regressed to those with greater correlations to concentration. This allowed for the number of latent variables used in the PLS models to be reduced to four, three, and three for Gd, Nd, and Sm, respectively. Finally, the genetic algorithm-filtered PLS models were used to predict the concentrations of the aerosol stream on a real-time dataset and resulted in a 73%, 18%, and 25% improvement in root mean squared error of prediction compared to the original PLS models developed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Non-cancellation of the parity anomaly in the strong-field regime of QED 2+1

Quantum fluctuations lead to an anomalous violation of parity symmetry in quantum electrodynamics for an even number of spatial dimensions. While the leading parity-odd electric current vanishes in vacuum, we uncover a non-cancellation of the anomaly for strong electric fields with distinct macroscopic signatures. We perform real-time lattice simulations with fully dynamical gauge fields and Wilson fermions in 2 + 1 space-time dimensions. In the static field limit, relevant at early times, we solve the problem analytically. Our results point out the fundamental role of quantum anomalies for strong-field phenomena, relevant for a wide range of condensed matter and high-energy applications, but also for the next generation of gauge theory quantum simulators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Decoupling copolymer, lipid and carbon nanotube interactions in hybrid, biomimetic vesicles

Bilayer vesicles that mimic a real biological cell can be tailored to carry out a specific function by manipulating the molecular composition of the amphiphiles. These bio-inspired and bio-mimetic structures are increasingly being employed for a number of applications from drug delivery to water purification and beyond. Complex hybrid bilayers are the key building blocks for fully synthetic vesicles that can mimic biological cell membranes, which often contain a wide variety of molecular species. While the assembly and morpholgy of pure phospholid bilayer vesicles is well understood, the functionality and structure dramaticlly changes when copolymer and/or carbon nanotube porins (CNTP) are added. The aim of this study is to understand how the collective molecular interactions within hybrid vesicles affect their nanoscale structure and properties. In situ small and wide angle X-ray scattering (SAXS/WAXS) and molecular dynamics simulations (MD) are used to investigate the morphological effect of molecular interactions between polybutadiene polyethylene oxide, lipids and carbon nanotubes (CNT) within the hybrid vesicle bilayer. Within the lipid/copolymer system, the hybrid bilayer morphology transitions from phase separated lipid and compressed copolymer at low copolymer loadings to a mixed bilayer where opposing lipids are mostly separated from the inner region. This transition begins between 60 wt% and 70 wt%, with full homogenization observed by 80 wt% copolymer. The incorporation of CNT into the hybrid vesicles increases the bilayer thickness and enhances the bilayer symmetry. Analysis of the WAXS and MD indicate that the CNT–dioleoyl interactions are much stronger than the CNT–polybutadiene.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Comprehensive Scheduling Framework using SP-ADMM for Residential Demand Response with Weather and Consumer Uncertainties

This paper presents a comprehensive scheduling framework for residential demand response (DR) programs considering both the day-ahead and real-time electricity markets. In the first stage, residential customers determine the operating status of their responsive devices such as heating, ventilation, and air conditioning (HVAC) systems and electric water heaters (EWHs), while the distribution system operator (DSO) computes the amount of electricity to be purchased in the day-ahead electricity market. In the second stage, the DSO purchases insufficient (or sells surplus) electricity in the real-time electricity market to maintain the supply-demand balance. Due to its computational complexity and data privacy issues, the proposed model cannot be directly solved in a centralized manner, especially with a large number of uncertain scenarios. Therefore, this paper proposes a combination of stochastic programming (SP) and the alternating direction method of multipliers (ADMM) algorithm, called SP-ADMM, to decompose the original model and then solve each sub-problem in a distributed manner while considering multiple uncertain scenarios. The simulation study is performed on the IEEE 33-bus system including 121 residential houses. Here, the results demonstrate the effectiveness of the proposed approach for large-scale residential DR applications under weather and consumer uncertainties.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Power Electronics Manufacturing Improvements for Heavy-Duty Fuel Cell Vehicles

The Marel Power Solutions project, funded by the U.S. Department of Energy under Award DE-SC0023801, focused on advancing manufacturing techniques for power electronics in heavy-duty fuel cell vehicles. The research aimed to enhance system efficiency, reduce costs, and support broader adoption of hydrogen fuel cell technology. Key areas of investigation included power topology, thermal modeling, system architecture, and accessibility through software tools. Key Accomplishments: 1. Power Topology: - Developed an interleaved boost converter with optimized phase count, leveraging Marel’s proprietary Power Stacks. - Achieved reduced parasitic inductance and resistance, enabling high efficiency in DC-DC converters. 2. Thermal Modeling: - Integrated innovative cooling systems into compact Silicon Carbide (SiC) modules. - Simulations demonstrated the ability to dissipate significant heat (up to 7.5 kW), ensuring device reliability under heavy loads. 3. System Architecture: - Utilized simulation tools to analyze the impact of various fuel cell and vehicle parameters on efficiency. - Highlighted the role of smaller, modular improvements, such as enhanced DC-DC converters, in achieving system-wide gains. 4. Accessibility: - Evaluated and implemented MATLAB/Simulink code generation tools for real-world hardware applications. - Demonstrated the potential for rapid prototyping of custom power systems with reduced development costs. Impact and Benefits: - Efficiency and Cost Reduction: Marel’s cooling technology enhances SiC die performance, reducing the number of dies required and overall system size. - Scalability and Flexibility: The innovations support tailored solutions for diverse applications, from mass transit to mining vehicles. - Sustainability: The research promotes the integration of electrification technologies, helping meet rising energy demands sustainably. Conclusion: The project’s outcomes advance the state of power electronics for hydrogen fuel cell vehicles, enabling more efficient, compact, and cost-effective solutions. These developments lay a foundation for future innovation, contributing to the broader adoption of clean energy technologies in transportation and other industries.

08 HYDROGEN↗
DOE OSTI2024

PV Generation and Load Forecasting for Adjuntas PR Community Microgrids

Existing frameworks to forecast time-series photovoltaic (PV) output power and consumer load for microgrid operations and controls assume a near-continuous availability of real-time input features from the field assets such as PV inverters, energy meters, and weather station. These incoming data points are used to periodically retrain models and update forecast snapshots over a moving horizon window, be it one hour-ahead, one-day ahead, or one-week ahead. However, such frameworks are not resilient to disruptions in data availability caused by losses in communications between the field sensors and data loggers. Hence, there is a need for programs that assume no availability of real-time microgrid asset data and still make reliable forecasts that can be used for decision-making. Such programs would be apt to function in extreme weather events such as hurricanes and would use lightweight recursive time-series models to independently forecast solar irradiance and ambient temperature, then compute PV power from those forecasts, as well as independently forecast consumer load. The codebase performs forecasting for the scenario of when the microgrid does not have a reliable access to forecasts or real-time observations of solar irradiance (I) and ambient temperature (AT) and load (Load) to be able to adequately forecast, in real-time, the PV power production or a business' load. In this case, using historical values of PV power and load, a univariate forecasting of generation and consumption are respectively made. The use-case in particular has two sub-scenarios: one, a normal 7-day ahead forecast where the unavailability of real-time data is assumed due to infrastructure issues such as loss of communication or sensor maintenance or service downtimes. Whereas a hurricane-caused unavailability of real-time data requires a second model trained specifically on historical hurricane days to be able to capture the extreme day behavior of generation in particular, and load if applicable. A gradient boosted regression tree comprises an ensemble of additive models that map between the input of historical values (be it irradiance, temperature, or load) and their corresponding output forecasts of a given horizon such that the individual learner predictions are summed up over the total number of such learners in the ensemble to produce an aggregate forecast. A weighting mechanism is applied to the training data in each iteration, where actual and forecast values are compared to penalize incorrect forecasts by increasing the weight and reducing it to reward correct forecasts. The code's benefits are that it: (a) accounts for a contingency where communication loss renders newly measured real-time data unavailable for model tuning and snapshot updates; (b) presents blind forecasting that recursively determines the next time-step value in a horizon using the forecast of the same attribute from a prior step; and (c) employs lightweight models that, once trained, can reliably generalize for different horizons, which make them suitable for enhancing the resilience of field microgrids prone to extreme events that encounter disruptions to data availability.

Sundararajan, Aditya [Oak Ridge National Laborator↗
DOE OSTI2025

Smart Contracts for Power Grid Applications Using the Advanced DLT Cyber Grid Guard Testbed

In this study is presented two power system applications with distributed ledger technology (DLT) and smart contracts (SC) that were assessed in a Cyber-Grid-Guard System (CGGS) advanced testbed, with protective relays, power meters, communication devices, DLT devices, synchronized time source, clock displays and real time simulator. This CGGS testbed was set in the Advanced Protection Lab, 252 lab space of the Grid Research Integration and Deployment Center (GRID- C), at Oak Ridge National Laboratory. In power grids, customer-owned distributed energy resources (DERs) are more frequent than in the past, and the numbers of points of interconnection (POI) with customer-owned DERs have increased. Disruptive operation from DERs presents a risk to grid operations, and protective relays located at the POI are used to isolate out-of-tolerance or poorly behaving of DERs. Ensuring the integrity of data from the relays at the POI, and DLT could enhance the security of the power grids. The first application is a SC to define and control the allowable total power factor (TPF) of the DER (wind farm) output, and the terms of the SC are implemented using DLT with a CGGS for a customer-owned DER. The TPF SC was implemented by the CGGS using DLT. The experimental model was performed with a real-time simulator using a CGGS and relay in-the-loop. The data collected from the CGGS were used to execute the TPF SC. The TPF limits were between +0.9 and +1.0, and the breakers’ operation in the POI was controlled by the relay using the SC. The events were collected from the real-time simulator, CGGS, and SEL 700GT relay to validate a successful application of the TPF SC using DLT. The second application is a SC to measure and control the allowable voltage service limits (VSL) by the CGGS using DLT. The tests were performed by using a real-time simulator, CGGS and relay in-the-loop. The data was collected from the CGGS that executed the SC. The main constraints were defined based on ANSI C84.1 service voltage limits, and the operation of the breakers in the POI. The events were collected from the CGGS, and SEL 700GT relay to assess a successful operation of the VSL SC using DLT.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

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

On the meaning of Berry force for unrestricted systems treated with mean-field electronic structure

We show that the Berry force as computed by an approximate, mean-field electronic structure can be meaningful if properly interpreted. In particular, for a model Hamiltonian representing a molecular system with an even number of electrons interacting via a two-body (Hubbard) interaction and a spin–orbit coupling, we show that a meaningful nonzero Berry force emerges whenever there is spin unrestriction—even though the Hamiltonian is real-valued and formally the on-diagonal single-surface Berry force must be zero. Moreover, if properly applied, this mean-field Berry force yields roughly the correct asymptotic motion for scattering through an avoided crossing. That being said, within the context of a ground-state calculation, several nuances do arise as far interpreting the Berry force correctly, and as a practical matter, the Berry force diverges near the Coulson–Fischer point (which can lead to numerical instabilities). Here, we do not address magnetic fields here.

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

Star-spots and magnetism: testing the activity paradigm in the Pleiades and M67

We measure star-spot filling fractions for 240 stars in the Pleiades and M67 open star clusters using APOGEE high-resolution H-band spectra. For this work, we developed a modified spectroscopic pipeline which solves for star-spot filling fraction and star-spot temperature contrast. We exclude binary stars, finding that the large majority of binaries in these clusters (80 per cent) can be identified from Gaia DR3 and APOGEE criteria – important for field star applications. Our data agree well with independent activity proxies, indicating that this technique recovers real star-spot signals. In the Pleiades, filling fractions saturate at a mean level of 0.248 ± 0.005 for active stars with a decline at slower rotation; we present fitting functions as a function of Rossby number. In M67, we recover low mean filling fractions of 0.030 ± 0.008 and 0.003 ± 0.002 for main sequence GK stars and evolved red giants, respectively, confirming that the technique does not produce spurious spot signals in inactive stars. Star-spots also modify the derived spectroscopic effective temperatures and convective overturn time-scales. Effective temperatures for active stars are offset from inactive ones by −109 ± 11 K, in agreement with the Pecaut & Mamajek empirical scale. Star-spot filling fractions at the level measured in active stars changes their inferred overturn time-scale, which biases the derived threshold for saturation. Finally, we identify a population of stars statistically discrepant from mean activity–Rossby relations and present evidence that these are genuine departures from a Rossby scaling. Our technique is applicable to the full APOGEE catalogue, with broad applications to stellar, galactic, and exoplanetary astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2020

A Study of Imbalance Levels Attributed to Photovoltaic Penetration in Distribution Systems

Abstract— The number of renewable sources, namely photovoltaics (PV), continues to grow within distribution systems. How placements impact a system’s performance remains a concern within both the research and the utility communities. Beyond system node voltages and real power losses, this work looks to industry standard imbalance metrics and the changes experienced on these values with varied levels of PV. Specifically, four imbalance metrics are presented and studied for an actual 906-bus distribution system. Multiple loading conditions and photovoltaic placement schemes are explored to compare the changes seen on these metrics across a variety of scenarios experienced by the network.

14 SOLAR ENERGY↗
DOE OSTI2022

Support Vector Machines for Classification of Direct Energy Deposition Standoff Distance for Improved Process Control

A critical factor in the implementation of direct energy deposition is the ability to maintain the standoff distance between the nozzle and the build surface, as this influences powder capture efficiency and overall part quality. Due to process-related variations, layer height may vary, causing unintended variation in standoff distance and poor build quality. While prior work has utilized contact probing to qualify standoff distance during processing, in situ methods for qualification of standoff distance are of major interest. The present work seeks to understand efficacy of image-based methods for classifying standoff distance variation in real-time using support vector machines (SVMs). It was hypothesized that the size of the melt pool and the amount of spatter will have significant correlations with deviations in the standoff distance; thus, SVMs were used on a dataset that is comprised of morphological features of melt pool size and image entropy. The SVM model was used to classify melt pool images into categories according to standoff distance variation from nominal. K-folds cross validation was used to find the optimal hyperparameters for the SVM model. To understand the impact of the selected features on the classification performance and inference speed, multiple models were trained with differing numbers of included features. Results for classification score, inference time, and image preprocessing/feature extraction from these data are reported. The present results show that the SVM model was able to predict the standoff distance classification with an accuracy of 97 percent and a speed of 0.122 s per image, making it a viable solution for real-time control of standoff distance.

Klesmith, Zoe↗
DOE OSTI2021

ATLAS toward the High Luminosity era: challenges on electronic systems

To maximize the physics reach, the Large Hadron Collider (LHC) plans to increase its instantaneous luminosity to $7.5×10^{34} cm^{-2} s^{-1}$ and deliver $3 - 4 ab^{-1}$ of data at a center-of-mass energy of 14 TeV. To profit from this high-luminosity LHC operation with a pile-up up to 200 inelastic collisions per bunch crossing, performance of the ATLAS detector needs to be maintained and in many systems, improved. The upgraded detector will have an unprecedented output data rate of up to 200 TB/s. Real-time processing of this large data volume in a short time period is extremely challenging. New sets of both front-end and back-end electronics are required for all sub-detectors. The large number of detector channels, huge volumes of input and output data, short time available to process and transmit data, harsh radiation environment, and the need of low power consumption all impose great challenges on the designs of electronic systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

FY2021 Isolated Grids and Grid-Connected Turbine Reference Systems

For individuals, businesses, and communities focused on building resilient electrical grid infrastructure, wind energy can provide an affordable, accessible, and compatible distributed energy resource option that also enhances the capabilities of local grid operations. However, there are technical barriers to realizing the market value and resilience benefits of distributed wind, and there is little to no ability to quantify those benefits so that stakeholders can compare grid investment options. The central aims of this report are: (1) to drive technology transfer of the methods and technologies developed under the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) project and (2) increase the number of referenceable case studies available to stakeholders interested in additional value-added capabilities of wind systems beyond bulk energy supply. We achieve this aim by applying three major methods developed under MIRACL to two real-world distributed wind reference systems. The two real-world distributed wind reference systems are the isolated grid of St. Mary’s, Alaska, and the two 10.5-megawatt (MW) front-of-the-meter wind turbine deployments owned and operated by Iowa Lakes Electric Cooperative (ILEC).

17 WIND ENERGY↗
DOE OSTI2022

Cyber risk assessment and investment optimization using game theory and ML-based anomaly detection and mitigation for wide-area control in smart grids

The electric power grid is increasingly becoming susceptible to cyber attacks that exploit vulnerabilities in the smart grid control, information, and physical layers. Successful cyber attacks can have catastrophic impacts on the social and economic well-being of any nation all over the globe. It has, thus, become imperative to secure the smart grid against such adversarial actions to ensure stable, secure, and reliable operation of the grid. The existing research and industry practices prove to be inadequate in terms of providing pragmatic and effective defense methodologies and measures for long-term cybersecurity planning and real-time cybersecurity for grid operation. For example, existing works lack models that incorporate uncertain behavior of cyber-attackers and pragmatic defense measures for cyber risk assessment and cybersecurity investment optimization which often provide unreliable and strictly qualitative solutions to these problems. At the same time, with the growing number of cyber incidents in the grid, there still exists a need to develop attack-resilient algorithms for wide-area monitoring, protection, and control (WAMPAC) applications like the wide-area voltage control systems (WAVCS) for Flexible AC Transmissions Systems (FACTS) that lack in scalable and feasible solutions from the cybersecurity perspective. This dissertation proposes novel models and methodologies for: (1) Cybersecurity planning, and (2) Cybersecurity for system operation. The cybersecurity planning is achieved through cyber risk assessment and cybersecurity resource investment optimization for long-term cybersecurity of the grid using game theory and attack-defense trees. Cybersecurity for system operation consists of development of cyber anomaly detection and mitigation algorithms for flexible AC transmission system (FACTS) controller-based wide-area voltage control systems (WAVCS) using machine learning (ML), and software defined networking-based moving target defense network routing for achieving real-time cyber-physical security for grid operations. This is followed by hardware-in-the-loop (HIL) implementation and evaluation of these attack prevention, detection, and mitigation algorithms and methodologies showcasing their feasibility in a close to real-world environment. For cybersecurity planning, a novel approach involving a combination of game theory and attack defense trees (ADT) for optimal cybersecurity resource allocation in the smart grid is proposed. This methodology involves modeling of the cyber-physical smart grid substations as ADTs, defining attacker costs, defense costs, and attack probabilities for attack access points. Using game theoretical formulation, optimal defense strategies for the defender of the system to invest cybersecurity resources in the grid are obtained. Additionally, a game-theoretic framework is developed for quantitative cyber-physical risk assessment of the grid under a dynamically changing cyber threat space and uncertain behavior of cyber attackers which is further used to optimize investments in the smart grid's cybersecurity resources. The attacker, defender, and the smart grid system are modeled while incorporating attacker-stochasticity and federal guidelines for smart grid cybersecurity. This allows quantification of threat, vulnerabilities, and attack impact of the grid for quantitative risk assessment. The defender's budget to invest in the security resources in the grid is optimized based on the strategies leading to minimum system risk. The evaluation of the proposed solutions highlight the feasibility for practical implementation of these methodologies and algorithms in the smart grid, while taking the federal requirements and guidelines for smart grid security into consideration. For achieving cybersecurity for system operation, attack prevention, detection, and mitigation algorithms and methodologies are developed specifically for FACTS-based WAVCS. Anomaly detection and mitigation in the WAVCS are achieved using algorithms based on machine learning which involves offline training and testing of ML models with CPS datasets incorporating physics-based features that allow accurate distinction between system faults and cyber attacks. For attack prevention, a methodology based on software defined network (SDN)-based moving target defense (MTD) network routing is proposed that enables prevention of Denial of Service (DoS) type attacks on the smart grid communication system. Subsequently, these methodologies and algorithms are implemented and evaluated on an HIL testbed that allows for real-time attack prevention, detection, and mitigation of emulated cyber attacks on the WAVCS in a close to real-world environment. The results show highly accurate and efficient performance of the implemented algorithms and methodologies with the smart grid system operating within the NERC's system operation limits even in the presence of DoS and data integrity cyber attacks. This work opens up future research opportunities in other directions such as (1) Expanding cybersecurity planning methodologies to real-time cyber contingency analysis with different game formulations; and (2) Applying the cybersecurity for system operation algorithms to broader categories of wide-area control applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Constrained non-negative matrix factorization enabling real-time insights of in situ and high-throughput experiments

Non-negative matrix factorization (NMF) is an appealing class of methods for performing unsupervised learning on streaming spectral data, particularly in time-sensitive applications such as in situ characterization of materials. These methods seek to decompose a dataset into a small number of components and weights that can compactly represent the underlying signal while effectively reconstructing the observations with minimal error. However, canonical NMF methods have no underlying requirement that the reconstruction uses components or weights that are representative of the true physical processes. In this work, we demonstrate how constraining a subset of the NMF weights or components as rigid priors, provided as known or assumed values, can provide significant improvement in revealing true underlying phenomena. We present a PyTorch-based method for efficiently applying constrained NMF and demonstrate its application to several synthetic examples. Our implementation allows an expert researcher-in-the-loop to provide and dynamically adjust the constraints during a live experiment involving streaming spectral data. Such interactive priors allow researchers to specify known or identified independent components, as well as functional expectations about the mixing or transitions between the components. We further demonstrate the application of this method to measured synchrotron x-ray total scattering data from in situ beamline experiments. In such a context, constrained NMF can result in a more interpretive and scientifically relevant decomposition than canonical NMF or other decomposition techniques. As a result, the details of the method are provided, along with general guidance for employing constrained NMF in the extraction of critical information and insights during time-sensitive experimental applications.

36 MATERIALS SCIENCE↗