Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Line Rating”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

HMI Data Driven Magnetohydrodynamic Model Predicted Active Region Photospheric Heating Rates: Their Scale Invariant, Flare Like Power Law Distributions, and Their Possible Association With Flares

A data driven, near photospheric, 3 D, non-force free magnetohydrodynamic model predicts time series of the complete current density, and the resistive heating rate Q at the photosphere in neutral line regions (NLRs) of 14 active regions (ARs). The model is driven by time series of the magnetic field B observed by the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory (SDO) satellite. Spurious Doppler periods due to SDO orbital motion are filtered out of the time series for B in every AR pixel. Errors in B due to these periods can be significant. The number of occurrences N(q) of values of Q > or = q for each AR time series is found to be a scale invariant power law distribution, N(Q) / Q−s, above an AR dependent threshold value of Q, where 0.3952 < or = s < or = 0.5298 with mean and standard deviation of 0.4678 and 0.0454, indicating little variation between ARs. Observations show that the number of occurrences N(E) of coronal flares with a total energy released > or = E obeys the same type of distribution, N(E) / E−S, above an AR dependent threshold value of E, with 0.38 < or approx. S < or approx. 0.60, also with little variation among ARs. Within error margins the ranges of s and S are nearly identical. This strong similarity between N(Q) and N(E) suggests a fundamental connection between the process that drives coronal flares and the process that drives photospheric NLR heating rates in ARs. In addition, results suggest it is plausible that spikes in Q, several orders of magnitude above background values, are correlated with times of the subsequent occurrence of M or X flares.

Goodman, Michael L.↗
DOE OSTI2015

Far-infrared and accretion luminosities of the present-day active galactic nuclei

We investigate the relation between star formation (SF) and black hole accretion luminosities, using a sample of 492 type-2 active galactic nuclei (AGNs) at z<0.22, which are detected in the far-infrared (FIR) surveys with AKARI and Herschel. We adopt FIR luminosities at 90 and 100 μm as SF luminosities, assuming the proposed linear proportionality of star formation rate with FIR luminosities. By estimating AGN luminosities from [O iii]λ5007 and [O i]λ6300 emission lines, we find a positive linear trend between FIR and AGN luminosities over a wide dynamical range. This result appears to be inconsistent with the recent reports that low-luminosity AGNs show essentially no correlation between FIR and X-ray luminosities, while the discrepancy is likely due to the Malmquist and sample selection biases. By analyzing the spectral energy distribution, we find that pure-AGN candidates, of which FIR radiation is thought to be AGN-dominated, show significantly low-SF activities. These AGNs hosted by low-SF galaxies are rare in our sample (∼1%). However, the low fraction of low-SF AGNs is possibly due to observational limitations since the recent FIR surveys are insufficient to examine the population of high-luminosity AGNs hosted by low-SF galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Study of the light scalar a 0 ( 980 ) through the decay D 0 → a 0 ( 980 ) − e + ν e with a 0 ( 980 ) − → η π −

Using 7.93 fb − 1 of e + e − collision data collected at a center-of-mass energy of 3.773 GeV with the BESIII detector, we present an analysis of the decay D 0 → η π − e + ν e . The branching fraction of the decay D 0 → a 0 ( 980 ) − e + ν e with a 0 ( 980 ) − → η π − is measured to be ( 0.86 ± 0.1 7 stat ± 0.0 5 syst ) × 10 − 4 . The decay dynamics of this process is studied with a single-pole parametrization of the hadronic form factor and the Flatté formula describing the a 0 ( 980 ) line shape in the differential decay rate. The product of the form factor f + a 0 ( 0 ) and the Cabibbo-Kobayashi-Maskawa matrix element | V c d | is determined for the first time with the result f + a 0 ( 0 ) | V c d | = 0.126 ± 0.01 3 stat ± 0.00 3 syst . Published by the American Physical Society 2025

Ablikim, M.↗
NASA NTRS1985

Multiline Transfer and the Dynamics of Stellar Winds

A Monte Carlo technique for treating multiline transfer in stellar winds is described. With a line list containing many thousands of transitions and with fairly realistic treatments of ionization, excitation and line formation, the resulting code allows the dynamic effects of overlapping lines the investigation of and provides the means to directly synthesize the complete spectrum of a star and its wind. It is found that the computed mass loss rate for data Puppis agrees with the observed rate. The synthesized spectrum of zeta Puppis also agrees with observational data. This confirms that line driving is the dominant acceleration mechanism in this star's wind.

Abbott, D. C.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
DOE OSTI2024

