Search NASA⌕ Search

SEARCH · Search NASA

Results for “timescale separation”

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 19 records

Resilience of the slow component in timescale separated synchronized oscillators

Physiological networks are usually made of a large number of biological oscillators evolving on a multitude of different timescales. Phase oscillators are particularly useful in the modelling of the synchronization dynamics of such systems. If the coupling is strong enough compared to the heterogeneity of the internal parameters, synchronized states might emerge where phase oscillators start to behave coherently. Here, we focus on the case where synchronized oscillators are divided into a fast and a slow component so that the two subsets evolve on separated timescales. We assess the resilience of the slow component by, first, reducing the dynamics of the fast one using Mori-Zwanzig formalism. Second, we evaluate the variance of the phase deviations when the oscillators in the two components are subject to noise with possibly distinct correlation times. From the general expression for the variance, we consider specific network structures and show how the noise transmission between the fast and slow components is affected. Interestingly, we find that oscillators that are among the most robust when there is only a single timescale, might become the most vulnerable when the system undergoes a timescale separation. We also find that layered networks seem to be insensitive to such timescale separations.

97 MATHEMATICS AND COMPUTING↗

Resilience of the slow component in timescale-separated synchronized oscillators

Physiological networks are usually made of a large number of biological oscillators evolving on a multitude of different timescales. Phase oscillators are particularly useful in the modelling of the synchronization dynamics of such systems. If the coupling is strong enough compared to the heterogeneity of the internal parameters, synchronized states might emerge where phase oscillators start to behave coherently. Here, we focus on the case where synchronized oscillators are divided into a fast and a slow component so that the two subsets evolve on separated timescales. We assess the resilience of the slow component by, first, reducing the dynamics of the fast one using Mori-Zwanzig formalism. Second, we evaluate the variance of the phase deviations when the oscillators in the two components are subject to noise with possibly distinct correlation times. From the general expression for the variance, we consider specific network structures and show how the noise transmission between the fast and slow components is affected. Interestingly, we find that oscillators that are among the most robust when there is only a single timescale, might become the most vulnerable when the system undergoes a timescale separation. We also find that layered networks seem to be insensitive to such timescale separations.

59 BASIC BIOLOGICAL SCIENCES↗

Distributed Barrier Certificates for Safe Operation of Inverter-Based Microgrids

Inverter-interfaced microgrids differ from the traditional power systems due to their lack of inertia. Vanishing timescale separation between voltage and frequency dynamics makes it critical that faster-timescale stabilizing control laws also guarantee by-construction the satisfaction of voltage limits during transients. In this article, we apply a barrier functions method to compute distributed active and reactive power setpoint control laws that certify satisfaction of voltage limits during transients. Using sum-of-squares optimization tools, we propose an algorithmic construction of these control laws. Numerical simulations are provided to illustrate the proposed method.

Kundu, Soumya↗

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗

Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches

Simulating optical spectra in the condensed phase remains a challenge for theory due to the need to capture spectral signatures arising from anharmonicity and dynamical effects, such as vibronic progressions and their induced asymmetry. As such, numerous simulation methods have been developed that invoke different approximations and vary in their ability to capture different physical regimes. In this work, we use several models of chromophores in the condensed phase and ab initio molecular dynamics simulations to rigorously assess the applicability of methods to simulate optical absorption spectra. Specifically, we focus on the ensemble scheme, which can address anharmonic potential energy surfaces but relies on the applicability of extreme nuclear-electronic timescale separation; the Franck-Condon method, which includes dynamical effects but only at the harmonic level; as well as the recently introduced ensemble zero-temperature Franck-Condon approach, which straddles these limits. We also devote particular attention to the performance of methods derived from a cumulant expansion of the energy gap fluctuations and test the ability to approximate the requisite time correlation functions using classical dynamics with quantum correction factors. These results provide insights as to when these methods are applicable and able to capture the features of condensed phase spectra qualitatively and, in some cases, quantitatively across a range of regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations

While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute–solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.

Chemistry↗

Voltage positioning using co-optimization of controllable grid assets in radial networks

