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

Diagonal isometric form for tensor network states in two dimensions

Isometric tensor network states (isoTNS) generalize the isometric form of the one-dimensional matrix product states to tensor networks in two and higher dimensions. Here, we introduce an alternative isometric form for isoTNS by incorporating auxiliary tensors to represent the orthogonality hypersurface. We implement the time evolving block decimation algorithm on this new isometric form and benchmark the method by computing ground states and the real time evolution of the transverse field Ising model in two dimensions on large square lattices of up to 1250 sites. Our results demonstrate that isoTNS can efficiently capture the entanglement structure of two-dimensional area law states. The short-time dynamics is also accurately reproduced even at the critical point. Our isoTNS formulation further allows for a natural extension to different lattice geometries, such as the honeycomb or kagome lattice.

2-dimensional systems↗
DOE OSTI2025

Anatomy of information scrambling and decoherence in the integrable Sachdev-Ye-Kitaev model

The growth of information scrambling, captured by out-of-time-order correlation functions (OTOCs), is a central indicator of the nature of many-body quantum dynamics. Here, we compute analytically the complete time dependence of the OTOC for an integrable Sachdev-Ye-Kitaev (SYK) model, 𝑁 Majoranas with random two-body interactions of infinite range, coupled to a Markovian bath at finite temperature. In the limit of no coupling to the bath, the time evolution of scrambling experiences different stages. For 𝑡≲√𝑁, after an initial polynomial growth, the OTOC approaches saturation in a power-law fashion with oscillations superimposed. At 𝑡∼√𝑁, the OTOC reverses trend and starts to decrease linearly in time. The reason for this linear decrease is that, despite being a subleading 1/𝑁 effect, the OTOC in this region is governed by the spectral form factor of the antisymmetric couplings of the SYK model. The linear decrease stops at 𝑡∼2⁢𝑁, the Heisenberg time, where saturation occurs. The effect of the environment is an overall exponential decay of the OTOC for times longer than the inverse of the coupling strength to the bath. Here, the oscillations at 𝑡≲√𝑁 indicate lack of thermalization—a desired feature for better performance of quantum information devices.

Open quantum systems & decoherence↗
DOE OSTI2026

AT2025ulz and S250818k: Zooming in with the Hubble Space Telescope

AT2025ulz is an optical/near-infrared transient discovered during follow-up of the candidate gravitational wave (GW) event S250818k. Its young age (≲1 day), rapid decline, and strong color evolution over the first 48 hr classify it as a potential kilonova candidate. In this work, we present the results of our observing campaign, carried out with the Gran Telescopio Canarias and the Hubble Space Telescope (HST). Although the early time evolution of AT2025ulz resembles some aspects of a kilonova, its rapid onset (∼3 hr after the GW trigger) and luminosity (a factor of ∼5 brighter than AT2017gfo in g band) are difficult to reproduce. Only a small subset of our kilonova models matches its multicolor light curve, and the inferred ejecta mass is uncomfortably large given the low chirp mass (≲0.87 M ⊙ ) of the GW candidate. HST observations place the transient within a nearby (z = 0.08489) spiral galaxy with ongoing star formation and measure a color (F336W − F160W ≈ 1.4 mag) that is too blue to match with a kilonova. Our data support the classification of AT2025ulz as a supernova, initially undergoing a shock-cooling phase and later entering its photospheric phase and spectroscopically identified via its broad absorption features.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

Dust Coagulation in Protoplanetary Accretion Disks

The time evolution of dust particles in circumstellar disk-like structures around protostars and young stellar objects is discussed. In particular, we consider the coagulation of grains due to collisional aggregation. The coagulation of the particles is calculated by solving numerically the non-linear Smoluchowski equation. The different physical processes leading to relative velocities between the grains are investigated. The relative velocities may be induced by Brownian motion, turbulence and drift motion. Starting from different regimes which can be identified during the grain growth we also discuss the evolution of dust opacities. These opacities are important for both the derivation of the circumstellar dust mass from submillimeter/millimeter continuum observations and the dynamical behavior of the disks. We present results of our numerical studies of the coagulation of dust grains in a turbulent protoplanetary accretion disk described by a time-dependent one-dimensional (radial) alpha-model. For several periods and disk radii, mass distributions of coagulated grains have been calculated. From these mass spectra, we determined the corresponding Rosseland mean dust opacities. The influence of grain opacity changes due to dust coagulation on the dynamical evolution of a protostellar disk is considered. Significant changes in the thermal structure of the protoplanetary nebula are observed. A 'gap' in the accretion disk forms at the very frontier of the coagulation, i.e., behind the sublimation boundary in the region between 1 and 5 AU.

Schmitt, W.↗
NASA NTRS1991

Evolution of magnetic flux ropes associated with flux transfer events and interplanetary magnetic clouds

