Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Development and Comparison of Purgeable Mercury Values in SRR Samples Measured by SRNL and Eurofins FGS

Savannah River Remediation (SRR) requested the development of mercury speciation capabilities at the Savannah River National Laboratory (SRNL) to support the Liquid Waste Operations at SRS. As part of that method development, SRR requested that SRNL Analytical Development (AD) compare their results with those obtained from their outside contract laboratory, Eurofins Frontier Global Sciences (FGS). This document reports on this method development work performed at SRNL as well as the comparative analyses conducted between the two laboratories. Development, optimization, and validation were undertaken at SRNL to produce a method for the species-specific analysis of dissolved gaseous mercury (DGM) in samples from SRR. It was determined that DGM present in SRR samples represented a combination of various soluble and volatile mercury species- and is thus referred to in this document as purgeable mercury. A method has been developed and analytically validated for the analysis of purgeable mercury in caustic nuclear waste. Using multivariate optimization, an efficient and sensitive method was produced. Total recovery of spiked purgeable mercury from deionized water was 99.4 ±10.6% and 91.6 ±14.9% from spiked high-activity tank waste samples. Calibration, by external calibration and standard addition calibration, maintained a mean accuracy of 101 ±5% and 99.9 ±1.7%, respectively. An interlaboratory comparison was performed between SRNL and Eurofins FGS in the co-analysis of three quarterly SRR Tank 50 samples. The statistical evaluations presented in this document support the following conclusions: for internal validation data, the two laboratories demonstrate biases of -4.3% for SRNL and 4.4% for Eurofins FGS, though neither are statistically significant at the 5% level. While the percent relative standard deviation for each of the laboratories is rather large (21.5% for SRNL and 18.4% for Eurofins FGS), there is no indication of a statistically significant difference between the two laboratories in the precision of these results. For the three sets of SRR quarterly Tank 50 samples, there is no indication of a statistically significant bias between the two laboratories at the 5% significance level. These results indicate that SRNL may provide data to the customer of statistically similar quality to Eurofins FGS.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗
DOE OSTI2019

Phenomenology of large scale structure in scalar-tensor theories: Joint prior covariance of w DE , Σ, and μ in Horndeski theories

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, w DE , and the phenomenology of large scale structure. They will allow us to constrain deviations from the Λ cold dark matter predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions Σ and μ, respectively. In this work, we derive the theoretical prior for the joint covariance of w DE , Σ and μ, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. Here, we identify several interesting features and trends in the distribution functions of w DE , Σ and μ, as well as in their covariances; we confirm the high degree of correlation between Σ and μ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly w DE , Σ and μ in a nonparametric way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

How Far Is Quasar UV/Optical Variability from a Damped Random Walk at Low Frequency?

Studies have shown that UV/optical light curves of quasars can be described using the prevalent damped random walk (DRW) model, also known as the Ornstein–Uhlenbeck process. A white noise power spectral density (PSD) is expected at low frequency in this model; however, a direct observational constraint to the low-frequency PSD slope is difficult due to the limited lengths of the light curves available. Meanwhile, quasars show scatter in their DRW parameters that is too large to be attributed to uncertainties in the measurements and dependence on the variation of known physical factors. In this work we present simulations showing that, if the low-frequency PSD deviates from the DRW, the red noise leakage can naturally produce large scatter in the variation parameters measured from simulated light curves. The steeper the low-frequency PSD slope, the larger scatter we expect. Based on observations of SDSS Stripe 82 quasars, we find that the low-frequency PSD slope should be no steeper than −1.3. The actual slope could be flatter, which consequently requires that the quasar variabilities should be influenced by other unknown factors. We speculate that the magnetic field and/or metallicity could be such additional factors.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