SEARCH · Engineering Papers
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.
The thermodynamic technique of Landau and Lifshitz originally developed for inter-species diffusion in a binary neutral gas mixture is extended to a quasi-neutral plasma with two ion species. It is shown that, while baro- and electro-diffusion coefficients depend on the choice of the thermodynamic system, prediction for the total diffusive mass flux is invariant.
Thermodynamic and diffusion properties of water + methane mixtures in a single liquid phase are studied using NVT molecular dynamics. An extensive comparison is reported for the thermal pressure coefficient, compressibilities, expansion coefficients, heat capacities, Joule-Thomson coefficient, zero frequency speed of sound, and diffusion coefficient at methane concentrations up to 15% in the temperature range of 298–650 K. The simulations reveal a complex concentration dependence of the thermodynamic properties of water + methane mixtures. The compressibilities, heat capacities, and diffusion coefficients decrease with increasing methane concentration, whereas values of the thermal expansion coefficients and speed of sound increase. Increasing methane concentration considerably retards the self-diffusion of both water and methane in the mixture. These effects are caused by changes in hydrogen bond network, solvation shell structure, and dynamics of water molecules induced by the solvation of methane at constant volume conditions.
The well-known classical nucleation theory (CNT) for the free energy barrier towards formation of a nucleus of critical size of the new stable phase within the parent metastable phase fails to take into account the influence of other metastable phases having density/order intermediate between the parent metastable phase and the final stable phase. This lacuna can be more serious than capillary approximation or spherical shape assumption made in CNT. This issue is particularly significant in ice nucleation because liquid water shows rich phase diagram consisting of two (high and low density) liquid phases in supercooled state. The explanations of thermodynamic and dynamic anomalies of supercooled water often invoke the possible influence of a liquid-liquid transition between two metastable liquid phases. To investigate both the role of thermodynamic anomalies and presence of distinct metastable liquid phases in supercooled water on ice nucleation, we employ density functional theoretical approach to find nucleation free energy barrier in different regions of phase diagram. The theory makes a number of striking predictions, such as a dramatic lowering of nucleation barrier due to presence of a metastable intermediate phase and crossover in the dependence of free energy barrier on temperature near liquid-liquid critical point. These predictions can be tested by computer simulations as well as by controlled experiments.
First-principles calculations were employed to investigate hypothetical complex hydrides (M,M{sup ′}){sub 4}FeH{sub 8} (M = Na, Li; M{sup ′}=Mg, Zn, Y, Al). Besides complex anion [FeH{sub 6}]{sup 4–}, these materials contain two H{sup –} anions, which raise the total anionic charge state from tetravalent to hexavalent, and thereby significantly increasing the number of combinations of countercations. We have determined that similar to complex hydrides (M,M{sup ′}){sub 2}FeH{sub 6} containing only [FeH{sub 6}]{sup 4–}, the thermodynamic stability is tuned by the average cation electronegativity. Thus, the chemical flexibility provided by incorporating H{sup –} enhances the tunability of thermodynamic stability, which will be beneficial in obtaining optimal stability for hydrogen storage materials.
We report a theoretical study of stochastic processes modeling the growth of first-order Markov copolymers, as well as the reversed reaction of depolymerization. These processes are ruled by kinetic equations describing both the attachment and detachment of monomers. Exact solutions are obtained for these kinetic equations in the steady regimes of multicomponent copolymerization and depolymerization. Thermodynamic equilibrium is identified as the state at which the growth velocity is vanishing on average and where detailed balance is satisfied. Away from equilibrium, the analytical expression of the thermodynamic entropy production is deduced in terms of the Shannon disorder per monomer in the copolymer sequence. The Mayo-Lewis equation is recovered in the fully irreversible growth regime. The theory also applies to Bernoullian chains in the case where the attachment and detachment rates only depend on the reacting monomer.
A novel thermodynamic integration (TI) scheme is presented to compute the crystal-liquid interfacial free energy (γ{sub cl}) from molecular dynamics simulation. The scheme is applied to a Lennard-Jones system. By using extremely short-ranged and impenetrable Gaussian flat walls to confine the liquid and crystal phases, we overcome hysteresis problems of previous TI schemes that stem from the translational movement of the crystal-liquid interface. Our technique is applied to compute γ{sub cl} for the (100), (110), and (111) orientation of the crystalline phase at three temperatures under coexistence conditions. For one case, namely, the (100) interface at the temperature T = 1.0 (in reduced units), we demonstrate that finite-size scaling in the framework of capillary wave theory can be used to estimate γ{sub cl} in the thermodynamic limit. Thereby, we show that our TI scheme is not associated with the suppression of capillary wave fluctuations.
A theory to describe structural fluctuations of protein induced by thermodynamic perturbations, pressure, temperature, and denaturant, is proposed. The theory is formulated based on the three methods in the statistical mechanics: the generalized Langevin theory, the linear response theory, and the three dimensional interaction site model (3D-RISM) theory. The theory clarifies how the change in thermodynamic conditions, or a macroscopic perturbation, induces the conformational fluctuation, which is a microscopic property. The theoretical results are applied, on the conceptual basis, to explain the experimental finding by Akasaka et al., concerning the NMR experiment which states that the conformational change induced by pressure corresponds to structural fluctuations occurring in the ambient condition. A method to evaluate the structural fluctuation induced by pressure is also suggested by means of the 3D-RISM and the site-site Kirkwood-Buff theories.
General methods for combining interactions between particles characterised by non-identical intermolecular potentials are investigated. The combination methods are tested by performing molecular dynamics simulations to determine the pressure, energy, isochoric and isobaric heat capacities, thermal expansion coefficient, isothermal compressibility, Joule-Thomson coefficient, and speed of sound of 10-5 + 12-6 Mie potential binary mixtures. In addition to the two non-identical Mie potentials, mixtures are also studied with non-identical intermolecular parameters. The combination methods are compared with results obtained by simply averaging the Mie exponents. When either the energy or size parameters are non-identical, very significant differences emerge in the thermodynamic properties predicted by the alternative combination methods. The isobaric heat capacity is the thermodynamic property that is most affected by the relative magnitude of the intermolecular potential parameters and the method for combining non-identical potentials. Either the arithmetic or geometric combination of potentials provides a simple and effective way of performing simulations involving mixtures of components characterised by non-identical intermolecular potentials, which is independent of their functional form.
A rigorous thermodynamic formulation of the geometric model for heterogeneous nucleation including line tension effect is missing till date due to the associated mathematical hurdles. In this work, we develop a novel thermodynamic formulation based on Classical Nucleation Theory (CNT), which is supposed to illustrate a systematic and a more plausible analysis for the heterogeneous nucleation on a planar surface including the line tension effect. The appreciable range of the critical microscopic contact angle (θ{sub c}), obtained from the generalized Young’s equation and the stability analysis, is θ{sub ∞} < θ{sub c} < θ′ for positive line tension and is θ{sub M} < θ{sub c} < θ{sub ∞} for negative line tension. θ{sub ∞} is the macroscopic contact angle, θ′ is the contact angle for which the Helmholtz free energy has the minimum value for the positive line tension, and θ{sub M} is the local minima of the nondimensional line tension effect for the negative line tension. The shape factor f, which is basically the dimensionless critical free energy barrier, becomes higher for lower values of θ{sub ∞} and higher values of θ{sub c} for positive line tension. The combined effect due to the presence of the triple line and the interfacial areas (f{sup L} + f{sup S}) in shape factor is always within (0, 3.2), resulting f in the range of (0, 1.7) for positive line tension. A formerly presumed appreciable range for θ{sub c}(0 < θ{sub c} < θ{sub ∞}) is found not to be true when the effect of negative line tension is considered for CNT. Estimation based on the property values of some real fluids confirms the relevance of the present analysis.
The compressibility, elastic anisotropy, and thermodynamic properties of the recently proposed tetragonal B{sub 4}C{sub 4} (t-B{sub 4}C{sub 4}) are investigated under high temperature and high pressure by using of first-principles calculations method. The elastic constants, bulk modulus, shear modulus, Young’s modulus, Vickers hardness, Pugh’s modulus ratio, and Poisson’s ratio for t-B{sub 4}C{sub 4} under various pressures are systematically explored, the obtained results indicate that t-B{sub 4}C{sub 4} is a stiffer material. The elastic anisotropies of t-B{sub 4}C{sub 4} are discussed in detail under pressure from 0 GPa to 100 GPa. The thermodynamic properties of t-B{sub 4}C{sub 4}, such as Debye temperature, heat capacity, and thermal expansion coefficient are investigated by the quasi-harmonic Debye model.
The structural, mechanical, electronic, and thermodynamic properties of Mg{sub 2}Sr with C14-type structure under pressures ranging from 0 to 40 GPa have been systematically studied within the framework of density functional theory in this work. The results at zero pressure were in good agreement with the available theoretical and experimental values. The pressure dependence of structure and elastic constants, elastic anisotropy, Cauchy pressure, melting points, and hardness was successfully calculated and discussed. In addition, the electronic density of states (DOSs) under various pressures were investigated. Debye temperature and the dependences of thermodynamic properties on temperature and pressure were also discussed in the present paper.
Thermodynamical properties like Entropy (S), Internal energy (E) and Helmholtz free energy (F) of liquid lanthanides using a variation principle based on the Gibbs-Bogoliubuv (GB) inequality with Percus Yevick hard sphere reference system have been reported in the present investigation. To describe electron-ion interaction we have used our newly constructed parameter free model potential along with Sarkar et al. local field correction function. Lastly, we conclude that our newly constructed model potential is capable to explain the thermodynamical properties of liquid lanthanides.
In this paper, we define and experimentally verify thermodynamic characteristics of the liquid-glass transition, taking into account a kinetic origin of the process. Using the density scaling law and the four-point measure of the dynamic heterogeneity of molecular dynamics of glass forming liquids, we investigate contributions of enthalpy, temperature, and density fluctuations to spatially heterogeneous molecular dynamics at the liquid-glass transition, finding an equation for the pressure coefficient of the glass transition temperature, dTg/dp. This equation combined with our previous formula for dTg/dp, derived solely from the density scaling criterion, implies a relationship among thermodynamic coefficients at Tg. Since this relationship and both the equations for dTg/dp are very well validated using experimental data at Tg, they are promising alternatives to the classical Prigogine-Defay ratio and both the Ehrenfest equations in case of the liquid-glass transition.
Here, a traditional single-junction solar panel cannot harness ground-scattered light (albedo reflectance, R A ), and also suffers from the fundamental sub-band-gap and the thermalization losses. In this paper, we explain how a “bifacial tandem” panel would dramatically reduce these losses, with corresponding improvement in thermodynamic performance. Our study predicts (i) the optimum combination of the band-gaps, empirically given by E opt g(t) ≈ E opt g(b) (2+R A )/3+(1–R A ) and the (ii) corresponding optimum normalized output power given by η* T(opt) ≈ R A (2η SJ(opt) ) + (1–R A )η DJ(opt) . Empirically, η* T(opt) interpolates between the thermodynamic efficiency limit of classical double-junction tandem cell (η DJ ) and twice that of a single-junction cell (η SJ ). We conclude by explaining how the fundamental loss mechanisms evolve with R A in a bifacial tandem cell.
The thermodynamic properties, including the density, volume expansion coefficient, ratio of specific heat to emissivity of intermetallic Ni{sub 7}Zr{sub 2} alloy, have been measured using the non-contact electrostatic levitation technique. These properties vary linearly with temperature at solid and liquid states, even down to the obtained maximum undercooling of 317 K. The enthalpy, glass transition, diffusion coefficient, shear viscosity, and surface tension were obtained by using molecular dynamics simulations. Ni{sub 7}Zr{sub 2} has a relatively poor glass forming ability, and the glass transition temperature is determined as 1026 K. The inter-diffusivity of Ni{sub 7}Zr{sub 2} alloy fitted by Vogel–Fulcher–Tammann law yields a fragility parameter of 8.49, which indicates the fragile nature of this alloy. Due to the competition of increased thermodynamic driving force and decreased atomic diffusion, the dendrite growth velocity of Ni{sub 7}Zr{sub 2} compound exhibits double-exponential relationship to the undercooling. The maximum growth velocity is predicted to be 0.45 m s{sup −1} at the undercooling of 335 K. Theoretical analysis reveals that the dendrite growth is a diffusion-controlled process and the atomic diffusion speed is only 2.0 m s{sup −1}.
We have set up and manually curated a dataset containing experimental information on the impact of amino acid substitutions in a protein on its thermal stability. It consists of a repository of experimentally measured melting temperatures (T{sub m}) and their changes upon point mutations (ΔT{sub m}) for proteins having a well-resolved x-ray structure. This high-quality dataset is designed for being used for the training or benchmarking of in silico thermal stability prediction methods. It also reports other experimentally measured thermodynamic quantities when available, i.e., the folding enthalpy (ΔH) and heat capacity (ΔC{sub P}) of the wild type proteins and their changes upon mutations (ΔΔH and ΔΔC{sub P}), as well as the change in folding free energy (ΔΔG) at a reference temperature. These data are analyzed in view of improving our insights into the correlation between thermal and thermodynamic stabilities, the asymmetry between the number of stabilizing and destabilizing mutations, and the difference in stabilization potential of thermostable versus mesostable proteins.
In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.
Coarse-grained models are a core computational tool in theoretical chemistry and biophysics. A judicious choice of a coarse-grained model can yield physical insights by isolating the essential degrees of freedom that dictate the thermodynamic properties of a complex, condensed-phase system. The reduced complexity of the model typically leads to lower computational costs and more efficient sampling compared with atomistic models. Designing “good” coarse-grained models is an art. Generally, the mapping from fine-grained configurations to coarse-grained configurations itself is not optimized in any way; instead, the energy function associated with the mapped configurations is. In this work, we explore the consequences of optimizing the coarse-grained representation alongside its potential energy function. We use a graph machine learning framework to embed atomic configurations into a low-dimensional space to produce efficient representations of the original molecular system. Because the representation we obtain is no longer directly interpretable as a real-space representation of the atomic coordinates, we also introduce an inversion process and an associated thermodynamic consistency relation that allows us to rigorously sample fine-grained configurations conditioned on the coarse-grained sampling. We illustrate that this technique is robust, recovering the first two moments of the distribution of several observables in proteins such as chignolin and alanine dipeptide.