Quarkonium transport in weakly and strongly coupled plasmas

We report on progress in the nonperturbative understanding of quarkonium dynamics inside a thermal plasma. The time evolution of small-size quarkonium is governed by two-point correlation functions of chromoelectric fields dressed with an adjoint Wilson line, known in this context as generalized gluon distributions (GGDs). The GGDs have been calculated in both weakly and strongly coupled plasmas by using perturbative and holographic methods. Strikingly, the results of our calculations for a strongly coupled plasma indicate that the quarkonium dissociation and recombination rates vanish in the transport descriptions that assume quarkonium undergoes Markovian dynamics. However, this does not imply that the dynamics is trivial. As a starting point to explore the phenomenological consequences of the result at strong coupling, we show a calculation of the Y(1 S ) formation probability in time-dependent perturbation theory. This is a first step towards the development of a transport formalism that includes non-Markovian effects, which, depending on how close the as of yet undetermined nonperturbative QCD result of the GGDs is to the strongly coupled N = 4 SYM result, could very well dominate over the Markovian ones in quark-gluon plasma produced at RHIC and the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Stellarator microinstabilities and turbulence at low magnetic shear

Gyrokinetic simulations of drift waves in low-magnetic-shear stellarators reveal that simulation domains comprised of multiple turns can be required to properly resolve critical mode structures important in saturation dynamics. Marginally stable eigenmodes important in saturation of ion temperature gradient modes and trapped electron modes in the Helically Symmetric Experiment (HSX) stellarator are observed to have two scales, with the envelope scale determined by the properties of the local magnetic shear and an inner scale determined by the interplay between the local shear and magnetic field-line curvature. Properly resolving these modes removes spurious growth rates that arise for extended modes in zero-magnetic-shear approximations, enabling use of a zero-magnetic-shear technique with smaller simulation domains and attendant cost savings. Analysis of subdominant modes in trapped electron mode (TEM)-driven turbulence reveals that the extended marginally stable modes play an important role in the nonlinear dynamics, and suggests that the properties induced by low magnetic shear may be exploited to provide another route for turbulence saturation.

Faber, B. J. (ORCID:000000034934400X)↗
NASA NTRS1989

Experimental study on two-dimensional free-flying robot satellite model

The experimental study on a control method for a free flying space robotic arm was treated by means of a two-dimensional laboratory model. The main target is to develop a new control method for trajectory tracking or target capturing, considering dynamical interaction between the manipulator arm and the base vehicle in space micro-gravity environment. In order to simulate the micro-gravity environment mechanically, a laboratory model of a robot satellite supported on air bearings was developed. The model comprises a base equipped with power and air supplies and a two-link manipulator arm. This model has relatively low gravitational and frictional disturbance in planar motion. An on-line resolved motion rate control scheme with vision feedback is developed for experimenting capture operations. The scheme utilizes the Generalized Jacobian Matrix. In experiment, the acceleration environment of the model is evaluated firstly, then target capture operations are examined. The manipulator can properly chase and capture both a standing target and a moving target in spite of complex satellite/manipulator dynamical interactions. The experimental results confirm the validity of the Generalized Jacobian Matrix concept and the proposed control method.

Umetani, Yoji↗
NASA NTRS1999

Cardiac interbeat interval dynamics from childhood to senescence : comparison of conventional and new measures based on fractals and chaos theory

BACKGROUND: New methods of R-R interval variability based on fractal scaling and nonlinear dynamics ("chaos theory") may give new insights into heart rate dynamics. The aims of this study were to (1) systematically characterize and quantify the effects of aging from early childhood to advanced age on 24-hour heart rate dynamics in healthy subjects; (2) compare age-related changes in conventional time- and frequency-domain measures with changes in newly derived measures based on fractal scaling and complexity (chaos) theory; and (3) further test the hypothesis that there is loss of complexity and altered fractal scaling of heart rate dynamics with advanced age. METHODS AND RESULTS: The relationship between age and cardiac interbeat (R-R) interval dynamics from childhood to senescence was studied in 114 healthy subjects (age range, 1 to 82 years) by measurement of the slope, beta, of the power-law regression line (log power-log frequency) of R-R interval variability (10(-4) to 10(-2) Hz), approximate entropy (ApEn), short-term (alpha(1)) and intermediate-term (alpha(2)) fractal scaling exponents obtained by detrended fluctuation analysis, and traditional time- and frequency-domain measures from 24-hour ECG recordings. Compared with young adults (<40 years old, n=29), children (<15 years old, n=27) showed similar complexity (ApEn) and fractal correlation properties (alpha(1), alpha(2), beta) of R-R interval dynamics despite lower spectral and time-domain measures. Progressive loss of complexity (decreased ApEn, r=-0.69, P<0.001) and alterations of long-term fractal-like heart rate behavior (increased alpha(2), r=0.63, decreased beta, r=-0.60, P<0.001 for both) were observed thereafter from middle age (40 to 60 years, n=29) to old age (>60 years, n=29). CONCLUSIONS: Cardiac interbeat interval dynamics change markedly from childhood to old age in healthy subjects. Children show complexity and fractal correlation properties of R-R interval time series comparable to those of young adults, despite lower overall heart rate variability. Healthy aging is associated with R-R interval dynamics showing higher regularity and altered fractal scaling consistent with a loss of complex variability.

