Reconciling results of MOCVD of a CNT composite with equilibrium thermodynamics

Composition and microstructure of the composite films can be tailored by controlling the CVD process parameters if an appropriate model can be suggested for quantitative prediction of growth. This is possible by applying equilibrium thermodynamics. A modification of such standard modeling procedure was required to account for the deposition of a hybrid film comprised of carbon nanotubes (CNTs), metallic iron (Fe), and magnetite (Fe3O4), a composite useful for energy storage. In contrast with such composite nature of the deposits obtained by inert-ambient CVD using Fe(acac)3 as precursor, equilibrium thermodynamic modeling with standard procedure predicts the deposition of only Fe3C and carbon, without any co-deposition of Fe and Fe3O4. A modification of the procedure comprising chemical reasoning is therefore proposed herein, which predicts simultaneous deposition of FeO1-x, Fe3C, Fe3O4 and C. At high temperatures and in a carbon-rich atmosphere, these convert to Fe3O4, Fe and C, in agreement with experimental CVD. Close quantitative agreement between the modified thermodynamic modeling and experiment validates the reliability of the modified procedure. Understanding of the chemical process through thermodynamic modeling provides potential for control of CVD process parameters to achieve desired hybrid growth. (C) 2016 Elsevier B.V. All rights reserved.

Dynamical horizon entropy and equilibrium thermodynamics of generalized gravity theories

We study the relation between the thermodynamics and field equations of generalized gravity theories on the dynamical trapping horizon with sphere symmetry. We assume the entropy of a dynamical horizon as the Noether charge associated with the Kodama vector and point out that it satisfies the second law when a Gibbs equation holds. We generalize two kinds of Gibbs equations to Gauss-Bonnet gravity on any trapping horizon. Based on the quasilocal gravitational energy found recently for f(R) gravity and scalar-tensor gravity in some special cases, we also build up the Gibbs equations, where the nonequilibrium entropy production, which is usually invoked to balance the energy conservation, is just absorbed into the modified Wald entropy in the Friedmann-Robertson-Walker spacetime with slowly varying horizon. Moreover, the equilibrium thermodynamic identity remains valid for f(R) gravity in a static spacetime. Our work provides an alternative treatment to reinterpret the nonequilibrium correction and supports the idea that the horizon thermodynamics is universal for generalized gravity theories.

Out-of-equilibrium thermodynamics of quantum optomechanical systems

We address the out-of-equilibrium thermodynamics of an isolated quantum system consisting of a cavity optomechanical device. We explore the dynamical response of the system when driven out of equilibrium by a sudden quench of the coupling parameter and compute analytically the full distribution of the work generated by the process. We consider linear and quadratic optomechanical coupling, where the cavity field is parametrically coupled to either the position or the square of the position of a mechanical oscillator, respectively. In the former case we find that the average work generated by the quench is zero, whilst the latter leads to a non-zero average value. Through fluctuations theorems we access the most relevant thermodynamical figures of merit, such as the free energy difference and the amount of irreversible work generated. We thus provide a full charac- terization of the out-of-equilibrium thermodynamics in the quantum regime for nonlinearly coupled bosonic modes. Our study is the first due step towards the construction and full quantum analysis of an optomechanical machine working fully out of equilibrium.

A characterization of conserved quantities in non-equilibrium thermodynamics

The well-known Noether theorem in Lagrangian and Hamiltonian mechanics associates symmetries in the evolution equations of a mechanical system with conserved quantities. In this work, we extend this classical idea to problems of non-equilibrium thermodynamics formulated within the GENERIC (General Equations for Non-Equilibrium Reversible-Irreversible Coupling) framework. The geometric meaning of symmetry is reviewed in this formal setting and then utilized to identify possible conserved quantities and the conditions that guarantee their strict conservation. Examples are provided that demonstrate the validity of the proposed definition in the context of finite and infinite dimensional thermoelastic problems.

Non-equilibrium thermodynamics of harmonically trapped bosons

