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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quenching time-scales in the IllustrisTNG simulation

ABSTRACT The time-scales for galaxy quenching offer clues to its underlying physical drivers. We investigate central galaxy quenching time-scales in the IllustrisTNG 100-1 simulation, their evolution over time, and the pre-quenching properties of galaxies that predict their quenching time-scales. Defining quenching duration τq as the time between crossing specific star formation rate (sSFR) thresholds, we find that ${\sim} 40{{\ \rm per\ cent}}$ of galaxies quench rapidly with τq < 1 Gyr, but a substantial tail of galaxies can take up to 10 Gyr to quench. Furthermore, 29 per cent of galaxies that left the star-forming main sequence (SFMS) more than 2 Gyr ago never fully quench by z = 0. While the median τq is fairly constant with epoch, the rate of galaxies leaving the SFMS increases steadily over cosmic time, with the rate of slow quenchers being dominant around z ∼ 2–0.7. Compared to fast quenchers (τq < 1 Gyr), slow-quenching galaxies (τq > 1 Gyr) were more massive, had more massive black holes, had larger stellar radii, and accreted gas with higher specific angular momentum (AM) prior to quenching. These properties evolve little by z = 0, except for the accreting gas AM for fast quenchers, which reaches the same high AM as the gas in slow quenchers. By z = 0, slow quenchers also have residual star formation in extended gas rings. Using the expected relationship between stellar age gradient and τq for inside-out quenching we find agreement with Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) integral field unit (IFU) observations. Our results suggest the accreting gas AM and potential well depth determine the quenching time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Magnetohydrodynamic Riemann problem and the structure of the magnetic reconnection layer

We present a complete solution for a set of magnetohydrodynamic (MHD) Riemann problems in which the upstream and downstream states have the same total pressure, and in which the normal component of the magnetic field is very small. These solutions are pertinent to subfast flows in the earth's magnetic tail and near the magnetopause. In a coplanar situation a family of solutions exists that depend on two parameters as well as on dissapation mechanisms. In the parallel case the transverse magnetic field either does not change direction or changes the direction twice by involving two intermediate shocks. In the antiparallel case an intermediate shock is always required, except when the solution consists of two switch-off shocks. In a noncoplanar case the solution is not self-similar as a function of x/t, but continues to evolve. At early times the evolution is similar to the coplanar case. In general two time-dependent intermediate shocks are required to rotate the magnetic fields. The velocity shear has a strong effect on the Riemann solution. In some cases no Riemann solution can exist because of the cavitation caused by the slow refraction waves. The calculated magnetopause structure resembles the observed structure for northward interplanetary magnetic field (IMF). However, for southward IMF, the MHD result shows the existence of a depletion layer, which is not supported by observations. We also show that on the magnetosheath side, the Walen relation, which is exact for a rotational discontinuity, can also be well satisfied by a slow shock, an intermediate shock, or the head of a slow rarefaction wave.

Wu, C. C.↗
NASA NTRS2008

Mergers of Non-spinning Black-hole Binaries: Gravitational Radiation Characteristics

We present a detailed descriptive analysis of the gravitational radiation from black-hole binary mergers of non-spinning black holes, based on numerical simulations of systems varying from equal-mass to a 6:1 mass ratio. Our primary goal is to present relatively complete information about the waveforms, including all the leading multipolar components, to interested researchers. In our analysis, we pursue the simplest physical description of the dominant features in the radiation, providing an interpretation of the waveforms in terms of an implicit rotating source. This interpretation applies uniformly to the full wavetrain, from inspiral through ringdown. We emphasize strong relationships among the l = m modes that persist through the full wavetrain. Exploring the structure of the waveforms in more detail, we conduct detailed analytic fitting of the late-time frequency evolution, identifying a key quantitative feature shared by the l = m modes among all mass-ratios. We identify relationships, with a simple interpretation in terms of the implicit rotating source, among the evolution of frequency and amplitude, which hold for the late-time radiation. These detailed relationships provide sufficient information about the late-time radiation to yield a predictive model for the late-time waveforms, an alternative to the common practice of modeling by a sum of quasinormal mode overtones. We demonstrate an application of this in a new effective-one-body-based analytic waveform model.

