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

In situ probing of interfacial roughness and transient phases during ceramic cold sintering process

The ceramic cold sintering process (CSP) offers an eco-friendly approach to producing fully dense ceramics at low temperatures. However, an incomplete mechanistic understanding hinders its optimization and widespread adoption. In this study, we analyze the microstructural and structural changes in ZnO, a model CSP system, using in situ synchrotron-based high-energy small-angle X-ray scattering and X-ray diffraction techniques. Our results reveal the time evolution of ZnO particles' surface area and roughness, reflecting the dissolution and re precipitation processes that enable densification. The in situ measurements supply valuable kinetic data for these stages of CSP. Alongside microstructural changes and densification, we observed the evolution of secondary phases representing reaction products between ZnO and acetic acid, the solvent used. The initial ZnO/solvent mixture's dominant secondary phase is attributed to zinc acetate, which is gradually replaced by a zinc soap-type structure during CSP. This structure has a large (≈ 21 Å) lattice parameter and is assumed to have a layered nature. The formation of this soap phase, which is retained in the sintered product as an intergranular component, appears to be a signature of successful cold sintering as it facilitates mass transport, leading to densification. Here, our study underscores the potential of in situ synchrotron characterization for revealing microstructural and phase-evolution details during CSP. These findings, which would be challenging to obtain through ex situ measurements, provide crucial data to guide and validate theoretical models, ultimately enhancing CSP's effectiveness and adoption.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Temperature-dependent dynamic fouling on superhydrophobic and slippery nonwetting copper surfaces

Bioinspired, superhydrophobic and slippery liquid infused surfaces that offer nonwetting characteristics have been explored in recent years for fouling mitigation. However, most of the studies are in the context of biofouling or under static immersion at ambient temperature conditions that are not reflective of the dynamic flow environment in practice. This article presents, for the first time, a systematic study of dynamic fouling of superhydrophobic (SHS) and slippery lubricant-infused surfaces (LIS) over a range of flow and temperature conditions. In this work, copper metallic surfaces were textured via electrodeposition or etching and further functionalized to achieve SHS and, additionally, infiltrated with a lubricant to fabricate LIS. The nonwetting surfaces were studied for their fouling behavior in a rotating Couette flow of a supersaturated calcium sulfate solution at different rotational speed and temperature. Fouling mineral mass accumulation on the different surfaces was measured as a function of time over a period of days using inductively coupled plasma mass spectroscopy and the fouled surfaces were investigated using scanning electron microscopy. Both SHS and LIS showed superior anti-scaling performance at all ranges of variables. An analytical Hill-Langmuir model is presented, for the first time, to describe the time evolution of scaling within 20% accuracy over the range of parameters studied. The study is the first to juxtapose two surface texturing methods, electrodeposition and etching, and two nonwetting surface types, SHS and LIS, subject to a common suite of experiments to elucidate fundamental understanding of mineral fouling on nonwetting surfaces.

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

Modeling gas release from a Bingham plastic slurry and deconvoluting measured data

Physics-based models for predicting the off-gassing characteristics of agitated slurries are presented. This approach decomposes the system into two separate but coupled slurry and headspace models. The physics driving bubble transport through the slurry and gas mixing within the headspace are discussed. An analytical expression for predicting the time evolution of the headspace concentration using first principles theory is also presented. Predictions from the numerical models, as well as expectations from the analytical solution, both agree with measured off gassing data for two different experimental operating conditions. After benchmarking the numerical model predictions against experimental data, the model was used to make predictions for gas release from a full-scale vessel and to perform sensitivity analyses to examine the sensitivity of gas release to parameters such as yield stress, consistency index, slurry density, bubble size, bubble concentration, and impeller speed.

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

Explainable artificial intelligence relates perovskite luminescence images to current-voltage metrics

