Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 379 records · Page 21

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗

Uncertainty quantification for Joule heating processes in fibrous pore-resolved media

Joule heating (JH) is an energy-efficient and sustainable technique for heating materials. Its application for industrial heating, particularly, has been gaining attention due to its potential for increasing the yield of various chemical products. The process involves the use of heating elements (materials that are highly conductive electrically and thermally) to heat up other materials or substances. These conductors, however, can exhbit varying degrees of uncertainty due to non-linearities in their temperature-dependent properties, which could result in variable material behavior. In this work, we carry out uncertainty quantification (UQ) at the pore scale to describe the uncertainty of such materials. In so doing, we applied the non-intrusive polynomial chaos expansion (PCE) technique to quantify the uncertainty within the system. The steady state Joule heating equation was solved numerically at the pore scale mimicking conditions within a heating chamber for propane dehydrogenation, and various electro-thermal profiles were obtained. We also examined the effect of the number of sampling points (20 – 100) and order of the PCE coefficients (2 – 5) on the accuracy of the temperature evaluations. The results were then benchmarked with the standard Monte Carlo (MC) method. The average temperature of the 4th-order global PCE showed good agreement with the MC results (which were positively skewed). Orders greater than 4 gave an underestimation of the temperatures while predictions for the peak temperature improved as the number of sampling points increased.

Fagbemi, Samuel [ORNL] (ORCID:0000000236995025)↗

Scaling CFB risers: Maintaining microstructure dynamics

Over the years there have been numerous approaches to develop a scaleup methodology for circulating fluidized beds. Many of these past approaches have relied on macroscopic definitions of similitude to define dimensionless groups such as Reynolds number, Froude number, density ratios, and dimensional ratios involving different length scales. These approaches have demonstrated various levels of success, but often fail when used outside the bounds of experimental datasets. The issue with these methods is that they rely on macroscopic similitude with the hopes that those conditions would provide microscopic similitude as it does in single phase systems where the Buckingham Pi techniques have given reasonable success. Two of the main reasons for failure are 1) changing fluidization regimes and 2) changing material properties that result in the material falling into a different Geldart group. Both are certain to change the microstructure (local instantaneous volume solids fraction) and changes to the microstructure affect the interphase transport properties of the reactor system, namely the mass transfer and heat transfer. Here, this work presents a new approach towards scaling based upon maintaining similitude at the microscale. It reviews a novel dimensionless flow regime map developed from analysis of the deterministic chaos parameters and higher order moments to ensure microstructure similitude. The work then examines four cases from the literature to assess the scaling approach used: (1) same particle-different scale risers for Geldart Group B particles, (2) Glicksman scaling of atmospheric circulating fluidized bed combustor, (3) same particle-different scale risers for Geldart Group A particles and (4) Glicksman scaling for pressurized gasification. General findings are that none of the scaling approaches produced directly applicable microstructure similar behavior based upon the chaotic and higher moment analysis. The positive take away is that the approach can define scale conditions that will provide similar microstructure that can be us to predict transport behavior.

42 ENGINEERING↗

Model correction and updating of a stochastic degradation model for failure prognostics of miter gates

Understanding the degradation of the quoin block is vital for failure prognostics in miter gates. Due to the complicated degradation mechanism, degradation models based on simplifications and assumptions cannot accurately describe the damage evolution. It is observed that small errors in a simplified degradation model can lead to a large discrepancy in the remaining useful life estimation attributed to error accumulation over time. Aiming to address this issue in failure prognostics, this paper presents a dynamic model correction framework for a simplified degradation model using strain measurements. In the proposed framework, a polynomial chaos expansion (PCE) model is employed to compensate the missing physics in a simplified stochastic degradation model. Here, a maximum likelihood estimation method is developed to estimate the uncertain parameters of the simplified physics-based degradation model along with the unknown PCE model parameters using strain measurements as the observables. The updated damage degradation model is then applied to failure prognostics of a miter gate. Results of a case study show that the proposed approach can effectively improve the accuracy of failure prognostics in miter gates.

42 ENGINEERING↗

Time domain probabilistic seismic risk analysis using ground motion prediction equations of Fourier amplitude spectra

