Search NASA⌕ Search

SEARCH · Search NASA

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

Gravity and count probabilities in an expanding universe

The time evolution of nonlinear clustering on large scales in cold dark matter, hot dark matter, and white noise models of the universe is investigated using N-body simulations performed with a tree code. Count probabilities in cubic cells are determined as functions of the cell size and the clustering state (redshift), and comparisons are made with various theoretical models. We isolate the features that appear to be the result of gravitational instability, those that depend on the initial conditions, and those that are likely a consequence of numerical limitations. More specifically, we study the development of skewness, kurtosis, and the fifth moment in relation to variance, the dependence of the void probability on time as well as on sparseness of sampling, and the overall shape of the count probability distribution. Implications of our results for theoretical and observational studies are discussed.

Bouchet, Francois R.↗
NASA NTRS1993

Analysis of the effects of mean local node-crossing time on the evolution of Sun-synchronous orbits

An investigation of the effect of mean local node-crossing time on the evolution of Sun-synchronous orbits was undertaken during Phase-A orbit analysis for the National Oceanic and Atmospheric Administration (NOAA) O,P,Q environmental spacecraft. That analysis added to the growing body of evidence that individual Sun-synchronous missions, at differing node-crossing times, experience nodal drift rates that can differ in both magnitude and direction. A Sun-synchronous orbit is obtained by means of a nodal drift rate approximating the 0.9856-degree-per-day apparent precession of the position of the mean Sun. This drift rate is achieved through the interaction of the orbital semimajor axis and inclination in Earth's geopotential field. Influencing perturbations include atmospheric drag and, most important, the effects of solar gravitation on inclination. The present analysis examines a series of Sun-synchronous orbits with mean local node-crossing times at 1-hour intervals from 6 a.m. to 6 p.m. It considers the fixed geometry of each orbital plane with respect to both the Sun and the diurnal atmospheric bulge, then analyzes the influence of these features upon the evolution of the semimajor axix and inclination and thus upon the rate of the nodal drift in the course of 1 year.

Jordan, P.↗
NASA NTRS1993

The dynamics of large-scale cyclogenesis over the North Pacific Ocean

Wintertime cases of persistent cyclonic flow abnormalities over the North Pacific were studied with the aim of identifying the primary mechanisms for persistent anomaly development. Processes influencing the time evolution of eddy potential enstrophy and potential vorticity anomalies were examined using intensive diagnostic analyses. Analysis indicates that the primary source region for large-scale development is located over the mid-Pacific jet exit region. Potential enstrophy analysis shows that eddy enstrophy increases result from local conversions between the large-scale eddy and the climatological-mean flow. The occurrence of a form of baroclinic 'self development' during large-scale intensification is suggested by potential vorticity inversions. Baroclinic and barotropic processes contribute positively to the large-scale developments. Results indicate that the primary mechanism for the developments is a large-scale instability of the three-dimensional time-mean flow, and suggest that nonmodel transient growth plays a significant role during development.

Black, Robert X.↗
DOE OSTI2024

Quarkonium transport in weakly and strongly coupled plasmas

We report on progress in the nonperturbative understanding of quarkonium dynamics inside a thermal plasma. The time evolution of small-size quarkonium is governed by two-point correlation functions of chromoelectric fields dressed with an adjoint Wilson line, known in this context as generalized gluon distributions (GGDs). The GGDs have been calculated in both weakly and strongly coupled plasmas by using perturbative and holographic methods. Strikingly, the results of our calculations for a strongly coupled plasma indicate that the quarkonium dissociation and recombination rates vanish in the transport descriptions that assume quarkonium undergoes Markovian dynamics. However, this does not imply that the dynamics is trivial. As a starting point to explore the phenomenological consequences of the result at strong coupling, we show a calculation of the Y(1 S ) formation probability in time-dependent perturbation theory. This is a first step towards the development of a transport formalism that includes non-Markovian effects, which, depending on how close the as of yet undetermined nonperturbative QCD result of the GGDs is to the strongly coupled N = 4 SYM result, could very well dominate over the Markovian ones in quark-gluon plasma produced at RHIC and the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Ozone transport during a cut-off low event studied in the frame of the TOASTE program