As the demand for low-cost, high-efficiency solar energy technologies grows, metal halide perovskite (MHP) solar cells have emerged as a promising candidate for next-generation photovoltaics due to their high power conversion efficiencies. However, their poor durability and issues with manufacturing consistency remain significant barriers to commercialization. In this work, we develop deep learning models to support materials characterization and provide insight into features and processes influencing performance. The models are trained using transfer learning of a pretrained model to predict relevant current-voltage (IV) metrics based on different combinations of input electroluminescence (EL) and photoluminescence (PL) images of MHP devices. We examine which image types are most informative in accurately predicting different IV metrics. Additionally, we use explainable artificial intelligence (XAI) techniques to provide insights into specific spatial features in the devices that drive differences in performance. We find that stabilized luminescence images (e.g. those collected after biasing the devices for at least 1 min) are better for predicting metrics of open-circuit voltage (by PL) and short-circuit current (by PL with EL), but that predicting fill factor and overall power output may use the time-evolution of EL images. Based on attribution masks generated by integrated gradients for each device performance metric, we further suggest different loss mechanisms associated with categories of large and small spatial defects. Overall, this case study highlights the potential applicability of XAI methodology for streamlining MHP device analysis and accelerating detailed understanding of the relationships between spatial defects and impacts on performance.

14 SOLAR ENERGY↗
DOE OSTI2025

A Comprehensive Chemistry Evaluation and Diagnostics Package for E3SM – ChemDyg Version 1.1.0

The Chemistry Evaluation and Diagnostics Package (ChemDyg) is an open-source tool designed for the Energy Exascale Earth System Model (E3SM) developed by the U.S. Department of Energy. ChemDyg facilitates routine evaluation, tailored development, and in-depth analysis of atmospheric chemistry through its modular architecture, allowing users to compare model outputs with observational data. Version 1.1.0 introduces a robust set of diagnostic capabilities, including climatology, time evolution of key tracers, diurnal and annual cycle analyses, and extensive budget diagnostics. These features help identify model discrepancies and enhance the representation of atmospheric chemistry in E3SM. Each self-contained diagnostic set includes dedicated scripts and documentation for ease of use. The interactive HTML output improves data accessibility, accelerating chemistry model development. Additionally, ChemDyg's flexible framework allows for customization, enabling users to create unique diagnostic sets for specific scientific contributions.

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

Development of spatially and temporally resolved electron density measurements for the assessment of radiation hydrodynamics simulations of laboratory X-ray photoionized plasmas

The photoionized plasma gas cell experiment is an established platform we use to make at-parameter (ξ > > 1 ergs cm s –1 ) measurements of plasma properties with application to high-energy astrophysical systems. We model the experiments with 1D radiation hydrodynamics simulations using the HELIOS-CR code to inform our understanding and assist in the interpretation of results. The simulations predict that the bulk of the plasma is in a quasi-uniform and hydrodynamically unperturbed state throughout the duration of the experiment. To evaluate this prediction, we introduced a photonic Doppler velocimetry (PDV) diagnostic to measure spatially and temporally resolved plasma electron density. The initial measurements were successful but had limitations that made model-data comparisons challenging. To address this, we re-designed the gas cell PDV diagnostic and doubled the number of measurement locations to sample across two thirds of the depth of the cell. We also present a comparison of the results from the upgraded PDV diagnostic to the HELIOS-CR simulations for the first time. As a result, the experimental data confirms the prediction of an unperturbed region in the bulk of the plasma but reveals discrepancies in the time evolution and spatial distribution of the simulated electron density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

High and Ultra-High Temperature Reaction Kinetics by Single Nanoparticle Mass Spectrometry

