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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Adaptive variational quantum minimally entangled typical thermal states for finite temperature simulations

Scalable quantum algorithms for the simulation of quantum many-body systems in thermal equilibrium are important for predicting properties of quantum matter at finite temperatures. Here we describe and benchmark a quantum computing version of the minimally entangled typical thermal states (METTS) algorithm for which we adopt an adaptive variational approach to perform the required quantum imaginary time evolution. The algorithm, which we name AVQMETTS, dynamically generates compact and problem-specific quantum circuits, which are suitable for noisy intermediate-scale quantum (NISQ) hardware. We benchmark AVQMETTS on statevector simulators and perform thermal energy calculations of integrable and nonintegrable quantum spin models in one and two dimensions and demonstrate an approximately linear system-size scaling of the circuit complexity. We further map out the finite-temperature phase transition line of the two-dimensional transverse field Ising model. Finally, we study the impact of noise on AVQMETTS calculations using a phenomenological noise model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Activation and residual dose rates in the NSLS-II storage ring

It is well known that components in high energy accelerators may get activated and expose workers during maintenance activities. While the potential for activation is higher in positive ion accelerators due to the prolific production of neutrons, certain regions and components in high energy electron accelerators too can be activated enough to create residual dose rates that require careful attention for ALARA compliance. In this report, the activation and residual dose rates from a few materials are studied from various beam loss scenarios ranging from chronic point losses to intermittent real loss scenarios obtained from the operational parameters of the NSLS-II storage ring. FLUKA Monte Carlo radiation transport code is used to model the geometry, the time profile of the beam losses and score the resulting time evolution of the remnant activity and the residual dose rates.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Deep Conservation: A latent dynamics model for exact satisfaction of physical conservation laws [Report]

This work proposes an approach for latent dynamics learning that exactly enforces physical conservation laws. The method comprises two steps. First, we compute a low-dimensional embedding of the high-dimensional dynamical-system state using deep convolutional autoencoders. This defines a low-dimensional nonlinear manifold on which the state is subsequently enforced to evolve. Second, we define a latent dynamics model that associates with a constrained optimization problem. Specifically, the objective function is defined as the sum of squares of conservation-law violations over control volumes in a finite-volume discretization of the problem; nonlinear equality constraints explicitly enforce conservation over prescribed subdomains of the problem. The resulting dynamics model—which can be considered as a projection-based reduced-order model—ensures that the time-evolution of the latent state exactly satisfies conservation laws over the prescribed subdomains. In contrast to existing methods for latent dynamics learning, this is the only method that both employs a nonlinear embedding and computes dynamics for the latent state that guarantee the satisfaction of prescribed physical properties. Numerical experiments on a benchmark advection problem illustrate the method's ability to significantly reduce the dimensionality while enforcing physical conservation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

QIS for Applied Quantum Field Theory

Precision predictions in HEP and non-equilibrium dynamics in the early universe require quantum computational technologies. Future quantum devices offer the possibility of performing these, complementing lattice quantum chromodynamics (QCD) studies on classical computers. We are investigating many aspects of the quantum simulation of quantum field theories. Specific topics include simulation-suitable formulations of gauge theories, state preparation, time evolution, extracting HEP-related observables from simulations, quantum error correction (QEC), and entanglement. We are exploring strategies for simulation of scattering in QFTs, and dynamics of the early Universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Quantum Quench Dynamics–Crossover Phenomena in Nonequilibrium Correlated Quantum Systems (Final Technical Report)

