Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble prediction”

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 469 records · Page 26

Using Data Assimilation Methods for Prediction of Solar Activity

The variability of solar magnetic activity known as the 11-­‐year solar cycles has the longest history of observations. These solar cycles dramatically affect conditions in the heliosphere and the Earth’s space environment. Our current understanding of the physical processes that make up global solar dynamics and the dynamo that generates the magnetic fields is sketchy, resulting in unrealistic descriptions in theoretical and numerical models of the solar cycles. The absence of long-­‐term observations of solar interior dynamics and photospheric magnetic fields hinders development of accurate dynamo models and their calibration. In such situations, mathematical data assimilation methods provide an optimal approach for combining the available observational data and their uncertainties with theoretical models in order to estimate the state of the solar dynamo and predict future cycles. In this presentation, we will discuss the implementation and performance of an Ensemble Kalman Filter data assimilation method based on the Parker migratory dynamo model complemented by the equation of magnetic helicity conservation and long-­‐term sunspot data series. This approach has allowed us to reproduce the general properties of the solar cycles and has already demonstrated a good predictive capability for the current cycle, 24. We will discuss further development of this approach, which includes a more sophisticated dynamo model, synoptic magnetogram data, and employs the DART Data Assimilation Research Testbed.

Solar Activity↗

Ab initio prediction of an order-disorder transition in Mg 2 Ge O 4 : Implication for the nature of super-Earth's mantles

Here we present an ab initio prediction of an order-disorder transition (ODT) from a I¯42d-type to Th 3 P 4 -type phase in the cation sublattices of Mg 2 GeO 4 , a post-post-perovskite phase. This uncommon type of prediction is achieved by carrying out a high-throughput sampling of atomic configurations in a 56-atom supercell followed by a Boltzmann ensemble statistics calculation. Mg 2 GeO 4 is a low-pressure analog of I¯42d-type Mg 2 SiO 4 , a predicted major planet-forming phase of super-Earths' mantles. Therefore, a similar ODT is anticipated in I¯42d-type Mg 2 SiO 4 as well, which should impact the internal structure and dynamics of these planets. Furthermore, the prediction of this Th 3 P 4 -type phase in Mg 2 GeO 4 further enhances the relationship between the crystal structures of Earth/planet-forming silicates and oxides at extreme pressures and those of rare-earth sesquisulfides at low pressures.

36 MATERIALS SCIENCE↗

Operational Dust Prediction

Over the last few years, numerical prediction of dust aerosol concentration has become prominent at several research and operational weather centres due to growing interest from diverse stakeholders, such as solar energy plant managers, health professionals, aviation and military authorities and policymakers. Dust prediction in numerical weather prediction-type models faces a number of challenges owing to the complexity of the system. At the centre of the problem is the vast range of scales required to fully account for all of the physical processes related to dust. Another limiting factor is the paucity of suitable dust observations available for model, evaluation and assimilation. This chapter discusses in detail numerical prediction of dust with examples from systems that are currently providing dust forecasts in near real-time or are part of international efforts to establish daily provision of dust forecasts based on multi-model ensembles. The various models are introduced and described along with an overview on the importance of dust prediction activities and a historical perspective. Assimilation and evaluation aspects in dust prediction are also discussed.

assimilation↗

Probabilistic Day-Ahead Forecasting Using an Analog Ensemble Approach for Wind Farm Grid Services

Wind resource assessment and wind power forecasting are used in research and industry to anticipate future power output at scales ranging from individual wind turbines to entire wind farms. Probabilistic day-ahead wind forecasting is useful for anticipating how a wind farm could potentially participate in the day-ahead market by providing upper and lower bounds for expected power generation, thus informing grid operators of its uncertainty. Understanding this uncertainty is part of a larger project focused on building a platform that combines efforts in weather forecasting, aerodynamic and economic modeling to create maximum value of a wind plant to better provide services to the grid. This effort is also known as the Atmosphere to Electrons to Grid (A2E2G) project. One method for producing a probabilistic forecast is through the analog ensemble approach (Delle Monache et al., 2011). This method leverages historical forecasts and their corresponding observations as a training data set from which future forecasts can be made. For some future forecast, the most similar historical forecasts (analogs) are identified on a regular time basis such as once per a 3-hour window. The most similar analogs, based on a metric such as root mean square error (RMSE), are recorded and their corresponding verifying observations are used as an ensemble member for this future forecast. Prior work in this area demonstrates improvements over raw Numerical Weather Prediction (NWP) forecasts and shows skill similar to techniques such as logistic regression and machine learning (Delle Monache et al., 2013; Alessandrini et al., 2015). Here, we take the High-Resolution Rapid Refresh model (HRRR) day-ahead forecast (0-36 hours) to create a probabilistic day-ahead forecast using an analog ensemble approach. The HRRR has an hourly temporal resolution, with a spatial resolution of 3 km. The 12 UTC HRRR model run is downloaded every day for one year from August 2019 - July 2020, with the first 11 months serving as a bank of analogs from which the forecasting algorithm can create a probabilistic forecast. Once downloaded, the original HRRR forecast is temporally interpolated to 5-minutes, aligning with both the temporal resolution of the observations as well as the timescale relevant for day-ahead power forecasts. The forecast is validated at the M2 tower at the Flatirons Campus of the National Renewable Energy Laboratory (NREL) at a typical wind turbine height of 80 m. Variables such as wind speed, wind direction, and turbulence intensity are incorporated into the probabilistic forecast model and weighted according to their relative importance to the forecast. Based on metrics such as mean bias error (MBE), mean absolute error (MAE), and root mean square error, the analog ensemble forecast outperforms the raw HRRR forecast during the testing period of July 2020. Figure 1 illustrates an example day-ahead forecast compared against the verifying observations. The general variability and ramps are captured throughout the day, with potential to further improve the analog ensemble model through machine learning techniques.

