Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene and carbon tetrachloride

Size: px
Start display at page:

Download "Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene and carbon tetrachloride"

Transcription

1 - 1 - PLMMP 2016, Kyiv Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene and carbon tetrachloride Jadran Vrabec Tatjana Janzen, Gabriela Guevara-Carrión, Yonny Mauricio Munoz Munoz () Paderborn University, Germany

2 - 2 - Motivation Diffusion coefficients are essential for modelling of complex systems and processes in science and engineering Problems: - Only few experimental data available, especially for multicomponent mixtures - Measurements are challenging and time consuming - Theoretical and empirical approaches often fail in predictions for strongly non-ideal mixtures Goal: - Establish precise predictive methods Prediction of diffusion coefficients by molecular simulation

3 - 3 - Molecular simulation Calculation of macroscopic behavior from molecular interactions - Molecular models based on force fields - Equilibrium molecular dynamics - Force calculations to get molecular trajectories static and dynamic fluid properties - Calculation of transport properties Green-Kubo formalism ms2: simulation tool for thermodynamic properties Glass et al., Comp. Phys. Commun. 185 (2014) 3302

4 ms2: simulation tool for thermodynamic properties Molecular dynamics / Monte Carlo Arbitrary mixtures of rigid molecules Grand equilibrium method (for VLE) Several classical ensembles All static properties (thermal, caloric, entropic) Gradual insertion for entropic properties Transport properties (Green-Kubo) Consistent FORTRAN90 code Object oriented All loops vectorized MD and MC parallelized 3D visualization interface Glass et al., Comp. Phys. Commun. 185 (2014)

5 - 5 - Description of diffusion (1) nn 1 Fick: JJ ii = ρρ mm DD iiii jj=1 xx jj coefficients from experiment Maxwell-Stefan: nn xx jj uu ii uu jj jj ii=1 Đ iiii = 1 RRRR μμ ii coefficients from molecular simulation xx ii uu ii = 1 RRRR Λ iiii μμ jj nn jj=1 Onsager reciprocal relation: Λ iiii = Λ jjii Binary mixture: Đ iiii = xx jj xx ii Λ iiii + xx ii xx jj Λ jjjj Λ iijj Λ jjii Fick Maxwell-Stefan

6 - 6 - Description of diffusion (2) Fick Maxwell-Stefan Binary mixture: DD = ĐГ Г = 1 + xx 1 dd ln γγ 1 ddxx 1 Thermodynamic factor Г Describes the thermodynamic non-ideality of a mixture Ideal mixture: Г = 1 Thermodynamic instability (phase separation): Г < 0 Can be calculated from G E models (e.g. Wilson, NRTL, UNIQUAC) Fitting of the model parameters to experimental VLE data or to simulation data

7 - 7 - Green-Kubo formalism (1) Microscopic fluctuations around equilibrium Description of non-equilibrium phenomena Transport coefficients from time dependent autocorrelation functions of corresponding flux DD ii, DD iijj, ηη, λλ e.g. self-diffusion NN DD ii = 1 dddd vv 3NN ii ii 0 ii (0) vv ii (tt)

8 - 8 - Green-Kubo formalism (2) Transport coefficients: DD ii, DD iijj, ηη, λλ Self-diffusion DD ii = 1 NN 3NN ii dddd vv ii 0 ii (0) vv ii (tt) Onsager coefficients mutual diffusion DD iiii Shear viscosity Thermal conductivity NN ii Λ iiii = 1 3NN dddd vv ii,kk 0 vv jj,ll tt 0 kk=1 NN jj ll=1 ηη = 1 VVkk BB TT dddd JJ xxxx pp (0) JJ xxxx pp (tt) 0 λλ = 1 VVkk BB TT 2 dddd JJ qq xx (0) JJ xx qq (tt) 0

9 - 9 - Molecular models Methanol CCl 4 Ethanol Acetone - Rigid molecules (united atom) - Lennard-Jones sites, point charges, dipole, quadrupole - Parameters optimized to saturated liquid density and vapor pressure, (self-diffusion) - Mixing behavior: predicted Toluene Benzene Cyclohexane

10 Molecular models: new models Benzene Toluene CCl 4

11 Studied mixtures (1) Methanol CCl 4 Toluol Ethanol Aceton Benzol Cyclohexan Molecular models - Rigid molecules (united atom) - Lennard-Jones sites, point charges, dipole, quadrupole - Parameters optimized to saturated liquid density and vapor pressure, (self-diffusion) - Mixing behavior: predicted 20 binary mixtures