Modeling of Fourier amplitude spectra (FAS) of seismic motions has gained much attention in engineering seismology. In the past few years, several ground motion prediction equations (GMPEs) and inter-frequency correlation structure of FAS have been established. Due to many preferable characteristics of FAS, probabilistic seismic hazard/risk analysis is rapidly changing from ergodic, spectrum acceleration Sa(T 0 )-based approach to non-ergodic, site-specific, FAS-based approach. This paper presents time domain intrusive framework for probabilistic seismic risk analysis using GMPE of FAS. Herein, methodology for time domain stochastic ground motion modeling based on GMPEs of FAS is presented in some detail. The simulated uncertain motions are modeled as a random process and represented by polynomial chaos Karhunen-Loève expansion. The random process excitations are further propagated into the uncertain structural system using Galerkin stochastic finite element method (SFEM). Probabilistic evolution of structural response is solved, and such solution is used to develop seismic risk for any damage state. The presented framework is illustrated through seismic risk analysis of a four-story building subjected to possible earthquakes from two strike slip faults. The influences of the epistemic uncertainties in source stress drop Δσ and site attenuation κ0 on seismic risk are investigated. The need for non-ergodic seismic risk analysis with source-specific and site specific characterizations is emphasized.

58 GEOSCIENCES↗

Scattering of radiofrequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radiofrequency (RF) electromagnetic waves that propagate through this region. It is important to know the spectral properties of these scattered RF waves, whether used for diagnostics or for heating and current drive. The spectral changes influence the interpretation of the obtained diagnostic data, and the current and heating profiles. A full-wave, three-dimensional (3-D) electromagnetic code ScaRF has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method used for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma by the anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can also be used for the study of the scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U and DIII-D. For the case of medium-sized tokamaks, there is experimental evidence that drift waves and rippling modes are present in the edge region. Hence, we have studied the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes. The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion method in conjunction with the Smolyak sparse-grid integration, which is known as the PCE-SG method. The PCE-SG method is proven to be accurate and more efficient compared with alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Role of Subsurface Flow on Evapotranspiration: A Global Sensitivity Analysis

Abstract Water resources are impacted by water‐energy balance fluxes at the land surface, most notably evapotranspiration ( ET ), the largest component of surface energy balance. While integrated hydrologic models show promise in quantifying the nonlinear dynamics at this interface, the model results are plagued by parametric uncertainties. Given the high computational demand of running multiple parameter spaces of these models, little is known about how these uncertainties propagate into land surface processes. In this work, we perform a global sensitivity analysis by computing Sobol and AMAE indices using a surrogate model constructed with a polynomial chaos expansion to assess the impacts of the subsurface physical properties on ET . We do so by modeling a semisynthetic test case. Our results show that the effects of vertical hydraulic conductivity, porosity, and the water retention curve parameter Van Genuchten α mainly control ET . However, we note that while evaporation ( E ) shows behavior similar to the ET with high sensitivity to the parameters controlling the flow in the unsaturated zone, transpiration ( T ) is very sensitive to the saturated zone parameters and groundwater flow, especially during periods without rain. Our results also show that heterogeneities in land cover impact the flow in the saturated zone. Furthermore, our work demonstrates that reduced‐order models can be developed to make integrated models more accessible for rigorous sensitivity analyses and calibration purposes.

Maina, Fadji Zaouna↗

Coherent Streamflow Variability in Monsoon Asia Over the Past Eight Centuries—Links to Oceanic Drivers

The Monsoon Asia region is home to ten of the world’s biggest rivers, supporting the lives of 1.7 billion people who rely on streamflow for water, energy, and food. Yet, a synoptic understanding of multi-centennial streamflow variability for this region is lacking. Here we produce the first large scale streamflow reconstruction over Monsoon Asia (63 stations in 16 countries), using a novel climate-informed dynamic algorithm that is skillful over 92% of the gauging stations. In this work, we show that streamflow in Monsoon Asia is spatially coherent, owing to common drivers from the Pacific, Indian, and Atlantic Oceans that exert their greatest influence over the Mekong and Chao Phraya basins. This work increases our understanding of streamflow variability over Monsoon Asia. We suggest that future water management in the region should be coordinated between basins, taking into account the states of the oceans.

54 ENVIRONMENTAL SCIENCES↗

