Engineering Papers⌕ Search

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.

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

At least 271 records · Page 15

Clear air turbulence as a response to meso- and synoptic-scale dynamic processes

The paper presents both the Roach equation for the rate of energy dissipation due to clear air turbulence and Richardson number tendencies in isentropic coordinates and examines the implications of these formulations to determine whether there is a dynamic interdependence between Ri and the nonturbulent deformation processes. The equation representing the ln(Ri) tendency is applied diagnostically to grids from an isentropic analysis of archived soundings. The evolution of the Richardson number fields over 12-hour time periods is examined using a mechanistic model. It is suggested that the application of the Roach equation for the turbulent dissipation rate should have a more restricted use. Analyses of the meso-alpha scale Richardson number and of the Richardson number tendency fields reveal a phase relationship consistent with the theoretical predictions.

Keller, John L.↗

Lyapunov exponent as a signature of dissipative many-body quantum chaos

A distinct feature of Hermitian quantum chaotic dynamics is the exponential increase of certain out-of-time-order correlation (OTOC) functions around the Ehrenfest time with a rate given by a Lyapunov exponent. Physically, the OTOCs describe the growth of quantum uncertainty that crucially depends on the nature of the quantum motion. Here, we employ the OTOC in order to provide a precise definition of dissipative quantum chaos. For this purpose, we compute analytically the Lyapunov exponent for the vectorized formulation of the large- q limit of a q -body Sachdev-Ye-Kitaev model coupled to a Markovian bath. These analytic results are confirmed by an explicit numerical calculation of the Lyapunov exponent for several values of q ≥ 4 based on the solutions of the Schwinger-Dyson and Bethe-Salpeter equations. We show that the Lyapunov exponent decreases monotonically as the coupling to the bath increases and eventually becomes negative at a critical value of the coupling signaling a transition to a dynamics which is no longer quantum chaotic. Therefore, a positive Lyapunov exponent is a defining feature of dissipative many-body quantum chaos. The observation of the breaking of the exponential growth for sufficiently strong coupling suggests that dissipative quantum chaos may require in certain cases a sufficiently weak coupling to the environment. Published by the American Physical Society 2024

García-García, Antonio M.↗

Analysis of Transition-Sensitized Turbulent Transport Equations

The dynamics of an ensemble of linear disturbances in boundary-layer flows at various Reynolds numbers is studied through an analysis of the transport equations for the mean disturbance kinetic energy and energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field even in proximity to the wall, the early stage transition regime fluctuations studied here are influenced cross the boundary layer by the solid boundary. The dominating dynamics in the disturbance kinetic energy and dissipation rate equations are described. These results are then used to formulate transition-sensitized turbulent transport equations, which are solved in a two-step process and applied to zero-pressure-gradient flow over a flat plate. Computed results are in good agreement with experimental data.

Rumsey, Christopher L.↗

Metriplectic four-bracket algorithm for constructing thermodynamically consistent dynamical systems

A unified thermodynamic algorithm is presented for constructing thermodynamically consistent dynamical systems, i.e., systems that have Hamiltonian and dissipative parts that conserve energy while producing entropy. The algorithm is based on the metriplectic 4-bracket given in Morrison and Updike [Phys. Rev. E 109, 045202 (2024)]. A feature of the unified thermodynamic algorithm is the force-flux relation 𝐉 𝛼 =−𝐿 𝛼⁢𝛽 𝛁(𝛿⁢𝐻⁡/𝛿⁢𝜉 𝛽 ) for phenomenological coefficients 𝐿 𝛼⁢𝛽 , Hamiltonian 𝐻, and dynamical variables 𝜉 𝛽 . The algorithm is applied to the Navier-Stokes-Fourier, the Cahn-Hilliard-Navier-Stokes, and Brenner-Navier-Stokes-Fourier systems, and significant generalizations of these systems are obtained.

Multiphase flows↗

