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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Search for direct production of electroweakinos in final states with one lepton, jets and missing transverse momentum in pp collisions at $\sqrt{s} $ = 13 TeV with the ATLAS detector

Searches for electroweak production of wino-like chargino pairs, ${\overset{\sim }{\chi}}_1^{+}{\overset{\sim }{\chi}}_1^{-}$, and of wino-like chargino and next-to-lightest neutralino, ${\overset{\sim }{\chi}}_1^{\pm }{\overset{\sim }{\chi}}_2^0$, are presented. The models explored assume that the charginos decay into a W boson and the lightest neutralino, ${\overset{\sim }{\chi}}_1^{\pm}\to {W}^{\pm }{\overset{\sim }{\chi}}_1^0$. The next-to-lightest neutralinos are degenerate in mass with the chargino and decay to ${\overset{\sim }{\chi}}_1^0$ and either a Z or a Higgs boson, ${\overset{\sim }{\chi}}_2^0\to Z{\overset{\sim }{\chi}}_1^0$ or $h{\overset{\sim }{\chi}}_1^0$. The searches exploit the presence of a single isolated lepton and missing transverse momentum from the W boson decay products and the lightest neutralinos, and the presence of jets from hadronically decaying Z or W bosons or from the Higgs boson decaying into a pair of b-quarks. The searches use 139 fb –1 of $\sqrt{s}$ = 13 TeV proton-proton collisions data collected by the ATLAS detector at the Large Hadron Collider between 2015 and 2018. No deviations from the Standard Model expectations are found, and 95% confidence level exclusion limits are set. Chargino masses ranging from 260 to 520 GeV are excluded for a massless ${\overset{\sim }{\chi}}_1^0$ in chargino pair production models. Degenerate chargino and next-to-lightest neutralino masses ranging from 260 to 420 GeV are excluded for a massless ${\overset{\sim }{\chi}}_1^0$ for ${\overset{\sim }{\chi}}_2^0\to Z{\overset{\sim }{\chi}}_1^0$. For decays through an on-shell Higgs boson and for mass-splitting between ${\overset{\sim }{\chi}}_1^{\pm }/{\overset{\sim }{\chi}}_2^0$ and ${\overset{\sim }{\chi}}_1^0$ as small as the Higgs boson mass, mass limits are improved by up to 40 GeV in the range of 200–260 GeV and 280–470 GeV compared to previous ATLAS constraints.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Wide-angle effects for peculiar velocities

ABSTRACT The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as ‘wide-angle’ effects in redshift-space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work, we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in autocorrelation and in cross-correlation with the density field. Using both linear theory and the Zel'dovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low-redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration space, Fourier space, and spherical analyses; providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Synthesis of ARM User Facility Surface Rainfall Datasets to Construct a Best Estimate Value Added Product (PrecipBE)

