Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time Evolution”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Efficient quantum algorithm for dissipative nonlinear differential equations

Significance Nonlinear differential equations appear in many domains and are notoriously difficult to solve. Whereas previous quantum algorithms for general nonlinear differential equations have complexity exponential in the evolution time, we give the first quantum algorithm for dissipative nonlinear differential equations that is efficient provided the dissipation is sufficiently strong relative to nonlinear and forcing terms and the solution does not decay too rapidly. We also establish a lower bound showing that differential equations with sufficiently weak dissipation have worst-case complexity exponential in time, giving an almost tight classification of the quantum complexity of simulating nonlinear dynamics. Furthermore, numerical results for the Burgers equation suggest that our algorithm may potentially address complex nonlinear phenomena even in regimes with weaker dissipation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Circuit complexity and functionality: A statistical thermodynamics perspective

Circuit complexity, defined as the minimum circuit size required for implementing a particular Boolean computation, is a foundational concept in computer science. Determining circuit complexity is believed to be a hard computational problem. Recently, in the context of black holes, circuit complexity has been promoted to a physical property, wherein the growth of complexity is reflected in the time evolution of the Einstein-Rosen bridge (“wormhole”) connecting the two sides of an anti-de Sitter “eternal” black hole. Here, we are motivated by an independent set of considerations and explore links between complexity and thermodynamics for functionally equivalent circuits, making the physics-inspired approach relevant to real computational problems, for which functionality is the key element of interest. In particular, our thermodynamic framework provides an alternative perspective on the obfuscation of programs of arbitrary length—an important problem in cryptography—as thermalization through recursive mixing of neighboring sections of a circuit, which can be viewed as the mixing of two containers with “gases of gates.” This recursive process equilibrates the average complexity and leads to the saturation of the circuit entropy, while preserving functionality of the overall circuit. The thermodynamic arguments hinge on ergodicity in the space of circuits which we conjecture is limited to disconnected ergodic sectors due to fragmentation. The notion of fragmentation has important implications for the problem of circuit obfuscation as it implies that there are circuits of same size and functionality that cannot be connected via a polynomial number of local moves. Furthermore, we argue that fragmentation is unavoidable unless the complexity classes NP and coNP coincide, a statement that implies the collapse of the polynomial hierarchy of computational complexity theory to its first level.

Science & Technology - Other Topics↗
DOE OSTI2020

Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation

Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Bayesian inference of electron density and ion temperature profiles from neutral beam and halo Balmer- α emission at Wendelstein 7-X

By employing Bayesian inference techniques, the full electron density profile from the plasma core to the edge of Wendelstein 7-X (W7-X) is inferred solely from neutral hydrogen beam and halo Balmer- α (H α ) emission data. The halo is a cloud of neutrals forming in the vicinity of the injected neutral beam due to multiple charge exchange reactions. W7-X is equipped with several neutral hydrogen beam heating sources and an H α spectroscopy system that views these sources from different angles and penetration depths in the plasma. As the beam and halo emission form complex spectra for each spatial point that are non-linearly dependent on the plasma density profile and other parameters, a complete model from the neutral beam injection and halo formation through to the spectroscopic measurements is required. The model is used here to infer electron density profiles for a range of common W7-X plasma scenarios. The inferred profiles show good agreement with profiles determined by the Thomson scattering and interferometry diagnostics across a broad range of absolute densities without any changes to the input or fitting parameters. The time evolution of the density profile in a discharge with continuous core density peaking is successfully reconstructed, demonstrating sufficient spatial resolution to infer strongly shaped profiles. Furthermore, it is shown as a proof of concept that the model is also able to infer the main ion temperature profile using the same data set.

forward modeling↗
DOE OSTI2024

Hamiltonian switching control of noisy bipartite qubit systems

Abstract We develop a Hamiltonian switching ansatz for bipartite control that is inspired by the quantum approximate optimization algorithm, to mitigate environmental noise on qubits. We demonstrate the control for a central spin coupled to bath spins via isotropic Heisenberg interactions, and then make physical applications to the protection of quantum gates performed on superconducting transmon qubits coupling to environmental two-level-systems (TLSs) through dipole-dipole interactions, as well as on such qubits coupled to both TLSs and a Lindblad bath. The control field is classical and acts only on the system qubits. We use reinforcement learning with policy gradient to optimize the Hamiltonian switching control protocols, using a fidelity objective for specific target quantum gates. We use this approach to demonstrate effective suppression of both coherent and dissipative noise, with numerical studies achieving target gate implementations with fidelities over 0.9999 (four nines) in the majority of our test cases and showing improvement beyond this to values of 0.999 999 999 (nine nines) upon a subsequent optimization by GRadient Ascent Pulse Engineering (GRAPE). We analyze how the control depth, total evolution time, number of environmental TLS, and choice of optimization method affect the fidelity achieved by the optimal protocols and reveal some critical behaviors of bipartite control of quantum gates.