Baker, John G.↗
NASA NTRS1976

Interstellar cloud evolution and the abundance of formaldehyde

The time scale for essentially complete conversion of C(+) to CO in interstellar clouds can be comparable to, or greater than, dynamical time scales for evolution, therefore suggesting steady state time independent abundances to be inappropriate. The solutions for the time-dependent carbon chemistry in dense clouds, with density not less than 500/cu cm, indicate that significant amounts of neutral carbon will be present throughout a cloud's lifetime. These nonequilibrium values of C I can explain the relatively large abundances observed for formaldehyde, isotopes of carbon monoxide, and other trace molecules

Langer, W. D.↗
NASA NTRS1983

Initial-value problems for Rossby waves in a shear flow with critical level

The time-dependent evolution of sheared Rossby waves starting from an initial disturbance is studied for the simple case in which the shear is uniform. The uniform-shear assumption allows explicit solutions to be obtained which are useful in addressing the issue of the long-time asymptotic approach to normal modes and in assessing the relative importance of viscosity, nonlinearity and time-dependence in the evolution of Rossby waves in the presence of critical layers.

Tung, K. K.↗
DOE OSTI2026

Quark flavor equilibration of the quark-gluon plasma

The early stage of a heavy-ion collision is marked by rapid entropy production and the transition from a gluon saturated initial condition to a plasma of quarks and gluons that evolves hydrodynamically. However, during the early times of the hydrodynamic evolution, the chemical composition of the QCD medium is still largely unknown. We present a study of quark chemical equilibration in the (Q)GP using a novel model of viscous hydrodynamic evolution in partial chemical equilibrium. Motivated by the success of gluon saturated initial condition models, we initialize the QCD medium as a completely gluon dominated state. Local quark production during the hydrodynamic phase is then simulated through the evolution of time-dependent fugacities for each independent quark flavor, with the timescales set as free parameters to compare different rates of equilibration. We present the results of complete heavy-ion collision simulations using this partial chemical equilibrium model, and show the effects on hadronic and electromagnetic observables. In particular, we show that the development of flow is sensitive to the equilibration timescale, providing an empirical way to probe the chemical equilibration of the QCD medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Asymptotic errors in adiabatic evolution

The adiabatic theorem in quantum mechanics implies that if a system is in a discrete eigenstate of a Hamiltonian and the Hamiltonian evolves in time arbitrarily slowly, the system will remain in the corresponding eigenstate of the evolved Hamiltonian. Understanding corrections to the adiabatic result that arise when the evolution of the Hamiltonian is slow—but not arbitrarily slow—has become increasingly important, especially since adiabatic evolution has been proposed as a method of state preparation in quantum computing. Here, this paper identifies two regimes, an adiabatic regime in which corrections are generically small and can depend on details of the evolution throughout the path, and a hyperadiabatic regime in which the error is given by a form similar to an asymptotic expansion in the inverse of the evolution time with the coefficients depending principally on the behavior at the endpoints. However, the error in this hyperadiabatic regime is neither given by a true asymptotic series nor solely dependent on the endpoints: the coefficients combine the contributions from both endpoints, with relative phase factors that depend on the average spectral gaps along the trajectory, multiplied by the evolution time. The central result of this paper is to identify a quantity, referred to as the typical error, which is obtained by appropriately averaging the error over evolution times that are small compared to the evolution time itself. This typical error is characterized by an asymptotic series and depends solely on the endpoints of the evolution, remaining independent of the details of the intermediate evolution.

adiabatic approximation↗
NASA NTRS1992

A first-order global model of Late Cenozoic climatic change: Orbital forcing as a pacemaker of the ice ages

