SEARCH · Search NASA

Results for “oscillation cases”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Effect of Hydrogen on Steady-State and Transient Combustion Instability Characteristics

This paper examines the effects of steady-state and transient hydrogen enrichment on thermoacoustic instability in a model gas turbine combustor. Combustion instability, a feedback loop between flame heat release rate oscillations and combustor acoustics, is characterized in a swirl-stabilized flame operated at a range of hydrogen–natural gas fuel blends and heat rates. Measurements of combustor chamber pressure fluctuations and CH* chemiluminescence imaging are used to characterize instability at a range of operating conditions. Steady-state tests show that both mixture heat rate and hydrogen content affect system stability. At a given heat rate, higher levels of hydrogen result in unstable combustion. As heat rate increases, instability occurs at lower concentrations of hydrogen in the fuel. Transient operation was tested in two directions—instability onset and decay—and two hydrogen-addition times—a short time of 1 ms and a longer time of 4 s. Results show that instability onset processes, through the transient addition of hydrogen, are highly repeatable regardless of the timescale of hydrogen addition. Certain instability decay processes are less repeatable, resulting in cases that do not fully transition from unstable to stable combustion despite similar changes in hydrogen fuel flow rate. Flame behavior before, during, and after the transient is characterized using high-speed CH* chemiluminescence imaging. Analysis of the high-speed images shows changes in flame stabilization and dynamics during the onset and decay processes. Finally, the results of this study can have implications for systems that experience variations in fuel composition, particularly in light of growing interest in hydrogen as a renewable fuel.

42 ENGINEERING↗

DOE OSTI2025

Ergodic Lagrangian dynamics in a superhero universe

We present a fictional scenario that, while undeniably whimsical, provides the foundation for a unique exercise in extended problem solving, physics analysis, and quantitative model development. Starting with the foundational premise of the Wild Cards shared-world superhero universe, we demonstrate how a variety of concepts appropriate to the advanced undergraduate level—ergodicity, functional analysis, Lagrangian mechanics, and the ever-important simplifying approximation—can be combined into a rich, coherent mathematical model. The goal of this case study is to develop a useful pedagogical exercise in exploring an open-ended research question that presents, at first glance, no clear path forward. Being both eclectic and lengthy, this exercise offers a unique way for students to apply their core physics and mathematics education. It is perhaps best used within a senior honors seminar or within a brief (e.g., January term) elective class.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Investigations of Long-Range Wakefield Effects in a TESLA-type Cryomodule at FAST

The preservation of low emittance of electron beams during transport in the accelerating structures of large facilities is an ongoing challenge. In the cases of the TESLA-type superconducting rf cavities currently used in the European X-ray Free-electron Laser (XFEL) and the under-construction Linac Coherent Light Source upgrade (LCLS-II), off-axis beam transport may result in emittance dilution due to transverse long-range wakefields (LRWs) and short-range wakefields (SRW)***. To investigate such effects, experiments were performed at the Fermilab Accelerator Science and Technology (FAST) facility with its unique configuration of two TESLA-type cavities after the photocathode rf gun followed by an 8-cavity cryomodule CM). We generated beam trajectory changes with the H/V125 corrector set located 4 m upstream of the cryomodule. At 125 pC/bunch, 50 bunches, 25-MeV input, and 100-MeV exit energy, we observed for the first time submacropulse position slews of up to 500 microns at locations ~3 m after the CM and a centroid oscillation at a difference frequency of 240 kHz further downstream. Both are emittance-dilution effects which we mitigated with selective upstream beam steering.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Numerical simulation of parametric resonance in point absorbers using a simplified model

Abstract Parametric resonance is a non‐linear phenomenon in which a system can oscillate at a frequency different from its exciting frequency. Some wave energy converters are prone to this phenomenon, which is usually detrimental to their performance. Here, a computationally efficient way of simulating parametric resonance in point absorbers is presented. The model is based on linear potential theory, so the wave forces are evaluated at the mean position of the body. However, the first‐order variation of the body's centres of gravity and buoyancy is taken into account. This gives essentially the same result as a more rigorous approach of keeping terms in the equation of motion up to second order in the body motions. The only difference from a linear model is the presence of non‐zero off‐diagonal elements in the mass matrix. The model is benchmarked against state‐of‐the‐art non‐linear Froude–Krylov and computational fluid dynamics models for free decay, regular wave, and focused wave group cases. It is shown that the simplified model is able to simulate parametric resonance in pitch to a reasonable accuracy even though no non‐linear wave forces are included. The simulation speed on a standard computer is up to two orders of magnitude faster than real time.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