The CO Emission in the Taffy Galaxies (UGC 12914/15) at 60 pc Resolution. I. The Battle for Star Formation in the Turbulent Taffy Bridge

We present Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of $0^{"}_{.}2$ (60 pc) of CO emission from the Taffy galaxies (UGC 12914/5). The observations are compared with narrowband Paα, mid-IR, radio continuum and X-ray imaging, plus optical spectroscopy. The galaxies have undergone a recent head-on collision, creating a massive gaseous bridge that is known to be highly turbulent. The bridge contains a complex web of narrow molecular filaments and clumps. The majority of the filaments are devoid of star formation, and fall significantly below the Kennicutt–Schmidt relationship for normal galaxies, especially for the numerous regions undetected in Paα emission. Within the loosely connected filaments and clumps of gas we find regions of high velocity dispersion that appear gravitationally unbound for a wide range of likely values of X CO . Like the "Firecracker" region in the Antennae system, they would require extremely high external dynamical or thermal pressure to stop them dissipating rapidly on short crossing timescales of 2–5 Myr. We suggest that the clouds may be transient structures within a highly turbulent multiphase medium that is strongly suppressing star formation. Despite the overall turbulence in the system, stars seem to have formed in compact hotspots within a kiloparsec-sized extragalactic H ii region, where the molecular gas has a lower velocity dispersion than elsewhere, and shows evidence for a collision with an ionized gas cloud. Like the shocked gas in the Stephan's Quintet group, the conditions in the Taffy bridge shows how difficult it is to form stars within a turbulent, multiphase, gas.

79 ASTRONOMY AND ASTROPHYSICS↗

A 3-dimensional Navier-Stokes-Euler code for blunt-body flow computations

The shock-layer flowfield is obtained with or without viscous and heat-conducting dissipations from the conservative laws of fluid dynamics equations using a shock-fitting implicity finite-difference technique. The governing equations are cast in curvilinear-orthogonal coordinates and transformed to the domain between the shock and the body. Another set of equations is used for the singular coordinate axis, which, together with a cone generator away from the stagnation point, encloses the computation domain. A time-dependent alternating direction implicit factorization technique is applied to integrate the equations with local-time increment until a steady solution is reached. The shock location is updated after the flowfield computation, but the wall conditions are implemented into the implicit procedure. Innovative procedures are introduced to define the initial flowfield, to treat both perfect and equilibrium gases, to advance the solution on a coarse-to-fine grid sequence, and to start viscous flow computations from their corresponding inviscid solutions. The results are obtained from a grid no greater than 28 by 18 by 7 and converged within 300 integration steps. They are of sufficient accuracy to start parabolized Navier-Stokes or Euler calculations beyond the nose region, to compare with flight and wind-tunnel data, and to evaluate conceptual designs of reentry spacecraft.

Li, C. P.↗

Magnetic blowout during collisionless reconnection

Analytic calculations and simulations of the steady-state structure of the dissipation region during collisionless reconnection are described. Ions are treated as particles and electrons as a dissipative fluid. A limiting reconnection rate results as the ions 'blowout' the magnetic field in the dissipation area. The implications for comprehending substorm dynamics are discussed.

Drake, J. F.↗

Research on Earth's rotation and the effect of atmospheric pressure on vertical deformation and sea level variability

The work done under NASA grant NAG5-485 included modelling the deformation of the earth caused by variations in atmospheric pressure. The amount of deformation near coasts is sensitive to the nature of the oceanic response to the pressure. The PSMSL (Permanent Service for Mean Sea Level) data suggest the response is inverted barometer at periods greater than a couple months. Green's functions were constructed to describe the perturbation of the geoid caused by atmospheric and oceanic loading and by the accompanying load-induced deformation. It was found that perturbation of up to 2 cm are possible. Ice mass balance data was used for continental glaciers to look at the glacial contributions to time-dependent changes in polar motion, the lod, the earth's gravitational field, the position of the earth's center-of-mass, and global sea level. It was found that there can be lateral, non-hydrostatic structure inside the fluid core caused by gravitational forcing from the mantle, from the inner core, or from topography at the core/mantle or inner core/outer core boundaries. The nutational and tidal response of a non-hydrostatic earth with a solid inner core was modeled. Monthly, global tide gauge data from PSMSL was used to look at the 18.6-year ocean tide, the 14-month pole tide, the oceanic response to pressure, the linear trend and inter-annual variability in the earth's gravity field, the global sea level rise, and the effects of post glacial rebound. The effects of mantle anelasticity on nutations, earth tides, and tidal variation in the lod was modeled. Results of this model can be used with Crustal Dynamics observations to look at the anelastic dissipation and dispersion at tidal periods. The effects of surface topography on various components of crustal deformation was also modeled, and numerical models were developed of post glacial rebound.

Wahr, John↗

A New K-epsilon Eddy Viscosity Model for High Reynolds Number Turbulent Flows: Model Development and Validation

A new k-epsilon eddy viscosity model, which consists of a new model dissipation rate equation and a new realizable eddy viscosity formulation, is proposed. The new model dissipation rate equation is based on the dynamic equation of the mean-square vorticity fluctuation at large turbulent Reynolds number. The new eddy viscosity formulation is based on the realizability constraints: the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. We find that the present model with a set of unified model coefficients can perform well for a variety of flows. The flows that are examined include: (1) rotating homogeneous shear flows; (2) boundary-free shear flows including a mixing layer, planar and round jets; (3) a channel flow, and flat plate boundary layers with and without a pressure gradient; and (4) backward facing step separated flows. The model predictions are compared with available experimental data. The results from the standard k-epsilon eddy viscosity model are also included for comparison. It is shown that the present model is a significant improvement over the standard k-epsilon eddy viscosity model.

Shih, T.-H.↗

Mapping protein dynamics at high spatial resolution with temperature-jump X-ray crystallography

Understanding and controlling protein motion at atomic resolution is a hallmark challenge for structural biologists and protein engineers because conformational dynamics are essential for complex functions such as enzyme catalysis and allosteric regulation. Time-resolved crystallography offers a window into protein motions, yet without a universal perturbation to initiate conformational changes the method has been limited in scope. Here we couple a solvent-based temperature jump with time-resolved crystallography to visualize structural motions in lysozyme, a dynamic enzyme. We observed widespread atomic vibrations on the nanosecond timescale, which evolve on the submillisecond timescale into localized structural fluctuations that are coupled to the active site. An orthogonal perturbation to the enzyme, inhibitor binding, altered these dynamics by blocking key motions that allow energy to dissipate from vibrations into functional movements linked to the catalytic cycle. Because temperature jump is a universal method for perturbing molecular motion, the method demonstrated here is broadly applicable for studying protein dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

Magnetic storm dynamics of the thermosphere

A theoretical study of the Dst component of magnetic storms is presented. The dynamic characteristics are found significantly different for Joule dissipation and electron precipitation, leading to the conclusion that the former is probably the predominant heat source for the upper thermosphere. Composition measurements on OGO-6, which reveal markedly different characteristics in N2, O and He, can be explained on the basis of energy advection and diffusive mass transport by thermospheric winds. Essential features in the F2-region response are explicable in terms of these dynamic processes. Electric field induced motions are estimated and it is concluded that resultant adiabatic heating could be significant.

Mayr, H. G.↗

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗

Relative permeability as a stationary process: Energy fluctuations in immiscible displacement

Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work, we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) there can be no net change for extensive measures of the system state over the time averaging interval; (2) the net energy inputs into the system are zero, meaning that the net rate of work done on the system must balance with the heat removed; and (3) there is no net work performed due to the contribution of internal energy fluctuations. Results are then evaluated based on direct numerical simulation. Dynamic connectivity is observed during steady-state flow, which is quantitatively assessed based the Euler characteristic. We show that even during steady-state flow at low capillary number ([Formula: see text]), typical flow processes will explore multiple connectivity states. The residence time for each connectivity state is captured based on the time-and-space average. The distribution for energy fluctuations is shown to be multi-modal and non-Gaussian when terms are considered independently. However, we demonstrate that their sum is zero. Given an appropriate choice of the thermodynamic driving force, we show that the conventional relative permeability relationship is sufficient to model the energy dissipation in systems with complex pore-scale dynamics that routinely alter the structure of fluid connected pathways.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ac losses in field-cooled type I superconducting cavities

