SEARCH · Search NASA

Results for “entropy viscosity”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

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 OSTI2020

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Effective viscosities in a hydrodynamically expanding boost-invariant QCD plasma

Background: The near-equilibrium properties of a QCD plasma can be encoded into transport coefficients such as bulk and shear viscosity. In QCD, the ratio of these transport coefficients to entropy density, ζ/$\textit{s}$ and η/$\textit{s}$, depends nontrivially on the plasma's temperature. This is unlike in conformal systems where they take constant values such as η/$\textit{s}$ = 1/(4π). Purpose: In this work, we show that in a (0 + 1)D boost-invariant fluid with no transverse expansion, a temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ can be described by an equivalent effective viscosity $\langleζ/s\rangle_{\text{eff}}$ or $\langleη/s\rangle_{\text{eff}}$. This effective viscosity combines the actual temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ with the temperature profile of the fluid. Further, we extend the concept of effective viscosity in systems with transverse expansion and discuss how effective viscosities can be used to identify families of ζ/$\textit{s(T)}$ and η /$\textit{s(T)}$ that lead to similar hydrodynamic evolution. Methods: The Navier-Stokes relativistic hydrodynamic equations are used to provide a first definition of effective viscosity, in (0 + 1)D and (1 + 1)D. In the (0 + 1)D case, the analysis is extended to Israel-Stewart-type second-order hydrodynamics to clarify the effect of higher-order hydrodynamics corrections on the effective viscosity. Results: In a boost-invariant fluid with no transverse expansion [(0 + 1)D], the effective viscosity is expressed as a simple integral of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ over temperature, with a weight determined by the speed of sound of the fluid. The result is general for any equation of state with a moderate temperature dependence of the speed of sound, including the QCD equation of state. This definition of effective viscosity can be used to identify infinite families of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ that produce essentially indistinguishable temperature profiles. In a boost-invariant cylindrical system [(1 + 1)D], a similar definition of effective viscosity is obtained in terms of characteristic trajectories in time and transverse direction. This leads to an infinite number of constraints on an infinite functional space for ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$. Realistic examples are presented by using a finite number of constraints on a finite functional space. Conclusions: The definition of effective viscosity in a (0 + 1)D system clarifies how infinite families of ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$ can result in nearly identical hydrodynamic temperature profiles. By extending the study to a boost-invariant cylindrical [(1 + 1)D] fluid, we identify an approximate but more general definition of effective viscosity that highlights the potential and limits of the concept of effective viscosity in fluids with limited symmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Empirical correlations between the function of entropy ( Z S ) and net artificial viscous work in a shock physics hydrocode

Entropy is a state variable that may be obtained from any thermodynamically complete equation of state (EOS). However, hydrocode calculations that output the entropy often contain numerical errors; this is not because of the EOS, but rather the solution techniques that are used in hydrocodes (especially Eulerian) such as convection, remapping, and artificial viscosity. Here, in this work, empirical correlations are investigated to reduce the errors in entropy without altering the solution techniques for the conservation of mass, momentum, and energy. Specifically, these correlations are developed for the function of entropy Z S , and they depend upon the net artificial viscous work, as determined via Sandia National Laboratories’ shock physics hydrocode CTH. These results are a continuation of a prior effort to implement the entropy-based CREST reactive burn model in CTH, and they are presented here to stimulate further interest from the shock physics community. Future work is planned to study higher-dimensional shock waves, shock wave interactions, and possible ties between the empirical correlations and a physical law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Dicarbonyl Electrolyte for High-Voltage Lithium Metal Batteries: Importance of the Entropy of Solvation in Bidentate Solvent

