Search NASA⌕ Search

SEARCH · Search NASA

Results for “Prediction Chaos”

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 55 records · Page 3

Efficient Probabilistic Prediction and Uncertainty Quantification of Tropical Cyclone–Driven Storm Tides and Inundation

Abstract This study proposes and assesses a methodology to obtain high-quality probabilistic predictions and uncertainty information of near-landfall tropical cyclone–driven (TC-driven) storm tide and inundation with limited time and resources. Forecasts of TC track, intensity, and size are perturbed according to quasi-random Korobov sequences of historical forecast errors with assumed Gaussian and uniform statistical distributions. These perturbations are run in an ensemble of hydrodynamic storm tide model simulations. The resulting set of maximum water surface elevations are dimensionality reduced using Karhunen–Loève expansions and then used as a training set to develop a polynomial chaos (PC) surrogate model from which global sensitivities and probabilistic predictions can be extracted. The maximum water surface elevation is extrapolated over dry points incorporating energy head loss with distance to properly train the surrogate for predicting inundation. We find that the surrogate constructed with third-order PCs using elastic net penalized regression with leave-one-out cross validation provides the most robust fit across training and test sets. Probabilistic predictions of maximum water surface elevation and inundation area by the surrogate model at 48-h lead time for three past U.S. landfalling hurricanes (Irma in 2017, Florence in 2018, and Laura in 2020) are found to be reliable when compared to best track hindcast simulation results, even when trained with as few as 19 samples. The maximum water surface elevation is most sensitive to perpendicular track-offset errors for all three storms. Laura is also highly sensitive to storm size and has the least reliable prediction. Significance Statement The purpose of this study is to develop and evaluate a methodology that can be used to provide high-quality probabilistic predictions of hurricane-induced storm tide and inundation with limited time and resources. This is important for emergency management purposes during or after the landfall of hurricanes. Our results show that sampling forecast errors using quasi-random sequences combined with machine learning techniques that fit polynomial functions to the data are well suited to this task. The polynomial functions also have the benefit of producing exact sensitivity indices of storm tide and inundation to the forecasted hurricane properties such as path, intensity, and size, which can be used for uncertainty estimation. The code implementing the presented methodology is publicly available on GitHub.

54 ENVIRONMENTAL SCIENCES↗

Thermophysical and Optical Properties of Semiconducting Ga2Te3 Melt

The majority of bulk semiconductor single crystals are presently grown from their melts. The thermophysical and optical properties of the melts provide a fundamental understanding of the melt structure and can be used to optimize the growth conditions to obtain higher quality crystals. In this paper, we report several thermophysical and optical properties for Ga2Te3 melts, such as electrical conductivity, viscosity, and optical transmission for temperatures ranging from the melting point up to approximately 990 C. The conductivity and viscosity of the melts are determined using the transient torque technique. The optical transmission of the melts is measured between the wavelengths of 300 and 2000 nm by an dual beam reversed-optics spectrophotometer. The measured properties are in good agreement with the published data. The conductivities indicate that the Ga2Te3 melt is semiconductor-like. The anomalous behavior in the measured properties are used as an indication of a structural transformation in the Ga2Te3 melt and discussed in terms of Eyring's and Bachinskii's predicted behaviors for homogeneous melts.

Li, Chao↗

Orbital stability analysis of hypothetical Earth-mass and Luna-mass moons in the Sagarmatha (HD 100777) star system

Abstract Six of the solar system planets have 150 confirmed moons (Earth—1, Mars—2, Jupiter—53, Saturn—53, Uranus—27, Neptune—14) and seven of them (Ganymede, Titan, Callisto, Io, the Moon, Europa, Triton) have masses >0.001 $M_\oplus$. However, no exomoons have yet been discovered despite the successful detection of ∼5000 exoplanets. We can infer, based on the solar system planets, that these exoplanets are capable of hosting one or more exomoons. In this paper, we study the possible existence of hypothetical Earth-mass and Luna-mass moons orbiting the Jupiter-mass planet, Laligurans (HD 100777b) in the Sagarmatha (HD 100777) star system by means of orbital stability. We apply long-term orbital integrations and the MEGNO (mean exponential growth of nearby orbits) chaos indicator to study the orbital stability of the moons and predict a phase-space region comprising periodic, chaotic, and unstable orbits. The phase spaces primarily constitute the moon's semimajor axis, which extends from the host planet's Roche radius to the Hill radius, and full range of eccentricity. Specific points are picked from three different regions of the MEGNO map and run as single-orbit integration for up to 10 billion periods of the innermost orbit. Furthermore, the lifetime and maximum eccentricity maps are generated from the direct integration to inspect the stable and unstable orbital configurations. The analyses of these maps, with the aid of time-series plots, show that both moons maintain stable orbits in the low-eccentricity regime and semimajor axis between the Roche limit and 28.4% of the Hill radius of the planet.

