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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Trapped-ion quantum simulation of collective neutrino oscillations

It is well known that the neutrino flavor in extreme astrophysical environments changes under the effect of three contributions: the vacuum oscillation, the interaction with the surrounding matter, and the collective oscillations due to interactions between different neutrinos. The latter adds a nonlinear contribution to the equations of motion, making the description of their dynamics complex. In this work we study various strategies to simulate the coherent collective oscillations of a system of N neutrinos in the two-flavor approximation using quantum computation. This was achieved by using a pair-neutrino decomposition designed to account for the fact that the flavor Hamiltonian, in the presence of the neutrino-neutrino term, presents an all-to-all interaction that makes the implementation of the evolution dependent on the qubit topology. We analyze the Trotter error caused by the decomposition demonstrating that the complexity of the implementation of time evolution scales polynomially with the number of neutrinos and that the noise from near-term quantum device simulation can be reduced by optimizing the quantum circuit decomposition and exploiting a full-qubit connectivity. We find that the gate complexity using second order Trotter-Suzuki formulas scales better with system size than with other decomposition methods such as quantum signal processing. In conclusion, we finally present the application and the results of our algorithm on a real quantum device based on trapped-ion qubits.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Parallel Time Integration for Constrained Optimization

The number of transistors in an average processor continues to increase, but individual clock speeds have plateaued. Those transistors are instead going into additional cores, increasing the number of different things that a processor can do at once and placing an emphasis on parallel computation. Many problems in scientific computing follow a time-evolution model, and it can be difficult to solve such problems in parallel across the temporal domain. The Multi-Grid Reduction In Time (MGRIT) algorithm, developed at Lawrence Livermore National Laboratory (LLNL), solves differential equations with a method designed specifically to take advantage of extreme numbers of processors by parallelizing across time. The Tri-diagonal MGRIT (TriMGRIT) algorithm, also developed at LLNL, is a generalization of MGRIT which enables parallel-in-time solving of a greater number of problems. Constrained optimization problems, in particular, may be solved in parallel using TriMGRIT. These consist of choosing a control function such that an objective functional is minimized, constrained by a differential-equation. We consider two such problems: applying torque to a pendulum to bring it to a gentle stop and moving a crowd of people from one distribution into another. We also perform some miscellaneous theoretical and practical research, including investigating the use of a line-search subroutine to refine intermediate TriMGRIT results and preliminary work on strategies for choosing operators for TriMGRIT to use.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exploiting Synergies Between Lensing and BAO Surveys for Improved Cosmological Constraints

Pinning down the nature of dark energy is one of the most pressing questions in modern physics and is one of the five key science drivers identified in the P5 report [P514]. Dark energy is thought to be either a cosmological constant with an equation of state w = P/ρ = −1 which remains constant at all times, a new type of fluid with an equation of state that varies with time (w ̸= constant), or dark energy might indicate a breakdown of Einstein’s Theory of General Relativity (GR). It is of critical importance to distinguish between these three scenarios. This can only be accomplished by ambitious and demanding measurements of both the expansion rate of the universe (to track the time evolution of dark energy) together with measurements of the rate at which cosmic structures, such as galaxies and clusters of galaxies, grow with time (the “growth rate”). Because of the challenging nature of such observations, and because no single probe simultaneously measures both expansion and growth, the Dark Energy Task Force [A+06] emphasized the importance of using multiple distinct methods to characterize dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evolution of Robustness in Growing Random Networks

Networks are widely used to model the interaction between individual dynamic systems. In many instances, the total number of units and interaction coupling are not fixed in time, and instead constantly evolve. In networks, this means that the number of nodes and edges both change over time. Various properties of coupled dynamic systems, such as their robustness against noise, essentially depend on the structure of the interaction network. Therefore, it is of considerable interest to predict how these properties are affected when the network grows as well as their relationship to the growth mechanism. Here, we focus on the time evolution of a network’s Kirchhoff index. We derive closed-form expressions for its variation in various scenarios, including the addition of both edges and nodes. For the latter case, we investigate the evolution where single nodes with one or two edges connecting to existing nodes are added recursively to a network. In both cases, we derive the relations between the properties of the nodes to which the new node connects along with the global evolution of network robustness. In particular, we show how different scalings of the Kirchhoff index can be obtained as a function of the number of nodes. We illustrate and confirm this theory via numerical simulations of randomly growing networks.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Evolution and dynamics of shear-layer structures in near-wall turbulence

