Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Non linear dynamics”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Geometric phase in coupled cluster theory

It has been well-established that the topography around conical intersections between excited electronic states is incorrectly described by coupled cluster and many other single reference theories (the intersections are “defective”). Despite this, we show both analytically and numerically that the geometric phase effect (GPE) is correctly reproduced upon traversing a path around a defective excited-state conical intersection (CI) in coupled cluster theory. The theoretical analysis is carried out by using a non-Hermitian generalization of the linear vibronic coupling approach. Interestingly, the approach qualitatively explains the characteristic (incorrect) shape of the defective CIs and CI seams. Moreover, the validity of the approach and the presence of the GPE indicate that defective CIs are local (and not global) artifacts. This implies that a sufficiently accurate coupled cluster method could predict nuclear dynamics, including geometric phase effects, as long as the nuclear wavepacket never gets too close to the conical intersections.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Using extended MHD to explore lasers as a trigger for x-pinches

X-pinches have been shown to be a source of extremely intense x-ray emissions useful for diagnosing plasma dynamics and imaging biological objects. The most striking feature of an x-pinch is the hotspot, the point source from where all the x rays come from. Unfortunately, the exact timing and location of the hotspot are still unpredictable. Since an x-pinch hotspot forms from instabilities (like an m = 0 mode), we will computationally explore whether we can use a high-power laser to control the timing and location of these instabilities. Our goal is to reduce the temporal and spatial jitter associated with the x-ray burst. Using an extended MHD (XMHD) code, we explore the non-relativistic instability generation using a current profile of a 250 kA linear transformer driver system and laser characteristics of SLAC's Matter in Extreme Conditions laboratory. Our results include both laser-penetration results using a boundary-defined electromagnetic (EM) wave and instability results from a power-deposition method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Kinetics and mechanism of light-induced phase separation in a mixed-halide perovskite

Halide ion phase separation is a barrier to the application of mixed-halide perovskites, whereby the presence of large populations of photo-generated or injected carriers causes undesirable changes in the local band gap. Here, we investigate the mechanism of phase separation in CsPbI x Br 3-x perovskite single crystals driven by light. The phase separation process and its dynamics are visualized at the nanometer scale at cryogenic temperatures using in-situ scanning transmission electron microscopy and cathodoluminescence. Combined with phase field modeling that accounts for the coupling between electronic and chemical driving forces, our observations point to a spinodal decomposition mechanism in which both the amplitude of composition fluctuation and the characteristic length scale grow non-linearly with time. Our findings provide microscopic insights that can assist in further engineering mixed-halide perovskites either for stability or for intentional programming of the local halide ion composition, opening pathways to a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Aging phenomena in the two-dimensional complex Ginzburg-Landau equation

The complex Ginzburg-Landau equation with additive noise is a stochastic partial differential equation that describes a remarkably wide range of physical systems which include coupled non-linear oscillators subject to external noise near a Hopf bifurcation instability and spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations or oscillatory chemical reactions. Here, we employ a finite-difference method to numerically solve the noisy complex Ginzburg-Landau equation on a two-dimensional domain with the goal to investigate its non-equilibrium dynamics when the system is quenched into the “defocusing spiral quadrant”. We observe slow coarsening dynamics as oppositely charged topological defects annihilate each other, and characterize the ensuing aging scaling behavior. We conclude that the physical aging features in this system are governed by non-universal aging scaling exponents. We also investigate systems with control parameters residing in the “focusing quadrant”, and identify slow aging kinetics in that regime as well. Finally, we provide heuristic criteria for the existence of slow coarsening dynamics and physical aging behavior in the complex Ginzburg-Landau equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Low-temperature spin dynamics in the TmFeO 3 orthoferrite with a non-Kramers ion

We investigate the magnetic dynamics of the orthorhombic perovskite TmFeO 3 at low temperatures, below the spin reorientation transition at TSR≈80 K, by means of time-of-flight neutron spectroscopy. We find that the magnetic excitation spectrum combines two emergent collective modes associated with different magnetic sublattices. The Fe subsystem orders below TN~632 K into a canted antiferromagnetic structure and exhibits sharp, high-energy magnon excitations. We describe them using linear spin-wave theory, and reveal a pronounced anisotropy between in- and out-of-plane exchange interactions, which was mainly neglected in previous reports on the spin dynamics in orthoferrites. At lower energies, we find two crystalline electrical field (CEF) excitations of Tm3+ ions at energies of ~2 and 5 meV. In contrast to the sister compound YbFeO 3 , where the Yb 3+ ions form quasi-one-dimensional chains along the c axis, the Tm excitations show dispersion along both directions in the (0KL) scattering plane. Analysis of the neutron scattering polarization factor reveals a longitudinal polarization of the 2 meV excitation. To evaluate the effect of the CEF on the Tm 3+ ions, we perform point-charge model calculations, and their results quantitatively capture the main features of Tm single-ion physics, such as energies, intensities, and polarization of the CEF transitions, and the type of magnetic anisotropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Structural basis of RNA conformational switching in the transcriptional regulator 7SK RNP