Gautam, Aayush (ORCID:0000000247257589)↗

Entanglement Structure of Non-Gaussian States and How to Measure It

Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. Here, we present a protocol that constrains quantum states using experimentally measured correlation functions. This method enables measurement of a quantum state’s entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol’s usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest nontrivial expansion quantitatively predicts early time thermalization dynamics, including signaling the onset of quantum chaos indicated by the entanglement Hamiltonian.

Fermi gases↗

The Influence of Planetary Mass on the Dynamical Lifetime of Planetary Systems

Recent numerical and analytic studies of planetary orbits have demonstrated the importance of resonances and chaos in destabilizing planetary systems. Newton's "clockwork" description of regular, predictable planetary orbits has been replaced by a view in which many systems can have long but finite lifetimes. This new knowledge has altered our perceptions of the later stages of planetary growth and of the stability of planetary systems. Stability criteria are inexact and time dependent. Most previous studies have focused on the effects in initial planetary orbits on the stability of the system. We are conducting an investigation which focuses on the dependence of stability criteria on planetary mass. Synthetic systems are created by increasing the masses of the planets in our Solar System or of the moons of a particular planet; these systems are then integrated until orbit crossing occurs. We have found that over some ranges, the time until orbit crossing varies to a good approximation as a power clothe factor by which the masses of the secondaries arc increased; some scatter occurs as a consequence of vie chaotic nature of orbital evolution. The slope of this power law varies substantially from system to system, and for moons it is mildly dependent on the inclusion of the planet's quadrupole moment in the gravitational potential.

Lissauer, J. J.↗

Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy

The increasing use of dynamic atomic force microscopy (d-AFM) for nanoscale materials characterization calls for a deeper understanding of the cantilever dynamics influencing scan stability, predictability, and image quality. Model development is critical to such understanding. Renormalization of the equations governing d- AFM provides a simple interpretation of cantilever dynamics as a single spring and mass system with frequency dependent cantilever stiffness and damping parameters. The renormalized model is sufficiently robust to predict the experimentally observed splitting of the free-space cantilever resonance into multiple resonances upon cantilever-sample contact. Central to the model is the representation of the cantilever sample interaction force as a polynomial expansion with coefficients F(sub ij) (i,j = 0, 1, 2) that account for the effective interaction stiffness parameter, the cantilever-to-sample energy transfer, and the amplitude of cantilever oscillation. Application of the Melnikov method to the model equation is shown to predict a homoclinic bifurcation of the Smale horseshoe type leading to a cascade of period doublings with increasing drive displacement amplitude culminating in chaos and loss of image quality. The threshold value of the drive displacement amplitude necessary to initiate subharmonic generation depends on the acoustic drive frequency, the effective damping coefficient, and the nonlinearity of the cantilever-sample interaction force. For parameter values leading to displacement amplitudes below threshold for homoclinic bifurcation other bifurcation scenarios can occur, some of which lead to chaos.

Cantrell, John H.↗

Mechanistic calculation of the effective silver diffusion coefficient in polycrystalline silicon carbide: Application to silver release in AGR-1 TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission products produced and not retained in the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110m Ag isotope. However, accurately determining the fission gas release rate requires knowing the diffusion coefficient through the SiC layer. In this study, we leverage atomistic calculations of Ag diffusivity in SiC bulk and grain boundaries (GBs) to develop a mesoscale effective Ag diffusion coefficient (D eff ) in SiC. Since GBs serve as pathways for Ag diffusion, D eff is defined as a function of temperature and microstructure variables. In particular, the size of SiC grains in the direction perpendicular to diffusion is shown to significantly affect Ag diffusion. The prediction of the mechanistic, mesoscale approach falls within one order of magnitude of empirical values. The temperature and microstructure-dependent effective Ag diffusivity in SiC is implemented in the fuel performance code Bison with a correction factor to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size on Ag release and improve Bison’s predictions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Light Diffraction by Large Amplitude Ultrasonic Waves in Liquids