We apply the framework of non-equilibrium quantum thermodynamics to the physics of quenched small-size bosonic quantum gases in a harmonic trap. By studying the temporal behaviour of the Loschmidt echo and of the atomic density profile within the trap, which are informative of the non-equilibrium physics and the correlations among the particles, we establish a link with the statistics of (irreversible) work done on the system. This highlights interesting connections between the degree of inter-particle entanglement and the non-equilibrium thermodynamics of the system.

Cosmological equations and thermodynamics on apparent horizon in thick braneworld

We derive the generalized Friedmann equation governing the cosmological evolution inside the thick brane model in the presence of two curvature correction terms: a four-dimensional scalar curvature from induced gravity on the brane, and a five-dimensional Gauss-Bonnet curvature term. We find two effective four-dimensional reductions of the generalized Friedmann equation in some limits and demonstrate that the reductions but not the generalized Friedmann equation can be rewritten as the first law of equilibrium thermodynamics on the apparent horizon of thick braneworld.

Deformation with coupled chemical diffusion

The deformation of rocks is commonly intimately associated with metamorphic reactions. This paper is a step towards understanding the behaviour of fully coupled, deforming, chemically reacting systems by considering a simple example of the problem comprising a single layer system with elastic-power law viscous constitutive behaviour where the deformation is controlled by the diffusion of a single chemical component that is produced during a metamorphic reaction. Analysis of the problem using the principles of non-equilibrium thermodynamics allows the energy dissipated by the chemical reaction-diffusion processes to be coupled with the energy dissipated during deformation of the layers. This leads to strain-rate softening behaviour and the resultant development of localised deformation which in turn nucleates buckles in the layer. All such diffusion processes, in leading to Herring-Nabarro, Coble or “pressure solution” behaviour, are capable of producing mechanical weakening through the development of a “chemical viscosity”, with the potential for instability in the deformation. For geologically realistic strain rates these chemical feed-back instabilities occur at the centimetre to micron scales, and so produce structures at these scales, as opposed to thermal feed-back instabilities that become important at the 100–1000 m scales.

Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

Entropic bounds for multi-scale and multi-physics coupling in earth sciences

The ability to understand and predict how thermal, hydrological,mechanical and chemical (THMC) processes interact is fundamental to many research initiatives and industrial applications. We present (1) a new Thermal– Hydrological–Mechanical–Chemical (THMC) coupling formulation, based on non-equilibrium thermodynamics; (2) show how THMC feedback is incorporated in the thermodynamic approach; (3) suggest a unifying thermodynamic framework for multi-scaling; and (4) formulate a new rationale for assessing upper and lower bounds of dissipation for THMC processes. The technique is based on deducing time and length scales suitable for separating processes using a macroscopic finite time thermodynamic approach. We show that if the time and length scales are suitably chosen, the calculation of entropic bounds can be used to describe three different types of material and process uncertainties: geometric uncertainties,stemming from the microstructure; process uncertainty, stemming from the correct derivation of the constitutive behavior; and uncertainties in time evolution, stemming from the path dependence of the time integration of the irreversible entropy production. Although the approach is specifically formulated here for THMC coupling we suggest that it has a much broader applicability. In a general sense it consists of finding the entropic bounds of the dissipation defined by the product of thermodynamic force times thermodynamic flux which in material sciences corresponds to generalized stress and generalized strain rates, respectively.

A strain-based failure criterion for pillar stability analysis

