Local equilibrium in liquid phase shock waves.

We have assessed the assumption of local thermodynamic equilibrium in a shock wave by comparing local thermodynamic data generated with nonequilibrium molecular dynamics (NEMD) simulations with results from corresponding equilibrium simulations. The shock had a Mach number approximately equal to 2 in a Lennard-Jones spline liquid. We found that the local equilibrium assumption holds perfectly well behind the wave front, and is a very good approximation in the front itself. This was supported by calculations of the excess entropy production in the shock front with four different methods that use the local equilibrium assumption in different ways. Two of the methods assume local equilibrium between excess thermodynamic variables by treating the shock as an interface in Gibbs's sense. The other two methods are based on the local equilibrium assumption in a continuous description of the shock front. We show for the shock studied in this work that all four methods give excess entropy productions that are in excellent agreement, with an average variance of 3.5% for the nonequilibrium molecular dynamics (NEMD) simulations. In addition, we solved the Navier-Stokes (N-S) equations numerically for the same shock wave using an equilibrium equation of state (EoS) based on a recently developed perturbation theory. The results for the density, pressure, and temperature profiles agree well with the profiles from the NEMD simulations. For instance, the shock waves generated in the two simulations travel with almost the same speed; the average absolute Mach-number deviation of the N-S simulations relative to NEMD is 2.6% in the investigated time interval.

Seebeck, Peltier, and Soret effects: On different formalisms for transport equations in thermogalvanic cells.

Thermogalvanic cells convert waste heat directly to electric work. There is an abundance of waste heat in the world and thermogalvanic cells may be underused. We discuss theoretical tools that can help us understand and therefore improve on cell performance. One theory is able to describe all aspects of the energy conversion: nonequilibrium thermodynamics. We recommend to use the theory with operationally defined, independent variables, as others have done before. These describe well-defined experiments. Three invariance criteria serve as a basis for any description: of local electroneutrality, entropy production invariance, and emf's independence of the frame of reference. Alternative formalisms, using different sets of variables, start with ionic or neutral components. We show that the heat flux is not the same in the two formalisms and derive a new relationship between the heat fluxes. The heat flux enters the definition of the Peltier coefficient and is essential for the understanding of the Peltier heat at the electrode interfaces and of the Seebeck coefficient of the cell. The Soret effect can occur independently of any Seebeck effect, but the Seebeck effect will be affected by the presence of a Soret effect. Common misunderstandings are pointed out. Peltier coefficients are needed for the interpretation and design of experiments.

Legendre-Fenchel transforms capture layering transitions in porous media.

We have investigated the state of a nanoconfined fluid in a slit pore in the canonical and isobaric ensembles. The systems were simulated with molecular dynamics simulations. The fluid has a transition to a close-packed structure when the height of the slit approaches the particle diameter. The Helmholtz energy is a non-convex function of the slit height if the number of particles does not exceed that of one monolayer. As a consequence, the Legendre transform cannot be applied to obtain the Gibbs energy. The Gibbs energy of a non-deformable slit pore can be transformed into the Helmholtz energy of a deformable slit pore using the Legendre-Fenchel transform. The Legendre-Fenchel transform corresponds to the Maxwell construction of equal areas.

Soret separation and thermo-osmosis in porous media.

When a temperature difference is applied over a porous medium soaked with a fluid mixture, two effects may be observed, a component separation (the Ludwig-Soret effect, thermodiffusion) and a pressure difference due to thermo-osmosis. In this work, we have studied both effects using non-equilibrium thermodynamics and molecular dynamics. We have derived expressions for the two characteristic parameters, the Soret coefficient and the thermo-osmotic coefficient in terms of phenomenological transport coefficients, and we show how they are related. Numerical values for these coefficients were obtained for a two-component fluid in a solid matrix where both fluid and solid are Lennard-Jones/spline particles. We found that both effects depend strongly on the porosity of the medium and weakly on the interactions between the fluid components and the matrix. The Soret coefficient depends strongly on whether the fluid is sampled from inside the porous medium or from bulk phases outside, which must be considered in experimental measurements using packed columns. If we use a methane/decane mixture in bulk as an example, our results for the Soret coefficient give that a temperature difference of 10 K will separate the mixture to about 49.5/50.5 and give no pressure difference. In a reservoir with 30% porosity, the separation will be 49.8/50.2, whereas the pressure difference will be about 15 bar. Thermo-osmotic pressures with this order or magnitude have been observed in frost-heave experiments.