Near-wall flow structures in turbulent shear flows are analyzed, with particular emphasis on the study of their space-time evolution and connection to turbulence production. The results are obtained from investigation of a database generated from direct numerical simulation of turbulent channel flow at a Reynolds number of 180 based on half-channel width and friction velocity. New light is shed on problems associated with conditional sampling techniques, together with methods to improve these techniques, for use both in physical and numerical experiments. The results clearly indicate that earlier conceptual models of the processes associated with near-wall turbulence production, based on flow visualization and probe measurements need to be modified. For instance, the development of asymmetry in the spanwise direction seems to be an important element in the evolution of near-wall structures in general, and for shear layers in particular. The inhibition of spanwise motion of the near-wall streaky pattern may be the primary reason for the ability of small longitudinal riblets to reduce turbulent skin friction below the value for a flat surface.

Johansson, Arne V.↗
NASA NTRS1986

An observational study of the final warming in the Southern Hemisphere stratosphere

The final warming in the Southern Hemisphere stratosphere during the period September 1-November 30 for the years 1978-1983 is studied. The change in zonally averaged temperatures during that period is almost the same for each year. However, there is substantial interannual variability in the time evolution of the warming as events of enhanced wave, mean-flow interaction associated with bursts of upward propagating planetary waves develop at different times in different years. Polar temperatures in the middle stratosphere can increase as much as about 20 K in 4 days during these events. Such dynamical processes can play a significant role in the interannual variability of stratospheric ozone and other constituents.

Farrara, John D.↗
DOE OSTI2022

A low-rank power iteration scheme for neutron transport criticality problems

Computing effective eigenvalues for neutron transport often requires a fine numerical resolution. Here, the main challenge of such computations is the high memory effort of classical solvers, which limits the accuracy of chosen discretizations. In this work, we derive a method for the computation of effective eigenvalues when the underlying solution has a low-rank structure. This is accomplished by utilizing dynamical low-rank approximation (DLRA), which is an efficient strategy to derive time evolution equations for low-rank solution representations. The main idea is to interpret the iterates of the classical inverse power iteration as pseudo-time steps and apply the DLRA concepts in this framework. In our numerical experiment, we demonstrate that our method significantly reduces memory requirements while achieving the desired accuracy. Analytic investigations show that the proposed iteration scheme inherits the convergence speed of the inverse power iteration, at least for a simplified setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Collective neutrino oscillations on a quantum computer with hybrid quantum-classical algorithm

We simulate the time evolution of collective neutrino oscillations in two-flavor settings on a quantum computer. We explore the generalization of Trotter-Suzuki approximation to time-dependent Hamiltonian dynamics. The trotterization steps are further optimized using the Cartan decomposition of two-qubit unitary gates U ϵ SU(4) in the minimum number of controlled-NOT (CNOT) gates making the algorithm more resilient to the hardware noise. As a result, a more efficient hybrid quantum-classical algorithm is also explored to solve the problem on noisy intermediate-scale quantum devices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1978

Motion of a curved vortex filament with decaying vortical core and axial velocity

The motion and decay of a curved vortex filament having large axial and circumferential velocity components in a three-dimensional stream are analyzed by using the method of matched asymptotic expansions of the incompressible Navier-Stokes equations. The small parameter is the square root of the ratio of the kinematic viscosity to the circulation. The outer region is analyzed by the classical Biot-Savart law, and its solution is matched to that of the inner region, where viscous effects are important. Equations describing the coupling between the inner vortex structure and the motion of the vortex filament as well as the time evolution of the inner vortex structure are obtained. Equations are derived for the motion of the vortex filament and for the change and decay in time and space of the leading-order circumferential and axial velocity and vorticity components. Solutions are constructed for these components in terms of initial data.

Callegari, A. J.↗
NASA NTRS1979

A thermal scanning study of coastal upwelling in Lake Superior

The use of a thermal scanner to monitor the time evolution of the thermal structure of the coastal waters in Lake Superior during an upwelling event is described. Mosaics of thermal imagery from ten different times are described. Qualitative descriptions of the imagery give insight into the upwelling event. Recommendations for future use of a thermal scanner to monitor an upwelling event are discussed.

Scarpace, F. L.↗
NASA NTRS1990

Second mode interactions in supersonic boundary layers