Methodology is presented for non-destructive, optically-detected single nanoparticle (NP) mass spectrometry, with the goal of extracting surface reaction kinetics for single NPs at high temperatures. Methods for determining the NP charge, mass, and temperature as a function of time are discussed, and the data are used to extract both the absolute kinetics for mass change, as well as the efficiencies of the surface processes that cause them. Factors that contribute to the uncertainties in absolute and relative mass determination, and in the resulting kinetic parameters, are discussed. The method allows the NP-to-NP variations in initial reactivity to be measured directly, along with the time evolution of reactivity resulting from NP structural/compositional changes that occur under reaction conditions. The strengths and limitations of single nanoparticle mass spectrometry as a high temperature surface kinetics tool are discussed in the context of sublimation and O2 oxidation kinetics experiments for single hafnium (Hf) NPs at temperatures ranging above 2400 K. The Hf oxidation kinetics are compared to analogous oxidation experiments for silicon, graphite, and carbon black NPs. In all four cases, the oxidation chemistry was dominated by processes that result in net mass loss, and the distinct mechanisms responsible are discussed. All four NPs also eventually passivated, i.e., the efficiencies for oxidative etching decreased by at least two orders of magnitude, relative to the initial efficiencies. Furthermore, the passivation mechanisms, which are quite different for carbon, compared to silicon or hafnium, are discussed. Carbon NP passivation is attributed to structural isomerization leading to fully coordinated, fullerene-like NP surfaces, while for silicon and hafnium, passivation results from delayed formation of an oxide layer, triggered by accumulation of oxygen in the NP sub-surface region.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Maximum bound principle preserving integrating factor Runge–Kutta methods for semilinear parabolic equations

A large class of semilinear parabolic equations satisfy the maximum bound principle (MBP) in the sense that the time-dependent solution preserves for any time a uniform pointwise bound imposed by its initial and boundary conditions. Here, the MBP plays a crucial role in understanding the physical meaning and the well-posedness of the mathematical model. Investigation on numerical algorithms with preservation of the MBP has attracted increasingly attentions in recent years, especially for the temporal discretizations, since the violation of MBP may lead to nonphysical solutions or even blow-ups of the algorithms. In this paper, we study high-order MBP-preserving time integration schemes by means of the integrating factor Runge-Kutta (IFRK) method. Beginning with the space-discrete system of semilinear parabolic equations, we present the IFRK method in general form and derive the sufficient conditions for the method to preserve the MBP. In particular, we show that the classic four-stage, fourth-order IFRK scheme is MBP preserving for some typical semilinear systems although not strong stability preserving, which can be instantly applied to the Allen-Cahn type of equations. To our best knowledge, this is the first time to present a fourth-order linear numerical method preserving the MBP. In addition, convergence of these numerical schemes is proved theoretically and verified numerically, as well as their efficiency by simulations of 2D and 3D long-time evolutional behaviors. Numerical experiments are also carried out for a model which is not a typical gradient flow as the Allen-Cahn type of equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Cluster dynamics modeling of irradiation growth in single crystal Zr

Irradiation growth of single-crystal Zr is modeled with a reduced set of cluster dynamics equations. Nucleation and growth of basal and prismatic dislocation loops are both accounted for in the model. Equations are developed for the time evolution of single point defects, a limited number of point defect clusters, interstitial and vacancy loop nucleation rates, as well as for the growth of vacancy loops on basal planes and interstitial loops on prismatic planes. Reduction of the usual infinite hierarchy of cluster dynamics equations to the simple set studied here is justified on the physical basis of the stability of small loops once nucleated. This simplified cluster dynamics model with a small number of adjustable parameters avoids the complexity of explicit representation of higher order point defect clusters. The model shows consistency with experimental observations of the following aspects: (1) the growth rates of Zr crystals along the a- and c-axes; (2) the onset dose for breakaway irradiation growth; (3) the saturation dislocation loop densities of vacancy and interstitial loops; (4) the effects of cold work and temperature on irradiation growth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Do atoms age?

Time evolution generically entangles a quantum state with environmental degrees of freedom. The resulting increase in entropy changes the properties of that quantum system leading to “aging”. It is interesting to ask if this familiar property also applies to simple, single particle quantum systems such as the decay of a radioactive particle. We propose a test of such aging in an ion clock setup where we probe for temporal changes to the energies of the electronic state of an ion containing a radioactive nucleus. Such effects are absent in standard quantum mechanics and this test is thus a potent null test for violations of quantum mechanics. As a proof of principle, we show that these effects exist in causal non-linear modifications of quantum mechanics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Stopping to reflect: Asymptotic static moving mirrors as quantum analogs of classical radiation

Radiation from an accelerating charge is a basic process that can serve as an intersection between classical and quantum physics. We present two exactly soluble electron trajectories that permit analysis of the radiation emitted, exploring its time evolution and spectrum by analogy with the moving mirror model of the dynamic Casimir effect. These classical solutions are finite energy, rectilinear (nonperiodic), asymptotically zero velocity worldlines with corresponding quantum analog beta Bogolyubov coefficients. One of them has an interesting connection to uniform acceleration and Leonardo da Vinci's water pitcher experiment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum state preparation with resolution refinement