Theory and simulation of shock waves: Entropy production and energy conversion.

We have considered a shock wave as a surface of discontinuity and computed the entropy production using nonequilibrium thermodynamics for surfaces. The results from this method, which we call the "Gibbs excess method" (GEM), were compared with results from three alternative methods, all based on the entropy balance in the shock-front region, but with different assumptions about local equilibrium. Nonequilibrium molecular dynamics (NEMD) simulations were used to simulate a thermal blast in a one-component gas consisting of particles interacting with the Lennard-Jones/spline potential. This provided data for the theoretical analysis. Two cases were studied, a weak shock with Mach number M≈2 and a strong shock with M≈6 and with a Prandtl number of the gas Pr≈1.4 in both cases. The four theoretical methods gave consistent results for the time-dependent surface excess entropy production for both Mach numbers. The internal energy was found to deviate only slightly from equilibrium values in the shock front. The pressure profile was found to be consistent with the Navier-Stokes equations. The entropy production in the weak and strong shocks were approximately proportional to the square of the Mach number and decayed with time at approximately the same relative rate. In both cases, some 97% of the total entropy production in the gas occurred in the shock wave. The GEM showed that most of the shock's kinetic energy was converted reversibly into enthalpy and entropy, and a small amount was dissipated as produced entropy. The shock waves traveled at almost constant speed, and we found that the overpressure determined from NEMD simulations agreed well with the Rankine-Hugoniot conditions for steady-state shocks.

Peltier effects in lithium-ion battery modeling.

A high battery temperature has been shown to be critical for lithium-ion batteries in terms of performance, degradation, and safety. Therefore, a precise knowledge of heat sources and sinks in the battery is essential. We have developed a thermal model for lithium-ion batteries, a model that includes terms not included before, namely, Peltier and Dufour heat effects. The model is derived using non-equilibrium thermodynamics for heterogeneous systems, the only theory which is able to describe in a systematic manner the coupling of heat, mass, and charge transport. The idea of this theory is to deal with surfaces as two-dimensional layers. All electrochemical processes in these layers are defined using excess variables, implying, for instance, that the surface has its own temperature. We show how the Peltier and Dufour heats affect a single cell and may produce an internal temperature rise of 8.5 K in a battery stack with 80 modules. The heat fluxes leaving the cell are also functions of these reversible heat effects. Most of the energy that is dissipated as heat occurs in the electrode surfaces and the electrolyte-filled separator. The analysis shows that better knowledge of experimental data on surface resistances, transport coefficients, and Dufour and Peltier heats is essential for further progress in thermal modeling of this important class of systems.

Adsorption of an Ideal Gas on a Small Spherical Adsorbent.

The ideal gas model is an important and useful model in classical thermodynamics. This remains so for small systems. Molecules in a gas can be adsorbed on the surface of a sphere. Both the free gas molecules and the adsorbed molecules may be modeled as ideal for low densities. The adsorption energy, Us, plays an important role in the analysis. For small adsorbents this energy depends on the curvature of the adsorbent. We model the adsorbent as a sphere with surface area Ω=4πR2, where R is the radius of the sphere. We calculate the partition function for a grand canonical ensemble of two-dimensional adsorbed phases. When connected with the nanothermodynamic framework this gives us the relevant thermodynamic variables for the adsorbed phase controlled by the temperature T, surface area Ω, and chemical potential µ. The dependence of intensive variables on size may then be systematically investigated starting from the simplest model, namely the ideal adsorbed phase. This dependence is a characteristic feature of small systems which is naturally expressed by the subdivision potential of nanothermodynamics. For surface problems, the nanothermodynamic approach is different, but equivalent to Gibbs' surface thermodynamics. It is however a general approach to the thermodynamics of small systems, and may therefore be applied to systems that do not have well defined surfaces. It is therefore desirable and useful to improve our basic understanding of nanothermodynamics.

Nanothermodynamic Description and Molecular Simulation of a Single-Phase Fluid in a Slit Pore.