An electrolyte system based on dicarbonyl solvents has been proposed and implemented in lithium metal batteries (LMBs). Among several dicarbonyl solvents explored, dimethyl dicarbonate (DMDC) emerged as the most promising electrolyte candidate. Interestingly, despite the structural similarities between DMDC and its monocarbonyl counterpart DMC, the DMDC electrolyte exhibits higher solvent coordination, even at the same lithium-to-carbonyl ratio. Further, this leads to the formation of a more robust solid-electrolyte interphase and enhances the anodic stability of the electrolyte. Notably, in various electrolyte systems, the solvation energy is typically enthalpy-driven. However, the bidentate nature of DMDC introduces a significant role for the entropy of solvation because the Li + -DMDC complexes release fewer free molecules. With the addition of non-solvating co-solvent 1,1,2,2-tetrafluoro-ethyl 2,2,3,3-tetrafluoro-propyl ether (TTE), which reduces the viscosity of the electrolyte, the DMDC-TTE system enables remarkably stable cycling of the lithium anode by manipulating the entropy of solvation. This introduces a new strategy for tailoring functional electrolyte systems.

25 ENERGY STORAGE↗

DOE OSTI2022

Glass length: Workability and meltability of glass as a function of glass melt fragility

Workability and meltability of molten glass are properties important for balancing glass formulation with the design and operations of glass melting and forming facilities. The working and melting “lengths” are the temperature intervals within which glass viscosity allows the melting and forming glass to be performed. They are related to glass melt structure through the glass melt fragility, which, by the augmented Adam-Gibbs equation, is a function of configuration entropy. Both working length and melting length are high for strong glasses and decrease as the melt fragility increases. In this work, analytical formulas are obtained for three-parameter viscosity models VFT, AM, and MYEGA. These models agree with each other and with the Adam-Gibbs equation on the relationship between workability and fragility of glasses but diverge at high-temperatures. The compositional effects on workability and meltability are mediated via the fragility as a function of glass composition expressed in terms of mass and mole component coefficients.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Consequences of sp 2 –sp 3 boron isomerization in supercooled liquid borates

In this study, time-resolved high-energy synchrotron x-ray total scattering measurements on supercooled molten lithium metaborate (LiBO 2 ) reveal an isomerization reaction involving conversion of trigonal sp 2 boron to tetrahedral sp 3 boron during quenching and glass formation. Van't Hoff analysis yields an accurate enthalpy change, ΔH = 21(1) kJ mol –1 boron, from which we develop an analytical model for the sp 3 isomer fraction and its contribution to configurational heat capacity (C p conf ) and entropy as a function of temperature and composition. Isomerization constitutes 40% of the total calorimetric C p conf at the glass transition for LiBO 2 and directly contributes to the observed rise in liquid fragility with the lithium content.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Liquid fragility maximum in lithium borate glass‐forming melts related to the local structure

Abstract The structure of liquid lithium pyroborate, Li 4 B 2 O 5 ( J = Li/B = 2), has been measured over a wide temperature range by high‐energy X‐ray diffraction, and compared to that of its glass and borate liquids of other compositions. The results indicate a gradual increase in tetrahedral boron fraction from 3(1)% to 6(1)% during cooling from T = 1271(15) to 721(8) K, consistent with the larger N 4 = 10(1)% found for the glass, and literature 11 B nuclear magnetic resonance measurements. van't Hoff analysis based on a simple boron isomerization reaction BØ 3 O 2 – ⇌ BØO 2 2– yields Δ H = 13(1) kJ mol –1 and Δ S = 40(1) J mol –1 K –1 for the boron coordination change from 4 to 3, which are, respectively, smaller and larger than found for singly charged isomers for J ≤ 1. With these, we extend our model for N 4 ( J , T ), nonbridging oxygen fraction f nbr ( J , T ), configurational heat capacity , and entropy S conf ( J , T ) contributions up to J = 3. A maximum is revealed in at J = 1, and shown semi‐quantitatively to lead to a corresponding maximum in fragility contribution, akin to that observed in the total fragilities by temperature‐modulated differential scanning calorimetry. Lithium is bound to 4.6(2) oxygen in the pyroborate liquid, with 2.7(1) bonds centered around 1.946(8) Å and 1.9(1) around 2.42(1) Å. In the glass, n LiO = 5.4(4), the increase being due to an increase in the number of short Li–O bonds.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A general metriplectic framework with application to dissipative extended magnetohydrodynamics

