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 163 records · Page 9

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗

Gate-Based Quantum Simulation of Gaussian Bosonic Circuits on Exponentially Many Modes

We introduce a framework for simulating, on an ( n + 1 )-qubit quantum computer, the action of a Gaussian bosonic (GB) circuit on a state over 2 n modes. Specifically, we encode the initial bosonic state’s expectation values over quadrature operators (and their covariance matrix) as an input qubit state. This is then evolved by a quantum circuit that effectively implements the symplectic propagators induced by the GB gates. We find families of GB circuits and initial states leading to efficient quantum simulations. For this purpose, we introduce a dictionary that maps between GB and qubit gates such that particle- (non-particle-) preserving GB gates lead to real- (imaginary-) time evolutions at the qubit level. For the special case of particle-preserving circuits, we present a bounded-error-quantum-polynomial time (BQP)-complete GB decision problem, indicating that GB evolutions of Gaussian states on exponentially many modes are as powerful as universal quantum computers. We also perform numerical simulations of an interferometer on ∼ 8 × 10 9 modes, illustrating the power of our framework. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical calculations of the collapse of nonrotating, magnetic gas clouds

Results of the first self-consistent numerical calculations of the dynamic collapse of a magnetized protostellar gas cloud are presented. Symmetry about an axis parallel to the initial magnetic field direction has been assumed, so that the calculations could be performed on a two-dimensional grid. Also, the cloud was taken to be nonrotating and isothermal, and the magnetic field was assumed to remain frozen in to the gas. As starting models for the calculations, gas spheres with uniform density and magnetic field were used. The time evolution of the clouds has been calculated for roughly two initial free-fall times, at which point the central density has increased by a factor of approximately 10,000 to 1,000,000. Several such calculations have been performed for different values of the cloud's initial thermal, magnetic, and gravitational energies. In virtually all cases it is found that, once a flattened core forms in the cloud, the central magnetic field strength, B, varies with gas density, rho, according to (d log B/d log rho) = 1/2. This behavior is independent of the initial energy ratios mentioned above. It is also found that the magnetic field is able to prevent completely the collapse of part of the outer envelope of the cloud.

Scott, E. H.↗

Multiparameter radar study of rainfall: Potential application to area-time integral studies

Multiparameter radars measure one or more additional parameters in addition to the coventional reflectivity factor. The combination of radar observations from a multiparameter radar is used to study the time evolution of rainstorms. A technique is presented to self-consistently compare the area-time integral (ATI) and rainfall volume estimates from convective storms, using two different measurements from a multiparameter radar. Rainfall volumes for the lifetime of individual storms are computed using the reflectivity at S band (10-cm wavelength) as well as one-way specific attenuation at X band (3-cm wavelength). Area-time integrals are computed by summing all areas in each radar snapshot having reflectivities (S band) in excess of a preselected threshold. The multiparameter radar data used in this study were acquired by the National Center for Atmospheric Research (NCAR) CP-2 radar during the Cooperative Huntsville Meteorological Experiment (COHMEX) and the Convection and Precipitation/Electrification Experiment (CaPE), respectively. ATI studies were accomplished in this work using multiparameter radar data acquired during the lifetime of six convective events that occurred in the COHMEX radar coverage area. A case study from the COMHEX field campaign (20 July 1986) was selected to depict the various stages in the evolution of a storm over which the ATI and rainfall volume computations were performed using multiparameter radar data. Another case study from the CaPE field campaign (12 August 1991) was used to demonstrate the evolution of a convective cell based on differential reflectivity observations.

Raghavan, R.↗

Anelastic response of the earth to a dip slip earthquake

The deformation induced by a vertical dip slip earthquake is examined using a variety of rheologic models. In this way the complications of dipping faults are avoided, and the phenomenon of transient peripheral warping is clearly revealed. A thrust fault dipping at 30 deg is investigated, and the important effects of dip and the existence of a slab on the asymmetry of strain pulses propagated into the overthrust and subducted lithosphere are demonstrated. One of the signal results of the study is the essential similarity of the strain patterns for Newtonian and non-Newtonian flow laws: the two rheologies give nearly identical strain field geometries. The principal difference between the two, which is readily observable, is in their time evolution. Relaxation in non-Newtonian rheologies tends to be initially fast, then slow at times that are late in comparison with relaxation in a Newtonian rheology. The possibility of simply recalling the time dependence of a Newtonian solution to obtain an approximate solution to a non-Newtonian problem is demonstrated.

