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 559 records · Page 31

Astrometric Calibration and Performance of the Dark Energy Spectroscopic Instrument Focal Plane

Abstract The Dark Energy Spectroscopic Instrument, consisting of 5020 robotic fiber positioners and associated systems on the Mayall telescope at Kitt Peak, Arizona, is carrying out a survey to measure the spectra of 40 million galaxies and quasars and produce the largest 3D map of the universe to date. The primary science goal is to use baryon acoustic oscillations to measure the expansion history of the universe and the time evolution of dark energy. A key function of the online control system is to position each fiber on a particular target in the focal plane with an accuracy of 11 μ m rms 2D. This paper describes the set of software programs used to perform this function along with the methods used to validate their performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Full Transport General Relativistic Radiation Magnetohydrodynamics for Nucleosynthesis in Collapsars

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above ${Y}_{{\rm{e}}}\sim 0.3$ and preventing third-peak r-process material from being synthesized. We analyze the time evolution of neutrino processes and electron fraction in the disk and present a simple one-dimensional model for the vertical structure that emerges. Finally, we compare our simulation to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

79 ASTRONOMY AND ASTROPHYSICS↗

The Twins Embedding of Type Ia Supernovae. I. The Diversity of Spectra at Maximum Light

In this work, we study the spectral diversity of Type Ia supernovae (SNe Ia) at maximum light using high signal-to-noise spectrophotometry of 173 SNe Ia from the Nearby Supernova Factory. We decompose the diversity of these spectra into different extrinsic and intrinsic components, and we construct a nonlinear parameterization of the intrinsic diversity of SNe Ia that preserves pairings of "twin"SNe Ia. We call this parameterization the "Twins Embedding."Our methodology naturally handles highly nonlinear variability in spectra, such as changes in the photosphere expansion velocity, and uses the full spectrum rather than being limited to specific spectral line strengths, ratios, or velocities. We find that the time evolution of SNe Ia near maximum light is remarkably similar, with 84.6% of the variance in common to all SNe Ia. After correcting for brightness and color, the intrinsic variability of SNe Ia is mostly restricted to specific spectral lines, and we find intrinsic dispersions as low as ~0.02 mag between 6600 and 7200 Å. With a nonlinear three-dimensional model plus one dimension for color, we can explain 89.2% of the intrinsic diversity in our sample of SNe Ia, which includes several different kinds of "peculiar"SNe Ia. A linear model requires seven dimensions to explain a comparable fraction of the intrinsic diversity. We show how a wide range of previously established indicators of diversity in SNe Ia can be recovered from the Twins Embedding. In a companion article, we discuss how these results can be applied to the standardization of SNe Ia for cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

A Machine Learning–Based Approach to Time-series Wave Identification in the Solar Wind

The Wind spacecraft has yielded several decades of high-resolution magnetic field data, a large fraction of which displays small-scale structures. In particular, the solar wind is full of wavelike fluctuations that appear in both the field magnitude and its components. The nature of these fluctuations can be tied to the properties of other structures in the solar wind, such as shocks, that have implications for the time evolution of the solar wind. As such, having a large collection of wave events would facilitate further study of the effects that these fluctuations have on solar wind evolution. Given the large volume of magnetic field data available, machine learning is the most practical approach to classifying the myriad small-scale structures observed. To this end, a subset of Wind data is labeled and used as a training set for a multibranch 1D convolutional neural network aimed at classifying circularly polarized wave modes. Using this algorithm, a preliminary statistical study of 1 yr of data is performed, yielding about 300,000 wave intervals out of about 5,000,000 solar wind intervals. The wave intervals come about more often in the fast solar wind and at higher temperatures, and the number of waves per day is highly periodic. This machine learning–based approach to wave detection has the potential to be a powerful, inexpensive way to catalog waves throughout decades of spacecraft data.

79 ASTRONOMY AND ASTROPHYSICS↗

An Exceptionally Powerful, Radiatively Driven Ultrafast Outflow in the Rapidly Accreting AGN RE J1034+396