Surface precipitation measurements are essential for Earth system model (ESM) evaluation and understanding cloud processes. An ever-growing need for robust, temporally evolving, and easy-to-use statistical datasets provides motivation for a baseline ground-based precipitation properties data product. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility operates an extensive suite of precipitation instruments with various sensitivities and operating mechanisms, which render the decision of which instrument to use based on one or more fixed thresholds challenging and prone to errors and bias. Using a long-term instrument inter-comparison from a unique per-precipitation event perspective, rather than instantaneous sample comparison, we demonstrate that ARM rainfall-measuring instruments are generally consistent with each other at the statistical level. Inter-instrument deviations at the single event level can be large, especially for specific rainfall event properties such as maximum precipitation rates. A machine-learning (ML) analysis using a random forest regressor indicates that in some cases, depending on instrument, local site climatology, and/or specific deployment configuration, certain atmospheric state variables influence the measured quantities in an unpredictable manner. Thus, a-priori weighting of different instruments does not necessarily lead to more accurate and less biased synthesis of instrument data. These results motivate the design of the ARM precipitation best-estimate (PrecipBE) value-added product, which incorporates all valid precipitation data while considering data quality and other instrument limitations. PrecipBE consists of time series and tabular statistics datasets in an easy-to-use and insightful per-precipitation event format. It provides a large set of precipitation event properties supplemented with ancillary data from ARM datasets that correspond to the detected precipitation events. We describe the PrecipBE algorithm and demonstrate its use via the examination of a single-day output as well as a long-term trend analysis of precipitation events at the ARM Southern Great Plains (SGP) site, covering more than 30 years of data. The trend analysis tentatively suggests a long-term temporal tendency for mainly shorter and less intense precipitation events at the SGP site, but a long-term increase in annual rainfall by more than 36 mm (5 %) per decade. This rainfall trend is catalyzed primarily by more extreme event properties of relatively rare, intense precipitation events, with event total and 1 min maximum precipitation rate at a 1 year timeframe increasing up to 5 mm and 9 mm h −1 (several percent) per decade, respectively. While the currently available PrecipBE datasets (at https://adc.arm.gov/discovery/, last access: 8 December 2025) cover rainfall from multiple ARM deployments up to March 2025, PrecipBE is planned to be expanded to include solid-phase precipitation and will soon become an operational product with a several-day lag from real-time. We invite the ARM user community to leverage this new product and welcome user feedback to further enhance the dataset.

Silber, Israel [Pacific Northwest National Laborat↗
DOE OSTI2022

\(J/\psi \) and \(\psi (2S)\) Production in Small Systems with PHENIX

The suppression of the ψ(2S) nuclear modification factor has been seen as a trademark signature of final-state effects in large collision systems for decades. In small systems, deviations of the nuclear modification from unity had been attributed to cold nuclear matter effects until the observation of strong differential suppression of the ψ(2S) state in p/d+A collisions, which suggests the presence of final-state effects. Here in this paper, we present results of J/ψ and ψ(2S) measurements in the dimuon decay channel for p + p, p+Al, and p+Au collision systems at $\sqrt{s_{NN}}$ = 200 GeV. Key results include the nuclear modification factors R pA as a function of centrality and rapidity. The measurements are compared with shadowing and transport model predictions, as well as complementary measurements at LHC energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1974

Vibrationally excited nitrogen in stable auroral red arcs and its effect on ionospheric recombination

The time-dependent continuity equations, including diffusion, were solved for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to stable auroral red (SAR) arcs. The results show that molecular nitrogen is excited vibrationally to the degree that the rate constant for the ionospheric loss process, O(+) + N2 yields NO(+) + N, is increased by as much as a factor of 7.6 at F2 region altitudes. It was found that deviations from the energetically equivalent Boltzmann distribution were large, causing the rate constant to be as much as 1.6 times the rate constant calculated for the Boltzmann distribution. These results indicate that SAR arc intensities as small as 58 R can produce noticeable increases in the ionosphere ion-atom interchange reaction rate and hence in the rate of loss of ionospheric electrons. It is suggested that the observed decrease of electron density in the F2 region in SAR arcs can probably be explained by enhanced reaction rates for ion-atom interchange between O(+) and N2 caused by vibrational excitation of molecular nitrogen by electron impact.

Newton, G. P.↗
NASA NTRS1979

A Comparison of the Aerodynamic Characteristics of Eight Sailwing Airfoil Sections

The operational characteristics of sailwings are discussed with emphasis placed on the importance of the trailing edge cable tension. The three dimensional aerodynamic characteristics were obtained from wind tunnel tests, and the results compared to determine the magnitude of the aerodynamic penalties paid for various structural simplifications. For the sectional thickness ratios, it is concluded that, while the basic double-membraned sailwing has exceptional aerodynamic performance, even superior for some applications to the conventional hardwing, any notable deviation from this configuration results in an unacceptably large performance penalty.

Maughmer, M. D.↗
NASA NTRS1985

Orientation and temperature dependence of some mechanical properties of the single-crystal nickel-base superalloy Rene N4. 3: Tension-compression anisotropy

Single crystal superalloy specimens with various crystallographic directions along their axes were tested in compression at room temperature, 650, 760, 870, and 980 deg C. These results are compared with the tensile behavior studied previously. The alloy, Rene N4, was developed for gas turbine engine blades and has the nominal composition 3.7 Al, 4.2 Ti, 4 Ta, 0.5 Nb, 6 W, 1.5 Mo 9 Cr. 7.5 Co, balance Ni, in weight percent. Slip trace analysis showed that primary cube slip occurred even at room temperature for the 111 specimens. With increasing test temperature more orientations exhibited primary cube slip, until at 870 deg C only the 100 and 011 specimens exhibited normal octahedral slip. The yield strength for octahedral slip was numerically analysed using a model proposed by Lall, Chin, and Pope to explain deviations from Schmid's Law in the yielding behavior of a single phase Gamma prime alloy, Ni3(Al, Nb). The Schmid's Law deviations in Rene N4 were found to be largely due to a tension-compression anisotropy. A second effect, which increases trength for orientations away from 001, was found to be small in Rene N4. Analysis of recently published data on the single crystal superalloy PWA 1480 yielded the same result.

Miner, R. V.↗
NASA NTRS1985

Trends in the Southern Hemisphere

The station data of monthly mean sea level pressure, surface air temperature as well as monthly mean upper air rawinsonde observations were used to study the trends in the southern winter (June, July and August) and summer (December, January and February). This study covers the period 1951 to 1981. Most of the surface data have continuous 31 year records, however, for some high latitude stations, both surface and upper air records started in 1957 or even later. The normalized anomaly for a given variable Z is defined as the departure from the 1951 to 1981 grand mean divided by the standard deviation during the same period. To eliminate the large fluctuations in the data, the 1-3-1 smoothing was applied to the time series of normalized anomalies. The trends of monthly mean temperatures and height from the surface to 200 mb level is documented.

Mo, K. C.↗
NASA NTRS1986

Compensation for Primary Reflector Wavefront Error

The object of the invention is to compensate for errors in a large telescope primary reflector by making certain compensating deviations in a smaller, auxiliary reflector of the telescope. At least one intermediate element forms an image of the primary surface onto the secondary surface, so each point on the secondary surface corresponds to a point on the primary surface. The secondary surface is formed with a deviation from an ideal secondary surface, with the piston distance of each point on the actual secondary surface equal to the piston distance of a corresponding piston on the actual primary surface from the ideal primary surface. It is found that this results in electromagnetic (e.g., light) rays which strike a deviating area of the actual primary surface being brought to the same focus as if the actual primary surface did not have a diviation from an ideal primary surface.

Meinel, A. B.↗
NASA NTRS1986

Lithologic mapping of silicate rocks using TIMS

Common rock-forming minerals have thermal infrared spectral features that are measured in the laboratory to infer composition. An airborne Daedalus scanner (TIMS) that collects six channels of thermal infrared radiance data (8 to 12 microns), may be used to measure these same features for rock identification. Previously, false-color composite pictures made from channels 1, 3, and 5 and emittance spectra for small areas on these images were used to make lithologic maps. Central wavelength, standard deviation, and amplitude of normal curves regressed on the emittance spectra are related to compositional information for crystalline igneous silicate rocks. As expected, the central wavelength varies systematically with silica content and with modal quartz content. Standard deviation is less sensitive to compositional changes, but large values may result from mixed admixture of vegetation. Compression of the six TIMS channels to three image channels made from the regressed parameters may be effective in improving geologic mapping from TIMS data, and these synthetic images may form a basis for the remote assessment of rock composition.

Gillespie, A. R.↗
NASA NTRS1999

The NASA/MSFC Global Reference Atmospheric Model: 1999 Version (GRAM-99)

The latest version of Global Reference Atmospheric Model (GRAM-99) is presented and discussed. GRAM-99 uses either (binary) Global Upper Air Climatic Atlas (GUACA) or (ASCII) Global Gridded Upper Air Statistics (GGUAS) CD-ROM data sets, for 0-27 km altitudes. As with earlier versions, GRAM-99 provides complete geographical and altitude coverage for each month of the year. GRAM-99 uses a specially-developed data set, based on Middle Atmosphere Program (MAP) data, for 20-120 km altitudes, and NASA's 1999 version Marshall Engineering Thermosphere (MET-99) model for heights above 90 km. Fairing techniques assure smooth transition in overlap height ranges (20-27 km and 90-120 km). GRAM-99 includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He and H). A variable-scale perturbation model provides both large-scale (wave) and small-scale (stochastic) deviations from mean values for thermodynamic variables and horizontal and vertical wind components. The small-scale perturbation model includes improvements in representing intermittency ("patchiness"). A major new feature is an option to substitute Range Reference Atmosphere (RRA) data for conventional GRAM climatology when a trajectory passes sufficiently near any RRA site. A complete user's guide for running the program, plus sample input and output, is provided. An example is provided for how to incorporate GRAM-99 as subroutines in other programs (e.g., trajectory codes).

Justus, C. G.↗
NASA NTRS2006

Visualization of a Capsule Entry Vehicle Reaction-Control System (RCS) Thruster

Planar laser-induced fluorescence (PLIF) was used to visualize the reaction control system (RCS) jet flow emanating from the aft-body of an Apollo-geometry capsule test article in the NASA Langley Research Center 31-Inch Mach 10 Air wind tunnel. The RCS jet was oriented normal to the aft surface of the model and had a nominal Mach number of 2.94. The composition of the jet gas by mass was 95% nitrogen (N2) and 5% nitric oxide (NO). The RCS jet flowrate varied between zero and 0.5 standard liters per minute and the angle of attack and tunnel stagnation pressure were also varied. PLIF was used to excite the NO molecules for flow visualization. These flow visualization images were processed to determine the trajectory and to quantify the flapping of the RCS jet. The spatial resolution of the jet trajectory measurement was about 1 mm and the single-shot precision of the measurement was estimated to be 0.02 mm in the far field of the jet plume. The jet flapping, measured by the standard deviation of the jet centerline position was as large as 0.9 mm, while the jet was 1.5-4 mm in diameter (full width at half maximum). Schlieren flow visualization images were obtained for comparison with the PLIF. Surface pressures were also measured and presented. Virtual Diagnostics Interface (VIDI) technology developed at NASA Langley was used to superimpose and visualize the data sets. The measurements demonstrate some of the capabilities of the PLIF method while providing a test case for computational fluid dynamics (CFD) validation.

Danehy, P. M.↗
NASA NTRS2016

Soil Moisture Active Passive Mission L4_SM Data Product Assessment (Version 2 Validated Release)

During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.

SMAP L4_SM↗
NASA NTRS2022

Determining the Primary Sources of Uncertainty in Retrieval of Marine Remote Sensing Reflectance From Satellite Ocean Color Sensors

Uncertainties in the retrieval of the remote sensing reflectance, R_rs, from Ocean Color (OC) satellite sensors have a strong impact on the performance of algorithms for the estimation of chlorophyll-a, mineral concentrations, and inherent optical properties (IOPs). The uncertainties are highest in the blue bands. The total radiance measured at the top of the atmosphere captures the instantaneous state of the atmosphere-ocean system: the in-water conditions, sky and Sun glint reflected from the wind-roughened ocean surface, as well as light scattered from molecules and aerosols in the atmosphere. Each of these components has associated uncertainties, and when combined with the additional uncertainties from the instrument noise and the atmospheric correction process, they contribute to the total uncertainty budget for the retrieved R_rs. We analyzed the contribution of each component uncertainties to the total R_rs uncertainties in SNPP-VIIRS level 2 products, taking advantage of the spectral differences between the components. We examined multiple scenes in the open ocean and coastal waters at spatial resolutions ranging from 2250 to 5250 m by comparing the retrieved R_rs to in situ measurements made at several AERONET-OC sites and at the MOBY site. It was shown that uncertainties associated with the molecular (Rayleigh) scattering play the most significant role, while the contributions of other components are usually smaller. Uncertainties in Rayleigh scattering are primarily attributed to the variability of Rayleigh optical thickness (ROT) with a standard deviation of approximately 1.5% of ROT, which can largely explain the frequency of negative R_rs retrievals as observed using the current standard atmospheric correction process employed by NASA. Variability of the sky light reflected from the ocean surface in some conditions also contributed to uncertainties in the blue; water variability proportional to R_rs had a very pronounced peak in the green at coastal sites.

remote sensing reflectance↗
NASA NTRS2016

Conjunction Assessment Late-Notice High-Interest Event Investigation: Space Weather Aspects

Late-notice events usually driven by large changes in primary (protected) object or secondary object state. Main parameter to represent size of state change is component position difference divided by associated standard deviation (epsilon divided by sigma) from covariance. Investigation determined actual frequency of large state changes, in both individual and combined states. Compared them to theoretically expected frequencies. Found that large changes ( (epsilon divided by sigma) is greater than 3) in individual object states occur much more frequently than theory dictates. Effect is less pronounced in radial components and in events with probability of collision (Pc) greater than 1 (sup -5) (1e-5). Found combined state matched much closer to theoretical expectation, especially for radial and cross-track. In-track is expected to be the most vulnerable to modeling errors, so not surprising that non-compliance largest in this component.

Conjunction↗
DOE OSTI2023

Rotor and wake aerodynamic analysis of the Hybrid-Lambda concept - an offshore low-specific-rating rotor concept

The low-specific-rating rotor concept Hybrid-Lambda introduces a blade design with nonuniform distribution of design variables (design tip speed ratio and axial induction) along the blade span to alleviate loads of the outboard section in strong winds. In this paper, we validate aerodynamic design calculations, which were carried out with the blade element momentum theory by comparing the results to free vortex wake investigations (FVW). Furthermore, we investigate the wake behaviour with FVW and large-eddy simulations. The results show good agreements between the blade element momentum theory and FVW for integrated rotor quantities (power and thrust). Small deviations are present when the gradients of axial induction along the span are large. The wake of the Hybrid-Lambda Rotor shows advantages in the near-wake region (up to 4 diameters [D] downstream), especially in the outer wake annulus and in low turbulence scenarios. For further downstream positions, the wake is comparable to that of the reference turbine.

17 WIND ENERGY↗
DOE OSTI2023

Frequency Recovery in Power Grids Using High-Performance Computing

Maintaining electric power system stability is paramount, especially in extreme contingencies involving unexpected outages of multiple generators or transmission lines that are typical during severe weather events. Such outages often lead to large supply-demand mismatches followed by subsequent system frequency deviations from their nominal value. The extent of frequency deviations is an important metric of system resilience, and its timely mitigation is a central goal of power system operation and control. This paper develops a novel nonlinear model predictive control (NMPC) method to minimize frequency deviations when the grid is affected by an unforeseen loss of multiple components. Our method is based on a novel multi-period alternating current optimal power flow (ACOPF) formulation that accurately models both nonlinear electric power flow physics and the primary and secondary frequency response of generator control mechanisms. We develop a distributed parallel Julia package for solving the large-scale nonlinear optimization problems that result from our NMPC method and thereby address realistic test instances on existing high-performance computing architectures. Our method demonstrates superior performance in terms of frequency recovery over existing industry practices, where generator levels are set based on the solution of single-period classical ACOPF models.

nonlinear model predictive control↗