Melosh, H. J.↗

Time-Variable Gravity Signal Due to Extratropic Pacific Water Mass Redistribution

Cox and Chao [2002] reported the detection of a large anomaly in the form of a positive "jump" in the time series of Earth's lowest-degree gravity harmonic J2, or the dynamic oblateness, during 1998. This prompted us to examine the mass redistribution in the global oceans. We report here a seesaw of the sea-surface height (SSH) in the extratropic north + south Pacific basins -- the leading (nonseasonal) EOF/PC mode in SSH derived from the 10-year TOPEX/Poseidon altimetry data in the extratropic Pacific region. The mode underwent a step-like jump with time evolution that match remarkably well with the observed J2 anomaly. However, the magnitude is several times too small to explain the observed J2, even if assuming the SSH jump was all mass-induced (as opposed to any steric effect which causes no time-variable gravity signal). If one accepts the notion that this extratropic Pacific seesaw is part of the geophysical process that produced the observed 1998 J2 anomaly, then this finding suggests strong geophysical connection of the interannual-to-decadal variation of J2 with the Pacific Decadal Oscillation (PDO), as the time series of the above EOF/PC mode is actually a formally defined PDO Index series.

Chao, Benjamin F.↗

An Analysis of the Orbital Distribution of Solid Rocket Motor Slag

The contribution made by orbiting solid rocket motors (SRMs) to the orbital debris environment is both potentially significant and insufficiently studied. A combination of rocket motor design and the mechanisms of the combustion process can lead to the emission of sufficiently large and numerous by-products to warrant assessment of their contribution to the orbital debris environment. These particles are formed during SRM tail-off, or the termination of burn, by the rapid expansion, dissemination, and solidification of the molten Al2O3 slag pool accumulated during the main burn phase of SRMs utilizing immersion-type nozzles. Though the usage of SRMs is low compared to the usage of liquid fueled motors, the propensity of SRMs to generate particles in the 100 m and larger size regime has caused concern regarding their contributing to the debris environment. Particle sizes as large as 1 cm have been witnessed in ground tests conducted under vacuum conditions and comparable sizes have been estimated via ground-based telescopic and in-situ observations of sub-orbital SRM tail-off events. Using sub-orbital and post recovery observations, a simplistic number-size-velocity distribution of slag from on-orbit SRM firings was postulated. In this paper we have developed more elaborate distributions and emission scenarios and modeled the resultant orbital population and its time evolution by incorporating a historical database of SRM launches, propellant masses, and likely location and time of particulate deposition. From this analysis a more comprehensive understanding has been obtained of the role of SRM ejecta in the orbital debris environment, indicating that SRM slag is a significant component of the current and future population.

Horstman, Matthew F.↗

Multiple energetic injections in a strong spike-like solar burst

An intense and fast spike-like solar burst was built up of short time scale structures superimposed on an underlying gradual emission, the time evolution of which shows remarkable proportionality between hard X-ray and microwave fluxes. The finer time structure were best defined at mm-microwaves. At the peak of the event, the finer structures repeat every 30x60ms. The more slowly varying component with a time scale of about 1 second was identified in microwave hard X-rays throughout the burst duration. It is suggested that X-ray fluxes might also be proportional to the repetition rate of basic units of energy injection (quasi-quantized). The relevant parameters of one primary energy release site are estimated both in the case where hard X-rays are produced primarily by thick-target bremsstrahlung, and when they are purely thermal. The relation of this figure to global energy considerations is discussed.

Kaufmann, P.↗

Multiple energetic injections in a strong spike-like solar burst

An intense and fast spike-like solar burst was built up of short time scale structures superimposed on an underlying gradual emission, the time evolution of which shows remarkable proportionality between hard X-ray and microwave fluxes. The finer time structures were best defined at mm-microwaves. At the peak of the event, the finer structures repeat every 30 x 60 ms. The more slowly varying component with a time scale of about 1 second was identified in microwave hard X-rays throughout the burst duration. It is suggested that X-ray fluxes might also be proportional to the repetition rate of basic units of energy injection (quasi-quantized). The relevant parameters of one primary energy release site are estimated both in the case where hard X-rays are produced primarily by thick-target bremsstrahlung, and when they are purely thermal. The relation of this figure to global energy considerations is discussed. Previously announced in STAR as N83-35983