Autocalibration of the E3SM Version 2 Atmosphere Model Using a PCA-Based Surrogate for Spatial Fields

Global Climate Model tuning (calibration) is a tedious and time-consuming process, with high-dimensional input and output fields. Experts typically tune by iteratively running climate simulations with hand-picked values of tuning parameters. Many, in both the statistical and climate literature, have proposed alternative calibration methods, but most are impractical or difficult to implement. We present a practical, robust, and rigorous calibration approach on the atmosphere-only model of the Department of Energy's Energy Exascale Earth System Model (E3SM) version 2. Our approach can be summarized into two main parts: (a) the training of a surrogate that predicts E3SM output in a fraction of the time compared to running E3SM, and (b) gradient-based parameter optimization. To train the surrogate, we generate a set of designed ensemble runs that span our input parameter space and use polynomial chaos expansions on a reduced output space to fit the E3SM output. We use this surrogate in an optimization scheme to identify values of the input parameters for which our model best matches gridded spatial fields of climate observations. To validate our choice of parameters, we run E3SMv2 with the optimal parameter values and compare prediction results to expertly-tuned simulations across 45 different output fields. This flexible, robust, and automated approach is straightforward to implement, and we demonstrate that the resulting model output matches present day climate observations as well or better than the corresponding output from expert tuned parameter values, while considering high-dimensional output and operating in a fraction of the time.

54 ENVIRONMENTAL SCIENCES↗

Quantum computing universal thermalization dynamics in a (2 + 1)D Lattice Gauge Theory

Simulating non-equilibrium phenomena in strongly-interacting quantum many-body systems, including thermalization, is a promising application of near-term and future quantum computation. By performing experiments on a digital quantum computer consisting of fully-connected optically-controlled trapped ions, we study the role of entanglement in the thermalization dynamics of a Z 2 lattice gauge theory in 2+1 spacetime dimensions. Using randomized-measurement protocols, we efficiently learn a classical approximation of non-equilibrium states that yields the gap-ratio distribution and the spectral form factor of the entanglement Hamiltonian. These observables exhibit universal early-time signals for quantum chaos, a prerequisite for thermalization. Our work, therefore, establishes quantum computers as robust tools for studying universal features of thermalization in complex many-body systems, including in gauge theories.

97 MATHEMATICS AND COMPUTING↗

Turnkey locking of quantum-dot lasers directly grown on Si

Ultralow-noise laser sources are crucial for a variety of applications, including microwave synthesizers, optical gyroscopes and the manipulation of quantum systems. Silicon photonics has emerged as a promising solution for high-coherence applications due to its ability to reduce the system size, weight, power consumption and cost. Semiconductor lasers based on self-injection locking have achieved fibre laser coherence, but typically require a high-quality-factor external cavity to suppress coherence collapse through frequency-selective feedback. Lasers based on external-cavity locking are a low-cost and turnkey operation option, but their coherence is generally inferior to self-injection locking lasers. In this work, we demonstrate quantum-dot lasers grown directly on Si that achieve self-injection-locking laser coherence under turnkey external-cavity locking. The high-performance quantum-dot laser offers a scalable and low-cost heteroepitaxial integration platform. Moreover, the chaos-free nature of the quantum-dot laser enables a 16 Hz Lorentzian linewidth under external-cavity locking using a low-quality-factor external cavity, and improves the frequency noise by an additional order of magnitude compared with conventional quantum-well lasers.

Dong, Bozhang↗

Axon-like active signal transmission

Any electrical signal propagating in a metallic conductor loses amplitude due to the natural resistance of the metal. Compensating for such losses presently requires repeatedly breaking the conductor and interposing amplifiers that consume and regenerate the signal. This century-old primitive severely constrains the design and performance of modern interconnect-dense chips. Here we present a fundamentally different primitive based on semi-stable edge of chaos (EOC), a long-theorized but experimentally elusive regime that underlies active (self-amplifying) transmission in biological axons. By electrically accessing the spin crossover in LaCoO 3 , we isolate semi-stable EOC, characterized by small-signal negative resistance and amplification of perturbations. In a metallic line atop a medium biased at EOC, a signal input at one end exits the other end amplified, without passing through a separate amplifying component. While superficially resembling superconductivity, active transmission offers controllably amplified time-varying small-signal propagation at normal temperature and pressure, but requires an electrically energized EOC medium. Operando thermal mapping reveals the mechanism of amplification—bias energy of the EOC medium, instead of fully dissipating as heat, is partly used to amplify signals in the metallic line, thereby enabling spatially continuous active transmission, which could transform the design and performance of complex electronic chips.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Extended harmonic mapping connects the equations in classical, statistical, fluid, quantum physics and general relativity

