Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Integral”

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 793 records · Page 44

Integrated Electrocoagulation, Ultrafiltration, Membrane Distillation, and Crystallization for Treating Produced Water

Produced water (PW) generated from hydraulic fracturing operations was treated using an integrated electrocoagulation, ultrafiltration, membrane distillation, and crystallization process (EC UF MDC). The aim was to determine the viability of this integrated process for maximizing water recovery. The results obtained here indicate that optimizing the various unit operations could lead to increased recovery of PW. Membrane fouling limits all membrane separation processes. A pretreatment step to suppress fouling is essential. Here, removal of total suspended solids (TSS) and total organic carbon (TOC) was achieved by electrocoagulation (EC) followed by ultrafiltration (UF). The hydrophobic membrane used in membrane distillation may be fouled by dissolved organic compounds. Reducing membrane fouling is essential to increase the long-term durability of the membrane distillation (MD) system. In addition, combining membrane distillation with crystallization (MDC) can help reduce scale formation. By inducing crystallization in the feed tank, scale formation on the MD membrane was suppressed. The integrated EC UF MDC process can impact Water Resources/Oil & Gas Companies. Conservation of surface and groundwater is possible by treating and reusing PW. Additionally, treating PW reduces the amount of PW disposed in Class II disposal wells and promotes more environmentally sustainable operations.

42 ENGINEERING↗

Integrating Radio-Over-Fiber Communication System and BOTDR Sensor System

In this paper, we propose and experimentally demonstrate for the first time, the integration of a radio-over-fiber (RoF) communication system and a Brillouin optical time-domain reflectometry (BOTDR) distributed sensor system using a single optical fiber link. In this proof-of-concept integrated system, the communication system is composed of three modulation formats of quadrature phase-shift keying (QPSK), 16-quadrature amplitude modulation (16-QAM) and 64-QAM, which are modulated onto an orthogonal frequency division multiplexing (OFDM) signal. Whereas, the BOTDR sensor system is used for strain and/or temperature monitoring over the fiber distance with a spatial resolution of 5 m using a 25 km single-mode silica fiber. The error vector magnitude (EVM) is analyzed in three modulation formats in the presence of various BOTDR input pump powers. Using QPSK modulation, optimized 18 dBm sensing and 10 dBm data power, the measured EVM values with and without bandpass filter are 3.5% and 14.5%, respectively. The proposed system demonstrates a low temperature measurement error (±0.49 °C at the end of 25 km) and acceptable EVM values, which were within the 3GPP requirements. The proposed integrated system can be effectively applied for practical applications, which significantly reduces the fiber infrastructure cost by effective usage of a single optical fiber link.

42 ENGINEERING↗

Integration of Decentralized Graph-Based Multi-Agent Reinforcement Learning with Digital Twin for Traffic Signal Optimization

Machine learning (ML) methods, particularly Reinforcement Learning (RL), have gained widespread attention for optimizing traffic signal control in intelligent transportation systems. However, existing ML approaches often exhibit limitations in scalability and adaptability, particularly within large traffic networks. This paper introduces an innovative solution by integrating decentralized graph-based multi-agent reinforcement learning (DGMARL) with a Digital Twin to enhance traffic signal optimization, targeting the reduction of traffic congestion and network-wide fuel consumption associated with vehicle stops and stop delays. In this approach, DGMARL agents are employed to learn traffic state patterns and make informed decisions regarding traffic signal control. The integration with a Digital Twin module further facilitates this process by simulating and replicating the real-time asymmetric traffic behaviors of a complex traffic network. The evaluation of this proposed methodology utilized PTV-Vissim, a traffic simulation software, which also serves as the simulation engine for the Digital Twin. The study focused on the Martin Luther King (MLK) Smart Corridor in Chattanooga, Tennessee, USA, by considering symmetric and asymmetric road layouts and traffic conditions. Comparative analysis against an actuated signal control baseline approach revealed significant improvements. Experiment results demonstrate a remarkable 55.38% reduction in Eco_PI, a developed performance measure capturing the cumulative impact of stops and penalized stop delays on fuel consumption, over a 24 h scenario. In a PM-peak-hour scenario, the average reduction in Eco_PI reached 38.94%, indicating the substantial improvement achieved in optimizing traffic flow and reducing fuel consumption during high-demand periods. These findings underscore the effectiveness of the integrated DGMARL and Digital Twin approach in optimizing traffic signals, contributing to a more sustainable and efficient traffic management system.