Kaufmann, P.↗

Plasmaspheric Mass Loss and Refilling as a Result of a Magnetic Storm

Using the sounding measurements from the radio plasma imager on IMAGE, and a plasma density inversion algorithm, we derive the plasma density profiles along the magnetic field in a few L-shells every 14 hours at magnetic local noon before, during, and after the March 31,2001 magnetic storm. An empirical model of the plasmaspheric plasma density distribution is derived as a reference using the measurements before the storm. During the storm, the equatorial plasma was substantially depleted in a range of L-shells. The flux tubes were refilled after the storm. The filling ratio, the equatorial plasma density normalized by its quiet-time value before the storm, is introduced to assess the time evolution of the depletion and refilling processes. The depletion, more than two thirds of the quiet time content, appeared to occur rather quickly after the storm onset, as determined by the limited temporal resolution of the measurements. The refilling proceeded, although more slowly than the depletion process, significantly faster than the theoretical prediction of a 3-day time scale. Dynamic structures are observed in situ and confirmed by the extreme ultraviolet imager (EUV) measurements.

Reinisch, B. W.↗

Accreting on the Edge: A Luminosity-Dependent Cyclotron Line in the Be/X-ray Binary 2S 1553-542 Accompanied by Accretion Regimes Transition

Accreting X-ray pulsars (XRPs) undergo luminous X-ray outbursts during which the luminosity-dependent spectral and timing features of the neutron star’s emission can be analyzed in detail, thus shedding light on the accretion regime at work. We took advantage of a monitoring campaign performed with NuSTAR, Swift/XRT, AstroSat and NICER, to follow the Be/X-ray Binary 2S 1553-542 along one of its rare outbursts and trace its spectral and timing evolution. We report the discovery of a luminosity-dependent cyclotron line energy for the first time in this source. The pulse profiles and pulsed fraction also show variability along the outburst, consistently with the interpretation that the source transitions from the sub-critical to the super-critical accretion regime, separated by a critical luminosity ofLcrit≈4×1037erg/s.

X-ray binary stars↗

Universal time scalings of sensitivity in Markovian quantum metrology

Assuming Markovian time evolution of a quantum sensing system, we study the general characterization of the optimal sensitivity scalings with time, under most general quantum control protocols. We allow the estimated parameter to influence both the Hamiltonian as well as the dissipative part of the quantum master equation and focus on the asymptotic-time along with the short-time sensitivity scalings. We find that via simple algebraic conditions (in terms of the Hamiltonian, the jump operators as well as their parameter derivatives), one can characterize the four classes of metrological models that represent: quadratic-linear, quadratic-quadratic, linear-linear, and linear-quadratic time scalings. We also investigate the relevant time scales on which the transition between the two regimes appears. Additionally, we provide universal numerical methods to obtain quantitative bounds on sensitivity that are the tightest that exist in the literature. Simplicity and universality of our results make it suitable for diverse applications in quantum metrology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An optimal smoothing algorithm for Viking spacecraft attitude determination

A fundamental requirement for ground analysis of onboard attitude control systems is the determination of the time evolution of spacecraft angular orientation and velocity, and of external torques acting on the spacecraft. However, at times, there are limitations on the data received on the ground, and essential but missing variables must be estimated. This paper contains a description of the design and implementation of an algorithm for angular position, velocity, and torque estimation. The algorithm has been used to analyze the performance of the Viking Orbiter attitude control system. Typical results of this analysis are presented.

Rodriguez, G.↗

Multiple transonic solutions and a new class of shock transitions in solar and stellar winds

The steady isothermal solar wind equations are shown to admit, under certain circumstances, mutliple transonic solutions when, for example, momentum deposition gives rise to multiplee critical points in the flow. These multiple solutions consist of a continuous solution and solutions which involve shock transitions between critical solutions. The ambiguity arising from the multiplicity of the solutions can be resolved by following the time evolution of a wind profile with one critical point. Results of the numerical integration of the time-dependent equations with momentum addition show that each of these multiple solutions is physically accessible and depends on the rate of change of momentum deposition. These results suggest that standing shocks are likely to be present in the inner solar wind flow.