Abstract One potential pathway to find an ultimate rule governing our universe is to hunt for a connection among the fundamental equations in physics. Recently, Ren et al. reported that the harmonic maps with potential introduced by Duan, named extended harmonic mapping (EHM), connect the equations of general relativity, chaos and quantum mechanics via a universal geodesic equation. The equation, expressed as Euler–Lagrange equations on the Riemannian manifold, was obtained from the principle of least action. Here, we further demonstrate that more than ten fundamental equations, including that of classical mechanics, fluid physics, statistical physics, astrophysics, quantum physics and general relativity, can be connected by the same universal geodesic equation. The connection sketches a family tree of the physics equations, and their intrinsic connections reflect an alternative ultimate rule of our universe, i.e. , the principle of least action on a Finsler manifold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of multipacting threshold sensitivity to the random distributions of the secondary electron yield parameters

The way multipacting develops, depends strongly on the secondary emission property of the surface material. The knowledge of secondary electron yield is crucial for accurate prediction of the multipacting threshold. Variations in secondary electron yield parameters from experimental measurements create uncertainty, stemming from handling and surface preparation, and these uncertainties significantly affect multipacting threshold predictions. Despite their significance, the previous studies on the multipacting phenomenon did not adequately address the effect of an assumed random distribution of the secondary emission parameters on the multipacting threshold. Therefore, this paper aims to provide a comprehensive statistical study on how the different random distributions of the secondary emission parameters and, as a result, the uncertainty in the secondary electron yield affect multipacting thresholds. We focus on three commonly used distributions, namely uniform, normal, and truncated normal distributions, to define the uncertainty of random inputs. We use the chaos polynomial expansion method to determine how much each of the random parameters contributes to the multipacting threshold uncertainty. Additionally, we calculate Sobol sensitivity indices to evaluate the impact of the individual parameters or groups of parameters on the model outputs and study how different random distributions of these parameters affected the Sobol index results.

physics↗

Quantum many-body physics from a gravitational lens

The past two decades have seen the emergence of remarkable interconnections among previously remotely related disciplines, such as condensed matter, nuclear physics, gravity and quantum information, fuelled both by experimental advances and by the new powerful theoretical methods offered by holographic duality. In this Review, we sample some recent developments in holographic duality in connection with quantum many-body dynamics. These include insights into strongly correlated phases without quasiparticles and their transport properties, quantum many-body chaos and the scrambling of quantum information. We also discuss recent progress in understanding the structure of holographic duality itself using quantum information, including a ‘local’ version of the duality, as well as the quantum error-correction interpretation of quantum many-body states with a gravity dual, and how such notions help to demonstrate the unitarity of black hole evaporation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting nonresonant pressure-driven MHD modes in equilibria with low magnetic shear

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q 0 . A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. Finally, this work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural inference of networked dynamical systems with universal differential equations

Networked dynamical systems are common throughout science in engineering; e.g., biological networks, reaction networks, power systems, and the like. For many such systems, nonlinearity drives populations of identical (or near-identical) units to exhibit a wide range of nontrivial behaviors, such as the emergence of coherent structures (e.g., waves and patterns) or otherwise notable dynamics (e.g., synchrony and chaos). Here, we seek to infer (i) the intrinsic physics of a base unit of a population, (ii) the underlying graphical structure shared between units, and (iii) the coupling physics of a given networked dynamical system given observations of nodal states. These tasks are formulated around the notion of the Universal Differential Equation, whereby unknown dynamical systems can be approximated with neural networks, mathematical terms known a priori (albeit with unknown parameterizations), or combinations of the two. We demonstrate the value of these inference tasks by investigating not only future state predictions but also the inference of system behavior on varied network topologies. The effectiveness and utility of these methods are shown with their application to canonical networked nonlinear coupled oscillators.

97 MATHEMATICS AND COMPUTING↗