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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Inference of Induced Fracture Geometries Using Fiber-Optic Distributed Strain Sensing in Hydraulic Fracture Test Site 2

Fiber-optic Distributed Strain Sensor observations in a vertical well allow rapid inference of the time evolution of nearby hydraulic fracture heights and their vertical aperture distributions during stimulation and depressurization. We describe a computationally efficient deconvolutional model based on an elastostatic kernel for fracture opening within an isotropic, homogeneous, linear elastic medium. Furthermore, an application of this model to observations from the Hydraulic Fracture Test Site 2 reveals that ultimate fracture heights are typically greater than those indicated by microseismic monitoring, and some hydraulic fractures continue their upward growth for several hours after high-pressure injections ceased.

58 GEOSCIENCES↗
DOE OSTI2020

Atomic electron tomography in three and four dimensions

Atomic electron tomography (AET) has become a powerful tool for atomic-scale structural characterization in three and four dimensions. It provides the ability to correlate structures and properties of materials at the single-atom level. With recent advances in data acquisition methods, iterative three-dimensional (3D) reconstruction algorithms, and post-processing methods, AET can now determine 3D atomic coordinates and chemical species with sub-Angstrom precision, and reveal their atomic-scale time evolution during dynamical processes. In this work, we review the recent experimental and algorithmic developments of AET and highlight several groundbreaking experiments, which include pinpointing the 3D atom positions and chemical order/disorder in technologically relevant materials and capturing how atoms rearrange during early nucleation at four-dimensional atomic resolution.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Simulating Magnetized Electron Cooling for EIC with JSPEC

We present a possible electron cooling configuration for the proposed Electron Ion Collider (EIC) facility, developed using a Nelder-Mead Simplex optimization procedure built into JSPEC, an electron cooling code developed at Jefferson Lab. We show the time evolution of the emittance of the ion beam in the presence of this cooler evaluated assuming the ion distribution remains Gaussian. We also show that bi-gaussian distributions emerge in simulations of ion macro-particles. We show how intra-beam scattering can be treated with a core-tail model in simulations of ion macro-particles. The Sirepo/JSPEC* and Sirepo/Jupyter** apps will be presented, with instructions enabling the community to reproduce our simulations.

Coleman, S. J.↗
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 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 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↗