Physics↗
DOE OSTI2022

Is Jet Re-orientation the Elusive Trigger for Star Formation Suppression in Radio Galaxies?

Abstract Jet re-orientation associated with the time evolution of radio quasars explains the formation of X-shaped radio galaxies and their preference for isolated environments. But since X-shaped radio galaxies are generally not found in dense environments (e.g., groups/clusters), the jet re-orientation phenomenon for radio galaxies in groups and clusters has been ignored. We take a closer look at the re-orientation of FRI jets with respect to FRII jets, and find that it may constitute the as-yet unidentified trigger for star formation suppression in radio galaxies. We show how the recently explored radio “red geyser” galaxies can be interpreted in this context and ultimately reveal a deeper understanding of why FRII radio galaxies are on one side of the star formation enhancement/suppression divide compared to FRI radio galaxies.

Astronomy & Astrophysics↗
DOE OSTI2025

Stability evaluation and mitigation strategies in advanced tokamaks using 3D MHD spectroscopy

Multi-modal, active 3D MHD spectroscopy is applied in high-performance advanced tokamak scenarios to study their stability time evolution, revealing an intriguing dependence on both and . A tailored applied 3D field provides a 3D plasma response to extract the growth rate of the least stable mode. The estimated growth rate finds a decrease in stability when the minimum in the safety factor (q) passes through 2.0 and reveals inherent risks of crossing an additional rational surface at integer q min , even above the usual q = 1 sawtooth condition. Based on this result, the potential scenario in which q min ~ 2 can be safely crossed during a more stable lower β N phase was investigated, and the improved stability of this scenario is confirmed by the estimated growth rate. This shows that 3D MHD spectroscopy can offer insights into strategies for improving stability by identifying the vulnerable aspects of such scenarios. In addition, the method highlights its potential for instability avoidance by enabling early detection of multiple modes, even before magnetic coils can measure them. The measured growth rate by the 3D MHD spectroscopy shows its reliability by exhibiting a correlation with the programmed rises in plasma beta across various high β N and high q min discharges. In addition, this method is successfully applied during rapidly evolving I p ramp-up phases, a key part of the scenario development. By achieving reasonable growth rate measurements at high-performance scenario developments, this technique contributes to the development of advanced diagnostic tools for tokamak scenario stability, which will help identify an effective pathway to stable, high-performance scenarios.

