Long-Term Atomistic Simulation of Heat and Mass Transport
|
|
- Dwight Thompson
- 5 years ago
- Views:
Transcription
1 Long-Term Atomistic Simulation of Heat and Mass Transport Kevin G. Wang 1, Mauricio Ponga 2, Michael Ortiz 2 1 Department of Aerospace and Ocean Engineering, Virginia Polytechnic Institute and State University 2 Division of Engineering and Applied Sciences, California Institute of Technology
2 Motivation Slow deformation-diffusion coupled processes hydrogen embrittlement in fuel cells initial after 10 cycles Li et al, JACS (2007) Mg nanowires collapse after 10 cycles of charging/discharging; temperature: 300K ~ 600K; time scale: seconds to hours!
3 Overview Long-term atomistic simulation of heat and mass transport - the essential difficulty: multiple disparate time scales atomic level rate-limiting processes: thermal vibrations, individual hops of atoms macroscopic processes of interest: hydrogen embrittlement, ductile fracture, fatigue, etc. time-scale gap: from molecular dynamics (femtoseconds) to to macroscopic (seconds, minutes, )! - modeling approach fine scale: non-equilibrium statistical thermodynamics coarse scale: discrete kinetic laws - previous work: Kulkarni et al ( HotQC ), JMPS (2008), Venturini et al, JMCE (2012), Ariza et al, IJF (2012)
4 Non-Equilibrium Thermodynamics Statistical thermodynamics for multi-species crystals - N particles, M species occupancy function 1, particle i is from species k for i = 1,, N, k = 1,, M, n ik = 0, oterwise - assumption: separation of time scales fine scale: thermal vibration (~fs) modeled coarse scale: global relaxation (>>fs) solved for probability: q, p, *n+ ~ p q, p, *n+, - maximum entropy principle: max p with S p = k B log p (information entropy) 1 f = 3N f q, p, *n+ p q, p, *n+ dqdp n O Γ NM subject to: H i = E e i i = 1,, N (atomic internal energy) nn ik = x ik = N k i = a 1,, N (atomic molar fraction) i<1 n ik
5
6 Non-Equilibrium Thermodynamics Meanfield approximation (Yeomans 1992) - given a class of trial Hamiltonian parameterized by π : 0i ( p, q, *n+; π ), the trial probability density function can be defined as p 0 q, p, n; π ; β, *γ+ = 1 Ξ 0 exp β T 0 + γ T n - for fixed Lagrange multipliers β, *γ+, the optimal probability of the trial space *p 0 + is obtained by maximizing the constrained entropy with respect to π, i.e. p 0 = argmax *π+ L, p 0 π, β, *γ+ - - theorem (bounding principle for the free entropy): Φ *β+, *γ+ max Φ( 0, β, *γ+) h 0 the maximum free entropy of the trial space provides a lower bound of the actual free entropy p 0 converges to p as 0
7 Non-Equilibrium Thermodynamics Meanfield model: an example - in practice there is always a trade-off between accuracy and computability - we define the trial Hamiltionian as 0i q, p ; q, p, m, ω = 1 p 2m i p 2 i + m 2 iω i q i 2 i q 2 i *π+ uncoupled harmonic oscillators N p 0 ( q, p, *n+) = 1 1 exp β Ξ i p 0 2m i p 2 i + m 2 iω i q i 2 i q 2 i i<1 Gaussian distribution p 0 = argmax L,p 0, β, *γ+ - p i, q i, ω i,*m i + meanfield parameters: q i (mean atomic position), p i (mean atomic momentum), ω i (frequency of oscillation), m i (mean atomic mass) i = p p 0 2 2m + V i(*q+) explicit formulations of S, Φ, F,...
8 Modal Validation Thermal expansion coefficient: Mg single crystals - at fixed, uniform temperature, maximizes grand-canonical free entropy (Φ) with respect to atomic frequency (ω i ) and mean atomic position (q i ) q = argmin q,*ω+ Φ( q, ω, T) α V T = 1 V 0 V T 1 V 0 V 1 V 0 T 1 T 0 (volumetric thermal expansion coefficient) - comparison with experiments temperature simulation experimental rel. error (K) ( 10 ;6 ) ( 10 ;6 data source ) (%) Goens & Schmid (1936) Raynor & Hume-Rothery (1939) Goens & Schmid (1936) Hanawalt & Frevel (1938) Raynor & Hume-Rothery (1939) Raynor & Hume-Rothery (1939) 14.53
9 Modal Validation Nickel-Palladium lattice parameter - problem setup: FCC random solution of nickel (Ni) and palladium (Pd) 365 lattice sites; periodic boundary condition temperature and molar fraction are prescribed x Ni varies from 0 (pure Pd) to 1 (pure Ni) equilibrium lattice parameter a obtained by relaxing *q+ and *ω+ to maximize grand canonical free entropy uniform molar fraction (x Ni, x Pd ), uniform temperature T
10 Detailed balance of energy - balance of energy at a particle u i = w i + μ T i x i + q i Discrete Kinetic Laws (w i = v h i A A α α α<1 is the external power) requires kinetic models for x i and q i - assume x i and q i can be divided into particle-to-particle fluxes of the form with - assume linear kinetics Q ij = κ ij (β i β j ) J ij = D ij (γ i γ j ) κ ij : pairwise conductivity coefficient D ij : pairwise diffusivity coefficient - question: how to determine κ ij and D ij?
11 A simple example Discrete Kinetic Laws - consider only pairwise flux for only the nearest neighbors (based on a perfect crystal lattice) J ij = D ij γ i γ j, if i and j are nearest neigbors 0, oterwise - require consistency with Fick s law in continuum mechanics (in a finite-difference sense) Fick s second law: x = D 0 (x γ) discrete diffusion model: x i = D ij (γ i γ j ) j Nei(i) assume x, γ C 2 (R 3 ), x i = D 0 x i γ i + O(a 2 ) iff D ij = D 0 2a 2 (x i + x j ) (for fcc and bcc lattices) - can be extended to more than one layers of neighbors - the same procedure can be applied to heat conduction
12 Summary Thermo-chemo-mechanical coupled analysis given m n;1, q n;1, p n;1, ω n;1, β n, *x n + [1] assume trial Hamiltonian H 0 q, p ; π with π i = (m i, q i, p i, ω i ) [2] obtain optimal π n by maximizing free entropy: max m i, p i, q i, *ω i + L,ρ 0, β, γ - [5] solve kinetic equations and obtain β n:1, *x n:1 + [3] extract quantities of interest e.g. S n, Φ n, E n, [4]
13 Model Validation Heat conduction in Si nanowires - experiment (D. Li et al. 2003) specimen: single-crystal Si (111) nanowires R = 11, 18.5, 28, 57.5 nm measurement: thermal conductivity main discovery: significant size effect HRTEM image of a 22 nm SiNW
14 Model Validation Computational model and simulations - consider a slice of the SiNWs - create amorphous layer by random perturbation - pair-wise thermal conductivity: thermal conductivity amorphous layer atomic temperature
15 Model Validation Hydrogen diffusion in Pd nanofilms - experiment (Y. Li and Y.-T. Cheng, 1996) specimen: Pd (111) thin films with thickness L = 22 nm, 46 nm, 135 nm method: electrochemical stripping step 1: introduce H to Pd by cathodic polarization (set ΔV = 0.83V) Pd + H 2 O + e ; Pd H ads + OH ; step 2: remove H from Pd-H by anodic stripping (switch ΔV to 0.5V) Pd H ads + OH ; Pd + H 2 O + e ; l electrolyte solution substrate: Au-coated Ni Pd film measurement: time history of discharge current at the surface of Pd film C l, t I t = FS J H (l, t) = FSD H = FSD HZ x l, t l l<0 N A V cell l l<0
16 Models and simulations Model Validation - simplification of thermodynamic model 1D, constant temperature, fixed Pd lattice linearization of Pd-H interatomic potential equilibrium condition γ i = log x i 1 x i + - discrete Fick s law x i = D ij (γ j γ i ) j Nei(i) C ij k B T x j + B i k B T j I H,j i x i - highlights fully atomistic t max = 1.0 sec, dt ~ 10-6 sec L = 46 nm L = 135 nm
17 26nm 5nm Model Verification Nanovoid growth under tension (coupled with QC) - problem setup objective: study the evolution of nanovoid in a single FCC crystal of Cu representative volume: 26nm x 26nm x 26nm (~ 10 6 atoms) initial void diameter: 2nm initial atomic zone: 5nm x 5nm x 5nm (12a 0 x 12a 0 x 12a 0 ) interatomic potential: EAM-Mishin (Y. Mishin et al. 2001) uniaxial and triaxial loading with prescribed boundary displacement - strain rate (1/s): 10 5, 10 7, 10 8, 10 9, atomic zone void QC approximation cut-view of initial FE mesh
18 Nanovoid Growth Uniaxial loading, T = 300K, ε = /s 1/6,112-1/6,112-1/6,112-1/6,112- ε = 6.6% (emission of dislocation) ε = 6.8% (leading Shockley partials) atomic temperature ε = 7.0% (trailing Shockley partials) *111+ dislocation core identified by CSD* * centro-symmetric deviation: CSD = N 2 i<1 q i + q i: N 2 a 0 2 before dislocation emission after dislocation emission
19 virial stress (GPa) Nanovoid Growth Tri-axial loading, T = 300K, ε = 10 5, 10 8, 10 9, /s stress-strain response 10 dislocation core colored by CSD (ε = /s) - virial stress (Allen and Tildesley, 1989) ς ij = 1 V Ω m k q i k q i k q j k q j k q i l q i k f j kl k Ω ε = 10 5, ε = 10 8, ε = 10 9, ε = ε, = HotQC, 10 10, dynamic ε = ε, = MD 10 10, MD ref engineering strain
20 Summary and Future Work Long-term atomistic simulation of heat and mass transport - non-equilibrium statistical thermodynamics model maximum entropy principle meanfield approximation - discrete kinetic laws discrete Fick s law discrete Fourier s law - coupling: a partitioned procedure - validations and applications Future work - computational methods for thermo-chemo-mechanical coupling - design and validation of discrete kinetic laws - validations and applications Reference - Venturini, Wang, Romero, Ariza, and Ortiz, JMPS (submitted) Acknowledgement: Army Research Lab (ARL)
of Long-Term Transport
Atomistic Modeling and Simulation of Long-Term Transport Phenomena in Nanomaterials K.G. Wang, M. and M. Ortiz Virginia Polytechnic Institute and State University Universidad de Sevilla California Institute
More informationLong-Term Atomistic Simulation of Hydrogen Diffusion in Metals
Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals K.G. Wang 1, M. and M. Ortiz 2 Universidad de Sevilla 1 Virginia Polytechnic Institute and State University 2 California Institute of Technology
More informationLong-Term Atomistic Simulation of Hydrogen Diffusion in Metals
Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals Kevin Wang (1), Pilar Ariza (2), and Michael Ortiz (3) (1) Department of Aerospace and Ocean Engineering Virginia Polytechnic Institute and
More informationIAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.
IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 8: Lattice Waves in 1D Monatomic Crystals Outline Overview of Lattice Vibrations so far Models for Vibrations in Discrete 1-D Lattice Example of Nearest
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationDiffusion in multicomponent solids. Anton Van der Ven Department of Materials Science and Engineering University of Michigan Ann Arbor, MI
Diffusion in multicomponent solids nton Van der Ven Department of Materials Science and Engineering University of Michigan nn rbor, MI Coarse graining time Diffusion in a crystal Two levels of time coarse
More informationSample geometry and the brittle-ductile behavior of edge cracks in 3D atomistic simulations by molecular dynamics
Sample geometry and the brittle-ductile behavior of edge cracks in 3D atomistic simulations by molecular dynamics V. Pelikán, P. Hora, A. Machová, R. Kolman, A. Uhnáková Institute of Thermomechanics AS
More informationAn Atomistic-based Cohesive Zone Model for Quasi-continua
An Atomistic-based Cohesive Zone Model for Quasi-continua By Xiaowei Zeng and Shaofan Li Department of Civil and Environmental Engineering, University of California, Berkeley, CA94720, USA Extended Abstract
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationKinetic Monte Carlo: from transition probabilities to transition rates
Kinetic Monte Carlo: from transition probabilities to transition rates With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More informationDIFFUSION IN SOLIDS. IE-114 Materials Science and General Chemistry Lecture-5
DIFFUSION IN SOLIDS IE-114 Materials Science and General Chemistry Lecture-5 Diffusion The mechanism by which matter is transported through matter. It is related to internal atomic movement. Atomic movement;
More informationWhy Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important
Why Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important The Silicon Age is Ending New materials and growth required Useful in design of Devices
More informationDensity Functional Modeling of Nanocrystalline Materials
Density Functional Modeling of Nanocrystalline Materials A new approach for modeling atomic scale properties in materials Peter Stefanovic Supervisor: Nikolas Provatas 70 / Part 1-7 February 007 Density
More informationThermodynamics of Damage: A Reliability Engineering Perspective
Thermodynamics of Damage: A Reliability Engineering Perspective Mohammad Modarres Presented at Department of Materials Science and Engineering Seminar UCLA, Los Angeles CA February 17, 217 Department of
More informationfiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2
More informationEntropy as an Indication of Damage in Engineering Materials
Entropy as an Indication of Damage in Engineering Materials Mohammad Modarres Presented at Entropy 2018: From Physics to Information Sciences and Geometry Barcelona, Spain, May 14-16, 2018 16 May 2018
More informationElastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals
PHYSICAL REVIEW B 70, 214113 (2004) Elastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals Sang Kyu Kwak and David A. Kofke Department of Chemical and Biological
More informationMorphological evolution of single-crystal ultrathin solid films
Western Kentucky University From the SelectedWorks of Mikhail Khenner March 29, 2010 Morphological evolution of single-crystal ultrathin solid films Mikhail Khenner, Western Kentucky University Available
More informationCalculation of thermal expansion coefficient of Fe3Al with the addition of transition metal elements Abstract Introduction Methodology
Calculation of thermal expansion coefficient of Fe Al with the addition of transition metal elements Tatiana Seletskaia, Leonid Muratov, Bernard Cooper, West Virginia University, Physics Departement, Hodges
More informationConcepts in Surface Physics
M.-C. Desjonqueres D. Spanjaard Concepts in Surface Physics Second Edition With 257 Figures Springer 1. Introduction................................. 1 2. Thermodynamical and Statistical Properties of
More informationEntropic methods to study the evolution of damage and degradation of materials
Entropic methods to study the evolution of damage and degradation of materials Mohammad Modarres Presented at 14th International Conference on Fracture Rhodes, Greece, June 18-23, 2017 20 JUNE 2017 Department
More informationReactive potentials and applications
1.021, 3.021, 10.333, 22.00 Introduction to Modeling and Simulation Spring 2011 Part I Continuum and particle methods Reactive potentials and applications Lecture 8 Markus J. Buehler Laboratory for Atomistic
More informationPhonons I - Crystal Vibrations (Kittel Ch. 4)
Phonons I - Crystal Vibrations (Kittel Ch. 4) Displacements of Atoms Positions of atoms in their perfect lattice positions are given by: R 0 (n 1, n 2, n 3 ) = n 10 x + n 20 y + n 30 z For simplicity here
More informationJournal of the Mechanics and Physics of Solids
Accepted Manuscript Diffusive molecular dynamics simulations of lithiation of silicon nanopillars J.P. Méndez, M. Ponga, M. Ortiz PII: S0022-5096(17)30880-3 DOI: 10.1016/j.jmps.2018.03.008 Reference: MPS
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationEE C245 ME C218 Introduction to MEMS Design Fall 2007
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 13: Material
More informationSteady-state diffusion is diffusion in which the concentration of the diffusing atoms at
Chapter 7 What is steady state diffusion? Steady-state diffusion is diffusion in which the concentration of the diffusing atoms at any point, x, and hence the concentration gradient at x, in the solid,
More information3.091 Introduction to Solid State Chemistry. Lecture Notes No. 5a ELASTIC BEHAVIOR OF SOLIDS
3.091 Introduction to Solid State Chemistry Lecture Notes No. 5a ELASTIC BEHAVIOR OF SOLIDS 1. INTRODUCTION Crystals are held together by interatomic or intermolecular bonds. The bonds can be covalent,
More informationPredicting the Dislocation Nucleation Rate as a Function of Temperature and Stress. Abstract
Journal of Materials Research, Invited Feature Article, in press (2011) Predicting the Dislocation Nucleation Rate as a Function of Temperature and Stress Seunghwa Ryu 1, Keonwook Kang 2, and Wei Cai 3
More informationPALLADIUM HYDROGEN System
The PALLADIUM HYDROGEN System F. A. LEWIS Department of Chemistry The Queen's University of Belfast Northern Ireland ACADEMIC PBESS 1967 LONDON NEW YORK CONTENTS Preface.......... v Chapter 1. General
More informationKinetics. Rate of change in response to thermodynamic forces
Kinetics Rate of change in response to thermodynamic forces Deviation from local equilibrium continuous change T heat flow temperature changes µ atom flow composition changes Deviation from global equilibrium
More informationModeling of slow time-scale behavior of fast molecular dynamic systems
Modeling of slow time-scale behavior of fast molecular dynamic systems ikun Tan, Amit Acharya, Kaushik Dayal Civil and Environmental Engineering, Carnegie Mellon University Abstract. This paper aims to
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More information3.021J / 1.021J / J / J / 22.00J Introduction to Modeling and Simulation Markus Buehler, Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 3.021J / 1.021J / 10.333J / 18.361J / 22.00J Introduction to Modeling and Simulation Markus Buehler, Spring 2008 For information about citing these materials or our
More informationT. Egami. Model System of Dense Random Packing (DRP)
Introduction to Metallic Glasses: How they are different/similar to other glasses T. Egami Model System of Dense Random Packing (DRP) Hard Sphere vs. Soft Sphere Glass transition Universal behavior History:
More informationIntroduction to atomic scale simulations
Powder Technology course Autumn semester 2017 Introduction to atomic scale simulations Peter M Derlet Condensed Matter Theory Paul Scherrer Institut peter.derlet@psi.ch Outline of lecture The big picture
More informationQuasi-Harmonic Theory of Thermal Expansion
Chapter 5 Quasi-Harmonic Theory of Thermal Expansion 5.1 Introduction The quasi-harmonic approximation is a computationally efficient method for evaluating thermal properties of materials. Planes and Manosa
More informationPREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures
PREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures Kevin G. Wang Virginia Tech Patrick Lea, Alex Main, Charbel Farhat Stanford University
More informationNUMERICAL SIMULATION OF CONCRETE EXPOSED TO HIGH TEMPERATURE DAMAGE AND EXPLOSIVE SPALLING
NUMERICAL SIMULATION OF CONCRETE EXPOSED TO HIGH TEMPERATURE DAMAGE AND EXPLOSIVE SPALLING Prof. Joško Ožbolt 1 Josipa Bošnjak 1, Goran Periškić 1, Akanshu Sharma 2 1 Institute of Construction Materials,
More informationMultiscale modeling and simulation of shock wave propagation
University of Iowa Iowa Research Online Theses and Dissertations Spring 2013 Multiscale modeling and simulation of shock wave propagation Cory Nelsen University of Iowa Copyright 2013 Cory Nelsen This
More informationHow to model chemical interactions II
1.021, 3.021, 10.333, 22.00 Introduction to Modeling and Simulation Spring 2011 Part I Continuum and particle methods How to model chemical interactions II Lecture 6 Markus J. Buehler Laboratory for Atomistic
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationApplications to simulations: Monte-Carlo
Applications to simulations: Monte-Carlo A.C. Maggs ESPCI, Paris June 2016 Summary Algorithms Faster/simpler codes Thermodynamics of Electric fields Partition function of electric field Thermal Casimir/Lifshitz
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationFrom Atoms to Materials: Predictive Theory and Simulations
From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology
More informationComputational Materials Science
Computational Materials Science 50 (2011) 1731 1743 Contents lists available at ScienceDirect Computational Materials Science journal homepage: www. elsevi er. com/ locate/ c ommatsci Stability and size-dependency
More informationStructure of Surfaces
Structure of Surfaces C Stepped surface Interference of two waves Bragg s law Path difference = AB+BC =2dsin ( =glancing angle) If, n =2dsin, constructive interference Ex) in a cubic lattice of unit cell
More informationUniversity of Ljubljana. Faculty of Mathematics and Physics. Department of Physics. High-entropy alloys. Author: Darja Gačnik
University of Ljubljana Faculty of Mathematics and Physics Department of Physics Seminar I 1 st year of Master study programme High-entropy alloys Author: Darja Gačnik Mentor: prof. Janez Dolinšek Ljubljana,
More informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More informationModule 16. Diffusion in solids II. Lecture 16. Diffusion in solids II
Module 16 Diffusion in solids II Lecture 16 Diffusion in solids II 1 NPTEL Phase II : IIT Kharagpur : Prof. R. N. Ghosh, Dept of Metallurgical and Materials Engineering Keywords: Micro mechanisms of diffusion,
More informationCoarse-graining molecular dynamics in crystalline solids
in crystalline solids Xiantao Li Department of Mathematics Penn State University. July 19, 2009. General framework Outline 1. Motivation and Problem setup 2. General methodology 3. The derivation of the
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationAtoms, electrons and Solids
Atoms, electrons and Solids Shell model of an atom negative electron orbiting a positive nucleus QM tells that to minimize total energy the electrons fill up shells. Each orbit in a shell has a specific
More informationSurface Physics Surface Diffusion. Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science
Surface Physics 008 8. Surface Diffusion Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science Random-Walk Motion Thermal motion of an adatom on an ideal crystal surface: - Thermal excitation the adatom
More informationUsing the Anharmonic Correclated Einstein Model for Calculation the Thermodynamic Parameters and Cumulants of Dopant Face Cubic Center Crystals
AASCIT Journal of Physics 015; 1(1: 1-5 Published online March 0, 015 (http://www.aascit.org/journal/physics Using the Anharmonic Correclated instein Model for Calculation the Thermodynamic Parameters
More informationI. Collective Behavior, From Particles to Fields
I. Collective Behavior, From Particles to Fields I.A Introduction The object of the first part of this course was to introduce the principles of statistical mechanics which provide a bridge between the
More informationCombined atomistic-continuum modeling of short-pulse laser melting and disintegration of metal films
Combined atomistic-continuum modeling of short-pulse laser melting and disintegration of metal films Dmitriy S. Ivanov and Leonid V. Zhigilei* Department of Materials Science and Engineering, University
More informationLattice Strain Due to an Atomic Vacancy
Int. J. Mol. Sci. 2009, 10, 2798-2808; doi:10.3390/ijms10062798 OPEN ACCESS Article International Journal of Molecular Sciences ISSN 1422-0067 www.mdpi.com/journal/ijms Lattice Strain Due to an Atomic
More informationLectures on: Introduction to and fundamentals of discrete dislocations and dislocation dynamics. Theoretical concepts and computational methods
Lectures on: Introduction to and fundamentals of discrete dislocations and dislocation dynamics. Theoretical concepts and computational methods Hussein M. Zbib School of Mechanical and Materials Engineering
More informationChapter 2. Theoretical background. 2.1 Itinerant ferromagnets and antiferromagnets
Chapter 2 Theoretical background The first part of this chapter gives an overview of the main static magnetic behavior of itinerant ferromagnetic and antiferromagnetic materials. The formation of the magnetic
More informationApplication to modeling brittle materials
1.01, 3.01, 10.333,.00 Introduction to Modeling and Simulation Spring 011 Part I Continuum and particle methods Application to modeling brittle materials Lecture 7 Markus J. Buehler Laboratory for Atomistic
More informationSimulation of Thermomechanical Couplings of Viscoelastic Materials
Simulation of Thermomechanical Couplings of Viscoelastic Materials Frank Neff 1, Thomas Miquel 2, Michael Johlitz 1, Alexander Lion 1 1 Institute of Mechanics Faculty for Aerospace Engineering Universität
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
More informationBridging to the Continuum Scale for Ferroelectric Applications
Bridging to the Continuum Scale for Ferroelectric Applications Shanfu Zheng and Alberto Cuitiño Mechanical and Aerospace Engineering, Rutgers University Alejandro Strachan Materials Engineering, Purdue
More informationSurface physics, Bravais lattice
Surface physics, Bravais lattice 1. Structure of the solid surface characterized by the (Bravais) lattice + space + point group lattice describes also the symmetry of the solid material vector directions
More informationHydrogen Embrittlement,, Microcracking, and Piezonuclear Fission Reactions at the Ni and Pd Electrodes of Electrolysis Cold. Fusion Experiments
Hydrogen Embrittlement,, Microcracking, and Piezonuclear Fission Reactions at the Ni and Pd Electrodes of Electrolysis Cold Fusion Experiments A. Carpinteri, O. Borla,, A. Goi, A. Manuello, D. Veneziano
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationM02M.1 Particle in a Cone
Part I Mechanics M02M.1 Particle in a Cone M02M.1 Particle in a Cone A small particle of mass m is constrained to slide, without friction, on the inside of a circular cone whose vertex is at the origin
More information1 Bulk Simulations in a HCP lattice
1 Bulk Simulations in a HCP lattice 1.1 Introduction This project is a continuation of two previous projects that studied the mechanism of melting in a FCC and BCC lattice. The current project studies
More informationN = N A Pb A Pb. = ln N Q v kt. = kt ln v N
5. Calculate the energy for vacancy formation in silver, given that the equilibrium number of vacancies at 800 C (1073 K) is 3.6 10 3 m 3. The atomic weight and density (at 800 C) for silver are, respectively,
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationA Quasicontinuum for Complex Crystals
A Quasicontinuum for Complex Crystals Ellad B. Tadmor Aerospace Engineering and Mechanics University of Minnesota Collaborators: Previous Work: Current Work: U. V. Waghmare, G. S. Smith, N. Bernstein,
More informationField Method of Simulation of Phase Transformations in Materials. Alex Umantsev Fayetteville State University, Fayetteville, NC
Field Method of Simulation of Phase Transformations in Materials Alex Umantsev Fayetteville State University, Fayetteville, NC What do we need to account for? Multi-phase states: thermodynamic systems
More informationConstitutive models: Incremental (Hypoelastic) Stress- Strain relations. and
Constitutive models: Incremental (Hypoelastic) Stress- Strain relations Example 5: an incremental relation based on hyperelasticity strain energy density function and 14.11.2007 1 Constitutive models:
More informationPhonons (Classical theory)
Phonons (Classical theory) (Read Kittel ch. 4) Classical theory. Consider propagation of elastic waves in cubic crystal, along [00], [0], or [] directions. Entire plane vibrates in phase in these directions
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 informationMolecular 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 informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationCombinatorial RF Magnetron Sputtering for Rapid Materials Discovery: Methodology and Applications
Combinatorial RF Magnetron Sputtering for Rapid Materials Discovery: Methodology and Applications Philip D. Rack,, Jason D. Fowlkes,, and Yuepeng Deng Department of Materials Science and Engineering University
More informationA MECHANISM FOR DUCTILE FRACTURE IN SHEAR. Viggo Tvergaard Department of Mechanical Engineering,
A MECHANISM FOR DUCTILE FRACTURE IN SHEAR Viggo Tvergaard Department of Mechanical Engineering, Solid Mechanics Technical University of Denmark, Bldg. 44, DK-8 Kgs. Lyngby, Denmark N.A. Fleck, J.W. Hutchinson
More information(a) What are the probabilities associated with finding the different allowed values of the z-component of the spin after time T?
1. Quantum Mechanics (Fall 2002) A Stern-Gerlach apparatus is adjusted so that the z-component of the spin of an electron (spin-1/2) transmitted through it is /2. A uniform magnetic field in the x-direction
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 informationATOMISTIC MODELING AND SIMULATION OF LONG-TERM TRANSPORT PHENOMENA IN NANOMATERIALS
IS Atomistic - Evolving modeling Microstructures and simulation in Solids of long-term transport phenomena in nanomaterials XIII International Conference on Computational Plasticity. Fundamentals and Applications
More informationBulk Structures of Crystals
Bulk Structures of Crystals 7 crystal systems can be further subdivided into 32 crystal classes... see Simon Garrett, "Introduction to Surface Analysis CEM924": http://www.cem.msu.edu/~cem924sg/lecturenotes.html
More informationStatistical 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 informationNanoscale 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 informationModèle de liaisons fortes au 4ème moment pour traiter l ordre-désordre dans les alliages
Modèle de liaisons fortes au 4ème moment pour traiter l ordre-désordre dans les alliages Jan Los, Christine Mottet, Guy Tréglia CINaM, Marseille Christine Goyhenex IPCMS, Strasbourg Outline Context Tight
More informationIn-class exercises. Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 8 Exercises due Mon March 19 Last correction at March 5, 2018, 8:48 am c 2017, James Sethna,
More informationMolecular Aggregation
Molecular Aggregation Structure Analysis and Molecular Simulation of Crystals and Liquids ANGELO GAVEZZOTTI University of Milano OXFORD UNIVERSITY PRESS Contents PART I FUNDAMENTALS 1 The molecule: structure,
More informationUnderstanding Phonon Dynamics via 1D Atomic Chains
Understanding Phonon Dynamics via 1D Atomic Chains Timothy S. Fisher Purdue University School of Mechanical Engineering, and Birck Nanotechnology Center tsfisher@purdue.edu Nanotechnology 501 Lecture Series
More informationChapter 19: Elastic Inclusions
Elastic Inclusions J.W. Morris Jr. January 8 Chapter 9: Elastic Inclusions CHAPTER 9: ELASTIC INCLUSIONS 9. Introduction 9. The Eshelby Cycle 3 9.3 The Macroscopic Strain 6 9.3. Boundary conditions; microscopic
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 5, April 14, 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationLecture contents. A few concepts from Quantum Mechanics. Tight-binding model Solid state physics review
Lecture contents A few concepts from Quantum Mechanics Particle in a well Two wells: QM perturbation theory Many wells (atoms) BAND formation Tight-binding model Solid state physics review Approximations
More informationPhysics of disordered materials. Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University
Physics of disordered materials Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University Course plan Familiarity with the basic description of disordered structures
More informationRatcheting and Rolling Contact Fatigue Crack Initiation Life of Rails under Service Loading. Wenyi YAN Monash University, Australia
Ratcheting and Rolling Contact Fatigue Crack Initiation Life of Rails under Service Loading Wenyi YAN Monash University, Australia Chung Lun PUN Peter Mutton Qianhua Kan Guozheng Kang Contents Introduction
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More information