42 ENGINEERING↗

A Causal Approach to Integrate Component Health Data into System Reliability Models

Two of the challenges of current plant reliability approaches are the ability to integrate plant health data, and to support decision making. Condition based data and diagnostic/prognostic information are in fact not considered into plant reliability models to inform system engineers on the most critical components. Currently, the propagation of quantitative health data from the component to the system level is a challenge given the diverse nature/structure of the data. On the other hand, plant reliability methods (which are typically based on fault-trees or reliability block diagrams) can effectively propagate data from the component to the system level, but values of failure rates or failure probabilities are an approximated integral representation of the past industry-wide operational experience, and it neglects the present component health status (e.g., diagnostic and condition-based data) and health projection (when available from prognostic data). Our first claim is that system reliability models should propagate health information from the component to the system/plant level in order to provide a quantitative snapshot of system/plant health and identify the most critical components. Our second claim is that component health should be informed solely by that specific component current and historical performance data and should not be an approximated integral representation of the past industry-wide operational experience. This paper is directly supporting these two claims by proposing a different approach to perform reliability modeling which relies on available component diagnostic, prognostic and condition-based data to measure component health, and it propagates this information through fault tree models. The propagation of health data from the component to the system level is performed not in terms of probability, but in terms of margins where margin is defined as the “distance” between the present actual status and an undesired event (e.g., failure or unacceptable performance). Through a cause-effect lens, while classical reliability models target the effect associated to a component performance, a margin-based approach focuses on the cause of an undesired component performance (i.e., component health). Hence, thinking of reliability in terms of margins implies decision making based on causal reasoning. We will show how fault tree models can be solved using a margin language and how this process can effectively assist system engineers to identify the most critical components.

97 MATHEMATICS AND COMPUTING↗

Energy-dependent and Energy-integrated Two-moment General-relativistic Neutrino Transport Simulations of a Hypermassive Neutron Star

Abstract We compare two-moment-based energy-dependent and three variants of energy-integrated neutrino transport general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. To study the impacts due to the choice of the neutrino transport schemes, we perform simulations with the same setups and input neutrino microphysics. We show that the main differences between energy-dependent and energy-integrated neutrino transport are found in the disk and ejecta properties, as well as in the neutrino signals. The properties of the disk surrounding the neutron star and the ejecta in energy-dependent transport are very different from the ones obtained using energy-integrated schemes. Specifically, in the energy-dependent case, the disk is more neutron-rich at early times and becomes geometrically thicker at later times. In addition, the ejecta is more massive and, on average, more neutron-rich in the energy-dependent simulations. Moreover, the average neutrino energies and luminosities are about 30% higher. Energy-dependent neutrino transport is necessary if one wants to better model the neutrino signals and matter outflows from neutron star merger remnants via numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

An Integral-based Technique to Accelerate the Monte Carlo Radiative Transfer Computation for Supernovae

We present an integral-based technique (IBT) algorithm to accelerate supernova (SN) radiative transfer calculations. The algorithm utilizes “integral packets,” which are calculated by the path integral of the Monte Carlo (MC) energy packets, to synthesize the observed spectropolarimetric signal at a given viewing direction in a 3D time-dependent radiative transfer program. Compared to the event-based technique (EBT) proposed by M. Bulla et al., our algorithm significantly reduces the computation time and increases the MC signal-to-noise ratio (S/N). Using a 1D spherical symmetric Type Ia SN ejecta model DDC10 and its derived 3D model, the IBT algorithm has successfully passed the verification of spherical symmetry and cross comparison on a 3D SN model with the direct-counting technique and EBT. Notably, with our algorithm implemented in the 3D MC radiative transfer code SEDONA, the computation time is faster than EBT by a factor of 10−30, and the S/N is better by a factor of 1.5−3, with the same number of MC quanta.

