Growing Microstructures using Phase-Field Crystal
|
|
- Neil Fox
- 6 years ago
- Views:
Transcription
1 Growing Microstructures using Phase-Field Crystal Stefan Bringuier 1 1 University of Arizona, Department of Materials Science and Engineering December 16, Overview The phase-field (PF) and phase-field-crystal (PFC) methods are relatively new approaches to modelling materials based on the variational principles of minimizing the free energy of a given system. More specifically, PFC is concerned with minimizing the free-energy functional by considering density fields that have periodic spatial variation 1. Thus for crystalline solids this inherently captures the periodic description. Furthermore such model naturally incorporates elastic/plastic deformations and captures varying grain boundary orientations. Computational investigation of microstructure evolution using PF to date has been well studied 2 4. The PF method makes use of field variables such as impurity concentration and temperature to evolve the system microstructure. This has been extensively used to study dendritic growth in metals and alloys as well as spinodal decomposition. The dynamics of the system are driven by the dissipative minimization of a phenomenological determined free-energy functional. One of the limitations concerning phase field modelling is that the formulation requires uniform fields that are in equilibrium. As a result of this properties such as elasticity, anisotropy, and grain orientations are not inherent in the physical description. Although this has been addressed via different methods, PF has yet been unable to handle diffusive phase transformations and anisotropic surface energy to name a few 1,5. More recently the PFC method has shown promise in addressing the short-comings of PF. The primary difference between PFC and PF is that the order parameter is refashioned as a temporally coarse but spatially described by an atomic probability density (APD) 5. This is done by choosing stefanb@ .arizona.edu 1
2 a free-energy functional which is minimized by an APD that is spatially periodic for crystalline solids and constant for liquids. It so happens that such description incorporates the necessary physics to describe things like surface tension, elasticity, anisotropy, dislocations etc. Additionally diffusive behaviour can be studied as a consequence of the time scale. In Fig. 1 an example from Provatas and Elder 6 is showing the ADP field of a partially solidified simulation. Many aspects such as orientation, liquid/solid interface, and dislocations are highlighted. Figure 1: PFC simulation showing ADP field and the various micro structure properties captured. 6 2 Proposed Research The importance of structural information such as grain boundary (GB) interfaces is crucial to understand the properties of a given material. For this reason computational methods such as molecular dynamics (MD) and first principle calculations are utilized to subject these interfaces to varies conditions. However, due to constraints from time and spatial scales it is computational difficult to grow crystals/polycrystals and furthermore interfaces using MD or first principle calculations. In some cases the use of mathematical geometric constructions make ease in building GB, an example of such approach are the coincident-site lattice model and Voronoi tessellation method. Although very useful, such methods are limiting and do not necessarily provide agreement with experimental microstructures. In addition even after considerable annealing and equilibration cycles the resulting interfaces might still reside in local minima. Thus determining the most suitable initial configuration requires some prior knowledge from 2
3 experimental or other computational methods. Phase-field-crystal may provide some value insight since it is possible to grow microstructures and interfaces.these structures can then be subjected to external fields or non-equilibrium events to extract properties 7. For example in Fig. 2 the dislocation glide velocity and crack propagation behaviours have been study. Figure 2: Examples of dislocation glide during shear and crack growth using the PFC method 7, 1 However there are situations where atomic length and time scales are needed and therefore MD might be better suited. An example of this is the effect of grain boundaries on ion bombardment. The benefit here is that PFC method has generated a blurred atomic picture but to some degree can be used to represent the atomic positions. The resulting microstructures from PFC can be used by MD simulations, this is a rare example where information across scales occurs in reverse. An example of a 3D generated polycrystalline BCC sample using PFC is shown in 3. Figure 3: 3D nucleation and growth of a BCC material where the purple indicates atomic coordinates with correct coordination and yellow otherwise. 8 3
4 The benefit to using PFC rather than other techniques such as Voronoi construction or coincident-stie lattice model is the correctness of the geometry generated by PFC. This is because the variational backbone accounts for minimizing free energy associated with interfaces. As a result mechanism associated with dislocation movement and interaction in polycrystalline materials can be observed. As a side note, the construction of phase diagrams is readily available since the free energy descriptions are used. 2.1 Background The first step is to determine the a free-energy functional which is minimized by a periodic function. If we consider the following equation: F = d r (ψ G2 ψ + u4 ) ψ4 (1) Here G is an operator by which when larger than zero has a single well and smaller than zero two wells. This is similar to other functionals such as Model A and Model B 6. The key point is that when ψ is a periodic function G can be chosen such that it is negative and when ψ is uniform it is positive. This type of problem was tackled by Swift and Hohenberg to study Rayleigh-Bénard convection of a fluid between a hot and cold plate. It was shown, now known as the Swift-Hohenberg model, that the free-energy functional for periodic systems is: F = [ ψ d r 2 ( ɛ + ( go 2 + 2) ) ] 2 ψ + ψ4 4 This equation sets much of the need framework for the formulation for PFC. I will not provide the full derivation on how to go from the Swift- Hohenberg model to PFC since it is still under investigation by the author. However I refer the reader to the book by First Steps As part of demonstrating capability to carry out the proposed work we have authored a python code based on GNU license libraries which simulates evolution of a non-convserved field (e.g. magnetization) using the Cahn- Allen or Model A. This code is based on the steps outlined in the textbook Phase-Field Methods in Materials Science and Engineering by Provatas and Elder 6. Here the phase-field variable is assumed to be a non-conserved quantity and the free energy functional is given by: (2) f(φ) = a(t ) + a 2 2 φ2 + a 4 4 φ4 + O(φ 6 ) (3) The dynamics of this system is of Langevin type and is given by: 4
5 ( φ t = M W o 2 φ f(φ, T ) ) φ Here M is related to the time scale change from disordered to order phase. Inserting Eq. 3 into Eq. 4 gives the following which describes the dynamics of a non-conserved field in the absence of an external field: (4) φ t = M ( W o 2 φ a 2 φ a 4 φ ) (5) Figure 4: Order to disordered phase transition at different timesteps using Chan-Allen/Model A approach. The source code and parameters used are provided in the appendix. In Fig. 4 microstructure evolution is shown at varies timesteps. We see that the system decomposes into two separate and distinguishable phases. This is an example spinodal decomposition. 3 Design of Simulation Although PFC has shown great promise to date it has been used on fairly simple systems. To better understand its capabilities the author suggest a possible investigation into understanding how to represent materials which are classified as chemically complex and mechanically brittle. The intial work would be to come up with a free-energy functional using a selected prototype system and extend then extend it to other materials. Once determined the relative dynamics can be worked out and PFC simulations can be regularly carried out. 5
6 References [1] K. R. Elder and M. Grant, Modeling elastic and plastic deformations in non-equilibrium processing using phase field crystals, arxiv:cond-mat/ , Jun [Online]. Available: [2] J. Warren and W. Boettinger, Prediction of dendritic growth and microsegregation patterns in a binary alloy using the phase-field method, Acta Metallurgica et Materialia, vol. 43, no. 2, pp , Feb [Online]. Available: [3] S.-L. Wang and R. F. Sekerka, Computation of the dendritic operating state at large supercoolings by the phase field model, Physical Review E, vol. 53, no. 4, pp , Apr [Online]. Available: [4] N. Provatas, P. J. A. Dantzig, B. Athreya, P. Chan, P. Stefanovic, N. Goldenfeld, and K. R. Elder, Using the phase-field crystal method in the multi-scale modeling of microstructure evolution, JOM, vol. 59, no. 7, pp , Jul [Online]. Available: [5] K. R. Elder, N. Provatas, J. Berry, P. Stefanovic, and M. Grant, Phasefield crystal modeling and classical density functional theory of freezing, Physical Review B, vol. 75, no. 6, p , Feb [Online]. Available: [6] N. Provatas and K. Elder, Phase-Field Methods in Materials Science and Engineering. John Wiley & Sons, Sep [7] P. Stefanovic, M. Haataja, and N. Provatas, Phase-field crystals with elastic interactions, Physical Review Letters, vol. 96, no. 22, p , Jun [Online]. Available: [8] W. J. Boettinger, J. A. Warren, C. Beckermann, and A. Karma, Phase-field simulation of solidification1, Annual Review of Materials Research, vol. 32, no. 1, pp , [Online]. Available: 4 Appendix 4.1 Model A Code 6
7 #! / usr / bin / python # ModelA Based on t he code in Phase F i e l d Methods in M a t e r i a l s Science # by N. Provatas and K. Elder. In order to have comparable performance to # o f numpexpr i s implemented. Also note t h a t numpy. r o l l ( ) which makes a c # Author Stefan Bringuier # 11/15/2013 # TODO Try p a r a l l i z a t i o n v i a MPI4py and PyCUDA import numpy a s np import numexpr a s nex import s c i p y as sp import mpi4py a s mpi import m a t p l o t l i b. pyplot as p l t import time class V a r i a b l e s ( ) : Class to d e f i n e v a r i a b l e s and arrays a s s o c i a t e d with Model A. def i n i t ( s e l f, kwargs ) : #D e f a u l t run v a l u e s s e l f. dx = 0. 8 s e l f. dy = 0. 8 s e l f. dt = 0. 1 s e l f. a 2 = 1.0 s e l f. a 4 = 1. 0 s e l f.w2 = 0.25 s e l f. tmax = s e l f. Nx = 256 s e l f. Ny = 256 s e l f.m = 1. 0 s e l f. cn = s e l f. dx2 in = 1. 0 / ( ( s e l f. dx 2 + s e l f. dy 2)/2) s e l f. a v e p s i = 0. 0 #Check f o r f i l e to read from #t r y : 7
8 # s e l f. ReadFile ( kwargs [ FileName ] ) #e x c e p t : # p r i n t No f i l e to read. Using D e f a u l t v a l u e s! s e l f. A l l o c a t e ( ) def ReadFile ( F i l e ) : #f i l e = open ( s t r ( FileName ), ro ) print read f i l e def A l l o c a t e ( s e l f ) : s e l f. PSI = np. z e r o s ( ( s e l f. Nx, s e l f. Ny), dtype=f l o a t ) s e l f. Grad2 = np. z e r o s ( ( s e l f. Nx, s e l f. Ny), dtype=f l o a t ) s e l f.rhs = np. z e r o s ( ( s e l f. Nx, s e l f. Ny), dtype=f l o a t ) s e l f. Grad2Mu = np. z e r o s ( ( s e l f. Nx, s e l f. Ny), dtype=f l o a t ) def D e a l l o c a t e ( s e l f ) : s e l f. PSI. d e l e t e ( ) s e l f. Grad2. d e l e t e ( ) s e l f.rhs. d e l e t e ( ) s e l f. Grad2Mu. d e l e t e ( ) class U t i l i t y : Set o f f u n c t i o n s t h a t do not need to be i n s t a n t i a t e def Print2DFile ( f i e l d, time ) : np. s a v e t x t ( f i e l d 2 D +time+. out, f i e l d, fmt= %23.16 f def RandNumArry( nx, ny ) : a = np. random. random sample ( ( nx, ny ) ) return def GaussArry ( nx, ny ) : mu = 0. 0 sigm = 0. 5 a =np. random. normal (mu, sigm, ( nx, ny ) ) return def NABALA2(A, dx2 in, num proc = 1 ) : 8
9 nex. set num threads ( num proc ) fwd = np. r o l l (A, 1, 1 ) up = np. r o l l (A, 1, 0 ) fwdup = np. r o l l ( fwd, 1, 0 ) bck = np. r o l l (A, 1,1) dwn = np. r o l l (A, 1,0) bckdwn = np. r o l l ( bck, 1,0) fwddwn = np. r o l l ( fwd, 1,0) bckup = np. r o l l ( bck, 1, 0 ) # See d e s c r i p t i o n in book f o r l a p l a c i a n d e r i v a t i o n GR2 = nex. e v a l u a t e ( ( fwd + bck + up + dwn) / \ ( fwdup + bckup + fwddwn + bckdwn ) / 4. 0 ( 3. 0 A dx2 in ) ) return def RightHS (V, num proc =1): W2 = V.W2 Grad2 = V. Grad2 a 2 = V. a 2 PSI = V. PSI a 4 = V. a 4 RHS = nex. e v a l u a t e ( W2 Grad2 a 2 PSI a 4 ( PSI 3) ) return RHS i f name == m a i n : v = V a r i a b l e s ( FileName = T e s t F i l e ) v. tmax = # Use u t i l i t y to g e n erate random phi f i e l d or # read in f i l e. 9
10 v. PSI = U t i l i t y. GaussArry ( v. Nx, v. Ny) #v. PSI = np. l o a d t x t ( input. c f g ). reshape ( v. Nx, v. Ny) threads = nex. d e t e c t n u m b e r o f c o r e s ( ) # Evolve s t r u c t u r e s t a r t = time. time ( ) for t in range ( v. tmax ) : v. Grad2 = U t i l i t y.nabala2( v. PSI, v. dx2 in, num proc=threads ) v.rhs = U t i l i t y. RightHS ( v, num proc=threads ) v. PSI += ( v. dt v.m v.rhs) i f t %5000 == 0 : t t = t /5000 p l t. imshow ( v. PSI ) p l t. c o l o r b a r ( ) p l t. s a v e f i g ( image +s t r ( t t )+. jpg ) p l t. c l f ( ) print Time :, time. time ( ) s t a r t 10
Density 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 informationPhase-field crystal modeling and classical density functional theory of freezing
Phase-field crystal modeling and classical density functional theory of freezing K. R. Elder, 1 Nikolas Provatas, 2 Joel Berry, 1,3 Peter Stefanovic, 2 and Martin Grant 3 1 Department of Physics, Oakland
More informationPhase Field Crystal (PFC) Model and Density Functional Theory (DFT) of Freezing
Phase Field Crystal (PFC) Model and Density Functional Theory (DFT) of Freezing Pyrite Project Meeting October 14 th 2010 Arvind Baskaran John Lowengrub Density functional Theory of Freezing [ Ramakrishnan
More informationAn In-Depth Examination of a Thermodynamic Framework for the Phase-Field Crystal Model
An In-Depth Examination of a Thermodynamic Framework for the Phase-Field Crystal Model by Victor W. L. Chan A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor
More informationAmplitude expansion of the binary phase-field-crystal model
Amplitude expansion of the binary phase-field-crystal model K. R. Elder Department of Physics, Oakland University, Rochester, Michigan 4809, USA Zhi-Feng Huang Department of Physics and Astronomy, Wayne
More informationErrata for SOLIDIFICATION (Second Edition, 2016)
Errata for SOLIDIFICATION (Second Edition, 2016) J. A. Dantzig and M. Rappaz September 6, 2017 Nomenclature There are several minor typographical errors in this table. Please download the corrected version
More informationSIMULATION OF DENDRITIC CRYSTAL GROWTH OF PURE Ni USING THE PHASE-FIELD MODEL
46 Rev. Adv. Mater. Sci. 33 (13) 46-5 Yu. Zhao and H. Hou SIMULATION OF DENDRITIC CRYSTAL GROWTH OF PURE Ni USING THE PHASE-FIELD MODEL Yuhong Zhao and Hua Hou College of Materials Science & Engineering,
More informationDresden, Dresden, Germany Published online: 09 Jan 2009.
This article was downloaded by: [SLUB Dresden] On: 11 December 2013, At: 04:59 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationRenormalization Group Approach to Multiscale Modelling in Materials Science
Journal of Statistical Physics, Vol. 125, Nos. 5/6, December 2006 C 2006 ) DOI: 10.1007/s10955-005-9013-7 Renormalization Group Approach to Multiscale Modelling in Materials Science Nigel Goldenfeld, 1
More informationPatterns, universality and computational algorithms
Patterns, universality and computational algorithms Nigel Goldenfeld Department of Physics University of Illinois at Urbana-Champaign Collaborators: Yoshi Oono, M. Mondello (Cell dynamical systems) Navot
More informationarxiv: v1 [cond-mat.mtrl-sci] 26 Oct 2012
A Phase Field Crystal study of Solute Trapping Harith Humadi and Jeffrey J. Hoyt Department of Materials Science and Engineering and Brockhouse Institute for Materials Research, McMaster University, 1280
More informationModeling structural transformations in binary alloys with phase field crystals
PHYSICAL REVIEW B 84,064104(2011) Modeling structural transformations in binary alloys with phase field crystals Michael Greenwood, 1,2 Nana Ofori-Opoku, 2 Jörg Rottler, 1 and Nikolas Provatas 2 1 Department
More informationSpinodal decomposition
Part V. Solid-solid transformations II Module 1 : Spinodal decomposition Part V. Solid-solid transformations II In this part, we discuss two transformations namely spinodal decomposition and ordering.
More informationSIMULATION OF POLYMER CRYSTAL GROWTH WITH VARIOUS MORPHOLOGIES USING A PHASE-FIELD MODEL
SIMULATION OF POLYMER CRYSTAL GROWTH WITH VARIOUS MORPHOLOGIES USING A PHASE-FIELD MODEL M. Asle Zaeem,*, S. Nouranian, Mark F. Horstemeyer Department of Materials Science and Engineering Missouri University
More informationErrata for SOLIDIFICATION (First Edition, 2009)
Errata for SOLIDIFICATION (First Edition, 29) J. A. Dantzig and M. Rappaz March 2, 217 Chapter 1: Overview Page 2, first line: the definition of r should have a square root. i.e., r = (ξ 2 + y 2 + z 2
More informationMODELLING MICROSTRUCTURE EVOLUTION IN MATERIALS SCIENCE
MODELLING MICROSTRUCTURE EVOLUTION IN MATERIALS SCIENCE PHASE-FIELD-CRYSTAL MODELS FOR MICROSTRUCTURAL EVOLUTION AND PHASE SELECTION IN MATERIALS SCIENCE By NANA OFORI-OPOKU B.Eng Materials Science and
More information4. STRUCTURES AND DEFECTS IN CRYSTALS
53 4. STRUCTURES AND DEFECTS IN CRYSTALS 4.1 Review and Introduction In the last three introductory chapters, we discussed the difference between variability and reliability and three types reliability
More informationMesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai
Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations
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 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 informationMULTISCALE MODELING OF SOLIDIFICATION: PHASE-FIELD METHODS TO ADAPTIVE MESH REFINEMENT
International Journal of Modern Physics B Vol. 19, No. 31 (2005) 4525 4565 c World Scientific Publishing Company20 December 2005 MULTISCALE MODELING OF SOLIDIFICATION: PHASE-FIELD METHODS TO ADAPTIVE MESH
More informationLAMMPS Simulation of a Microgravity Shear Cell 299r Progress Report 2 Taiyo Wilson. Introduction:
Introduction: After our initial simulations outlined in the first progress report, we determined that a shear cell with a bidisperse composition of particles would be a promising system to explore. However,
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 6 () 77 745 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Multiphase
More informationPHASE-FIELD MODELS FOR MICROSTRUCTURE EVOLUTION
Annu. Rev. Mater. Res. 2002. 32:113 40 doi: 10.1146/annurev.matsci.32.112001.132041 Copyright c 2002 by Annual Reviews. All rights reserved PHASE-FIELD MODELS FOR MICROSTRUCTURE EVOLUTION Long-Qing Chen
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationA simple model for block-copolymer layer morphology formation
A simple model for block-copolymer layer morphology formation A. J. Wagner and A. Croll, Department of Physics, North Dakota State University, Fargo, ND 58102 Prato Workshop, February 12, 2016 A big thanks
More informationUniversity of Bath DOI: / Publication date: Document Version Peer reviewed version. Link to publication
Citation for published version: Kundin, J, Emmerich, H & Zimmer, J 2011, 'Mathematical concepts for the micromechanical modelling of dislocation dynamics with a phase-field approach' Philosophical Magazine,
More informationAdaptive Mesh Refinement Computation of Solidification Microstructures Using Dynamic Data Structures
Journal of Computational Physics 148, 265 290 (1999) Article ID jcph.1998.6122, available online at http://www.idealibrary.com on Adaptive Mesh Refinement Computation of Solidification Microstructures
More informationINTERFACE THERMODYNAMICS WITH APPLICATIONS TO ATOMISTIC SIMULATIONS
INTERFACE THERMODYNAMICS WITH APPLICATIONS TO ATOMISTIC SIMULATIONS by Timofey Frolov A Dissertation Submitted to the Graduate Faculty of George Mason University in Partial Fulfillment of The Requirements
More informationStudy of Ni/Al Interface Diffusion by Molecular Dynamics Simulation*
Engineering, 2011, 3, 227-232 doi:10.4236/eng.2011.33026 Published Online March 2011 (http://www.scirp.org/journal/eng) Study of Ni/Al Interface Diffusion by Molecular Dynamics Simulation* Chunguang Zhang,
More informationVIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition
VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium
More informationMultiscale modelling of microstructure formation in polymer casting
TECHNISCHE MECHANIK, 30, 1-3, (2010), 259 268 submitted: October 31, 2009 Multiscale modelling of microstructure formation in polymer casting T. Pitkänen, S. Majaniemi, T. Ala-Nissila A data bank approach
More informationA STOCHASTIC PHASE-FIELD MODEL DETERMINED FROM MOLECULAR DYNAMICS
A STOCHASTIC PHASE-FIELD MODEL DETERMINED FROM MOLECULAR DYNAMICS ERIK VON SCHWERIN AND ANDERS SZEPESSY Abstract. The dynamics of dendritic growth of a crystal in an undercooled melt is determined by macroscopic
More informationc 2006 by Badrinarayan P. Athreya. All rights reserved
c 2006 by Badrinarayan P. Athreya. All rights reserved PHASE-FIELDS AND THE RENORMALIZATION GROUP: A CONTINUUM APPROACH TO MULTISCALE MODELING OF MATERIALS BY BADRINARAYAN P. ATHREYA B. E., University
More informationThermodynamically Consistent Algorithms for the Solution of Phase-Field Models. Dissertation by Philippe Antoine Vignal Atherton
Thermodynamically Consistent Algorithms for the Solution of Phase-Field Models Dissertation by Philippe Antoine Vignal Atherton In Partial Fulfillment of the Requirements For the Degree of Doctor of Philosophy
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationMagnetic islands modelled by a phase-field-crystal approach
Eur. Phys. J. B (018) 91: 55 https://doi.org/10.1140/epjb/e018-80543-9 Regular Article THE EUROPEAN PHYSICAL JOURNAL B Magnetic islands modelled by a phase-field-crystal approach Niloufar Faghihi 1,a,
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 informationSpatiotemporal Chaos in Rayleigh-Bénard Convection
Spatiotemporal Chaos in Rayleigh-Bénard Convection Michael Cross California Institute of Technology Beijing Normal University June 2006 Janet Scheel, Keng-Hwee Chiam, Mark Paul Henry Greenside, Anand Jayaraman
More informationAnalysis of the Efficiency PETSc and PETIGA Libraries in Solving the Problem of Crystal Growth
Analysis of the Efficiency PETSc and PETIGA Libraries in Solving the Problem of Crystal Growth Ilya Starodumov 1, Evgeny Pavlyuk 1, Leonid Klyuev 2, Maxim Kovalenko 3, and Anton Medyankin 1 1 Ural Federal
More informationPhase-field modeling of the dendrite orientation transition in Al-Zn alloys
Phase-field modeling of the dendrite orientation transition in Al-Zn alloys Jonathan Friedli, Paolo Di Napoli, Michel Rappaz and Jonathan A. Dantzig Laboratoire de simulation des matériaux, Institut des
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationCoarsening process in the 2d voter model
Alessandro Tartaglia (LPTHE) Coarsening in the 2d voter model May 8, 2015 1 / 34 Coarsening process in the 2d voter model Alessandro Tartaglia LPTHE, Université Pierre et Marie Curie alessandro.tartaglia91@gmail.com
More informationEquilibria in Materials
009 fall Advanced Physical Metallurgy Phase Equilibria in Materials 10. 13. 009 Eun Soo Park Office: 33-316 Telephone: 880-71 Email: espark@snu.ac.kr Office hours: by an appointment 1 Contents for previous
More informationA Diffuse Interface Model of Grain Boundary Faceting. F. Abdeljawad, 1 D. L. Medlin, 2 J. A. Zimmerman, 2 K. Hattar, 1 and S. M.
A Diffuse Interface Model of Grain Boundary Faceting F. Abdeljawad, 1 D. L. Medlin, 2 J. A. Zimmerman, 2 K. Hattar, 1 and S. M. Foiles 1 1) Sandia National Laboratories, Albuquerque, New Mexico 87185,
More informationTopological Defects inside a Topological Band Insulator
Topological Defects inside a Topological Band Insulator Ashvin Vishwanath UC Berkeley Refs: Ran, Zhang A.V., Nature Physics 5, 289 (2009). Hosur, Ryu, AV arxiv: 0908.2691 Part 1: Outline A toy model of
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 informationRheology. What is rheology? From the root work rheo- Current: flow. Greek: rhein, to flow (river) Like rheostat flow of current
Rheology What is rheology? From the root work rheo- Current: flow Greek: rhein, to flow (river) Like rheostat flow of current Rheology What physical properties control deformation? - Rock type - Temperature
More informationA sharp diffuse interface tracking method for approximating evolving interfaces
A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex Overview Introduction Phase field models Double well and double obstacle
More informationPHASE FIELD MODELS, ADAPTIVE MESH REFINEMENT AND LEVEL SETS FOR SOLIDIFICATION PROBLEMS
PHASE FIELD MODELS, ADAPTIVE MESH REFINEMENT AND LEVEL SETS FOR SOLIDIFICATION PROBLEMS Department of Physics University of Illinois at Urbana-Champaign COWORKERS AND COLLABORATORS Phase-field calculations
More informationPhase-Field Models. Mathis Plapp
Phase-Field Models Mathis Plapp Physique de la Matière Condensée, École Polytechnique, CNRS, 91128 Palaiseau, France Abstract. Phase-field models have become popular in recent years to describe a host
More informationLAMMPS Simulation of a Microgravity Shear Cell 299r Progress Report Taiyo Wilson. Units/Parameters:
Units/Parameters: In our simulations, we chose to express quantities in terms of three fundamental values: m (particle mass), d (particle diameter), and τ (timestep, which is equivalent to (g/d)^0.5, where
More informationOn the nonlinear anelastic behaviour of AHSS
Journal of Physics: Conference Series PAPER OPEN ACCESS On the nonlinear anelastic behaviour of AHSS To cite this article: A Torkabadi et al 2016 J. Phys.: Conf. Ser. 734 032100 Related content - Essential
More informationUniversity of Michigan, Ann Arbor, Michigan, 48109, USA. University of Michigan, Ann Arbor, Michigan, 48109, USA. (Dated: December 11, 2012) Abstract
Modeling Classical Density Functional Theory and the Phase-Field Crystal Method using a Rational Function to Describe the Two-body Direct Correlation Function N. Pisutha-Arnond,, V. W. L. Chan,, M. Iyer,,
More informationDendritic solidification of binary alloys
University of Ljubljana Faculty of Mathematics and Physics Department of Physics Seminar I b - 3rd year, 2nd cycle Dendritic solidification of binary alloys Author: Tadej Dobravec Mentor: prof. dr. Božidar
More informationPapers Cited >1000X GOOGLE SCHOLAR
Papers Cited >1000X GOOGLE SCHOLAR March 2019 Citations 60861 15529 h-index 111 57 i10-index 425 206 1. Title: Formation of glasses from liquids and biopolymers Source: Science, 1995 sciencemag.org Abstract
More informationPhase field approach to heterogeneous crystal nucleation in alloys
Phase field approach to heterogeneous crystal nucleation in alloys James A. Warren, 1 Tamás Pusztai, 2 László Környei, 2 and László Gránásy 3 1 Metallurgy Division, National Institute of Standards and
More informationREVIEW ARTICLE. Phase-field-crystal models for condensed matter dynamics on atomic length and diffusive time scales: an overview
Advances in Physics Vol. 61, No. 6, November December 2012, 665 743 REVIEW ARTICLE Phase-field-crystal models for condensed matter dynamics on atomic length and diffusive time scales: an overview Heike
More informationRelaxation of a Strained Elastic Film on a Viscous Layer
Mat. Res. Soc. Symp. Proc. Vol. 695 Materials Research Society Relaxation of a Strained Elastic Film on a Viscous Layer R. Huang 1, H. Yin, J. Liang 3, K. D. Hobart 4, J. C. Sturm, and Z. Suo 3 1 Department
More informationPhysics 541: Condensed Matter Physics
Physics 541: Condensed Matter Physics In-class Midterm Exam Wednesday, October 26, 2011 / 14:00 15:20 / CCIS 4-285 Student s Name: Instructions There are 23 questions. You should attempt all of them. Mark
More informationPHASE-FIELD MODELLING OF NONEQUILIBRIUM PARTITIONING DURING RAPID SOLIDIFICATION IN A NON-DILUTE BINARY ALLOY. Denis Danilov, Britta Nestler
First published in: DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 15, Number 4, August 2006 pp. 1035 1047 PHASE-FIELD MODELLING OF NONEQUILIBRIUM PARTITIONING DURING
More informationAdaptive 3D finite elements with high aspect ratio for dendritic growth of a binary alloy including fluid flow induced by shrinkage
Copyright c 2007 Tech Science Press FDMP, vol.3, no.1, pp.49-64, 2007 Adaptive 3D finite elements with high aspect ratio for dendritic growth of a binary alloy including fluid flow induced by shrinkage
More informationFluid Model of Crystal Plasticity - Mathematical Properties and Computer Simulations
Fluid Model of Crystal Plasticity - Mathematical Properties and Computer Simulations Piotr Minakowski An extended abstract of the PhD dissertation 1 Introduction The main purpose of this dissertation is
More informationLecture 15: Phase Field Modelling
Materials Science & Metallurgy Master of Philosophy, Materials Modelling, Course MP6, Kinetics and Microstructure Modelling, H. K. D. H. Bhadeshia Lecture 15: Phase Field Modelling Introduction Imagine
More informationSTRESS UPDATE ALGORITHM FOR NON-ASSOCIATED FLOW METAL PLASTICITY
STRESS UPDATE ALGORITHM FOR NON-ASSOCIATED FLOW METAL PLASTICITY Mohsen Safaei 1, a, Wim De Waele 1,b 1 Laboratorium Soete, Department of Mechanical Construction and Production, Ghent University, Technologiepark
More informationSupersolid. Yixing Fu 2013/3/14
Supersolid Yixing Fu 2013/3/14 Supersolid "Superfluid-like solid" Intensively discussed by theorists in 1970s Still uncertain 3/14/2013 Yixing Fu - Physics 141A 2 Outline The pre-history of supersolid
More informationPhenomenological Theories of Nucleation
Chapter 1 Phenomenological Theories of Nucleation c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 16 September 2012 1.1 Introduction These chapters discuss the problems of nucleation, spinodal
More informationSTRONG CONFIGURATIONAL DEPENDENCE OF ELASTIC PROPERTIES OF A CU-ZR BINARY MODEL METALLIC GLASS
Chapter 3 STRONG CONFIGURATIONAL DEPENDENCE OF ELASTIC PROPERTIES OF A CU-ZR BINARY MODEL METALLIC GLASS We report the strong dependence of elastic properties on configurational changes in a Cu-Zr binary
More informationLecture 6 Surface Diffusion Driven by Surface Energy
Lecture 6 Surface Diffusion Driven by Surface Energy Examples Flattening a surface. Spherodizing. Rayleigh instability. Grain boundary grooving. Sintering Self-assembled quantum dots Atomic Flux and Surface
More informationThermal properties of Engineering Materials
Thermal properties of Engineering Materials Engineering materials are important in everyday life because of their versatile structural properties. Other than these properties, they do play an important
More informationAn Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation
An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University
More informationConvective Heat and Mass Transfer Prof. A. W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay
Convective Heat and Mass Transfer Prof. A. W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay Module No.# 01 Lecture No. # 41 Natural Convection BLs So far we have considered
More informationc 2009 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 47 No 3 pp 69 88 c 009 Society for Industrial and Applied Mathematics AN ENERGY-STABLE AND CONVERGENT FINITE-DIFFERENCE SCHEME FOR THE PHASE FIELD CRYSTAL EQUATION S M WISE C WANG
More informationContact melting and the structure of binary eutectic near the eutectic point
Contact melting and the structure of binary eutectic near the eutectic point 1, ystrenko O.V. and Kartuzov V.V. 1 1 Frantsevich Institute for material sсience problems, Kiev, Ukraine, ogolyubov Institute
More informationHIGHER ORDER OPERATOR SPLITTING FOURIER SPECTRAL METHODS FOR THE ALLEN CAHN EQUATION
J. KSIAM ol., No., 6, 7 http://dx.doi.org/.94/jksiam.7.. HIGHER ORDER OPERATOR SPLITTING FOURIER SPECTRAL METHODS FOR THE ALLEN CAHN EQUATION JAEMIN SHIN, HYUN GEUN LEE, AND JUNE-YUB LEE 3 INSTITUTE OF
More informationPattern Formation and Spatiotemporal Chaos
Pattern Formation and Spatiotemporal Chaos - Chennai, 2004 1 Pattern Formation and Spatiotemporal Chaos Insights from Large Scale Numerical Simulations of Rayleigh-Bénard convection Collaborators: Mark
More informationPhase Field Approach to Heterogeneous Crystal Nucleation in Alloys
Phase Field Approach to Heterogeneous Crystal Nucleation in Alloys James A. Warren 1, Tamás Pusztai, László Környei, and László Gránásy 3 1 Metallurgy Division, National Institute of Standards and Technology,
More information(Crystal) Nucleation: The language
Why crystallization requires supercooling (Crystal) Nucleation: The language 2r 1. Transferring N particles from liquid to crystal yields energy. Crystal nucleus Δµ: thermodynamic driving force N is proportional
More informationHeat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection
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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationWhat we should know about mechanics of materials
What we should know about mechanics of materials 0 John Maloney Van Vliet Group / Laboratory for Material Chemomechanics Department of Materials Science and Engineering Massachusetts Institute of Technology
More informationMOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM
MOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM MOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM By YUNFEI BAI, B. Eng.
More informationSupporting Material. A Mechanistic Collective Cell Model for Epithelial Colony Growth and Contact Inhibition
Supporting Material A Mechanistic Collective Cell Model for Epithelial Colony Growth and Contact Inhibition S. Aland 1,, H. Hatzikirou 1,2, J. Lowengrub,, A. Voigt 1,2,4 1 Department of Mathematics, TU
More informationMicrostructural Randomness and Scaling in Mechanics of Materials. Martin Ostoja-Starzewski. University of Illinois at Urbana-Champaign
Microstructural Randomness and Scaling in Mechanics of Materials Martin Ostoja-Starzewski University of Illinois at Urbana-Champaign Contents Preface ix 1. Randomness versus determinism ix 2. Randomness
More informationL UNAM Université, ONIRIS, CNRS, GEPEA, UMR6144, Nantes, F-44322, France.
ing Heating During Batch Processing of Liquid Sample in a Single Mode Cavity S. Curet 1, F. Bellincanta Begnini, O. Rouaud 1, L. Boillereau 1 1 L UNAM Université, ONIRIS, CNRS, GEPEA, UMR6144, Nantes,
More informationAbstract. The phase field model is used to compute numerically the temporal evolution
Computations of dendrites in 3-D and comparison with microgravity experiments Y. B. Altundas and G. Caginalp University of Pittsburgh Pittsburgh, PA 15260 ybast@pitt.edu caginalp@pitt.edu Abstract. The
More informationPhase-field modeling of binary alloy solidification with coupled heat and solute diffusion
PHYSICAL REVIEW E 69, 051607 (004) Phase-field modeling of binary alloy solidification with coupled heat and solute diffusion J. C. Ramirez, 1 C. Beckermann, 1, * A. Karma, and H.-J. Diepers,3 1 Department
More informationPhase-Field Modeling of Nucleation in Solid-State Phase Transformations
JOM, Vol. 66, No. 8, 2014 DOI: 10.1007/s11837-014-1033-9 Ó 2014 The Minerals, Metals & Materials Society Phase-Field Modeling of Nucleation in Solid-State Phase Transformations TAE WOOK HEO 1,3 and LONG-QING
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationPart I : Preliminaries (Thermodynamics and Kinetics) Module 1: Gibbs free energy and equilibria
Part I : Preliminaries (Thermodynamics and Kinetics) Module 1: Gibbs free energy and equilibria Introduction Phase transformations are ubiquitous; every one of us have kept water in the refrigerator to
More informationGrowth of Secondary Dendrite Arms of Fe C Alloy during Transient Directional Solidification by Phase-field Method
, pp. 43 436 Growth of Secondary Dendrite Arms of Fe C Alloy during Transient Directional Solidification by Phase-field Method Yu XIE, 1) Hongbiao DONG 1) * and Jonathan DANTZIG,3) 1) Department of Engineering,
More informationDendritic Growth with Fluid Flow in Pure Materials
Dendritic Growth with Fluid Flow in Pure Materials Jun-Ho Jeong 1, Jonathan A. Dantzig 1 and Nigel Goldenfeld 2 1 Department of Mechanical and Industrial Engineering 2 Department of Physics University
More informationSTUDY ON CORROSION P HENOMENA OF STEELS IN PB-BI FLOW
11 th International Conference on Nuclear Engineering Tokyo, JAPAN, April -3, 3 ICONE11-36375 STUDY ON CORROSION P HENOMENA OF STEELS IN PB-BI FLOW Yingxia Qi Research Laboratory for Nuclear Reactors Tokyo
More informationChapter 10. Lesson Starter. Why did you not smell the odor of the vapor immediately? Explain this event in terms of the motion of molecules.
Preview Lesson Starter Objectives The Kinetic-Molecular Theory of Gases The Kinetic-Molecular Theory and the Nature of Gases Deviations of Real Gases from Ideal Behavior Section 1 The Kinetic-Molecular
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationNumerical modelling of ice & ice-structure interactions
Numerical modelling as defined in Encyclopædia Britannica: A computer generated description of a mathematical system to represent the behaviour of a real or proposed system that uses a set of equations
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationComputer simulation of spinodal decomposition in constrained films
Acta Materialia XX (2003) XXX XXX www.actamat-journals.com Computer simulation of spinodal decomposition in constrained films D.J. Seol a,b,, S.Y. Hu b, Y.L. Li b, J. Shen c, K.H. Oh a, L.Q. Chen b a School
More informationAnalytical Methods for Materials
Analytical Methods for Materials Lesson 11 Crystallography and Crystal Structures, Part 3 Suggested Reading Chapter 6 in Waseda Chapter 1 in F.D. Bloss, Crystallography and Crystal Chemistry: An Introduction,
More information