We have described for the first time the thermodynamic state of a highly confined single-phase and single-component fluid in a slit pore using Hill's thermodynamics of small systems. Hill's theory has been named nanothermodynamics. We started by constructing an ensemble of slit pores for controlled temperature, volume, surface area, and chemical potential. We have presented the integral and differential properties according to Hill, and used them to define the disjoining pressure on the new basis. We identified all thermodynamic pressures by their mechanical counterparts in a consistent manner, and have given evidence that the identification holds true using molecular simulations. We computed the entropy and energy densities, and found in agreement with the literature, that the structures at the wall are of an energetic, not entropic nature. We have shown that the subdivision potential is unequal to zero for small wall surface areas. We have showed how Hill's method can be used to find new Maxwell relations of a confined fluid, in addition to a scaling relation, which applies when the walls are far enough apart. By this expansion of nanothermodynamics, we have set the stage for further developments of the thermodynamics of confined fluids, a field that is central in nanotechnology.

Small size effects in open and closed systems: What can we learn from ideal gases about systems with interacting particles?

Small systems have higher surface area-to-volume ratios than macroscopic systems. The thermodynamics of small systems therefore deviates from the description of classical thermodynamics. One consequence of this is that properties of small systems can be dependent on the system's ensemble. By comparing the properties in grand canonical (open) and canonical (closed) systems, we investigate how a small number of particles can induce an ensemble dependence. Emphasis is placed on the insight that can be gained by investigating ideal gases. The ensemble equivalence of small ideal gas systems is investigated by deriving the properties analytically, while the ensemble equivalence of small systems with particles interacting via the Lennard-Jones or the Weeks-Chandler-Andersen potential is investigated through Monte Carlo simulations. For all the investigated small systems, we find clear differences between the properties in open and closed systems. For systems with interacting particles, the difference between the pressure contribution to the internal energy, and the difference between the chemical potential contribution to the internal energy, are both increasing with the number density. The difference in chemical potential is, with the exception of the density dependence, qualitatively described by the analytic formula derived for an ideal gas system. The difference in pressure, however, is not captured by the ideal gas model. For the difference between the properties in the open and closed systems, the response of increasing the particles' excluded volume is similar to the response of increasing the repulsive forces on the system walls. This indicates that the magnitude of the difference between the properties in open and closed systems is related to the restricted movement of the particles in the system. The work presented in this paper gives insight into the mechanisms behind ensemble in-equivalence in small systems, and illustrates how a simple statistical mechanical model, such as the ideal gas, can be a useful tool in these investigations.

Fluctuation-Dissipation Theorems for Multiphase Flow in Porous Media.

A thermodynamic description of porous media must handle the size- and shape-dependence of media properties, in particular on the nano-scale. Such dependencies are typically due to the presence of immiscible phases, contact areas and contact lines. We propose a way to obtain average densities suitable for integration on the course-grained scale, by applying Hill's thermodynamics of small systems to the subsystems of the medium. We argue that the average densities of the porous medium, when defined in a proper way, obey the Gibbs equation. All contributions are additive or weakly coupled. From the Gibbs equation and the balance equations, we then derive the entropy production in the standard way, for transport of multi-phase fluids in a non-deformable, porous medium exposed to differences in boundary pressures, temperatures, and chemical potentials. Linear relations between thermodynamic fluxes and forces follow for the control volume. Fluctuation-dissipation theorems are formulated for the first time, for the fluctuating contributions to fluxes in the porous medium. These give an added possibility for determination of the Onsager conductivity matrix for transport through porous media. Practical possibilities are discussed.

A Legendre-Fenchel Transform for Molecular Stretching Energies.

Single-molecular polymers can be used to analyze to what extent thermodynamics applies when the size of the system is drastically reduced. We have recently verified using molecular-dynamics simulations that isometric and isotensional stretching of a small polymer result in Helmholtz and Gibbs stretching energies, which are not related to a Legendre transform, as they are for sufficiently long polymers. This disparity has also been observed experimentally. Using molecular dynamics simulations of polyethylene-oxide, we document for the first time that the Helmholtz and Gibbs stretching energies can be related by a Legendre-Fenchel transform. This opens up a possibility to apply this transform to other systems which are small in Hill's sense.

Entropy Production beyond the Thermodynamic Limit from Single-Molecule Stretching Simulations.