Habbal, S. R.↗

Time Accurate Unsteady Pressure Loads Simulated for the Space Launch System at a Wind Tunnel Condition

Using the Fully Unstructured Three-Dimensional (FUN3D) computational fluid dynamics code, an unsteady, time-accurate flow field about a Space Launch System configuration was simulated at a transonic wind tunnel condition (Mach = 0.9). Delayed detached eddy simulation combined with Reynolds Averaged Naiver-Stokes and a Spallart-Almaras turbulence model were employed for the simulation. Second order accurate time evolution scheme was used to simulate the flow field, with a minimum of 0.2 seconds of simulated time to as much as 1.4 seconds. Data was collected at 480 pressure taps at locations, 139 of which matched a 3% wind tunnel model, tested in the Transonic Dynamic Tunnel (TDT) facility at NASA Langley Research Center. Comparisons between computation and experiment showed agreement within 5% in terms of location for peak RMS levels, and 20% for frequency and magnitude of power spectral densities. Grid resolution and time step sensitivity studies were performed to identify methods for improved accuracy comparisons to wind tunnel data. With limited computational resources, accurate trends for reduced vibratory loads on the vehicle were observed. Exploratory methods such as determining minimized computed errors based on CFL number and sub-iterations, as well as evaluating frequency content of the unsteady pressures and evaluation of oscillatory shock structures were used in this study to enhance computational efficiency and solution accuracy. These techniques enabled development of a set of best practices, for the evaluation of future flight vehicle designs in terms of vibratory loads.

Alter, Stephen J.↗

The Chaotic Long-Term X-Ray Variability of 4U 1705-44

The low-mass X-ray binary 4U1705-44 exhibits dramatic long-term X-ray time variability with a time-scale of several hundred days. The All-SkyMonitor (ASM) aboard the Rossi X-ray Timing Explorer (RXTE) and the Japanese Monitor of All-sky X-ray Image (MAXI) aboard the International Space Station together have continuously observed the source from 1995 December through 2014 May. The combined ASM-MAXI data provide a continuous time series over 50 times the length of the time-scale of interest. Topological analysis can help us identify 'fingerprints' in the phase space of a system unique to its equations of motion. The Birman-Williams theorem postulates that if such fingerprints are the same between two systems, then their equations of motion must be closely related. The phase-space embedding of the source light curve shows a strong resemblance to the double-welled non-linear Duffing oscillator. We explore a range of parameters for which the Duffing oscillator closely mirrors the time evolution of 4U1705-44. We extract low period, unstable periodic orbits from the 4U1705-44 and Duffing time series and compare their topological information. The Duffing and 4U1705-44 topological properties are identical, providing strong evidence that they share the same underlying template. This suggests that we can look to the Duffing equation to help guide the development of a physical model to describe the long-term X-ray variability of this and other similarly behaved X-ray binary systems.

Rossi X-ray Timing Explorer (RXTE)↗

Physics and evolution of constant opening angle jets using a quasi-one-dimensional magnetohydrodynamic model

We discuss the significance of the assumptions of infinite conductivity and time independence in the context of an ideal MHD model for constant opening angle jets. The model is developed by projecting the MHD equations onto the jet axis. We find that for initially sub-Alfvenic flows (i.e., flows emanating from active galactic nuclei and neutron stars) wind-type solutions exist only when the field lines at the origin are wound up in a direction opposite to the direction of rotation of the compact source. We discuss the possibility that the time evolution of these outflows may be a cycle between breeze- and wind-type solutions as a result of continuous changes in the boundary conditions at the origin due to accretion. We propose that such cycles may explain the apparent one-sideness of some jets, especially the ones for which we cannot use arguments of relativistic beaming. We examine the dependence of the wind-type solutions on the following parameters describing the outflow at the origin: the degree of winding of the field lines, the value of the gas pressure, the polytropic index, the strength of the magnetic field, the value of the rotational velocity, the gravitational potential of the compact object, and the injection velocity. We compare results with results obtained previously, and discuss briefly the qualitative features and physical interpretation of the solutions for outflows emanating from neutron stars and protostars.

Koupelis, Theodoros↗