The nonlinear evolution of a two-dimensional second mode unstable wave in a Mach 4.5 wall-bounded flow is computed by solving the full time-dependent compressible Navier-Stokes equations. A highly accurate solution is obtained using spectral collocation methods. It is shown that departure from linearity first occurs in the critical layer due to the cubic nonlinearities in the momentum equation. This is a direct result of the large density perturbations in this regime. Time evolution studies of the growth rate as a function of normal distance from the plate suggests that the mode is evolving toward a nonlinear saturated state, and that this problem is possibly amenable to standard weakly nonlinear perturbation methods.

Erlebacher, Gordon↗
DOE OSTI2025

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗
NASA NTRS2004

Effect of gravity on the caloric stimulation of the inner ear

Robert Barany won the 1914 Nobel Prize in medicine for his convection hypothesis for caloric stimulation. Microgravity caloric tests aboard the 1983 SpaceLab 1 mission produced nystagmus results that contradicted the basic premise of Barany's convection theory. In this paper, we present a fluid structural analysis of the caloric stimulation of the lateral semicircular canal. Direct numerical simulations indicate that on earth, natural convection is the dominant mechanism for endolymphatic flow. However, in the microgravity environment of orbiting spacecraft, where buoyancy effects are mitigated, an expansive convection becomes the sole mechanism for producing endolymph motion and cupular displacement. Transient 1 g and microgravity case studies are presented to delineate the different dynamic behaviors of the 1 g and microgravity endolymphatic flows. The associated fluid-structural interactions are also analyzed based on the time evolution of cupular displacements.

Labyrinth/anatomy & histology/physiology↗
NASA NTRS1986

The X-ray background and the evolution of quasars

The contribution of QSOs and active galaxies to the 2-10-keV X-ray background is calculated theoretically, assuming that the X-ray luminosity of each GSO has a power-law time evolution (dL/dt = AL exp alpha) and applying the continuity equation to derive the X-ray luminosity function at arbitrary redshift. The observed X-ray background is shown to require alpha greater than 1.2, ruling out pure luminosity evolution (alpha = 1).

Tucker, W. H.↗
DOE OSTI2022

Population balance models for polymer upcycling: signatures of the mechanism in the molecular weight evolution

Chemical and catalytic upcycling processes could help realize a circular plastics economy, but current models for testing mechanistic hypotheses and designing catalysts remain primitive. This work shows how proposed catalytic mechanisms can be incorporated into population balance models to predict the time evolution of molecular weight distributions. We develop models for homogeneous and heterogeneous catalysts, including catalysts that cut at chain ends and catalysts that cut at random locations. For heterogeneous catalysts, we illustrate the effect of adsorption constants that depend on polymer chain length. Here, we discuss ongoing efforts and challenges in measuring and modeling the time evolving molecular weight distributions in polymer upcycling processes.

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

X-ray spectral and flux variability of the microquasar GRS 1758-258 on timescales from weeks to years

We present the spectral and timing evolution of the persistent black hole X-ray binary GRS 1758-258 based on almost 12 years of observations using the Rossi X-ray Timing Explorer Proportional Counter Array. While the source was predominantly found in the hard state during this time, it entered the thermally dominated soft state seven times. In the soft state GRS 1758-258 shows a strong decline in flux above 3 keV rather than the pivoting flux around 10 keV more commonly shown by black hole transients. In its 3–20 keV hardness intensity diagram, GRS 1758-258 shows a hysteresis of hard and soft state fluxes typical for transient sources in outburst. The RXTE-PCA and RXTE-ASM long-term light curves do not show any orbital modulations in the range of 2–30 d. However, in the dynamic power spectra significant peaks drift between 18.47 and 18.04 d for the PCA data, while less significant signatures between 19 d and 20 d are seen for the ASM data as well as for the Swift/BAT data. We discuss different models for the hysteresis behavior during state transitions as well as possibilities for the origin of the long term variation in the context of a warped accretion disk.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Multiple steady states for characteristic initial value problems

The time dependent, isentropic, quasi-one-dimensional equations of gas dynamics and other model equations are considered under the constraint of characteristic boundary conditions. Analysis of the time evolution shows how different initial data may lead to different steady states and how seemingly anamolous behavior of the solution may be resolved. Numerical experimentation using time consistent explicit algorithms verifies the conclusions of the analysis. The use of implicit schemes with very large time steps leads to erroneous results.

Salas, M. D.↗