We introduce a method called resolution refinement that allows one to bootstrap eigenstate preparation on a quantum computer. We first prepare an eigenstate of a low-resolution Hamiltonian using any method of choice. The eigenstate is then lifted to higher resolution and adiabatically evolved to produce the corresponding eigenstate of a higher-fidelity Hamiltonian. We give examples of resolution refinement applied to both single-particle basis states as well as a spatial lattice grid. For basis refinement, we compute few-body ground states of the Busch model for interacting particles in a harmonic trap in one dimension. For lattice refinement, we compute Hartree-Fock nuclear states for a central Woods-Saxon potential in three dimensions, and we compute bound states and continuum states in a multi-species Hubbard model of fermions in one dimension. In all cases, the method is efficient and requires an adiabatic evolution time that scales with the inverse of the energy gap times the square root of the system size. We show that this very favorable scaling arises from the fact that resolution refinement does not make large changes to the structure or energies of the low-energy eigenstates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Quasi-periodic travelling gravity–capillary waves

We present a numerical study of spatially quasi-periodic travelling waves on the surface of an ideal fluid of infinite depth. This is a generalization of the classic Wilton ripple problem to the case when the ratio of wavenumbers satisfying the dispersion relation is irrational. We propose a conformal mapping formulation of the water wave equations that employs a quasi-periodic variant of the Hilbert transform to compute the normal velocity of the fluid from its velocity potential on the free surface. We develop a Fourier pseudo-spectral discretization of the travelling water wave equations in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on the torus. This leads to an overdetermined nonlinear least-squares problem that we solve using a variant of the Levenberg–Marquardt method. We investigate various properties of quasi-periodic travelling waves, including Fourier resonances, time evolution in conformal space on the torus, asymmetric wave crests, capillary wave patterns that change from one gravity wave trough to the next without repeating and the dependence of wave speed and surface tension on the amplitude parameters that describe a two-parameter family of waves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2026

Learning Latent Representations to Bridge Coarse-Grained and Atomistic Resolutions in Polymer Simulations

We present a machine-learning-based framework for learning reduced-order representations of polymer chain conformations across coarse-grained (CG) and united-atom (UA) fidelities. By employing linear singular value decomposition and nonlinear autoencoders, we compress high-dimensional polymer configurations into latent spaces with minimal loss of structural accuracy. Crucially, we demonstrate a near-perfect linear mapping between CG and UA latent spaces, enabling an efficient super-resolution back-mapping procedure that reconstructs high-fidelity UA configurations from CG simulations. While minor structural inaccuracies occur, they are effectively corrected through a brief molecular dynamics relaxation, forming a practical hybrid machine learning−physics scheme. This approach establishes the key structural prerequisites for accelerated polymer dynamics simulations: a compact and accurate latent encoding of polymer chain conformations and a validated multi-fidelity mapping that permits reconstruction of UA structures from CG configurations. The extension of this framework to explicit time evolution within the latent space, enabling dynamics to be propagated at CG fidelity and decoded to UA resolution only when required, represents a natural and well-motivated direction for future work.

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

A Multireference Quantum Krylov Algorithm for Strongly Correlated Electrons

We introduce a multireference selected quantum Krylov (MRSQK) algorithm suitable for quantum simulation of many-body problems. MRSQK is a low-cost alternative to the quantum phase estimation algorithm that generates a target state as a linear combination of non-orthogonal Krylov basis states. This basis is constructed from a set of reference states via real-time evolution; thus, avoiding the numerical optimization of parameters. An efficient algorithm for the evaluation of the off-diagonal matrix elements of the overlap and Hamiltonian matrices is discussed and a selection procedure is introduced to identify a basis of orthogonal references that ameliorates the linear dependency problem. Preliminary benchmarks on linear H 6 , H 8 , and BeH 2 indicate that MRSQK can predict the energy of these systems accurately using very compact Krylov bases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