Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large deviations”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Elucid—exploring the local universe with the reconstructed initial density field. I. Hamiltonian Markov chain Monte Carlo method with particle mesh dynamics

Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input nonlinear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle-mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. These constrained simulations accurately reproduce both the amplitudes and phases of the input simulations at various z. Using a PM model with a grid cell size of 0.75 h {sup –1} Mpc and 40 time steps in the HMC can recover more than half of the phase information down to a scale k ∼ 0.85 h Mpc{sup –1} at high z and to k ∼ 3.4 h Mpc{sup –1} at z = 0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method is limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Effect of noncircularity of experimental beam on CMB parameter estimation

Measurement of Cosmic Microwave Background (CMB) anisotropies has been playing a lead role in precision cosmology by providing some of the tightest constrains on cosmological models and parameters. However, precision can only be meaningful when all major systematic effects are taken into account. Non-circular beams in CMB experiments can cause large systematic deviation in the angular power spectrum, not only by modifying the measurement at a given multipole, but also introducing coupling between different multipoles through a deterministic bias matrix. Here we add a mechanism for emulating the effect of a full bias matrix to the PLANCK likelihood code through the parameter estimation code SCoPE. We show that if the angular power spectrum was measured with a non-circular beam, the assumption of circular Gaussian beam or considering only the diagonal part of the bias matrix can lead to huge error in parameter estimation. We demonstrate that, at least for elliptical Gaussian beams, use of scalar beam window functions obtained via Monte Carlo simulations starting from a fiducial spectrum, as implemented in PLANCK analyses for example, leads to only few percent of sigma deviation of the best-fit parameters. However, we notice more significant differences in the posterior distributions for some of the parameters, which would in turn lead to incorrect errorbars. These differences can be reduced, so that the errorbars match within few percent, by adding an iterative reanalysis step, where the beam window function would be recomputed using the best-fit spectrum estimated in the first step.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Probing the electroweak phase transition at the LHC

We study the correlation between the value of the triple Higgs coupling and the nature of the electroweak phase transition. We use an effective potential approach, including higher order, non-renormalizable terms coming from integrating out new physics. We show that if only the dimension six operators are considered, large positive deviations of the triple Higgs coupling from its Standard Model (SM) value are predicted in the regions of parameter space consistent with a strong first order electroweak phase transition (SFOEPT). We also show that at higher orders sizable and negative deviations of the triple Higgs coupling may be obtained, and the sign of the corrections tends to be correlated with the order of the phase transition. We also consider a singlet extension of the SM, which allows us to establish the connection with the effective field theory (EFT) approach and analyze the limits of its validity. Furthermore, we study how to probe the triple Higgs coupling from the double Higgs production at the LHC. Lastly, we show that selective cuts in the invariant mass of the two Higgs bosons should be used, to maximize the sensitivity for values of the triple Higgs coupling significantly different from the Standard Model one.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Angular Momentum Removal by Neutron and γ-Ray Emissions during Fission Fragment Decays

In this work, we investigate the angular momentum removal from fission fragments (FFs) through neutron and γ-ray emission, finding that about half the neutrons are emitted with angular momenta ≥ 1.5 ℏ and that the change in angular momentum after the emission of neutrons and statistical γ rays is significant, contradicting usual assumptions. Per fission event, in our simulations, the neutron and statistical γ-ray emissions change the spin of the fragment by 3.5 – 5 ℏ , with a large standard deviation comparable to the average value. Such wide angular momentum removal distributions can hide any underlying correlations in the fission fragment initial spin values. Within our model, we reproduce data on spin measurements from discrete transitions after neutron emissions, especially in the case of light FFs. The agreement further improves for the heavy fragments if one removes from the analysis the events that would produce isomeric states. Finally, we show that while in our model the initial FF spins do not follow a sawtoothlike behavior observed in recent measurements, the average FF spin computed after neutron and statistical γ emissions exhibits a shape that resembles a sawtooth. This suggests that the average FF spin measured after statistical emissions is not necessarily connected with the scission mechanism as previously implied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Observation of Medium-Induced Yield Enhancement and Acoplanarity Broadening of Low-𝑝T Jets from Measurements in 𝑝⁢𝑝 and Central Pb-Pb Collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The ALICE Collaboration reports the measurement of semi-inclusive distributions of charged-particle jets recoiling from a high transverse momentum (high p T ) hadron trigger in proton-proton and central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. A data-driven statistical method is used to mitigate the large uncorrelated background in central Pb-Pb collisions. Recoil jet distributions are reported for jet resolution parameter R = 0.2, 0.4, and 0.5 in the range 7 < p T.jet < 140 GeV=c and trigger-recoil jet azimuthal separation π/2 < Δφ < π. The measurements exhibit a marked medium-induced jet yield enhancement at low p T and at large azimuthal deviation from Δφ ~ π. The enhancement is characterized by its dependence on Δφ, which has a slope that differs from zero by 4.7σ. Comparisons to model calculations incorporating different formulations of jet quenching are reported. These comparisons indicate that the observed yield enhancement arises from the response of the QGP medium to jet propagation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

