Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium potential”

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

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

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

Structure of the electrical double layer at the ice–water interface

The surface of ice in contact with water contains sites that undergo deprotonation and protonation and can act as adsorption sites for aqueous ions. Therefore, an electrical double layer should form at this interface and existing models for describing the electrical double layer at metal oxide–water interfaces should be able to be modified to describe the surface charge, surface potential, and ionic occupancy at the ice–water interface. I used a surface complexation model along with literature measurements of the zeta potential of ice in brines of various strength and pH to constrain equilibrium constants. I then made predictions of ion site occupancy, surface charge density, and partitioning of counterions between the Stern and diffuse layers. The equilibrium constant for cation adsorption is more than 5 orders of magnitude larger than the other constants, indicating that this reaction dominates even at low salinity. Deprotonated OH sites are predicted to be slightly more abundant than dangling O sites, consistent with previous work. Surface charge densities are on the order of ±0.001 C/m2 and are always negative at the moderate pH values of interest to atmospheric and geophysical applications (6–9). In this pH range, over 99% of the counterions are contained in the Stern layer. This suggests that diffuse layer polarization will not occur because the ionic concentrations in the diffuse layer are nearly identical to those in the bulk electrolyte and that electrical conduction and polarization in the Stern layer will be negligible due to reduced ion mobility.

Daigle, Hugh (ORCID:0000000260628321)↗
DOE OSTI2025

Overcoming time and complexity limitations in molecular dynamics investigations of equilibrium melting

Abstract A hybrid Monte-Carlo molecular-dynamics method for determining solidus and liquidus compositions in multicomponent systems is presented that overcomes both the time limitations in conventional molecular dynamics that prevent the evolution of distinct solid and liquid compositions via diffusion and the complexity challenge that prevents use of thermodynamic assessment in systems of many components. This hybrid method is validated in the Cu–Ni system against an independent assessment of solidus and liquidus compositions based on the regular solution model. Strategies for efficient mapping of different phase diagrams, based on the thermodynamic parameter T 0 , the temperature at which two phases of the composition X 0 have equal free energies, are presented and then demonstrated for the copper-nickel fully miscible system and the gold–silicon eutectic system. A calculation of the solidus and liquidus sampled during the equilibrium melting of equiatomic CrMnFeCoNi is performed, indicating that this method has potential to be extended to the study of many component alloys.

Au-Si↗
DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

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

Spin-polarized atomic nitrogen and the 7Sigma + u state of N2

The first self-consistent field (SCF) calculation of the repulsive 7Sigma + u state of molecular nitrogen is presented. This calculation is used with attractive dispersion to provide a model potential for electron spin-polarized atomic nitrogen N appropriate in the range of 2.5-20A(O). Potential parameters obtained from the calculation are used to provide estimates of equilibrium thermodynamic properties of spin-polarized N in the quantum theorem of corresponding states framework. Consideration is given to increased stability of spin-polarized N with respect to single electron spin-flip processes, as compared to spin-polarized H.

Ferrante, R. F.↗
NASA NTRS1987

Determination of the masses of elliptical galaxies and clusters of galaxies with AXAF

The procedures to be used in the determination of the mass and mass distribution of elliptical galaxies and clusters of galaxies, using spatially resolved X-ray spectra of galaxies and clusters of galaxies to be obtained by the imaging spectrometers on AXAF, are described. The procedures utilize the unique AXAF features of large band width, large field of view, high angular resolution, long exposure times, and the availability of sophisticated AXAF instrumentation. The techniques described involve the measurement of the surface brightness and temperature profiles of the hot gas trapped in the potential wells of these systems. Then, by using the equation of hydrostatic equilibrium, which relates the total mass inside a given radius to the size, absolute temperature, density gradient, and temperature profile, it will be possible to infer the mass distribution in these objects to large physical scales.

Mushotzky, Richard F.↗
NASA NTRS1989

Properties of possible polar rings around Neptune

The near-polar-position rings whose possible existence and character around Neptune are presently considered correspond to the equilibrium configuration of a test-particle orbital plane in the gravitational potential of Neptune and Triton. At these equilibria, of which two whose invariable-plane nodal line is approximately perpendicular to the nodal line of Triton on the same plane are dynamically stable, the orbit of the test particle maintains a fixed orientation with respect to the invariable plane and the node of triton on that plane. The stable equilibria are noted to correspond to a maximum of the total energy for a given semimajor axis; their stability is found to remain unperturbed despite collisions among particles.

Borderies, Nicole↗
NASA NTRS1992

Tidal deceleration of the moon's mean motion