12 Studied mixtures (2) Groups according to deviation of thermodynamic factor from ideal behavior Group I (less 10 %) Group II (10 %...45 %) Group III (>60 %) Ethanol + Methanol Benzene + Toluene Benzene + CCl 4 Cyclohexane + CCl 4 Toluene + CCl 4 Methanol + Acetone Ethanol + Acetone Acetone + Benzene Acetone + Toluene Acetone + CCl 4 Benzene + Cyclohexane Cyclohexane + Toluene Methanol + Benzene Methanol + Toluene Methanol + CCl 4 Ethanol + Benzene Ethanol + Cyclohexane Ethanol + Toluene Ethanol + CCl 4 All mixtures studied at T = K and p = 0.1 MPa Acetone + Cyclohexane

13 Predictive phenomenological models

14 Group I: Benzene + Toluene Self-diffusion Diffusion coefficients Shear viscosity Thermal conductivity

15 Group I: Fick diffusion Ethanol + Methanol Benzene + CCl 4 Cyclohexane + CCl 4 Toluene + CCl 4

16 Group II: Acetone + Benzene Self-diffusion Diffusion coefficients Shear viscosity Thermal conductivity

17 Group II: Fick diffusion Methanol + Acetone Ethanol + Acetone Acetone + Toluene Acetone + CCl 4 Benzene + Cyclohexane Cyclohexane + Toluene

18 Group III: Methanol + Toluene Self-diffusion Diffusion coefficients Shear viscosity Thermal conductivity

19 Group III: Fick diffusion Methanol + Benzene Methanol + CCl 4 Ethanol + Benzene Ethanol + Cyclohexane Ethanol + Toluene Acetone + Cyclohexane

20 Transport properties: average deviation Mean average deviation for all mixtures Fick diffusion: 16 % (calculated with Wilson model)

21 Radial distribution functions G I: Benzene + Toluene G II: Acetone + Benzene G III: Methanol + Toluene

22 Radial distribution functions - Hydrogen bonding in mixtures with ethanol and methanol - Strong non-idealities in mixtures with one polar and one non-polar component - RDF show that extensive clustering occurs in mixtures of group III - Micro-heterogeneity is present in these mixtures clusters can be quite large, but no phase separation - Even stronger thermodynamic non-idealities lead to liquid-liquid equilibrium e.g. Methanol + Cyclohexane

23 Methanol + Benzene x 1 = 0.1 mol/mol

24 Methanol + Benzene x 1 = 0.1 mol/mol

25 Summary Method: - Prediction of transport coefficient by molecular simulation - Equilibrium MD, Green-Kubo formalism Results: - Mutual diffusion in different binary liquid mixtures - Further transport properties - Good predictons of all kind of mixtures Outlook: - Diffusion coefficients in mixtures with LLE - Diffusion coefficients in ternary mixtures

Mutual diffusion in the ternary mixture of water + methanol + ethanol: Experiments and Molecular Simulation

Mutual diffusion in the ternary mixture of water + methanol + ethanol: Experiments and Molecular Simulation - 1 - Mutual diffusion in the ternary mixture of water + methanol + ethanol: Experiments and Molecular Simulation Tatjana Janzen, Gabriela Guevara-Carrión, Jadran Vrabec University of Paderborn, Germany

More information

Atomistic molecular simulations for engineering applications: methods, tools and results. Jadran Vrabec

Atomistic molecular simulations for engineering applications: methods, tools and results. Jadran Vrabec Atomistic molecular simulations for engineering applications: methods, tools and results Jadran Vrabec Motivation Simulation methods vary in their level of detail The more detail, the more predictive power

More information

Lennard-Jones force field parameters for cyclic alkanes from cyclopropane to cyclohexane

Lennard-Jones force field parameters for cyclic alkanes from cyclopropane to cyclohexane Lennard-Jones force field parameters for cyclic alkanes from cyclopropane to cyclohexane Y. Mauricio Muñoz-Muñoz a,b, Gabriela Guevara-Carrion b, Mario Llano-Restrepo a, Jadran Vrabec* b (a) School of

More information

Molecular Modeling and Simulation of Phase Equilibria for Chemical Engineering

Molecular Modeling and Simulation of Phase Equilibria for Chemical Engineering InPROMT 2012, Berlin, 16. November 2012 DFG Transregio CRC 63 Molecular Modeling and Simulation of Phase Equilibria for Chemical Engineering Hans Hasse 1, Martin Horsch 1, Jadran Vrabec 2 1 Laboratory

More information

Molecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties

Molecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties 12 th HLRS Results and Review Workshop Molecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties High Performance Computing Center Stuttgart (HLRS), October 8,

More information

Effects of concentration and temperature on the coupled heat and mass transport in liquid mixtures

Effects of concentration and temperature on the coupled heat and mass transport in liquid mixtures University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Papers in Thermal Mechanics Chemical and Biomolecular Engineering Research and Publications 3-26-2001 Effects of concentration

More information

Heat, Work, and the First Law of Thermodynamics. Chapter 18 of Essential University Physics, Richard Wolfson, 3 rd Edition