General equations for conservative yet dissipative (entropy producing) extended magnetohydrodynamics are derived from two-fluid theory. Keeping all terms generates unusual cross-effects, such as thermophoresis and a current viscosity that mixes with the usual velocity viscosity. While the Poisson bracket of the ideal version of this model has already been discovered, we determine its metriplectic counterpart that describes the dissipation. This is done using a new and general thermodynamic point of view to derive dissipative brackets, a means of derivation that is natural for understanding and creating dissipative dynamics without appealing to underlying kinetic theory orderings. Finally, the formalism is used to study dissipation in the Lagrangian variable picture where, in the context of extended magnetohydrodynamics, non-local dissipative brackets naturally emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Transverse momentum fluctuations and their correlation with elliptic flow in nuclear collision

In this article, we propose observables v 0 and v 0 (p T ) which quantify the relative fluctuations in the total transverse momentum at fixed multiplicity. We first study the factorization of the fixed multiplicity momentum dependent two particle correlation function into a product of v 0 ($p^{a}_{T}$) and v 0 ($p^{b}_{T}$) within realistic hydrodynamic simulations. Then we present computations of v 0 (p T ) for different particle types. We determine the relation between the integrated v 0 and previously measured observables, and compare results from a hybrid hydrodynamics based model to experimental data. The effects of bulk viscosity and an initial pre-equilibrium stage on the results are quantified. We find that v 0 is strongly correlated with the initial state entropy per elliptic area, S/A. Using this result, we explain how the observed correlations between the elliptic flow and the transverse momentum (both in simulations and experiment) reflect the initial state correlations between 1/A and ellipticity ε 2 at fixed multiplicity. Furthermore, we argue that the systematic experimental study of v 0 , with the same sophistication as used for the other v n , can contribute significantly to our understanding of quark gluon plasma properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Dual modulation of the anion-driven thermodynamic properties of aqueous choline halide-based deep eutectic solvents

Deep eutectic solvents (DESs) are considered tunable solvents because their specific properties can be achieved based on the choice of components and their relative concentrations in a mixture. In this article, we investigate the influence of the variation in halide ions (F − , Cl − , Br − , I − ) of choline salts used on the thermodynamic and physicochemical properties of choline halide-based DESs. Our findings show that the density of choline halide-based DESs decreases nonlinearly with an increasing mole fraction of water, following a trend based on the size of the halides, with choline iodide showing the highest density. Temperature-dependent density data reveal that the thermal expansion coefficient decreases slightly with increasing water content, indicating a more stable volume at a higher mole fraction of water. The excess molar volume (V E ) of the DES mixtures exhibits complex behavior depending on the choline halide used, with both negative and positive V E values observed across different water mole fractions. These variations are linked to the hydrogen bonding interactions between the DES components and water molecules. In addition, viscosities decrease with increasing water content, suggesting the disruption of hydrogen bonding networks and enhanced mobility of the ions, which contributes to the observed increase in conductivity. The excess molar Gibbs energies, enthalpies, and entropies of activation have also been determined.

Choline halide↗

DOE OSTI2023

Effects of Methylating Imidazolium-Based Ionic Liquids on Viscosity: New Insights from the Compensated Arrhenius Formalism

Methylation of the C(2) carbon on imidazolium-based room temperature ionic liquids (RTILs) causes an unexpected increase in viscosity when paired with the anion bis(trifluoromethylsulfonamide) [Tf 2 N] – , but the viscosity decreases when the methylated imidazolium is paired with a tetracyanoborate [B(CN) 4 ] – anion. Herein, this paper investigates these different observations in viscosity using the compensated Arrheniusformalism (CAF) for fluidity (inverse viscosity), which assumes fluidity to be a thermally activated process. CAF activation energies are determined for imidazolium [Tf 2 N] – and methylated imidazolium [Tf 2 N] – and compared to imidazolium [B(CN) 4 ] – and methylated imidazolium [B(CN) 4 ] – . The results show that the activation energy increases with methylation for [Tf 2 N] – , but it decreases with methylation for [B(CN) 4 ] – . The CAF results also yield information concerning the entropy of activation, which are compared for the two systems

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗

DOE OSTI2025

Shear and bulk viscosity for a pure glue theory using an effective matrix model

