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 721 records · Page 40

ITeM: Independent temporal motifs to summarize and compare temporal networks

We report networks are a fundamental and flexible way of representing various complex systems. Many domains such as communication, citation, procurement, biology, social media, and transportation can be modeled as a set of entities and their relationships. Temporal networks are a specialization of general networks where every relationship occurs at a discrete time. The temporal evolution of such networks is as important to understand as the structure of the entities and relationships. We present the Independent Temporal Motif (ITeM) to characterize temporal graphs from different domains. ITeMs can be used to model the structure and the evolution of the graph. In contrast to existing work, ITeMs are edge-disjoint directed motifs that measure the temporal evolution of ordered edges within the motif. For a given temporal graph, we produce a feature vector of ITeM frequencies and the time it takes to form the ITeM instances. We apply this distribution to measure the similarity of temporal graphs. We show that ITeM has higher accuracy than other motif frequency-based approaches. We define various ITeM-based metrics that reveal salient properties of a temporal network. We also present importance sampling as a method to efficiently estimate the ITeM counts. We present a distributed implementation of the ITeM discovery algorithm using Apache Spark and GraphFrame. We evaluate our approach on both synthetic and real temporal networks.

97 MATHEMATICS AND COMPUTING↗

Adaptive Control for Load-Following of Boiling Water Reactors Part I: Linear Systems and Fully-Observable Dynamics

Automation control is a key strategy to improve the economic competitiveness of nuclear power plants. Not only does it help reduce operational costs, but it also extends the value proposition of these plants to nontraditional markets, including unattended operations in remote villages and space. However, the dynamics of the operating environments of nuclear reactors are subject to changes over time, and there are no widely adopted methods to ensure that the automation strategy will remain effective over the extended durations required for these applications. Adaptive control is a discipline that offers the possibility to accommodate such changes online. However, it relies on mathematical assumptions that must be respected to ensure robustness and reliability. In this work, we derive an adaptive control formulation for linear systems in which all states are observable and apply it to an instance of load-follow operation for Boiling Water Reactors. We assumed uncertainty in two factors: the temperature coefficient and the control rod worth, both of which are affected over time by the evolution of the nuclear reactor core environment. With an arbitrary penalty factor of 5, we found that the mean absolute and integral time absolute errors can be reduced by more than 90%, underscoring the strength of adaptive control. To extend the application to more challenges, different uncertainties and load-follow trajectories, as well as new formulations that include non-linearity and partial observability, are currently being developed.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Plateau de Bure High-z Blue Sequence Survey 2 (PHIBSS2): Search for Secondary Sources, CO Luminosity Functions in the Field, and the Evolution of Molecular Gas Density through Cosmic Time

We report on the results of a search for serendipitous sources in CO emission in 110 cubes targeting CO(2 − 1), CO(3 − 2), and CO(6 − 5) at z ∼ 1–2 from the second Plateau de Bure High-z Blue Sequence Survey (PHIBSS2). The PHIBSS2 observations were part of a 4 yr legacy program at the IRAM Plateau de Bure Interferometer aimed at studying early galaxy evolution from the perspective of molecular gas reservoirs. We present a catalog of 67 candidate secondary sources from this search, with 45 of the 110 data cubes showing sources in addition to the primary target that appear to be field detections, unrelated to the central sources. This catalog includes redshifts, line widths, and fluxes, as well as an estimation of their reliability based on their false-positive probability. We perform a search in the 3D Hubble Space Telescope/CANDELS catalogs for the secondary CO detections and tentatively find that ∼64% of these have optical counterparts, which we use to constrain their redshifts. Finally, we use our catalog of candidate CO detections to derive the CO(2 − 1), CO(3 − 2), CO(4 − 3), CO(5 − 4), and CO(6 − 5) luminosity functions over a range of redshifts, as well as the molecular gas mass density evolution. Despite the different methodology, these results are in very good agreement with previous observational constraints derived from blind searches in deep fields. They provide an example of the type of “deep-field” science that can be carried out with targeted observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Revisiting shear stress tensor evolution: Nonresistive magnetohydrodynamics with momentum-dependent relaxation time