Single-molecular systems are a test bed to analyze to what extent thermodynamics applies when the size of the system is drastically reduced. Isometric and isotensional single-molecule stretching experiments and their theoretical interpretations have shown the lack of a thermodynamic limit at those scales and the nonequivalence between their corresponding statistical ensembles. This disparity between thermodynamic results obtained in both experimental protocols can also be observed in entropy production, as previous theoretical results have shown. In this work, we present results from molecular dynamics simulations of stretching of a typical polymer, polyethylene-oxide, where this framework is applied to obtain friction coefficients associated with stretching at the two different statistical ensembles for two different system sizes, from which the entropy production follows. In the smallest system, they are different up to a factor of 2, and for the bigger system, the difference is smaller, as predicted. In this way, we provide numerical evidence that a thermodynamic description is still meaningful for the case of single-molecule stretching.

When Thermodynamic Properties of Adsorbed Films Depend on Size: Fundamental Theory and Case Study.

Small system properties are known to depend on geometric variables in ways that are insignificant for macroscopic systems. Small system considerations are therefore usually added to the conventional description as needed. This paper presents a thermodynamic analysis of adsorbed films of any size in a systematic and general way within the framework of Hill's nanothermodynamics. Hill showed how to deal with size and shape as variables in a systematic manner. By doing this, the common thermodynamic equations for adsorption are changed. We derived the governing thermodynamic relations characteristic of adsorption in small systems, and point out the important distinctions between these and the corresponding conventional relations for macroscopic systems. We present operational versions of the relations specialized for adsorption of gas on colloid particles, and we applied them to analyze molecular simulation data. As an illustration of their use, we report results for CO2 adsorbed on graphite spheres. We focus on the spreading pressure, and the entropy and enthalpy of adsorption, and show how the intensive properties are affected by the size of the surface, a feature specific to small systems. The subdivision potential of the film is presented for the first time, as a measure of the film's smallness. For the system chosen, it contributes with a substantial part to the film enthalpy. This work can be considered an extension and application of the nanothermodynamic theory developed by Hill. It provides a foundation for future thermodynamic analyses of size- and shape-dependent adsorbed film systems, alternative to that presented by Gibbs.

Two-Phase Equilibrium Conditions in Nanopores.

It is known that thermodynamic properties of a system change upon confinement. To know how, is important for modelling of porous media. We propose to use Hill's systematic thermodynamic analysis of confined systems to describe two-phase equilibrium in a nanopore. The integral pressure, as defined by the compression energy of a small volume, is then central. We show that the integral pressure is constant along a slit pore with a liquid and vapor in equilibrium, when Young and Young-Laplace's laws apply. The integral pressure of a bulk fluid in a slit pore at mechanical equilibrium can be understood as the average tangential pressure inside the pore. The pressure at mechanical equilibrium, now named differential pressure, is the average of the trace of the mechanical pressure tensor divided by three as before. Using molecular dynamics simulations, we computed the integral and differential pressures, p ^ and p, respectively, analysing the data with a growing-core methodology. The value of the bulk pressure was confirmed by Gibbs ensemble Monte Carlo simulations. The pressure difference times the volume, V, is the subdivision potential of Hill, ( p - p ^ ) V = ϵ . The combined simulation results confirm that the integral pressure is constant along the pore, and that ϵ / V scales with the inverse pore width. This scaling law will be useful for prediction of thermodynamic properties of confined systems in more complicated geometries.

Gibbs Ensemble Monte Carlo Simulation of Fluids in Confinement: Relation between the Differential and Integral Pressures.

The accurate description of the behavior of fluids in nanoporous materials is of great importance for numerous industrial applications. Recently, a new approach was reported to calculate the pressure of nanoconfined fluids. In this approach, two different pressures are defined to take into account the smallness of the system: the so-called differential and the integral pressures. Here, the effect of several factors contributing to the confinement of fluids in nanopores are investigated using the definitions of the differential and integral pressures. Monte Carlo (MC) simulations are performed in a variation of the Gibbs ensemble to study the effect of the pore geometry, fluid-wall interactions, and differential pressure of the bulk fluid phase. It is shown that the differential and integral pressure are different for small pores and become equal as the pore size increases. The ratio of the driving forces for mass transport in the bulk and in the confined fluid is also studied. It is found that, for small pore sizes (i.e., < 5 σ fluid ), the ratio of the two driving forces considerably deviates from 1.

The thermal boundary resistance at semiconductor interfaces: a critical appraisal of the Onsager vs. Kapitza formalisms.