7SK non-coding RNA (7SK) negatively regulates RNA polymerase II (RNA Pol II) elongation by inhibiting positive transcription elongation factor b (P-TEFb), and its ribonucleoprotein complex (RNP) is hijacked by HIV-1 for viral transcription and replication. Methylphosphate capping enzyme (MePCE) and La-related protein 7 (Larp7) constitutively associate with 7SK to form a core RNP, while P-TEFb and other proteins dynamically assemble to form different complexes. Here, we present the cryo-EM structures of 7SK core RNP formed with two 7SK conformations, circular and linear, and uncover a common RNA-dependent MePCE-Larp7 complex. Further, together with NMR, biochemical, and cellular data, these structures reveal the mechanism of MePCE catalytic inactivation in the core RNP, unexpected interactions between Larp7 and RNA that facilitate a role as an RNP chaperone, and that MePCE-7SK-Larp7 core RNP serves as a scaffold for switching between different 7SK conformations essential for RNP assembly and regulation of P-TEFb sequestration and release.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Scaling of ionic conductance in a fluctuating single-layer nanoporous membrane

Abstract Single-layer membranes have emerged as promising candidates for applications requiring high transport rates due to their low resistance to molecular transport. Owing to their atomically thin structure, these membranes experience significant microscopic fluctuations, emphasizing the need to explore their impact on ion transport processes. In this study, we investigate the effects of membrane fluctuations on the elementary scaling behavior of ion conductance $$G$$ G as a function of ion concentration $$c_{0}$$ c 0 , represented as $$G = \beta c_{0}^{\alpha }$$ G = β c 0 α , using molecular dynamics simulations. Our findings reveal that membrane fluctuations not only alter the conductance coefficient $$\beta$$ β but also the power-law exponent $$\alpha$$ α . We identify two distinct frequency regimes of membrane fluctuations, GHz-scale and THz-scale fluctuations, and examine their roles in conductance scaling. Furthermore, we demonstrate that the alteration of conductance scaling arises from the non-linearity between ion conductance and membrane shape. This work provides a fundamental understanding of ion transport in fluctuating membranes.

42 ENGINEERING↗
DOE OSTI2021

Stability and pulsation of the first dark stars

ABSTRACT The first bright objects to form in the Universe might not have been ‘ordinary’ fusion-powered stars, but ‘dark stars’ (DSs) powered by the annihilation of dark matter (DM) in the form of weakly interacting massive particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of the order of M⊙ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, \mathrm{ M}_\odot$, fueled by WIMPS weighing 100 GeV. We improve the previous implementation of the DM energy source into the stellar evolution code mesa. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \!1000 \, \mathrm{ M}_{\odot }$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code gyre. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, \mathrm{ M}_\odot$ are excited by the κ − γ and γ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Spin dynamics and exchange interaction in orthoferrite TbFeO 3 with non-Kramers rare-earth ion

The low-temperature spin dynamics of the orthorhombic TbFeO 3 perovskite has been studied. It has been found that the inelastic neutron scattering (INS) spectrum of the investigated compound contains two modes corresponding to different sublattices. It is shown that the iron subsystem orders antiferromagnetically at T N =632 K and exhibits the high-energy magnon dispersion. The magnetic dynamics of this subsystem has been described within the linear spin wave theory and the in-plane and out-of-plane exchange anisotropy has been demonstrated. The approach proposed previously to describe the magnon dispersion in the TmFeO 3 compound has been used. Three levels of the nondispersive crystal electric field corresponding to Tb 3+ ions have been found in the energy region below 40 meV at about 17, 26, and 35 meV. The behavior of the magnetic correlation length of the terbium subsystem has been determined by studying the diffuse scattering at different temperatures. The evolution of this subsystem has been numerically described within the point charge model. In conclusion, it is shown that the numerical data agree satisfactorily with the experiment and with the general concept of the single-ion approximation applied to the rare-earth subsystem of orthorhombic perovskites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Planck 2018 results