NASA Discipline Cardiopulmonary↗
NASA NTRS1985

Pitch change thresholds as influenced by angular rate and duration of the apparent horizon's displacement

The displacement threshold (DT) during voluntary visual tracking of a line and during fixation of a stable spot located at the initial stimulus position was determined to help in the modeling pilot manual dynamics during the nose-up 'flare' maneuver. Forty-six observers made paired comparison, forced choice judgments of the maximal downward displacements of a horizontal line moving at various speeds for the durations of 0.25, 0.50, and 1.00 sec. The results indicated that the DT is significantly lower, the confidence is higher, and the number of errors and guesses is smaller when a fixation dot is present; moreover, increasing both the stimulus velocity and duration significantly improves the accuracy of displacement discrimination and increases mean confidence.

Haines, R. F.↗
DOE OSTI2014

ALMA imaging of gas and dust in a galaxy protocluster at redshift 5.3: [C II] emission in 'typical' galaxies and dusty starbursts ≈1 billion years after the big bang

We report interferometric imaging of [C II]({sup 2} P {sub 3/2}→{sup 2} P {sub 1/2}) and OH({sup 2}Π{sub 1/2} J = 3/2→1/2) emission toward the center of the galaxy protocluster associated with the z = 5.3 submillimeter galaxy (SMG) AzTEC-3, using the Atacama Large (sub)Millimeter Array (ALMA). We detect strong [C II], OH, and rest-frame 157.7 μm continuum emission toward the SMG. The [C II]({sup 2} P {sub 3/2}→{sup 2} P {sub 1/2}) emission is distributed over a scale of 3.9 kpc, implying a dynamical mass of 9.7 × 10{sup 10} M {sub ☉}, and a star formation rate (SFR) surface density of Σ{sub SFR} = 530 M {sub ☉} yr{sup –1} kpc{sup –2}. This suggests that AzTEC-3 forms stars at Σ{sub SFR} approaching the Eddington limit for radiation pressure supported disks. We find that the OH emission is slightly blueshifted relative to the [C II] line, which may indicate a molecular outflow associated with the peak phase of the starburst. We also detect and dynamically resolve [C II]({sup 2} P {sub 3/2}→{sup 2} P {sub 1/2}) emission over a scale of 7.5 kpc toward a triplet of Lyman-break galaxies with moderate UV-based SFRs in the protocluster at ∼95 kpc projected distance from the SMG. These galaxies are not detected in the continuum, suggesting far-infrared SFRs of <18-54 M {sub ☉} yr{sup –1}, consistent with a UV-based estimate of 22 M {sub ☉} yr{sup –1}. The spectral energy distribution of these galaxies is inconsistent with nearby spiral and starburst galaxies, but resembles those of dwarf galaxies. This is consistent with expectations for young starbursts without significant older stellar populations. This suggests that these galaxies are significantly metal-enriched, but not heavily dust-obscured, 'normal' star-forming galaxies at z > 5, showing that ALMA can detect the interstellar medium in 'typical' galaxies in the very early universe.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

An Assessment of Helium Evolution from Helium-Saturated Propellant Depressurization in Space