With increasing penetration of solar PV, some distribution feeders are experiencing highly variable net-load flows and even reverse flows. To optimize distribution systems under such conditions, the scheduling of mechanical devices, such as OLTCs and capacitor banks, needs to take into account forecasted solar PV and actual grid conditions. However, these legacy switching assets are operated on a daily or hourly timescale, due to the wear and tear associated with mechanical switching, which makes them unsuitable for real-time control. Therefore, there is a natural timescale- separation between these slower mechanical assets and the responsive nature of inverter-based resources. In this paper, we present a network admissible convex formulation for holistically scheduling controllable grid assets to position voltage optimally against solar PV. An optimal hourly schedule is presented that utilizes mechanical resources to position the predicted voltages close to nominal values, while minimizing the use of inverter-based resources (i.e., DERs), making them available for control at a faster time-scale (after the uncertainty reveals itself). A convex, inner approximation of the OPF problem is adapted to a mixed-integer linear program that minimizes voltage deviations from nominal (i.e., maximizes voltage margins). Here, the resulting OPF solution respects all the network constraints and is, hence, robust against modeling simplifications. Simulation based analysis on IEEE distribution feeders validates the approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Injection, confinement, and diagnosis of electrons and positrons in a permanent magnet dipole trap

Abstract Prerequisites for the goal of studying long-lived, magnetically confined, electron–positron pair plasmas in the laboratory include the injection of both species into the trap, long trapping times, and suitable diagnostic methods. Here we report recent progress on these tasks achieved in a simple dipole trap based on a supported permanent magnet. For the injection of electrons, both an $$\textbf{E}\times B$$ E × B drift technique (of a $$\sim $$ ∼ 2– $$\upmu $$ μ A, 6-eV beam) and “edge injection” (from a filament emitting a few mA and biased to some tens of volts) have been demonstrated; the former is suitable for low-density beams with smaller spatial and velocity spreads, while the latter employs fluctuations arising from collective behavior. To diagnose the edge-injected electrons, image potentials and currents induced on a wall probe, the magnet case, and wall electrodes were measured. Confinement of drift-injected positrons, measured experimentally, exhibited at least two well-separated timescales. Simulations reproduced this qualitatively, using a simple model of elastic collisions with residual background gas, and point to small adjustments for increasing trapping times. In a major upgrade to diagnostic capabilities, 25 bismuth germanate detectors, placed in three reentrant ports, are able to localize annihilation gammas, which will be used in future experiments to distinguish between different loss channels. Graphical abstract

von der Linden, J. (ORCID:0000000327477813)↗

Model reduction for a power grid model

We examine the complexity of constructing reduced order models for subsets of the variables needed to represent the state of the power grid. In particular, we apply model reduction techniques to the DeMarco-Zheng power grid model. We show that due to the oscillating nature of the solutions and the absence of timescale separation between resolved and unresolved variables, the construction of accurate reduced models becomes highly non-trivial because one has to account for long memory effects. In addition, we show that a reduced model that includes even a short memory is drastically better than a memoryless model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Dynamics of the Uranian Rings

Some of the problems of the shepherding satellite model of Goldreich ant tremaine are discussed. The following topics are studied: (1) optical depths of the all the observed narrow rings; (2) satellite and ring separation timescales; (3) ring edge sharpness; (4) shock formation in narrow rings; (5) the existence of small satellites near the Uranian rings; and (6) the apse and node alignments of the eccentric and inclined rings.

Dermott, S. F.↗

Multiscale Temporal Variability of the Global Air‐Sea CO 2 Flux Anomaly

Abstract The global air‐sea CO 2 flux (F) impacts and is impacted by a plethora of climate‐related processes operating at multiple time scales. In bulk mass transfer formulations, F is driven by physico‐ and bio‐chemical factors such as the air‐sea partial pressure difference (∆pCO 2 ), gas transfer velocity, sea surface temperature, and salinity–all varying at multiple time scales. To de‐convolve the impact of these factors on variability in F at different time scales, time‐resolved estimates of F were computed using a global data set assembled between 1988 and 2015. The F anomalies were defined as temporal deviations from the 28‐year time‐averaged value. Spectral analysis revealed four dominant timescales of variability in F–subseasonal, seasonal, interannual, and decadal with relative amplitude differences varying across regions. A second‐order Taylor series expansion was then conducted along these four timescales to separate drivers across differing regions. The analysis showed that on subseasonal timescales, wind speed variability explains some 66% of the global F anomaly and is the dominant driver. On seasonal, interannual, and decadal timescales, the ∆pCO 2 effect controlled by the ∆pCO 2 anomaly, explained much of the F anomaly. On decadal timescales, the F anomaly was almost entirely governed by the ∆pCO 2 effect with large contributions from high latitudes. The main drivers across timescales also dominate the regional F anomaly, particularly in the mid‐high latitude regions. Finally, the driver of the ∆pCO 2 effect was closely connected with the relative strength of atmospheric pCO 2 and the nonthermal component of oceanic pCO 2 anomaly associated with dissolved inorganic carbon and alkalinity.