We report on the implications for cosmic inflation of the 2018 release of the Planck cosmic microwave background (CMB) anisotropy measurements. The results are fully consistent with those reported using the data from the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be n s = 0.9649 ± 0.0042 at 68% CL. We find no evidence for a scale dependence of n s , either as a running or as a running of the running. The Universe is found to be consistent with spatial flatness with a precision of 0.4% at 95% CL by combining Planck with a compilation of baryon acoustic oscillation data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, r0.002 < 0.10, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain r 0.002 < 0.056. In the framework of standard single-field inflationary models with Einstein gravity, these results imply that: (a) the predictions of slow-roll models with a concave potential, V"(Φ) < 0, are increasingly favoured by the data; and (b) based on two different methods for reconstructing the inflaton potential, we find no evidence for dynamics beyond slow roll. Three different methods for the non-parametric reconstruction of the primordial power spectrum consistently confirm a pure power law in the range of comoving scales 0.005 Mpc -1 ≲ k ≲ 0.2 Mpc -1 . A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectra. For the case of oscillatory features that are logarithmic or linear in k, this result is further strengthened by a new combined analysis including the Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions for the cosmological fluctuations. In correlated, mixed adiabatic and isocurvature models, the non-adiabatic contribution to the observed CMB temperature variance is constrained to 1.3%, 1.7%, and 1.7% at 95% CL for cold dark matter, neutrino density, and neutrino velocity, respectively. Planck power spectra plus lensing set constraints on the amplitude of compensated cold dark matter-baryon isocurvature perturbations that are consistent with current complementary measurements. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadupolar modulation of the primordial fluctuations. However, the polarization data do not support physical models for a scale-dependent dipolar modulation. All these findings support the key predictions of the standard single-field inflationary models, which will be further tested by future cosmological observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Time-dependent saturation and physics-based nonlinear model of cross-beam energy transfer

The nonlinear physics of cross-beam energy transfer (CBET) for multi-speckled laser beams is examined using large-scale particle-in-cell simulations for a range of laser and plasma conditions relevant to indirect-drive inertial confinement fusion (ICF) experiments. The time-dependent growth and saturation of CBET involve complex, nonlinear ion and electron dynamics, including ion trapping-induced enhancement and detuning, ion acoustic wave (IAW) nonlinearity, oblique forward stimulated Raman scattering (FSRS), and backward stimulated Brillouin scattering (BSBS) in a CBET-amplified seed beam. Ion-trapping-induced detuning of CBET is captured in the kinetic linear response by a new δf-Gaussian-mixture algorithm, enabling an accurate characterization of trapping-induced non-Maxwellian distributions. Ion trapping induces nonlinear processes, such as changes to the IAW dispersion and nonlinearities (e.g., bowing and self-focusing), which, together with pump depletion, FSRS, and BSBS, determine the time-dependent nature and level of CBET gain as the system approaches a steady state. Using VPIC simulations at intensities at and above the onset threshold for ion trapping and the insight from the time-dependent saturation analyses, we construct a nonlinear CBET model from local laser and plasma conditions that predicts the CBET gain and the energy deposition into the plasma. This model is intended to provide a more accurate, physics-based description of CBET saturation over a wide range of conditions encountered in ICF hohlraums compared with linear CBET gain models with ad hoc saturation clamps often used in laser ray-based methods in multi-physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of the Grazing Incidence Geometry on X-ray Photon Correlation Spectroscopy Measurements

X-ray photon correlation spectroscopy (XPCS) is a versatile tool to measure dynamics on the nanometer to micrometer scale in bulk samples. XPCS has also been applied in grazing incidence (GI) geometry to examine the dynamics of surface layers. However, considering GI scattering experiments more universally, the GI geometry leads to a superposition of signals due to reflection and refraction effects, also known from the distorted-wave Born approximation (DWBA). In this paper, the impact of these reflection and refraction effects on the correlation analysis is determined experimentally by measuring grazing incidence transmission XPCS (GT-XPCS) and grazing incidence XPCS (GI-XPCS) simultaneously for a thin film sample, showing non-equilibrium dynamics. Here the results of the GI and GT geometry comparisons are combined within the framework of the standardly applied, simplified DWBA. These calculations allow identifying the main contributions of the detected signal from the leading scattering terms along the out-of-plane direction q z , which dominate the measured intensity pattern on the detector. In combination with the calculation of the non-linear effect of refraction in GTSAXS and GISAXS, it is possible to identify experimental conditions that can be chosen to run experiments and data analysis as close as possible to transmission XPCS and to explain which limitations for data interpretations are observed. Consequently, the beam exposure can be significantly reduced by using GI geometry only. Calculations of experimental settings prior to experiments are detailed to determine suitable q z regions for a variety of material systems measured in bulk-sensitive GI-XPCS experiments, allowing us to determine the scaling behavior of typical decay times as a function of q that is comparable to the scaling behavior obtained in distortion-free GT-XPCS or transmission XPCS experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗
DOE OSTI2020

Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics

