Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-equilibrium thermodynamics”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗
DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
NASA NTRS2019

Hybrid Simulation of Solar-Wind-Like Turbulence

We present 2.5D hybrid simulations of the spectral and thermodynamic evolution of an initial state of magnetic field and plasma variables that in many ways represents solar wind fluctuations. In accordance with Helios near-Sun high-speed stream observations, we start with Alfvénic fluctuations along a mean magnetic field in which the fluctuations in the magnitude of the magnetic field are minimized. Since fluctuations in the radial flow speed are the dominant free energy in the observed fluctuations, we include a field-aligned v∥(k⊥) with an k(exp −1) spectrum of velocity fluctuations to drive the turbulent evolution. The flow rapidly distorts the Alfvénic fluctuations, yielding spectra (determined by spacecraft-like cuts) transverse to the field that become comparable to the k∥ fluctuations, as in spacecraft observations. The initial near constancy of the magnetic field is lost during the evolution; we show this also takes place observationally. We find some evolution in the anisotropy of the thermal fluctuations, consistent with expectations based on Helios data. We present 2D spectra of the fluctuations, showing the evolution of the power spectrum and cross-helicity. Despite simplifying assumptions, many aspects of simulations and observations agree. The greatly faster evolution in the simulations is at least in part due to the small scales being simulated, but also to the non-equilibrium initial conditions and the relatively low overall Alfvénicity of the initial fluctuations.

Roberts, D. Aaron↗
DOE OSTI2026

A modified scheil approach for nucleation-dependent solidification pathways

As-solidified microstructures of near-eutectic alloys often contain multiple primary phases that are not expected from equilibrium phase diagrams. Such microstructures are caused by cooling-rate-dependent solidification pathways, a factor not captured by the Scheil–Gulliver model or variations thereof. Here, we present a model and algorithm that incorporate the critical nucleation undercooling for each solid phase into the Scheil–Gulliver model. We hypothesize that the non-equilibrium microstructure formation is primarily governed by a nucleation-competition mechanism. This mechanism accounts for both stable/metastable phase selection and primary-phase formation within eutectic regions driven by asymmetric nucleation barriers. The model is validated against a hypereutectic Al-Fe alloy, where it successfully reproduces the observed microstructural constituents, revealing the key dependencies of solidification microstructure on nucleation kinetics. Applicability to multicomponent systems is demonstrated through a hypereutectic Al–Fe–Si ternary alloy, where the model successfully predicts divorced eutectic microstructures and the associated oscillatory solidification pathways along univariant lines. As a result, the proposed framework establishes a nucleation-dependent computational approach for interpreting and predicting solidification microstructures.

Alloy design↗
DOE OSTI2024

Flame kinetics at scramjet-engine-relevant conditions: Role of prompt dissociation of weakly-bound radicals

Combustion in high-speed ram-based propulsion engines occurs under distinct thermodynamic conditions of high reactant temperatures (greater than 1000 K) and relatively low pressures (<5 atm). There is a lack of fundamental flame measurements at such conditions that result in adiabatic flame temperatures (T ad ) exceeding 2500 K. In this work, we have measured laminar flame speeds of oxygen-enriched CH 4 /oxidizer mixtures at sub-atmospheric conditions to probe kinetics at high T ad using the isobaric spherically expanding flame approach. Simulations with recent kinetic models revealed increasing differences between data and model predictions with increasing T ad , reaching up to 25 %. Kinetic analyses reveal that at the thermodynamic conditions in these O 2 -enriched flames, i.e., lower pressures and higher T ad , the effects of HCO prompt dissociation are accentuated. In addition to HCO, the prompt dissociations of CH 2 OH and C 2 H 5 are also considered. Here, the prompt dissociations of all three radicals were evaluated and their effects considered in flame speed simulations. Reaction path analysis for the present flames revealed that approximately half of the reaction flux for HCO formation undergoes prompt dissociation to H + CO. Furthermore, these analyses also revealed that the pathways and sensitive reactions are similar between oxygen-enriched fuel/oxidizer mixtures and preheated fuel/air mixtures, if both have similar T ad . Thus, flames of oxygen-enriched mixtures could be a surrogate to probe the flame chemistry of highly preheated mixtures at relatively low pressures that are often encountered in ram-based propulsion engine combustors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

