A sharp diffuse interface tracking method for approximating evolving interfaces
|
|
- Thomasine McCarthy
- 5 years ago
- Views:
Transcription
1 A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex
2 Overview Introduction Phase field models Double well and double obstacle potentials Applications of front tracking using the sharp diffuse interface method. Allen-Cahn equation (motion by mean curvature) Diffusion induced grain boundary motion (forced MBMC) Cahn-Hilliard equation (surface diffusion) Void electro migration Epitaxial growth Conclusions
3 Notation We are interested in tracking an evolving interface Γ t n Γ t evolving interface n unit normal V normal velocity κ mean curvature (negative when Γ t is a circle) s arc length d(x, t) signed distance function to Γ t (n = d, V = d t Γt, κ = d Γt )
4 Mathematical approaches Direct method parametric methods of differential geometry Γ t is parameterised by X(s, t) : [0, 1] [0, T ] IR 2 ( ) n = X s X s, V = X t X s X s, κ = 1 Xs X s X s Level set method Osher & Sethian 88 Γ t is taken to be the zero level set of a scalar ω. ( ) n = ω ω, V = ω t ω, κ = ω ω Phase field method Γ t is approximated by a diffuse transition layer in which an order parameter changes from 1 to +1. s
5 Phase field model To model motion by mean curvature V = κ Cahn & Allen 79 approximated Γ t by a diffuse transition layer. In this layer an order parameter ϕ changes from 1 to 1. The idea is based on the free energy ε F = 2 ϕ W (ϕ) (ε small) ε W ( ) is a double well potential with equal minima at ±1. The associated gradient flow of the free energy is εϕ t = ε ϕ + 1 ε W (ϕ) (Allen-Cahn equation)
6 Phase field model Setting P (ϕ) = εϕ t ε ϕ 1 ε W (ϕ) we have P (ψ(d/ε)) = ψ (d t d) 1 ε (ψ W (ψ)). Let ψ = ψ(z) be the unique solution of ψ + W (ψ) = 0 z IR, ψ(0) = 0, ψ (z) > 0, ψ(z) ±1 as z ± Then - (d t d) Γt = 0 - ψ ( d ε ) is only O(1) in an ε neighbourhood of Γ t - ψ ( d ε ) 0 away from Γ t So ψ( d ) is close to a solution of the Allen-Cahn equation. ε
7 Double Well Potential ϕ = 0 ϕ 1 ϕ W (ϕ) = 1 4 (1 ϕ2 ) 2 The diffuse interface is not well defined. Interface profile is ψ = tanh( d ε ) width O(ε log(ε) ).
8 Double Obstacle Potential ϕ < 1 ϕ = 1 ϕ = W (ϕ) = { 1 2 (1 ϕ2 ) ϕ 1 ϕ > 1 ϕ = 0 The diffuse interface is sharply defined. Interface profile is ψ = sin( d ε ), d ε < π 2 width O(ε).
9 Double Obstacle Potential Using the double obstacle potential W (ϕ) = { 1 2 (1 ϕ2 ) for ϕ 1 for ϕ > 1 the Allen-Cahn equation is now a parabolic variational inequality ε(ϕ t, η ϕ) + ε( ϕ, η ϕ) 1 (ϕ, η ϕ) ε η K = {η H 1 (Ω) : η 1}. Elliott & Chen 93, Nochetto, Paolini & Verdi 93 This problem approximates mean curvature flow.
10 Considering the gradient flow of the energy functional F = ε ϕ ε W (ϕ) + C Wfϕ Double Obstacle Potential yields ε(ϕ t, η ϕ) + ε( ϕ, η ϕ) ( ) ϕ ε C Wf, η ϕ η K C W = W (ϕ) = π 4 This converges to forced mean curvature flow V = κ + f. A finite element approximation of this model yields (for θ = 0, 1) ε t (ϕn ϕ n 1, χ ϕ n ) + ε( ϕ n θ, χ ϕ n ) ( 1 ε ϕn θ C W f n, χ ϕ n ) χ K h Convergence analysis Nochetto & Verdi 96, 97.
11 Double Obstacle Potential Phase field models can naturally deal with topological changes 1 t=0 1 t= t= f = 40 These computations and many others in this talk were obtained using the adaptive finite element toolbox ALBERT (Schmidt & Siebert 00).
12 Anisotropic Allen-Cahn with Forcing Considering the gradient flow of the energy functional F = εa( ϕ) + 1 W (ϕ) + fu ε yields an anisotropic Allen-Cahn equation with forcing ( ) ϕ ε(ϕ t, η ϕ) + ε(a ( ϕ), η ϕ) ε f, η ϕ η K A( ϕ) = γ2 2 (θ) ϕ 2 with tan θ = ϕ y ϕ x. Wheeler & McFadden 96, Bellettini & Paolini 96, Elliott & Schätzle 96 Asymptotics yield anisotropic mean curvature flow V = γ(γ + γ )κ + γf This model can be discretized using a standard explicit in time finite element method.
13 Anisotropic Allen-Cahn with Forcing γ = 1 γ = cos 4θ 1 t= t= f =
14 Diffusion Induced grain boundary motion (DIGM) When a film of Fe is immersed in a vapour containing Zn atoms the Zn diffuses into the Fe along the grain boundaries causing them to move. VAPOUR z y x A A+V B VAPOUR
15 Two dimensional reductions: Thick film VAPOUR Diffusion Induced grain boundary motion (DIGM) Assume y-independence A A+V B The grain boundaries span Ω VAPOUR Two dimensional reductions: Thin film C+V A A+V C A+V B Assume z-independence Rapid diffusion of solute into the film.
16 Thick film uni-directional motion Cahn, Fife & Penrose 97 proposed a phase field model for DIGM using the energy functional 1 F = 2ε c2 + ε ϕ W (ϕ) + ϕc2 ( ) ε c denotes the concentration of solute atoms in the material. The first term takes into account the chemical energy The second and third represent the interfacial energy The last term approximates the elastic energy caused by a difference in stress across the interface.
17 Thick film uni-directional motion Considering the gradient flow of the energy yields the Allen-Cahn equation with forcing ( ) 1 ε(ϕ t, η ϕ) + ε( ϕ, η ϕ) ε ϕ c2, η ϕ and a diffusion equation for c εc t = (D(ϕ) v). η K Here D(ϕ) = 1 ϕ 2 and v = c + εf ϕ. Deckelnick & Elliott 01 Existence with v = c. Deckelnick, Elliott & Styles 01 Finite element approximation, stability and convergence (v = c and the forcing term c 2 replaced by c). Fife, Cahn & Elliott 01 Formal asymptotics yield V = κ + c 2, c ss = V c
18 Thin film bi-directional motion C.M. Elliott & V. Styles Using the ideas of Fife & Wang 02, yields the following model for thin film bi-directional DIGM V = κ + [c] Γ, { 1 (x, t) s.t. x Γ( t) for some t t c(x, t) = 0 otherwise [c] Γ denotes the jump in concentration across Γ. We derived corresponding phase field and level set models with forcing terms that mimic the jump in concentration across the interfacial region. For the the phase field model we set F[c, ϕ] = c(x + δ + ϕ n, t) c(x δ ϕ n, t) x s.t. ϕ < 1. Garcke & Styles 04 Multi-order parameter phase field system for bi-directional DIGM in the presence of triple junctions.
19 Thin film bi-directional motion bi-directional motion bi-directional motion with triple junctions
20 DIGM with elasticity H. Garcke R. Nürnberg & V. Styles We study a phase field model for DIGM based on the free energy 1 F = + ε 2ε 2c2 2 ϕ ε W (ϕ)+ 1 (E(u) Ē : C(E(u) Ē)) 2 The last term models the elastic energy with u displacement vector E(u) = 1 2 ( u + ( u)t ) strain tensor (linear elasticity) Ē(x, ϕ, c) = Ē 1 (x)ϕ + Ē 2 (x)c + Ē 3 (x) stress free strain C = C (1 ϕ)c+ elasticity tensor (C ± constant tensors) C C +
21 DIGM with elasticity Considering the gradient flow of the free energy gives ( ) ϕ ε t, η ϕ ε ( ϕ, (η ϕ)) + 1 (ϕ, η ϕ) ε 1 2 (E(u) Ē : C ϕ (E(u) Ē)(η ϕ)) +(E(u) Ē : CĒ ϕ (η ϕ)) η K (E(u) Ē : CE(ξ)) = ( C(E(u) Ē), E(ξ) ) ξ (H 1 (Ω)) 2 The concentration c of solute satisfies c(x, t) = { 1 (x, t) s.t. ϕ(x, t) < 1 for some t t 0 otherwise We consider a finite element approximation of this model.
22 The following computations model DIGM in systems with cubic anisotropy. We consider two grains separated by a grain boundary. The two crystal lattices differ by a rotation. We set Ē = c(x) and C = C (1 ϕ)c+. We take C R to be a given elasticity tensor with cubic anisotropy with respect to a given reference orientation of the lattice. C ± are computed by setting for i, j, k, l = 1, 2 where C ± ijkl = 2 p,q,r,s=1 d ± ip d± jq d± kr d± ls CR pqrs ( ) cos θ sin θ D = sin θ cos θ This models a tilt boundary with rotation axis parallel to the x 2 -axis.
23 DIGM with elasticity θ = π 8 θ = 0 θ = π 8 θ = π 6
24 DIGM with elasticity t = 0, 4, 8, θ = π 8 θ = 0 θ = π 6 θ = 0
25 The Cahn-Hilliard equation The original motivation lay in the modelling of phase separation and spinodal decomposition in binary alloys. Order parameter u measuring the amount of one component satisfies the fourth order equation u t = ( ε u + 1 ε W (u)). Second order splitting with double obstacle potential, Blowey, Elliott 94 u t = w, ε( u, η u) (w + 1 u, η u) η K. ε chemical potential w is the functional derivative of F(u) := (ε u ε W (u)). Ω
26 Degenerate Cahn-Hilliard Introducing a concentration u dependent mobility M(u) in the diffusion equation yields:- εu t = (M(u) w), w = 4 π ( ε u + 1 ε W (u)) Taking M(u) := π 4 (1 u2 ) yields a degenerate system. Cahn, Elliott & Novick-Cohen 94 Formal asymptotics yield a geometric evolution law V = κ ss (surface diffusion), w approximates κ. Elliott & Garcke 96 Existence Barrett, Blowey & Garcke Finite element approximation, stability and convergence.
27 Void Electro Migration J.W. Barrett, R. Nürnberg & V.Styles The motion of voids in interconnect lines on microelectronic circuits is modelled by V = (κ + αϕ) ss ϕ is the electric potential associated with an applied electric field α > 0 is related to various physical constants. ϕ n = 0 ϕ n = 1 ϕ = 0 ϕ Γ = 0 ϕ n = 1 ϕ n = 0
28 A sharp diffuse interface approximation of this model is Void Electro Migration εu t = (M(u) (w + αϕ)) M(u) = 1 u 2 ε( u, χ u) (w + 1 u, χ u) ε χ K (c(u) ϕ) = 0 c(u) = 1 (1 u) 2 Barrett, Nürnberg & Styles 02 F.E.A., stability and convergence α = 0 Cahn-Hilliard
29 Void Electro Migration α = 40π α = 120π
30 Epitaxial growth Modern technology of growing single crystals that inherit atomic structure from substrates. It produces almost defect-free high quality materials. Consider the simple geometry of a single island Ω + on a terrace Ω. Set Γ t to be the boundary of the island. Ω Ω +
31 Epitaxial growth The following Burton-Cabrera-Frank type model was studied by Bansch, Hausser, Lakkis, Li & Voigt, 04 The adatom diffusion on the terrace and the island is described by t ρ± D ρ ± = F τ ρ ± ρ (ρ + ) denotes the adatom density on the terrace (island) D diffusion constant, F deposition flux rate, On the island boundary Γ(t) we have 1 τ desorbtion rate. q := D ρ n + V ρ = k (ρ κ) flux from terrace q + := D ρ + n V ρ + = k + (ρ + κ) flux from island k k + > 0 model the Ehrlich-Schwoebel effect. On Γ(t) we have V = κ ss + k + (ρ + κ) + k (ρ κ).
32 Epitaxial growth We study a sharp diffuse interface model of this problem εu t = (M(u) w) + 1 ε M(u)(k+ (ρ + w) + k (ρ w)) ε( u, χ u) (w + 1 u, χ u) ε χ K t t ( H (u)ρ ) (H (u) ρ ) + 1 ε M(u)k (ρ w) = H (u)(f τ ρ ) ( H + (u)ρ +) (H + (u) ρ + ) + 1 ε M(u)k+ (ρ + w) = H + (u)(f τ ρ + ) Here M(u) = 1 u 2, H + (u) = 1 2 (1 + u), H (u) = 1 (1 u). 2 We consider a finite element approximation to the weak formulation of this model.
33 Growth of a single circular island with no desorbtion. t=0, Ω 0 Ω
34 adatom density on the island ρ + adatom density on the terrace ρ
35 Conclusions We approximate interface evolution using phase field models with a double obstacle potential. The method is widely applicable and can be used for second and fourth order problems. It is frequently physically based as the gradient flow of a free energy functional. It can naturally deal with topological changes. Alternative methods used to study interface evolution are direct methods or level set methods.
Modelling 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 informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationPhase-field modeling of step dynamics. University of California, Irvine, CA Caesar Research Center, Friedensplatz 16, 53111, Bonn, Germany.
Phase-field modeling of step dynamics ABSTRACT J.S. Lowengrub 1, Zhengzheng Hu 1, S.M. Wise 1, J.S. Kim 1, A. Voigt 2 1 Dept. Mathematics, University of California, Irvine, CA 92697 2 The Crystal Growth
More informationUniversität Regensburg Mathematik
Universität Regensburg Mathematik On the stable discretization of strongly anisotropic phase field models with applications to crystal growth John W. Barrett, Harald Garcke and Robert Nürnberg Preprint
More informationPhase-field modeling of nanoscale island dynamics
Title of Publication Edited by TMS (The Minerals, Metals & Materials Society), Year Phase-field modeling of nanoscale island dynamics Zhengzheng Hu, Shuwang Li, John Lowengrub, Steven Wise 2, Axel Voigt
More informationLECTURE 4: GEOMETRIC PROBLEMS
LECTURE 4: GEOMETRIC PROBLEMS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and Related
More informationwith deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics
12-3-2009 Modeling and Asymptotics for a general non-homogeneous Cahn-Hilliard equation with deterministic and noise terms D.C. Antonopoulou (Joint with G. Karali and G. Kossioris) Department of Applied
More informationQuantum Dots and Dislocations: Dynamics of Materials Defects
Quantum Dots and Dislocations: Dynamics of Materials Defects with N. Fusco, G. Leoni, M. Morini, A. Pratelli, B. Zwicknagl Irene Fonseca Department of Mathematical Sciences Center for Nonlinear Analysis
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 informationWeak solutions for the Cahn-Hilliard equation with degenerate mobility
Archive for Rational Mechanics and Analysis manuscript No. (will be inserted by the editor) Shibin Dai Qiang Du Weak solutions for the Cahn-Hilliard equation with degenerate mobility Abstract In this paper,
More informationSolid-State Dewetting: From Flat to Curved Substrates
Solid-State Dewetting: From Flat to Curved Substrates Yan Wang with Prof. Weizhu Bao, Prof. David J. Srolovitz, Prof. Wei Jiang and Dr. Quan Zhao IMS, 2 May, 2018 Wang Yan (CSRC) Solid-State Dewetting
More informationMultiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models
Multiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models Russel Caflisch 1, Mark Gyure 2, Bo Li 4, Stan Osher 1, Christian Ratsch 1,2, David Shao 1 and Dimitri Vvedensky 3 1 UCLA,
More informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationLecture 1 Modeling and simulation for the growth of thin films
Lecture 1 Modeling and simulation for the growth of thin films Russel Caflisch Mathematics Department Materials Science and Engineering Department UCLA & IPAM www.math.ucla.edu/~material 1 Outline Epitaxial
More informationA finite element level set method for anisotropic mean curvature flow with space dependent weight
A finite element level set method for anisotropic mean curvature flow with space dependent weight Klaus Deckelnick and Gerhard Dziuk Centre for Mathematical Analysis and Its Applications, School of Mathematical
More informationCoarsening Dynamics for the Cahn-Hilliard Equation. Jaudat Rashed
Coarsening Dynamics for the Cahn-Hilliard Equation Jaudat Rashed Coarsening Dynamics for the Cahn-Hilliard Equation Research Thesis Submitted in partial fulfillment of the requirements for the degree of
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 informationarxiv: v1 [math.oc] 9 Dec 2013
arxiv:1312.2356v1 [math.oc] 9 Dec 2013 Multi-material phase field approach to structural topology optimization Luise Blank, M. Hassan Farshbaf-Shaker, Harald Garcke, Christoph Rupprecht, Vanessa Styles
More informationContinuum Methods 1. John Lowengrub Dept Math UC Irvine
Continuum Methods 1 John Lowengrub Dept Math UC Irvine Membranes, Cells, Tissues complex micro-structured soft matter Micro scale Cell: ~10 micron Sub cell (genes, large proteins): nanometer Macro scale
More informationSolid-State Dewetting: Modeling & Numerics
Solid-State Dewetting: Modeling & Numerics Wei Jiang (ö ) Collaborators: Weizhu Bao (NUS, Singapore), Yan Wang (CSRC, Beijing), Quan Zhao (NUS, Singapore), Tiezheng Qian (HKUST, HongKong), David J. Srolovitz
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationStability Analysis of Stationary Solutions for the Cahn Hilliard Equation
Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Peter Howard, Texas A&M University University of Louisville, Oct. 19, 2007 References d = 1: Commun. Math. Phys. 269 (2007) 765
More informationCOMPARISON OF DIFFERENT NUMERICAL SCHEMES FOR THE CAHN HILLIARD EQUATION
J KSIAM Vol17, No3, 197 207, 2013 http://dxdoiorg/1012941/jksiam201317197 COMPARISON OF DIFFERENT NUMERICAL SCHEMES FOR THE CAHN HILLIARD EQUATION SEUNGGYU LEE 1, CHAEYOUNG LEE 1, HYUN GEUN LEE 2, AND
More informationStep Bunching in Epitaxial Growth with Elasticity Effects
Step Bunching in Epitaxial Growth with Elasticity Effects Tao Luo Department of Mathematics The Hong Kong University of Science and Technology joint work with Yang Xiang, Aaron Yip 05 Jan 2017 Tao Luo
More informationPHASE-FIELD APPROACH TO CRYSTAL GROWTH
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2010 Czech Technical University in Prague, August 30 - September 4, 2010 pp. 18 25 PHASE-FIELD APPROACH TO CRYSTAL GROWTH HUNG HOANG DIEU
More informationWell-posedness and asymptotic analysis for a Penrose-Fife type phase-field system
0-0 Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system Buona positura e analisi asintotica per un sistema di phase field di tipo Penrose-Fife Salò, 3-5 Luglio 2003 Riccarda
More informationNonsmooth Schur Newton Methods and Applications
Nonsmooth Schur Newton Methods and Applications C. Gräser, R. Kornhuber, and U. Sack IMA Annual Program Year Workshop Numerical Solutions of Partial Differential Equations: Fast Solution Techniques November
More informationModelling and Numerical Treatment of Diffuse-Interface Models with Applications in Epitaxial Growth
Modelling and Numerical Treatment of Diffuse-Interface Models with Applications in Epitaxial Growth Dissertation zur Erlangung des Doktorgrades (Dr. rer. nat.) der Mathematisch-Naturwissenschaftlichen
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationApplication of Elastic Flows to Cell Motility. Björn Stinner
Mathematics Institute Centre for Scientific Computing Application of Elastic Flows to Cell Motility Björn Stinner (joint work with C Venkataraman and C Elliott) FBP - Theory and Applications 2012 Cell
More informationFrom nonlocal to local Cahn-Hilliard equation. Stefano Melchionna Helene Ranetbauer Lara Trussardi. Uni Wien (Austria) September 18, 2018
From nonlocal to local Cahn-Hilliard equation Stefano Melchionna Helene Ranetbauer Lara Trussardi Uni Wien (Austria) September 18, 2018 SFB P D ME S. Melchionna, H. Ranetbauer, L.Trussardi From nonlocal
More informationSingular Diffusion Equations With Nonuniform Driving Force. Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009
Singular Diffusion Equations With Nonuniform Driving Force Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009 1 Contents 0. Introduction 1. Typical Problems 2. Variational Characterization
More informationOn the Parametric Finite Element Approximation of Evolving Hypersurfaces in R 3
On the Parametric Finite Element Approximation of Evolving Hypersurfaces in R 3 John W. Barrett Harald Garcke Robert Nürnberg Abstract We present a variational formulation of motion by minus the Laplacian
More informationTHE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS
THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS Alain Miranville Université de Poitiers, France Collaborators : L. Cherfils, G. Gilardi, G.R. Goldstein, G. Schimperna,
More informationSurface stress-induced self-assembly
Surface stress-induced self-assembly Zhigang Suo Work with Wei Lu Solid surface Liquid surface Stretched solid surface Energy per surface atom depends on stretch Stretched liquid surface Energy per surface
More informationFINITE ELEMENT METHOD FOR CONSERVED PHASE FIELD MODELS: SOLID STATE PHASE TRANSFORMATIONS MOHSEN ASLE ZAEEM
FINITE ELEMENT METHOD FOR CONSERVED PHASE FIELD MODELS: SOLID STATE PHASE TRANSFORMATIONS By MOHSEN ASLE ZAEEM A dissertation submitted in partial fulfillment of the requirements for the degree of DOCTOR
More informationModeling and Simulation for Solid-State Dewetting Problems
Modeling and Simulation for Solid-State Dewetting Problems Weizhu Bao Department of Mathematics National University of Singapore Email: matbaowz@nus.edu.sg URL: http://www.math.nus.edu.sg/~bao Collaborators:
More informationAn unfortunate misprint
An unfortunate misprint Robert L. Jerrard Department of Mathematics University of Toronto March 22, 21, BIRS Robert L. Jerrard (Toronto ) An unfortunate misprint March 22, 21, BIRS 1 / 12 Consider the
More informationLevel-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation
Level-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation Bo Li Dept. of Math & NSF Center for Theoretical Biological Physics UC San Diego, USA Collaborators: Li-Tien Cheng
More informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More informationPhase Field Formulation for Microstructure Evolution in Oxide Ceramics
Phase Field Formulation for Microstructure Evolution in Oxide Ceramics Corning Problem Presenter: Leslie Button (Corning) Problem Participants: Daniel Anderson 1, Chris Breward 2, Yanping Chen 3, Brendan
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationFission of a longitudinal strain solitary wave in a delaminated bar
Fission of a longitudinal strain solitary wave in a delaminated bar Karima Khusnutdinova Department of Mathematical Sciences, Loughborough University, UK K.Khusnutdinova@lboro.ac.uk and G.V. Dreiden, A.M.
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationSupplementary Material
1 2 3 Topological defects in confined populations of spindle-shaped cells by G. Duclos et al. Supplementary Material 4 5 6 7 8 9 10 11 12 13 Supplementary Note 1: Characteristic time associated with the
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
More informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
More informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
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 informationTHE TOPOLOGICAL DESIGN OF MATERIALS WITH SPECIFIED THERMAL EXPANSION USING A LEVEL SET-BASED PARAMETERIZATION METHOD
11th. World Congress on Computational Mechanics (WCCM XI) 5th. European Conference on Computational Mechanics (ECCM V) 6th. European Conference on Computational Fluid Dynamics (ECFD VI) July 20-25, 2014,
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationOptimal controls for gradient systems associated with grain boundary motions
Optimal controls for gradient systems associated with grain boundary motions Speaker: Shirakawa, Ken (Chiba Univ., Japan) Based on jointworks with: Yamazaki, Noriaki (Kanagawa Univ., Japan) Kenmochi, Nobuyuki
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 informationUniversität Regensburg Mathematik
Universität Regensburg Mathematik Parametric approximation of Willmore flow and related geometric evolution equations John W. Barrett, Harald Garcke and Robert Nürnberg Preprint Nr. 22/2007 Parametric
More informationGiordano Tierra Chica
NUMERICAL METHODS FOR SOLVING THE CAHN-HILLIARD EQUATION AND ITS APPLICABILITY TO MIXTURES OF NEMATIC-ISOTROPIC FLOWS WITH ANCHORING EFFECTS Giordano Tierra Chica Department of Mathematics Temple University
More informationUNCONDITIONALLY ENERGY STABLE SCHEME FOR CAHN-MORRAL EQUATION 3
UNCONDITIONALLY ENERGY STABLE TIME STEPPING SCHEME FOR CAHN-MORRAL EQUATION: APPLICATION TO MULTI-COMPONENT SPINODAL DECOMPOSITION AND OPTIMAL SPACE TILING R. TAVAKOLI Abstract. An unconditionally energy
More informationDynamics and stochastic limit for phase transition problems. problems with noise. Dimitra Antonopoulou IACM-FORTH
Dynamics and stochastic limit for phase transition problems with noise IACM-FORTH International Conference on Applied Mathematics, September 16-20, 2013, ACMAC, Heraklion, Greece Table of contents 1 1.1
More informationEVOLUTION OF ELASTIC THIN FILMS WITH CURVATURE REGULARIZATION VIA MINIMIZING MOVEMENTS
EVOLUTION OF ELASTIC THIN FILMS WITH CURVATURE REGULARIZATION VIA MINIMIZING MOVEMENTS PAOLO PIOVANO Abstract. The evolution equation, with curvature regularization, that models the motion of a two-dimensional
More informationA Differential Variational Inequality Approach for the Simulation of Heterogeneous Materials
A Differential Variational Inequality Approach for the Simulation of Heterogeneous Materials Lei Wang 1, Jungho Lee 1, Mihai Anitescu 1, Anter El Azab 2, Lois Curfman McInnes 1, Todd Munson 1, Barry Smith
More informationAnisotropic multi-phase-field model: Interfaces and junctions
First published in: PHYSICAL REVIEW E VOLUME 57, UMBER 3 MARCH 1998 Anisotropic multi-phase-field model: Interfaces and junctions B. estler 1 and A. A. Wheeler 2 1 Foundry Institute, University of Aachen,
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 informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationSurfaces in spheres: biharmonic and CMC
Surfaces in spheres: biharmonic and CMC E. Loubeau (joint with C. Oniciuc) Université de Bretagne Occidentale, France May 2014 Biharmonic maps Definition Bienergy functional at φ : (M, g) (N, h) E 2 (φ)
More informationKinetic cross-coupling between non conserved and conserved fields in phase field models
Kinetic cross-coupling between non conserved and conserved fields in phase field models Guillaume Boussinot 1,2,3, Efim A. Brener 1 1 P. Gruenberg Institut, Forschungszentrum Juelich, Germany 2 Max-Planck-Institut
More informationResidual Stress analysis
Residual Stress analysis Informations: strains Macro elastic strain tensor (I kind) Crystal anisotropic strains (II kind) Fe Cu C Macro and micro stresses Applied macro stresses Residual Stress/Strain
More informationA New Approach for the Numerical Solution of Diffusion Equations with Variable and Degenerate Mobility
A New Approach for the Numerical Solution of Diffusion Equations with Variable and Degenerate Mobility Hector D. Ceniceros Department of Mathematics, University of California Santa Barbara, CA 93106 Carlos
More informationCS 468. Differential Geometry for Computer Science. Lecture 17 Surface Deformation
CS 468 Differential Geometry for Computer Science Lecture 17 Surface Deformation Outline Fundamental theorem of surface geometry. Some terminology: embeddings, isometries, deformations. Curvature flows
More informationA generalized MBO diffusion generated motion for constrained harmonic maps
A generalized MBO diffusion generated motion for constrained harmonic maps Dong Wang Department of Mathematics, University of Utah Joint work with Braxton Osting (U. Utah) Workshop on Modeling and Simulation
More informationThe Phase Field Method
The Phase Field Method To simulate microstructural evolutions in materials Nele Moelans Group meeting 8 June 2004 Outline Introduction Phase Field equations Phase Field simulations Grain growth Diffusion
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 informationCoupled second order singular perturbations for phase transitions
Coupled second order singular perturbations for phase transitions CMU 06/09/11 Ana Cristina Barroso, Margarida Baía, Milena Chermisi, JM Introduction Let Ω R d with Lipschitz boundary ( container ) and
More informationRisk Averse Shape Optimization
Risk Averse Shape Optimization Sergio Conti 2 Martin Pach 1 Martin Rumpf 2 Rüdiger Schultz 1 1 Department of Mathematics University Duisburg-Essen 2 Rheinische Friedrich-Wilhelms-Universität Bonn Workshop
More informationSurface Energy, Surface Tension & Shape of Crystals
Surface Energy, Surface Tension & Shape of Crystals Shape of Crystals Let us start with a few observations: Crystals (which are well grown ) have facets Under certain conditions of growth we may observe
More informationSurface phase separation and flow in a simple model of drops and vesicles
Surface phase separation and flow in a simple model of drops and vesicles Tutorial Lecture 4 John Lowengrub Department of Mathematics University of California at Irvine Joint with J.-J. Xu (UCI), S. Li
More informationVariational analysis of a striped-pattern model
Variational analysis of a striped-pattern model Mark Peletier TU Eindhoven Show three things: - variational analysis provides additional tools (for singular limit analysis) - show the tools on an example
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationCML Lecture Series Primer on Homogenization
C ecture Series Primer on Homogenization Bing C. Chen The University of ichigan, Department of echanical Engineering Computational echanics aboratory utline odel problem: asymptotic analysis for a uni-directional
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationUniversität Regensburg Mathematik
Universität Regensburg Mathematik Optimal control of Allen-Cahn systems Luise Blank, M. Hassan Farshbaf-Shaker, Claudia Hecht Josef Michl and Christoph Rupprecht Preprint Nr. 21/2013 Optimal control of
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationHARALD GARCKE AND AMY NOVICK{COHEN modied chemical potentials which are prescribed as functional derivatives of the free energy functional E(u) = Z (u
A SINGULAR LIMIT FOR A SYSTEM OF DEGENERATE CAHN{HILLIARD EQUATIONS HARALD GARCKE AND AMY NOVICK{COHEN Abstract. A singular limit is considered for a system of Cahn-Hilliard equations with a degenerate
More informationSurface finite elements for parabolic equations
1 Surface finite elements for parabolic equations G. Dziuk Abteilung für Angewandte Mathematik University of Freiburg Hermann-Herder-Straße 1 D 7914 Freiburg i. Br. Germany C. M. Elliott Department of
More informationsystem CWI, Amsterdam May 21, 2008 Dynamic Analysis Seminar Vrije Universiteit
CWI, Amsterdam heijster@cwi.nl May 21, 2008 Dynamic Analysis Seminar Vrije Universiteit Joint work: A. Doelman (CWI/UvA), T.J. Kaper (BU), K. Promislow (MSU) Outline 1 2 3 4 Outline 1 2 3 4 Paradigm U
More informationA Post-treatement of the homogenization method in shape optimization
8 th World Congress on Structural and Multidisciplinary Optimization June 1-5, 2009, Lisbon, Portugal A Post-treatement of the homogenization method in shape optimization Olivier Pantz, Karim Trabelsi
More informationBlack Hole Universe with Rotation Chan Park KAIST
Black Hole Universe with Rotation 2016.01.02. Chan Park KAIST Motivation FLRW cosmology Highly symmetric assumption : spatial homogeneity and isotropy Metric ds 2 = dt 2 + a 2 t a t : scale factor Friedmann
More informationHomogenization of Lattice Systems
ANDREA BRAIDES Homogenization of Lattice Systems Ginzburg-Landau equations, Dislocations and Homogenization May 23, 2011, Ile de Ré From discrete to continuous energies Discrete system: with discrete variables
More informationRigorous derivation of a sixth order thin film equation
Rigorous derivation of a sixth order thin film equation Boris Muha Department of Mathematics, Faculty of Science, University of Zagreb Waves in Flows, Prague, August 2018 Joint work with M. Bukal, University
More informationModeling and computational studies of surface shape dynamics of ultrathin single-crystal films with anisotropic surface energy, part I
Modeling and computational studies of surface shape dynamics of ultrathin single-crystal films with anisotropic surface energy, part I Mikhail Khenner Department of Mathematics, Western Kentucky University
More informationA MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION
1 A MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique The most recent results were obtained in collaboration with F. de Gournay, F. Jouve, O. Pantz, A.-M. Toader.
More informationModulation-Feature based Textured Image Segmentation Using Curve Evolution
Modulation-Feature based Textured Image Segmentation Using Curve Evolution Iasonas Kokkinos, Giorgos Evangelopoulos and Petros Maragos Computer Vision, Speech Communication and Signal Processing Group
More informationGeometric Evolution Laws for Thin Crystalline Films: Modeling and Numerics
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 6, No. 3, pp. 433-482 Commun. Comput. Phys. September 2009 REVIEW ARTICLE Geometric Evolution Laws for Thin Crystalline Films: Modeling and Numerics Bo Li 1,
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 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 informationStability of Travelling Wave Solutions for Coupled Surface and Grain Boundary Motion
Stability of Travelling Wave Solutions for Coupled Surface and Grain Boundary Motion Margaret Beck Zhenguo Pan Brian Wetton May 14, 2010 Abstract We investigate the spectral stability of the travelling
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Engineering Acoustics Session 1aEA: Thermoacoustics I 1aEA7. On discontinuity
More informationTHE EXISTENCE OF GLOBAL ATTRACTOR FOR A SIXTH-ORDER PHASE-FIELD EQUATION IN H k SPACE
U.P.B. Sci. Bull., Series A, Vol. 8, Iss. 2, 218 ISSN 1223-727 THE EXISTENCE OF GLOBAL ATTRACTOR FOR A SIXTH-ORDER PHASE-FIELD EQUATION IN H k SPACE Xi Bao 1, Ning Duan 2, Xiaopeng Zhao 3 In this paper,
More informationGeneralized Phase Field Models with Anisotropy and Non-Local Potentials
Generalized Phase Field Models with Anisotropy and Non-Local Potentials Gunduz Caginalp University of Pittsburgh caginalp@pitt.edu December 6, 2011 Gunduz Caginalp (Institute) Generalized Phase Field December
More information