Strain-based failure criteria have several advantages over stress-based failure criteria: they can account for elastic and inelastic strains, they utilise direct, observables effects instead of inferred effects (strain gauges vs. stress estimates), and model complete stress-strain curves including pre-peak, non-linear elasticity and post-peak strain weakening. In this study, a strain-based failure criterion derived from thermodynamic first principles utilising the concepts of continuum damage mechanics is presented. Furthermore, implementation of this failure criterion into a finite-element simulation is demonstrated and applied to the stability of underground mining coal pillars. In numerical studies, pillar strength is usually expressed in terms of critical stresses or stress-based failure criteria where scaling with pillar width and height is common. Previous publications have employed the finite-element method for pillar stability analysis using stress-based failure criterion such as Mohr-Coulomb and Hoek-Brown or stress-based scalar damage models. A novel constitutive material model, which takes into consideration anisotropy as well as elastic strain and damage as state variables has been developed and is presented in this paper. The damage threshold and its evolution are strain-controlled, and coupling of the state variables is achieved through the damage-induced degradation of the elasticity tensor. This material model is implemented into the finite-element software ABAQUS and can be applied to 3D problems. Initial results show that this new material model is capable of describing the non-linear behaviour of geomaterials commonly observed before peak strength is reached as well as post-peak strain softening. Furthermore, it is demonstrated that the model can account for directional dependency of failure behaviour (i.e. anisotropy) and has the potential to be expanded to environmental controls like temperature or moisture.

Theoretical and computational approaches on heterogeneous nucleation

Nucleation is the first step of a first order phase transition. A new phase is always sprung up in nucleation phenomena. The two main categories of nucleation are homogeneous nucleation, where the new phase is formed in a uniform substance, and heterogeneous nucleation, when nucleation occurs on a pre-existing surface. In this thesis the main attention is paid on heterogeneous nucleation. This thesis wields the nucleation phenomena from two theoretical perspectives: the classical nucleation theory and the statistical mechanical approach. The formulation of the classical nucleation theory relies on equilibrium thermodynamics and use of macroscopically determined quantities to describe the properties of small nuclei, sometimes consisting of just a few molecules. The statistical mechanical approach is based on interactions between single molecules, and does not bear the same assumptions as the classical theory. This work gathers up the present theoretical knowledge of heterogeneous nucleation and utilizes it in computational model studies. A new exact molecular approach on heterogeneous nucleation was introduced and tested by Monte Carlo simulations. The results obtained from the molecular simulations were interpreted by means of the concepts of the classical nucleation theory. Numerical calculations were carried out for a variety of substances nucleating on different substances. The classical theory of heterogeneous nucleation was employed in calculations of one-component nucleation of water on newsprint paper, Teflon and cellulose film, and binary nucleation of water-n-propanol and water-sulphuric acid mixtures on silver nanoparticles. The results were compared with experimental results. The molecular simulation studies involved homogeneous nucleation of argon and heterogeneous nucleation of argon on a planar platinum surface. It was found out that the use of a microscopical contact angle as a fitting parameter in calculations based on the classical theory of heterogeneous nucleation leads to a fair agreement between the theoretical predictions and experimental results. In the presented cases the microscopical angle was found to be always smaller than the contact angle obtained from macroscopical measurements. Furthermore, molecular Monte Carlo simulations revealed that the concept of the geometrical contact parameter in heterogeneous nucleation calculations can work surprisingly well even for very small clusters.

Thermocapillary Effect in a Thermal Rheological Jet

The process of die swell in polymer jets is an important feature within polymer processing and can be explained through a study of its rheological effects. The existence of a thermocapillary effect, driven by the gradient of its surface tension, should be considered when examining a thermal jet that has a non-uniform temperature distribution on its free surface, as in various polymer processings. Both the rheological effect and thermocapillary effect on die swell can be studied numerically through a finite element method as used on a two-dimensional and unsteady model, in which a Coleman-Noll second-order fluid model is employed. The results show that the expanding angle depends on both the rheological property of the fluid and the pressure at the vessel exit. Although both the thermocapillary and the rheological effects contribute to the cross-section expansion of the fluid jet, the latter is more important in determining the expansion.

Oscillatory instability of Rayleigh-Marangoni-Bénard convection in two-layer liquid systems