Large Solar-Rejection Filter

NASA utilized Image Intensified Video Cameras for ATV data acquisition from a jet flying at 12.8 km. Afterwards the video was digitized and then analyzed with a modified commercial software package, Image Systems Trackeye. Astrometric results were limited by saturation, plate scale, and imposed linear plate solution based on field reference stars. Time-dependent fragment angular trajectories, velocities, accelerations, and luminosities were derived in each video segment. It was evident that individual fragments behave differently. Photometric accuracy was insufficient to confidently assess correlations between luminosity and fragment spatial behavior (velocity, deceleration). Use of high resolution digital video cameras in future should remedy this shortcoming. A scenically accurate description of matter interpreted as a substance made up of corpuscular constituents was established during the course of the 19th century. In this description, atoms--the building blocks of the matter--form molecules. The properties of the molecules were described by chemistry or thermodynamics depending on what characteristics of the matter were investigated. In both theories, the molecules can dissociate to atoms when the kinetic energies of the atoms exceed the strength of the chemical bonds. The number of atoms is always preserved in a closed system. This is not true, however, when the matter takes up much higher energies at relativistic scales. New particles can be produced at the expense of the kinetic energy. The number of particles is no longer preserved. There are other conserved quantities, however, these quantities, the charge, baryon number, lepton number, are associated with particles that are considered elementary today. The properties and behavior of these elementary particles is the subject of Particle Physics or High Energy Physics. Practice Page: A crewmember touching a positively charged surface was thought to be galvanically isolated from the vehicle ground analogous to a bird on a high voltage power wire. Recent analysis confirms that positive floating potentials, ionospheric currents to the EVA suit, can be hazardous. The analysis is wrong in that the ionospheric plasma itself can close the circuit. Parametric analysis of very low voltage exposures (2 to 15 volts) could cause pain and/or involuntary muscle tetani or spinal cord shock. NASA worked with the Naval Health Research Center Detachment Directed Energy Bioeffects Laboratory to examine the affects electrical hazards could have on extravehicular activity using two models. The results of the two computational models were combined to predict areas of the body in which neurons of different diameters would be excited. They predicted that physiologically active current could be conducted across the crew member causing catastrophic hazards. Future work to analyze additional current paths was proposed. The FUSE spectrum of BB Dor, observed in a high state, is modeled with an accretion disk with a very low inclination (possibly lower than 10 degrees). Assuming an average WD mass of 0.8 solar mass leads to a distance of the order of approximately 650pc, consistent with the extremely low galactic reddening in its direction, and a mass accretion rate of 10 (exp -9) solar mass a year. The spectrum presents some broad and deep silicon and sulfur absorption lines, indicating that these elements are over-abundant: silicon is 3 times solar, and sulfur is 20 times solar. The FUSE spectrum of BB Dor, observed in a high state, is modeled with an accretion disk with a very low inclination (possibly lower than 10 degrees). Assuming an average WD mass of 0.8 solar mass leads to a distance of the order of approximately 650pc, consistent with the extremely low galactic reddening in its direction, and a mass accretion rate of 10 (exp -9) solar mass a year. The spectrum presents some broad and deep silicon and sulfur absorption lines, indicating that these elements are over-abundant: silicon is times solar, and sulfur is 20 times solar. From its previously observed optical emission lines, P831-57 (WD 0334 6400 or Ret 1 in A Catalog and Atlas of Cataclysmic Variables: Living Edition) has been suspected to contain an accretion disk associated with a companion star in orbit around a subdwarf star with a temperature T is greater than 21,000K. P831-57 has therefore been classified as a nova-like. However, our present observations show it to be a DA + dMe binary. The analysis of its Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum (continuum and lines) reveal an average mass white dwarf (Log(g) approximately equals 7.8 plus or minus 0.1)with a temperature T approximately equals 37,000 plus or minus 500K, an extremely low projected rotational velocity, and a distance of about 115 plus or minus 5pc. The photosphere contains C, N, Si, and S (at about 1% of solar abundances). The dMe star is seen as a flux excess in near-infrared photometry and appears to show occasional flaring of about one magnitude as seen in the Harvard plates. There is no evidence of periodic variability in the spectroscopic or photometric data. We find no evidence of a an accretion disk, instead we find evidence of wind accretion as the stellar carbon abundance (N(C)/N(H) = 2.5 x 10 (exp -6) is about ten times larger than predicted by radiative levitation for such a gravity and temperature. The power needs and solutions for the space exploration and lunar mobility program are discussed. Long term missions in space and on the lunar surface require high energy batteries. Rechargeable batteries for mobility systems and portable utility pallet are needed for successful exploration missions. Nanomaterial usage increases the energy density of the cells apart from increasing the power density. The symptoms and threats from acute mountain sickness (AMS) are discussed. The underlying assumptions concerning spacecraft atmosphere mean there is a potential risk to astronauts. The baseline worst case scenario is about 25% based on direct ascent to 8.0 psia with 32% O2. The practical mitigation approach is a staged depressurization scheme. Due to the uncertainty about potential AMS risk the flight surgeon should prepare, focused research should proceed, and current analytical efforts should continue. This slide presentation reviews the threat astronauts face from acute mountain sickness (AMS). The incidence of AMS is highly variable and the risk to astronauts is unclear. However, the potential risk is higher than current EAAs suggest and precautions should be implemented. Due to the uncertain nature of the risk flight surgeons should prepare, and focused research should proceed. Quasiparticle tunneling can be understood at low temperatures using a non-equilibrium kinetic theory. Low-temperature tunneling can be suppressed using existing techniques: island-lead gap engineering, lead cooling with SIN junctions, and quasiparticle traps. The slide presentation reviews the science goals and mission planning for the proposed Jupiter System Observer (JSO). The spacecraft would accommodate 9 instruments making the instrument payload more capable than any previously flown to Jupiter. The would-be mission includes study programs for Jupiter's satellites, interior, magnetosphere and atmosphere. The power needs and solutions for the space exploration and lunar mobility program are discussed. Long term missions in space and on the lunar surface require high energy batteries. Rechargeable batteries for mobility systems and portable utility pallet are needed for successful exploration missions. Nanomaterial usage increases the energy density of the cells apart from increasing the power density. This slide presentation reviews the power requirements for the space exploration and the lunar surface mobility programs. It includes information about the specifications for high energy batteries and the power requirements for lunar rovers,

Roberts, William↗
DOE OSTI2024

MaterialsMap: A CALPHAD-based tool to design composition pathways through feasibility map for desired dissimilar materials, demonstrated with resistance spot welding joining of Ag-Al-Cu

Assembly of dissimilar metals can be achieved by different methods, for example, casting, welding, and additive manufacturing (AM). However, undesired phases formed in liquid-phase assembling processes due to solute segregation during solidification diminish mechanical and other properties of the processed parts. In the present work, an open-source software named MaterialsMap, has been developed based on the CALculation of Phase Diagrams (CALPHAD) approach. The primary objective of MaterialsMap is to facilitate the design of an optimal composition pathway for assembling dissimilar alloys with liquid-phases based on the formation of desired and undesired phases along the pathway. In MaterialsMap, equilibrium thermodynamic calculations are used to predict equilibrium phases formed at slow cooling rate, while Scheil-Gulliver simulations are employed to predict non-equilibrium phases formed during rapid cooling. Additionally, by combining these two simulations, MaterialsMap offers a thorough guide for understanding phase formation in various manufacturing processes, assisting users in making informed decisions during material selection and production. As a demonstration of this approach, a compositional pathway was designed from pure Al to pure Cu through Ag using MaterialsMap. The design was experimentally verified using resistance spot welding (RSW).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design of novel refractory equiatomic multi-principal elemental alloys based on Mo-Nb-Ti system for Gen IV reactor applications

Excellent irradiation damage resistance demonstrated by multi-principal elemental alloys (MPEAs) has sparked significant interest among researchers, prompting exploration into their vast compositional space, to validate their suitability for nuclear applications. A combined approach of thermodynamic and empirical parameters calculations alongside CALPHAD (CALculation of PHAse Diagrams) for phase formation predictions enable high-throughput material selection for sophisticated applications like nuclear, overcoming laborious and time-consuming experiments. Key thermodynamic and empirical parameters for eight novel equiatomic MPEAs, based on seven low thermal neutron cross section refractory elements, for predicting phase formation were calculated, and equilibrium and non-equilibrium simulations in CALPHAD were employed to comprehensively model the systems. Pseudo binary phase diagram simulations showed that Zr, V or equiatomic CrV additions to the base MoNbTi alloy (MoNbTi-Zr, MoNbTi-V and MoNbTi-CrV alloys) favor the formation of isomorphous body-centered cubic (BCC) phase at high temperatures, while Cr, Al, equiatomic ZrV, or equiatomic CrAl additions (MoNbTi-Cr, MoNbTi-Al, MoNbTi-ZrV or MoNbTi-CrAl alloys) limit the solubility of them. Equilibrium CALPHAD simulations at 750 oC were consistent with XRD results on MoNbTi, MoNbTiZr and MoNbTiCr alloys, and partially for others. Notably, elemental segregation observed in the backscattered electron (BSE) scanning electron microscopy (SEM) images of the alloys was accurately simulated through non-equilibrium Scheil solidification calculations in CALPHAD, further verified by experiments. The precipitation of TiCr2 Laves phase in Cr containing MoNbTiCr and MoNbTiCrAl was accurately predicted while discrepancies were noted in MoNbTiCrV. The equilibrium simulations also provided insights into phase compositions at specific temperatures offering a pathway for tailoring the desired microstructure and properties of these systems. Empirical parameters calculations successfully predicted random solid solution in the base MoNbTi alloy, and with an exception in MoNbTiV and MoNbTiAl, predicted intermetallic precipitation in the rest, especially, Laves phase precipitation in Cr containing alloys.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Pseudo-equilibrium theory for extrinsic doping control of the topological semimetal Cd 3 As 2

The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Out of time order correlation of the Hubbard model with random local disorder

The out-of-time-order correlator (OTOC) serves as a powerful tool for investigating quantum information spreading and chaos in complex systems. We present a method employing non-equilibrium dynamical mean-field theory and coherent potential approximation combined with diagrammatic perturbation on the Schwinger–Keldysh contour to calculate the OTOC for correlated fermionic systems subjected to both random disorder and electron interaction. Furthermore, our key finding is that random disorder enhances the OTOC decay in the Hubbard model for the metallic phase in the weakly interacting limit. However, the current limitation of our perturbative solver restricts the applicability to weak interaction regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗
DOE OSTI2026

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Non-equilibrium Dynamics of the Quark Gluon Plasma (Final Report)

During this funding period I continued my work on far-from-equilibrium formulations of hydrodynamics and bottomonium suppression in the quark-gluon plasma. I published 23 peer-reviewed articles and 9 proceedings contributions. The publications produced during this period have thus far received 705 citations. I delivered 22 talks, seminars, and colloquia about my work, most of which were invited. In addition to the work on hydrodynamics and heavy quarkonia, I also worked on high-order perturbative calculations of N=4 supersymmetric Yang-Mills thermodynamics and contributed to the efforts of the MUSES collaboration whose focus is developing a modular equation-of-state framework for use in modeling neutron stars and neutron star mergers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Nonlinear Dynamics and Nucleation Kinetics in Near-Critical Liquids

The objective of our study is to model the nonlinear behavior of a near-critical liquid following a rapid change of the temperature and/or other thermodynamic parameters (pressure, external electric or gravitational field). The thermodynamic critical point is manifested by large, strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, concentration in binary solutions) in the critical range of scales. The largest critical length scale is the correlation radius r(sub c). According to the scaling theory, r(sub c) increases as r(sub c) = r(sub 0)epsilon(exp -alpha) when the nondimensional distance epsilon = (T - T(sub c))/T(sub c) to the critical point decreases. The normal gravity alters the nature of correlated long-range fluctuations when one reaches epsilon approximately equal to 10(exp -5), and correspondingly the relaxation time, tau(r(sub c)), is approximately equal to 10(exp -3) seconds; this time is short when compared to the typical experimental time. Close to the critical point, a rapid, relatively small temperature change may perturb the thermodynamic equilibrium on many scales. The critical fluctuations have a hierarchical structure, and the relaxation involves many length and time scales. Above the critical point, in the one-phase region, we consider the relaxation of the liquid following a sudden temperature change that simultaneously violates the equilibrium on many scales. Below T(sub c), a non-equilibrium state may include a distribution of small scale phase droplets; we consider the relaxation of such a droplet following a temperature change that has made the phase of the matrix stable.