In this project we have developed analytical and numerical tools to study the interplay between interactions, dimensionality, fluctuations and disorder in quantum impurity systems. We have studied the transient dynamics following a quench, i.e. the sudden change in a control variable of the model Hamiltonian, which offers a unique way to examine the system’s spatial, temporal and energetic structure. Specifically, we have investigated how crossover scales, i.e. given by Kondo couplings, are affected by disorder. We have addressed a wide range of phenomena, such as the effects of static disorder on dynamical screening of magnetic impurities, and the influence of a fully time-fluctuating bath environment on the dephasing of system qubits. We have found that many of these systems exhibit fascinating counterintuitive behavior, such as e.g. the possibility of asymptotic unitary time evolution in specifically conditioned system-bath Hamiltonians. From a technical point of view, we have developed a supersymmetric formalism to analyze scaling functions of disordered quantum impurity spin chains, we generalized the strong-coupling real-space renormalization group method to probe the quench dynamics of higher-dimensional quantum impurity systems, and we devised a master equation approach that included the effects of a fully fluctuating environment. This way, we have addressed a very broad spectrum of transient phenomena in quantum impurity systems, with disorder ranging from the quenched to the annealed limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

CMDV (CM)4 Project - University of Washington contribution. Final report

The original goal of this project was to diagnose and improve CLUBB, the turbulence and cloud fraction parameterization used in DOE’s E3SM model. For this purpose, we planned to use large-eddy simulation (LES) and ARM observations from the NE Pacific Ocean and the SGP site, in coordination with ARM’s LASSO program. It was determined that for these cloud regimes, many of the turbulent ‘moment closures’ which underlie CLUBB’s mathematical formulation are inconsistent with LES, which is an appropriate benchmark for testing this. The research assistant found that these closures could be more accurately formulated using machine learning using the LES as a training dataset. However the resulting parameterization proved to quickly drift away from physical plausibility. A novel machine-learning based boundary layer parameterization called MARBLE based on matching the time evolution of a parameterized cloud-topped boundary layer to reanalysis was then developed and published.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Table-Top High Photon Energy Sources for Chemical Dynamics Investigations

Time-resolved spectroscopies using high-energy photons in the vacuum ultraviolet (VUV) to the X-ray region of the electromagnetic spectrum, have proven to be powerful probes of chemical dynamics. These high-energy photons can access valence and core orbitals of molecules and materials, providing key information on molecular and electronic structure and their time evolution. This report details the development of table-top sources of extreme ultraviolet (XUV) and VUV pulses at Sandia National Laboratories for use in studies of gas phase chemical dynamics. Femtosecond duration XUV pulses are produced using laser-driven high harmonic generation and their detected range span ~40-140 eV photon energies. These pulses are used in conjunction with ultraviolet pulses in a pump-probe scheme to study excited state dynamics of gas phase molecules. VUV pulses at 7.75 eV are generated using a four-wave-mixing scheme driven by 800 nm and 266 nm pulses in an argon-filled hollow-core fiber.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Data-Driven Specification of Earth's Magnetic Field

As our society is increasingly dependent on technologies in space, predicting space weather becomes ever more important. We have come to rely on technologies such as cellphones, GPS, the Internet, and other commercial and military assets. In order to protect these systems a high priority is to develop better understanding of the harmful conditions of space weather encountered in the near-Earth radiation environment which are mostly dictated by the dynamics of energetic charged particles. The main objective of our proposal was to improve the specification and forecast of the particle distribution in the Earth’s inner magnetosphere using data assimilation together with machine learning (ML) algorithms in our ring current-atmosphere interactions model with self-consistent magnetic and electric fields (RAM-SCBE). Specifically, we explored the application of machine learning algorithms to emulate and improve, through infusion of observational data, the magnetic field (SCBE) of our ring current model. Our work advanced towards the characterization of the RAM-SCBE particle flux field and its potentially different simulation scenarios. We explored how the model behaves in a number of simulation scenarios in order to better understand how the model reacts to perturbations in the solar wind. We were able to capture this behavior through the simulation of an ensemble of RAM-SCB model solutions. The ensemble showed the important correlations between the particle fluxes and the magnetic field. We were able to start the training of a neural network (Autoencoder) to identify the dominant features of the magnetic field, and the main correlations (linear and non-linear) between the particle fluxes and the magnetic field. Our next steps consists of training a neural network to capture the time evolution of the dominant features and obtain a consistent surrogate model for the RAM-SCB.