Environmental Sciences & Ecology↗

Dissociation-induced depletion of high-energy reactant molecules as a mechanism for pressure-dependent rate constants for bimolecular reactions

In 1922, Lindemann proposed the now-well-known mechanism for pressure-dependent rate constants for unimolecular reactions: reactant molecules with sufficiently high energies dissociate more quickly than collisions can reestablish the Boltzmann distribution of the internal energies of the molecule during its dissociation at low pressures – yielding pressure-dependent rate constants for unimolecular reactions due to the preferential depletion of the high energy states capable of dissociation. In the last century, incredible progress has been made in achieving a far greater understanding of and quantitative predictions for unimolecular and association reactions. In the modern era, pressure-dependent phenomenological rate constants are now nearly universally used to describe the rates of unimolecular and associative reactions in phenomenological kinetic modeling. However, there is a second, more indirect, implication of Lindemann’s mechanism that relates to how these dissociation-induced non-equilibrium distributions impact bimolecular reactions, including non-associative bimolecular reactions – which are generally not considered to have pressure-dependent rate constants. Yet, as we show herein, the same high energy states depleted due to dissociation would otherwise react most rapidly in high-activation-energy bimolecular reactions – yielding a mechanism for pressure-dependent rate constants for bimolecular reactions (including non-associative reactions). Here, we present results from a case study for CH 2 O dissociation, isomerization, and bimolecular reaction with O2 to explore this question. Results from our master equation calculations indicate that the effect of dissociation-induced non-equilibrium distributions on bimolecular reactions can be substantial – even when chemical timescales are well separated from internal energy relaxational timescales (i.e. when the traditional rate constant description would be thought to apply). This effect is found to be more pronounced – and more complex – for bimolecular reactions involving molecular entities whose chemical timescales are merged with the internal energy relaxational timescales. Lastly, we present some ideas for discussion regarding what should be considered as “chemical species” in phenomenological kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Noise-induced stabilization of dynamical states with broken time-reversal symmetry

Under a high frequency drive, Josephson junctions demonstrate Shapiro steps of quantized voltage. These are dynamically stabilized states in which the phase across the junction locks to the external drive. We explore the stochastic switching between two symmetric steps at $\frac{hw}{2e}$ and –$\frac{hw}{2e}$. Surprisingly, the switching rate exhibits a pronounced nonmonotonicity as a function of temperature, violating the general expectation that transitions should become faster with temperature. As a result, we explain this behavior by realizing that the system retains memory of the dynamic state from which it is switching, thereby breaking the conventional simplifying assumptions about separations of timescales.

36 MATERIALS SCIENCE↗

Properties of Nonlinear Dynamo Waves

Dynamo theory offers the most promising explanation of the generation of the sun's magnetic cycle. Mean field electrodynamics has provided the platform for linear and nonlinear models of solar dynamos. However, the nonlinearities included are (necessarily) arbitrarily imposed in these models. This paper conducts a systematic survey of the role of nonlinearities in the dynamo process, by considering the behaviour of dynamo waves in the nonlinear regime. It is demonstrated that only by considering realistic nonlinearities that are non-local in space and time can modulation of the basic dynamo wave he achieved. Moreover, this modulation is greatest when there is a large separation of timescales provided by including a low magnetic Prandtl number in the equation for the velocity perturbations.

Tobias, S. M.↗

OGLE-2016-BLG-1227L: A Wide-separation Planet from a Very Short-timescale Microlensing Event

