A Monte Carlo Simulation of Locally Nonchaotic Barrier

Total Page:16

File Type:pdf, Size:1020Kb

A Monte Carlo Simulation of Locally Nonchaotic Barrier Production of Useful Work in an Isothermal Cycle: A Monte Carlo Simulation of a Locally Nonchaotic Energy Barrier Y. Qiao,1,2,* Z. Shang1 1 Program of Materials Science and Engineering, University of California – San Diego, La Jolla, CA 92093, U.S.A. 2 Department of Structural Engineering, University of California – San Diego, La Jolla, CA 92093-0085, U.S.A. * Corresponding author (Email: [email protected]; phone: 858-534-3388) Abstract: A Monte Carlo simulation is performed on a billiard-type model system, which contains a locally nonchaotic energy barrier. Without extensive particle collision across the energy barrier, the system steady state is nonequilibrium; that is, the particles follow a non-Boltzmann distribution. Remarkably, as the energy barrier varies in an isothermal cycle, the total produced work is greater than the total consumed work, because of the asymmetry in the cross-influence of the thermally correlated thermodynamic driving forces. Such a phenomenon cannot be explained by the second law of thermodynamics. Similar anomalous effects may be achieved if the barrier is switchable or asymmetric. In essence, the energy barrier is a spontaneously nonequilibrium dimension. It is fundamentally different from Maxwell’s demon, unrelated to the physical nature of information. KEYWORDS: Nonequilibrium; Nonchaotic; The second law of thermodynamics; Energy barrier; Monte Carlo simulation 1. Introduction Randomness is a fundamental concept in statistical mechanics [1,2]. An ergodic and chaotic system always reaches thermodynamic equilibrium, while a nonergodic or nonchaotic system may not [e.g., 3-5]. The latter should be analyzed in the framework of nonequilibrium stochastic thermodynamics [6]. 1 Figure 1. (A) In a vertical 푦 − 푧 plane in a gravitational field (푔), if 푧̂ ≪ , at the steady state the particle density does not follow the Maxwell-Boltzmann distribution. (B) The billiard-type model system, wherein elastic particles randomly move in the horizontal 푥 − 푦 dimension in between a upper “plateau” and a lower “plain”, across a narrow transition step. If 푧̂ ≪ , the steady-state ⁄ particle density ratio ( = ) is inherently nonequilibrium. One example is shown in Figure 1(A), in which billiard-like elastic particles randomly move in a vertical 푦 − 푧 plane in a gravitational field (푔). The 푦 dimension is infinitely large. The upper boundary (푧 = 푧̂) is a specular wall, with 푧̂ being the plane height. The lower boundary (푧 = 0) is a diffusive wall at constant temperature (푇), where the reflected particle velocity (푣) is governed by the two-dimensional (2D) Maxwell-Boltzmann distribution [7]: 푝(푣) = ∙ exp − , with 푝 being the probability density, 푚 the particle mass, 퐾 = 푘 푇 the average particle kinetic energy (퐾), and 푘 the Boltzmann constant. If the mean free path of the particles () is small compared to the plane height (푧̂), the system is ergodic and chaotic. At 푧 = 푧̂, the local particle density, , is equal to [8], where is the particle density ratio, is the local ̂ particle density at the bottom of the plane, =푒 is the Boltzmann factor, and 훽 = . If ≫ 푧̂, the characteristic interaction time 푡̅ = ⁄푣̅ is much longer than the characteristic event duration 푡̅ =2푧̂⁄푣̅ , where 푣̅ is the average particle velocity and 푣̅ is the average 푧-dimension particle velocity component (푣). On the time scale of 푡̅ , few particle collisions take place and thus, the particle trajectories tend to be nonchaotic. The probability for a particle to travel from 푧 =0 to 푧 = 푧̂ is mainly determined by 푣, relatively unrelated to the 푦- dimension momentum. As a first-order approximation, the particle density ratio at 푧̂ can be 2 assessed as = ∫ 푝 (푣 )d푣 = , where 푝 (푣 ) = exp − is the 2D Maxwell- ̂ Boltzmann distribution of 푣. Clearly, the nonchaotic is different from . Figure 2. (A) In a system containing a spontaneously nonequilibrium dimension (SND), the steady-state particle distribution is intrinsically in a non-Boltzmann form; i.e., ≠ . (B) In an ergodic and chaotic system, at thermodynamic equilibrium, = . Figure 2(A) depicts two large ergodic and chaotic areas (푖 and 푗) that are separated by a nonchaotic and/or nonergodic barrier. A large number of elastic particles randomly move in the system. The crossing ratio of the barrier, = ⁄, does not equal to the Boltzmann factor, 훿 = ∆ 푒 , where ∆퐸 is the potential difference between the two areas, and and are the probabilities for the particles to cross the barrier from area 푖 to 푗 and from area 푗 to 푖, respectively. Hereafter, such a barrier will be referred to as a spontaneously nonequilibrium dimension (SND). Without any specific knowledge of the system microstate, SND offers a mechanism to reach a ⁄ nonequilibrium steady state; i.e., ≠ 훿, where = , with and being the steady-state particle densities in areas 푖 and 푗, respectively. In comparison, Figure 2(B) shows the same two large areas connected through a regular open gate, across which 훿 = 훿; at thermodynamic equilibrium, = 훿 [8]. In the current study, we investigate a SND-based system that satisfies the following six conditions: 1) The system is divided by the SND into two large areas (푖 and 푗). 2) Areas 푖 and 푗 are dominated by two different thermodynamic driving forces (퐹 and 퐹), respectively. 3) 퐹 and 퐹 are thermally correlated through particle diffusion; the energy change associated with their cross- 3 influence (퐹 = and 퐹 = ) is from heat, where 푥 and 푥 are the conjugate variables of 퐹 and 퐹, respectively. The operations of 퐹 and 퐹 are 4) reversible and 5) independent of each other, and 6) do not rely on temperature variation. 2. Monte Carlo Simulation 2.1 The billiard-type model system Figure 1(B) shows a SND-based billiard-type model system. A large number of elastic particles randomly move in the horizontal 푥푦 dimension. A uniform gravitational field (푔) is along the out-of-plane direction, −푧. The central area is higher, which will be referred to as “plateau”. The surrounding lower area will be referred to as “plain”. The plateau and the plain are separated by a transition step, which imposes an energy barrier to the particle motion from the plain to the plateau. The plateau height (푧̂) can be changed by the lifting force on the plateau, 퐹. The total particle number 푁 = 푁 + 푁 is constant, with subscripts “P” and “G” indicating the plain and the plateau, respectively. The plain area (퐴) can be adjusted by moving the outer system boundary; the area of the plateau (퐴) is fixed. The thermodynamic driving forces under investigation are the lifting force (퐹) and the in-plane pressure of the plain (푃), with the conjugate variables being 푧̂ and −퐴, respectively. It is assumed that i) the particle motion is frictionless; ii) the changes of 푧̂ and 퐴 are reversible; iii) the transition step is smooth, i.e., as the particles move across it, no energy is dissipated; and iv) the environment is a perfect heat reservoir of a constant average particle kinetic energy (퐾 = 푘푇). Moreover, v) we choose to study a system wherein 푧̂ is much smaller than the plain/plateau size; if 푔 =0, 푧̂ has little influence on the particle distribution. The model system meets the six conditions listed in the last paragraph of the introductory section: The plateau and the plain are separated by the energy barrier; they are dominated by 퐹 and 푃, respectively; 퐹 is directional and has a different governing equation from 푃; a variation in 퐴 or 푧̂ would cause a particle redistribution across the transition step, resulting in an exchange of thermal energy with the environment; 퐴 and 푧̂ can be adjusted reversibly and independently, and 퐾 is maintained constant by the thermal bath. 4 When 푧̂ ≫ , the system is ergodic and chaotic. If 푧̂ ≪ , the transition step becomes a locally nonchaotic SND. Notice that the in-plane pressure (푃), the plain area (퐴), the lifting force (퐹), and the plateau height (푧̂) can be measured from the surface of the system. In the following discussion, the system state will be defined by these macroscopic variables; no specific microstate will be directly observed. 2.2 Governing equations With Condition (v) in Section 2.1, the model system may be viewed as two large areas (the plain and the plateau) bordering each other across the line of transition step. The in-plane pressure of the plain is governed by [7] 푃퐴 = 푁퐾 (1) The lifting force on the plateau (퐹) contains two components: the particle weight 퐹 = 푚푔푁, and the centrifugal force (퐹) caused by the particles changing direction in the transition step. Denote the characteristic time for a particle to change its direction in the transition step by ∆푡̃. At the steady state, during ∆푡̃, on average 퐿(푣∆푡̃) particles pass through the transition step, where 푣 represents the average 푧̃-dimension velocity component of these particles, and 퐿 is the plateau circumference. The average centrifugal force per particle is 푚푣⁄∆푡̃. Thus, 퐹 is a fraction of 퐹, where 퐾 = 푚푣 and 퐷 is the plateau size. When 퐷 is much larger than 푧 = 퐾⁄(푚푔), 퐹 is negligible compared to 퐹, since 퐾 is on the same scale as 퐾. Under this condition, the lifting force may be calculated as 퐹 = 푚푔푁 (2) 2.3 Setup for the Monte Carlo simulation Nonequilibrium stochastic processes are often analyzed through Monte Carlo (MC) simulation. The program code of our MC simulation is partly based on [9] and can be downloaded from [10]. The simulation setup is illustrated in Figure 3. A relatively large number of billiard-like elastic particles freely move in a square box in the 푥 − 푦 plane.
Recommended publications
  • Interconversion of Hydrated Protons at the Interface Between Liquid Water and Platinum
    Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is © the Owner Societies 2019 Supporting Information: Interconversion of Hydrated Protons at the Interface between Liquid Water and Platinum Peter S. Rice, Yu Mao, Chenxi Guo and P. Hu* School of Chemistry and Chemical Engineering, The Queen’s University of Belfast, Belfast BT9 5AG, N. Ireland Email: [email protected] Contents Page S1. Ab Initio Molecular Dynamics (AIMD) based Umbrella Sampling (US) ................. 14 S2. Unit Cell Conversion ................................................................................................ 17 S3. Atomic Density Profile ............................................................................................. 17 S4. Radial Pair Distribution Function, (g(r)) ................................................................. 18 S5. Coordination Number (CN) ..................................................................................... 18 S6. Charge Density Difference ....................................................................................... 19 S7. Self-Diffusion Coefficient (D0) .................................................................................. 21 S8. Dipole Orientation Analysis ..................................................................................... 22 13 S1. Ab Initio Molecular Dynamics (AIMD) based Umbrella Sampling (US) x To study the interchange of species at the Pt(111)/H2n+1On interface, MD simulations based on umbrella sampling were used to investigate
    [Show full text]
  • Thermodynamic Temperature
    Thermodynamic temperature Thermodynamic temperature is the absolute measure 1 Overview of temperature and is one of the principal parameters of thermodynamics. Temperature is a measure of the random submicroscopic Thermodynamic temperature is defined by the third law motions and vibrations of the particle constituents of of thermodynamics in which the theoretically lowest tem- matter. These motions comprise the internal energy of perature is the null or zero point. At this point, absolute a substance. More specifically, the thermodynamic tem- zero, the particle constituents of matter have minimal perature of any bulk quantity of matter is the measure motion and can become no colder.[1][2] In the quantum- of the average kinetic energy per classical (i.e., non- mechanical description, matter at absolute zero is in its quantum) degree of freedom of its constituent particles. ground state, which is its state of lowest energy. Thermo- “Translational motions” are almost always in the classical dynamic temperature is often also called absolute tem- regime. Translational motions are ordinary, whole-body perature, for two reasons: one, proposed by Kelvin, that movements in three-dimensional space in which particles it does not depend on the properties of a particular mate- move about and exchange energy in collisions. Figure 1 rial; two that it refers to an absolute zero according to the below shows translational motion in gases; Figure 4 be- properties of the ideal gas. low shows translational motion in solids. Thermodynamic temperature’s null point, absolute zero, is the temperature The International System of Units specifies a particular at which the particle constituents of matter are as close as scale for thermodynamic temperature.
    [Show full text]
  • Ideal Gas with a Varying (Negative Absolute) Temperature
    Ideal gas with a varying (negative absolute) temperature: An alternative to dark energy? Subhajit Saha1, Anindita Mondal2 and Christian Corda3 February 5, 2018 1 Department of Mathematics, Panihati Mahavidyalaya, Kolkata 700110, West Bengal, India, e-mail: [email protected]. 2Department of Astrophysics and Cosmology, S. N. Bose National Centre for Basic Sciences, JD Block, Sector III, Salt Lake City, Kolkata 700106, India, e-mail: [email protected]. 3 Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), P.O. Box 55134-441, Maragha, Iran, e-mail: [email protected]. Abstract The present work is an attempt to investigate whether the evolutionary history of the Universe from the offset of inflation can be described by assuming the cosmic fluid to be an ideal gas with a specific gas constant but a varying negative absolute temperature (NAT). The motivation of this work is to search for an alternative to the "exotic" and "supernatural" dark energy (DE). In fact, the NAT works as an “effective quintessence” and there is need to deal neither with exotic matter like DE nor with modified gravity theories. For the sake of completeness, we release some clarifications on NATs in Section 3 of the paper. Keywords: Negative absolute temperature; Cosmic accel- eration; Ideal gas law; Dark energy arXiv:1802.00742v1 [physics.gen-ph] 26 Jan 2018 PACS Numbers: 98.80.-k 1 Introduction Since the intriguing discovery of the late time cosmic acceleration in 1998, which was expected by the Indian physicist B. Sidharth [1] before the cosmological observations of the High-Z Supernova Search Team [2] and of the Supernova 1 Cosmology Project [3], there have been rigorous attempts to incorporate this unexpected observations into standard cosmology.
    [Show full text]
  • Selvi Afstudeerverslag
    Eindhoven University of Technology MASTER On the reservoir model for CO2-laser amplification Selvi, S.C. Award date: 2018 Link to publication Disclaimer This document contains a student thesis (bachelor's or master's), as authored by a student at Eindhoven University of Technology. Student theses are made available in the TU/e repository upon obtaining the required degree. The grade received is not published on the document as presented in the repository. The required complexity or quality of research of student theses may vary by program, and the required minimum study period may vary in duration. General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain On the reservoir model for CO2-laser amplification S.C. Selvi November 15, 2018 2 Abstract In many technological applications, such as extreme ultraviolet lithography, high-intensity CO2- laser pulses are used, which are created in a CO2-laser amplification system. Research concerned with CO2-laser amplification tries to describe the dynamics of such amplifications systems to be able to make predictions about its characteristics and performance. The system under consider- ation is a single plasma cell of the amplification system having a pre-mixture of CO2-CO-N2-He and infra-red radiation with a wavelength of 10:58 µm and 10:206 µm.
    [Show full text]
  • A Spin-Half System As a Working Substance of a Heat Engine: Exploring Its finite-Time Thermodynamic Quantities
    A spin-half system as a working substance of a heat engine: exploring its finite-time thermodynamic quantities A Dissertation Submitted to the Department of Physics Addis Ababa University in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy (in Physics) By Tolasa Adugna Dima Advisor Dr. Mulugeta Bekele (Assoc. Professor), Addis Ababa University June 27, 2017 Declaration I hereby declare that except where specific reference is made to the work of others,the contents of this dissertation are original and have not been submitted in whole or in part for consideration for any other degree or qualification in this, or any other University. This dissertation is the result of my own work and that all sources of material used for the dissertation have been duly acknowledged. Name: Tolasa Adugna Dima Signature: Date: This PhD dissertation has been submitted for examination with my approval as uni- versity advisor. Advisor’s Name: Dr. Mulugeta Bekele Signature: Date: Approved by the Examination Committee Signature Date External Examiner: Prof. Purushotam D. Gujrati Internal Examiner: Dr. Lemi Demeyu Thesis Advisor: Dr. Mulugeta Bekele Chairman: Dr. Deribe Hirpo i This work is dedicated solely to Hawi B. Mosissa. ii Acknowledgements First and above all, I praise Waaqayyo, the almighty for providing me the opportunity to do this work. Next, my special thank goes to my advisor Dr. Mulugeta Bekele for his unreserved support throughout my study. From his friendly approach and fa- therhood advice, I have learned and shared lots of experiences. My sincere thanks also goes to Dr. Lemi Demeyu, course instructor and ex-head of physics Department, AAU, for his encouragement, valuable comments and support during my PhD study.
    [Show full text]
  • David Elsing
    Department of Physics and Astronomy University of Heidelberg Bachelor Thesis in Physics submitted by David Elsing born in Freiburg im Breisgau (Germany) 2018 Plasma screening effects in laser-generated plasmas This Bachelor Thesis has been carried out by David Elsing at the Max Planck Institute for Nuclear Physics in Heidelberg under the supervision of PD Dr. Adriana P´alffy-Bußand Dr. Yuanbin Wu Abstract Plasma screening effects affect the rate at which nuclear reactions occur in hot and dense astrophysical plasmas. Their direct measurement in the laboratory has not succeeded so far and is vital for our understanding of nuclear processes occurring in stellar nucle- osynthesis. A promising route for such measurements involve petawatt laser facilities capable of producing hot and dense plasmas in the laboratory. In this work, we consider and compare theoretically four plasma screening models originally introduced for astro- physical plasmas, and apply them to nuclear reactions in laser-generated plasmas for a large range of densities and temperatures. We consider two astrophysically relevant nuclear reactions, namely 13C(α; n) 16O, which is one of the important helium burning processes in advanced stellar phases as well as one of the interesting neutron sources for the slow neutron caption process of nucleosynthesis, and 7Li(d; α)α, which might be the key to explain the lithium production problem of the standard big bang nucle- osynthesis. Our results anticipate experiments on nuclear reactions in plasmas being rendered possible at high power laser facilities such as the Extreme Light Infrastructure under construction today, and shedding light on the validity of plasma screening models.
    [Show full text]
  • The Boltzmann Distribution
    THE BOLTZMANN DISTRIBUTION ZHENGQU WAN Abstract. This paper introduces some of the basic concepts in statistical mechanics. It focuses how energy is distributed to different states of a physical system, i.e. under certain hypothesis, it obeys the Boltzmann distribution. I will demonstrate three ways that the Boltzmann distribution will arise. Contents 1. The Motivating Problem 1 2. Boltzmann Distribution Arises from the Principle of Indifference 4 3. Validity of the Principle of Indifference 7 4. Boltzmann Distribution Arises as the Maximal Entropy Distribution 9 5. Application to Particles in the Box 12 6. Application to Gas Particles under Gravity 13 7. Boltzmann Distribution Arises as the Equilibrium Distribution 13 Acknowledgments 14 References 14 1. The Motivating Problem Consider the following game. Suppose there are N people and everyone has some cash. For simplicity, each person's cash is a natural number. Now imagine the following process. At each second we randomly uniformly pick a person A. If A has positive amount of cash, then we randomly uniformly pick a person B(it could be the same person) and transfer a dollar from A to B. Keep this game running for a sufficiently long period of time, will the distribution of wealth among this group of people reach an equilibrium. If so, what will the equilibrium distribution be? The answer to the above question is affirmative. Let Xi denote the wealth of the i-th person. Let u denote the average wealth so that Nu is a positive integer. The game can be modeled by a Markov chain where the state space is n X S = (X1; :::XN ): Xi 2 Z>0; Xi = nu i=1 Theorem 1.1.
    [Show full text]
  • Partition Function Estimation in Computational Protein Design with Continuous-Label Markov Random fields
    Partition function estimation in computational protein design with continuous-label Markov random fields Aditya Mukund A thesis submitted to the Department of Computer Science for honors Duke University Durham, North Carolina 2017 Contents 1 Abstract 2 2 Computational Structure-Based Protein Design 3 2.1 Conformation Optimization and the GMEC Model . .3 2.2 Computing Ensemble Properties of Proteins . .5 3 From Protein Design to Graphical Models 8 3.1 Constructing the MRF . .8 3.2 A Toy Example . .9 4 Mean-Field Approximations 12 4.1 Variational Optimization . 12 4.2 The Mean Field Algorithm . 14 5 Tree-Reweighted Belief Propagation 17 5.1 The Marginal Polytope . 17 5.2 Jensen's Inequality . 19 5.3 A Message-Passing Algorithm . 20 6 Reproducing Kernel Hilbert Spaces 25 6.1 Constructing the RKHS . 25 6.2 RKHS Representations of Functions . 27 6.3 The Mean Map . 28 7 Results and Discussion 30 7.1 Bounds on log Z .................................. 30 7.2 HIV Envelope Glycoprotein gp120 . 37 7.3 Discussion . 41 8 Acknowledgements 43 1 1 Abstract Proteins perform a variety of biological tasks, and drive many of the dynamic processes that make life possible. Computational structure-based protein design (CSPD) involves com- puting optimal sequences of amino acids with respect to particular backbones, or folds, in order to produce proteins with novel functions. In particular, it is crucial to be able to accurately model protein-protein interfaces (PPIs) in order to realize desired functionalities. Accurate modeling of PPIs raises two significant considerations. First, incorporating contin- uous side-chain flexibility in the design process has been shown to significantly improve the quality of designs.
    [Show full text]
  • Introduction to the Transition State Theory 29
    DOI: 10.5772/intechopen.78705 ProvisionalChapter chapter 2 Introduction toto thethe TransitionTransition StateState TheoryTheory Petr PtáPtáček,ček, FrantiFrantišekšek ŠŠoukaloukal andand TomáTomášš OpravilOpravil Additional information isis availableavailable atat thethe endend ofof thethe chapterchapter http://dx.doi.org/10.5772/intechopen.78705 Abstract The transition state theory (TST), which is also known as theory of absolute reaction rates (ART) and the theory of activated state (complex), is essentially a refined version of crude collision theory, which treats the reacting molecules as the rigid spheres without any internal degree of freedom. The theory explains the rate of chemical reaction assuming a special type of chemical equilibrium (quasi-equilibrium) between the reactants and acti- vated state (transition state complex). This special molecule decomposes to form the products of reaction. The rate of this reaction is then equal to the rate of decomposition of activated complex. This chapter also explains the limitation of TST theory and deals with the kinetics isotope effect. Keywords: transition state theory, theory of absolute reaction rates, activated complex, equilibrium constant, reaction rate, kinetics, activation energy, frequency factor, partition function, thermodynamics of activated state, kinetics isotope effect 1. Introduction “Two major findings have, however, stood the test of time. First, there was the unmistakable evidence that molecular collisions play the all-important role in communicating the activation energy to the molecules which are to be transformed and, second, there was the recognition of the activation energy (derivable for the Arrhenius Law of the temperature dependence) as a dominant, though by no means the sole, factor determining reactivity in general.”—C.N.
    [Show full text]
  • Jarzynski Equality and Crooks Relation for Local Models of Air–Sea Interaction
    Earth Syst. Dynam., 12, 689–708, 2021 https://doi.org/10.5194/esd-12-689-2021 © Author(s) 2021. This work is distributed under the Creative Commons Attribution 4.0 License. Jarzynski equality and Crooks relation for local models of air–sea interaction Achim Wirth1 and Florian Lemarié2 1LEGI, Univ. Grenoble Alpes, CNRS, Grenoble INP, 38000 Grenoble, France 2LJK, Univ. Grenoble Alpes, Inria, CNRS, Grenoble INP, 38000 Grenoble, France Correspondence: Achim Wirth ([email protected]) Received: 26 October 2020 – Discussion started: 15 December 2020 Revised: 30 March 2021 – Accepted: 23 April 2021 – Published: 4 June 2021 Abstract. We show that the most prominent of the work theorems, the Jarzynski equality and the Crooks rela- tion, can be applied to the momentum transfer at the air–sea interface using a hierarchy of local models. In the more idealized models, with and without a Coriolis force, the variability is provided from Gaussian white noise which modifies the shear between the atmosphere and the ocean. The dynamics is Gaussian, and the Jarzynski equality and Crooks relation can be obtained analytically solving stochastic differential equations. The more involved model consists of interacting atmospheric and oceanic boundary layers, where only the dependence on the vertical direction is resolved, the turbulence is modeled through standard turbulent models and the stochas- ticity comes from a randomized drag coefficient. It is integrated numerically and can give rise to a non-Gaussian dynamics. Also in this case the Jarzynski equality allows for calculating a dynamic beta βD of the turbulent −1 fluctuations (the equivalent of the thermodynamic beta β D (kBT ) in thermal fluctuations).
    [Show full text]
  • Isobaric Heat Capacity of Supercritical Fluids
    ISOBARIC HEAT CAPACITY OF SUPERCRITICAL FLUIDS: EXPERIMENTAL MEASUREMENTS AND MODELING A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy by Mitchell Price Ishmael May 2017 © 2017 Mitchell Price Ishmael ALL RIGHTS RESERVED ISOBARIC HEAT CAPACITY OF SUPERCRITICAL FLUIDS: EXPERIMENTAL MEASUREMENTS AND MODELING Mitchell Price Ishmael, Ph. D. Cornell University 2017 Fluids at high temperature and pressure are commonplace across engineering applications. These fluids are sometimes the product of a harsh environment, as in subsurface operations such as oil and gas development or geothermal energy production, or at other times, exploited for their desirable transport and thermodynamic properties in, for example, chemical extractions and material synthesis. At these elevated conditions, many single and multi-component fluids exist in a supercritical state, where small changes in temperature and pressure cause significant property variations without induction of a phase change. A fluid’s isobaric heat capacity (Cp), a property important to most thermal processes, fluctuates strongly in the supercritical region, ranging from a few times the ideal gas state Cp to infinity exactly at the critical point for pure fluids. Accurately capturing the thermophysical behavior of supercritical pure fluids and mixtures is necessary for optimal design and efficient operation of numerous engineering processes that utilize these substances. The main goal of this research was to increase understanding of isobaric heat capacity property changes in supercritical fluids. The first step in a threefold approach was to design and construct a calorimeter for precise and accurate measurements of Cp in the supercritical region.
    [Show full text]
  • Monte Carlo Method Elena Akhmatskaya (BCAM) January 24 2013 Introduction: Questions to Answer
    Extensions, Variations and Applications of Monte Carlo Method Elena Akhmatskaya (BCAM) January 24 2013 Introduction: Questions to Answer Why so many different Monte Carlo based algorithms are around? • MC is one of the most universal algorithms but not always the most efficient way to treat a problem • The effective, good MC algorithms are in high demand What are good MC algorithms? • Accurate ü Create a Markov chain ü Obey the detailed balance condition ü Ergodic • Fast ü Identify the most important samples ü Identification process is fast • The ideas often come from application The purpose of this course is to introduce some good MC algorithms and show the ways of constructing efficient Monte Carlo based algorithms. Introduction: Topics to Cover • Monte Carlo Applications to Molecular Systems • Hamiltonian Dynamics • Numerical Methods for Hamiltonian Problems • Hamiltonian (Hybrid) Monte Carlo (HMC) • HMC in Practice and Theory Monte Carlo Applications to Molecular Systems Introduction Outline • Metropolis Algorithm • Simulated Annealing • Parallel Tempering • Hybrid Monte Carlo • References Monte Carlo Applications to Molecular Systems (I) Some well known Monte Carlo based algorithms were initially designed as applications to molecular systems: Metropolis algorithm, Simulated Annealing, Hybrid Monte Carlo, Parallel Tempering. • Metropolis algorithm [1] was the first MCMC algorithm ü It is associated with a second computer called MANIAC built in Los Alamos under the direction of a physicist and a mathematician Nicolas Metropolis in early 1952 ü The method was designed for simulation of a liquid in equilibrium with its gas phase ü The idea of the method published in 1953 in the Journal of Chemical Physics was to create a Markov chain (a random walk in thermodynamic phase space) of molecular configurations whose appearance is proportional to their statistical probability ü Metropolis algorithm: in the year 2000 was ranked as one of the “10 algorithms with the greatest influence on the development and practice of science and engineering in the 20th century”[2].
    [Show full text]