The secular change in the mean motion of the moon, n, caused by the tidal dissipation in the ocean and solid earth is due primarily to the effect of the diurnal and semidiurnal tides. The long-period ocean tides produce an increase in n, but the effects are only 1 percent of the diurnal and semidiurnal ocean tides. In this investigation, expressions for these effects are obtained by developing the tidal potential in the ecliptic reference system. The computation of the amplitude of equilibrium tide and the phase corrections is also discussed. The averaged tidal deceleration of the moon's mean motion, n, from the most recent satellite ocean tide solutions is -25.25 +/- 0.4 arcseconds/sq century. The value for n inferred from the satellite-determined ocean-tide solution is in good agreement with the value obtained from the analysis of 20 years of lunar laser-ranging observations.

Cheng, M. K.↗
NASA NTRS1997

Kinetics and Thermochemistry of the Br((sup 2)P3/2) + NO2 Association Reaction

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the Br((sup 2)P3/2) + NO2 association reaction as a function of temperature (259-432 K) pressure (12.5 - 700 Torr), and buffer gas identity (He, Ar, H2, N2, CO2, CF4, SF6). The reaction is found to be in the falloff regime between third and second order over the entire range of conditions investigated. At temperatures below 350 K, the association reaction is found to be irreversible on the time scale of the experiment (approximately 30 ms). At higher temperatures reversible addition is observed, allowing equilibrium constants for BrNO2 formation and dissociation to be determined. Second- and third-law analyses of the equilibrium data are in only fair agreement and lead to the following thermochemical parameters for the association reaction: Delta-H(298) = 19.6 +/- 1.7 kcal/mol, Delta-H(0) = -18.6 +/- 2.0 kcal/mol, Delta-S(298) = 29.3 +/- 4.2 cal/mol/K, Delta-H(sub f)(sub 298)(BrNO2) = 17.0 +/-1.8 kcal/mol(uncertainties are 2 sigma estimates of absolute accuracy). The value for Delta-H(0) determined in this study has been employed to calculate k(sub 0)(sup SC), the low-pressure third-order rate coefficient in the strong collision limit, by using the method of Troe; calculated values of k(sub 0)(sup SC) are inconsistent with experimental results unless Delta-H(0) is assigned a value near the lower limit derived from analysis of the high-temperature approach to equilibrium data, i.e. delta-H(0) approximately equals -16.6 kcal/mol. A potential source of systematic error in the calculation of both k(sub 0)(sup SC) and the absolute entropy of BrNO2 results from the complete lack of knowledge of the energies and degeneracies of the electronic states of BrNO3. The procedure developed by Troe and co-workers has been employed to extrapolate experimental falloff curves to the low- and high-pressure limits. Derived values for k(sub 0)(M,298K) in units of 10(exp -31) cm(exp 6)/sq molecule/s range from 2.75 for M = He to 6.54 for M = CO2; 2 sigma uncertainties are estimated to be +/- 20%. Values for k(sub 0)(N2,T) in units of 10(exp -31) cm(exp 6)/sq molecule/s are 5.73 at 259 K, 4.61 at 298 K, and 3.21 at 346 K; the observed temperature dependence for k(sub 0)(N2,T) is consistent with the theoretical temperature dependence for Beta(sub c)k(sub 0)(sup SC). Values for k(sub infinity)(T) in units of 10(exp -11) cu cm/molecule/s are 2.86 at 259 K, 3.22 at 298 K, and 3.73 at 346 K; 2 sigma uncertainties are estimated to be a factor of 2. Approximate falloff parameters in a convenient format for atmospheric modeling are also derived.

Kreutter, K. D.↗
NASA NTRS2003

Searching for Alien Life Having Unearthly Biochemistry

The search for alien life in the solar system should include exploring unearth-like environments for life having an unearthly biochemistry. We expect alien life to conform to the same basic chemical and ecological constraints as terrestrial life, since inorganic chemistry and the laws of ecosystems appear to be universal. Astrobiologists usually assume alien life will use familiar terrestrial biochemistry and therefore hope to find alien life by searching near water or by supplying hydrocarbons. The assumption that alien life is likely to be based on carbon and water is traditional and plausible. It justifies high priority for missions to search for alien life on Mars and Europa, but it unduly restricts the search for alien life. Terrestrial carbon-water biochemistry is not possible on most of the bodies of our solar system, but all alien life is not necessarily based on terrestrial biochemistry. If alien life has a separate origin from Earth life, and if can survive in an environment extremely different from Earth's, then alien life may have unearthly biochemistry. There may be other solvents than water that support alien life and other elements than carbon that form complex life enabling chain molecules. Rather than making the exploration-restricting assumption that all life requires carbon, water, and terrestrial biochemistry, we should make the exploration-friendly assumption that indigenous, environmentally adapted, alien life forms might flourish using unearthly biochemistry in many places in the solar system. Alien life might be found wherever there is free energy and a physical/chemical system capable of using that energy to build living structures. Alien life may be discovered by the detection of some general non-equilibrium chemistry rather than of terrestrial biochemistry. We should explore all the potential abodes of life in the solar system, including those where life based on terrestrial biochemistry can not exist.