This study aims to develop second-order relativistic viscous magnetohydrodynamics (MHD) derived from kinetic theory within an extended relaxation time approximation (momentum/energy dependent) for the collision kernel. The investigation involves a detailed examination of shear stress tensor evolution equations and associated transport coefficients. The Boltzmann equation is solved using a Chapman-Enskog-like gradient expansion for a charge-conserved conformal system, incorporating a momentum-dependent relaxation time. The derived relativistic nonresistive, viscous second-order MHD equations for the shear stress tensor reveal significant modifications in the coupling with dissipative charge current and magnetic field due to the momentum dependence of the relaxation time. By utilizing a power law parametrization to quantify the momentum dependence of the relaxation time, the anisotropic magnetic field-dependent shear coefficients in the Navier-Stokes limit have been investigated. The resulting viscous coefficients are seen to be sensitive to the momentum dependence of the relaxation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermal Behavior During the Selective Laser Melting Process of Ti-6Al-4V Powder in the Point Exposure Scan Pattern

Currently, there are two main scan patterns, including the continuous exposure scan pattern and the point exposure scan pattern, during the selective laser melting (SLM) process. The point exposure scan pattern allows a three-dimensional (3-D) printer to build finer detail features as a static molten pool that is more stable than a dynamic one. However, there has been limited theoretical research on the thermal behavior characteristics during the process in the point exposure scan pattern. Therefore, in this study, the simulation of thermal behavior during SLM of Ti-6Al-4V powder in the point exposure scan pattern was performed. The temperature evolution behavior of different positions and the effects of exposure time on the temperature evolution behavior of different positions, temperature distributions, and dimensions of the molten pool were investigated. The results showed that the direct exposure position and unexposed position had significantly different temperature evolution behaviors under a given condition, and the changed exposure time had the greatest influence on the direct exposure position compared with unexposed position. Moreover, the thermal accumulation effect of a former exposure point on a later one decreased with increasing exposure time. In addition, with the increase of exposure time, the maximum temperature of the molten pool was enhanced and the surface morphology of the molten pool changed from an approximate ellipse to an approximate circle. Besides, the molten pool dimensions were found to increase with exposure time, which indicated that the exposure time played an important role in the stability of the molten pool and the metallurgical bonding in the process. Furthermore, the dimensions of the molten pool and metallurgical bonding in the cross-sectional view were obtained through experiments. Good agreement was obtained when comparing the simulated results with the experimental ones.

36 MATERIALS SCIENCE↗

Efficient Berry phase calculation via adaptive variational quantum computing approach

We present an adaptive variational quantum algorithm to estimate the Berry phase accumulated by a nondegenerate ground state under cyclic, adiabatic evolution of a time-dependent Hamiltonian. Our method leverages cyclic adiabatic evolution of the Hamiltonian and employs adaptive variational quantum algorithms for state preparation and evolution, optimizing circuit efficiency while maintaining high accuracy. We benchmark our approach on dimerized Fermi–Hubbard chains with four sites, demonstrating precise Berry phase simulations in both noninteracting and interacting regimes. Our results show that circuit depths reach up to 106 layers for noninteracting systems and increase to 279 layers for interacting systems due to added complexity. In addition, we demonstrate the robustness of our scheme across a wide range of parameters governing adiabatic evolution and variational algorithms. These findings highlight the potential of adaptive variational quantum algorithms for advancing quantum simulations of topological materials and computing geometric phases in strongly correlated systems.