Helium evolution from the transfer of helium-saturated propellant in space is quantified to assess its impacts from creating two-phase gas/liquid flow from the supply tank, gas injection into the receiving tank, and liquid discharge from the receiving tank. Propellant transfer takes place between two similar tanks whose maximum storage capacity is approximately 2.55 cubic meters each. The maximum on-orbit propellants transfer capability is 9000 lbm (fuel and oxidizer). The transfer line is approximately 1.27 cm in diameter and 6096 cm in length and comprised of the fluid interconnect system (FICS), the orbiter propellant transfer system (OPTS), and the International Space Station (ISS) propulsion module (ISSPM). The propellant transfer rate begins at approximately 11 liter per minute (lpm) and subsequently drops to approximately 0.5 lpm. The tank nominal operating pressure is approximately 1827 kPa (absolute). The line pressure drops for Monomethy1hydrazine (MMH) and Nitrogen tetroxide (NTO) at 11.3 lpm are approximately 202 kPa and 302 kPa, respectively. The pressure-drop results are based on a single-phase flow. The receiving tank is required to vent from approximately 1827 kPa to a lower pressure to affect propellant transfer. These pressure-drop scenarios cause the helium-saturated propellants to release excess helium. For tank ullage venting, the maximum volumes of helium evolved at tank pressure are approximately 0.5 ft3 for MMH and 2 ft3 for NTO. In microgravity environment, due to lack of body force, the helium evolution from a liquid body acts to propel it, which influences its fluid dynamics. For propellant transfer, the volume fractions of helium evolved at line pressure are 0.1% by volume for MMH and 0.6 % by volume for NTO at 11.3 lpm. The void fraction of helium evolved varies as an approximate second order power function of flow rate.

Nguyen, Bich N.↗
NASA NTRS1985

On the use of the cross section concept as applied to pulsed CO2 laser dynamics

The 'cross sections' which are commonly used in combination with the molecular vibrational level densities to describe induced transition rates are time dependent in a pulsed discharge. This greatly affects the relative rates of buildup of radiation at the various lines of a molecular gain medium in a nondispersive cavity, and the consequences for line selection by injection in a TEA-CO2 laser are discussed.

Flamant, P. H.↗
DOE OSTI2016

Mechanical behavior and microstructure of Ti{sub 20}Hf{sub 20}Zr{sub 20}Ta{sub 20}Nb{sub 20} high-entropy alloy loaded under quasi-static and dynamic compression conditions

The microstructure and the mechanical behavior of equimolar Ti{sub 20}Hf{sub 20}Zr{sub 20}Ta{sub 20}Nb{sub 20} high-entropy alloy in a wide range of initial strain rates between ~ 10{sup −3} s{sup −1} and ~ 3.4 × 10{sup 3} s{sup −1} were studied. A significant increment in the yield strength with increasing strain rate was observed. The yield strength at ~ 3.4 × 10{sup 3} s{sup −1} was about 40% higher than that measured at ~ 10{sup −3} s{sup −1}. Analysis by electron backscatter diffraction shows that in the low strain rate regime (up to ~ 10 s{sup −1}) the deformation occurs mainly in evenly distributed bands, while in the dynamic regime the deformation is strongly localized in macroscopic shear bands accompanied by softening even after the onset of yielding. The Kernel Average Misorientation technique reveals a high level of lattice rotation within these bands that also carries intense shear. In addition, X-ray diffraction line profile analysis indicates that the sharp increase in the flow stress is mostly related to an increase of the dislocation density. - Highlights: • Strain rate effect on the plastic behavior of Ti{sub 20}Hf{sub 20}Zr{sub 20}Ta{sub 20}Nb{sub 20} HEA was studied. • Low strain rate regime was characterized by a continuous hardening. • At high strain rates softening occurred shortly after the onset of yielding. • Intense strain localization in shear bands occurred in the high strain rate regime. • High dislocation density explained the upturn of flow stress at high strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Electromagnetic effects on dynamics of high-beta filamentary structures

The impacts of the electromagnetic effects on blob dynamics are considered. Electromagnetic BOUT++ simulations on seeded high-beta blobs demonstrate that inhomogeneity of magnetic curvature or plasma pressure along the filament leads to bending of the blob filaments and the magnetic field lines due to increased propagation time of plasma current (Alfvén time). The bending motion can enhance heat exchange between the plasma facing materials and the inner scrape-off layer (SOL) region. The effects of sheath boundary conditions on the part of the blob away from the boundary are also diminished by the increased Alfvén time. Using linear analysis and BOUT++ simulations, it is found that electromagnetic effects in high temperature and high density plasmas reduce the growth rate of resistive drift wave instability when resistivity drops below a certain value. The blobs temperature decreases in the course of its motion through the SOL and so the blob can switch from the electromagnetic to the electrostatic regime where resistive drift waves become important again.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Enhanced visible and near-infrared capabilities of the JET mirror-linked divertor spectroscopy system