Quantum computers can be used to simulate nonlinear non-Hamiltonian classical dynamics on phase space by using the generalized Koopman–von Neumann formulation of classical mechanics. The Koopman–von Neumann formulation implies that the conservation of the probability distribution function on phase space, as expressed by the Liouville equation, can be recast as an equivalent Schrödinger equation on Hilbert space with a Hermitian Hamiltonian operator and a unitary propagator. This Schrödinger equation is linear in the momenta because it derives from a constrained Hamiltonian system with twice the classical phase-space dimension. A quantum computer with finite resources can be used to simulate a finite-dimensional approximation of this unitary evolution operator. Quantum simulation of classical dynamics is exponentially more efficient than a deterministic Eulerian discretization of the Liouville equation if the Koopman–von Neumann Hamiltonian is sparse. Utilizing quantum walk techniques for state preparation and amplitude estimation for the calculation of observables leads to a quadratic improvement over classical probabilistic Monte Carlo algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Coherently coupled quantum oscillators for quantum reservoir computing

Here, we analyze the properties of a quantum system composed of two coherently coupled quantum oscillators and show through simulations that it fulfills the two properties required for reservoir computing: non-linearity and fading memory. We first show that the basis states of this system apply a set of nonlinear transformations on the input signals and thus can implement neurons. We then show that the system exhibits a fading memory that can be controlled by its dissipation rates. Finally we show that a strong coupling between the oscillators is important in order to ensure complex dynamics and to populate a number of basis state neurons that is exponential in the number of physical devices.

Dudas, Julien↗
DOE OSTI2025

Scrambling Dynamics with Imperfections in a Solvable Model

We study how probes of quantum scrambling dynamics respond to two kinds of imperfections -- unequal forward and backward evolutions and decoherence -- in a solvable Brownian circuit model. We calculate a ``renormalized'' out-of-time-order correlator (ROTOC) in the model with $N$ qubits, and we show that the circuit-averaged ROTOC is controlled by an effective probability distribution in operator weight space which obeys a system of $N$ non-linear equations of motion. These equations can be easily solved numerically for large system sizes which are beyond the reach of exact methods. Moreover, for an operator initially concentrated on weight one $w_0=1$, we provide an exact solution to the equations in the thermodynamic limit of many qubits that is valid for all times, all non-vanishing perturbation strengths $p\gtrsim 1/\sqrt{N}$, and all decoherence strengths. We also show that a generic initial condition $w_0 >1$ leads to a metastable state that eventually collapses to the $w_0=1$ case after a lifetime $\sim \log(N/w_0)$. Our results highlight situations where it is still possible to extract the unperturbed chaos exponent even in the presence of imperfections, and we comment on the applications of our results to existing experiments with nuclear spins and to future scrambling experiments.

FOS: Physical sciences↗
DOE OSTI2022

Addressing Uncertainty in MultiSector Dynamics Research

MultiSector Dynamics (MSD) research explores the dynamics and co-evolutionary pathways of human and Earth systems and the emerging interdependent sectors of energy, water, agriculture, and transportation among others. The interactions between these sectors are central to MSD science, as they capture how processes and feedbacks across Earth, environmental, infrastructure, and societal systems shape transitions, socioeconomic risks, and the provision of services. By definition, MSD research requires deep integration across diverse scientific disciplines, ranging from the natural to the social sciences and engineering. All these disciplines apply a variety of numerical simulation models to study and understand their underlying systems of focus. The utility of these models hinges on the fidelity with which they represent the real systems and their ability to produce novel insights about systems and their interactions. The coupled human-natural systems typically represented are shaped by a multitude of interdependent human and natural processes which, when modeled, translate to highly complex, non-linear, interacting behaviors. This ever-increasing complexity massively expands the uncertainty space of a model and can be found in model inputs, processes and parameters. This is further amplified when several models are combined to answer multisectoral questions, as additional uncertainty regarding coupling relationships and interactions is introduced.

IM3, Uncertainty and Sensitivity Analysis, multise↗