We critically readdress the definition of thermal boundary resistance at an interface between two semiconductors. By means of atomistic simulations we provide evidence that the widely used Kapitza formalism predicts thermal boundary resistance values in good agreement with the more rigorous Onsager non-equilibrium thermodynamics picture. The latter is, however, better suited to provide physical insight on interface thermal rectification phenomena. We identify the factors that determine the temperature profile across the interface and the source of interface thermal rectification. To this end we perform non-equilibrium molecular dynamics computational experiments on a Si-Ge system with a graded compositional interface of varying thickness, considering thermal bias of different sign.

Fluid-Fluid Interfaces of Multi-Component Mixtures in Local Equilibrium.

We derive in a new way that the intensive properties of a fluid-fluid Gibbs interface are independent of the location of the dividing surface. When the system is out of global equilibrium, this finding is not trivial: In a one-component fluid, it can be used to obtain the interface temperature from the surface tension. In other words, the surface equation of state can serve as a thermometer for the liquid-vapor interface in a one-component fluid. In a multi-component fluid, one needs the surface tension and the relative adsorptions to obtain the interface temperature and chemical potentials. A consistent set of thermodynamic properties of multi-component surfaces are presented. They can be used to construct fluid-fluid boundary conditions during transport. These boundary conditions have a bearing on all thermodynamic modeling on transport related to phase transitions.

Ensemble distribution for immiscible two-phase flow in porous media.

We construct an ensemble distribution to describe steady immiscible two-phase flow of two incompressible fluids in a porous medium. The system is found to be ergodic. The distribution is used to compute macroscopic flow parameters. In particular, we find an expression for the overall mobility of the system from the ensemble distribution. The entropy production at the scale of the porous medium is shown to give the expected product of the average flow and its driving force, obtained from a black-box description. We test numerically some of the central theoretical results.

Size and shape effects on the thermodynamic properties of nanoscale volumes of water.

Small systems are known to deviate from the classical thermodynamic description, among other things due to their large surface area to volume ratio compared to corresponding big systems. As a consequence, extensive thermodynamic properties are no longer proportional to the volume, but are instead higher order functions of size and shape. We investigate such functions for second moments of probability distributions of fluctuating properties in the grand-canonical ensemble, focusing specifically on the volume and surface terms of Hadwiger's theorem, explained in Klain, Mathematika, 1995, 42, 329-339. We resolve the shape dependence of the surface term and show, using Hill's nanothermodynamics [Hill, J. Chem. Phys., 1962, 36, 3182], that the surface satisfies the thermodynamics of a flat surface as described by Gibbs [Gibbs, The Scientific Papers of J. Willard Gibbs, Volume 1, Thermodynamics, Ox Bow Press, Woodbridge, Connecticut, 1993]. The Small System Method (SSM), first derived by Schnell et al. [Schnell et al., J. Phys. Chem. B, 2011, 115, 10911], is extended and used to analyze simulation data on small systems of water. We simulate water as an example to illustrate the method, using TIP4P/2005 and other models, and compute the isothermal compressibility and thermodynamic factor. We are able to retrieve the experimental value of the bulk phase compressibility within 2%, and show that the compressibility of nanosized volumes increases by up to a factor of two as the number of molecules in the volume decreases. The value for a tetrahedron, cube, sphere, polygon, etc. can be predicted from the same scaling law, as long as second order effects (nook and corner effects) are negligible. Lastly, we propose a general formula for finite reservoir correction to fluctuations in subvolumes.

Heat transport through a solid-solid junction: the interface as an autonomous thermodynamic system.

We perform computational experiments using nonequilibrium molecular dynamics simulations, showing that the interface between two solid materials can be described as an autonomous thermodynamic system. We verify the local equilibrium and give support to the Gibbs description of the interface also away from the global equilibrium. In doing so, we reconcile the common formulation of the thermal boundary resistance as the ratio between the temperature discontinuity at the interface and the heat flux with a more rigorous derivation from nonequilibrium thermodynamics. We also show that thermal boundary resistance of a junction between two pure solid materials can be regarded as an interface property, depending solely on the interface temperature, as implicitly assumed in some widely used continuum models, such as the acoustic mismatch model. Thermal rectification can be understood on the basis of different interface temperatures for the two flow directions.