numerical weather prediction↗

Attribution of NAO Predictive Skill Beyond 2 Weeks in Boreal Winter

Weeks 3–6 averaged winter North Atlantic Oscillation (NAO) predictive skill in a state-of-the-art coupled climate prediction system is attributed to two principle sources: upper and lower boundary conditions linked to the stratosphere and El Niño-Southern Oscillation (ENSO), respectively. A 20-member ensemble of 45-day reforecasts over 1999–2015 is utilized, together with uninitialized simulations with the atmospheric component of the prediction system forced with observed radiative forcing and lower boundary conditions. NAO forecast skill for lead times out to 6 weeks is higher following extreme stratospheric polar vortex conditions (weak and strong vortex events) compared to neutral states. Enhanced weeks 3-6 NAO predictive skill for weak vortex events results primarily from stratospheric downward coupling to the troposphere, while enhanced skill for strong vortex events can be partly attributed to lower boundary forcing related to the ENSO phenomenon. Implications for forecast system development and improvement are discussed.

54 ENVIRONMENTAL SCIENCES↗

Causes of Long-Term Drought in the United States Great Plains

The United States Great Plains (USGP) experienced a number of multi-year droughts during the last century, most notably the droughts of the 1930s and 1950s. This study examines the causes of such droughts using ensembles of long term (1930-1999) simulations carried out with the NASA Seasonal-to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed sea surface temperatures (SSTs). The results show that the model produces long-term (multi-year) variations in the USGP precipitation that are similar to those observed. A correlative analysis suggests that the ensemble mean low frequency (time scales longer than about 6 years) rainfall variations in the USGP are linked to a pan-Pacific pattern of SST variability that is the leading empirical orthogonal function (EOF) in the low frequency SST data. The link between the SST and the Great Plains precipitation is confirmed in idealized AGCM simulations, in which the model is forced by the 2 polarities of the pan-Pacific SST pattern. The idealized simulations further show that it is primarily the tropical part of the SST anomalies that influence the USGP. As such, the USGP tend to have above normal precipitation when the tropical Pacific SSTs are above normal, while there is a tendency for drought when the tropical SSTs are cold. The upper tropospheric response to the pan-Pacific SST EOF shows a global-scale pattern with a strong wave response in the Pacific and a substantial zonally-symmetric component in which USGP pluvial (drought) conditions are associated with reduced (enhanced) heights throughout the extra-tropics. The potential predictability of rainfall in the USGP associated with SSTs is rather modest, with on average about 1/3 of the total low frequency rainfall variance forced by SST anomalies. Further idealized experiments with climatological SST, suggest that the remaining low frequency variance in the USGP precipitation is the result of interactions with soil moisture. In particular, simulations with soil moisture feedback show a six-fold increase in the variance in annual USGP precipitation compared with simulations in which the soil feedback is excluded. In addition to increasing variance, the interactions with the soil introduce year-to-year memory in the hydrological cycle that is consistent with a red noise process, in which the low frequencies in the deep soil are the result of integrating a net forcing (precipitation-evaporation-runoff) that is white noise on interannual time scales. As such, the role of low frequency SST variability is to introduce a bias to the net forcing on the soil moisture that drives the random process preferentially to either wet or dry conditions.

Schubert, Siegfried D.↗

Predicting Emission Source Terms in a Reduced-Order Fire Spread Model—Part 1: Particulate Emissions