Jones, Harry↗
NASA NTRS2015

Open-Source Conceptual Sizing Models for the Hyperloop Passenger Pod

Hyperloop is a new mode of transportation proposed as an alternative to California's high speed rail project, with the intended benefits of higher performance at lower overall costs. It consists of a passenger pod traveling through a tube under a light vacuum and suspended on air bearings. The pod travels up to transonic speeds resulting in a 35 minute travel time between the intended route from Los Angeles and San Francisco. Of the two variants outlined, the smaller system includes a 1.1 meter tall passenger capsule traveling through a 2.2 meter tube at 700 miles per hour. The passenger pod features water-based heat exchangers as well as an on-board compression system that reduces the aerodynamic drag as it moves through the tube. Although the original proposal looks very promising, it assumes that tube and pod dimensions are independently sizable without fully acknowledging the constraints of the compressor system on the pod geometry. This work focuses on the aerodynamic and thermodynamic interactions between the two largest systems; the tube and the pod. Using open-source toolsets, a new sizing method is developed based on one-dimensional thermodynamic relationships that accounts for the strong interactions between these sub-systems. These additional considerations require a tube nearly twice the size originally considered and limit the maximum pod travel speed to about 620 miles per hour. Although the results indicate that Hyperloop will need to be larger and slightly slower than originally intended, the estimated travel time only increases by approximately five minutes, so the overall performance is not dramatically affected. In addition, the proposed on-board heat exchanger is not an ideal solution to achieve reasonable equilibrium air temperatures within the tube. Removal of this subsystem represents a potential reduction in weight, energy requirements and complexity of the pod. In light of these finding, the core concept still remains a compelling possibility, although additional engineering and economic analyses are markedly necessary before a more complete design can be developed.

design analysis↗
NASA NTRS2017

A Library of ATMO Forward Model Transmission Spectra for Hot Jupiter Exoplanets

We present a grid of forward model transmission spectra, adopting an isothermal temperature-pressure profile, alongside corresponding equilibrium chemical abundances for 117 observationally significant hot exoplanets (equilibrium temperatures of 547-2710 K). This model grid has been developed using a 1D radiative-convective-chemical equilibrium model termed ATMO, with up-to-date high-temperature opacities. We present an interpretation of observations of 10 exoplanets, including best-fitting parameters and X(exp 2) maps. In agreement with previous works, we find a continuum from clear to hazy/cloudy atmospheres for this sample of hot Jupiters. The data for all the 10 planets are consistent with subsolar to solar C/O ratio, 0.005 to 10 times solar metallicity and water rather than methane-dominated infrared spectra. We then explore the range of simulated atmospheric spectra for different exoplanets, based on characteristics such as temperature, metallicity, C/O ratio, haziness and cloudiness. We find a transition value for the metallicity between 10 and 50 times solar, which leads to substantial changes in the transmission spectra. We also find a transition value of C/O ratio, from water to carbon species dominated infrared spectra, as found by previous works, revealing a temperature dependence of this transition point ranging from approximately 0.56 to approximately 1-1.3 for equilibrium temperatures from approximately 900 to approximately 2600 K. We highlight the potential of the spectral features of HCN and C2H2 to constrain the metallicities and C/O ratios of planets, using James Webb Space Telescope (JWST) observations. Finally, our entire grid (approximately 460 000 simulations) is publicly available and can be used directly with the JWST simulator PandExo for planning observations.

techniques: spectroscopic – planets and satellit↗
NASA NTRS1990

Application of the embedded atom method to Pb and Be

We have derived the embedding energy functional and two-body potential of the embedded atom method using decreasing exponentials for both the electron density and the two body potential. The embedding function was obtained from the equation of state given by Rose et al. (1984). Because of the form of the embedding function, the equilibrium lattice constant, cohesive energy, and bulk modulus are automatically satisfied. The two parameters phi(e) and gamma of the two-body potential were determined by fitting to shear modulus and the single vacancy formation energy. Contributions of up to the third nearest neighbors were included in the evaluation of the charge density rho and the two-body potential phi. The stability and anisotropy of each structure were estimated and compared with the available experimental data.

Karimi, M.↗
NASA NTRS1977

Radio and EUV observations of a coronal hole