Mootz, Martin [Ames Laboratory (AMES), Ames, IA (U↗

The Most Metal-poor Stars in the Inner Bulge

The bulge is the oldest component of the Milky Way. Since numerous simulations of Milky Way formation have predicted that the oldest stars at a given metallicity are found on tightly bound orbits, the Galaxy’s oldest stars are likely metal-poor stars in the inner bulge with small apocenters (i.e., R {sub apo} ≲ 4 kpc). In the past, stars with these properties have been impossible to find due to extreme reddening and extinction along the line of sight to the inner bulge. We have used the mid-infrared metal-poor star selection of Schlaufman & Casey (2014) on Spitzer/Galactic Legacy Infrared Mid-Plane Survey Extraordinaire data to overcome these problems and target candidate inner bulge metal-poor giants for moderate-resolution spectroscopy with Anglo-Australian Telescope/AAOmega. We used those data to select three confirmed metal-poor giants ([Fe/H] = −3.15, −2.56, −2.03) for follow-up high-resolution Magellan/Magellan Inamori Kyocera Echelle spectroscopy. A comprehensive orbit analysis using Gaia DR2 astrometry and our measured radial velocities confirms that these stars are tightly bound inner bulge stars. We determine the elemental abundances of each star and find high titanium and iron-peak abundances relative to iron in our most metal-poor star. We propose that the distinct abundance signature we detect is a product of nucleosynthesis in the Chandrasekhar-mass thermonuclear supernova of a CO white dwarf accreting from a helium star with a delay time of about 10 Myr. Even though chemical evolution is expected to occur quickly in the bulge, the intense star formation in the core of the nascent Milky Way was apparently able to produce at least one Chandrasekhar-mass thermonuclear supernova progenitor before chemical evolution advanced beyond [Fe/H]∼−3.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbulence-driven flow dynamics in general axisymmetric toroidal geometry

The present work gives the equations governing the generation of toroidally axisymmetric flows by turbulent Reynolds and Maxwell stresses in finite aspect ratio, general cross section tokamak plasmas. Inclusion of the divergence-free flow constraint in lowest order changes the time scale for evolution of the poloidal flows driven by turbulence by substantial factors. In the pedestal region for present day machines, comparing to earlier cylindrical models, the time scale evaluated using a large aspect ratio circular cross section model can be two orders of magnitude longer while the present, general geometry result can be about one order of magnitude longer. Inclusion of gyroviscosity in the calculation shows that the only lowest order radial velocity fluctuations that enter the problem are those due to fluctuating E×B flows. Toroidal and poloidal flow effects on the toroidally axisymmetric flows are inextricably coupled due to the neoclassical poloidal viscosity. Accordingly, the physics is inherently three dimensional and measurements of all three velocity components are required to obtain the information needed to quantitatively test the theory. The parallel and angular momentum equations for the lowest order, toroidally axisymmetric flows look like radial transport equations when the turbulence is included. The turbulence terms provide the radial transport fluxes. In standard neoclassical theory, the parallel flow equation is local on each flux surface; there is no radial derivative term. However, adding turbulence gives a way, in principle, for radial transport to lead to poloidal flows that deviate from the neoclassical prediction. As a result, the inclusion of the Maxwell stress provides a mechanism for MHD fluctuations to alter the toroidally axisymmetric flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of powder layer thickness on the microstructural development of additively manufactured SS316

By Increasing the powder layer thickness we could increase material deposition rates and reduce build times during laser powder bed fusion (LPBF). However, increasing the powder layer thickness could also negatively affect the local melting and subsequent microstructure of the printed material. This study systematically investigated, for the first time, the defect evolution and microstructural changes of LPBF stainless steel 316 (SS316) processed by using various powder layer thicknesses of 40, 80, and 120 μm while also varying the laser dwell time for the pulsed laser. Through image analysis, the pore density was found to increase with increasing powder layer thickness, particularly when the dwell time was too high or too low. The maximum achievable density (>99%) was found over a wide processing window when using the smallest powder layer thickness of 40 μm between dwell times of 60 and 100 μs. Increasing the powder layer thicknesses to 80 and 120 μm resulted in maximum densities of 98.7 and 96.8%, respectively, but the window for acceptable laser dwell times that could achieve these densities narrowed considerably. A microstructural analysis of the melt pools was performed to measure the melt pool depths and widths, both of which increased with increasing dwell time. However, increasing the powder layer thickness did not affect the melt pool depth and only minorly affected the width. A sub-grain cellular structure distinguished the melt pool boundaries. The cell size increased with increasing dwell time and decreasing powder layer thickness. Moreover, the cell size was used to calculate a cooling rate that had a magnitude of 10 7 K/s and increased with increasing powder layer thickness.

36 MATERIALS SCIENCE↗

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗

Similarities behind the high- and low- α disc: small intrinsic abundance scatter and migrating stars

The detailed age-chemical abundance relations of stars measure time-dependent chemical evolution. These trends offer strong empirical constraints on nucleosynthetic processes, as well as the homogeneity of star-forming gas. Characterizing chemical abundances of stars across the Milky Way over time has been made possible very recently, thanks to surveys like Gaia , APOGEE, and Kepler . Studies of the low- α disc have shown that individual elements have unique age–abundance trends and the intrinsic dispersion around these relations is small. In this study, we examine and compare the age distribution of stars across both the high and low- α disc and quantify the intrinsic dispersion of 16 elements around their age–abundance relations at [Fe/H] = 0 using APOGEE DR16. Here, we examine the age–metallicity relation and visualize the temporal and spatial distribution of disc stars in small chemical cells. We find: (1) the high- α disc has shallower age–abundance relations compared to the low- α disc, but similar median intrinsic dispersions of ~0.03 dex; (2) turnover points in the age-[Fe/H] relations across radius for both the high- and low- α disc. The former constrains the mechanisms that set similar intrinsic dispersions, regardless of differences in the enrichment history, for stars in both disc, and the latter indicates the presence of radial migration in both disc. Our study is accompanied by an age catalogue for 64 317 stars in APOGEE derived using THE CANNON with a median uncertainty of 1.5 Gyr (26 per cent; APO-CAN stars), and a red clump catalogue of 22 031 stars with a contamination rate of 2.7 per cent.

79 ASTRONOMY AND ASTROPHYSICS↗

Operation above the Greenwald density limit in high performance DIII-D negative triangularity discharges

The density limit in strongly-shaped negative triangularity (NT) discharges is studied experimentally in the DIII-D tokamak. Record-high Greenwald fractions f G are obtained, using gas puff injection only, with values up to near 2, where f G is defined as the ratio of the line-averaged density over n G = I p /(π a 2 ), with I p [MA] the plasma current and a[m] the plasma minor radius. A clear higher operational limit with higher auxiliary power is also demonstrated, with the ohmic density limit about two times lower than with additional neutral beam injection heating. The evolution of the electron density, temperature and pressure profiles are analyzed as well. The core density can be up to twice the Greenwald density and keeps increasing, while the value at the separatrix remains essentially constant and slightly below n G . The edge temperature gradient collapses to near zero and NT plasmas are shown to be resilient to such profiles in terms of disruptivity. We also present the time evolution of the inverse electron pressure scale length with the value at the last closed flux surface (LCFS) decreasing below the value at the normalized radius 0.9 near the density limit, demonstrating the clear drop of confinement starting from the edge. This inverse scale length ‘collapse’ at the LCFS also defines well the characteristic behavior of the kinetic profiles approaching a density limit.

density limit↗

The ultraviolet upturn in field luminous red galaxies at 0.3 < z < 0.7

ABSTRACT We derive the evolution of the ultraviolet upturn colour from a sample of field luminous red galaxies at 0.3 < z < 0.7 with −24 < Mr < −21.5. No individual objects are securely detected, so we stack several hundred galaxies within absolute magnitude and redshift intervals. We find that the colour of the ultraviolet upturn (in observed NUV − i which is approximately equivalent to the classical FUV − V at the redshifts of our targets) does not change strongly with redshift to z = 0.7. This behaviour is similar to that observed in cluster ellipticals over this same mass range and at similar redshifts, and we speculate that the processes involved in the origin of the UV upturn are the same. The observations are most consistent with spectral synthesis models containing a fraction of a helium rich stellar population with abundances between 37 per cent and 42 per cent, although we cannot formally exclude a contribution due to residual star formation at the $\sim 0.5\, \mathrm{ per\,cent}$ level (however, this appears unlikely for cluster galaxies that are believed to be more quenched). This suggests that the ultraviolet upturn is a primordial characteristic of early-type galaxies at all redshifts and that an unexpected nucleosynthesis channel may lead to nearly complete chemical evolution at early times.

79 ASTRONOMY AND ASTROPHYSICS↗

A Possible Tidal Disruption Event Candidate in the Black Hole Binary System of OJ 287

The BL Lacertae OJ 287 is a supermassive black hole binary system with complex physics of its irregular flares. During 2016 October–2017 April, a surprising outburst in the X-ray, UV, and optical bands was detected, while no variability was seen in the γ-ray light curve. During the outburst, the X-ray light curves were dominated by the soft X-rays, whose peak in luminosity was ~10 46 erg s -1 —more than 10 times higher than the mean level before the outburst—and a "softer-when-brighter" phenomenon was exhibited. These above phenomena have been reported by some previous works. The hardness ratio showed negligible evolution with flare time and soft X-ray luminosity. Critically, the luminosity of the soft X-rays decayed following a power law of $t$ -5/3 , which occurs in most tidal disruption events (TDEs), and a similar trend can be seen in the UV and optical bands during the soft X-ray declining period. Helium and oxygen narrow emission lines are strengthened prominently in the optical spectra of post outburst epochs, that could be attributed to the surrounding gas caused by TDE. In conclusion, we discuss three possible origins of the event, including the jet's precession, the aftereffects of the black hole–disk impaction, and the TDE. Our results show that the TDE is the more likely scenario to explain the outburst.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Spectral Evolution of Hot White Dwarf Stars. II. Time-dependent Simulations of Element Transport in Evolving White Dwarfs with STELUM

White dwarf stars are subject to various element transport mechanisms that can cause their surface composition to change radically as they cool, a phenomenon known as spectral evolution. In this paper, we undertake a comprehensive theoretical investigation of the spectral evolution of white dwarfs. First, we introduce STELUM, a new implementation of the stellar evolutionary code developed at the Université de Montréal. We provide a thorough description of the physical content and numerical techniques of the code, covering the treatment of both stellar evolution and chemical transport. Then, we present two state-of-the-art numerical simulations of element transport in evolving white dwarfs. Atomic diffusion, convective mixing, and mass loss are considered simultaneously as time-dependent diffusive processes and are fully coupled to the cooling. We first model the PG 1159-DO-DB-DQ evolutionary channel: a helium-, carbon-, and oxygen-rich PG 1159 star transforms into a pure-helium DB white dwarf due to gravitational settling and then into a helium-dominated, carbon-polluted DQ white dwarf through convective dredge-up. We also compute for the first time the full DO-DA-DC evolutionary channel: a helium-rich DO white dwarf harboring residual hydrogen becomes a pure-hydrogen DA star through the float-up process and then a helium-dominated, hydrogen-bearing DC star due to convective mixing. We demonstrate that our results are in excellent agreement with available empirical constraints. In particular, our DO-DA-DC simulation perfectly reproduces the lower branch of the bifurcation observed in the Gaia color–magnitude diagram, which can therefore be interpreted as a signature of spectral evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Attosecond time delay in atomic photoionization: Angular-dependent transition from dipole to quadrupole and spin-flip dynamics

The angular distribution of attosecond Wigner time delay has been investigated including both quadrupole and relativistic effects using the relativistic-random-phase approximation that is based on the Dirac equation and includes significant aspects of many-body correlations. Here the results show a dramatic evolution of the time delay from (essentially nonrelativistic) dipole to quadrupole and spin-flip dynamics as a function of angle, thereby providing a venue for studying these interactions at the attosecond level.

74 ATOMIC AND MOLECULAR PHYSICS↗

Arctic Sea Ice in CMIP6

Abstract We examine CMIP6 simulations of Arctic sea‐ice area and volume. We find that CMIP6 models produce a wide spread of mean Arctic sea‐ice area, capturing the observational estimate within the multimodel ensemble spread. The CMIP6 multimodel ensemble mean provides a more realistic estimate of the sensitivity of September Arctic sea‐ice area to a given amount of anthropogenic CO 2 emissions and to a given amount of global warming, compared with earlier CMIP experiments. Still, most CMIP6 models fail to simulate at the same time a plausible evolution of sea‐ice area and of global mean surface temperature. In the vast majority of the available CMIP6 simulations, the Arctic Ocean becomes practically sea‐ice free (sea‐ice area <1 × 10 6 km 2 ) in September for the first time before the Year 2050 in each of the four emission scenarios SSP1‐1.9, SSP1‐2.6, SSP2‐4.5, and SSP5‐8.5 examined here.

Notz, Dirk↗