We report the analysis of ∼1 Ms of XMM-Newton observations of the rapidly accreting active galactic nucleus RE J1034+396. The 0.3–9 keV EPIC-pn spectra are well described by a model consisting of steep continuum emission from the corona accompanied by relativistically blurred reflection from a highly ionized accretion disk. The source is known to exhibit strong excess soft X-ray emission, which we show is well represented by thermal disk photons Comptonized by a warm plasma spanning the inner accretion flow. Additionally, the EPIC-pn data provide compelling evidence (ΔC ∼ 60 for four additional parameters) for the presence of an ultrafast outflow (UFO) with a line-of-sight velocity v/c = 0.307$^{+0.001}_{-0.005}$, and an emission signature consistent with reflection of the corona from modestly ionized, outflowing gas. The simultaneous 0.5–2.5 keV RGS spectra show clear absorption lines. Modelling of these data confirms the presence of the UFO and constrains its equivalent hydrogen column density, log N H /(atom cm −2 ) = 21.7$^{+0.1}_{-0.2}$. The RGS data also reveal at least two warm absorber components with a modest outflow velocity (1680$^{+40}_{-50}$ km s -1 ). The measured properties and time evolution of the UFO in RE J1034+396 suggest that it is formed from collisionally ionized plasma, launched from the disk surface and accelerated by radiation pressure. The high terminal velocity and substantial absorbing column density imply that the outflow carries sufficient momentum and energy to transform its environment, being capable of driving out essentially all dust and gas it interacts with along the line of sight, even if the AGN were initially surrounded by a Compton-thick absorber.

79 ASTRONOMY AND ASTROPHYSICS↗

Variable X-Ray Reverberation in the Rapidly Accreting Active Galactic Nucleus Ark 564: The Response of the Soft Excess to the Changing Geometry of the Inner Accretion Flow

X-ray reverberation, which exploits the time delays between variability in different energy bands as a function of Fourier frequency, probes the structure of the inner accretion disks and X-ray coronae of active galactic nuclei. We present a systematic X-ray spectroscopic and reverberation study of the high-Eddington-ratio narrow-line Seyfert 1 galaxy Ark 564, using over 900 ks of XMM-Newton and NuSTAR observations spanning 13 yr. The time-averaged spectra can be well described by a model consisting of a coronal continuum, relativistic disk reflection, warm Comptonization, and three warm absorbers. Leveraging the high X-ray brightness of Ark 564, we are able to resolve the time evolution of the spectra and contemporaneous reverberation lags. The soft-band lag relative to the continuum increases with the X-ray flux, while Fe Kα lags are detected in only a subset of epochs and do not correlate with soft lags. Models based on a lamppost corona and reflection from a standard thin disk can broadly reproduce the observed lag-energy spectra of low-flux epochs; however, additional reverberation from the warm Comptonized atmosphere is required to explain the soft lags observed in high-flux epochs. A vertically puffed-up inner disk and a variable, vertically extended corona can better explain the observed evolution of the lags and covariance spectra. Our study underscores the importance of multiepoch, multiband analyses for a comprehensive understanding the inner accretion disk and corona.

Active galactic nuclei↗

Young Stellar Objects, Accretion Disks, and Their Variability with Rubin Observatory LSST

Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multiband filters option available in Rubin LSST, in particular the u, g, r, i filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g., field crowdedness, massive neighbors, metallicity) by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Optimization Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide–Fast–Deep observing mode that will be used to survey the sky over the full duration of the main survey (≈10 yr). The combination of these two modes will be vital to investigate the connection between the inner-disk dynamics and longer-term eruptive variability behaviors, such as those observed on EX Lupi–type objects.

79 ASTRONOMY AND ASTROPHYSICS↗

Multimessenger Characterization of Markarian 501 during Historically Low X-Ray and γ-Ray Activity

We study the broadband emission of Mrk 501 using multiwavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi's Large Area Telescope (LAT), NuSTAR, Swift, GASP-WEBT, and the Owens Valley Radio Observatory. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wave bands, with the highest occurring at X-rays and very-high-energy (VHE) γ-rays. A significant correlation (>3σ) between X-rays and VHE γ-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between the Swift X-Ray Telescope and Fermi-LAT. We additionally find correlations between high-energy γ-rays and radio, with the radio lagging by more than 100 days, placing the γ-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE γ-rays from mid-2017 to mid-2019 with a stable VHE flux (>0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2 yr long low state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED toward the low state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic Error Cancellation in Analog Quantum Simulation

Analog quantum simulation is a promising path towards solving classically intractable problems in many-body physics on near-term quantum devices. However, the presence of noise limits the size of the system and the length of time that can be simulated. In our work, we consider an error model in which the actual Hamiltonian of the simulator differs from the target Hamiltonian we want to simulate by small local perturbations, which are assumed to be random and unbiased. We analyze the error accumulated in observables in this setting and show that, due to stochastic error cancellation, with high probability the error scales as the square root of the number of qubits instead of linearly. We explore the concentration phenomenon of this error as well as its implications for local observables in the thermodynamic limit. Moreover, we show that stochastic error cancellation also manifests in the fidelity between the target state at the end of time-evolution and the actual state we obtain in the presence of noise. This indicates that, to reach a certain fidelity, more noise can be tolerated than implied by the worst-case bound if the noise comes from many statistically independent sources.

