SEARCH · Engineering Papers

Results for “dissipative dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Effective dissipation rate in a Liouvillian-graph picture of high-temperature quantum hydrodynamics

At high temperature, generic strongly interacting spin systems are expected to display hydrodynamics: local transport of conserved quantities, governed by classical partial differential equations like the diffusion equation. I argue that the emergence of this dissipative long-wavelength dynamics from the system's unitary microscopic dynamics is controlled by the structure of the Liouvillian graph of the system's Hamiltonian, that is, the graph induced on Pauli strings by commutation with that Hamiltonian. The Liouvillian graph decomposes naturally into subgraphs of Pauli strings of constant diameter, and the coherent dynamics of these subgraphs determines the rate at which operator weight spreads to long operators. In conclusion, this argument provides a quantitative theory of the emergence of a dissipative effective dynamics from unitary microscopic dynamics; it also leads to an effective model with Hilbert space dimension linear in system size and exponential in the UV cutoff for diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1990

Dissipative controller designs for second-order dynamic systems

The passivity theorem may be used to design robust controllers for structures with positive transfer functions. This result is extended to more general configurations using dissipative system theory. A stability theorem for robust, model-independent controllers of structures which lack collocated rate sensors and actuators is given. The theory is illustrated for non-square systems and systems with displacement sensors.

Morris, K. A.↗

DOE OSTI2026

Dissipation Enables Strongly Detuning-Dependent Interference in Pulsed Dynamical Decoupling

One of the defining features of pulsed dynamical decoupling is its suppression of a driven qubit's sensitivity to static detuning errors between the drive field and qubit resonance. In this paper, we show that dissipation, in the form of excited-state decay, can change this behavior entirely, producing an interference signal with a strong detuning dependence. This signal arises from decay during driven evolution and relies on coherence retained by the qubit after such a decay event. We develop analytical and numerical models that capture the underlying mechanism and observe this same dissipation-induced detuning dependence experimentally in both a free-space strontium atom interferometer and a superconducting transmon qubit system. We also use this dissipation-induced detuning dependence as the basis for a new spectroscopic technique called Dissipative Carr-Purcell Spectroscopy (DCPS) and compare it with a traditional Ramsey sequence. Our results establish a regime of pulsed dynamical decoupling in which dissipation reshapes, rather than merely degrades, coherent control, and we expect these dynamics to be relevant to a wide range of quantum systems.