As superconductors are cooled below their critical temperature, stray magnetic flux can become trapped in regions that remain normal. The presence of trapped flux facilitates dissipation of ac current in a superconductor, leading to losses in superconducting elements of microwave devices. In type II superconductors, dissipation is well-understood in terms of the dynamics of vortices hosting a single flux quantum. In contrast, the ac response of type I superconductors with trapped flux has not received much attention. Building on Andreev's early work (Andreev 1967 Sov. Phys. JETP 24 1019), here we show theoretically that the dominant dissipation mechanism is the absorption of the ac field at the exposed surfaces of the normal regions, while the deformation of the superconducting/normal interfaces is unimportant. We use the developed theory to estimate the degradation of the quality factors in field-cooled cavities, and we satisfactorily compare these theoretical estimates to the measured field dependence of the quality factors of two aluminum cavities. We also identify a regime in which the dissipated power depends weakly on the Ginzburg-Landau parameter; this makes it possible to apply our findings to cavities made of other materials, such as niobium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Molecular Picture of the Transient Nature of Lipid Rafts

In biological membranes, lipid rafts are now thought to be transient and nanoscopic. However, the mechanism responsible for these nanoscopic assemblies remains poorly understood, even in the case of model membranes. Therefore, it has proven extremely challenging to probe the physicochemical properties of lipid rafts at the molecular level. Here, we use all-atom molecular dynamics (MD) simulations and inelastic X-ray scattering (IXS), an intrinsically nanoscale technique, to directly probe the energy transfer and collective short-wavelength dynamics (phonons) of biologically relevant model membranes. We show that the nanoscale propagation of stress in lipid rafts takes place in the form of collective motions made up of longitudinal (compression waves) and transverse (shear waves) molecular vibrations. Importantly, we provide a molecular picture for the so-called van der Waals mediated "force from lipid" [Anishkin, A. et al. Proc. Natl. Acad. Sci. U.S.A. 2014, 111, 7898], a key parameter for the ionic channel mechano-transduction and the mechanism for the lipid transfer of molecular level stress [Aponte-Santamaría, C. et al. J. Am. Chem. Soc. 2017, 139, 13588]. Specifically, we describe how lipid rafts are formed and maintained through the propagation of molecular stress, lipid raft rattling dynamics, and a relaxation process. Eventually, the rafts dissipate through the self-diffusion of lipids making up the rafts. We also show that the molecular stress and viscoelastic properties of transient lipid rafts can be modulated through the use of hydrophobic biomolecules such as melatonin and tryptophan. Overall, the herein proposed mechanism describing the molecular interactions for the formation and dissolution of lipid rafts may offer insights as to how lipid rafts enable biological function.

36 MATERIALS SCIENCE↗

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Evolution of the lunar orbit with temperature- and frequency-dependent dissipation

The importance of the thermal evolution of the earth and the moon for the dynamical history of the earth-moon system is investigated using a model of coupled thermal-dynamical evolution of this system, in which the sun's gravity is ignored and only the gravitational interaction of the earth-moon system is included, along with oceanic dissipation. Results suggest that most of the lunar angular momentum was received from the earth via solid tides within the first few 100 million years after the formation of the earth and the moon, indicating that oceanic dissipation plays a minor role. It is suggested that solid tides, rather than ocean tides, might be most important to the lunar-orbit evolution.

Ross, M. N.↗