At nonzero temperatures, the deconfining phase transition can be analyzed using an effective matrix model to characterize the change in holonomy. The model includes gluons and two-dimensional ghost fields in the adjoint representation, or “teens.” As ghosts, the teen fields are responsible for the decrease of the pressure as 𝑇 →𝑇 𝑑 , with 𝑇 𝑑 the transition temperature for deconfinement. Using the solution of this matrix model for a large number of colors, the parameters of the teen fields are adjusted so that the expectation value of the Polyakov loop is close to the values from the lattice. The shear, 𝜂, and bulk, 𝜁, viscosities are computed at nonzero holonomy to leading logarithmic order in weak coupling. In the pure glue theory, the value of the Polyakov loop is relatively large in the deconfined phase, ≈1/2 at 𝑇 𝑑 . Consequently, if 𝑠 is the entropy density, while 𝜂/𝑠 decreases as 𝑇 →𝑇 𝑑 , it is still well above the conformal bound. In contrast, 𝜁/𝑠 is largest at 𝑇 𝑑 , comparable to 𝜂/𝑠, then falls off rapidly with increasing temperature and is negligible by ∼2⁢𝑇 𝑑 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Deep glassy state dynamic data challenge glass models: Elastic models

The idea of an “ideal” glass transition temperature has persisted at least since the work from Kauzmann when it was observed that the entropy of glass-forming liquids extrapolated to below that of the crystal, thus suggesting the need for a phase transition at a finite or non-zero absolute temperature. This thermodynamic paradox was also found to be related to the observations of a diverging of the extrapolated viscosity or relaxation times at a temperature near to this ideal glass transition. Recently, however, we have carried out experiments using both an ancient amber material and an ultra-stable amorphous fluoropolymer that challenge the ideas of the divergence of the viscosity or relaxation times at this ideal glass transition. In the present manuscript we have evaluated two theories of the glass transition that are based on ideas related to elasticity of the amorphous glass-forming material. We find that the models from both J. Dyre and his group and of K.S. Schweizer and his group not only show non-diverging behavior but are also in some agreement with the new data, though still showing somewhat slower dynamics than those observed in the experiments. The work shows that the data are good enough to distinguish between theories and it is suggested that other mechanisms may be needed to fully describe the non-diverging responses of the ultra-stable glasses.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Stochastic relativistic viscous hydrodynamics from the Metropolis algorithm

We propose an algorithm for simulating stochastic relativistic fluid dynamics based on Metropolis updates. Each step of the algorithm begins with an update based on ideal hydrodynamics. This is followed by proposing random (spatial) momentum transfers between fluid cells, keeping the total energy fixed. These proposals are then accepted or rejected using the change in entropy as a statistical weight. The algorithm reproduces relativistic viscous hydrodynamics in the “density frame,” which is a formulation of viscous hydrodynamics we review and clarify here. This formulation is first order in time and requires no auxiliary dynamical fields such as Π 𝜇⁢𝜈 . The only parameters are the shear and bulk viscosities and the equation of state. Here, by adopting the 3+1 split of general relativity, we extend the Metropolis algorithm to general space-time coordinates, such as Bjorken coordinates, which are commonly used to simulate heavy-ion collisions.

Hydrodynamic noise↗

DOE OSTI2023

Late-time post-merger modeling of a compact binary: effects of relativity, r-process heating, and treatment of transport

Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first 10 s after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximations of axisymmetry and modeling turbulent angular momentum transport by a viscosity, previous simulations often make other simplifications related to the treatment of the equation of state and turbulent transport effects. In this paper, we test the effect of these modeling choices. By evolving with the same viscosity the exact post-merger initial configuration previously evolved in Newtonian viscous hydrodynamics, we find that the Newtonian treatment provides a good estimate of the disk ejecta mass but underestimates the outflow velocity. We find that the inclusion of heavy nuclei causes a notable increase in ejecta mass. An approximate inclusion of r-process effects has a comparatively smaller effect, except for its designed effect on the composition. Diffusion of composition and entropy, modeling turbulent transport effects, has the overall effect of reducing ejecta mass and giving it a speed with lower average and more tightly-peaked distribution. Also, we find significant acceleration of outflow even at distances beyond 10 000 km, so that thermal wind velocities only asymptote beyond this radius and at higher values than often reported.

79 ASTRONOMY AND ASTROPHYSICS↗