Light diffraction from ultrasound, which can be used to investigate nonlinear acoustic phenomena in liquids, is reported for wave amplitudes larger than that typically reported in the literature. Large amplitude waves result in waveform distortion due to the nonlinearity of the medium that generates harmonics and produces asymmetries in the light diffraction pattern. For standing waves with amplitudes above a threshold value, subharmonics are generated in addition to the harmonics and produce additional diffraction orders of the incident light. With increasing drive amplitude above the threshold a cascade of period-doubling subharmonics are generated, terminating in a region characterized by a random, incoherent (chaotic) diffraction pattern. To explain the experimental results a toy model is introduced, which is derived from traveling wave solutions of the nonlinear wave equation corresponding to the fundamental and second harmonic standing waves. The toy model reduces the nonlinear partial differential equation to a mathematically more tractable nonlinear ordinary differential equation. The model predicts the experimentally observed cascade of period-doubling subharmonics terminating in chaos that occurs with increasing drive amplitudes above the threshold value. The calculated threshold amplitude is consistent with the value estimated from the experimental data.

Adler, Laszlo↗

Data-Driven Exploration of Climate Attractor Manifolds For Long-Term Predictability

Focal Area: This white paper responds to Focal Area 3. We seek to gain insight into decadal-scale climate predictability by applying novel manifold-finding probabilistic AI techniques to the complex data produced by Earth System models (ESMs) such as E3SM. The associated portfolio of research activities leverages DOE’s asset mix of HPC platforms, climate expertise, climate simulation codes, and AI expertise. Science Challenge: Climate and climate models are dynamical systems exhibiting properties that are interpretable through chaos theory. The theory contains an important concept that is relevant to multi-decade-scale climate prediction: a chaotic attractor. While the space containing all the possible states of the Earth’s atmosphere and ocean, the possible weather, is large, the realized states tend to stay near the smaller-dimensioned attractor. This behavior is responsible for the “order behind the irregularity” [1] of climate phenomena. Climate change can be thought of as a change in the properties of the attractor, and predicting the climate over years to decades is equivalent to predicting how those properties will change. To date, the attractor has been a useful conceptual tool, but has not been amenable to direct characterization. A new development is the advent of efficient high-dimensional manifold-finding probabilistic AI techniques, which permit a data-driven characterization of the ESM attractor and its probability distribution over weather states. Such a characterization would result in a natural dimensional reduction — a “non-linear Principal Components Analysis (PCA) adapted to climate simulation data” — leading to important advances in scenario-based long-term climate prediction, long-term prediction of water cycle extremes, ESM verification, inter-model comparison, and process model development.

54 ENVIRONMENTAL SCIENCES↗

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben D Phillips↗

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben Phillips↗

Sensitivity Analysis of Irradiated Fueled Experiments using the MOOSE Framework [Slides]

Modeling and simulation (M&S) methods are able to predict uncertainties in experimental parameters (e.g., power and fission density) during irradiation. A shortfall exists in predicting how sensitive some of the parameters will behave during the experimental process. Sensitivity and Uncertainty Quantification (SUQ) is critical in support of qualification and licensing reactor fuels. The application of a method to quantify the uncertainty in these experiments is critical to the prediction of their performance. In this work, we propose the use of a polynomial chaos expansion (PCE) method to quantify the sensitive parameters in these simulations and, in an extension, their experimental surrogates. We propose to perform M&S using PCE uncertainty quantification on a previously irradiated fueled experiment in order to provide a validation case for Griffin and expand its use as a verification and validation (V&V) tool for experiments with a neutronics component. Griffin is an advanced, deterministic neutronics analysis code built using the MOOSE (multiphysics object-oriented simulation environment) framework which can provide state-of-the-art neutronic analysis on M&S of experiments. We will use the stochastic tools module (STM) in MOOSE to provide PCE uncertainty quantification on the proposed experimental setup. Idaho National Laboratory (INL) does not yet have an in-house developed code with V&V approval for experiments performed on-site; this work would provide a necessary addition of support for experiments performed at INL. The Nuclear Regulatory Commission (NRC) has explicitly requested uncertainties in calculated values such as fuel power and burnup, and the development of this capability would benefit the relationship between INL and the NRC.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessment of Arctic and Antarctic Sea Ice Predictability in CMIP5 Decadal Hindcasts

This paper examines the ability of coupled global climate models to predict decadal variability of Arctic and Antarctic sea ice. We analyze decadal hindcasts/predictions of 11 Coupled Model Intercomparison Project Phase 5 (CMIP5) models. Decadal hindcasts exhibit a large multimodel spread in the simulated sea ice extent, with some models deviating significantly from the observations as the predicted ice extent quickly drifts away from the initial constraint. The anomaly correlation analysis between the decadal hindcast and observed sea ice suggests that in the Arctic, for most models, the areas showing significant predictive skill become broader associated with increasing lead times. This area expansion is largely because nearly all the models are capable of predicting the observed decreasing Arctic sea ice cover. Sea ice extent in the North Pacific has better predictive skill than that in the North Atlantic (particularly at a lead time of 3-7 years), but there is a reemerging predictive skill in the North Atlantic at a lead time of 6-8 years. In contrast to the Arctic, Antarctic sea ice decadal hindcasts do not show broad predictive skill at any timescales, and there is no obvious improvement linking the areal extent of significant predictive skill to lead time increase. This might be because nearly all the models predict a retreating Antarctic sea ice cover, opposite to the observations. For the Arctic, the predictive skill of the multi-model ensemble mean outperforms most models and the persistence prediction at longer timescales, which is not the case for the Antarctic. Overall, for the Arctic, initialized decadal hindcasts show improved predictive skill compared to uninitialized simulations, although this improvement is not present in the Antarctic.