The development of a theory of the evolution of the climate of the earth over millions of years can be subdivided into three fundamental, nested, problems: (1) to establish by equilibrium climate models (e.g., general circulation models) the diagnostic relations, valid at any time, between the fast-response climate variables (i.e., the 'weather statistics') and both the prescribed external radiative forcing and the prescribed distribution of the slow response variables (e.g., the ice sheets and shelves, the deep ocean state, and the atmospheric CO2 concentration); (2) to construct, by an essentially inductive process, a model of the time-dependent evolution of the slow-response climatic variables over time scales longer than the damping times of these variables but shorter than the time scale of tectonic changes in the boundary conditions (e.g., altered geography and elevation of the continents, slow outgassing, and weathering) and ultra-slow astronomical changes such as in the solar radiative output; and (3) to determine the nature of these ultra-slow processes and their effects on the evolution of the equilibrium state of the climatic system about which the above time-dependent variations occur. All three problems are discussed in the context of the theory of the Quaternary climate, which will be incomplete unless it is embedded in a more general theory for the fuller Cenozoic that can accommodate the onset of the ice-age fluctuations. We construct a simple mathematical model for the Late Cenozoic climatic changes based on the hypothesis that forced and free variations of the concentration of atmospheric greenhouse gases (notably CO2), coupled with changes in the deep ocean state and ice mass, under the additional 'pacemaking' influence of earth-orbital forcing, are primary determinants of the climate state over this period. Our goal is to illustrate how a single model governing both very long term variations and higher frequency oscillatory variations in the Pleistocene can be formulated with relatively few adjustable parameters.

Saltzman, Barry↗
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Dust particles from comets and asteroids collected at the Earth's orbit: Parent-daughter relationships

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not well known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of Asteroids in the main belt. It is believed that one may combine lab analysis of the physics and chemistry of captured particles with orbital data in order to identify comet and asteroid parent bodies. It is possible to use the collected orbits of the dust to connect with its source in two ways. One is to consider the long time orbit evolution of the dust under Poynting-Robertson drag. The other is to look at the prompt orbit change of dust from comets onto trajectories that intersect the earth's orbit. In order to characterize the orbits of dust particles evolved over a long period of time, a study of its orbital evolution was undertaken. Various parameters associated with these dust orbits as they cross the Earth's orbit were considered in order to see if one may discriminate between particles evolved from comets and asteroids. The method was to calculate by a numerical procedure the orbits of dust particles after they left their parent bodies. It appears that as the particles pass the Earth's orbit, asteroidal grains and cometary grains can be differentiated on the basis of their measured orbital eccentricities even after much planetary perturbation. Broad parent daughter associations can be made on this basis from measurement of their trajectories intercepted in earth orbit.

Jackson, A. A.↗
NASA NTRS1980

Fluid mechanics simulation of fog formation associated with polluted atmosphere produced by energy related fuel combustion

It is noted that large quantities of atmospheric aerosols with composition SO4(-2), NO3(-1), and NH4(+1) have been detected in highly industrialized areas. Most aerosol products come from energy-related fuel combustion. Fluid mechanics simulation of both microphysical and macrophysical processes is considered in studying the time dependent evolution of the saturation spectra of condensation nuclei associated with polluted and clean atmospheres during the time periods of advection fog formation. The results demonstrate that the condensation nuclei associated with a polluted atmosphere provide more favorable conditions than condensation nuclei associated with a clean atmosphere to produce dense advection fog, and that attaining a certain degree of supersaturation is not necessarily required for the formation of advection fog having condensation nuclei associated with a polluted atmosphere.

Hung, R. J.↗
NASA NTRS1991

High-precision timing of millisecond pulsars. II - Astrometry, orbital evolution, and eclipses of PSR 1957 + 20