Spacecraft observations suggest that flux transfer events and interplanetary magnetic clouds may be associated with magnetic flux ropes which are magnetic flux tubes containing helical magnetic field lines. In the magnetic flux ropes, the azimuthal magnetic field is superposed on the axial field. The time evolution of a localized magnetic flux rope is studied. A two-dimensional compressible MHD simulation code with a cylindrical symmetry is developed to study the wave modes associated with the evolution of flux ropes. It is found that in the initial phase both the fast magnetosonic wave and the Alfven wave are developed in the flux rope. After this initial phase, the Alfven wave becomes the dominant wave mode for the evolution of the magnetic flux rope and the radial expansion velocity of the flux rope is found to be negligible. Numerical results further show that even for a large initial azimuthal component of the magnetic field, the propagation velocity along the axial direction of the flux rope remains the Alfven velocity. It is also found that the localized magnetic flux rope tends to evolve into two separate magnetic ropes propagating in opposite directions. The simulation results are used to study the evolution of magnetic flux ropes associated with flux transfer events observed at the earth's dayside magnetopause and magnetic clouds in the interplanetary space.

Wei, C. Q.↗
DOE OSTI2026

Kinetic Model of HoxEFU reduction by NADH [SWR-26-087]

This repository is used to release code generated for manuscripts on the Photosynthetic Energy Transduction core program. This code simulates the reduction of HoxEFU by NADH. The electro transfer rate constants for the simulation are specified in the .csv files. The two .csv files correspond tot he two models described in Dawson et al. Cell. Rep. Phys. Sci. 2026. The code utilizes a chemical master equation, a set of differential equations, defining the time evolution of the oxidation and reduction kinetics of NAD+, NADH, a FMN flavin, and a set of iron sulfur clusters. The kinetics of HoxEFU reduction by NADH are evaluated by numerical integration of the chemical master equation using a variable-time-step Runge-Kutta algorithm.