42 ENGINEERING↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Characterize the Rings in Protoplanetary Disks Formed by an Eccentric Instability [Slides]

Typical long-term time evolution of gas surface density in a disk undergoing eccentric mode instability (EMI). In a disk with no initial substructures or external perturbers, EMI generates a growing one arm spiral (top right). After EMI saturates, the disk continues to evolve with a long-lived coherent eccentric pattern. If the disk has strong wave damping, the disk is then left with circular rings. Synthesized dust continuum images of some EMI simulations from this IC project show some characteristics similar to the actual PPD images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

AMPPED LDRD: Advanced Models of the Physics & Phenomena of Electrical Discharge.

This Laboratory Directed Research and Development (LDRD) project aimed at addressing an unmet need for physics-based models that enable prediction of electrical breakdown in the presence of high voltage gradients and real-world environments. This project focused on understanding coupled plasma-surface feedback phenomena critical to electrical discharge processes and design of breakdown voltages, breakdown times, and breakdown paths. New secondary electron emission and coupled plasma-surface measurement methods and models were developed. Highly reproducible discharge experiments were developed to provide current-voltage data as a function of time, and methods were developed to condition electrodes for predictable discharge. To assess plasma-surface feedback mechanisms, new surface physics capabilities, including measurement systems for local and global photoemission yield and ion-induced secondary electron emission yield were developed. Discharge models were developed to incorporate these yields for predicting discharge behavior to understand influence of surface physics on discharge current and time evolution, and methods to accelerate these simulations were demonstrated. A key finding was that the initiation process of a non-vacuum discharge appears highly dependent on photoemission at the cathode. Advanced diagnostics were developed, including spatial electrical field imaging, spectroscopy of plasmas, and high-speed imaging of discharge event s in metal and dielectric environments. These diagnostics and improved electrical models are a step toward enabling quantitative design margins, assurance of safety architectures, and implementation of novel materials and technologies for wide ranging energy and electronics applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Assessing the Influence of Background State and Climate Variability on Tropical Cyclones Using Initialized Ensembles and Mesh Refinement in E3SM

The scientific goal of this project is to assess the influence of the background state and horizontal resolution of climate models on the simulation of extreme events like intense precipitation, tropical cyclones, atmospheric rivers, and heatwaves. The main results are: - Initialized ensembles of 4-week forecasts using a climate model (E3SM) show that dynamic (flow) and thermodynamic (moisture) biases asymptote to their climatological values at different rates—dynamic bias grows rapidly and approaches its climatological value within about 2 weeks whereas the thermodynamic bias takes about 4 weeks. - Moisture bias has a bigger impact on extreme events like tropical cyclones and intense precipitation than flow bias, and this impact is non-monotonic over bias evolution time. - Climate models are capable of simulating the flow features associated with unprecedented extreme events like the Western North America heatwave of June 2021, but moisture biases can weaken the surface manifestation of extreme heatwaves. Reducing model bias may be more important than increasing model resolution in improving the fidelity of heatwave simulations in climate models. - Climate models are capable of simulating weather patterns associated with “power droughts”, i.e., periods with low solar and wind power generation, provided model biases in the simulation of the background state can be corrected. - Regional mesh refinement can improve the simulation of extreme weather events but model biases may still persist.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Analysis of Neural Networks as Random Dynamical Systems

In this report we present our findings and outcomes of the NNRDS (analysis of Neural Networks as Random Dynamical Systems) project. The work is largely motivated by the analogy of a large class of neural networks (NNs) with a discretized ordinary differential equation (ODE) schemes. Namely, residual NNs, or ResNets, can be viewed as a discretization of neural ODEs (NODEs) where the NN depth plays the role of the time evolution. We employ several legacy tools from ODE theory, such as stiffness, nonlocality, autonomicity, to enable regularization of ResNets thus improving their generalization capabilities. Furthermore, armed with NN analysis tools borrowed from the ODE theory, we are able to efficiently augment NN predictions with uncertainty overcoming wellknown dimensionality challenges and adding a degree of trust towards NN predictions. Finally, we have developed a Python library QUiNN (Quantification of Uncertainties in Neural Networks) that incorporates improved-architecture ResNets, besides classical feed-forward NNs, and contains wrappers to PyTorch NN models enabling several major classes of uncertainty quantification methods for NNs. Besides synthetic problems, we demonstrate the methods on datasets from climate modeling and materials science.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