Patashinski, Alexander Z.↗
NASA NTRS2003

A Generalized Multi-Phase Framework for Modeling Cavitation in Cryogenic Fluids

A generalized multi-phase formulation for cavitation in fluids operating at temperatures elevated relative to their critical temperatures is presented. The thermal effects and the accompanying property variations due to phase change are modeled rigorously. Thermal equilibrium is assumed and fluid thermodynamic properties are specified along the saturation line using the NIST-12 databank. Fundamental changes in the physical characteristics of the cavity when thermal effects become pronounced are identified; the cavity becomes more porous, the interface less distinct, and has increased entrainment when temperature variations are present. Quantitative estimates of temperature and pressure depressions in both liquid nitrogen and liquid hydrogen were computed and compared with experimental data of Hord for hydrofoils. Excellent estimates of the leading edge temperature and pressure depression were obtained while the comparisons in the cavity closure region were reasonable. Liquid nitrogen cavities were consistently found to be in thermal equilibrium while liquid hydrogen cavities exhibited small, but distinct, non-equilibrium effects.

Dorney, Dan↗
NASA NTRS2010

An Overview of High-Resolution, Non-Dispersive, Imaging Spectrometers for High-Energy Photons

High-resolution x-ray spectroscopy has become a powerful tool for studying the evolving universe. The grating spectrometers on the XMM and Chandra satellites initiated a new era in x-ray astronomy. Despite their successes, there is still need for instrumentation that can provide higher spectral resolution with high throughput in the Fe-K band and for extended sources. What is needed is a non-dispersive imaging spectrometer - essentially a 14-bit x-ray color camera. And a requirement for a nondispersive spectrometer designed to provide eV-scale spectral resolution is a temperature below 0.1 K. The required spectral resolution and the constraints of thermodynamics and engineering dictate the temperature regime nearly independently of the details of the sensor or the read-out technology. Low-temperature spectrometers can be divided into two classes - - equilibrium and non-equilibrium. In the equilibrium devices, or calorimeters, the energy is deposited in an isolated thermal mass and the resulting increase in temperature is measured. In the non-equilibrium devices, the absorbed energy produces quantized excitations that are counted to determine the energy. The two approaches have different strong points, and within each class a variety of optimizations have been pursued. I will present the basic fundamentals of operation and the details of the most successful device designs to date. I will also discuss how the measurement priorities (resolution, energy band, count rate) influence the optimal choice of detector technology.