The guiding center and gyrokinetic theory of magnetized particle motion is extended to the regime of large electric field gradients perpendicular to the magnetic field. A gradient in the electric field directly modifies the oscillation frequency and causes the Larmor orbits to deform from circular to elliptical trajectories. In order to retain a good adiabatic invariant, there can only be strong dependence on a single coordinate at lowest order, so that resonances do not generate chaotic motion that destroys the invariant. When the gradient across magnetic flux surfaces is dominant, the guiding center drift velocity becomes anisotropic in response to external forces and additional curvature drifts must be included. The electric polarization density remains gyrotropic, but both the polarization and magnetization are modified by the change in gyrofrequency. The theory can be applied to shear flows that are even stronger than those observed in the edge transport barrier of a high-performance tokamak (H-mode) pedestal, even if the toroidal field is as small as or even smaller than the poloidal field. Yet, the theory retains a mathematical form that is similar to the standard case and can readily be implemented within existing simulation tools.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Understanding Low Level Jets in the US Atlantic Offshore

Low level jets (LLJs) in the atmosphere exhibit a local windspeed maximum in the boundary layer, with positive shear beneath the jet and negative shear above the jet. Wind turbines tend to experience increased loads with varying wake recovery characteristics in the presence of an LLJ, therefore understanding the mechanism and impact of LLJs is crucial to wind energy development. The US Mid-Atlantic offshore region is a huge potential wind energy resource, yet LLJs in this area are poorly understood. In particular, the coastal offshore environment does not exhibit the same diurnal cycle that leads to strong LLJs in the well-characterized Great Plains region. In this study, we use the Weather Research and Forecasting model (WRF) plus lidar buoy data to identify case studies for LLJ events in 2020 in the New York Bight. A reduced order model is presented to explain the onset of these events based on the competing effects of baroclinicity and eddy diffusivity. Finally, using the macro-scale WRF, we drive a micro-scale large eddy simulation (LES) to generate a more detailed characterization of the Marine Boundary Layer during an LLJ. Gaining a better understanding of LLJs and their impacts on offshore wind in the mid-Atlantic is crucial for a transition toward renewable energy.

atmospheric boundary layer↗

DOE OSTI2023

Electron-field instability: Excitation of electron plasma waves by an electric field

Electric fields are commonplace in plasmas and affect transport by driving currents and, in some cases, instabilities. The necessary condition for instability in collisionless plasmas is commonly understood to be described by the Penrose criterion, which quantifies a sufficient relative drift between different populations of particles that must be present for wave amplification via inverse Landau damping. For example, electric fields generate drifts between electrons and ions that can excite the ion-acoustic instability. Here, we use particle-in-cell simulations and linear stability analysis to show that the electric field can drive a fundamentally different type of kinetic instability, named the electron-field instability. This instability excites electron plasma waves with wavelengths ≳30λ De , has a growth rate that is proportional to the electric field strength, and does not require a relative drift between electrons and ions. The Penrose criterion does not apply when accounting for the electric field. Furthermore, the large value of the observed frequency, near the electron plasma frequency, further distinguishes it from the standard ion-acoustic instability, which oscillates near the ion plasma frequency. The ubiquity of macroscopic electric fields in quasineutral plasmas suggests that this instability is possible in a host of systems, including low-temperature and space plasmas. In fact, damping from neutral collisions in such systems is often not enough to completely damp the instability, adding to the robustness of the instability across plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Mountain waves can impact wind power generation

