Atomistic molecular simulations for engineering applications: methods, tools and results. Jadran Vrabec
|
|
- Maximillian Daniel
- 5 years ago
- Views:
Transcription
1 Atomistic molecular simulations for engineering applications: methods, tools and results Jadran Vrabec
2 Motivation Simulation methods vary in their level of detail The more detail, the more predictive power Quantum chemical methods scale with O(M 3 ) O(M 7 ) [M: basis functions] Force field methods are favorable with respect to scaling and thermodynamics
3 Force fields Thermodynamics Classical models for molecular interactions Parameters have a physical interpretation Contain all thermodynamic properties Static: thermal, caloric, entropic Dynamic: viscosity, diffusion, thermal conductivity, Surface properties, e.g. surface tension Iso-Butane Straightforwardly applicable to mixtures Excellent predictive power Technical accuracy Directly applicable for the study of fluids in geometries, e.g. wetting, adsorption, zeolites, in processes, e.g. condensation, flow,
4 Sampling force fields Molecular dynamics Monte Carlo
5 ms2: simulation tool for thermodynamic properties Deublein et al., Comp. Phys. Commun. 182 (2011) 2350 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
6 Equation of state for CO 2 (Span and Wagner, 1996) T = K, p = MPa F R T 0 Res,,, τ =T c / T δ = ρ / ρ c Ideal part: Residual part: i i i 4, ln a a a ln a ln 1 exp n 7 34 Re s t i di ti di ci, ai ai exp i 1 i 8 39 ti di 2 2 ai exp i( i) i( i) i bi 2 2 ai exp C i( 1) D i( 1) i 40
7 Derivatives of the Massieu-Planck potential F /(RT ) nm n nm / T m 1 NVT A total of nine independent properties sampled per NVT simulation with ms2 Any equilibrium thermodynamic property = combination of nm s
8 Cyclohexane Rigid 6CLJ united-atom model by Merker et al., Fluid Phase Equilib. 315 (2012) 77 Density / mol m -3
9 Cyclohexane potential energy Density / mol m -3 Present simulation data Span and Wagner, Int. J. Thermophys. 24 (2003) 41 Penoncello et. al., Int. J. Thermophys. 16 (1995) 519
10 Cyclohexane pressure Density / mol m -3 Present simulation data Span and Wagner, Int. J. Thermophys. 24 (2003) 41 Penoncello et. al., Int. J. Thermophys. 16 (1995) 519
11 Cyclohexane isochoric heat capacity Density / mol m -3 Present simulation data Span and Wagner, Int. J. Thermophys. 24 (2003) 41 Penoncello et. al., Int. J. Thermophys. 16 (1995) 519
12 Ongoing project with ms2 9 independent thermodynamic data types from one NVT simulation Generation of an extensive dataset in an automatized fashion Parallel execution of (~80 simulation runs) x (9 data points) Fit of a fundamental EOS in terms of F res to these data EOS may serve for the optimization of the force field EOS may be the starting point for technical EOS fitting Poster T21: Thol, Rutkai, Span, Vrabec: Molecular simulation of thermodynamic properties and an equation of state for the Lennard-Jones model fluid
13 Present status and outlook for ms2 An efficient molecular simulation tool for thermodynamic properties of homogeneous fluids New features Ewald summation for ionic systems Pair correlations functions MPI/OpenMP hybrid parallelization Next development steps Integer arithmetics Internal molecular degrees of freedom Execution on graphics processing units (GPUs)
14 Force fields for Hydrazine and two derivates 600 Hydrazine T / K Hydrazine Dimethylhydrazine Monomethylhydrazine Monomethylhydrazine Simulation, this work r / mol/l Experimental data from the literature Simulation, Gutowski et al Dimethylhydrazine
15 VLE and gas solubility of systems containing Hydrazine, Monomethylhydrazine and Dimethylhydrazine Elts et al., Fluid Phase Equilib (2012) Water + 1 bar H i / MPa H i / MPa Hydrazine Monomethylhydrazine Stickstoff Nitrogen N2 Argon Ar Stickstoff Nitrogen Argon T / K 380 Simulation, this work + Experimental data (literature) Peng-Robinson EOS x,y H2O / mol/mol H i / MPa Stickstoff Nitrogen Carbon Kohlenstoffmonoxid monoxide Argon Dimethylhydrazinee T / K
16 VLE and gas solubility of CO 2 -expanded liqiuds Cyclohexane (C 6 H 12 ) 6 LJ sites Carbon dioxide (CO 2 ) 3 LJ sites + point quadrupole Oxygen (O 2 ) 2 LJ sites + point quadrupole Cyclohexanol (C 6 H 10 O) 7 LJ sites + point dipole + CH - O Cyclohexanol (C 6 H 12 O) 7 LJ sites + 3 point charges Poster T37: Merker, Hsieh, Lin, Hasse, Vrabec: Fluid phase equilibria for the oxidation of cyclohexane in carbon dioxide expanded liquids from experiment, molecular simulation, Peng-Robinson EOS and COSMO-SAC + H
17 Speed of sound of Hexamethylsiloxane (MM) deviation [%] K EOS: Colonna et al., pressure [MPa] Poster T12: Dubberke, Vrabec: Speed of sound of siloxanes as workings fluids in Organic Rankine Cycles
18 Acetone in N 2 and O 2 under extreme conditions C -30 C 30 C 0 C 50 C 70 C Acetone Presentation Wed 10:40: Windmann, Köster, Vra.: Study on vapor-liquid equilibria of nitrogen + acetone and oxygen + acetone with a focus on the extended critical region p / MPa x N2 / mol/mol 90 C C
19 Massively parallel molecular dynamics code: ls1 Identical force field types as with ms2 Additionally, Tersoff potential for solids Large systems, long time scales Concurrency in space, not in time Flow Nucleation
20 Scaling tests on Cray XE6 (Hermit) Peak performance: Pflops operations / s Racks: 38 with 96 nodes each Nodes: 3552 Cores: (2 sockets each with 16 cores / node) Processor: AMD 2.3 GHz #12 (November 2011)
21 Rechenzeit / s Execution time / s Scaling of ls1 on Cray XE6 (Hermit) Computational effort ~ O(N 1 ) cores Molecule number Molekülzahl time steps
22
23 Octree load balancing strategy of ls1 Hierarchical subdivision of the simulation volume through recursive bisection Choice of planes for bisection that lead to an equal load in the subvolumes Designed for strongly inhomogeneous molecular systems Capable to deal with rapidly changing inhomogeneity
24 Scenarios for scaling tests of ls1 Bulk: homogeneous liquid (Ethylene oxide) Domain decomposition trivial Slab: liquid slab surrounded by vapor in equilibrium (Argon) Topology for domain decomposition simple Droplet: liquid droplet surrounded by vapor in equilibrium (Argon) Topology for domain decomposition more complex Bulk Slab Droplet
25 Scaling of ls1 on Cray XE6 (Hermit) Strong Scaling 2 22 = molecules Execution time / s 1000 Bulk Slab Film Tropfen Droplet Cores time steps
26 Scaling of ls1 on Cray XE6 (Hermit) Strong scaling 2 26 = molecules Execution time / s Film Bulk Tropfen Droplet Slab Bulk Cores time steps
27 VLE of Nitrogen + Ethane Stoll et al., AIChE J. 49 (2003) ξ = Simulation simulation Peng-Robinson equation of state Experiment experimental (lit.) data pressure [MPa] K 200 K mole fraction (N 2 )
28 Direct simulation of the LLE Mixture of 60 mol-% Nitrogen and 40 mol-% Ethane Identical molecular model as before 20,000 molecules Molecular dynamics with ls1 Canonical ensemble (NVT) Initial configuration with randomly dispersed components
29 LLE of N 2 + C 2 H 6 after 48 ns mole fraction 1,0 0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0, box length [nm] N2 C2H6 Temperature: 128 K Pressure: 11.0 MPa Average over 500,000 time steps
30 LLE temperature dependence of N 2 + C 2 H temperature [K] p = 1.8 to 4.0 MPa, depending on exp. data Simulation, this work Experimental data, literature mole fraction of N 2 [mol/mol]
31 LLE pressure dependence of N 2 + C 2 H 6 Poster T25: Eckelsbach, Vrabec: Prediction of liquid-liquid equilibria of nitrogen + ethane with a molecular model that was adjusted to vapor-liquid equilibria pressure [MPa] T = 127 to 129 K, depending on exp. data + Simulation, this work Experimental data, literature mole fraction of N 2 [mol/mol]
32 Summary The computational effort of classical force field methods scales linearly with the molecule number Classical force fields contain the thermodynamic properties adequately Molecular dynamics simulations of inhomogeneous fluids may efficiently use Petaflop machines The spectrum of possible applications is very wide
Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene and carbon tetrachloride
- 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,
More informationMolecular 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 information1 Introduction. Stefan Eckelsbach 1, Thorsten Windmann 1, Ekaterina Elts 1, and Jadran Vrabec 1
Simulation of liquid-liquid equilibria with molecular models optimized to vapor-liquid equilibria and model development for Hydrazine and two of its derivatives Stefan Eckelsbach 1, Thorsten Windmann 1,
More informationMD simulation of methane in nanochannels
MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall
More informationMolecular 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 informationSOLUBILITY OF CO 2 IN BRANCHED ALKANES IN ORDER TO EXTEND THE PPR78 MODEL TO SUCH SYSTEMS
SOLUBILITY OF CO IN BRANCHED ALKANES IN ORDER TO EXTEND THE PPR78 MODEL TO SUCH SYSTEMS Fabrice MUTELET, Stéphane VITU and Jean-Noël JAUBERT (*) Institut National Polytechnique de Lorraine, Ecole Nationale
More informationTHE DETERMINATION OF CRITICAL FLOW FACTORS FOR NATURAL GAS MIXTURES
A REPORT ON THE DETERMINATION OF CRITICAL FLOW FACTORS FOR NATURAL GAS MIXTURES Part 1: The Thermodynamic Properties of the Components of Natural Gas and of Natural Gas Mixtures in the Ideal-Gas State
More informationSurface property corrected modification of the classical nucleation theory
CCP5 Annual Meeting Surface property corrected modification of the classical nucleation theory Sheffield Hallam University, September 15, 2010 Martin Horsch, Hans Hasse, and Jadran Vrabec The critical
More informationScalable, 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 informationDependence of the surface tension on curvature
Dependence of the surface tension on curvature rigorously determined from the density profiles of nanodroplets Athens, st September M. T. Horsch,,, S. Eckelsbach, H. Hasse, G. Jackson, E. A. Müller, G.
More informationA First Course on Kinetics and Reaction Engineering Example 1.2
Example 1.2 Problem Purpose This example shows how to determine when it is permissible to choose a basis for your calculations. It also illustrates how to use reaction progress variables and the initial
More informationMolecular models for cyclic alkanes and ethyl acetate as well as surface tension data from molecular simulation
Molecular models for cyclic alkanes and ethyl acetate as well as surface tension data from molecular simulation Stefan Eckelsbach 1, Tatjana Janzen 1, Andreas Köster 1, Svetlana Miroshnichenko 1, Yonny
More informationComparison of the GERG-2008 and Peng-Robinson Equations of State for Natural Gas Mixtures
RESEARCH ARTICLE OPEN ACCESS Comparison of the and Peng-Robinson Equations of State for Natural Gas Mixtures L. F. Baladão*, R. P. Soares**, P. R. B. Fernandes*** * (Virtual Laboratory for Properties Prediction,
More informationPrediction of surface tension of binary mixtures with the parachor method
Prediction of surface tension of binary mixtures with the parachor method Tomáš Němec 1,a Institute of Thermomechanics ASCR, v.v.i., Dolejškova, 18 Praha 8, Czech Republic Abstract. The parachor method
More informationStructure-Property Relationships of Porous Materials for Carbon Dioxide Separation and Capture
Supporting Information Structure-Property Relationships of Porous Materials for Carbon Dioxide Separation and Capture Christopher E. Wilmer, 1 Omar K. Farha, 2 Youn-Sang Bae, 3,a Joseph T. Hupp, 2 and
More informationAn 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 informationIdeal Gas Behavior. NC State University
Chemistry 331 Lecture 6 Ideal Gas Behavior NC State University Macroscopic variables P, T Pressure is a force per unit area (P= F/A) The force arises from the change in momentum as particles hit an object
More informationCE 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 informationPhase 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 informationThermodynamic 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 informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
More informationSurface analysis algorithms in the mardyn program and the ls1 project
Surface analysis algorithms in the mardyn program and the ls1 project Stuttgart, 15 th December 1 M. T. Horsch Surface tension The virial route Bakker-Buff equation: γ R 2 out in dz z Normal pressure decays
More informationFluid Phase Behavior from Molecular Simulation: Hydrazine, Monomethylhydrazine, Dimethylhydrazine and Binary Mixtures Containing These Compounds
Fluid Phase Behavior from Molecular Simulation: Hydrazine, Monomethylhydrazine, Dimethylhydrazine and Binary Mixtures Containing These Compounds Ekaterina Elts 1, Thorsten Windmann 1, Daniel Staak 2, Jadran
More informationLarge-scale MD simulation of heterogeneous systems with ls1 mardyn
Large-scale MD simulation of heterogeneous systems with ls1 mardyn M. T. Horsch, R. Srivastava, S. J. Werth, C. Niethammer, C. W. Glass, W. Eckhardt, A. Heinecke, N. Tchipev, H.-J. Bungartz, S. Eckelsbach,
More informationUnderstanding 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 informationMolecular models for the hydrogen age: Hydrogen, nitrogen, oxygen, argon and water
Molecular models for the hydrogen age: Hydrogen, nitrogen, oxygen, argon and water Andreas Köster, Monika Thol, and Jadran Vrabec, Lehrstuhl für Thermodynamik und Energietechnik, Universität Paderborn,
More informationMutual 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 informationSurface tension of the two center Lennard-Jones plus quadrupole model fluid
Surface tension of the two center Lennard-Jones plus quadrupole model fluid Stephan Werth, Martin Horsch 1, Hans Hasse Laboratory of Engineering Thermodynamics, Department of Mechanical and Process Engineering,
More informationPhase 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 informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationIntermolecular Model Potentials and Virial Coefficients from Acoustic Data
JASEM ISSN 1119-8362 All rights reserved Full-text Available Online at https://www.ajol.info/index.php/jasem http://www.bioline.org.br/ja J. Appl. Sci. Environ. Manage. Vol.22 (2) 246-251. February 2018
More informationLiquid state isomorphism, Rosenfeld-Tarazona temperature scaling and Riemannian thermodynamic geometry
Liquid state isomorphism, Rosenfeld-Tarazona temperature scaling and Riemannian thermodynamic geometry Peter Mausbach Technical University of Cologne, 50678 Köln/Germany Andreas Köster and Jadran Vrabec
More informationIntroduction 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 informationPhase 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 informationADSORPTION IN MICROPOROUS MATERIALS: ANALYTICAL EQUATIONS FOR TYPE I ISOTHERMS AT HIGH PRESSURE
ADSORPTION IN MICROPOROUS MATERIALS: ANALYTICAL EQUATIONS FOR TYPE I ISOTHERMS AT HIGH PRESSURE A. L. MYERS Department of Chemical and Biomolecular Engineering University of Pennsylvania, Philadelphia
More informationThermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationThermodynamics at Small Scales. Signe Kjelstrup, Sondre Kvalvåg Schnell, Jean-Marc Simon, Thijs Vlugt, Dick Bedeaux
Thermodynamics at Small Scales Signe Kjelstrup, Sondre Kvalvåg Schnell, Jean-Marc Simon, Thijs Vlugt, Dick Bedeaux At small scales: 1. The systematic method of Hill 2. Small systems: surface energies 3.
More informationCOSMO-RS Theory. The Basics
Theory The Basics From µ to properties Property µ 1 µ 2 activity coefficient vapor pressure Infinite dilution Gas phase Pure compound Pure bulk compound Partition coefficient Phase 1 Phase 2 Liquid-liquid
More informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states
More informationImperfect Gases. NC State University
Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
More informationEngineering 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 informationChemical Potential. Combining the First and Second Laws for a closed system, Considering (extensive properties)
Chemical Potential Combining the First and Second Laws for a closed system, Considering (extensive properties) du = TdS pdv Hence For an open system, that is, one that can gain or lose mass, U will also
More informationMolecular 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 informationIntroduction. Monday, January 6, 14
Introduction 1 Introduction Why to use a simulation Some examples of questions we can address 2 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling Calculate
More informationChemistry 2000 Lecture 11: Chemical equilibrium
Chemistry 2000 Lecture 11: Chemical equilibrium Marc R. Roussel February 4, 2019 Marc R. Roussel Chemical equilibrium February 4, 2019 1 / 27 Equilibrium and free energy Thermodynamic criterion for equilibrium
More informationSchool of Chemical Engineering, Universidad del Valle, Ciudad Universitaria Melendez, Building 336, Apartado 25360, Cali, Colombia
Vapor-phase chemical equilibrium and combined chemical and vapor-liquid equilibrium for the ternary system ethylene + water + ethanol from reaction-ensemble and reactive Gibbs-ensemble molecular simulations
More informationa) ion-ion attractions b) London dispersion forces c) hydrogen bonding forces d) dipole-dipole attractions
Asgn #48: Intermolecular Forces Name Dec. 13, 2016 1. The intermolecular forces that are most significant in accounting for the high boiling point of liquid water relative to other substances of similar
More information591 TFLOPS Multi-TRILLION Particles Simulation on SuperMUC
International Supercomputing Conference 2013 591 TFLOPS Multi-TRILLION Particles Simulation on SuperMUC W. Eckhardt TUM, A. Heinecke TUM, R. Bader LRZ, M. Brehm LRZ, N. Hammer LRZ, H. Huber LRZ, H.-G.
More informationMolecular Simulation towards Efficient and Representative Subsurface Reservoirs Modeling. Dissertation by Ahmad Kadoura
Molecular Simulation towards Efficient and Representative Subsurface Reservoirs Modeling Dissertation by Ahmad Kadoura In Partial Fulfillment of the Requirements For the Degree of Doctor of Philosophy
More informationSOFTWARE INTELIGENT PACKAGE FOR PHASE EQULIBRIA (PHEQ) IN SYSTEMS APPLIED IN CHEMISTRY AND CHEMICAL ENGINEERING
SOFTWARE INTELIGENT PACKAGE FOR PHASE EQULIBRIA (PHEQ) IN SYSTEMS APPLIED IN CHEMISTRY AND CHEMICAL ENGINEERING Prof. Dr. Dan GEANĂ University Politechnica Bucharest Abstract The importance and role of
More informationAccurate thermodynamic-property models for CO 2 -rich mixtures
Available online at www.sciencedirect.com Energy Procedia 37 (2013 ) 2914 2922 GHGT-11 Accurate thermodynamic-property models for CO 2 -rich mixtures Roland Span*, Johannes Gernert, Andreas Jäger Thermodynamics,
More informationMolecular model for carbon dioxide optimized to vapor-liquid equilibria. Abstract
XXX Molecular model for carbon dioxide optimized to vapor-liquid equilibria Thorsten Merker and Cemal Engin Laboratory of Engineering Thermodynamics, University of Kaiserslautern, 67633 Kaiserslautern,
More informationMODELING OF PHASE EQUILIBRIA FOR BINARY AND TERNARY MIXTURES OF CARBON DIOXIDE, HYDROGEN AND METHANOL
MODELING OF PHASE EQUILIBRIA FOR BINARY AND TERNARY MIXTURES OF CARBON DIOXIDE, HYDROGEN AND METHANOL Neil R. Foster *, Keivan Bezanehtak, Fariba Dehghani School of Chemical Engineering and Industrial
More informationThermophysical properties of dry and humid air. by molecular simulation including dew point. calculations with the Mollier ensemble
Thermophysical properties of dry and humid air by molecular simulation including dew point calculations with the Mollier ensemble Bernhard Eckl 1, Thorsten Schnabel 1, Jadran Vrabec 2, Martin Wendland
More informationEQUATION OF STATE DEVELOPMENT
EQUATION OF STATE DEVELOPMENT I. Nieuwoudt* & M du Rand Institute for Thermal Separation Technology, Department of Chemical Engineering, University of Stellenbosch, Private bag X1, Matieland, 760, South
More informationBinary Hard-Sphere Mixtures Within Spherical Pores
Journal of the Korean Physical Society, Vol. 35, No. 4, October 1999, pp. 350 354 Binary Hard-Sphere Mixtures Within Spherical Pores Soon-Chul Kim Department of Physics, Andong National University, Andong
More informationMOLECULAR DYNAMICS SIMULATION OF VAPOR BUBBLE NUCLEATION ON A SOLID SURFACE. Tatsuto Kimura and Shigeo Maruyama
MOLECULAR DYNAMICS SIMULATION OF VAPOR BUBBLE NUCLEATION ON A SOLID SURFACE Tatsuto Kimura and Shigeo Maruyama * Department of Mechanical Engineering, The University of Tokyo, 7-- Hongo, Bunkyo-ku, Tokyo
More informationRigorous calculation of LNG flow reliefs using the GERG-2004 equation of state
Rigorous calculation of LNG reliefs using the GERG-2004 equation of state Luigi Raimondi Process Simulation Services www.xpsimworld.com Via Galvani 105, 20025 Legnano (MI) - Italy The design of process
More informationSystematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids
Systematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids Cameron F. Abrams Department of Chemical and Biological Engineering Drexel University Philadelphia, PA USA 9 June
More informationADSORPTION AND INTERFACIAL PROPERTIES OF FLUIDS FROM MOLECULAR SIMULATION
ADSORPTION AND INTERFACIAL PROPERTIES OF FLUIDS FROM MOLECULAR SIMULATION by Xiongce Zhao B.S. in Chemical Engineering, Tsinghua University, 1993 M.S. in Chemical Engineering, Tsinghua University, 1998
More informationarxiv: v1 [physics.comp-ph] 25 Jul 2015
To appear in Molecular Physics Vol. 11, No. 17, 014, 1 13 arxiv:07.07130v1 [physics.comp-ph] 5 Jul 0 Long range correction for multi-site Lennard-Jones models and planar interfaces Stephan Werth 1, Gabor
More informationReview Article. VLE Properties from ISM Equation of State: Application to Pure and Mixture
Review Article PHYSICAL CHEMISTRY RESEARCH Published by the Iranian Chemical Society www.physchemres.org info@physchemres.org Phys. Chem. Res., Vol., No. 1, 16-, March 015. DOI: 10.06/pcr.015.798 VLE Properties
More informationthe expansion for the Helmholtz energy derived in Appendix A, part 2, the expression for the surface tension becomes: σ = ( a + ½ k(ρ) ρ 2 x ) dx
Motivation The surface tension plays a major role in interfacial phenomena. It is the fundamental quantity that determines the pressure change across a surface due to curvature. This in turn is the basis
More informationMOLECULAR DYNAMICS STUDY OF THE NUCLEATION OF BUBBLE
CAV2:sessionA.5 MOLECULAR DYNAMICS STUDY OF THE NUCLEATION OF BUBBLE Takashi Tokumasu, Kenjiro Kamijo, Mamoru Oike and Yoichiro Matsumoto 2 Tohoku University, Sendai, Miyagi, 98-8577, Japan 2 The University
More informationSteady-State Molecular Diffusion
Steady-State Molecular Diffusion This part is an application to the general differential equation of mass transfer. The objective is to solve the differential equation of mass transfer under steady state
More informationCHAPTER 1. MEASURES OF ATMOSPHERIC COMPOSITION
1 CHAPTER 1. MEASURES OF ATMOSPHERIC COMPOSITION The objective of atmospheric chemistry is to understand the factors that control the concentrations of chemical species in the atmosphere. In this book
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationIMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER
IMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER ABSTRACT Hong Peng 1, Anh V Nguyen 2 and Greg R Birkett* School of Chemical Engineering, The University of Queensland Brisbane,
More informationA MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE
A MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE Shigeo Maruyama and Tatsuto Kimura Department of Mechanical Engineering The University of Tokyo 7-- Hongo, Bunkyo-ku, Tokyo -866,
More informationOFB Chapter 6 Condensed Phases and Phase Transitions
OFB Chapter 6 Condensed Phases and Phase Transitions 6-1 Intermolecular Forces: Why Condensed Phases Exist 6- The Kinetic Theory of Liquids and Solids 6-3 Phase Equilibrium 6-4 Phase Transitions 6-5 Phase
More informationCE 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 informationEEC 503 Spring 2009 REVIEW 1
EEC 503 Spring 2009 REVIEW 1 1. Why are chemical reactions important to energy, environmental and process engineering? Name as many reasons as you can think of. 2. What is a chemical reaction? 3. What
More informationevidyarthi.in Thermodynamics Q 1.
SUBJECTIVE PROBLEMS: Q 1. Thermodynamics The enthalpy for the following reaction ( H o ) at 25 o C are given below: (i) 1/2 H 2 (g) + 1/2 O 2 (g) OH(g) 10.06 kcal (ii) H 2 (g) 2H(g) 104.18 kcal (iii) O
More informationIndex to Tables in SI Units
Index to Tables in SI Units Table A-1 Atomic or Molecular Weights and Critical Properties of Selected Elements and Compounds 926 Table A-2 Properties of Saturated Water (Liquid Vapor): Temperature Table
More informationMolecular Dynamics Simulation on Permeation of Acetone/Nitrogen Mixed Gas
Molecular Dynamics Simulation on Permeation of Acetone/Nitrogen Mixed Gas through Al 2 O 3 Microporous Membranes Kailiang Yin 1,2, Duanjun Xu 2 and Jing Zhong 1 1. Department of Chemical Engineering, Jiangsu
More informationMolecular modeling and simulation in chemistry
Molecular modeling and simulation in chemistry 1/20 Elements of modeling: Nuclei + electrons + photons: QED Nuclei + electrons: Schrödinger equation Auxiliary interaction centers (sites) Atoms classical
More informationDescription of HFO-1234yf with BACKONE equation of state
Description of HFO-1234yf with BACKONE equation of state Ngoc Anh Lai a,*, Jadran Vrabec b, Gabriele Raabe c, Johann Fischer d, Martin Wendland d a Heat Engineering Department, Institute of Heat Engineering
More informationFinal Exam Chem 260 (12/20/99) Name (printed) Signature
Final Exam Chem 260 (12/20/99) Name (printed) Signature 1. For the reaction, (14 pts total) C(graph) + CO 2 (g) 2 CO(g) (a) Write the equilibrium constant expression: (2 pts) K = P CO 2 P CO2 (b) Using
More informationTable of Contents Preface PART I. MATHEMATICS
Table of Contents Preface... v List of Recipes (Algorithms and Heuristics)... xi PART I. MATHEMATICS... 1 Chapter 1. The Game... 3 1.1 Systematic Problem Solving... 3 1.2 Artificial Intelligence and the
More informationFUNDAMENTALS of Thermodynamics
SOLUTION MANUAL SI UNIT PROBLEMS CHAPTER 15 SONNTAG BORGNAKKE VAN WYLEN FUNDAMENTALS of Thermodynamics Sixth Edition CONTENT SUBSECTION PROB NO. Correspondence table Concept-Study Guide Problems 1-20 Equilibrium
More informationThe Monte Carlo Method in Condensed Matter Physics
The Monte Carlo Method in Condensed Matter Physics Edited by K. Binder With Contributions by A. Baumgärtner K. Binder A.N. Burkitt D.M. Ceperley H. De Raedt A.M. Ferrenberg D.W. Heermann H.J. Herrmann
More informationPORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES
PORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES Atichat Wongkoblap*, Worapot Intomya, Warangkhana Somrup, Sorod Charoensuk, Supunnee Junpirom and Chaiyot Tangsathitkulchai School of Chemical
More informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationNon-equilibrium molecular dynamics simulation study of the behavior of hydrocarbon-isomers in silicalite
Fluid Phase Equilibria 194 197 (2002) 309 317 Non-equilibrium molecular dynamics simulation study of the behavior of hydrocarbon-isomers in silicalite S. Furukawa a,b,, C. McCabe a,c, T. Nitta b, P.T.
More informationOpening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations
Supporting Information for: Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations D. Fairen-Jimenez *, S. A. Moggach, M. T. Wharmby, P. A. Wright, S. Parsons and T. Düren
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationMolecular simulation of the surface tension of 33 multi-site models for real fluids
Molecular simulation of the surface tension of 33 multi-site models for real fluids Stephan Werth a, Martin Horsch a,, Hans Hasse a a University of Kaiserslautern, Laboratory of Engineering Thermodynamics,
More informationSurface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water
Surface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water Gabriel V. Lau, Ian J. Ford, Patricia A. Hunt, Erich A. Müller, and George Jackson Citation: The Journal of
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationA MOLECULAR DYNAMICS SIMULATION OF WATER DROPLET IN CONTACT WITH A PLATINUM SURFACE
The 6th ASME-JSME Thermal Engineering Joint Conference March 16-2, 23 TED-AJ3-183 A MOLECULAR DYNAMICS SIMULATION OF WATER DROPLET IN CONTACT WITH A PLATINUM SURFACE Tatsuto KIMURA Department of Mechanical
More informationDevelopment 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 informationarxiv: v1 [physics.comp-ph] 21 Apr 2009
arxiv:0904.3191v1 [physics.comp-ph] 21 Apr 2009 Molecular Modeling and Simulation of Thermophysical Properties: Application to Pure Substances and Mixtures Bernhard Eckl, Martin Horsch, Jadran Vrabec,
More informationSupplementary material to Simultaneous description of bulk and interfacial properties of fluids by the Mie potential
To appear in Molecular Physics Vol. 00, No. 00, Month 2016, 1 11 Supplementary material to Simultaneous description of bul and interfacial properties of fluids by the Mie potential Stephan Werth 1, Katrin
More informationEquations of State. Equations of State (EoS)
Equations of State (EoS) Equations of State From molecular considerations, identify which intermolecular interactions are significant (including estimating relative strengths of dipole moments, polarizability,
More information2012 AP CHEMISTRY FREE-RESPONSE QUESTIONS
01 AP CHEMISTRY FREE-RESPONSE QUESTIONS. A sample of a pure, gaseous hydrocarbon is introduced into a previously evacuated rigid 1.00 L vessel. The pressure of the gas is 0.00 atm at a temperature of 17C.
More informationModelling of Adsorption and Diffusion in Dual-Porosity Materials: Applications to Shale Gas
Modelling of Adsorption and Diffusion in Dual-Porosity Materials: Applications to Shale Gas Martin Lísal Institute of Chemical Process Fundamentals, CAS, Prague, Czech Republic Faculty of Science, J. E.
More information