Kilbourne, Caroline↗
NASA NTRS2011

Non-Detection of L-Band Line Emission from the Exoplanet HD 189733b

We attempt to confirm bright non-local thermodynamic equilibrium (non-LTE) emission from the exoplanet HD 189733b at 3.25 microns, as recently reported by Swain et al. based on observations at low spectral resolving power (lambda/delta lambda approximately equals 30). Non-LTE emission lines from gas in an exoplanet atmosphere will not be significantly broadened by collisions, so the measured emission intensity per resolution element must be substantially brighter when observed at high spectral resolving power. We observed the planet before, during, and after a secondary eclipse event at a resolving power lambda/delta lambda equals 27,000 using the NIRSPEC spectrometer on the Keck II telescope. Our spectra cover a spectral window near the peak found by Swain et al., and we compare emission cases that could account for the magnitude and wavelength dependence of the Swain et al. result with our final spectral residuals. To model the expected line emission, we use a general non-equilibrium formulation to synthesize emission features from all plausible molecules that emit in this spectral region. In every case, we detect no line emission to a high degree of confidence. After considering possible explanations for the Swain et at., results and the disparity with our own data, we conclude that an astrophysical source for the putative non-LTE emission is unlikely. We note that the wavelength dependence of the signal seen by Swain et al. closely matches the 2(sub 1/2) band of water vapor at 300 K, and we suggest that an imperfect correction for telluric water is the source of the feature claimed by Swain et al.

Mandell, Avi M.↗