NuSTAR broad-band X-ray observation campaign of energetic pulsar wind nebulae in synergy with VERITAS, HAWC and Fermi gamma-ray telescopes

We report recent progress on the on-going NuSTAR observation campaign of 8 TeV-detected middle-aged pulsar wind nebulae (PWNe). This campaign constitutes a major part of our NuSTAR study of some of the most energetic TeV sources in our Galaxy detected by VERITAS and HAWC. NuSTAR is the first focusing X-ray telescope operating above 10 keV in space with sub-arcminute angular resolution. Broad-band X-ray imaging and spectroscopy data, obtained by NuSTAR, allow us to probe sub-PeV electron populations through detecting synchrotron X-ray radiation. Our targets include PeVatron candidates detected by HAWC, the Boomerang nebula, PWNe crushed by supernova remnant shocks (or relic PWNe) and G0.9‚0.1 in the Galactic Center. Combined with Fermi-LAT data and available TeV data, we aim to provide a complete, multi-wavelength view of a diverse class of middle-aged (~10–100 kyrs old) PWNe. Our NuSTAR analysis detected hard X-ray emission from the Eel and Boomerang PWNe and characterized their broad-band X-ray spectra most accurately. We plan to apply both time-evolution and multi-zone PWNe models to multi-wavelength spectral energy distribution (SED) data over the radio, X-ray, GeV and TeV bands. In this proceeding, we will review our observation campaign and discuss the preliminary results for some PWNe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Three-pion effects in $K^0-\bar{K}^0$ mixing

The rate of mixing between a neutral kaon and an anti-kaon ($K^0$-$\bar{K}^0$) is given, in part, by a long-range matrix element, defined with two insertions of the weak Hamiltonian separated by physical, Minkowski time evolution. For physical quark masses, the kaon mass lies above the two- and three-pion thresholds and, as a result, this long-range matrix element receives contributions from intermediate on-shell $2\pi$ and $3\pi$ states. These contributions cannot easily be captured in a finite Euclidean spacetime, meaning that such matrix elements are not directly accessible via lattice QCD. In this talk, we present a strategy for combining quantities that can be extracted in numerical lattice QCD calculations in order to reproduce the physical, infinite-volume long-range amplitude for $K^0-\bar{K}^0$. The key novelty relative to published work is that we fully include the effects of three-particle states that were previously neglected. The strategy is built on existing formalism for long-range matrix elements with two-particle intermediate states, together with the relativistic-field-theory finite-volume formalism for extracting three-hadron weak decays.

Jackura, Andrew↗
DOE OSTI2023

Improved Fermion Hamiltonians for Quantum Simulation

We developed a Hamiltonian inspired by ASQTAD and highly improved staggered quark (HISQ) actions and show how these Hamiltonians can be used for quantum simulations. Gate costs for the time evolution of these improved Hamiltonians are provided as well as a demonstration of the reduction of lattice spacing errors using the 1+1d lattice Schwinger model.

Gustafson, Erik↗
DOE OSTI2021

Quantum Algorithm for Simulating Hamiltonian Dynamics with an Off-diagonal Series Expansion

We propose an efficient quantum algorithm for simulating the dynamics of general Hamiltonian systems. Our technique is based on a power series expansion of the time-evolution operator in its off-diagonal terms. The expansion decouples the dynamics due to the diagonal component of the Hamiltonian from the dynamics generated by its off-diagonal part, which we encode using the linear combination of unitaries technique. Our method has an optimal dependence on the desired precision and, as we illustrate, generally requires considerably fewer resources than the current state-of-the-art. We provide an analysis of resource costs for several sample models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