Construction and sharp consistency estimates for atomistic/continuum coupling methods: a 2D model problem
|
|
- Claribel Robertson
- 5 years ago
- Views:
Transcription
1 Construction and sharp consistency estimates for atomistic/continuum coupling methods: a 2D model problem Lei Zhang with: Christoph Ortner Mathematical Institute University of Oxford
2 Point Defects in Crystal [Wikipedia]
3 determines their strength. You CT simulations using the fixed-parameter SW potential, will discuss the crystallographic two different TB schemes [22,23], and our hybrid aspectsscheme. of plastic deformation in Atomistic Com direct diagonalization would be approximately an or more detail in your A4 Material To be able to perform several fully quantum-mechanical Failure course. of magnitude slower. Improving on the accuracy of simulations we used a 64 Si atom cubic cell. As the cell is classical potentials is crucial to simulate this exten Figure 2. the exact TB forces to be used so small, here we compute!"#a $b a) Time: 2.98 fitpsby direct diagonalization b) Time: ps defect. LOTF simulations on this system reveal in the LOTF of the 3.02 periodic topology for the free energy minimum at 900 K wh system Hamiltonian. We note that the large difference!, differs from that obtained in the same conditions (and between the results obtained by fitting the scheme on the 0 Kunder [26]) from the SW model, with a square of ato Single crystals usually deform very easily %30 4 two different quantum models islowdue to the different values of stress. Since they have a nearly adjacent to an antiphase defect (red atom formed predictions of these models in accurate quantitative ideal structure, once thecomdeformation begins at Fig. 4), and a different defect migration pathway. Wh one point, it can continue the crystal. (111) putations, reproduced by the present method. This throughout emthe TB stress levels for plastic deformation are model used [23] may still not capture ev phasizes scheme can at best be (112)the fact that the present Typical in the range between 1 and 10 MPa,relevant much too feature of the phenomenon under study, th expected to reproduce the results of the QM model that it low for most engineering purposes. results indicate the need to enforce electronic-struct is given [24], but can in no way improve its accuracy. Most engineering alloys are precision. used in Indeed, in dislocation dynamics the inaccur To further elucidate how the scheme adaptsform to the local polycrystalline (Figure 3). This means that canonical phase space sampling provided by class 930#0 environment, we report in Fig. 3 the the up of a great eachtime pieceevolution of metal isofmade potentials is regarded as the main source of discrepa 205 Korsunsky; lecture notes Figure 3. number of single crystals, or grains, each having Csanyi etal;jphys.:condens.matter(2005) ARTICLEparameters IN PRESS Cijk at a dimerized Si[100] surface at room between theory and experiment [26]. al. / J. Mech. Phys. Solids 55 (2007) Dobson & Luskin (grain boundary) c) Time: 3.15Inps Time: 3.47 (crackremain tip inps temperature. the bulk these angled) parameters atsilicon) 4 M=Numbers of atom (edge dislocation) all times close to their!1=3 equilibrium value (left). However, on the reconstructed surface, the equilibrium Current computationa angle on the lower side of the buckled dimers (inset) is Empirical potent lowered to almost 90". The corresponding parameters learn this by switching to a value close to zero, and Semi-empirical p flip back and forth between zero and!1=3 as the buckling direction varies with time (right). Quantum mecha Moving to a problem where a fully quantum approach Nonconvex potential would be practically unfeasible, we simulate the gliding motion at 900 K of an opposing pair of 30" partial Can be modeled by c dislocations in Si, using a 4536 atom unit cell. We flag Csanyi et al; PRL 93; 2004 for QM Gavini; treatment all atoms within 7.0 A of Miller,Tadmor the dangling(2003) PhD Thesis ctron-density around the vacancy on (1 0 0) plane; (b) contours of groundfig. 4 (color). Equilibrium structure of the Use Si 30"ato par AtC method: (dislocation in silicon) Figurebonds 6. Four snapshots of theatoms) opening (111)[1 1 0] crack systembdry) simulated using LOTF. Atoms in (undercoordinated formed during dislocancy on (1 1 1) plane. (crack & grain (vacancy in aluminium) dislocation kink from a 900 K hybrid simulation (only tion motion,quantum corresponding to two disjoint QM of red were treated mechanically using the abregions initio SIESTA package.left The quantum region glide plane atoms are shown). One undercoordinated atomby ( freedom elsewhere ariation of the vacancy formation with sample size.the diffusing kink of each atoms,energy which follow Defects in Crystals
4 Scales in Materials Modelling Current computational limits: Empirical potentials (EAM): atoms Semi-empirical potentials (TB): 10 6 atoms Almost ab initio (KS-DFT): 10 4 atoms Schrödinger Equation: a few electrons Density of atoms in typical materials: atoms/mm 3 Goal: QM MM FEM Rountree et al, Annu.Rev.Mater.Res.2002
5 u Application: Crack at Grain Boundary Applications of Coarse-Graining [Iglesias/Leiva, Acta Mater. 06] complicated defects, defect interaction, crack growth, growth,... 8m atoms, void15000 repatoms Application: Shear in Twinned Aluminium B Langwallner (Oxford University) [Iglesia & Leiva, 2006], [Miller et al, 1998], The Quasicontinuum Method [QC software tutorial, Miller & Tadmor]
6 Atomistic Mechanics (0T statics) Atomistic body: N atoms at positions y = (y n ) N n=1 Rd N Total energy of configuration y: min Ea tot (y) := E a (y) + P a (y) E a = interaction potential, P a = potential of external frcs
7 Atomistic Mechanics (0T statics) Atomistic body: N atoms at positions y = (y n ) N n=1 Rd N Total energy of configuration y: min Ea tot (y) := E a (y) + P a (y) E a = interaction potential, P a = potential of external frcs Crystallisation:[Theil, Theil/Harris,... ]
8 The Quasicontinuum Idea (a) (b) Molecular statics problem: y a argmin E tot a (Y ) Coarse grained problem: y h argmin E tot a (Y h ) where mesh T h resolves the defect, and Y h = Y P 1 (T h ).
9 Atomistic Stored Energy: E a (y) = Cauchy Born Approximation x L Cauchy Born Stored Energy: E c (y) = W ( y) dv, Ω V ( y(x + r) y(x); r R ) where W (F) = V ({Fr; r R}).
10 Atomistic Stored Energy: E a (y) = Cauchy Born Approximation x L Cauchy Born Stored Energy: E c (y) = W ( y) dv, Ω V ( y(x + r) y(x); r R ) where W (F) = V ({Fr; r R}). Theorem: [Similar to result by E/Ming; 2007] Let y a argmin Ea tot be sufficiently smooth globally, then there exists y c argmin Ec tot such that y a y c L 2 C ( 3 y a L y a 2 ) L 4 If there are no defects, then the Cauchy Born model is a highly accurate continuum approximation.
11 Atomistic/Continuum Coupling: First Attempt (a) (b) (c) E a (y h ) E ac (y h ) := ω x V x + W ( y h ) dx x L Ω c a
12 Atomistic/Continuum Coupling: First Attempt Fails the patch test: δe a (y F ) = 0 and δe c (y F ) = 0, but δe qce (y F ) 0! (a) (b) (c) E a (y h ) E ac (y h ) := ω x V x + W ( y h ) dx x L Ω c a Error: x Component x 10 3 Error: y Component
13 Alternative Approaches 1 Energy-based coupling: interface correction 2 Force-based coupling: FeAt: Kohlhoff, Schmauder, Gumbsch (1989, 1991) Dead-load GF removal: Shenoy, Miller, Rodney, Tadmor, Phillips, Ortiz (1999) AtC: Parks, Gunzburger, Fish, Badia, Bochev, Lehoucq, et al. (2007,... ) CADD: Shilkrot & Curtin & Miller (2002,... ) Stress-based coupling: Makridakis/Ortner/Süli (2011)... 3 Blending methods: Belytschko & Xiao (2004) Klein & Zimmerman (2006) Parks, Gunzburger, Fish, Badia, Bochev, Lehoucq, et al. (2008)...
14 A Priori Error Analysis Framework: Let y a argmin Ea tot, y ac argmin Eac tot, then... (standard ideas)... [y a y ac ] L 2 CONSISTENCY STABILITY = δe a(y a ) δe ac (y a ) H 1 inf u L 2 =1 δ 2 E ac (y a )u, u 3 Steps: 1 CONSISTENCY: δe ac (y) δe a (y), u 2 y L a 2 (Ω c) 2 STABILITY: δ 2 E ac (y)u, u C stab u 2 L 2 3 REGULARITY: bounds on 2 y a, note that r a for defects
15 Consistency Numerical Analysis Literature on Energy-Based Coupling: 1D, NN, variational analysis: [Blanc/LeBris/Legoll, 2005] 1D, Error estimates for QCE and QNL: [Dobson/Luskin, 2009], [Ming/Yang, 2010], [Ortner, 2010], [Ortner/Wang, 2011] 1D, Sharp Stability Analysis, Linear Regime: [Dobson/Luskin/Ortner, 2010] 1D, finite range: [Li/Luskin, preprint] 1D, EAM potentials: [Li/Luskin, 2011] 1D, blending methods: [Luskin/van Koten, preprint] 2D, pair interactions, defects: [Ortner/Shapeev, preprint] 2D, first order consistency for general finite-range interactions: [Ortner, preprint] 2D, multi-body interactions: [Ortner/Zhang, preprint] Analysis for Force-based Coupling: Dobson & Luskin (2008); Ming (2009); Dobson & Luskin & Ortner (2009, 2010, 2010); Makridakis & Ortner & Süli (2010, preprint); Dobson & Ortner & Shapeev (preprint), Lu & Ming (manuscript) Analysis for Multi-Lattices: Dobson &Elliott & Luskin & Tadmor (2007); Abdulle & Lin & Shapeev (preprint)
16 Patch Test Consistency implies First-order Consistency? Suppose the interface potentials Ṽx are fitted numerically, s.t. E ac becomes patch test consistent: δe ac (y F ) = 0 F R d d Does this automatically imply that E ac is first-order consistent? Theorem: [First-order Consistency] Suppose δe ac passes the patch test, V finite range multi-body potential + technical conditions + d = 1; or d = 2, Ω a connected; or d = 3, flat interface, homogeneous Ṽx then δeac (y) δe a (y), u h h 2 y L 2 (Ω c) u h L 2 If E ac is also stable then this would imply y a y ac L 2 h 2 y a L 2 (Ω c)
17 Consistency in 1d no ghost forces = consistent in a negative Sobolev space. A/C method E ac has no ghost forces if F ac (y F ; x) = E ac(y) y(x) = 0 F > 0 x L. Theorem: [sharp estimate in 1d] Assume that the A/C energy E ac has no ghost forces, then there exists a constant C depending only on M 2 (µ) and M 3 (µ), such that sup δe ac (y) δe a (y), u C ( y l p (I) + y l p (C) + y 2 u R L l (C)). 2p u l p =1
18 Consistency in 1d no ghost forces = consistent in a negative Sobolev space. 1 A/C stress function δe ac(y), u = x L Σac(y; x)u (x), where F ac(y; x) = Σ ac(y; x + 1) Σ ac(y; x). Similarly Σ a(y; x), Σ c(y; x). 2 2nd order error for Σ c(y; x) (by symmetry, V(b)=V(-b)) Σa(y; x) Σ c(y; x) C ( y l (x+r) + y 2 l (x R)). 3 Σ ac(y; x) = Σ a(y; x) for x A; Σ ac(y; x) = Σ c(y; x) for x C, 4 F ac(y F ; x) = 0 Σ ac(y F ; x) = Σ a(y F ; x) = Σ c(y F ; x), 5 Σ ac and Σ a are Lipschitz continous when V is smooth Σac(y; x) Σ a(y; x) Σac(y; x) Σ ac(y F ; x) + Σa(y F ; x) Σ a(y; x) C F y l (I ± ) Σac(y; x) Σ a(y; x) C y l ((x R) I) 6 use δe ac(y) δe a(y), u = x L (Σac(y; x) Σa(y; x))u (x) to complete the proof.
19 Construction of General A/C Schemes E ac (y h ) = V x + Ṽ x + Vx c x L a x L i x L c Construct Ṽ s.t. δe ac(y F ) = 0 for all F R d d. General Construction: [Shimokawa et al, 2004; E/Lu/Yang, 2006] Ṽ x = V ( g x,r ; r R ) g x,r = s R x C x,r,s g s Find C x,r,s s.t. δe ac (y F ) = 0 geometric conditions only! F Explicit constructions for flat interfaces: [E/Lu/Yang, 2006] Explicit constructions for 2D general interface: [Ortner/Zhang; preprint] In general: compute C x,r,s numerically in preprocessing 2d, NN, multibody potential, triagular lattice
20 Construction of General A/C Schemes General Construction: E ac (y h ) = V x + Ṽ x + Vx c x L a x L i x L c Construct Ṽ s.t. δe ac(y F ) = 0 for all F R d d. Ṽ x = V ( g x,r ; r R ) g x,r = s R x C x,r,s g s Find C x,r,s s.t. δe ac (y F ) = 0 geometric conditions only! F 2. Patch Test Consistency 0 = δe ac(y F ), u = V F,r C x,r,sd su x L r R s R = (C x as,r,sv F,r C x,r,sd r V F,r )u(x) x L r R s R 1. Local Energy Consistency Ṽ (y F ) = V (y F ) r = s R x C x,r,ss. (a) (C x s,r,sv F,r C x,r,sv F,r ) = 0. (b) r R s R Solve (a) + (b) + B.C. in L a and L c to obtain C x,r,s for x L i. unknowns: L I R 2, eqns: 2 L I R.
21 Construction of General A/C Schemes E ac (y h ) = V x + Ṽ x + Vx c x L a x L i x L c Construct Ṽ s.t. δe ac(y F ) = 0 for all F R d d. General Construction: Flat interface Interface with corner C x,r,r for NN interaction, multibody potential, one-sided construction. 1. works for general interface in 2d 2. preprocessing for longer interaction range
22 Consisency of The Schemes E ac (y h ) = V x + Ṽ x + Vx c x L a x L i x L c Construct Ṽ s.t. δe ac (y F ) = 0 for all F R d d. Theorem: [Ortner/Zhang] There exists a constant C that depends only on M 2 (V ) and M 3 (V ) such that δe ac (y) δe a (y) U 1,p C ( 3 y l p (Ω c) + 2 y 2 l 2p (Ω c) + 2 y l p (Ω I )).
23 Stress Function in 2D Define atomistic and continuum stress, δe a (y), u = T T T Σ a (y; T ) : T u δe c (y), u = T Σ c (y; T ) : T u T T index conventions: δe a(y), u 6 = V x,j D j u(x) x L j=1 ( 1 = T 2 T T 6 ) V xt,j,j a j : T u j=1 Since D j u(x) = 1 2 ( T x,j u a j + Tx,j 1 u a j ).
24 Stress Function in 2D Define atomistic and continuum stress, δe a (y), u = T T δe c (y), u = T T T Σ a (y; T ) : T u T Σ c (y; T ) : T u δe a(y), u 6 = V x,j D ju(x) (D ju(x) = D j+3u(x + a j)) x L j=1 1 6 = T V xt 2 T,j,j a j : T u T T j=1 }{{} =Σ 1 a (y;t ) 1 6 1( = T VxT,j,j V xt 2 T 2,j+3,j+3) a j : T u T T j=1 }{{} =Σ 2 a (y;t ) index conventions:
25 Stress Function in 2D Define atomistic and continuum stress, Similarly, δec(y), u δe a (y), u = T T δe c (y), u = T T = T W ( T y) : T u(t ) T T 1 6 = T j V ( T ya) a j 2 T T T j=1 }{{} =Σ 1 c (y;t )= W ( T y) : T u 1 6 1( ) = T VT,j + V TT,j,j aj : T u 2 T 2 T T j=1 }{{} =Σ 2 c (y;t ) T Σ a (y; T ) : T u T Σ c (y; T ) : T u index conventions:
26 Stress Function in 2D Define atomistic and continuum stress, δe a (y), u = T T T Σ a (y; T ) : T u δe c (y), u = T Σ c (y; T ) : T u T T index conventions: stress functions are not unique Σ 1 a(y; T ) W ( T y) = O( D 2 y ) Σ 2 a(y; T ) Σ 2 c(y; T ) = O( D 2 y 2 + D 3 y ) Theorem: [Consistency for CB energy] δea(y) δe c(y) U 1,p C( 3 y l p + 2 y 2 l 2p ).
27 Discrete Divergence Free Tensor Field The stress function is not unque, up to a divergence free tensor field, namely, σ such that σ : T u = 0 T T Characterization of div-free field, [Arnold/Falk, Polthier/Preuß] Lemma: v piecewise constant vector on T, v divergence free, i.e., v T u = 0 T T iff a function w N 1 (T ), such that v = J w, where J is the counter-clockwise rotation by π/2. N 1 (T ) is Crouzeix Raviart finite element space. For basis ζ f N 1(T ) associated with f = T 1 T 2, J ζ f T u = Jn 1 T1 u + Jn 2 T2 u = u 1 u 2 + u 2 u 1 = 0 T T
28 Characterization of stress function If an A/C energy E ac satisfies patch test consistency, 0 = δe ac (y F ), u = T T T Σ ac (y F ; T ) : T u Lemma: a function ψ(f, T ) N 1 (T ) 2, such that Σ ac (y F ; T ) = W (F) + J ψ(f; T ) For general deformation y, deformation gradient average for patch ω f = (T 1 T 2 ), F f (y) = y dx, Corrector function: ˆψ(y; ) = f F ψ( F f (y); m f ) ζf Define the modified stress function, Σ ac (y; T ) := Σ ac (y; T ) J ˆψ(y; T ), for T T. Σ ac(y F ; T ) = W (F) = Σ a(y F ; T ) Σ ac(y; T ) = Σ a(y; T ), T Ω A and Σac(y; T ) = Σ c(y; T ), T Ω C ω f
29 Proof of Consistency δe ac (y) δe a (y), u = T T (Σ ac (y; T ) Σ a (y; T )) : u = T T ( Σ ac (y; T ) Σ a (y; T )) : u 1 T Ω A, Σac(y; T ) = Σ a(y; T ), 2 T Ω C, Σac(y; T ) Σ a(y; T ) = Σ c(y; T ) Σ a(y; T ), 2nd order consistency 3 T Ω I, Let y T = y y(t ), we have, Σac(y; T ) Σ a(y; T ) Σac(y; T ) W ( y(t )) + W ( y(t )) Σ a(y; T ) = Σac(y; T ) Σac(y T ; T ) + Σ a(y T ; T ) Σ a(y; T ) C y(t ) y l C D 2 y δe ac(y) δe a (y) U 1,p C ( 3 y lp (Ω c) + 2 y 2 l 2p (Ω c) + 2 y lp (Ω I )).
30 Numerical Experiment Test Problem: single vacancy in the triangular lattice, Vacancy defect in 269,100 atom cell atomic spacings Eac Ea / Ea E0 QCE B-QCE (DoF) 1 (DoF) 2 GR-AC yac ya 1,2/ ya id 1,2 QCE B-QCE (DoF) 1/2 GR-AC (DoF) yac ya 1, / ya x 1, QCE B-QCE (DoF) 1 GR-AC # DoFs # DoFs 10 3 # DoFs Error in Energy H 1 Error W 1, Error
31 Outlook on A/C Methods Summary Ghost force removal Patch test consistency consistency Construction of practical energy-based a/c methods Sharp consistency error estimates through stress based formulation
32 Outlook on A/C Methods Summary Ghost force removal Patch test consistency consistency Construction of practical energy-based a/c methods Sharp consistency error estimates through stress based formulation Open Problems Consistency in 3D: proof or counterexample Stability A/C coupling at surfaces Implementation, benchmarks, applications
33 Outlook on A/C Methods Summary Ghost force removal Patch test consistency consistency Construction of practical energy-based a/c methods Sharp consistency error estimates through stress based formulation Open Problems Consistency in 3D: proof or counterexample Stability A/C coupling at surfaces Implementation, benchmarks, applications Major Open Problems A/C methods for multi-lattices A/C methods for Coulomb interaction A/C methods for electronic structure models (done only for insulators)
Construction and Analysis of Consistent Atomistic/Continuum Coupling Methods
Construction and Analysis of Consistent Atomistic/Continuum Coupling Methods Lei Zhang Department of Mathematics & Institute of Natural Sciences Shanghai Jiao Tong University, China with Christoph Ortner
More informationarxiv: v1 [math.na] 11 Dec 2011
FORMULATION AND OPTIMIZATION OF THE ENERGY-BASED BLENDED QUASICONTINUUM METHOD arxiv:1112.2377v1 [math.na] 11 Dec 2011 M. LUSKIN, C. ORTNER, AND B. VAN KOTEN Abstract. We formulate an energy-based atomistic-to-continuum
More informationOn atomistic-to-continuum couplings without ghost forces
On atomistic-to-continuum couplings without ghost forces Dimitrios Mitsoudis ACMAC Archimedes Center for Modeling, Analysis & Computation Department of Applied Mathematics, University of Crete & Institute
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 informationarxiv: v1 [math.na] 30 Jun 2014
ATOMISTIC/CONTINUUM BLENDING WITH GHOST FORCE CORRECTION CHRISTOPH ORTNER AND L. ZHANG arxiv:1407.0053v1 [math.na] 30 Jun 2014 Abstract. We combine the ideas of atomistic/continuum energy blending and
More informationOvercoming Temporal and Spatial Multiscale Challenges In Materials Modeling And Computing
Overcoming Temporal and Spatial Multiscale Challenges In Materials Modeling And Computing Mitchell Luskin School of Mathematics University of Minnesota May 16, 2014 A. Binder (UMN), M. Dobson (Mass), W.-K.
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 informationFINITE ELEMENT ANALYSIS OF CAUCHY BORN APPROXIMATIONS TO ATOMISTIC MODELS
FINITE ELEMENT ANALYSIS OF CAUCHY BORN APPROXIMATIONS TO ATOMISTIC MODELS CHARALAMBOS MAKRIDAKIS AND ENDRE SÜLI ABSTRACT This paper is devoted to a new finite element consistency analysis of Cauchy Born
More informationto appear in International Journal for Multiscale Computational Engineering
to appear in International Journal for Multiscale Computational Engineering GOAL-ORIENTED ATOMISTIC-CONTINUUM ADAPTIVITY FOR THE QUASICONTINUUM APPROXIMATION MARCEL ARNDT AND MITCHELL LUSKIN Abstract.
More informationError Control for Molecular Statics Problems
International Journal for Multiscale Computational Engineering, 4(5&6)647 662(26) Error Control for Molecular Statics Problems Serge Prudhomme, Paul T. Bauman & J. Tinsley Oden Institute for Computational
More informationElectronic-structure calculations at macroscopic scales
Electronic-structure calculations at macroscopic scales M. Ortiz California Institute of Technology In collaboration with: K. Bhattacharya, V. Gavini (Caltech), J. Knap (LLNL) BAMC, Bristol, March, 2007
More informationError Control for Molecular Statics Problems
Error Control for Molecular Statics Problems Serge Prudhomme 1, Paul T. Bauman 2, and J. Tinsley Oden 3 Institute for Computational Engineering and Sciences The University of Texas at Austin Austin, Texas
More informationModeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA
HHHHM Modeling Materials Continuum, Atomistic and Multiscale Techniques ELLAD B. TADMOR University of Minnesota, USA RONALD E. MILLER Carleton University, Canada gg CAMBRIDGE ^0 UNIVERSITY PRESS Preface
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 28 Jun 2001
arxiv:cond-mat/665v [cond-mat.mtrl-sci] 28 Jun 2 Matching Conditions in -Continuum Modeling of Materials Weinan E and Zhongyi Huang 2 Department of Mathematics and PACM, Princeton University and School
More informationBridging methods for coupling atomistic and continuum models
Bridging methods for coupling atomistic and continuum models Santiago Badia 12, Pavel Bochev 1, Max Gunzburger 3, Richard Lehoucq 1, and Michael Parks 1 1 Sandia National Laboratories, Computational Mathematics
More information5 Multiscale Modeling and Simulation Methods
5 Multiscale Modeling and Simulation Methods This chapter is dedicated to a discussion of multiscale modeling and simulation methods. These approaches are aimed to provide a seamless bridge between atomistic
More informationarxiv: v3 [math.na] 8 Aug 2011
Consistent Energy-Based Atomistic/Continuum Coupling for Two-Body Potentials in One and Two Dimensions arxiv:1010.051v3 [math.na] 8 Aug 011 Alexander V. Shapeev October 5, 018 Abstract This paper addresses
More informationAb initio Berechungen für Datenbanken
J Ab initio Berechungen für Datenbanken Jörg Neugebauer University of Paderborn Lehrstuhl Computational Materials Science Computational Materials Science Group CMS Group Scaling Problem in Modeling length
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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationc 2010 Society for Industrial and Applied Mathematics
MULTISCALE MODEL. SIMUL. Vol. 8, No., pp. 571 590 c 010 Society for Industrial and Applied Mathematics A QUADRATURE-RULE TYPE APPROXIMATION TO THE QUASI-CONTINUUM METHOD MAX GUNZBURGER AND YANZHI ZHANG
More informationCOMPUTATIONAL MODELING OF SHAPE MEMORY MATERIALS
COMPUTATIONAL MODELING OF SHAPE MEMORY MATERIALS Jan Valdman Institute of Information Theory and Automation, Czech Academy of Sciences (Prague) based on joint works with Martin Kružík and Miroslav Frost
More informationANALYSIS OF A ONE-DIMENSIONAL NONLOCAL QUASICONTINUUM METHOD
ANALYSIS OF A ONE-DIMENSIONAL NONLOCAL QUASICONTINUUM METHOD PINGBING MING AND JERRY ZHIJIAN YANG Abstract. The accuracy of the quasicontinuum method is analyzed using a series of models with increasing
More informationAdaptive modelling in atomistic-to-continuum multiscale methods
Journal of the Serbian Society for Computational Mechanics / Vol. 6 / No. 1, 2012 / pp. 169-198 (UDC: 530.145.6) Adaptive modelling in atomistic-to-continuum multiscale methods E. Marenic 1,2, J. Soric
More informationADAPTIVE ATOMISTIC-CONTINUUM MODELING OF DEFECT INTERACTION WITH THE DEBDM
Journal for Multiscale Computational Engineering, 11 (6): 505 525 (2013) ADAPTIVE ATOMISTIC-CONTINUUM MODELING OF DEFECT INTERACTION WITH THE DEBDM Philip Moseley, 1 Jay Oswald, 2, & Ted Belytschko 1 1
More informationarxiv: v3 [math.na] 2 Jul 2010
ITERATIVE METHODS OR THE ORCE-BASED QUASICONTINUUM APPROXIMATION: ANALYSIS O A 1D MODEL PROBLEM M. DOBSON, M. LUSKIN, AND C. ORTNER arxiv:0910.2013v3 [math.na] 2 Jul 2010 Abstract. orce-based atomistic-continuum
More informationJournal of the Mechanics and Physics of Solids
Journal of the Mechanics and Physics of Solids 59 (2011) 775 786 Contents lists available at ScienceDirect Journal of the Mechanics and Physics of Solids journal homepage: www.elsevier.com/locate/jmps?
More informationFirst-Principles Quantum Simulations of Nanoindentation
Chapter 10 First-Principles Quantum Simulations of Nanoindentation Qing Peng Additional information is available at the end of the chapter http://dx.doi.org/10.5772/48190 1. Introduction Modern nanotechnology
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationA finite deformation membrane based on inter-atomic potentials for single atomic layer films Application to the mechanics of carbon nanotubes
A finite deformation membrane based on inter-atomic potentials for single atomic layer films Application to the mechanics of carbon nanotubes Marino Arroyo and Ted Belytschko Department of Mechanical Engineering
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 informationMultiscale Modeling with Extended Bridging Domain Method
Multiscale Modeling with Extended Bridging Domain Method Mei Xu, Robert Gracie, and Ted Belytschko Northwestern University 2145 Sherman Ave, Evanston, IL, 60208 1 Introduction This Chapter describes the
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationA Coupled Discrete/Continuous Method for Computing Lattices. Application to a Masonry-Like Structure
Author manuscript, published in "International Journal of Solids and Structures 48, 21 (2011) 3091-3098" DOI : 10.1016/j.ijsolstr.2011.07.002 A Coupled Discrete/Continuous Method for Computing Lattices.
More informationSTRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationA New Multigrid Method for Molecular Mechanics
A New Multigrid Method for Molecular Mechanics Pingbing Ming k mpb@lsec.cc.ac.cn http//lsec.cc.ac.cn/ mpb Joint with Jingrun Chen (ICMSEC) &W. E (Princeton University) Supported by 973, NSFC Ecole des
More informationA surface Cauchy Born model for nanoscale materials
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 26; 68:172 195 Published online 9 May 26 in Wiley InterScience (www.interscience.wiley.com). DOI: 1.12/nme.1754 A surface
More informationMulti-paradigm modeling of fracture of a silicon single crystal under mode II shear loading
Journal of Algorithms & Computational Technology Vol. 2 No. 2 203 Multi-paradigm modeling of fracture of a silicon single crystal under mode II shear loading Markus J. Buehler 1, *, Alan Cohen 2, and Dipanjan
More informationc 2012 Society for Industrial and Applied Mathematics
MULTISCALE MODEL. SIMUL. Vol. 10, No. 3, pp. 744 765 c 2012 Society for Industrial and Applied Mathematics THE SPECTRUM OF THE FORCE-BASED QUASICONTINUUM OPERATOR FOR A HOMOGENEOUS PERIODIC CHAIN M. DOBSON,
More informationA Multilattice Quasicontinuum for Phase Transforming Materials: Cascading Cauchy Born Kinematics
Journal of Computer-Aided Materials Design manuscript No. (will be inserted by the editor) A Multilattice Quasicontinuum for Phase Transforming Materials: Cascading Cauchy Born Kinematics M. Dobson R.
More informationLarge-scale real-space electronic structure calculations
Large-scale real-space electronic structure calculations YIP: Quasi-continuum reduction of field theories: A route to seamlessly bridge quantum and atomistic length-scales with continuum Grant no: FA9550-13-1-0113
More informationc 2018 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 40, No. 4, pp. A2087 A2119 c 2018 Society for Industrial and Applied Mathematics A POSTERIORI ERROR ESTIMATION AND ADAPTIVE ALGORITHM FOR ATOMISTIC/CONTINUUM COUPLING IN TWO DIMENSIONS
More informationAugust 2005 June 2010 Assistant Professor
Xiantao Li Curriculum Vita 219C McAllister Building Department of Mathematics Penn State University University Park, PA 16802 (814)8639081 xxl12@psu.edu www.personal.psu.edu/xxl12 Education Ph.D. 2002,
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
More informationAn Efficient Multigrid Method for Molecular Mechanics Modeling in Atomic Solids
Commun. Comput. Phys. doi: 10.4208/cicp.270910.131110a Vol. 10, No. 1, pp. 70-89 July 2011 An Efficient Multigrid Method for Molecular Mechanics Modeling in Atomic Solids Jingrun Chen 1 and Pingbing Ming
More informationMolecular Dynamics Simulation of Fracture of Graphene
Molecular Dynamics Simulation of Fracture of Graphene Dewapriya M. A. N. 1, Rajapakse R. K. N. D. 1,*, Srikantha Phani A. 2 1 School of Engineering Science, Simon Fraser University, Burnaby, BC, Canada
More informationSurface Cauchy-Born analysis of surface stress effects on metallic nanowires
Surface Cauchy-Born analysis of surface stress effects on metallic nanowires Harold S. Park 1, * and Patrick A. Klein 2 1 Department of Civil and Environmental Engineering, Vanderbilt University, Nashville,
More informationThe atomic-scale finite element method
Comput. Methods Appl. Mech. Engrg. 193 (2004) 1849 1864 www.elsevier.com/locate/cma The atomic-scale finite element method B. Liu, Y. Huang a, *, H. Jiang a,s.qu a, K.C. Hwang b a Department of Mechanical
More informationConcurrently coupled atomistic and XFEM models for dislocations and cracks. Robert Gracie and Ted Belytschko, SUMMARY
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng (2008) Published online in Wiley InterScience (www.interscience.wiley.com)..2488 Concurrently coupled atomistic and
More informationAccuracy and transferability of GAP models for tungsten
Accuracy and transferability of GAP models for tungsten Wojciech J. Szlachta Albert P. Bartók Gábor Csányi Engineering Laboratory University of Cambridge 5 November 214 Motivation Number of atoms 1 1 2
More informationKinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys
Kinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys, Scott Dunham Department of Electrical Engineering Multiscale Modeling Hierarchy Configuration energies and transition
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 informationCoupling atomistic and continuum modelling of magnetism
Coupling atomistic and continuum modelling of magnetism M. Poluektov 1,2 G. Kreiss 2 O. Eriksson 3 1 University of Warwick WMG International Institute for Nanocomposites Manufacturing 2 Uppsala University
More informationMULTIRESOLUTION MOLECULAR MECHANICS: DYNAMICS, ADAPTIVITY, AND IMPLEMENTATION. Emre Biyikli. B.S. in Mechanical Engineering, Koc University, 2007
MULTIRESOLUTION MOLECULAR MECHANICS: DYNAMICS, ADAPTIVITY, AND IMPLEMENTATION by Emre Biyikli B.S. in Mechanical Engineering, Koc University, 2007 M.S. in Mechanical Engineering, Koc University, 2009 Submitted
More informationMultiscale Simulations
Multiscale Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia nstitute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu DFT-MD coupling Topics First-Principles Molecular
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
More informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More informationSmall-Scale Effect on the Static Deflection of a Clamped Graphene Sheet
Copyright 05 Tech Science Press CMC, vol.8, no., pp.03-7, 05 Small-Scale Effect on the Static Deflection of a Clamped Graphene Sheet G. Q. Xie, J. P. Wang, Q. L. Zhang Abstract: Small-scale effect on the
More informationConcurrent AtC coupling based on a blend of the continuum stress and the atomistic force
Concurrent AtC coupling based on a blend of the continuum stress and the atomistic force Jacob Fish a, Mohan A. Nuggehally b, Mark S. Shephard c and Catalin R. Picu d Scientific Computation Research Center
More informationElasticity Constants of Clay Minerals Using Molecular Mechanics Simulations
Elasticity Constants of Clay Minerals Using Molecular Mechanics Simulations Jin-ming Xu, Cheng-liang Wu and Da-yong Huang Abstract The purpose of this paper is to obtain the elasticity constants (including
More informationIntroduction to Solving the Time- Dependent Schrödinger Equation. Tom Penfold
Introduction to Solving the Time- Dependent Schrödinger Equation Tom Penfold Outline 1 Introduction to Solving the Time-Dependent Schrödinger Equation What problems are we trying to solve? How can we use
More informationTemperature-related Cauchy Born rule for multiscale modeling of crystalline solids
Computational Materials Science 37 (2006) 374 379 www.elsevier.com/locate/commatsci Temperature-related Cauchy Born rule for multiscale modeling of crystalline solids Shaoping Xiao a, *, Weixuan Yang b
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 informationScanning Tunneling Microscopy. how does STM work? the quantum mechanical picture example of images how can we understand what we see?
Scanning Tunneling Microscopy how does STM work? the quantum mechanical picture example of images how can we understand what we see? Observation of adatom diffusion with a field ion microscope Scanning
More informationEnd forming of thin-walled tubes
Journal of Materials Processing Technology 177 (2006) 183 187 End forming of thin-walled tubes M.L. Alves a, B.P.P. Almeida b, P.A.R. Rosa b, P.A.F. Martins b, a Escola Superior de Tecnologia e Gestão
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationSingular integro-differential equations for a new model of fracture w. curvature-dependent surface tension
Singular integro-differential equations for a new model of fracture with a curvature-dependent surface tension Department of Mathematics Kansas State University January 15, 2015 Work supported by Simons
More informationCleavage Planes of Icosahedral Quasicrystals: A Molecular Dynamics Study
Cleavage Planes of Icosahedral Quasicrystals: A Molecular Dynamics Study Frohmut Rösch 1, Christoph Rudhart 1, Peter Gumbsch 2,3, and Hans-Rainer Trebin 1 1 Universität Stuttgart, Institut für Theoretische
More informationMOLECULAR MODELLING OF STRESSES AND DEFORMATIONS IN NANOSTRUCTURED MATERIALS
Int. J. Appl. Math. Comput. Sci., 2004, Vol. 14, No. 4, 541 548 MOLECULAR MODELLING OF STRESSES AND DEFORMATIONS IN NANOSTRUCTURED MATERIALS GWIDON SZEFER Institute of Structural Mechanics, AGH University
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 informationSupplementary Figures
Fracture Strength (GPa) Supplementary Figures a b 10 R=0.88 mm 1 0.1 Gordon et al Zhu et al Tang et al im et al 5 7 6 4 This work 5 50 500 Si Nanowire Diameter (nm) Supplementary Figure 1: (a) TEM image
More informationLong-Term Atomistic Simulation of Heat and Mass Transport
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
More informationSingle Crystal Gradient Plasticity Part I
Chair for Continuum Mechanics Institute of Engineering Mechanics (Prof. Böhlke) Department of Mechanical Engineering S. Wulfinghoff, T. Böhlke, E. Bayerschen Single Crystal Gradient Plasticity Part I Chair
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
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 informationGLOBAL WELL-POSEDNESS FOR NONLINEAR NONLOCAL CAUCHY PROBLEMS ARISING IN ELASTICITY
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 55, pp. 1 7. ISSN: 1072-6691. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu GLOBAL WELL-POSEDNESS FO NONLINEA NONLOCAL
More informationVariational coarse graining of lattice systems
Variational coarse graining of lattice systems Andrea Braides (Roma Tor Vergata, Italy) Workshop Development and Analysis of Multiscale Methods IMA, Minneapolis, November 5, 2008 Analysis of complex lattice
More informationConcrete Fracture Prediction Using Virtual Internal Bond Model with Modified Morse Functional Potential
Concrete Fracture Prediction Using Virtual Internal Bond Model with Modified Morse Functional Potential Kyoungsoo Park, Glaucio H. Paulino and Jeffery R. Roesler Department of Civil and Environmental Engineering,
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
More informationJournal of Computational Physics
ournal of Computational Physics 9 () 397 3987 Contents lists available at ScienceDirect ournal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A multiscale coupling method for the
More informationAbsorbing Boundary Conditions for Molecular Dynamics and Multiscale Modeling
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Chemical and Biomolecular Engineering Publications and Other Works Chemical and Biomolecular Engineering 10-24-2007 Absorbing
More informationMulti-Scale Modeling of Physical Phenomena: Adaptive Control of Models
Multi-Scale Modeling of Physical Phenomena: Adaptive Control of Models J. Tinsley Oden, 1 Serge Prudhomme, 2 Albert Romkes, 3 and Paul Bauman 4 Institute for Computational Engineering and Sciences The
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationPROBLEM OF CRACK UNDER QUASIBRITTLE FRACTURE V.A. KOVTUNENKO
PROBLEM OF CRACK UNDER QUASIBRITTLE FRACTURE V.A. KOVTUNENKO Overview: 1. Motivation 1.1. Evolutionary problem of crack propagation 1.2. Stationary problem of crack equilibrium 1.3. Interaction (contact+cohesion)
More informationAtomistic simulations on the mobility of di- and tri-interstitials in Si
Atomistic simulations on the mobility of di- and tri-interstitials in Si related publications (since 2001): Posselt, M., Gao, F., Zwicker, D., Atomistic study of the migration of di- and tri-interstitials
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
More informationEAM. ReaxFF. PROBLEM B Fracture of a single crystal of silicon
PROBLEM B Fracture of a single crystal of silicon This problem set utilizes a new simulation method based on Computational Materials Design Facility (CMDF) to model fracture in a brittle material, silicon.
More informationStudies of Bimaterial Interface Fracture with Peridynamics Fang Wang 1, Lisheng Liu 2, *, Qiwen Liu 1, Zhenyu Zhang 1, Lin Su 1 & Dan Xue 1
International Power, Electronics and Materials Engineering Conference (IPEMEC 2015) Studies of Bimaterial Interface Fracture with Peridynamics Fang Wang 1, Lisheng Liu 2, *, Qiwen Liu 1, Zhenyu Zhang 1,
More informationModeling Transport in Heusler-based Spin Devices
Modeling Transport in Heusler-based Spin Devices Gautam Shine (Stanford) S. Manipatruni, A. Chaudhry, D. E. Nikonov, I. A. Young (Intel) Electronic Structure Extended Hückel theory Application to Heusler
More informationAM 106/206: Applied Algebra Madhu Sudan 1. Lecture Notes 12
AM 106/206: Applied Algebra Madhu Sudan 1 Lecture Notes 12 October 19 2016 Reading: Gallian Chs. 27 & 28 Most of the applications of group theory to the physical sciences are through the study of the symmetry
More informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
More information