SU-E-T-458: Determining Threshold-Of-Failure for Dead Pixel Rows in EPID-Based Dosimetry

Purpose: A pixel correction map is applied to all EPID-based applications on the TrueBeam (Varian Medical Systems, Palo Alto, CA). When dead pixels are detected, an interpolative smoothing algorithm is applied using neighboring-pixel information to supplement missing-pixel information. The vendor suggests that when the number of dead pixels exceeds 70,000, the panel should be replaced. It is common for entire detector rows to be dead, as well as their neighboring rows. Approximately 70 rows can be dead before the panel reaches this threshold. This study determines the number of neighboring dead-pixel rows that would create a large enough deviation in measured fluence to cause failures in portal dosimetry (PD). Methods: Four clinical two-arc VMAT plans were generated using Eclipse's AXB algorithm and PD plans were created using the PDIP algorithm. These plans were chosen to represent those commonly encountered in the clinic: prostate, lung, abdomen, and neck treatments. During each iteration of this study, an increasing number of dead-pixel rows are artificially applied to the correction map and a fluence QA is performed using the EPID (corrected with this map). To provide a worst-case-scenario, the dead-pixel rows are chosen so that they present artifacts in the highfluence region of the field. Results: For all eight arc-fields deemed acceptable via a 3%/3mm gamma analysis (pass rate greater than 99%), VMAT QA yielded identical results with a 5 pixel-width dead zone. When 10 dead lines were present, half of the fields had pass rates below the 99% pass rate. With increasing dead rows, the pass rates were reduced substantially. Conclusion: While the vendor still suggests to request service at the point where 70,000 dead rows are measured (as recommended by the vendor), the authors suggest that service should be requested when there are greater than 5 consecutive dead rows.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Anomalous dips in reflection spectra of optical polymers deposited on plasmonic metals

Abstract We have studied reflection spectra of dye-doped and undoped polymers deposited onto Ag and Au substrates and found anomalous dips in the UV spectral range. On top of Ag substrates, the λ ∼ 375 nm dips were observed in undoped PMMA, PVP, and PS polymers as well as PMMA doped with Rh590 and HITC laser dyes. In silver-based samples, the spectral positions of the observed reflection dips were close to singularities in the refractive indexes of surface plasmon polaritons (SPPs) propagating at the interface between Ag and polymer. The latter singularities can tentatively explain the λ ∼ 375 nm reflection dips, if the scattering of Ag and polymeric films is large enough to launch SPP without any prism or grating. The dips observed in reflection of Rh590:PMMA and HITC:PMMA on top of Au, were more pronounced than those on Ag, broader, shifted to shorter wavelengths, and their spectral positions had large standard deviations. Furthermore, no anomalous dips in gold-based samples were observed in the reflection spectra of undoped PMMA, PVP, and PS polymers, and a modest singularity in the SPP refractive index, predicted theoretically at λ ∼ 500 nm, cannot explain the dips in the UV reflection spectra observed experimentally. It appears likely that the origin of the reflection dips on top of Au substrates is different from that on top of Ag substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimal Distribution System Voltage Regulation using State Estimation and DER Grid-Support Functions

Increasing solar energy penetrations may create challenges for distribution system operations because production variability can lead to large voltage deviations or protection system miscoordination. Instituting advanced management systems on distribution systems is one promising method for combating these challenges by intelligently controlling distribution assets to regulate voltage and ensure protection safety margins. While it is generally not the case today, greater deployment of power system sensors and interoperable distributed energy resources (DER)e.g., photovoltaic (PV) inverters, energy storage systems (ESS), electric vehicles (EVs)will enable situational awareness, control, and optimization of distribution systems. In this work, a control system was created which measures power system parameters to estimate the status of a feeder, forecasts the distribution state over a short-term horizon, and issues optimal set point commands to distribution-connected equipment to regulate voltage and protect the system. This two-year project integrated multiple research innovations into a management system designed to safely allow PV penetrations of 50% or greater. The integrated software was demonstrated through extensive real-time (RT) and power hardware-in-the-loop studies and a field demonstration on a live power system with a 684 kVA PV system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Intermediate-Term Thermal Aging Effect Evaluation for Grade 92 and 316L at The LWR Relevant Temperature