Analog quantum simulation↗

Cloudy-sky contributions to the direct aerosol effect

The radiative forcing of the aerosol–radiation interaction can be decomposed into clear-sky and cloudy-sky portions. Two sets of multi-model simulations within Aerosol Comparisons between Observations and Models (AeroCom), combined with observational methods, and the time evolution of aerosol emissions over the industrial era show that the contribution from cloudy-sky regions is likely weak. A mean of the simulations considered is 0.01±0.1 W m-2. Multivariate data analysis of results from AeroCom Phase II shows that many factors influence the strength of the cloudy-sky contribution to the forcing of the aerosol–radiation interaction. Overall, single-scattering albedo of anthropogenic aerosols and the interaction of aerosols with the short-wave cloud radiative effects are found to be important factors. A more dedicated focus on the contribution from the cloud-free and cloud-covered sky fraction, respectively, to the aerosol–radiation interaction will benefit the quantification of the radiative forcing and its uncertainty range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The 4DEnVar-based weakly coupled land data assimilation system for E3SM version 2

Abstract. A new weakly coupled land data assimilation (WCLDA) system based on the four-dimensional ensemble variational (4DEnVar) method is developed and applied to the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2). The dimension-reduced projection four-dimensional variational (DRP-4DVar) method is employed to implement 4DVar using the ensemble technique instead of the adjoint technique. With an interest in providing initial conditions for decadal climate predictions, monthly mean anomalies of soil moisture and temperature from the Global Land Data Assimilation System (GLDAS) reanalysis from 1980 to 2016 are assimilated into the land component of E3SMv2 within the coupled modeling framework with a 1-month assimilation window. The coupled assimilation experiment is evaluated using multiple metrics, including the cost function, assimilation efficiency index, correlation, root-mean-square error (RMSE), and bias, and compared with a control simulation without land data assimilation. The WCLDA system yields improved simulation of soil moisture and temperature compared with the control simulation, with improvements found throughout the soil layers and in many regions of the global land. In terms of both soil moisture and temperature, the assimilation experiment outperforms the control simulation with reduced RMSE and higher temporal correlation in many regions, especially in South America, central Africa, Australia, and large parts of Eurasia. Furthermore, significant improvements are also found in reproducing the time evolution of the 2012 US Midwest drought, highlighting the crucial role of land surface in drought lifecycle. The WCLDA system is intended to be a foundational resource for research to investigate land-derived climate predictability.

58 GEOSCIENCES↗

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗

Intense windstorms in the northeastern United States

Abstract. Windstorms are a major natural hazard in many countries. The objective of this study is to identify and characterize intense windstorms during the last 4 decades in the US Northeast and determine both the sources of cyclones responsible for these events and the manner in which those cyclones differ from the cyclone climatology. The windstorm detection is based on the spatial extent of locally extreme wind speeds at 100 m height from the ERA5 reanalysis database. During the top 10 windstorms, wind speeds exceed their local 99.9th percentile over at least one-third of land-based ERA5 grid cells in this high-population-density region of the USA. Maximum sustained wind speeds at 100 m during these windstorms range from 26 to over 43 ms−1, with wind speed return periods exceeding 6.5 to 106 years (considering the top 5 % of grid cells during each storm). Property damage associated with these storms, with inflation adjusted to January 2020, ranges from USD 24 million to over USD 29 billion. Two of these windstorms are linked to decaying tropical cyclones, three are Alberta clippers, and the remaining storms are Colorado lows. Two of the 10 re-intensified off the east coast, leading to development of nor'easters. These windstorms followed frequently observed cyclone tracks but exhibit maximum intensities as measured using 700 hPa relative vorticity and mean sea level pressure that is 5–10 times the mean values for cyclones that followed similar tracks over this 40-year period. The time evolution of wind speeds and concurrent precipitation for those windstorms that occurred after the year 2000 exhibit good agreement with in situ ground-based and remote sensing observations, plus storm damage reports, indicating that the ERA5 reanalysis data have a high degree of fidelity for large, damaging windstorms such as these. A larger pool of the top 50 largest windstorms exhibit evidence of only weak serial clustering, which is in contrast to the relatively strong serial clustering of windstorms in Europe.

58 GEOSCIENCES↗

Quantum Simulation of Field Theories Without State Preparation