79 ASTRONOMY AND ASTROPHYSICS↗

Practical Effects of Integrating Temperature with Strang Split Reactions

For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.

79 ASTRONOMY AND ASTROPHYSICS↗

A Fully Explicit Integrator for Modeling Astrophysical Reactive Flows

Simulating complex astrophysical reacting flows is computationally expensive—reactions are stiff and typically require implicit integration methods. The reaction update is often the most expensive part of a simulation, which motivates the exploration of more economical methods. In this research note, we investigate how the explicit Runge–Kutta–Chebyshev (RKC) method performs compared to an implicit method when applied to astrophysical reactive flows. These integrators are applied to simulations of X-ray bursts arising from unstable thermonuclear burning of accreted fuel on the surface of neutron stars. We show that the RKC method performs with similar accuracy to our traditional implicit integrator, but is more computationally efficient when run on CPUs.

GPU computing↗

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗

A Fortran–Python interface for integrating machine learning parameterization into earth system models

Abstract. Parameterizations in earth system models (ESMs) are subject to biases and uncertainties arising from subjective empirical assumptions and incomplete understanding of the underlying physical processes. Recently, the growing representational capability of machine learning (ML) in solving complex problems has spawned immense interests in climate science applications. Specifically, ML-based parameterizations have been developed to represent convection, radiation, and microphysics processes in ESMs by learning from observations or high-resolution simulations, which have the potential to improve the accuracies and alleviate the uncertainties. Previous works have developed some surrogate models for these processes using ML. These surrogate models need to be coupled with the dynamical core of ESMs to investigate the effectiveness and their performance in a coupled system. In this study, we present a novel Fortran–Python interface designed to seamlessly integrate ML parameterizations into ESMs. This interface showcases high versatility by supporting popular ML frameworks like PyTorch, TensorFlow, and scikit-learn. We demonstrate the interface's modularity and reusability through two cases: an ML trigger function for convection parameterization and an ML wildfire model. We conduct a comprehensive evaluation of memory usage and computational overhead resulting from the integration of Python codes into the Fortran ESMs. By leveraging this flexible interface, ML parameterizations can be effectively developed, tested, and integrated into ESMs.

54 ENVIRONMENTAL SCIENCES↗

A Fortran-Python Interface for Integrating Machine Learning Parameterization into Earth System Models

Parameterizations in Earth System Models (ESMs) are subject to biases and uncertainties arising from subjective empirical assumptions and incomplete understanding of the underlying physical processes. Recently, the growing representational capability of machine learning (ML) in solving complex problems has spawned immense interests in climate science applications. Specifically, ML-based parameterizations have been developed to represent convection, radiation and microphysics processes in ESMs by learning from observations or high-resolution simulations, which have the potential to improve the accuracies and alleviate the uncertainties. Previous works have developed some surrogate models for these processes using ML. These surrogate models need to be coupled with the dynamical core of ESMs to investigate the effectiveness and their performance in a coupled system. In this study, we present a novel Fortran-Python interface designed to seamlessly integrate ML parameterizations into ESMs. This interface showcases high versatility by supporting popular ML frameworks like PyTorch, TensorFlow, and Scikit-learn. We demonstrate the interface's modularity and reusability through two cases: a ML trigger function for convection parameterization and a ML wildfire model. We conduct a comprehensive evaluation of memory usage and computational overhead resulting from the integration of Python codes into the Fortran ESMs. By leveraging this flexible interface, ML parameterizations can be effectively developed, tested, and integrated into ESMs.

54 ENVIRONMENTAL SCIENCES↗

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗

Method of authenticating integrated circuits using optical characteristics of physically unclonable functions