A simple, easy-to-evaluate, surrogate model was developed for predicting the particle emission source term in wildfire simulations. In creating this model, we conceptualized wildfire as a series of flamelets, and using this concept of flamelets, we developed a one-dimensional model to represent the structure of these flamelets which then could be used to simulate the evolution of a single flamelet. A previously developed soot model was executed within this flamelet simulation which could produce a particle size distribution. Executing this flamelet simulation 1200 times with varying conditions created a data set of emitted particle size distributions to which simple rational equations could be tuned to predict a particle emission factor, mean particle size, and standard deviation of particle sizes. These surrogate models (the rational equation) were implemented into a reduced-order fire spread model, QUIC-Fire. Using QUIC-Fire, an ensemble of simulations were executed for grassland fires, southeast U.S. conifer forests, and western mountain conifer forests. Resulting emission factors from this ensemble were compared against field data for these fire classes with promising results. Also shown is a predicted averaged resulting particle size distribution with the bulk of particles produced to be on the order of 1 μm in size.

54 ENVIRONMENTAL SCIENCES↗

Mass Conservation and Positivity Preservation with Ensemble-type Kalman Filter Algorithms

Maintaining conservative physical laws numerically has long been recognized as being important in the development of numerical weather prediction (NWP) models. In the broader context of data assimilation, concerted efforts to maintain conservation laws numerically and to understand the significance of doing so have begun only recently. In order to enforce physically based conservation laws of total mass and positivity in the ensemble Kalman filter, we incorporate constraints to ensure that the filter ensemble members and the ensemble mean conserve mass and remain nonnegative through measurement updates. We show that the analysis steps of ensemble transform Kalman filter (ETKF) algorithm and ensemble Kalman filter algorithm (EnKF) can conserve the mass integral, but do not preserve positivity. Further, if localization is applied or if negative values are simply set to zero, then the total mass is not conserved either. In order to ensure mass conservation, a projection matrix that corrects for localization effects is constructed. In order to maintain both mass conservation and positivity preservation through the analysis step, we construct a data assimilation algorithms based on quadratic programming and ensemble Kalman filtering. Mass and positivity are both preserved by formulating the filter update as a set of quadratic programming problems that incorporate constraints. Some simple numerical experiments indicate that this approach can have a significant positive impact on the posterior ensemble distribution, giving results that are more physically plausible both for individual ensemble members and for the ensemble mean. The results show clear improvements in both analyses and forecasts, particularly in the presence of localized features. Behavior of the algorithm is also tested in presence of model error.

EnKR↗

Numerical and experimental investigation of nonswirling and swirling confined jets

An investigation of the influence of large scale structures on the flow development for coaxial jets with sudden expansion (with and without swirl) is presented. Both an experimental study and numerical predictions were performed for a configuration corresponding to that considered by Johnson and Bennett and Roback and Johnson. The effects of large scale structures on the swirling and nonswirling flows are documented, in particular their influence on turbulence modeling and the numerical simulation. The ensemble-averaged, time-dependent Navier-Stokes equations are solved by an LBI procedure to predict the turbulent flow field. Effects of artificial dissipation and placement of the upstream boundary in the numerical computation are also discussed.

Brondum, D. C.↗

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab initio studies on structural and thermodynamic properties of magnetic Fe

The present work systematically investigates the total energy, phonon spectra, and thermodynamic properties of different polymorphs of pure Fe, i.e., FCC, BCC, and HCP, with the ab initio approach, considering various magnetic configurations. In general, the calculated energy vs. volume curves and phonon spectra agree well with previous calculations and the experimental data. In addition, their thermodynamic properties are estimated by the quasiharmonic approximation (QHA). Specifically, a superposition approach based on the latest Zentropy theory was utilized to predict magnetic transition temperatures and thermodynamic properties of pure Fe. With the ensemble of the partition function considering the multiplicity of each magnetic microstate, the current work successfully reproduced the Curie/Néel temperature and the Schottky anomaly of heat capacity in FCC, BCC, and HCP Fe purely based on the ab initio input, which exhibits good agreement with the experimental data and CALPHAD modeling.

36 MATERIALS SCIENCE↗

Electrocatalytic Hydrogen Evolution at Full Atomic Utilization over ITO-Supported Sub-nano-Pt n Clusters: High, Size-Dependent Activity Controlled by Fluxional Pt Hydride Species

A combination of density functional theory (DFT) and experiments with atomically size-selected Pt n clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Pt n size on the electrocatalytic activity of Pt n (n = 1, 4, 7, 8) for the hydrogen evolution reaction (HER). Activity is found to be negligible for isolated Pt atoms on ITO, increasing rapidly with Pt n size, such that Pt 7 /ITO and Pt 8 /ITO have roughly double the activity per Pt atom compared to atoms in the surface layer of polycrystalline Pt. Both DFT and experiment find that hydrogen under-potential deposition (H upd ) results in Pt n /ITO (n = 4, 7, 8) adsorbing ~2 H atoms/Pt atom at the HER threshold potential, equal to ca. double the Hupd observed for Pt bulk or nanoparticles. Here, the cluster catalysts under electrocatalytic conditions are hence best described as a Pt hydride compound, significantly departing from a metallic Pt cluster. The exception is Pt 1 /ITO, where H adsorption at the HER threshold potential is energetically unfavorable. Theory combines global optimization with grand canonical approaches for the influence of potential, uncovering that several metastable structures contribute to HER, changing with the applied potential. It is hence critical to include reactions of the ensemble of energetically accessible Pt n H x /ITO structures to correctly predict the activity vs. Pt n size and applied potential. For the small clusters, spillover of H ads from the clusters to the ITO support is significant, resulting in a competing channel for loss of H ads , particularly at slow potential scan rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Amplification of Subtropical Extreme Precipitation in a Warming Climate