We propose an algorithm for computing real-time observables using a quantum processor while avoiding the need to prepare the full quantum state. This reduction in quantum resources is achieved by classically sampling configurations in imaginary-time using standard lattice field theory. These configurations can then be passed to quantum processor for time-evolution. This method encounters a signal-to-noise problem which we characterize, and we demonstrate the application of standard lattice QCD methods to mitigate it.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

From Atom to Engine: Understanding Fundamental Effects of Structure on Combustion Using Tandem Experiment and Computation

The ability to tie structural features to a fuel candidate’s ignition properties provides a path for rational design of advantaged fuels that facilitate higher-efficiency, lower-emitting combustion in engines. The time evolution of a fuel’s radical population directly correlates with global combustion parameters, and because the generated radicals depend on the initial structure and reaction pathways of the fuel molecule, it follows that the molecular structure of a fuel has a direct impact on this cascade, and therefore salient ignition-tied fuel properties. This work aims to understand the impact of structural features in the context of global combustion properties for direct application to emerging engine technologies. Specifically, we highlight the ability to identify potential fuel blendstocks for specific engine strategies by examining and elucidating the radical cascade tied to ignition. A representative subset of alcohols possessing a variety of structural features (e.g. branching, chain length) was selected for study using a robust combination of experimental and computational methods. The approach begins with identifying short-lived radicals and likely abstraction sites upon pyrolysis in a microreactor coupled with photoionization mass spectrometry for direct detection of product species. In tandem, electronic structure calculations were performed to calculate the relevant potential energy surfaces for fuel pyrolysis as well as the subsequent O2 addition-isomerization reactions, which govern ignition behavior. Furthermore, the theoretical energetic and ro-vibrational data are used to calculate rate constants for direct ignition modeling to explore how structural effects impact NTC behavior. The microreactor data combined with electronic structure theory calculations provide a well-informed picture of the reactions relevant to ignition for each fuel in this study, as well as the most likely radicals generated under engine conditions. This thorough fundamental picture of fuel decomposition and ignition chemistry is used to correlate the observed ignition behavior to molecular structural effects.

advanced biofuels↗

CHEMNODE: CHEMICAL KINETICS SOLVER APPROACH BASED ON NEURAL ORDINARY DIFFERENTIAL EQUATIONS

This software contains an algorithm to model and predict the time evolution of chemical kinetics in reacting flow simulations using neural ordinary differential equations (NODEs). The approach works by using artificial neural networks to predict the chemical source terms. It achieves this by using forward-mode automatic differentiation and the Levenberg-Marquardt algorithm to adjust the neural network parameters, such that the discrepancies between the actual and predicted species profiles are minimized.

OWOYELE, OPEOLUWA↗

Molecular Dynamics Study of Creep Deformation in Nickel-based Superalloy

The present study employs molecular dynamics simulations of Ni-based superalloy to investigate the creep behavior under uniaxial compression test. Dislocation dynamics is analyzed for the nickel-based single crystal superalloy with the presence of void and with varying the distribution of gamma-prime phase The results show that multiple-void systems are more prone to yield than single-void systems and single-void systems are more prone to yield than the system withoutvoid. From the simulations, it has been determined that the creep mechanism in Ni/Ni3Al is subject to change on the applied stress depending on the distribution of gamma-prime phases change. Dislocation behavior is also observed in addition to the time evolution of dislocation nucleation and growth to quantify the propensity and mechanism of creep in nickel-based single crystal superalloy. The studies confirmed that dislocation length prior to plastic deformation decreased as the number of voids increased. In addition, the deformation behavior for finegrained nickel-based polycrystalline superalloy upon grain size variation is also investigated. Our findings show that both the movement mechanisms of the grains and grain boundaries and the relationship between grain size and yield strength are highly affected by the grain size of the materials. All the databases have then been classified to accurately predict the target class vi machine learning using the data mining tool WEKA.

creep, nickel-based superalloy, molecular dynamics↗

Stout Smearing on a Quantum Computer

Smearing of gauge-field configurations in lattice field theory improves the results of lattice simulations by suppressing high energy modes from correlation functions. In quantum simulations, high kinetic energy eigenstates are introduced when the time evolution operator is approximated such as Trotterization. While improved Trotter product formulae exist to reduce the errors, they have diminishing accuracy returns with respect to resource costs. Therefore having an algorithm that has fewer resources than an improved Trotter formula is desirable. In this work I develop a representation agnostic method for quantum smearing and show that it reduces the coupling to high energy modes in the discrete nonabelian gauge theory $D_4$

Gustafson, Erik J.↗