Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative systems”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The coset construction for non-equilibrium systems

We propose a systematic coset construction of non-equilibrium effective field theories (EFTs) governing the long-distance and late-time dynamics of relativistic, finite-temperature condensed matter systems. Our non-equilibrium coset construction makes significant advances beyond more standard coset constructions in that it takes advantage of recently-developed techniques, which allow the formulation of non-equilibrium effective actions that account for quantum and thermal fluctuations as well as dissipation. Because these systems exist at finite temperature, the EFTs live on the closed-time-path of the Schwinger-Keldysh contour. Since the coset construction and the non-equilibrium effective actions may be unfamiliar to many readers, we include brief introductions to these topics in an effort to make this paper self-contained. To demonstrate the legitimacy of this coset construction, we successfully reproduce the known EFTs for fluids and superfluids at finite temperature. Then, to demonstrate its utility, we construct novel EFTs for solids, supersolids, and four phases of liquid crystals, all at finite temperature. We thereby combine the non-equilibrium effective action and the coset construction to create a powerful tool that can be used to study many-body systems out of thermal equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Data-driven reduced-order models for port-Hamiltonian systems with operator inference

Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Colloidal superionic conductors

Nanoparticles with highly asymmetric sizes and charges that self-assemble into crystals via electrostatics may exhibit behaviors reminiscent of those of metals or superionic materials. Here, we use coarse-grained molecular simulations with underdamped Langevin dynamics to explore how a binary charged colloidal crystal reacts to an external electric field. As the field strength increases, we find transitions from insulator (ionic state), to superionic (conductive state), to laning, to complete melting (liquid state). In the superionic state, the resistivity decreases with increasing temperature, which is contrary to metals, yet the increment decreases as the electric field becomes stronger. Additionally, we verify that the dissipation of the system and the fluctuation of charge currents obey recently developed thermodynamic uncertainty relation. Our results describe charge transport mechanisms in colloidal superionic conductors.

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

Detailed Structural and Chemical Analysis of Amorphous Compounds in Superconducting Qubit Systems

With massive improvements in device coherence times and gate fidelity over the past two decades, superconducting quantum devices have emerged as a leading technology platform for next generation quantum computing. While many of these improvements have been driven through optimized device designs and geometries, the presence of defects, impurities, interfaces, and surfaces in the constituent materials continue to limit performance and present a critical barrier in achieving scalable quantum systems with long coherence times. As part of the Superconducting Quantum Materials and Systems (SQMS) center, we have applied state-of-the-art materials characterization techniques, such as scanning/transmission electron microscopy, secondary ion mass spectrometry, and atom probe tomography, to identify structural defects and chemical inhomogeneities in superconducting qubits. In this review, we discuss the findings of these efforts as they relate to the inherent surfaces, interfaces, impurities, and defects in superconducting qubits. Through this comprehensive investigation, we provide insight into possible sources of dissipation in these systems. Further, we supplement these findings with potential strategies for mitigating these undesired features in an effort to systematically and intelligently improve the coherence times of superconducting qubits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
DOE OSTI2025

Quantum simulation of Lindbladian dynamics via repeated interactions

The Lindblad equation generalizes the Schrödinger equation to quantum systems that undergo dissipative dynamics. The quantum simulation of Lindbladian dynamics is therefore non-unitary, preventing a naive application of state-of-the-art quantum algorithms. Here, we make use of an approximate correspondence between Lindbladian dynamics and evolution based on repeated interaction (RI) CPTP maps to write down a Hamiltonian formulation of the Lindblad dynamics and derive a rigorous error bound on the master equation. Specifically, we show that the number of interactions needed to simulate the Liouvillian within error e scales in most physical scenarios as . This is significant because the error in the Lindbladian approximation to the dynamics is not explicitly bounded in existing quantum algorithms for open system simulations. We then provide quantum algorithms to simulate RI maps using an iterative qubitization approach and Trotter–Suzuki formulas, and specifically show that for iterative qubitization the number of operations needed to simulate the dynamics (for a fixed value of ?) scales as in the limit where a0 (the coefficient 1-norm for the system and bath Hamiltonians) asymptotically dominates over the corresponding factor for the interaction Hamiltonian, which is often the case in weak coupling. This scaling would appear to be optimal if the complexity of ? is not considered, which underscores the importance of considering the error in the Liouvillian that we reveal in this work.

Quantum Computing↗
DOE OSTI2020

Stochastic equation of motion approach to fermionic dissipative dynamics. I. Formalism

In this work, we establish formally exact stochastic equation of motion (SEOM) theory to describe the dissipative dynamics of fermionic open systems. The construction of the SEOM is based on a stochastic decoupling of the dissipative interaction between the system and fermionic environment, and the influence of environmental fluctuations on the reduced system dynamics is characterized by stochastic Grassmann fields. Meanwhile, numerical realization of the time-dependent Grassmann fields has remained a long-standing challenge. To solve this problem, we propose a minimal auxiliary space (MAS) mapping scheme with which the stochastic Grassmann fields are represented by conventional c-number fields along with a set of pseudo-levels. This eventually leads to a numerically feasible MAS-SEOM method. The important properties of the MAS-SEOM are analyzed by making connection to the well-established time-dependent perturbation theory and the hierarchical equations of motion theory. The MAS-SEOM method provides a potentially promising approach for the accurate and efficient simulation of fermionic open systems at ultra-low temperatures.