Current ideas and modeling assumptions concerning the structure of the transition region and corona are summarized. Observations of a specific coronal hole are reported which were made with the EUV spectroheliometer aboard Skylab and with high-resolution radio telescopes in Australia and West Germany. Theoretical EUV line intensities and radio brightness temperatures are calculated, and an attempt is made to analyze the observations on the basis of standard one-dimensional plane-parallel hydrostatic-equilibrium models for the transition region and corona. It is found that any given standard model predicts either higher radio brightnesses or lower EUV line intensities than are observed. Several potential problem areas are examined in order to resolve this conflict. It is concluded that a model involving departures from hydrostatic equilibrium, varying conductive flux, mass outflow, as well as effects of thermal diffusion and nonequilibrium ionization should be investigated in more detail.

Dulk, G. A.↗
DOE OSTI2020

General Method for Determining Redox Potentials without Electrolyte

A novel method to determine redox potentials without electrolyte is presented here. The method is based on a new ability to determine the dissociation constant, $K°_d$, for ion pairs formed between any radical anion and any inert electrolyte counterion. These dissociation constants can be used to determine relative shifts of redox potential as a function of electrolyte concentration, connecting referenced potentials determined with electrochemistry (with 0.1 M electrolyte) to electrolyte-free values. Pulse radiolysis created radical anions enabling determination of equilibrium constants for electron transfer between anions of donor and acceptor molecules as a function of electrolyte concentration in THF. The measurements determined “composite equilibrium constants”, $K_{eqC}$, which contain information about the dissociation constant for the electrolyte cations, X + , with the radical anions of both the donor, $K°_d$(D –• ,X + ) and the acceptor, $K°_d$(A –• ,X + ). Dissociation constants were obtained for a selection of radical anions with tetrabutylammonium (TBA + ). The electrolyte was found to shift the reduction potentials of small molecules 1-methylpyrene and trans-stilbene by close to +130 mV whereas oligo-fluorenes and polyfluorenes experienced shifts of only (+25 ± 6) mV due to charge delocalization weakening the ion pair. These shifts for reduction of aromatic hydrocarbon molecules are smaller than shifts of +232 and +451 mV seen previously for benzophenone radical anion with TBA + and Na + respectively where the charge on the radical anion is localized largely on one C=O bond, thus forming a more tightly bound ion pair.

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

Polymorphic Metastability in Colloidal Semiconductor Nanocrystals

Metastable polymorphs of inorganic solids often possess material properties not present in the corresponding thermodynamic polymorphs, making them targets for the development of new functional materials. In contrast with isolating metastable bulk materials, syntheses of metastable polymorphs on the nanoscale are aided by fast non-equilibrium reaction kinetics and the favorable thermodynamic influence of surface energies, giving rise to greater ease of access to metastable high-temperature polymorphs and, in some cases, new polymorphs that do not exist in the bulk. The syntheses of metastable semiconductor nanocrystals are of interest for their potentially unique optoelectronic and physicochemical properties. However, in many material systems, synthesizing nanocrystalline products away from thermodynamic equilibrium in a predictable manner remains an outstanding challenge. This review outlines direct synthetic methodologies that have been developed to enable control over the nucleation and growth of metastable polymorphs of semiconductor nanocrystals by tailoring reaction conditions, precursor kinetics, ligand and surface effects, and other synthetic levers. The case studies reviewed herein expound on the direct syntheses of metastable ZnSe, Cu 2 SnSe 3 , CuInSe 2 , Ag 2 Se, and AgInSe 2 nanocrystals, and although there remain numerous examples of metastable nanocrystal syntheses outside of these metal chalcogenide systems, the concepts discussed are of general utility to the field of metastable nanocrystal syntheses as a whole. Explicit examples in which new functional properties are afforded by metastable polymorphs of the aforementioned material systems are presented within the context of applications for solar cells, photonics, and optical sensing. Finally, the factors that affect the kinetic persistence of metastable nanocrystalline polymorphs are discussed at length for these material systems.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Collapse, equilibrium, and fragmentation of rotating, adiabatic clouds

Numerical calculations of the collapse of adiabatic clouds from uniform density and rotation initial conditions show that when restricted to axisymmetry, the clouds form either near-equilibrium spheroids or rings. Rings form in the collapse of low thermal energy clouds and have a ratio of rotational kinetic energy to the absolute value of gravitational potential energy greater than approximately 0.43. When the axisymmetric constraint is removed and an initial m = 2 density variation is introduced, clouds either collapse to form near-equilibrium ellipsoids or else fragment into binary systems through a bar phase. Ellipsoids form in the collapse of high thermal energy clouds and have a rotational kinetic energy/absolute value of gravitational potential energy ratio less than approximately 0.27. The results are consistent with the critical values of the rotational kinetic energy/absolute value of gravitational potential energy ratio for instabilities in Maclaurin spheroids, and suggest that protostellar clouds may undergo a dynamic fragmentation in the nonisothermal collapse regime.

Boss, A. P.↗