Life extension of the existing nuclear reactors imposes accumulated damages, such as higher fluences and longer period of corrosion, to structural materials, which would result in significant challenges to the traditional reactor materials such as type 304 and 316 stainless steels. Advanced alloys with superior radiation resistance will increase safety margins, design flexibility, and economics for not only the life extension of the existing fleet but also new builds with advanced reactor designs. The Electric Power Research Institute (EPRI) teamed up with Department of Energy (DOE) on the Advanced Radiation Resistant Materials (ARRM) program, aiming to develop and test degradation resistant alloys from current commercial alloy specifications by 2021 to a new advanced alloy with superior degradation resistance in light water reactor (LWR)-relevant environments by 2024. Based on a comprehensive microstructure and property screening performed in Phase-1 of the ARRM program, a total of five alloys, together with 316L and X-750 as references, were down-selected for Phase-2 neutron irradiation studies. Because thermal aging could exert a synergistic effect on neutron irradiation due to the low neutron damage rate on the order of 10 –7 displacements per atom per second (dpa/s), Grade 92 (one of the five down-selected alloys) and two heats of 316L were selected in this task to study the effect of aging at 350°C for 12.5–12.7 kh on microstructure and mechanical properties. The aging time is approximately corresponding to 5 dpa neutron irradiation, which can be served as a reference for the 5-dpa-irradiated samples of the alloys to help understand the independent neutron damage influence on microstructure and mechanical properties. Optical microscopy and scanning electron microscopy were used for microstructural characterization. Hardness, tensile, Charpy impact toughness, and fracture toughness in the ductile regime of the aged samples were examined. The aging of Grade 92 led to the formation of many Laves phase in sizes of ~100–200 nm and resulted in some reduction in hardness and yield/ultimate tensile strength with some increases in uniform and total plastic elongations, which may have helped the enhancement of Charpy impact toughness, e.g., ~4 J increase in upper-shelf energy and 20.5°C reduction in ductile-brittle transition temperature compared with the unaged condition. The fracture toughness of the aged Grade 92 showed decent toughness of ~303 MPa√m (K Jq ) with ~87 tearing modulus at 22°C, which reduced to ~241 MPa√m (~20% reduction) with ~75 tearing modulus at 300°C. The two heats of 316L are differentiated by their amounts of δ-ferrite, i.e., ~1 vol% in heat T1103 and ~4 vol% in heat N5B8. The two heats of 316L were subjected to 15% cold work (CW) as this condition is often used in nuclear reactors to help trapping radiation-induced defects and thus postpone the steady state swelling stage. The aging did not result in any noticeable microstructure changes, except for possible segregations at grain boundaries that need to be further investigated. The aging resulted in slight reduction in hardness of 316L-T1103 and slight increase in hardness of 316L-N5B8, which is consistent with their tensile testing results. The 15%CW might have introduced some inhomogeneity, resulting in large standard deviations in hardness, which were reduced after the aging. The aged 316L followed the same trend as the unaged condition for the strength and elongation results from the tensile testing. Minimum elongations appeared within ~300–600°C, above which the elongations seem to increase. Unlike the negligible or minor changes in microstructure, hardness and tensile results of the two heats of 316L, their fracture toughness showed noticeable difference. The aged 316L-T1103 showed good toughness of ~350 MPa√m with ~80 tearing modulus at 22°C, which reduced to ~271 MPa√m (~23% reduction) with ~83 tearing modulus at 300°C. In contrast, the aged 316L-N5B8 showed lower toughness of ~244 MPa√m with ~53 tearing modulus at 22°C, which reduced to ~180 MPa√m (~26% reduction) with ~31 tearing modulus at 300°C. The preliminary results indicate that the presence of high volume of δ-ferrite would noticeable impair fracture toughness although it may not influence hardness and tensile properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Open Source Fault-tolerant Grid Frequency Measurement for Solar Inverters