Han, Lu (ORCID:0000000339002225)↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
DOE OSTI2024

Transient anisotropic kernel for probabilistic learning on manifolds

PLoM (Probabilistic Learning on Manifolds) is a method introduced in 2016 for handling small training datasets by projecting an Itô equation from a stochastic dissipative Hamiltonian dynamical system, acting as the MCMC generator, for which the KDE-estimated probability measure with the training dataset is the invariant measure. PLoM performs a projection on a reduced-order vector basis related to the training dataset, using the diffusion maps (DMAPS) basis constructed with a time-independent isotropic kernel. In this paper, we propose a new ISDE projection vector basis built from a transient anisotropic kernel, providing an alternative to the DMAPS basis to improve statistical surrogates for stochastic manifolds with heterogeneous data. The construction ensures that for times near the initial time, the DMAPS basis coincides with the transient basis. For larger times, the differences between the two bases are characterized by the angle of their spanned vector subspaces. The optimal instant yielding the optimal transient basis is determined using an estimation of mutual information from Information Theory, which is normalized by the entropy estimation to account for the effects of the number of realizations used in the estimations. Consequently, this new vector basis better represents statistical dependencies in the learned probability measure for any dimension. Three applications with varying levels of statistical complexity and data heterogeneity validate the proposed theory, showing that the transient anisotropic kernel improves the learned probability measure.

Diffusion maps↗
DOE OSTI2026

Instabilities and phase transitions in architected metamaterials: a gradient-enhanced continuum approach

Architected metamaterials such as foams and lattices exhibit a wide range of properties governed by microstructural instabilities and emerging phase transitions. Their macroscopic response–including energy dissipation during impact, large recoverable deformations, morphing between configurations, and auxetic behavior–remains difficult to capture with conventional continuum models, which often rely on discrete approaches that limit scalability. In this work, we propose a nonlocal continuum formulation that captures both stable and unstable responses of elastic architected metamaterials. The framework extends anisotropic hyperelasticity by introducing nonlocal variables and internal length scales reflective of microstructural features. Local polyconvex free-energy models are systematically augmented with two families of non-(poly)convex energies, enabling both metastable and bistable responses. Implementation in a finite element framework enables solution using a hybrid monolithic–staggered strategy. Simulations capture densification fronts, forward and reverse transitions, hysteresis loops, imperfection sensitivity, and globally coordinated auxetic modes. Overall, this framework provides a robust foundation for accelerated modeling of instability-driven phenomena in architected metamaterials, while enabling extensions to anisotropic, dissipative, and active systems as well as integration with data-driven and machine learning approaches.

42 ENGINEERING↗
DOE OSTI2021

Numerical investigation of active control of tearing mode by magnetic coils and the role of Δ'

Magnetic feedback stabilization of the tearing mode (TM) is numerically investigated, utilizing the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681) for toroidal tokamak equilibria. With control coil configurations assumed in this study, magnetic feedback partially or fully stabilizes the TM, with either vanishing or finite equilibrium pressure. The best control is achieved by the combination of internal active coils and internal poloidal sensors. The internal and external tearing indices are evaluated for the close-loop system, based on the MARS-F computed mode eigenvalue and eigenfunction, respectively. In the absence of the favorable curvature effect, these two indices are real-valued and quantitatively agree well with each other. For the equilibrium with finite pressure gradient at the mode rational surface, the favorable average curvature effect becomes important and the close-loop tearing index also becomes complex-valued, partly due to interaction of the feedback system with the dissipative wall eddy current response. Isolating the inner layer and outer region response to magnetic feedback, with either proportional or proportional-derivative actions, allows to establish that feedback stabilization of the TM occurs mainly due to modification of the behavior of the external ideal solution, further confirming the analytic result reported in He et al 2021 Phys. Plasmas 28 012504.

Physics↗
DOE OSTI2024

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗
DOE OSTI2026

Generating entangled steady states in multistable open quantum systems via initial state control

Entanglement underpins the power of quantum technologies, yet it is fragile and typically destroyed by dissipation. Paradoxically, the same dissipation, when carefully engineered, can drive a system toward robust entangled steady states. However, this engineering task is nontrivial, as dissipative many-body systems are complex, particularly when they support multiple steady states. Here, we derive analytic expressions that predict how the steady state of a system evolving under a Lindblad equation depends on the initial state, without requiring integration of the dynamics. These results extend Refs. [V. V. Albert and L. Jiang, Phys. Rev. A 89, 022118 (2014); V. V. Albert et al., Phys. Rev. X 6, 041031 (2016)], showing that while the steady-state manifold is determined by the Liouvillian kernel, the weights within it depend on both the Liouvillian and the initial state. We identify a special class of Liouvillians for which the steady state depends only on the initial overlap with the kernel. Our framework provides analytical insight and a computationally efficient tool for predicting steady states in open quantum systems. As an application, we propose schemes to generate metrologically useful entangled steady states in spin ensembles via balanced collective decay.