Yang, S. M. [Princeton Plasma Physics Laboratory (↗
DOE OSTI2022

Statistical properties of sites visited by independent random walks

The set of visited sites and the number of visited sites are two basic properties of the random walk trajectory. Here, we consider two independent random walks on hyper-cubic lattices and study ordering probabilities associated with these characteristics. The first is the probability that during the time interval (0, t), the number of sites visited by a walker never exceeds that of another walker. The second is the probability that the sites visited by a walker remain a subset of the sites visited by another walker. Using numerical simulations, we investigate the leading asymptotic behaviors of the ordering probabilities in spatial dimensions d = 1, 2, 3, 4. We also study the time evolution of the number of ties between the number of visited sites. We show analytically that the average number of ties increases as a 1 ln t with a 1 = 0.970 508 in one dimension and as (ln t) 2 in two dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantum Davidson algorithm for excited states

Abstract Excited state properties play a pivotal role in various chemical and physical phenomena, such as charge separation and light emission. However, the primary focus of most existing quantum algorithms has been the ground state, as seen in quantum phase estimation and the variational quantum eigensolver (VQE). Although VQE-type methods have been extended to explore excited states, these methods grapple with optimization challenges. In contrast, the quantum Krylov subspace (QKS) method has been introduced to address both ground and excited states, positioning itself as a cost-effective alternative to quantum phase estimation. However, conventional QKS methodologies depend on a pre-generated subspace through real or imaginary-time evolutions. This subspace is inherently expansive and can be plagued with issues like slow convergence or numerical instabilities, often leading to relatively deep circuits. Our research presents an economic QKS algorithm, which we term the quantum Davidson (QDavidson) algorithm. This innovation hinges on the iterative expansion of the Krylov subspace and the incorporation of a pre-conditioner within the Davidson framework. By using the residues of eigenstates to expand the Krylov subspace, we manage to formulate a compact subspace that aligns closely with the exact solutions. This iterative subspace expansion paves the way for a more rapid convergence in comparison to other QKS techniques, such as the quantum Lanczos. Using quantum simulators, we employ the novel QDavidson algorithm to delve into the excited state properties of various systems, spanning from the Heisenberg spin model to real molecules. Compared to the existing QKS methods, the QDavidson algorithm not only converges swiftly but also demands a significantly shallower circuit. This efficiency establishes the QDavidson method as a pragmatic tool for elucidating both ground and excited state properties on quantum computing platforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Quantum simulation of boson-related Hamiltonians: techniques, effective Hamiltonian construction, and error analysis

Elementary quantum mechanics proposes that a closed physical system consistently evolves in a reversible manner. However, control and readout necessitate the coupling of the quantum system to the external environment, subjecting it to relaxation and decoherence. Consequently, system-environment interactions are indispensable for simulating physically significant theories. A broad spectrum of physical systems in condensed-matter and high-energy physics, vibrational spectroscopy, and circuit and cavity QED necessitates the incorporation of bosonic degrees of freedom, such as phonons, photons, and gluons, into optimized fermion algorithms for near-future quantum simulations. In particular, when a quantum system is surrounded by an external environment, its basic physics can usually be simplified to a spin or fermionic system interacting with bosonic modes. Nevertheless, troublesome factors such as the magnitude of the bosonic degrees of freedom typically complicate the direct quantum simulation of these interacting models, necessitating the consideration of a comprehensive plan. This strategy should specifically include a suitable fermion/boson-to-qubit mapping scheme to encode sufficiently large yet manageable bosonic modes, and a method for truncating and/or downfolding the Hamiltonian to the defined subspace for performing an approximate but highly accurate simulation, guided by rigorous error analysis. In this pedagogical tutorial review, we aim to provide such an exhaustive strategy, focusing on encoding and simulating certain bosonic-related model Hamiltonians, inclusive of their static properties and time evolutions. Specifically, we emphasize two aspects: (1) the discussion of recently developed quantum algorithms for these interacting models and the construction of effective Hamiltonians, and (2) a detailed analysis regarding a tightened error bound for truncating the bosonic modes for a class of fermion-boson interacting Hamiltonians.

bosonic Hamiltonian↗
DOE OSTI2020

Comparison of stellar populations in simulated and real post-starburst galaxies in MaNGA

ABSTRACT Recent integral field spectroscopic (IFS) surveys have revealed radial gradients in the optical spectral indices of post-starburst (PSB) galaxies, which can be used to constrain their formation histories. We study the spectral indices of post-processed mock IFS datacubes of binary merger simulations, carefully matched to the properties of the MaNGA IFS survey, with a variety of black hole (BH) feedback models, progenitor galaxies, orbits, and mass ratios. Based on our simulation sample, we find that only major mergers on prograde–prograde or retrograde–prograde orbits in combination with a mechanical BH feedback model can form galaxies with weak enough ongoing star formation, and therefore absent H α emission, to be selected by traditional PSB selection methods. We find strong fluctuations in nebular emission line strengths, even within the PSB phase, suggesting that H α selected PSBs are only a subsample of the underlying population. The global PSB population can be more robustly identified using stellar continuum-based approaches. The difficulty in reproducing the very young PSBs in simulations potentially indicates that new sub-resolution star formation recipes are required to properly model the process of star formation quenching. In our simulations, we find that the starburst peaks at the same time at all radii, but is stronger and more prolonged in the inner regions. This results in a strong time evolution in the radial gradients of the spectral indices that can be used to estimate the age of the starburst without reliance on detailed star formation histories from spectral synthesis models.

Zheng, Yirui↗
DOE OSTI2023

New constraints on cosmological modified gravity theories from anisotropic three-point correlation functions of BOSS DR12 galaxies

Here we report a new test of modified gravity theories using the large-scale structure of the Universe. This paper is the first attempt to (1) apply a joint analysis of the anisotropic components of galaxy two- and three-point correlation functions (2 and 3PCFs) to actual galaxy data and (2) constrain the non-linear effects of degenerate higher-order scalar-tensor (DHOST) theories on cosmological scales. Applying this analysis to the Baryon Oscillation Spectroscopic Survey (BOSS) data release 12, we obtain the lower bounds of –1.655 < ξ t and –0.504 < ξ s at the $95{{\ \rm per\ cent}}$ confidence level on the parameters characterizing the time evolution of the tidal and shift terms of the second-order velocity field. These constraints are consistent with GR predictions of ξ t = 15/1144 and ξs = 0. Moreover, they represent a 35-fold and 20-fold improvement, respectively, over the joint analysis with only the isotropic 3PCF. We ensure the validity of our results by investigating various quantities, including theoretical models of the 3PCF, window function corrections, cumulative S/N, Fisher matrices, and statistical scattering effects of mock simulation data. We also find statistically significant discrepancies between the BOSS data and the Patchy mocks for the 3PCF measurement. Finally, we package all of our 3PCF analysis codes under the name hitomi and make them publicly available so that readers can reproduce all the results of this paper and easily apply them to ongoing future galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Schoeffler, K. M. (ORCID:0000000276447381)↗
DOE OSTI2023

Measuring the streaming motion in the Milky Way disc with Gaia EDR3+

We map the 3D kinematics of the Galactic disc out to 3.5 kpc from the Sun, and within 0.75 kpc from the mid-plane of the Galaxy. To this end, we combine high-quality astrometry from Gaia EDR3, with heliocentric line-of-sight velocities from Gaia DR2, and spectroscopic surveys including APOGEE, GALAH, and LAMOST. We construct an axisymmetric model for the mean velocity field, and subtract this on a star-by-star basis to obtain the residual velocity field in the Galactocentric components (V Φ , V R , V z ), and V los . The velocity residuals are quantified using the power spectrum, and we find that the peak power (A/[km s –1 ]) in the mid-plane (|z| < 0.25 kpc) is (A Φ , A R , A Z , A los ) = (4.2,8.5,2.6,4.6), at 0.25 < |z|/[kpc] < 0.5, is (A Φ , A R , A Z , A los ) = (4.0,7.9,3.6,5.3), and at 0.5 < |z|/[kpc] < 0.75, is (A Φ , A R , A Z , A los ) = (1.9,6.9,5.2,6.4). Our results provide a sophisticated measurement of the streaming motion in the disc and in the individual components. We find that streaming is most significant in VR, and at all heights (|Z|) probed, but is also non-negligible in other components. Additionally, we find that patterns in velocity field overlap spatially with models for spiral arms in the Galaxy. Our simulations show that phase-mixing of disrupting spiral arms can generate such residuals in the velocity field, where the radial component is dominant, just as in real data. Furthermore, we also find that with time evolution, both the amplitude and physical scale of the residual motion decrease.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Beyond – ΛCDM constraints from the full shape clustering measurements from BOSS and eBOSS

We analyse the full shape of anisotropic clustering measurements from the extended Baryon Oscillation Spectroscopic Survey quasar sample together with the combined galaxy sample from the Baryon Oscillation Spectroscopic Survey. We obtain constraints on the cosmological parameters independent of the Hubble parameter h for the extensions of the Lambda cold dark matter (ΛCDM) models, focusing on cosmologies with free dark energy equation of state parameter w. We combine the clustering constraints with those from the latest cosmic microwave background data from Planck to obtain joint constraints for these cosmologies for w and the additional extension parameters – its time evolution w a , the physical curvature density ω K and the neutrino mass sum Σm ν . Our joint constraints are consistent with a flat ΛCDM cosmological model within 68 per cent confidence limits. We demonstrate that the Planck data are able to place tight constraints on the clustering amplitude today, σ 12 , in cosmologies with varying w and present the first constraints for the clustering amplitude for such cosmologies, which is found to be slightly higher than the ΛCDM value. Additionally, we show that when we vary w and allow for non-flat cosmologies and the physical curvature density is used, Planck prefers a curved universe at 4σ significance, which is ~2σ higher than when using the relative curvature density Ω K . Finally, when w is varied freely, clustering provides only a modest improvement (of 0.021 eV) on the upper limit of Σm ν .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Spatially quasi-periodic water waves of finite depth

We present a numerical study of spatially quasi-periodic gravity-capillary waves of finite depth in both the initial value problem and travelling wave settings. We adopt a quasi-periodic conformal mapping formulation of the Euler equations, where one-dimensional quasi-periodic functions are represented by periodic functions on a higher-dimensional torus. We compute the time evolution of free surface waves in the presence of a background flow and a quasi-periodic bottom boundary and observe the formation of quasi-periodic patterns on the free surface. Two types of quasi-periodic travelling waves are computed: small-amplitude waves bifurcating from the zero-amplitude solution and larger-amplitude waves bifurcating from finite-amplitude periodic travelling waves. We derive weakly nonlinear approximations of the first type and investigate the associated small-divisor problem. We find that waves of the second type exhibit striking nonlinear behaviour, e.g. the peaks and troughs are shifted non-periodically from the corresponding periodic waves due to the activation of quasi-periodic modes.

Science & Technology - Other Topics↗