A study of ozone transfer to the troposphere has been performed during two phases of the evolution of a cut-off low using both ozone vertical profiles and objective analysis of the ECMWF to compute potential vorticity distributions and air mass trajectories. Ozone profiles were measured by a ground based lidar system at the Observatoire de Haute Provence (OHP, 43 deg 55 N, 5 deg 42 E). A stratospheric ozone transport into the troposphere has been observed during a tropopause fold which occurred at the beginning of the cut-off low formation and during the erosion phase of the cut-off low. From the estimate of the maximum ozone content transferred to the troposphere, both mechanisms have the same order of magnitude of influence on the ozone flux to the troposphere. On a time scale of a few days, the correlation is very good between the potential vorticity and the ozone time evolution in the vicinity of the upper level frontal system.

Ancellet, G.↗
DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Time dependent simulation of cosmic-ray shocks including Alfven transport

Time evolution of plane, cosmic-ray modified shocks was simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven wave energy transport equations. Nonlinear feedback from the cosmic-rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the time scale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic-rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
DOE OSTI2024

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2023

DPOD2020: A DORIS Extension of the ITRF2020 for Precise Orbit Determination

As one of the tracking systems used to determine orbits of the altimeter mission satellites (such as TOPEX/Poseidon, Envisat, Jason-1/-2/-3, CryoSat-2, Saral/Altika, Sentinel-3A/-3B, HY-2A/C/D, Jason-CS/Sentinel-6A, SWOT), DORIS (Doppler Orbitography Radiopositionning Integrated by Satellite) allows to determine positions and velocities of tracking stations that define a stable reference for the estimation of the precise orbits and thus are fundamental for the quality of the altimeter data and derived mean sea level products. Due to the time evolution of the DORIS ground network, some stations included in the 2020 realization of the International Terrestrial Reference Frame (ITRF2020) have been decommissioned and since 2021.0 a few new stations were added to the tracking network. Therefore, to satisfy operational requirements for POD (Precise Orbit Determination) and routine delivery of geodetic products, the International DORIS Service (IDS) regularly updates the DPOD (DORIS terrestrial reference frame for Precise Orbit Determination). The DPOD solutions include mean positions and velocities of all the DORIS stations since 1993.0 derived from the stacking of the latest IDS weekly combined series aligned to the current ITRF. In this paper, we first present the stacking process of the DPOD2020 version 1.0. Then, we address the validation procedure of the DPO2020 including comparison with ITRF2020 and POD tests. For eighty percent of all the time segments of all the DORIS stations, the station position differences between DPOD2020 version 1.0 and ITRF2020 are smaller than ten millimeters. The major position differences between these two solutions are associated with the DORIS sites either localized in the South Atlantic Anomaly region or with time spans smaller than one year. Compared to DPOD2014, the DPOD2020 shows reduction of the main statistics of the DORIS-to-DORIS tie residuals (differences between the estimated and measured ties). The POD tests showed similar results for DPOD2020 and DPOD2014 for most of the altimetric satellites (TOPEX/Poseidon, Jason-1, CryoSat-2, Jason-3). In addition, we observed better POD results with DPOD2020 for the latest altimetric satellite Sentinel-3A as well as a slight degradation for Jason-2. That degradation was fully explained by slightly worst results for the stations localized in the South Atlantic Anomaly region.

DORIS↗
NASA NTRS2004

Torque Balances on the Taylor Cylinders in the Geomagnetic Data Assimilation

In this presentation we report on our continuing effort in geomagnetic data assimilation, aiming at understanding and predicting geomagnetic secular variation on decadal time scales. In particular, we focus on the effect of the torque balances on the cylindrical surfaces in the core co-axial with the Earth's rotation axis (the Taylor cylinders) on the time evolution of assimilated solutions. We use our MoSST core dynamics,model and observed geomagnetic field at the Earth's surface derived via Comprehensive Field Model (CFM) for the geomagnetic data assimilation. In our earlier studies, a model solution is selected randomly from our numerical database. It is then assimilated with the observations such that the poloidal field possesses the same field tomography on the core-mantel boundary (CMB) continued downward from surface observations. This tomography change is assumed to be effective through out the outer core. While this approach allows rapid convergence between model solutions and the observations, it also generates sevee numerical instabilities: the delicate balance between weak fluid inertia and the magnetic torques on the Taylor cylinders are completely altered. Consequently, the assimilated solution diverges quickly (in approximately 10% of the magnetic free-decay time in the core). To improve the assimilation, we propose a partial penetration of the assimilation from the CMB: The full-scale modification at the CMB decreases linearly and vanish at an interior radius r(sub a). We shall examine from our assimilation tests possible relationships between the convergence rate of the model solutions to observations and the cut-off radius r(sub a). A better assimilation shall serve our nudging tests in near future.