The oscillatory behaviour of the Rayleigh-Marangoni-Bénard convective instability (R-M-B instability) regarding two combinations of two-layer fluid systems has been investigated theoretically and numerically. For the two-layer system of Silicone oil (10cSt) over Fluorinert (FC70), both linear instability analysis and 2D numerical simulation show that the instability of the system depends strongly on the depth ratio Hr = H1/H2 of the two-layer liquid. The oscillatory regime at the onset of R-M-B convection enlarges with reducing Γ = Ra/Ma values. In the two-layer system of Silicone oil (2cSt) over water, it loses its stability and onsets to steady convection at first, then the steady convection bifurcates to oscillatory convection with increasing Rayleigh number Ra. This behaviour was found through numerical simulation above the onset of steady convection in the case of r = 2.9, ε=(Ra-Ruc)/Rac = 1.0, and Hr = 0.5. Our findings are different from the previous study of the Rayleigh-Benard instability and show the strong effects of the thermocapillary force at the interface on the time-dependent oscillations at or after the onset of convection. We propose a secondary oscillatory instability mechanism to explain the experimental observation of Degen et al. [Phys. Rev. E, 57 (1998), 6647-6659].

Part I. Influence of membrane morphology on ion-exchange and charge propagation in composite polyelectrolyte electrode coatings. Part II. Dynamic surface tension measurements of polarization relaxation at mercury pool electrodes by the method of Wilhelmy.

Part I. Novel composite polyelectrolyte materials were developed that exhibit desirable charge propagation and ion-retention properties. The morphology of electrode coatings cast from these materials was shown to be more important for its electrochemical behavior than its chemical composition.

Part II. The Wilhelmy plate technique for measuring dynamic surface tension was extended to electrified liquid-liquid interphases. The dynamical response of the aqueous NaF-mercury electrified interphase was examined by concomitant measurement of surface tension, current, and applied electrostatic potential. Observations of the surface tension response to linear sweep voltammetry and to step function perturbations in the applied electrostatic potential (e.g., chronotensiometry) provided strong evidence that relaxation processes proceed for time-periods that are at least an order of magnitude longer than the time periods necessary to establish diffusion equilibrium. The dynamical response of the surface tension is analyzed within the context of non-equilibrium thermodynamics and a kinetic model that requires three simultaneous first order processes.

Site-Specific Isotopes in Small Organic Molecules

Stable isotope geochemistry is a valuable toolkit for addressing a broad range of problems in the geosciences. Recent technical advances provide information that was previously unattainable or provide unprecedented precision and accuracy. Two such techniques are site-specific stable isotope mass spectrometry and clumped isotope thermometry. In this thesis, I use site-specific isotope and clumped isotope data to explore natural gas development and carbonate reaction kinetics. In the first chapter, I develop an equilibrium thermodynamics model to calculate equilibrium constants for isotope exchange reactions in small organic molecules. This equilibrium data provides a framework for interpreting the more complex data in the later chapters. In the second chapter, I demonstrate a method for measuring site-specific carbon isotopes in propane using high-resolution gas source mass spectrometry. This method relies on the characteristic fragments created during electron ionization, in which I measure the relative isotopic enrichment of separate parts of the molecule. My technique will be applied to a range of organic compounds in the future. For the third chapter, I use this technique to explore diffusion, mixing, and other natural processes in natural gas basins. As time progresses and the mixture matures, different components like kerogen and oil contribute to the propane in a natural gas sample. Each component imparts a distinct fingerprint on the site-specific isotope distribution within propane that I can observe to understand the source composition and maturation of the basin. Finally, in Chapter Four, I study the reaction kinetics of clumped isotopes in aragonite. Despite its frequent use as a clumped isotope thermometer, the aragonite blocking temperature is not known. Using laboratory heating experiments, I determine that the aragonite clumped isotope thermometer has a blocking temperature of 50-100°C. I compare this result to natural samples from the San Juan Islands that exhibit a maximum clumped isotope temperature that matches this blocking temperature. This thesis presents a framework for measuring site-specific carbon isotopes in organic molecules and new constraints on aragonite reaction kinetics. This study represents the foundation of a future generation of geochemical tools for the study of complex geologic systems.