Abstract. Mountains can modify the weather downstream of the terrain. In particular, when stably stratified air ascends a mountain barrier, buoyancy perturbations develop. These perturbations can trigger mountain waves downstream of the mountains that can reach deep into the atmospheric boundary layer where wind turbines operate. Several such cases of mountain waves occurred during the Second Wind Forecast Improvement Project (WFIP2) in the Columbia River basin in the lee of the Cascade Range bounding the states of Washington and Oregon in the Pacific Northwest of the United States. Signals from the mountain waves appear in boundary layer sodar and lidar observations as well as in nacelle wind speeds and power observations from wind plants. Weather Research and Forecasting (WRF) model simulations also produce mountain waves and are compared to satellite, lidar, and sodar observations. Simulated mountain wave wavelengths and wave propagation speeds (group velocities) are analyzed using the fast Fourier transform. We found that not all mountain waves exhibit the same speed and conclude that the speed of propagation, magnitudes of wind speeds, or wavelengths are important parameters for forecasters to recognize the risk for mountain waves and associated large drops or surges in power. When analyzing wind farm power output and nacelle wind speeds, we found that even small oscillations in wind speed caused by mountain waves can induce oscillations between full-rated power of a wind farm and half of the power output, depending on the position of the mountain wave's crests and troughs. For the wind plant analyzed in this paper, mountain-wave-induced fluctuations translate to approximately 11 % of the total wind farm output being influenced by mountain waves. Oscillations in measured wind speeds agree well with WRF simulations in timing and magnitude. We conclude that mountain waves can impact wind turbine and wind farm power output and, therefore, should be considered in complex terrain when designing, building, and forecasting for wind farms.

17 WIND ENERGY↗

DOE OSTI2024

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus ( T 2 ) and associated systems on the plane ( R 2 ), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

General relativistic pulsations of ultra-massive ZZ Ceti stars

ABSTRACT Ultra-massive white dwarf stars are currently being discovered at a considerable rate, thanks to surveys such as the Gaia space mission. These dense and compact stellar remnants likely play a major role in Type Ia supernova explosions. It is possible to probe the interiors of ultra-massive white dwarfs through asteroseismology. In the case of the most massive white dwarfs, general relativity could affect their structure and pulsations substantially. In this work, we present results of relativistic pulsation calculations employing relativistic ultra-massive ONe-core white dwarf models with hydrogen-rich atmospheres and masses ranging from 1.29 to $1.369 \ \mathrm{M}_{\odot }$ with the aim of assessing the impact of general relativity on the adiabatic gravity (g)-mode period spectrum of very high mass ZZ Ceti stars. Employing the relativistic Cowling approximation for the pulsation analysis, we find that the critical buoyancy (Brunt–Väisälä) and acoustic (Lamb) frequencies are larger for the relativistic case, compared to the Newtonian case, due to the relativistic white dwarf models having smaller radii and higher gravities for a fixed stellar mass. In addition, the g-mode periods are shorter in the relativistic case than those in the Newtonian computations, with relative differences of up to ∼$50$ per cent for the highest mass models ($1.369 \ \mathrm{M}_{\odot }$) and for effective temperatures typical of the ZZ Ceti instability strip. Hence, the effects of general relativity on the structure, evolution, and pulsations of white dwarfs with masses larger than ∼$1.29 \ \mathrm{M}_{\odot }$ cannot be ignored in the asteroseismological analysis of ultra-massive ZZ Ceti stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Gate-Based Quantum Simulation of Gaussian Bosonic Circuits on Exponentially Many Modes

We introduce a framework for simulating, on an ( n + 1 )-qubit quantum computer, the action of a Gaussian bosonic (GB) circuit on a state over 2 n modes. Specifically, we encode the initial bosonic state’s expectation values over quadrature operators (and their covariance matrix) as an input qubit state. This is then evolved by a quantum circuit that effectively implements the symplectic propagators induced by the GB gates. We find families of GB circuits and initial states leading to efficient quantum simulations. For this purpose, we introduce a dictionary that maps between GB and qubit gates such that particle- (non-particle-) preserving GB gates lead to real- (imaginary-) time evolutions at the qubit level. For the special case of particle-preserving circuits, we present a bounded-error-quantum-polynomial time (BQP)-complete GB decision problem, indicating that GB evolutions of Gaussian states on exponentially many modes are as powerful as universal quantum computers. We also perform numerical simulations of an interferometer on ∼ 8 × 10 9 modes, illustrating the power of our framework. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Investigation of oscillations above the acoustic cut-off frequency for a sample of Kepler stars. II. Identifying temporal frequency shifts