Kuang, Weijia↗
DOE OSTI2025

Control simulations of many-body quantum systems by a synergism of discrete real-time learning and optimal control theory

We present a self-consistent algorithm for optimal control simulations of many-body quantum systems. The algorithm features a two-step synergism that combines discrete real-time machine learning (DRTL) with Quantum Optimal Control Theory (QOCT) using the time-dependent Schrödinger equation. Specifically, in step (1), DRTL is employed to identify a compact working space (i.e., the important portion of the Hilbert space) for the time evolution of the many-body quantum system in the presence of a control field (i.e., the initial or previously updated field), and in step (2), QOCT utilizes the DRTL-determined working space to find a newly updated control field for a chosen objective. Steps 1 and 2 are iterated until a self-consistent control objective value is reached such that the resulting optimal control field yields the same targeted objective value when the corresponding working space is systematically enlarged. Furthermore, to demonstrate this two-step self-consistent DRTL-QOCT synergistic algorithm, we perform optimal control simulations of strongly interacting 1D as well as 2D Heisenberg spin systems. In both scenarios, only a single spin (at the left end site for 1D and the upper left corner site for 2D) is driven by the time-dependent control fields to create an excitation at the opposite site as the target. It is found that, starting from all spin-down zero excitation states, the synergistic method is able to identify working spaces and convergence of the desired controlled dynamics with just a few iterations of the overall algorithm. In the cases studied, the dimensionality of the working space scales only quasi-linearly with the number of spins.

Artificial neural networks↗
DOE OSTI2025

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗
NASA NTRSDate not supplied

Impact of DORISx Post Seismic Deformation and Seasonal Corrections on the DPOD2020 Solution

As one of the tracking systems used for the altimeter missions (such as TOPEX/Poseidon, Envisat, Jason-1/-2/-3, CryoSat-2, Saral/AltiKa, Sentinel-3A/-3B, HY-2A/C/D, Jason-CS/Sentinel-6A, SWOT), the position of the DORIS tracking stations provides a fundamental reference for the estimation of the precise orbits and so, by extension is fundamental for the quality of the altimeter data analysis and derived products. Due to the time evolution of the DORIS ground network, stations included in ITRF2020 may have been decommissioned and stations that were added to the tracking network after 2021.0 are, by definition, not part of ITRF2020. Therefore, to satisfy operational requirements for precise orbit determination and routine delivery of geodetic products, the International DORIS Service maintains the DORIS terrestrial reference frame for Precise Orbit Determination (DPOD) solutions. Since 2016, the DPOD products are computed by the IDS Combination Center (CC) as a DORIS cumulative position and linear velocity solution aligned to the latest available ITRF and using the latest IDS weekly combined series. Mid-2023, the IDS CC started the computation of the second version of the DPOD2020 based on the IDS 20 weekly solutions from 1993.0 to 2023.0. With that new solution, the IDS CC will include Post Seismic Deformation (PSD) corrections for some DORIS sites such as Socorro Island. These DPOD PSD corrections are determined from the analysis of DORIS coordinate time series only. The second version will also include optional annual and semi-annual station corrections. The presentation will describe the analysis used for extracting the PSD and seasonal DORIS station corrections, compare such corrections to those derived for ITRF2020, and evaluate the impact these new corrections have on station positioning and satellite POD.

Guilhem Moreaux↗
DOE OSTI2026

Efficient Measurement-Driven Eigenenergy Estimation with Classical Shadows

Quantum algorithms exploiting real-time evolution under a target Hamiltonian have demonstrated remarkable efficiency in extracting key spectral information. However, the broader potential of these methods, particularly beyond ground-state calculations, is underexplored. In this work, we introduce the framework of multiobservable dynamic mode decomposition (MODMD), which combines the observable dynamic mode decomposition (DMD), a measurement-driven eigensolver tailored for near-term implementation, with classical shadow tomography. MODMD leverages random scrambling in the classical shadow technique to construct, with exponentially reduced resource requirements, a signal subspace that encodes rich spectral information. Notably, we replace typical Hadamard-test circuits with a protocol designed to predict low-rank observables, thereby broadening the use of classical shadow tomography for predicting many low-rank observables. We establish theoretical guarantees on the spectral approximation from MODMD, taking into account distinct sources of error. In the ideal case, we prove that the spectral error scales as exp (−Δ⁢𝐸⁢𝑡 max ), where Δ⁢𝐸 is the Hamiltonian spectral gap and 𝑡 max is the maximal simulation time. This analysis provides a rigorous justification of the rapid convergence observed across simulations. To demonstrate the utility of our framework, we consider its application to fundamental tasks, such as determining the low-lying, i.e., ground or excited, energies of representative many-body systems. Our work paves the path for efficient designs of measurement-driven algorithms on near-term and early fault-tolerant quantum devices.