The Discrete Fourier transform (DFT) based measurement algorithms are one of the most common measurement algorithms for grid parameter estimation such as rms, phase angle, frequency. Over the past few years, many DFT based algorithms have been developed to enhance its measurement accuracy under steady-state and/or dynamic grid conditions. For example, an adaptive band-pass filter utilizing exponential modulation filter has been proposed to reduce measurement errors at the presence of large frequency deviations. Measurement accuracy of different algorithms including FIR filter, extended Kalman filtering (EKF), and enhanced DFT method have been compared in detail under different grid conditions. Two artificial signals that have 90-degree phase difference were constructed by the Clarke transformation to address the frequency spectrum leakage of DFT. A multi-module approach was developed to enhance both steady-state and dynamic measurement accuracies, in which each module was developed to eliminate some specific errors. Besides DFT-based measurement algorithms, some signal model-based algorithms have been developed to further improve the accuracy under dynamic conditions. However, a key drawback of the state-of-the-art algorithms is that they cannot perform measurements accurately during system transient faults. In the Blue Cut Fire event, there was a phase angle jump of about 26 degrees in the voltage waveform during the transient fault. The phase angle jump fault will cause waveform discontinuity, and these algorithms will fail to provide reliable measurements during this period because they typically assume the waveform to be measured is continuous, no matter what method (DFT, PLL, EKF, FIR, or Taylor WLS) is used for estimation. In fact, the measurement errors during the system transient faults like phase-jump is not required in the IEEE Standard. As a result, although a measurement instrument can pass the strict IEEE Standard, it could still be the source of the problem in the future if we have similar system transient faults, which could happen again. Therefore, developing the fault-tolerant measurement technology is the key to solve the problem.

14 SOLAR ENERGY↗
DOE OSTI2022

Residuals-based distributionally robust optimization with covariate information

We consider data-driven approaches that integrate a machine learning prediction model within distributionally robust optimization (DRO) given limited joint observations of uncertain parameters and covariates. Our framework is flexible in the sense that it can accommodate a variety of regression setups and DRO ambiguity sets. We investigate asymptotic and finite sample properties of solutions obtained using Wasserstein, sample robust optimization, and phi-divergence-based ambiguity sets within our DRO formulations, and explore cross-validation approaches for sizing these ambiguity sets. Through numerical experiments, we validate our theoretical results, study the effectiveness of our approaches for sizing ambiguity sets, and illustrate the benefits of our DRO formulations in the limited data regime even when the prediction model is misspecified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Signatures of Horndeski gravity on the dark matter bispectrum

We present a detailed study of second-order matter perturbations for the general Horndeski class of models. Being the most general scalar-tensor theory having second-order equations of motion, it includes many known gravity and dark energy theories and General Relativity with a cosmological constant as a specific case. This enables us to estimate the leading order dark matter bispectrum generated at late-times by gravitational instability. We parametrize the evolution of the first and second-order equations of motion as proposed by Bellini and Sawicki (2014), where the free functions of the theory are assumed to be proportional to the dark energy density. We show that it is unnatural to have large ∼> 10% (∼> 1%) deviations of the bispectrum introducing even larger ∼ 30% (∼ 5%) deviations in the linear growth rate. Considering that measurements of the linear growth rate have much higher signal-to-noise than bispectrum measurements, this indicates that for Horndeski models which reproduce the expansion history and the linear growth rate as predicted by GR the dark matter bispectrum kernel can be effectively modelled as the standard GR one. On the other hand, an observation of a large bispectrum deviation that can not be explained in terms of bias would imply either that the evolution of perturbations is strongly different than the evolution predicted by GR or that the theory of gravity is exotic (e.g., breaks the weak equivalence principle) and/or fine-tuned.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Rapidly rotating neutron stars with a massive scalar field—structure and universal relations

We construct rapidly rotating neutron star models in scalar-tensor theories with a massive scalar field. The fact that the scalar field has nonzero mass leads to very interesting results since the allowed range of values of the coupling parameters is significantly broadened. Deviations from pure general relativity can be very large for values of the parameters that are in agreement with the observations. We found that the rapid rotation can magnify the differences several times compared to the static case. The universal relations between the normalized moment of inertia and quadrupole moment are also investigated both for the slowly and rapidly rotating cases. The results show that these relations are still EOS independent up to a large extend and the deviations from pure general relativity can be large. This places the massive scalar-tensor theories amongst the few alternative theories of gravity that can be tested via the universal I -Love- Q relations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tracer Particles for Core-collapse Supernova Nucleosynthesis: The Advantages of Moving Backward