We present the analysis of the microlensing event OGLE-2016-BLG-1227. The light curve of this short-duration event appears to be a single-lens event affected by severe finite-source effects. Analysis of the light curve based on a single-lens single-source (1L1S) model yields very small values of the event timescale, t{sub E}∼3.5 days, and the angular Einstein radius, θ{sub E}∼0.009 mas, making the lens a candidate of a free-floating planet. Close inspection reveals that the 1L1S solution leaves small residuals with an amplitude of ΔI ≲ 0.03 mag. We find that the residuals are explained by the existence of an additional widely separated heavier lens component, indicating that the lens is a wide-separation planetary system rather than a free-floating planet. From Bayesian analysis, it is estimated that the planet has a mass of M{sub p}=0.79{sub −0.39}{sup +1.30} M{sub J} and it is orbiting a low-mass host star with a mass of M{sub host}=0.10{sub −0.05}{sup +0.17} M{sub ⊙} located with a projected separation of a{sub ⊥}=3.4{sub −1.0}{sup +2.1} au. The planetary system is located in the Galactic bulge with a line-of-sight separation from the source star of D{sub LS}=1.21{sub −0.63}{sup +0.96} kpc. The event shows that there is a range of deviations in the signatures of host stars for apparently isolated planetary lensing events and that it is possible to identify a host even when a deviation is subtle.

79 ASTRONOMY AND ASTROPHYSICS↗

Ultrafast Triplet Pair Separation and Triplet Trapping following Singlet Fission in Amorphous Pentacene Films

In this work, ultrafast infrared and electronic spectroscopy are used to examine the dynamics of triplet pair separation following singlet fission in amorphous and crystalline films of the model singlet fission chromophore, 6,13-bis(triisopropylsilylethynyl)pentacene (TIPS-Pn). Direct probing of correlated triplet pair intermediates through their unique vibrational frequencies, infrared electronic transitions, and indirectly through their visible triplet absorptions reveals that triplet pair separation occurs on similar picosecond timescales in both amorphous and crystalline films despite their markedly different average intermolecular coupling strengths. Although triplet pair separation occurs on similar timescales in both environments, measurements of diffusion-controlled triplet-triplet annihilation reveal that the diffusivity of triplet excited states is an order of magnitude lower in amorphous films. The data reveal the presence of sparse triplet traps in the amorphous environment that inhibit the transport of triplet excitons in comparison to crystalline films. These observations inform recent efforts to develop disordered and polymeric singlet fission sensitizers that contain amorphous regions. In particular, the data suggest that it may be possible to nanostructure amorphous or polymeric singlet fission sensitizers to allow ultrafast triplet pair separation and harvesting in photovoltaic and light emitting applications despite their low triplet exciton diffusivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Timescales of Soil Moisture Anomalies: Results from Two GCMs

Soil moisture anomalies dissipate over timescales that may span weeks to months. Characterizing the geographical and seasonal variations in these timescales can have important practical benefit; significant soil moisture "memory" allows long-lead forecasts of soil moisture, which have been found in recent studies to be essential for useful Ion--lead forecasts of precipitation in many regions. In this talk, we will present and compare the soil moisture timescales derived in two separate general circulation model (GCM) studies. Both studies employ multiple ensembles of short-term climate simulations. Timescales at a given point are effectively estimated by determining how quickly the soil moisture distribution generated in one ensemble of simulations (characterized by a unique set of initial soil moisture conditions) approaches that produced by another ensemble (characterized by a different set of initial soil moisture conditions). The talk will include a discussion of why the timescales produced by the two GCMs differ in some regions, and it will also describe the impact of soil moisture memory on simulated precipitation.

Koster, Randal D.↗

Global gyrokinetic simulations of ASDEX Upgrade up to the transport timescale with GENE–Tango

An accurate description of turbulence up to the transport timescale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport timescales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collisions, realistic geometries, toroidal rotation and electromagnetic effects for the first time. This approach overcomes gyrokinetic codes' limitations and enables high-fidelity profile calculations in experimentally relevant plasma conditions, significantly reducing the computational cost. We present numerical results of GENE–Tango for two ASDEX Upgrade discharges, one of which exhibits a pronounced peaking of the ion temperature profile not reproduced by TGLF–ASTRA. We show that GENE–Tango can correctly capture the ion temperature peaking observed in the experiment. By retaining different physical effects in the GENE simulations, e.g., collisions, toroidal rotation and electromagnetic effects, we show that the ion temperature profile's peaking can be linked to electromagnetic effects of submarginal (stable) KBM modes. Furthermore, based on these results, the expected GENE–Tango speedup for the ITER standard scenario is larger than two orders of magnitude compared to a single gyrokinetic simulation up to the transport timescale, possibly making first-principles ITER simulations feasible on current computing resources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