A method and apparatus for reading unique identifiers of an integrated circuit. The unique identifiers may be physically unclonable functions (PUFs), formed by high energy ions implanted into semiconductor material of the integrated circuit. The method may include electrically or optically stimulating each of the PUFs and sensing with an optical sensor optical characteristics of resulting light emitted from the PUFs. Then the method may include comparing values associated with the optical characteristics of the PUFs with groups of stored values in a circuit database. Each of the groups of stored values may be associated with optical characteristics of PUFs of a known authentic circuit. The method may then include the controller providing verification of authenticity of the integrated circuit when each of the values associated with the optical characteristics of the PUFs match the stored values of at least one of the groups in the circuit database.

42 ENGINEERING↗

INTEGRATION OF FLEX EQUIPMENT AND OPERATOR ACTIONS IN PLANT FORCE-ON-FORCE MODELS WITH DYNAMIC RISK ASSESSMENT

The overall operation and maintenance cost to protect nuclear power plants accounts for approximately 7% of the total cost of power generation, with labor accounting for half of this cost. In the current research, from interaction with utilities and other stakeholders, it was determined that physical security forces account for nearly 20% of the entire workforce at several nuclear power plants. Labor costs continue to rise in the U.S., so any measures to reduce the cost of operating a nuclear power plant will need to include a reduction in labor. The physical security pathway within the DOE’s Light Water Reactor Sustainability program aims to lower the cost of physical security through directed research into modeling and simulation, application of advanced sensors or deployment of advanced weapons. This report presents a modeling and simulation framework for integrating Diverse and Flexible Mitigation Capability (FLEX) portable equipment performance with Force on Force models of a plant’s physical security posture. The generic framework is described in detail, followed by a case study of modeling an adversarial attack aimed at causing a radiological release by sabotaging the plant’s power supply and its ultimate heat sink capabilities at a hypothetical nuclear power plant. Two different FLEX deployment strategies, series and parallel, are modeled with distinct timelines. The results of the adversarial attack modeled in a commercial Force on Force tool are integrated with the FLEX deployment model in INL’s dynamic modeling tool EMRALD. Monte Carlo simulation is used to model the distribution of the timeline in FLEX deployment strategies. The results demonstrate that, even in the extreme case of a successful adversarial attack, deployment of FLEX equipment can result in a significantly high likelihood of preventing radiological release. The modeling and simulation framework integrating FLEX equipment with Force on Force models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security.

97 MATHEMATICS AND COMPUTING↗

An Assessment of Impacts on Electric End-Use Load Profile of a Typical Residential Building from a Ground Source Heat Pump Systems Integrated with Underground Thermal Energy Storage

Buildings are the primary users of electricity, and a considerable amount of electric energy consumed in residential and commercial buildings is for meeting thermal demands (e.g., space heating/cooling, water heating, etc.). Therefore, integrating thermal energy storage with a building’s energy system has the potential to reshape the electric load profile of the building and mitigate the mismatch between the renewable power supply and the thermal demands of buildings. A novel ground-source heat pump system (GSHP) integrated with underground thermal energy storage (UTES) has been developed to enable flexible electric demand at buildings while meeting their thermal demands with an energy-efficient GSHP system. This paper assessed the potential impacts of the proposed system on the peak electrical demand and annual electricity consumption in a typical residential building at various climate zones in the United States. The results show that by replacing the conventional HVAC system (air conditioner and natural gas furnace) with the UTES integrated GSHP system, summer peak electric demand can be reduced by 27% to 50% depending on the weather. Although the winter peak electrical demand could be increased by up to 70% if the proposed system replaces natural gas furnace, the existing electricity supply system for single-family homes in most regions of the United States have adequate capacity for the replacement. While eliminating the natural gas consumption for space heating, the proposed system consumes more electricity on an annual basis in heating-dominated regions. However, it consumes less electricity in cooling dominated regions. The resulting variance in annual electricity consumption in various climate zones is within 15% compared with the baseline system.

Shi, Liang↗

Resilient buildings for fire-adapted landscapes: EE and flexible loads integrated with solar and storage microgrids