Projected precipitation changes in a warming climate vary considerably, spatially, and between intensities. The changes can be greater or less than the ~7% K -1 Clausius-Clapeyron (CC) prediction, owing to dynamic effects. Using two global-climate-model large ensembles, we quantify the dynamically induced changes to precipitation extremes from the present (1996–2005) to late-21st-century (2071–2080) climates, as a function of recurrence interval, focusing on the subtropics. We separate non-CC changes into a term proportional to the present-day vertical-velocity spatial pattern (i.e., an amplification or damping thereof by a constant factor) and a residual. The amplitude term varies with recurrence interval, approximately canceling (doubling) CC for moderate (large) extremes, increasing precipitation variability. Contrastingly, the residual is quasi-uniform across recurrence intervals but spatially heterogeneous, weakening extremes over dry zones. This residual may be related to Hadley cell expansion, although this explanation is insufficient to explain many features, and other possible mechanisms are discussed.

54 ENVIRONMENTAL SCIENCES↗

Two dimensional turbulence in inviscid fluids or guiding center plasmas

Analytic theory for two-dimensional turbulent equilibria for the inviscid Navier-Stokes equations is examined mathematically. Application of the technique to electrostatic guiding center plasma is discussed. A good fit is demonstrated for the approach to a predicted energy per Fourier mode obtained from a two-temperature canonical ensemble. Negative as well as positive temperature regimes are explored. Fluctuations about the mean energy per mode also compare well with theory. In the regime of alpha less than zero, beta greater than zero, with the minimum value of alpha plus beta times k squared near zero, contour plots of the stream function reveal macroscopic vortex structures similar to those seen previously in discrete vortex simulations. Eulerian direct interaction equations, which can be used to follow the approach to inviscid equilibrium, are derived.

Seyler, C. E., Jr.↗

MHD solution of interplanetary disturbances generated by simulated velocity perturbations

An MHD time-dependent numerical simulation, restricted to the solar equatorial plane, is used to demonstrate the interplanetary disturbances caused by several simplified coronal holes. Each 'hole' is assumed to have a configuration such that the higher solar wind velocity produced within their longitudinal extent is Gaussian over a 7-day period at the inner boundary (0.3 AU) of the numerical simulation. A second, twin coronal hole is assumed to rotate on the solar disk behind its predecessor. It is shown that the first coronal hole-produced interplanetary shock ensemble is overtaken by the second ensemble because of the higher velocity, lower density environment into which the latter propagates. A number of features predicted by MHD similarity theory are confirmed by the numerical simulation. These features include (1) strong azimuthal magnetic and plasma density compression, accompanied with average temperature depression, at the contact surface between forward and reverse shock ensembles, and (2) increasing spatial separation distance between forward and reverse shocks.

Dryer, M.↗

Stationarity of magnetohydrodynamic fluctuations in the solar wind

Solar wind research and studies of charged particle propagation often assume that the interplanetary magnetic field represents a stationary random process. The extent to which ensemble averages of the solar wind magnetic fields follow the asymptotic behavior predicted by the ergodic theorem was investigated. Several time periods, including a span of nearly two years, are analyzed. Data intervals which span many solar rotations satisfy the conditions of weak stationarity if the effects of solar rotation are included in the asymptotic analysis. Shorter intervals which include a small integral number of interplanetary sectors also satisfy weak stationarity. The results are illustrated using magnetometer data from the ISEE-3, Voyager and IMP spacecraft.

Matthaeus, W. H.↗

On the Causes of the 1930s Dust Bowl

During the 1930s the United States experienced one of the most devastating droughts of the last century. The drought affected almost 2/3 of the country and parts of Mexico and Canada and was infamous for the numerous dust storms that lead to the characterization of much of the Great Plains as the "Dust Bowl". Results from an ensemble of 100-year simulations with the NASA Seasonal- to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed SSTs show that the model reproduces the basic features of the 1930s drought. In t h s study we exploit this realism in the model simulation to examine in more detail the role of the SST and soil moisture in the development and maintenance the 1930s drought. Results will be presented from simulations in which the SST anomalies are confined to the separate ocean basins, as well as from runs in which soil moisture feedback is turn off.

Schubert, Siegfried↗