Pulse time-of-arrival measurements have been acquired for PSR 1957 + 20 on 61 days over a 2.7 yr interval, including the times of 44 eclipse disappearance or reappearance events. The timing data provide a reliable measurement of the pulsar's proper motion, 29 +/- 3 mas/yr, at position angle 216 +/- 6 deg; this direction is in excellent agreement with that suggested by the cometlike shape of a surrounding optical emission nebula. The rate of change of the orbital period has been measured: a surprisingly large (-3.9 +/- 0.9) x 10 exp -11, implying that significant orbital evolution must take place on a time scale of 30 Myr. Finally, the present observations extend the eclipse measurements to higher frequencies and provide better statistics for these highly variable phenomena. The eclipse data are used to compute a column-density profile of free electrons in the outer parts of the eclipsing region, and thereby provide an important glimpse into the material being ablated from the companion star.

Ryba, M. F.↗
DOE OSTI2026

Universal progression of structure and dynamics in colloidal nanocrystal gels during salt-accelerated aging

Controlling the structure and function of colloidal gels requires a detailed understanding of how the various components govern network formation and aging. In particular, molecular additives like salts are widely used to tune interparticle interactions, yet their influence on gelation pathways in complex systems such as colloidal nanocrystal gels remains inadequately understood. Here, we investigate how noncoordinating salts modulate the evolution of gels formed using chemically linked tin-doped indium oxide nanocrystals. Through combined structural, dynamic, and kinetic analyses, we demonstrate that increasing salt concentration accelerates gelation. When rescaled by salt-dependent characteristic times, the evolution collapses onto universal trajectories, revealing a time-salt superposition principle. The universality extends across length scales, suggesting a consistent salt-dependent mechanism that controls both local structuring and macroscopic network formation. This observed salt modulation of structure and dynamics provides a predictive basis for controlling the kinetics of nonequilibrium nanocrystal gel assembly, enhancing the rational design of functional nanomaterials with tunable properties.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Advection fog formation in a polluted atmosphere

Large quantities of atmospheric aerosols with compositions SO4(2-), NO3(-) and NH4(+) have been detected in highly industrialized areas. The major portions of aerosol products are the results of energy-related fuel combustion. Both microphysical and macrophysical processes are considered in investigating the time-dependent evolution of the saturation spectra of condensation nuclei associated with both polluted and clean atmospheres during the time periods of advection fog formation. The results show that the condensation nuclei associated with a polluted atmosphere provide more favorable conditions than condensation nuclei associated with a clean atmosphere to produce dense advection fog, and that attaining a certain degree of supersaturation is not necessarily required for the formation of advection fog with condensation nuclei associated with a polluted atmosphere for monodisperse distribution.

Hung, R. J.↗
DOE OSTI2022

Neural-network learning of SPOD latent dynamics

Here, we aim to reconstruct the latent space dynamics of high dimensional, quasi-stationary systems using model order reduction via the spectral proper orthogonal decomposition (SPOD). The proposed method is based on three fundamental steps: in the first, once that the mean flow field has been subtracted from the realizations (also referred to as snapshots), we compress the data from a high-dimensional representation to a lower dimensional one by constructing the SPOD latent space; in the second, we build the time-dependent coefficients by projecting the snapshots containing the fluctuations onto the SPOD basis and we learn their evolution in time with the aid of recurrent neural networks; in the third, we reconstruct the high-dimensional data from the learnt lower -dimensional representation. The proposed method is demonstrated on two different test cases, namely, a compressible jet flow, and a geophysical problem known as the Madden-Julian Oscillation. An extensive comparison between SPOD and the equivalent POD-based counterpart is provided and differences between the two approaches are highlighted. The numerical results suggest that the proposed model is able to provide low rank predictions of complex statistically stationary data and to provide insights into the evolution of phenomena characterized by specific range of frequencies. The comparison between POD and SPOD surrogate strategies highlights the need for further work on the characterization of the interplay of error between data reduction techniques and neural network forecasts.

97 MATHEMATICS AND COMPUTING↗