Energy efficiency (EE) and flexible loads can be part of a resilient buildings packagewhen they are combined with solar and storage in a clean energy microgrid to reduce the carbonfootprint of buildings and enable resilience to extreme events. Recent large wildfires and an emerging understanding of the fire-adapted landscapes in the American West underscore the urgency of work towards scaling and commercializing these systems. For example, a recent power shutoff in Northern California (October 9-12, 2019) resulted in 738k customers disconnected at the peak of the outage and is emblematic of challenges to come. Our paper reports on insights from a clean energy microgrid deployment pilot that integrates a 50 kWAC PV array, a 109 kW / 174 kWh battery system, switchgear to safely isolate from the regional power system, and communicating controllers for HVAC and refrigeration. The project will be commissioned in May 2020 and is sited at a critical infrastructure site in rural Northern California – in this case a gasoline station with convenience store. Our experience and results shed light on capabilities of integrated microgrids to provide value to customers during blue sky conditions and resilience during black sky days with high fire risk, and what opportunities and barriers exist for scaling these integrated microgrid systems in the near term. We use a simulation model to estimate how EE and flexibility can extend the run time of solar and storage, improving the reliability of power at critical sites.

Alstone, Peter↗

Life Cycle Assessment of Pasture-based Agrivoltaic Systems: Emissions and Energy Use of Integrated Rabbit Production

Agrivoltaic systems, which deliberately maximize the utility of a single parcel of land for both solar photovoltaic (PV) electricity production and agriculture, have been demonstrated as a viable technology that can ameliorate competing land uses and meet growing energy and food demands efficiently. The goal of this study is to assess the environmental impacts of a novel pasture-based agrivoltaic concept: co-farming rabbits and solar PV. A life cycle assessment (LCA) quantified the impacts of 1) the integrated agrivoltaic concept in comparison to conventional practices including 2) separate rabbit farming and PV production and 3) separate rabbit farming and conventional electricity production. The impact assessment methods employed to determine the environmental impacts were IPCC 2013 global warming potential 100a V1.03 and fossil energy demand V1.11. The results indicate that the pasture-based agrivoltaic system produces the least amount of greenhouse gas emissions (3.8 million kg CO₂ equivalent) and demands the least amount of fossil energy (46 million MJ) per functional unit of cumulative MWh output of electricity and cumulative kg of meat over 30 years in comparison to the two other scenarios under study. The pasture-based agrivoltaic system features a dual synergy that consequently produces 69.3% less emissions and demands 82.9% less fossil energy compared to non-integrated production. The potential for agrivoltaic systems to significantly reduce environmental impacts revealed by this LCA, which was limited to emissions and fossil energy demand, demonstrates that integrated solar and pasture-based agricultural systems are superior to conventional practices in terms of their comparatively lower emission and energy intensity. These findings provide empirical support for increased agrivoltaic system development more broadly.

14 SOLAR ENERGY↗

Probabilistic Look-ahead Contingency Analysis Integration with Commercial Tool and Practical Data

This paper presents an initial effort of integrating a smart sampling-based probabilistic look-ahead contingency analysis algorithm with General Electric (GE) Grid Solutions’ commercial energy management system (EMS) tool as a proof-of-concept for a seamless research tool integration using real world large-scale grid data. With the increasing impact of random forces such as variable generation and load, their stochastic behaviors cannot be ignored. However, the current practices are still dominated by deterministic tools. They are becoming increasingly inadequate for the future grid. The developed look-ahead contingency analysis algorithm incorporates forecast errors of variable energy and load to address the challenges brought by the increasing uncertainty of power system. The algorithm can reveal the potential violations caused by the variance of variable energy and load that are not normally detected by traditional deterministic approaches. To test its performance under practical environments ( real data with real commercial tool), significant efforts have been made to prepare test cases, modify GE EMS tool, and adapt an extreme value distribution algorithm to analyze the GE EMS’s violation-only outputs. The test results clearly demonstrate the effectiveness of the developed algorithm as new transformer violations that were not previously detected have been identified. This performance provides better situational awareness to engineers for their decision-making process under uncertainty. Moreover, with the discussion of computational performance and future work, this paper has shown a clear path for integrating the probabilistic algorithm with commercial tools to make us better equipped for the changing power system.

Modeling and simulation of power systems, constrai↗