DeRose, Kenneth [Northwestern U.; Fermilab] (ORCID↗

DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1990

On the mathematical modeling of the Reynolds stress's equations

By considering the Reynolds stress equations as a possible descriptor of complex turbulent fields, pressure-velocity interaction and turbulence dissipation are studied as two of the main physical contributions to Reynolds stress balancing in turbulent flow fields. It is proven that the pressure interaction term contains turbulence generation elements. However, the usual 'return to isotropy' element appears more weakly than in the standard models. In addition, convection-like elements are discovered mathematically, but there is no mathematical evidence that the pressure fluctuations contribute to the turbulent transport mechanism. Calculations of some simple one-dimensional fields indicate that this extra convection, rather than the turbulent transport, is needed mathematically. Similarly, an expression for the turbulence dissipation is developed. The end result is a dynamic equation for the dissipation tensor which is based on the tensorial length scales.

Lin, Avi↗

NASA NTRS2023

Ames Contrail Simulation Model: Modeling Aviation Induced Contrails and the Computation of Contrail Radiative Forcing Using Air Traffic Data

NASA is developing traffic flow management strategies to reduce the impact of aviation on climate and improve efficiency of aircraft routes in the presence of limited airport and airspace capacity constraints, which are partly due to convective weather and natural calamities. Development of Strategic Traffic Flow Management (STFM) requires models of aircraft emissions and contrails, and models of their effect on climate. There is extensive literature on the modeling of contrails and the Radiative Forcing (RF) associated with contrails. This report captures the latest NASA developments in this research area and describes the integrated modeling, analysis, and software development to support STFM. Aircraft contrails are long, thin and often linear clouds triggered by aircraft engine exhausts in the high-altitude ice-saturated atmosphere. Contrails, similar to that of natural cirrus clouds, can impact global climate by reflecting shortwave radiation and trapping longwave radiation. Recent studies from the Intergovernmental Panel on Climate Change (IPCC) have shown that aircraft contrails are estimated to have greater impact on global warming than aircraft CO 2 emissions. The Ames Contrail Simulation Model (ACSM) presented in this report simulates the full life cycle of aircraft contrails, including their formation, dynamic evolution, and dissipation, and calculates the associated RF with actual meteorological and air traffic data. ACSM combines models from previous studies that focus on contrail formation and persistence based on the Schmidt-Appleman theoretical criteria, and it incorporates elements of cloud dynamics, microphysics, and climate modeling found in other surveyed contrail models, while also making adjustment for improved computational efficiency. In addition, ACSM is integrated with NASA's state-of-the-art flight simulation software for rapid assessment of aircraft contrail impacts. Applications include assessments of long-term global climate impact resulting from aviation-induced contrails and the design of optimal contrail-mitigation aircraft operation strategies.

aviation contrails↗

DOE OSTI2022

Reconciling and Improving Formulations for Thermodynamics and Conservation Principles in Earth System Models (ESMs)

Abstract This paper provides a comprehensive derivation of the total energy equations for the atmospheric components of Earth System Models (ESMs). The assumptions and approximations made in this derivation are motivated and discussed. In particular, it is emphasized that closing the energy budget is conceptually challenging and hard to achieve in practice without resorting to ad hoc fixers. As a concrete example, the energy budget terms are diagnosed in a realistic climate simulation using a global atmosphere model. The largest total energy errors in this example are spurious dynamical core energy dissipation, thermodynamic inconsistencies (e.g., coupling parameterizations with the host model) and missing processes/terms associated with falling precipitation and evaporation (e.g., enthalpy flux between components). The latter two errors are not, in general, reduced by increasing horizontal resolution. They are due to incomplete thermodynamic and dynamic formulations. Future research directions are proposed to reconcile and improve thermodynamics formulations and conservation principles.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Josephson dynamics in two-dimensional ring-shaped condensates

We investigate Josephson transport in a fully closed, two-dimensional superfluid circuit formed by a ring-shaped 87 Rb Bose-Einstein condensate that contains two optical barriers acting as movable weak links. Translating these barriers at controlled speeds imposes a steady bias current, enabling direct mapping of the current-chemical-potential (𝐼−Δ⁢𝜇) characteristics. For narrow junctions (𝑤 ≈ 1µ⁢m), the circuit exhibits a pronounced dc branch that terminates at a critical current 𝐼 𝑐 = 9⁢(1) × 10 3 s −1 ; above this threshold, the system switches to an ac, resistive regime. Classical-field simulations that include the moving barriers quantitatively reproduce both the nonlinear 𝐼−Δ⁢𝜇 curve and the measured 𝐼 𝑐 , validating the underlying microscopic picture. Analysis of the ensuing phase dynamics shows that dissipation is mediated by the nucleation and traversal of vortex-antivortex pairs through the junctions, while the bulk condensate remains globally phase locked—direct evidence of the ring's topological constraint enforcing quantized circulation. These results establish a cold-atom analog of a superconducting quantum interference device in which Josephson dynamics can be resolved at the single-vortex level, providing a versatile platform for atomtronic circuit elements, nonreciprocal Josephson devices, and on-chip Sagnac interferometers for multiaxis rotation sensing.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1999

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

NASA NTRS1990

The route to chaos for the Kuramoto-Sivashinsky equation

The results of extensive numerical experiments of the spatially periodic initial value problem for the Kuramoto-Sivashinsky equation. This paper is concerned with the asymptotic nonlinear dynamics at the dissipation parameter decreases and spatio-temporal chaos sets in. To this end the initial condition is taken to be the same for all numerical experiments (a single sine wave is used) and the large time evolution of the system is followed numerically. Numerous computations were performed to establish the existence of windows, in parameter space, in which the solution has the following characteristics as the viscosity is decreased: a steady fully modal attractor to a steady bimodal attractor to another steady fully modal attractor to a steady trimodal attractor to a periodic attractor, to another steady fully modal attractor, to another periodic attractor, to a steady tetramodal attractor, to another periodic attractor having a full sequence of period-doublings (in parameter space) to chaos. Numerous solutions are presented which provide conclusive evidence of the period-doubling cascades which precede chaos for this infinite-dimensional dynamical system. These results permit a computation of the length of subwindows which in turn provide an estimate for their successive ratios as the cascade develops. A calculation based on the numerical results is also presented to show that the period doubling sequences found here for the Kuramoto-Sivashinsky equation, are in complete agreement with Feigenbaum's universal constant of 4,669201609... . Some preliminary work shows several other windows following the first chaotic one including periodic, chaotic, and a steady octamodal window; however, the windows shrink significantly in size to enable concrete quantitative conclusions to be made.

Papageorgiou, Demetrios T.↗

NASA NTRS1991

The route to chaos for the Kuramoto-Sivashinsky equation

The results of extensive numerical experiments of the spatially periodic initial value problem for the Kuramoto-Sivashinsky equation. This paper is concerned with the asymptotic nonlinear dynamics at the dissipation parameter decreases and spatio-temporal chaos sets in. To this end the initial condition is taken to be the same for all numerical experiments (a single sine wave is used) and the large time evolution of the system is followed numerically. Numerous computations were performed to establish the existence of windows, in parameter space, in which the solution has the following characteristics as the viscosity is decreased: a steady fully modal attractor to a steady bimodal attractor to another steady fully modal attractor to a steady trimodal attractor to a periodic attractor, to another steady fully modal attractor, to another periodic attractor, to a steady tetramodal attractor, to another periodic attractor having a full sequence of period-doublings (in parameter space) to chaos. Numerous solutions are presented which provide conclusive evidence of the period-doubling cascades which precede chaos for this infinite-dimensional dynamical system. These results permit a computation of the length of subwindows which in turn provide an estimate for their successive ratios as the cascade develops. A calculation based on the numerical results is also presented to show that the period doubling sequences found here for the Kuramoto-Sivashinsky equation, are in complete agreement with Feigenbaum's universal constant of 4,669201609 .... Some preliminary work shows several other windows following the first chaotic one including periodic, chaotic, and a steady octamodal window; however, the windows shrink significantly in size to enable concrete quantitative conclusions to be made.

Papageorgiou, Demetrios T.↗

DOE OSTI2025

Electron Energy-Loss Spectroscopy and Differential Phase Contrast Imaging with Active Decision in Multimodal Electron Microscopy: Isotopic detection at the atomic scale

Isotopic engineering provides a powerful route to control phonon behavior in crystalline solids, enabling fundamental studies of lattice dynamics and heat transport at the atomic scale. Here, we directly visualize isotope-dependent phonon propagation in epitaxial Cr 2 O 3 using aberration-corrected scanning transmission electron microscopy (STEM) combined with monochromated, high-energy-resolution electron energy-loss spectroscopy (EELS). Guided by ab initio phonon calculations, we demonstrate that optical phonon modes above 70 meV are predominantly oxygen-derived and exhibit measurable redshifts upon substitution of natural 16 O by enriched 18 O. Spatially resolved vibrational spectrum imaging reveals isotope-enriched tracer layers within Cr 2 O 3 thin films, correlating isotope concentration with phonon intensity variations and vibrational energy shifts. At the nanometer and atomic scales, vibrational EELS mapping uncovers coherent phonon propagation across isotopic interfaces, consistent with theoretical phonon density of states and dispersion relations. These results establish vibrational EELS as a quantitative probe for isotope-dependent phonon transport in materials, opening new possibilities for studying energy dissipation and lattice dynamics.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Phase Contrast X‐Ray Imaging of the Collapse of an Engineered Void in Single‐Crystal HMX

Abstract Imaging the collapse of a single void that creates a hot spot initiation site in an otherwise defect‐free explosive is challenging given the spatial and temporal scales involved in explosive systems. This work presents our attempt to examine a single hot spot mode (void collapse) in single‐crystal octahydro‐l,3,5,7‐tetranitro‐l,3,5,7‐tetrazocine (HMX) embedded in Sylgard. The hot spot heating mechanisms involved with pore collapse include adiabatic heating, jetting, and viscoplastic dissipation. Quantifying the dynamics of a pore collapse is a crucial step to understanding which mechanisms dominate during ignition events. Our experiments were conducted with a single‐stage, light‐gas gun at Argonne National Laboratory's Advanced Photon Source, applying the phase contrast imaging technique while collecting high‐speed video. The details of HMX single crystal production, defect (pore) engineering, and sample construction, along with experimental results are presented here. These results demonstrate that detailed collapse dynamics can be obtained from homogeneous, single‐crystal explosives with this approach. Qualitative comparisons are made with simulation data which show good agreement in the transition between a quasi‐symmetric pore collapse and an asymmetric collapse with jetting across the pore as measured with normalized pore area and pore circularity.

42 ENGINEERING↗

DOE OSTI2023

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗

DOE OSTI2021

On applications of quantum computing to plasma simulations

Quantum computing is gaining increased attention as a potential way to speed up simulations of physical systems, and it is also of interest to apply it to simulations of classical plasmas. However, quantum information science is traditionally aimed at modeling linear Hamiltonian systems of a particular form that is found in quantum mechanics, so extending the existing results to plasma applications remains a challenge. Here, we report a preliminary exploration of the long-term opportunities and likely obstacles in this area. First, we show that many plasma-wave problems are naturally representable in a quantumlike form and thus are naturally fit for quantum computers. Second, we consider more general plasma problems that include non-Hermitian dynamics (instabilities, irreversible dissipation) and nonlinearities. We show that by extending the configuration space, such systems can also be represented in a quantumlike form and thus can be simulated with quantum computers too, albeit that requires more computational resources compared to the first case. Third, we outline potential applications of hybrid quantum–classical computers, which include analysis of global eigenmodes and also an alternative approach to nonlinear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Efficient optimization method for finding minimum energy paths of magnetic transitions

Here, efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden–Fletcher–Goldfarb–Shanno (LBFGS) and conjugate gradient algorithms in conjunction with orthogonal spin optimization (OSO) approach is shown to greatly outperform the previously used velocity projection and dissipative Landau–Lifschitz dynamics optimization methods. The implementation makes use of energy weighted springs for the distribution of the discretization points along the path and this is found to improve performance significantly. The various methods are applied to several test problems using a Heisenberg-type Hamiltonian, extended in some cases to include Dzyaloshinskii–Moriya and exchange interactions beyond nearest neighbours. Minimum energy paths are found for magnetization reversals in a nano-island, collapse of skyrmions in two-dimensional layers and annihilation of a chiral bobber near the surface of a three-dimensional magnet. The LBFGS-OSO method is found to outperform the dynamics based approaches by up to a factor of 8 in some cases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1989

Summary of activities

Four space projects' activities are summarized. Design work on the Mars Penetrator Project, in cooperation with the NASA Ames Research Center, is being continued. Efforts are focused on the drilling mechanism which must penetrate the martian subsurface soil to collect and retrieve an uncontaminated sample. The new design consists of a rotary-percussive drill mechanism. This mechanism is optimum for dry drilling, necessary to avoid contamination of the soil sample, in many different soil types. The need for a small, relatively inexpensive device to study the chemical structure of this martian soil was also established. The egg design was chosen for its low cost compared to other systems and potential for a large number of eggs to be deployed on the martian surface. The design process included analysis of the dynamics of reentry, dissipation of heat during reentry, impact with the surface, access to undisturbed soil samples, and ability to gather samples from the soil at three depths. The egg consists of the reentry systems, soil probe lifting system, soil probe package, gas chromatograph, transmitter, and battery power supply. The egg must function only once, but is designed to withstand one martian year. The Mars Mole is designed as a rover-based device which penetrates the martian soil to a depth of up to 10 m, obtains a sample of soil, and returns it to the surface for analysis. The mole was designed to meet the following specifications: (1) weight less than 10 kg; (2) size less than 20 x 20 x 30 cm; (3) power less than 100 W; (4) ability to obtain a sample of at least 5 cc; (5) ability to penetrate fine, loose sand; and (6) need to obtain at least one sample. The space station umbilical connector project is a device which provides the translational motion of the connectors on the Space Station Freedom to allow engagement for power and data transfer. The design is capable of delivering a 20 lb force within the necessary tolerances and will operate reliably in the space environment.

Source record↗

NASA NTRS1995

Is chemical heating a major cause of the mesosphere inversion layer?

A region of thermal enhancement of the mesosphere has been detected on numerous occasions by in situ measurements, remote sensing from space, and lidar techniques. The source of these 'temperature inversion layers' has been attributed in the literature to the dissipation relating to dynamical forcing by gravity wave or tidal activity. However, evidence that gravity wave breaking can produce the inversion layer with amplitude as large as that observed in lidar measurements has been limited to results of numerical modeling. An alternative source for the production of the thermal inversion layer in the mesosphere is the direct deposition of heat by exothermic chemical reactions. Two-dimensional modeling combining a comprehensive model of the mesosphere photochemistry with the dynamical transport of long-lived species shows that the region from 80 to 95 km may be heated as much as 3 to 10 K/d during the night and half this rate during the day. Given the uncertainties in our understanding of the dynamics and chemistry for the mesopause region, separating the two sources by passive observations of the mesosphere thermal structure looks to be difficult. Therefore we have considered an active means for producing a mesopause thermal layer, namely the release of ozone into the upper mesosphere from a rocket payload. The induced effects would include artificial enhancements of the OH and Na airglow intensities as well as the mesopause thermal structure. The advantages of the rocket release of ozone is that detection of these effects by ground-based imaging, radar, and lidar systems and comparison of these effects with model predictions would help quantify the partition of the artificial inversion layer production into sources of dynamical and chemical forcing.

Meriwether, John W.↗