Heat, Work, and the First Law of Thermodynamics. Chapter 18 of Essential University Physics, Richard Wolfson, 3 rd Edition Heat, Work, and the First Law of Thermodynamics Chapter 18 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 Different ways to increase the internal energy of system: 2 Joule s apparatus

More information

Molecular Dynamic Simulation Study of the Volume Transition of PNIPAAm Hydrogels

Molecular Dynamic Simulation Study of the Volume Transition of PNIPAAm Hydrogels Molecular Dynamic Simulation Study of the Volume Transition of PNIPAAm Hydrogels Jonathan Walter 1, Jadran Vrabec 2, Hans Hasse 1 1 Laboratory of Engineering, University of Kaiserslautern, Germany 2 and

More information

Thermodynamics, Design, Simulation and Benchmarking of Biofuel Processes

Thermodynamics, Design, Simulation and Benchmarking of Biofuel Processes Thermodynamics, Design, Simulation and Benchmarking of Biofuel Processes Mauro Torli Philip Loldrup Fosbøl Georgios Kontogeorgis SYNFERON project work packages Commercial syngas fermentation technologies:

More information

EXPERIMENT 1 DETERMINATION OF GAS DIFFUSION COEFFICIENT

EXPERIMENT 1 DETERMINATION OF GAS DIFFUSION COEFFICIENT EXPERIMENT 1 DETERMINATION OF GAS DIFFUSION COEFFICIENT Objective: The objective of this experiment is to calculate diffusion coefficient of a volatile organic compound in air by means of Chapman Enskog

More information

Diffusion modeling for Dip-pen Nanolithography Apoorv Kulkarni Graduate student, Michigan Technological University

Diffusion modeling for Dip-pen Nanolithography Apoorv Kulkarni Graduate student, Michigan Technological University Diffusion modeling for Dip-pen Nanolithography Apoorv Kulkarni Graduate student, Michigan Technological University Abstract The diffusion model for the dip pen nanolithography is similar to spreading an

More information

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Review for Exam2 11. 13. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University

More information

Contents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21

Contents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21 Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic

More information

Name: Discussion Section:

Name: Discussion Section: CBE 141: Chemical Engineering Thermodynamics, Spring 2018, UC Berkeley Midterm 2 March 22, 2018 Time: 80 minutes, closed-book and closed-notes, one-sided 8 ½ x 11 equation sheet allowed Please show all

More information

Prediction of methanol content in natural gas with the GC-PR-CPA model Hajiw, Martha; Chapoy, Antonin; Coquelet, Christophe; Lauermann, Gerhard

Prediction of methanol content in natural gas with the GC-PR-CPA model Hajiw, Martha; Chapoy, Antonin; Coquelet, Christophe; Lauermann, Gerhard Heriot-Watt University Heriot-Watt University Research Gateway Prediction of methanol content in natural gas with the GC-PR-CPA model Hajiw, Martha; Chapoy, Antonin; Coquelet, Christophe; Lauermann, Gerhard

More information

Angular Momentum, Electromagnetic Waves

Angular Momentum, Electromagnetic Waves Angular Momentum, Electromagnetic Waves Lecture33: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay As before, we keep in view the four Maxwell s equations for all our discussions.

More information

7.3 The Jacobi and Gauss-Seidel Iterative Methods

7.3 The Jacobi and Gauss-Seidel Iterative Methods 7.3 The Jacobi and Gauss-Seidel Iterative Methods 1 The Jacobi Method Two assumptions made on Jacobi Method: 1.The system given by aa 11 xx 1 + aa 12 xx 2 + aa 1nn xx nn = bb 1 aa 21 xx 1 + aa 22 xx 2

More information

MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES

MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES V.V. Korenkov 1,3, a, E.G. Nikonov 1, b, M. Popovičová 2, с 1 Joint Institute for Nuclear

More information

Lecture 3. STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher

Lecture 3. STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher Lecture 3 STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher Previous lectures What is machine learning? Objectives of machine learning Supervised and

More information

Review for Exam Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa

Review for Exam Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa 57:020 Fluids Mechanics Fall2013 1 Review for Exam3 12. 11. 2013 Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa 57:020 Fluids Mechanics Fall2013 2 Chapter

More information

Worksheets for GCSE Mathematics. Algebraic Expressions. Mr Black 's Maths Resources for Teachers GCSE 1-9. Algebra

Worksheets for GCSE Mathematics. Algebraic Expressions. Mr Black 's Maths Resources for Teachers GCSE 1-9. Algebra Worksheets for GCSE Mathematics Algebraic Expressions Mr Black 's Maths Resources for Teachers GCSE 1-9 Algebra Algebraic Expressions Worksheets Contents Differentiated Independent Learning Worksheets

More information

(1) Correspondence of the density matrix to traditional method