ABSTRACT Stellar oscillations with frequencies greater than the acoustic cut-off frequency are not trapped within the star’s interior. Geometric interference between these high-frequency waves produces a peak-like structure in the power spectrum of the star, known as the pseudo-modes. In this study, we aim to analyse 51 stars previously identified (from a sample of 91 Kepler stars) to contain statistically significant evidence for pseudo-modes, to determine if the pseudo-mode frequencies vary in time. We obtained these frequency variations and the uncertainty by the resampled periodogram approach, where periodogram realisations were created from successive, overlapping 90-d time segments, and a cross-correlation function was generated between them. The computed shifts were compared against temporal p-mode frequency shifts and a stellar magnetic activity proxy, $S_{\mathrm{ph}}$. We found that for 16 per cent of the stars pseudo-mode frequency shifts were significantly anticorrelated with p-mode shifts, as is the case for the Sun. However, we also found pseudo-mode and p-mode shifts to be significantly correlated in-phase in 8 per cent of our sample. The remainder showed no significant correlation. We also searched for trends between the maximum variation of pseudo-mode frequency shifts and stellar parameters. Whilst we found no correlation between the maximum variation and stellar photospheric and chromospheric proxies, we found frequency variations to have a greater amplitude for cooler and slower rotating stars, which is in opposition to p-mode frequency shift variation. Temporal changes to pseudo-mode frequencies may be attributable to variations in stellar magnetic activity, and therefore can contribute to our understanding of stellar activity cycles.

Millson, L. J. (ORCID:0009000342542676)↗

DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

On the Regionality of Moist Kelvin Waves and the MJO: The Critical Role of the Background Zonal Flow

Abstract A global model with superparameterized physics is used to shed light on the observed regionality of convectively coupled Kelvin waves and the Madden‐Julian Oscillation (MJO). A series of aquaplanet simulations over zonally uniform sea‐surface temperatures is performed, in which the axisymmetric structure of the background zonal flow is altered through nudging, while maintaining a quasi‐fixed rainfall climatology. Results show that nudging at the equator to match profiles typical of the Indo‐Pacific or eastern Pacific sectors yields eastward‐moving tropical rain spectra typical of those sectors. Two different mechanistic pathways are identified as being responsible for this mean‐flow dependence, in addition to Doppler shifting effects. The first is through shifts of the Rossby wave critical line in the subtropical upper troposphere that affect the lateral forcing of Kelvin‐mode circulations at the equator by eastward and equatorward‐propagating eddies impinging on the tropics from higher latitudes. The second is through changes in the strength of the mean cyclonic shear in the lower tropical troposphere that affect the degree to which intraseasonal fluctuations in Kelvin‐mode zonal winds modulate the activity of higher‐frequency equatorial Rossby‐type eddies. In cases where the mean low‐level cyclonic shear is enhanced, the strength of this modulation, referred to as “shear‐induced eddy modulation” or SIEM, is also seen to be enhanced, such that MJO‐like modes of variability are rendered either unstable or near neutral, depending on the strength of the shear.

Tulich, S. N.↗

DOE OSTI2021

Using Disorder to Overcome Disorder: A Mechanism for Frequency and Phase Synchronization of Diode Laser Arrays

Noise and disorder are known, in certain circumstances and for certain systems, to improve the level of coherence over that of the noise-free system. Examples include cases in which disorder enhances response to periodic signals, and those where it suppresses chaotic behavior. We report a new type of disorder-enhancing mechanism, observed in a model that describes the dynamics of external cavity-coupled semiconductor laser arrays, where disorder of one type mitigates (and overcomes) the desynchronization effects due to a different disorder source. Here, we demonstrate stabilization of dynamical states due to frequency locking and subsequently frequency locking-induced phase locking. Here, we have reduced the equations to a potential model that illustrates the mechanism behind the misalignment-induced frequency and phase synchronization.

42 ENGINEERING↗

DOE OSTI2021

Kinetic Measurements in Heterogeneous Catalysis