The mirror-linked divertor spectroscopy diagnostic on JET has been upgraded with a new visible and near-infrared grating and filtered spectroscopy system. New capabilities include extended near-infrared coverage up to 1875 nm, capturing the hydrogen Paschen series, as well as a 2 kHz frame rate filtered imaging camera system for fast measurements of impurity (Be II) and deuterium Dα, Dβ, Dγ line emission in the outer divertor. The expanded system provides unique capabilities for studying spatially resolved divertor plasma dynamics at near-ELM resolved timescales as well as a test bed for feasibility assessment of near-infrared spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Solid-State Fault Current Limiter Development: Design and Testing Update of a 15kV SSCL Power Stack

The Solid-State Fault Current Limiter (SSCL) is a promising technology that can be applied to utility power delivery systems to address the problem of increasing fault currents associated with load growth. As demand continues to grow, more power is added to utility system either by increasing generator capacity or by adding distributed generators, resulting in higher available fault currents, often beyond the capabilities of the present infrastructure. The SSCL is power-electronics based equipment designed to work with the present utility system to address this problem. The SSCL monitors the line current and dynamically inserts additional impedance into the line in the event of a fault being detected. The SSCL is based on a modular design and can be configured for 5kV through 69kV systems at nominal current ratings of 1000A to 4000A. Results and Findings This report provides the final test results on the development of 15kV class SSCL single phase power stack. The scope of work included the design of the modular standard building block sub-assemblies, the design and manufacture of the power stack and the testing of the power stack for the key functional tests of continuous current capability and fault current limiting action. Challenges and Objectives Solid-State Current Limiter technology impacts a wide spectrum of utility engineering and operating personnel. It addresses the problems associated with load growth both at Transmission and Distribution class networks. The design concept is pioneering in terms of developing the most efficient and compact power electronics equipment for utility use. The initial test results of the standard building blocks are promising. The independent laboratory tests of the power stack are promising. However the complete 3 phase system needs rigorous testing for performance and reliability. Applications, Values, and Use The SSCL is an intelligent power-electronics device which is modular in design and can provide current limiting or current interrupting capabilities. It can be applied to variety of applications from distribution class to transmission class power delivery grids and networks. It can also be applied to single major commercial and industrial loads and distributed generator supplies. The active switching of devices can be further utilized for protection of substation transformers. The stress on the system can be reduced substantially improving the life of the power system. It minimizes the voltage sag by speedy elimination of heavy fault currents and promises to be an important element of the utility power system. DOE Perspective This development effort is now focused on a 15kV system. This project will help mitigate the challenges of increasing available fault current. DOE has made a major contribution in providing a cost effective SSCL designed to integrate seamlessly into the Transmission and Distribution networks of today and the future. Approach SSCL development program for a 69kV SSCL was initiated which included the use of the Super GTO advanced semiconductor device which won the 2007 R&D100 Award. In the beginning, steps were identified to accomplish the economically viable design of a 69kV class Solid State Current Limiter that is extremely reliable, cost effective, and compact enough to be applied in urban transmission. The prime thrust in design and development was to encompass the 1000A and the 3000A ratings and provide a modular design to cover the wide range of applications. The focus of the project was then shifted to a 15kV class SSCL. The specifications for the 15kV power stack are reviewed. The design changes integrated into the 15kV power stack are discussed. In this Technical Update the complete project is summarized followed by a detailed test report. The power stack independent high voltage laboratory test requirements and results are presented. Keywords Solid State Current Limiter, SSCL, Fault Current Limiter, Fault Current Controller, Power electronics controller, Intelligent power-electronics Device, IED

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A Kinetic Monte Carlo Approach for Simulating Cascading Transmission Line Failure

In this work, cascading transmission line failures are studied through a dynamical model of the power system operating under fixed conditions. The power grid is modeled as a stochastic dynamical system where first-principles electromechanical dynamics are excited by small Gaussian disturbances in demand and generation around a specified operating point. In this context, a single line failure is interpreted in a large deviation context as a first escape event across a surface in phase space defined by line security constraints. The resulting system of stochastic differential equations admits a transverse decomposition of the drift, which leads to considerable simplification in evaluating the quasipotential (rate function) and, consequently, computation of exit rates. Tractable expressions for the rate of transmission line failure in a restricted network are derived from large deviation theory arguments and validated against numerical simulations. Extensions to realistic settings are considered, and individual line failure models are aggregated into a Markov model of cascading failure inspired by chemical kinetics. Cascades are generated by traversing a graph composed of weighted edges representing transitions to degraded network topologies. Numerical results indicate that the Markov model can produce cascades with qualitative power-law properties similar to those observed in empirical cascades.

Kinetic Monte Carlo, rare events, Cascading Failur↗