(1) Correspondence of the density matrix to traditional method (1) Correspondence of the density matrix to traditional method New method (with the density matrix) Traditional method (from thermal physics courses) ZZ = TTTT ρρ = EE ρρ EE = dddd xx ρρ xx ii FF = UU

More information

Impact of Chamber Pressure on Sputtered Particle Energy

Impact of Chamber Pressure on Sputtered Particle Energy Wilmert De Bosscher Chief Technology Officer +32 9381 6177 wilmert.debosscher@soleras.com Impact of Chamber Pressure on Sputtered Particle Energy Tampa, October 18 th, 2017 Background Why Sputtering at

More information

Self-wrapping of an Ouzo Drop Induced by Evaporation on a Superamphiphobic Surface Supplementary Material

Self-wrapping of an Ouzo Drop Induced by Evaporation on a Superamphiphobic Surface Supplementary Material Electronic Supplementary Material (ESI) for Soft Matter. This journal is The Royal Society of Chemistry 2017 Self-wrapping of an Ouzo Drop Induced by Evaporation on a Superamphiphobic Surface Supplementary

More information

CE 530 Molecular Simulation

CE 530 Molecular Simulation 1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy

More information

Review of Last Class 1

Review of Last Class 1 Review of Last Class 1 X-Ray diffraction of crystals: the Bragg formulation Condition of diffraction peak: 2dd sin θθ = nnλλ Review of Last Class 2 X-Ray diffraction of crystals: the Von Laue formulation

More information

Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany

Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics

More information

A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions

A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin A Posteriori DG using Non-Polynomial Basis 1 A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin Department of Mathematics, UC Berkeley;

More information

Molecular modelling and simulation of hydrogen bonding fluids

Molecular modelling and simulation of hydrogen bonding fluids Molecular modelling and simulation of hydrogen bonding fluids Martin Horsch, Alptekin Celik and Hans Hasse Lehrstuhl für Thermodynamik (LTD) Technische Universität Kaiserslautern Rostock, 27 th March 13

More information

Materials & Advanced Manufacturing (M&AM)