This contribution is about the experimental determination of the rate of a heterogeneous catalytic reaction and the analysis of kinetic data. In this case, the reaction rate can be defined as the frequency at which the closed sequence of elementary steps transforming reactants into products is occurring. However, the rate of chemical reaction is not directly observed; rather, one records the rate of substance change. The rate of chemical reaction is calculated based on the rate of substance change and assumed stoichiometry of the reaction. In order to make this an intensive quantity (i.e., independent of the volume in which the reaction is performed, this frequency is divided by the latter or by a quantity proportional to the latter). As the focus of what follows is on heterogeneously catalyzed reactions this can be the catalyst mass, the catalyst surface area or the total number of active sites present in the reaction volume confined by the walls of a chemical reactor. For example, the rate may be reported as one of the following: Net rate of consumption/production of component i mol m-3 s-1 Net specific rate of consumption/production of component i mol kg-1 s-1 Rate of substance change (per unit volume of catalyst) mol m-3 s-1 Specific reaction rate (per unit mass of catalyst) mol kg-1 s-1 Net rate of consumption/production of component i mol m-2 s-1 The turnover frequency (TOF) is a characteristic originally introduced by Boudart [1] is expressed as TOF=R/G_tot , where R is the steady-state rate of reactant consumption or product generation and G_totis the total areal density of active sites [mol m-2]. The latter is typically taken from experimental chemisorption data obtained at low temperature. The advancement of a chemical reaction leads to changes in the amounts of reactants and products which, for a stoichiometric single reaction, are connected by the stoichiometric coefficients. In many cases, for the sake of simplicity, we will assume that the investigated reaction proceeds according to a single reaction pathway, so that the net production rate of a reaction component is directly proportional to the reaction rate and that the latter equals the net production rate of any of the involved components when divided by the appropriate stoichiometric coefficient. For complex stoichiometric reactions the net production rates of the involved components are linear combinations of the reaction rates. Different goals for kinetic measurements can be formulated: Catalyst testing, i.e., obtaining the kinetic dependences for the development and discrimination of efficient catalytic materials. Precise kinetic characterization of active materials via high-throughput screening is a part of such activity. Detailed kinetics, i.e., revealing the detailed mechanism of complex catalytic reaction via systematic kinetic studies, both steady-state and non-steady-state Industrial kinetics, i.e., obtaining data and relationships for describing and predicting the behaviour of catalytic reactors and processes at industrial scale Finally, mathematical modeling and analysis where kinetic measurements provide data for understanding complex kinetic phenomena, e.g., oscillations, non-linear self-organization, etc.

36 MATERIALS SCIENCE↗

DOE OSTI2026

SPT-3G D1: Axion early dark energy with CMB experiments and DESI observations

We present the most up-to-date constraints on axion early dark energy (AEDE) from cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) measurements. In particular, we assess the impact of data from ground-based CMB experiments, the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT)—both with and without Planck—on constraints on AEDE. We also highlight the impact that BAO information from the Dark Energy Spectroscopic Instrument (DESI) has on these constraints. From CMB data alone, we do not find statistically significant evidence for the presence of AEDE, and we find only moderate reduction in the Hubble tension. From the latest SPT data alone, we find the maximal fractional contribution of AEDE to the cosmic energy budget is f EDE < 0.12 at 95% confidence level (CL), and the Hubble tension between the SPT and Supernovae, H0, for the Equation of State of Dark Energy (SH0ES) results is reduced to the 2.3 σ level. When combining the latest SPT, ACT, and Planck datasets, we find f EDE < 0.070 at 95% CL and the Hubble tension at the 3.6 σ level. In contrast, adding DESI data to the CMB datasets results in mild preference for AEDE and, in some cases, non-negligible reduction in the Hubble tension. From SPT + DESI , we find f EDE = 0.08 1 - 0.052 + 0.037 at 68% CL, and the Hubble tension reduces to 1.5 σ . From the combination of DESI with all three CMB experiments, we get f EDE = 0.05 5 - 0.047 + 0.024 at 68% CL and a weak preference for AEDE over Λ CDM . This data combination, in turn, reduces the Hubble tension to 2.6 σ . We highlight that this shift in parameters when adding the DESI dataset is a manifestation of the discrepancy currently present between DESI and CMB experiments in the concordance model Λ CDM .

Khalife, A. R. [Paris, Inst. Astrophys.] (ORCID:00↗

DOE OSTI2020

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