Dissipative dynamics↗
DOE OSTI2022

Polariton enhanced free charge carrier generation in donor–acceptor cavity systems by a second-hybridization mechanism

Cavity quantum electrodynamics has been studied as a potential approach to modify free charge carrier generation in donor–acceptor heterojunctions because of the delocalization and controllable energy level properties of hybridized light–matter states known as polaritons. However, in many experimental systems, cavity coupling decreases charge separation. In this work, we theoretically study the quantum dynamics of a coherent and dissipative donor–acceptor cavity system, to investigate the dynamical mechanism and further discover the conditions under which polaritons may enhance free charge carrier generation. We use open quantum system methods based on single-pulse pumping to find that polaritons have the potential to connect excitonic states and charge separated states, further enhancing free charge generation on an ultrafast timescale of several hundred femtoseconds. The mechanism involves polaritons with optimal energy levels that allow the exciton to overcome the high Coulomb barrier induced by electron–hole attraction. Moreover, we propose that a second-hybridization between a polariton state and dark states with similar energy enables the formation of the hybrid charge separated states that are optically active. These two mechanisms lead to a maximum of 50% enhancement of free charge carrier generation on a short timescale. However, our simulation reveals that on the longer timescale of picoseconds, internal conversion and cavity loss dominate and suppress free charge carrier generation, reproducing the experimental results. Thus, our work shows that polaritons can affect the charge separation mechanism and promote free charge carrier generation efficiency, but predominantly on a short timescale after photoexcitation.

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

Effective field theory for quasicrystals and phasons dynamics

We build an effective field theory (EFT) for quasicrystals -- aperiodic incommensurate lattice structures -- at finite temperature, entirely based on symmetry arguments and a well-define action principle. By means of Schwinger-Keldysh techniques, we derive the full dissipative dynamics of the system and we recover the experimentally observed diffusion-to-propagation crossover of the phason mode. From a symmetry point of view, the diffusive nature of the phason at long wavelengths is due to the fact that the internal translations, or phason shifts, are symmetries of the system with no associated Noether currents. The latter feature is compatible with the EFT description only because of the presence of dissipation (finite temperature) and the lack of periodic order. Finally, we comment on the similarities with certain homogeneous holographic models and we formally derive the universal relation between the pinning frequency of the phonons and the damping and diffusion constant of the phason.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the violation of the zeroth law of turbulence in space plasmas

The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. This behaviour, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the dissipation – no matter how small – can thermalise the energy input. Using direct numerical simulations and theoretical arguments, we show that in strongly magnetised plasma turbulence such as that recently observed by the Parker Solar Probe spacecraft, the zeroth law is routinely violated. Namely, when such turbulence is ‘imbalanced’ – when the large-scale energy input is dominated by Alfvénic perturbations propagating in one direction (the most common situation in space plasmas) – nonlinear conservation laws imply the existence of a ‘barrier’ at scales near the ion gyroradius. This causes energy to build up over time at large scales. The resulting magnetic-energy spectra bear a strong resemblance to those observed in situ , exhibiting a sharp, steep kinetic transition range above and around the ion-Larmor scale, with flattening at yet smaller scales. The effect thus offers a possible solution to the decade-long puzzle of the position and variability of ion-kinetic spectral breaks in plasma turbulence. The existence of the ‘barrier’ also suggests that, how a plasma is forced at large scales (the imbalance) may have a crucial influence on thermodynamic properties such as the ion-to-electron heating ratio.

Physics↗
DOE OSTI2022

Nonadiabatic transition paths from quantum jump trajectories

We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. Here, we detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.

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

TOF-SIMS analysis of decoherence sources in superconducting qubits

Superconducting qubits have emerged as a potentially foundational platform technology for addressing complex computational problems deemed intractable with classical computing. Despite recent advances enabling multiqubit designs that exhibit coherence lifetimes on the order of hundreds of μs, material quality and interfacial structures continue to curb device performance. Two-level system defects in the thin superconducting film and adjacent dielectric regions introduce stochastic noise and dissipate electromagnetic energy at the cryogenic operating temperatures. In this study, we utilize time-of-flight secondary ion mass spectrometry to understand the role specific fabrication procedures play in introducing such dissipation mechanisms in these complex systems. We interrogated Nb thin films and transmon qubit structures fabricated through slight modifications in the processing and vacuum conditions. We find that when the Nb film is sputtered onto the Si substrate, oxide and silicide regions are generated at various interfaces. We also observe that impurity species, such as niobium hydrides and carbides, are incorporated within the niobium layer during the subsequent lithographic patterning steps. The formation of these resistive compounds likely impacts the superconducting properties of the Nb thin film. Additionally, we observe the presence of halogen species distributed throughout the patterned thin films. We conclude by hypothesizing the source of such impurities in these structures in an effort to intelligently fabricate superconducting qubits and extend coherence times moving forward.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