Dahl, Peter [National Laboratory of the Rockies (N↗
DOE OSTI2026

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗
NASA NTRS1993

Obliquity variation in a Mars climate evolution model

The existence of layered terrain in both polar regions of Mars is strong evidence supporting a cyclic variation in climate. It has been suggested that periods of net deposition have alternated with periods of net erosion in creating the layered structure that is seen today. The cause for this cyclic climatic behavior is variation in the annually averaged latitudinal distribution of solar insolation in response to obliquity cycles. For Mars, obliquity variation leads to major climatological excursion due to the condensation and sublimation of the major atmospheric constituent, CO2. The atmosphere will collapse into the polar caps, or existing caps will rapidly sublimate into the atmosphere, dependent upon the polar surface heat balance and the direction of the change in obliquity. It has been argued that variations in the obliquity of Mars cause substantial departures from the current climatological values of the surface pressure and the amount of CO2 stored in both the planetary regolith and polar caps. In this new work we have modified the Haberle et al. model to incorporate variable obliquity by allowing the polar and equatorial insolation to become functions of obliquity, which we assume to vary sinusoidally in time. As obliquity varies in the model, there can be discontinuities in the time evolution of the model equilibrium values for surface pressure, regolith, and polar cap storage. The time constant, tau r, for the regolith to find equilibrium with the climate is estimated--depending on the depth, thermal conductivity, and porosity of the regolith--between 10(exp 4) and 10(exp 6) yr. Thus, using 2000-yr timesteps to move smoothly through the 0.1250 m.y. obliquity cycles, we have an atmosphere/regolith system that cannot be assumed in equilibrium. We have dealt with this problem by limiting the rate at which CO2, can move between the atmosphere and regolith, mimicking the diffusive nature and effects of the temperature and pressure waves, by setting the time rate of change of regolith storage proportional to the difference between equilibrium storage and current storage.

Tyler, D.↗
NASA NTRS1989

The radiative signature of the shock acceleration process in active galactic nuclei

A time-dependent model in which plasma enters a shock region and undergoes first-order Fermi acceleration is considered. The time evolution of the particle distribution and the associated synchrotron emission has been calculated for an energy-dependent diffusion coefficient. When the emitting region is not spatially resolved, the spectrum is a power law with an exponential turnover in frequency. With increasing time, the power law and the turnover extend to higher frequencies. To test the model, observations at three or more frequencies should be made of time-resolved flux increases from the optically thin nonthermal emission region in variable active galactic nuclei. Existing data that can be used to test the model are scarce but offer qualitative support.

Bregman, Joel N.↗
NASA NTRS1999

Drop Ejection From an Oscillating Rod

The dynamics of a drop of a Newtonian liquid that is pendant from or sessile on a solid rod that is forced to undergo time-periodic oscillations along its axis is studied theoretically. The free boundary problem governing the time evolution of the shape of the drop and the flow field inside it is solved by a method of lines using a finite element algorithm incorporating an adaptive mesh. When the forcing amplitude is small, the drop approaches a limit cycle at large times and undergoes steady oscillations thereafter. However, drop breakup is the consequence if the forcing amplitude exceeds a critical value. Over a wide range of amplitudes above this critical value, drop ejection from the rod occurs during the second oscillation period from the commencement of rod motion. Remarkably, the shape of the interface at breakup and the volume of the primary drop formed are insensitive to changes in forcing amplitude. The interface shape at times close to and at breakup is a multi-valued function of distance measured along the rod axis and hence cannot be described by recently popularized one-dimensional approximations. The computations show that drop ejection occurs without the formation of a long neck. Therefore, this method of drop formation holds promise of preventing formation of undesirable satellite droplets.

Wilkes, E. D.↗
NASA NTRS2011

GRB 091127: The Cooling Break Race on Magnetic Fuel

Using high-quality, broad-band afterglow data for GRB 091127, we investigate the validity of the synchrotron fireball model for gamma-ray bursts, and infer physical parameters of the ultra-relativistic outflow. Methods. We used multi-wavelength (NIR to X-ray) follow-up observations obtained with GROND simultaneously in the g' r' t' i' z' JH filters and the XRT onboard the Swift satellite in the 0.3 to 10 keY energy range. The resulting afterglow light curve is of excellent accuracy with relative photometric errors as low as 1 %, and the spectral energy distribution (SED) is well-sampled over 5 decades in energy. These data present one of the most comprehensive observing campaigns for a single GRB afterglow and allow us to test several proposed emission models and outflow characteristics in unprecedented detail. Results. Both the multi-color light curve and the broad-band SED of the afterglow of GRB 091127 show evidence of a cooling break moving from high to lower energies. The early light curve is well described by a broken power-law, where the initial decay in the optical/NlR wavelength range is considerably flatter than at X-rays. Detailed fitting of the time-resolved SED shows that the break is very smooth with a sharpness index of 2.2 +/- 0.2, and evolves towards lower frequencies as a power-law with index -1.23 +/- 0.06. These are the first accurate and contemporaneous measurements of both the sharpness of the spectral break and its time evolution. Conclusions. The measured evolution of the cooling break (V(sub c) varies as t(sup -1.2) is not consistent with the predictions of the standard model, wherein V(sub c) varies as t(sup -05) is expected. A possible explanation for the observed behavior is a time dependence of the microphysical parameters, in particular the fraction of the total energy in the magnetic field epsilon(sub Beta). This conclusion provides further evidence that the standard fireball model is too simplistic, and time-dependent micro-physical parameters may be required to model the growing number of well-sampled afterglow light curves.

Filgas, R.↗
NASA NTRS2014

Timing and Flux Evolution of the Galactic Center Magnetar SGR J1745-2900

We present the X-ray timing and spectral evolution of the Galactic Center magnetar SGR J1745-2900 (SGR*) for the first 4 months post-discovery using data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift observatories. Our timing analysis reveals a large increase in the magnetar spin-down rate by a factor of 2.6 plus or minus 0.07 over our data span. We further show that the change in spin evolution was likely coincident with a bright X-ray burst observed in 2013 June by Swift, and if so, there was no accompanying discontinuity in the frequency. We find that the source 3 to 10 kiloelectronvolt flux has declined monotonically by a factor of approximately 2 over an 80-day period post-outburst accompanied by an approximately 20 percent decrease in the source's blackbody temperature, although there is evidence for both flux and kiloteslas having leveled off. We argue that the torque variations are likely to be magnetospheric in nature and will dominate over any dynamical signatures of orbital motion around Sgr A*.

Stars: magnetic field↗
NASA NTRS1982

Time-dependent solutions of multimode convection equations

Truncated modal equations are used to study the time evolution of thermal convection. In the Boussinesq approximation these nonlinear equations are obtained by expanding the fluctuating velocity and temperature fields in a finite set of planforms of the horizontal coordinates. Numerical studies dealing with two or three modes with triad interactions are discussed. Rich time dependence was found in these cases: periodic and aperiodic solutions can be obtained, along with various steady solutions. Three-mode solutions reproduce the qualitative appearance of spoke-pattern convection as observed in experiments at high Prandtl numbers. Though the values of the periods of the time-dependent solutions do not agree with those of the experiments, their variation with Rayleigh number compares favorably. Except at the highest Rayleigh number considered (10,000,000), the theoretical Nusselt numbers agree well with experiment.

Toomre, J.↗
DOE OSTI2024

Inferring fusion nuclear burnwidths with low gain photomultiplier impulse response functions

When an inertial confinement fusion implosion is compressed, it maintains thermonuclear density and temperatures for a very short time scale, about 100 ps. The Gamma Reaction History diagnostic measures the time evolution of the fusion burn, but its temporal resolution is limited by the use of a photomultiplier tube (PMT) to amplify the photon signal. Multichannel plate-based PMTs have a fast (~120 ps) full-width at half-max impulse response function (IRF), but the time scale is similar to the incoming physics signal. An analysis routine is used to remove the effect of the PMT IRF and infer the incident fusion burnwidth. With the National Ignition Facility achieving ignition and creating much brighter signals, the PMTs are run at gains three orders of magnitude lower than nominal operation. Calibration at these settings shows the PMT IRFs get ~15% wider. Taking the gain-dependent IRF can affect the inferred nuclear burnwidths by up to ~15%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of their predictions.

Speziale, Charles G.↗
NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of thier predictions.

Speziale, Charles G.↗
NASA NTRS1990

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of their predictions.

Speziale, Charles G.↗