sea ice↗

Instabilities in a three-level coherently pumped laser

A theory for a coherently pumped, homogeneously broadened laser is developed which predicts instability at excitations 1.6 times threshold. The system exhibits a period-doubling sequence, chaos, and a period-three window.

Ryan, J. C.↗

Intercomparison of Satellite Derived Gravity Time Series with Inferred Gravity Time Series from TOPEX/POSEIDON Sea Surface Heights and Climatological Model Output

The upcoming GRACE mission promises to open a window on details of the global mass budget that will have remarkable clarity, but it will not directly answer the question of what the state of the Earth's mass budget is over the critical last quarter of the 20th century. To address that problem we must draw upon existing technologies such as SLR, DORIS, and GPS, and climate modeling runs in order to improve our understanding. Analysis of long-period geopotential changes based on SLR and DORIS tracking has shown that addition of post 1996 satellite tracking data has a significant impact on the recovered zonal rates and long-period tides. Interannual effects such as those causing the post 1996 anomalies must be better characterized before refined estimates of the decadal period changes in the geopotential can be derived from the historical database of satellite tracking. A possible cause of this anomaly is variations in ocean mass distribution, perhaps associated with the recent large El Nino/La Nina. In this study, a low-degree spherical harmonic gravity time series derived from satellite tracking is compared with a TOPEX/POSEIDON-derived sea surface height time series. Corrections for atmospheric mass effects, continental hydrology, snowfall accumulation, and ocean steric model predictions will be considered.

Cox, C.↗

Investigation of several proxies to estimate sulfuric acid concentration in volcanic plume conditions

Sulfuric acid (H2SO4) is commonly accepted as a key precursor for atmospheric new particle formation (NPF). However, direct measurements of [H2SO4] remain challenging, thus preventing the determination of this important quantity, and, consequently, a complete understanding of its contribution to the NPF process. Several proxies have been developed to bridge the gaps, but their ability to predict [H2SO4] in very specific conditions such as those encountered in volcanic plumes (including in particular high sulphur dioxide mixing ratios) has not been evaluated so far. In this context, the main objective of the present study was to develop new proxies for daytime [H2SO4] in volcanic plume conditions and compare their performance to that of the proxies available in the literature. In specific, the data collected at Maïdo during the OCTAVE 2018 campaign, in the volcanic eruption plume of the Piton de la Fournaise, were first used to derive seven proxies based on the knowledge of sulphur dioxide (SO2) mixing ratio, global radiation, condensation sink (CS) and relative humidity (RH). In three of the seven proxies (F1–F3), all variables were given equal weight in the prediction of [H2SO4], while adjusted powers were allowed for the different variables in the other four proxies (A1–A4). Proxies A1–A4 were overall found to perform better compared to F1–F3, with, in specific, improved predictive ability for [H2SO4] > 2 × 108 per cu.cm. The CS was observed to play an important role in regulating [H2SO4], while, in contrast, the inclusion of RH did not improve the predictions. A last expression accounting for an additional sink term related to cluster formation, S1, was also tested and showed a very good predictive ability over the whole range of measured [H2SO4]. The newly developed proxies were in a second step further evaluated using airborne measurements performed in the passive degassing plume of Etna during the STRAP 2016 campaign. Increased correlations between observed and predicted [H2SO4] were obtained when the dependence of predicted [H2SO4] over CS was the lowest, and when the dependence over [SO2] was concurrently the highest. The best predictions were finally retrieved by the simple formulation of F2 (in which [SO2] and radiation alone were assumed to explain the variations of [H2SO4] with equal contributions), with a pre factor adapted to the STRAP data. All in all, our results illustrate the fairly good capacity of the proxy available in the literature to describe [H2SO4] in volcanic plume conditions, but highlight at the same time the benefit of the newly developed proxies for the prediction of the highest concentrations ([H2SO4] > 2–3 × 108 per cu.cm). Also, the contrasting behaviours of the new proxies in the two investigated datasets indicate that in volcanic plumes like in other environments, the relevance of a proxy can be affected by changes in environmental conditions, and that location specific coefficients do logically improve the predictions.

new particle formation↗