Abstract After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multidimensional simulations. However, obtaining nucleosynthesis yields based on such multidimensional core-collapse supernova simulations is not straightforward. Frequently, tracer particles are employed. Tracer particles may be tracked in situ during the simulation, but often they are reconstructed in a post-processing step based on the information saved during the hydrodynamic simulation. Reconstruction can be done in a number of ways, and here we compare the approaches of backward and forward integration of the equations of motion to the results based on inline particle trajectories. We find that both methods agree reasonably well with the inline results for isotopes for which a large number of particles contribute. However, for rarer isotopes that are produced only by a small number of particle trajectories, deviations can be large. For our setup, we find that backward integration leads to better agreement with the inline particles by more accurately reproducing the conditions following freeze-out from nuclear statistical equilibrium, because the establishment of nuclear statistical equilibrium erases the need for detailed trajectories at earlier times. Based on our results, if inline tracers are unavailable, we recommend backward reconstruction to the point when nuclear statistical equilibrium was last applied, with an interval between simulation snapshots of at most 1 ms for nucleosynthesis post-processing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Transfer Factors for the FRMAC Assessment Manual and Turbo FRMAC to Improve Radiological Dose Assessment

The Turbo FRMAC analysis tool is used to perform complex calculations to quickly evaluate radiological consequences and aid in decision making during an emergency response by assessing impacts to the public, workers, and the food supply. Turbo FRMAC calculations are based on methods established by the Federal Radiological Monitoring and Assessment Center (FRMAC). To be able to assess impacts, input data called transfer factors that describe radionuclide uptake by local plants and animals are required. During the code application exercises in late 2016, identifying, finding, and validating nonstandard transfer factors proved to be time consuming and diverted the teams’ activity away from other critical tasks. As a result, a task was undertaken to dramatically expand the list of available transfer factors (food and non-food) and incorporate these factors into Turbo FRMAC. This will ultimately result in improved efficiency of the assessment team to perform calculations during times when the FRMAC is activated and provide more defensible, vetted data from which to calculate results. As a result of these and subsequent exercises, transfer factors were needed for the following items: bell peppers, Christmas tree, deer, flowers, fresh cucumbers, tomatoes, grapefruit, lichens, mushrooms, oranges, snap beans, squash, strawberries, sugarcane for sugar and seed, sweet corn, tea, tobacco, tree bark, and watermelon. To expand the applicability of the information tables, generic transfer data were also provided for common categories and recommendations were made for expanding the list of chemicals based on chemical similarity. The data presented in this report were compiled from recent literature with most of the data encompassing the period from 2000 to 2018. The following radionuclides were targeted during the literature search: elements associated with reactor accidents or nuclear detonations (Sr 89/90, Cs 134/137, Ce 141/144, Ru 103/106, I 129/131/133, Pu 238/239, Am 241, Zr 95, and Nb 95) and elements associated with industrial accidents or dirty bombs (Ir 192 and Co 60). Data from other elements were evaluated if they were identified during the literature search. In addition, reports were evaluated that were recommended by the research consultants. For each plant or animal transfer factor, the goal was to determine the geometric mean, the geometric standard deviation, the minimum, the maximum, and the number of measurements used. If only one measurement was available that was presented as the mean. Concentration ratios with large geometric standard deviations (GSDs) were generally the result of a paucity of measurements or a few disparate measurements. Measurement disparity was observed for data from different soil types. The geometric mean (GM) is a good reference value for planning and responses purposes, but the location-specific concentrations are unlikely to be similar to the model results. Recommendations for further work include: segregating the data to reflect the influence of soil type, developing approaches to incorporate animal data based on aggregated transfer measurements, and including data for foliar deposition on plants.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

TRANSIENT BEAM LOADING EFFECTS IN RF SYSTEMS IN JLEIC

The pulsed electron bunch trains generated from the Continuous Electron Beam Accelerator Facility (CEBAF) linac to inject into the proposed Jefferson Lab Electron Ion Collider (JLEIC) e-ring will produce transient beam loading effects in the Superconducting Radio Frequency (SRF) systems that, if not mitigated, could cause unacceptably large beam energy deviation in the injection capture, or exceed the energy acceptance of CEBAF’s recirculating arcs. In the electron storage ring, the beam abort or ion clearing gaps or uneven bucket filling can cause large beam phase transients in the (S)RF cavity control systems and even beam loss due to Robinson instability. We have first analysed the beam stability criteria in steady state and estimated the transient effect in Feedforward and Feedback RF controls. Initial analytical models for these effects are shown for the design of the JLEIC e-ring from 3GeV to 12GeV.

Wang, Haipeng↗
DOE OSTI2015

Frequency Responsive Demand in U.S. Western Power System Model

Frequency control from the demand side presents a novel and viable way to complement control by traditional generators. This paper presents simulation tests of a previously developed load controller in a large interconnected system model. The controller was implemented as a user defined model in PowerWorld software. WECC system model was used. The simulation tests include comparison of situation with and without control, sensitivities to controller location, control gain parameters, and type of end-use loads. The performance is measured in terms of frequency and intertie power flow deviations after a large generator outage.

Power system dynamics, power system control↗