Materials & Advanced Manufacturing (M&AM) Modeling of Shear Thickening Fluids for Analysis of Energy Absorption Under Impulse Loading Alyssa Bennett (University of Michigan) Nick Vlahopoulos, PhD (University of Michigan) Weiran Jiang, PhD (Research

More information

Thermodynamic Cycles

Thermodynamic Cycles Thermodynamic Cycles Content Thermodynamic Cycles Carnot Cycle Otto Cycle Rankine Cycle Refrigeration Cycle Thermodynamic Cycles Carnot Cycle Derivation of the Carnot Cycle Efficiency Otto Cycle Otto Cycle

More information

Gradient expansion formalism for generic spin torques

Gradient expansion formalism for generic spin torques Gradient expansion formalism for generic spin torques Atsuo Shitade RIKEN Center for Emergent Matter Science Atsuo Shitade, arxiv:1708.03424. Outline 1. Spintronics a. Magnetoresistance and spin torques

More information

Coarse-Graining via the Mori-Zwanzig Formalism

Coarse-Graining via the Mori-Zwanzig Formalism Coarse-Graining via the Mori-Zwanzig Formalism George Em Karniadakis & Zhen Li http://www.pnnl.gov/computing/cm4/ Supported by DOE ASCR Mori-Zwanzig can be used to: 1. Derive coarse-grained equations 2.

More information

Lecture 11. Kernel Methods

Lecture 11. Kernel Methods Lecture 11. Kernel Methods COMP90051 Statistical Machine Learning Semester 2, 2017 Lecturer: Andrey Kan Copyright: University of Melbourne This lecture The kernel trick Efficient computation of a dot product

More information

J. Chem. Thermodynamics

J. Chem. Thermodynamics J. Chem. Thermodynamics 43 (2011) 1192 1203 Contents lists available at ScienceDirect J. Chem. Thermodynamics journal homepage: www.elsevier.com/locate/jct Form of multicomponent Fickian diffusion coefficients

More information

Engineering Molecular Models: Efficient Parameterization Procedure and Cyclohexanol as Case Study

Engineering Molecular Models: Efficient Parameterization Procedure and Cyclohexanol as Case Study Engineering Molecular Models: Efficient Parameterization Procedure and Cyclohexanol as Case Study Thorsten Merker, Jadran Vrabec, and Hans Hasse Laboratory of Engineering Thermodynamics, University of

More information

Charge carrier density in metals and semiconductors

Charge carrier density in metals and semiconductors Charge carrier density in metals and semiconductors 1. Introduction The Hall Effect Particles must overlap for the permutation symmetry to be relevant. We saw examples of this in the exchange energy in

More information

Part I.

Part I. Part I bblee@unimp . Introduction to Mass Transfer and Diffusion 2. Molecular Diffusion in Gasses 3. Molecular Diffusion in Liquids Part I 4. Molecular Diffusion in Biological Solutions and Gels 5. Molecular

More information

CH2351 Chemical Engineering Thermodynamics II Unit I, II Phase Equilibria. Dr. M. Subramanian

CH2351 Chemical Engineering Thermodynamics II Unit I, II   Phase Equilibria.   Dr. M. Subramanian CH2351 Chemical Engineering Thermodynamics II Unit I, II Phase Equilibria Dr. M. Subramanian Associate Professor Department of Chemical Engineering Sri Sivasubramaniya Nadar College of Engineering Kalavakkam

More information

Scalable, performant, and resilient large-scale applications of molecular process engineering

Scalable, performant, and resilient large-scale applications of molecular process engineering Scalable, performant, and resilient large-scale applications of molecular process engineering M. Horsch,1 P. Gralka,2 C. Niethammer,3 N. Tchipev,4 J. Vrabec,5 H. Hasse1 1 University of Kaiserslautern,

More information

Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons

Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons Gang Chen Massachusetts Institute of Technology OXFORD UNIVERSITY PRESS 2005 Contents Foreword,

More information

1.3 Molecular Level Presentation

1.3 Molecular Level Presentation 1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of

More information

Molecular Dynamics Simulation of Methanol-Water Mixture

Molecular Dynamics Simulation of Methanol-Water Mixture Molecular Dynamics Simulation of Methanol-Water Mixture Palazzo Mancini, Mara Cantoni University of Urbino Carlo Bo Abstract In this study some properties of the methanol-water mixture such as diffusivity,

More information

Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane

Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.6, pp 2975-2979, Oct-Dec 2013 Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane +

More information

CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions

CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions The aim of solution theory is to express the properties of liquid mixture in terms of intermolecular forces and liquid structure. The

More information

Jasmin Smajic1, Christian Hafner2, Jürg Leuthold2, March 23, 2015

Jasmin Smajic1, Christian Hafner2, Jürg Leuthold2, March 23, 2015 Jasmin Smajic, Christian Hafner 2, Jürg Leuthold 2, March 23, 205 Time Domain Finite Element Method (TD FEM): Continuous and Discontinuous Galerkin (DG-FEM) HSR - University of Applied Sciences of Eastern

More information

424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393

424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393 Index After-effect function, 368, 369 anthropic principle, 232 assumptions nature of, 242 autocorrelation function, 292 average, 18 definition of, 17 ensemble, see ensemble average ideal,23 operational,

More information

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Review for Exam3 12. 9. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University

More information

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore

More information

Radiation. Lecture40: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay

Radiation. Lecture40: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Radiation Zone Approximation We had seen that the expression for the vector potential for a localized cuent distribution is given by AA (xx, tt) = μμ 4ππ ee iiiiii dd xx eeiiii xx xx xx xx JJ (xx ) In

More information

Lecture No. 5. For all weighted residual methods. For all (Bubnov) Galerkin methods. Summary of Conventional Galerkin Method

Lecture No. 5. For all weighted residual methods. For all (Bubnov) Galerkin methods. Summary of Conventional Galerkin Method Lecture No. 5 LL(uu) pp(xx) = 0 in ΩΩ SS EE (uu) = gg EE on ΓΓ EE SS NN (uu) = gg NN on ΓΓ NN For all weighted residual methods NN uu aaaaaa = uu BB + αα ii φφ ii For all (Bubnov) Galerkin methods ii=1

More information

Exam 2 Fall 2015

Exam 2 Fall 2015 1 95.144 Exam 2 Fall 2015 Section instructor Section number Last/First name Last 3 Digits of Student ID Number: Show all work. Show all formulas used for each problem prior to substitution of numbers.

More information

PH4211 Statistical Mechanics Brian Cowan

PH4211 Statistical Mechanics Brian Cowan PH4211 Statistical Mechanics Brian Cowan Contents 1 The Methodology of Statistical Mechanics 1.1 Terminology and Methodology 1.1.1 Approaches to the subject 1.1.2 Description of states 1.1.3 Extensivity

More information

Introduction Statistical Thermodynamics. Monday, January 6, 14

Introduction Statistical Thermodynamics. Monday, January 6, 14 Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can

More information

Development and Evaluation of a New Equation of State for Polar Fluids: Pure Components and Mixtures

Development and Evaluation of a New Equation of State for Polar Fluids: Pure Components and Mixtures Development and Evaluation of a New Equation of State for Polar Fluids: Pure Components and Mixtures PhD thesis by Eirini Karakatsani National Center for Scientific Research Demokritos, Aghia Paraskevi

More information

Revision : Thermodynamics

Revision : Thermodynamics Revision : Thermodynamics Formula sheet Formula sheet Formula sheet Thermodynamics key facts (1/9) Heat is an energy [measured in JJ] which flows from high to low temperature When two bodies are in thermal

More information

2018 Fall Semester Quiz 4 For General Chemistry I (CH101)

2018 Fall Semester Quiz 4 For General Chemistry I (CH101) 2018 Fa Semester Quiz 4 For General Chemistry I (CH101) Date: November 5 (Mon), Time: 19:00 ~ 19:45 Professor Name Class Student I.D. Number Name Solution Avogadro constant (N A): 6.022 x 10 23 mol -1

More information

Molecular Driving Forces

Molecular Driving Forces Molecular Driving Forces Statistical Thermodynamics in Chemistry and Biology SUBGfittingen 7 At 216 513 073 / / Ken A. Dill Sarina Bromberg With the assistance of Dirk Stigter on the Electrostatics chapters

More information

Creeping Movement across a Long Strike-Slip Fault in a Half Space of Linear Viscoelastic Material Representing the Lithosphere-Asthenosphere System

Creeping Movement across a Long Strike-Slip Fault in a Half Space of Linear Viscoelastic Material Representing the Lithosphere-Asthenosphere System Frontiers in Science 214, 4(2): 21-28 DOI: 1.5923/j.fs.21442.1 Creeping Movement across a Long Strike-Slip Fault in a Half Space of Linear Viscoelastic Material Representing the Lithosphere-Asthenosphere

More information

Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses

Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses A normal liquid will crystallize at a melting temperature T m as it is cooled via a first-order phase transition (see figure above).

More information

CHEMISTRY - TRO 4E CH.11 - LIQUIDS, SOLIDS & INTERMOLECULAR FORCES

CHEMISTRY - TRO 4E CH.11 - LIQUIDS, SOLIDS & INTERMOLECULAR FORCES !! www.clutchprep.com CONCEPT: INTERMOLECULAR FORCES When looking at a molecular substance such as H 2 O you will discover two types of electrostatic forces at work: forces exist within a molecule and

More information

G : Statistical Mechanics

G : Statistical Mechanics G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble

More information

Modified Raoult's Law and Excess Gibbs Free Energy

Modified Raoult's Law and Excess Gibbs Free Energy ACTIVITY MODELS 1 Modified Raoult's Law and Excess Gibbs Free Energy Equilibrium criteria: f V i = L f i For vapor phase: f V i = y i i P For liquid phase, we may use an activity coefficient ( i ), giving:

More information

Interparticle interaction

Interparticle interaction Colloid chemistry for pharmacists Interparticle interaction Levente Novák, István Bányai Dep. of Colloid- and Environmental Chemistry www.kolloid.unideb.hu/ 2. lecture 1 Characterization of colloids colloidal

More information

NON-EQUILIBRIUM THERMODYNAMICS

NON-EQUILIBRIUM THERMODYNAMICS NON-EQUILIBRIUM THERMODYNAMICS S. R. DE GROOT Professor of Theoretical Physics University of Amsterdam, The Netherlands E MAZUR Professor of Theoretical Physics University of Leiden, The Netherlands DOVER

More information

Properties of Solutions. Review

Properties of Solutions. Review Properties of Solutions Review Matter Pure substance Mixture of substances compound element homogeneous heterogeneous Solution Definitions A solution is a homogeneous mixture of two or more substances.

More information

An Introduction to Two Phase Molecular Dynamics Simulation

An Introduction to Two Phase Molecular Dynamics Simulation An Introduction to Two Phase Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 19, 2016 date last updated: April

More information

THERMODYNAMICS THERMOSTATISTICS AND AN INTRODUCTION TO SECOND EDITION. University of Pennsylvania

THERMODYNAMICS THERMOSTATISTICS AND AN INTRODUCTION TO SECOND EDITION. University of Pennsylvania THERMODYNAMICS AND AN INTRODUCTION TO THERMOSTATISTICS SECOND EDITION HERBERT B. University of Pennsylvania CALLEN JOHN WILEY & SONS New York Chichester Brisbane Toronto Singapore CONTENTS PART I GENERAL

More information

A Generalized Correlation for Pool Boiling Heat Transfer Coefficient Based on Corresponding State Rule for Pure Compounds and Binary Mixtures

A Generalized Correlation for Pool Boiling Heat Transfer Coefficient Based on Corresponding State Rule for Pure Compounds and Binary Mixtures A Generalized Correlation for Pool Boiling Heat Transfer Coefficient Based on Corresponding State Rule for Pure Compounds and Binary Mixtures HASAN QABAZARD and MAHMOOD MOSHFEGHIAN 1 Petroleum Research

More information

M.5 Modeling the Effect of Functional Responses

M.5 Modeling the Effect of Functional Responses M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,

More information

Computer simulations as concrete models for student reasoning

Computer simulations as concrete models for student reasoning Computer simulations as concrete models for student reasoning Jan Tobochnik Department of Physics Kalamazoo College Kalamazoo MI 49006 In many thermal physics courses, students become preoccupied with

More information

Time-Dependent Statistical Mechanics 1. Introduction

Time-Dependent Statistical Mechanics 1. Introduction Time-Dependent Statistical Mechanics 1. Introduction c Hans C. Andersen Announcements September 24, 2009 Lecture 1 9/22/09 1 Topics of concern in the course We shall be concerned with the time dependent

More information

BAE 820 Physical Principles of Environmental Systems

BAE 820 Physical Principles of Environmental Systems BAE 820 Physical Principles of Environmental Systems Estimation of diffusion Coefficient Dr. Zifei Liu Diffusion mass transfer Diffusion mass transfer refers to mass in transit due to a species concentration

More information

Variations. ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra

Variations. ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra Variations ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra Last Time Probability Density Functions Normal Distribution Expectation / Expectation of a function Independence Uncorrelated

More information

Quantum state measurement

Quantum state measurement Quantum state measurement Introduction The rotation properties of light fields of spin are described by the 3 3 representation of the 0 0 SO(3) group, with the generators JJ ii we found in class, for instance

More information

Advanced data analysis

Advanced data analysis Advanced data analysis Akisato Kimura ( 木村昭悟 ) NTT Communication Science Laboratories E-mail: akisato@ieee.org Advanced data analysis 1. Introduction (Aug 20) 2. Dimensionality reduction (Aug 20,21) PCA,

More information

Thermodynamic behaviour of mixtures containing CO 2. A molecular simulation study

Thermodynamic behaviour of mixtures containing CO 2. A molecular simulation study Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues

More information

Fluctuating Hydrodynamics for Transport in Neutral and Charged Fluid Mixtures

Fluctuating Hydrodynamics for Transport in Neutral and Charged Fluid Mixtures Fluctuating Hydrodynamics for Transport in Neutral and Charged Fluid Mixtures Alejandro L. Garcia San Jose State University Lawrence Berkeley Nat. Lab. Hydrodynamic Fluctuations in Soft Matter Simulations

More information

Jacco Snoeijer PHYSICS OF FLUIDS

Jacco Snoeijer PHYSICS OF FLUIDS Jacco Snoeijer PHYSICS OF FLUIDS dynamics dynamics freezing dynamics freezing microscopics of capillarity Menu 1. surface tension: thermodynamics & microscopics 2. wetting (statics): thermodynamics & microscopics

More information

Analog Circuits Part 1 Circuit Theory

Analog Circuits Part 1 Circuit Theory Introductory Medical Device Prototyping Analog Circuits Part 1 Circuit Theory, http://saliterman.umn.edu/ Department of Biomedical Engineering, University of Minnesota Concepts to be Covered Circuit Theory

More information

CE 530 Molecular Simulation

CE 530 Molecular Simulation 1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops

More information

Free energy dependence along the coexistence curve

Free energy dependence along the coexistence curve Free energy dependence along the coexistence curve In a system where two phases (e.g. liquid and gas) are in equilibrium the Gibbs energy is G = GG l + GG gg, where GG l and GG gg are the Gibbs energies

More information

Statistical Mechanics

Statistical Mechanics Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2

More information

Colloidal Dispersions

Colloidal Dispersions Physical Pharmacy Lecture 5 Colloidal Dispersions Assistant Lecturer in Pharmaceutics Overview Dispersed Systems Classification Colloidal Systems Properties of Colloids Optical Properties Kinetic Properties

More information

The simultaneous prediction of vapor-liquid equilibrium and excess enthalpy. Kwon, Jung Hun. Thermodynamics and properties lab.

The simultaneous prediction of vapor-liquid equilibrium and excess enthalpy. Kwon, Jung Hun. Thermodynamics and properties lab. The simultaneous prediction of vapor-liquid equilibrium and excess enthalpy Kwon, Jung Hun. 2 Contents 1 A comparison of cubic EOS mixing rules for the simultaneous description of excess enthalpies and

More information

Overall: 75 ECTS: 7.0

Overall: 75 ECTS: 7.0 Course: Chemical Engineering Thermodynamics Language: English Lecturer: Prof. dr. sc. Marko Rogošić TEACHING WEEKLY SEMESTER Lectures 3 45 Laboratory 1 15 Seminar 1 15 Overall: 75 ECTS: 7.0 PURPOSE: Within

More information

Phase equilibrium studies of impure CO 2 systems to underpin developments of CCS technologies

Phase equilibrium studies of impure CO 2 systems to underpin developments of CCS technologies Phase equilibrium studies of impure CO 2 systems to underpin developments of CCS technologies Jie Ke, Martyn Poliakoff and Michael W. George School of Chemistry The University of Nottingham 1 September,

More information

Transport phenomena in a temperature gradient studied by NEMD. A chemical reaction and a phase transition

Transport phenomena in a temperature gradient studied by NEMD. A chemical reaction and a phase transition Jing Xu Transport phenomena in a temperature gradient studied by NEMD. A chemical reaction and a phase transition Thesis for the degree of Philosophiae Doctor Trondheim, October 2009 Norwegian University

More information

Constitutive Relations of Stress and Strain in Stochastic Finite Element Method

Constitutive Relations of Stress and Strain in Stochastic Finite Element Method American Journal of Computational and Applied Mathematics 15, 5(6): 164-173 DOI: 1.593/j.ajcam.1556. Constitutive Relations of Stress and Strain in Stochastic Finite Element Method Drakos Stefanos International

More information

Fall Possibly Useful Information: 1 atm = lb/in 2 = kpa. 1 atm = 101,325 N/m 2 = 760 mmhg. 1 atm = 101,325 Pa = 1.

Fall Possibly Useful Information: 1 atm = lb/in 2 = kpa. 1 atm = 101,325 N/m 2 = 760 mmhg. 1 atm = 101,325 Pa = 1. Chemistry 122 (Tyvoll) Fall 2005 PRACTICE EXAMINATION I Possibly Useful Information: 1 atm = 14.70 lb/in 2 = 101.325 kpa 1 atm = 101,325 N/m 2 = 760 mmg 1 atm = 101,325 Pa = 1.01325 bar 1 atm = 1013.25

More information

Linear relation between Heisenberg exchange and interfacial Dzyaloshinskii Moriya interaction in metal films

Linear relation between Heisenberg exchange and interfacial Dzyaloshinskii Moriya interaction in metal films Linear relation between Heisenberg exchange and interfacial Dzyaloshinskii Moriya interaction in metal films Hans T. Nembach, Justin M. Shaw, Mathias Weiler*, Emilie Jué and Thomas J. Silva Electromagnetics

More information

Simple lattice-gas model for water

Simple lattice-gas model for water Simple lattice-gas model for water A. Ciach, W. T. Góźdź and A. Perera Motivation H-bonds and packing of molecules Construction of the lattice model Results in mean-field approximation Water molecule and

More information

THE PROPERTIES OF GASES AND LIQUIDS

THE PROPERTIES OF GASES AND LIQUIDS THE PROPERTIES OF GASES AND LIQUIDS Bruce E. Poling University of Toledo John M. Prausnitz University of California at Berkeley John P. O'Connell University of Virginia Fifth Edition McGRAW-HILL New York

More information

CH 240 Chemical Engineering Thermodynamics Spring 2007

CH 240 Chemical Engineering Thermodynamics Spring 2007 CH 240 Chemical Engineering Thermodynamics Spring 2007 Instructor: Nitash P. Balsara, nbalsara@berkeley.edu Graduate Assistant: Paul Albertus, albertus@berkeley.edu Course Description Covers classical

More information

TECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES

TECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES COMPUTERS AND STRUCTURES, INC., FEBRUARY 2016 TECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES Introduction This technical note

More information

General Chemistry A

General Chemistry A General Chemistry 1140 - A May 6, 2004 (6 Pages, 43 Parts) Name Each of the 40 multiple choice questions counts 2 point. Give the letter of the correct answer. 1. 2. Crystalline solids differ from amorphous

More information

Development new correlations for NRTL Parameters in Polymer Solutions

Development new correlations for NRTL Parameters in Polymer Solutions Rhodes, Greece, August 0-, 008 Development new correlations for NRTL Parameters in Polymer Solutions A. Saatchi, M. Edalat* Oil and Gas Center of Excellence, Department of Chemical Engineering University

More information

Molecular simulation of adsorption from dilute solutions

Molecular simulation of adsorption from dilute solutions Vol. 52 No. 3/2005 685 689 on-line at: www.actabp.pl Molecular simulation of adsorption from dilute solutions Werner Billes Rupert Tscheliessnig and Johann Fischer Institut für Verfahrens- und Energietechnik

More information

Simulation of Ethanol Dehydration Using Cyclohexane as an Entrainer

Simulation of Ethanol Dehydration Using Cyclohexane as an Entrainer Simulation of Ethanol Dehydration Using Cyclohexane as an Entrainer 2013 년 9 월 23 일 ( 월 ) 공주대학교화학공학부조정호 Limit of Distillation by Azeotrope 1.0 Ethanol / Water System Distillation range is restricted by

More information

Outline of the Course

Outline of the Course Outline of the Course 1) Review and Definitions 2) Molecules and their Energies 3) 1 st Law of Thermodynamics Conservation of Energy. 4) 2 nd Law of Thermodynamics Ever-Increasing Entropy. 5) Gibbs Free

More information