quantum algorithms & computation↗
NASA NTRS1995

The Dripping Handrail Model: Transient Chaos in Accretion Systems

We define and study a simple dynamical model for accretion systems, the "dripping handrail" (DHR). The time evolution of this spatially extended system is a mixture of periodic and apparently random (but actually deterministic) behavior. The nature of this mixture depends on the values of its physical parameters - the accretion rate, diffusion coefficient, and density threshold. The aperiodic component is a special kind of deterministic chaos called transient chaos. The model can simultaneously exhibit both the quasiperiodic oscillations (QPO) and very low frequency noise (VLFN) that characterize the power spectra of fluctuations of several classes of accretion systems in astronomy. For this reason, our model may be relevant to many such astrophysical systems, including binary stars with accretion onto a compact object - white dwarf, neutron star, or black hole - as well as active galactic nuclei. We describe the systematics of the DHR's temporal behavior, by exploring its physical parameter space using several diagnostics: power spectra, wavelet "scalegrams," and Lyapunov exponents. In addition, we note that for large accretion rates the DHR has periodic modes; the effective pulse shapes for these modes - evaluated by folding the time series at the known period - bear a resemblance to the similarly- determined shapes for some x-ray pulsars. The pulsing observed in some of these systems may be such periodic-mode accretion, and not due to pure rotation as in the standard pulsar model.

Young, Karl↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Adaptive variational quantum dynamics simulations with compressed circuits and fewer measurements

The adaptive variational quantum dynamics simulation (AVQDS) method performs real-time evolution of quantum states using automatically generated parametrized quantum circuits that often contain substantially fewer gates than Trotter circuits. Here we report an improved version of the method, which we call AVQDS(T), by porting the tiling efficient trial circuits with rotations implemented simultaneously technique. The algorithm adaptively adds layers of disjoint unitary gates to the ansatz circuit so as to keep the McLachlan distance, a measure of the accuracy of the variational dynamics, below a fixed threshold. Here we perform benchmark noiseless AVQDS(T) simulations of quench dynamics in local spin models and compare with an alternative adaptive variational approach on quantum resource requirement. Quantum dynamical simulations implementing realistic noise channels are also reported. Finally, we propose a way to substantially alleviate the measurement overhead of AVQDS(T) while maintaining high accuracy by synergistically integrating quantum circuit calculations on quantum processing units with classical calculations using, e.g., tensor networks to evaluate the quantum geometric tensor. We showcase that this approach enables AVQDS(T) to deliver more accurate results than simulations using a fixed ansatz of comparable final depth for a significant time duration with fewer quantum resources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Digital quantum magnetism on a trapped-ion quantum computer

Digital quantum matter—realized when discrete quantum gates approximate continuous time evolution—is susceptible to heating into chaotic, structureless states. If digitization errors are adequately suppressed, a long-lived transient regime of approximately energy-conserving dynamics can be observed on gate-based quantum computers. Conservation of energy, in turn, enables the exploration of a wide variety of complex behaviours observed in equilibrium systems, ranging from the non-trivial microscopic origins of thermalization itself to the stabilization of effective models hosting exotic emergent properties. Here we use Quantinuum’s H2 quantum computer to simulate digitized dynamics of the quantum Ising model, suppressing digitization errors well enough to observe thermalization on timescales that severely challenge classical simulation methods. Relaxation of an inhomogeneous state reveals an emergent hydrodynamics owing to approximate energy conservation and we compute the associated diffusion constant. By reprogramming our simulations to take place on a triangular lattice with periodic boundary conditions, we observe thermalization consistent with emergent gauge and topological constraints resulting from lattice frustration. Furthermore, our results were enabled by continued advances in two-qubit gate quality (native partial entangler fidelities of 99.94(1)%) and establish digital quantum computers as powerful tools for studying (effectively) continuous-time dynamics.

Information theory and computation↗