MSE 561, Atomic Modeling in Material Science Assignment 2
|
|
- Edwin Mason
- 6 years ago
- Views:
Transcription
1 Department of Material Science and Engineering, University of Pennsylvania MSE 561, Atomic Modeling in Material Science Assignment 2 1. Construction of the Two-dimensional Block As derived in Assignment 1, rmin for the Lennard-Jones potential is, Yang Lu So the block we construct is a square that contains 400 atoms with equal spacing 3.822Å. The subroutine named SETUP is shown below, SUBROUTINE SETUP(M) INTEGER I, J, N REAL RMIN REAL M(400,2) This block is stored in M, which contains coordinates of 400 atoms OPEN (1, FILE='BLOCK') RMIN = DO I=1, 20 FORTRAN I CODES represents ARE the i th DELETED atom in y direction IN THIS COPY DO J=1, 20 CONTACT J represents ME IF the YOU j th atom NEED in x direction MORE INFO. N=(I-1)*20+J Current number of atom whose coordinates are calculated M(N,1)=(J-1)*RMIN x coordinate is calculated using j M(N,2)=(I-1)*RMIN y coordinate is calculated using i WRITE(1,*) M(N,1),M(N,2) CLOSE(1) Finally, this block is shown in Fig.1 (next page). Each atom is donated by a red circle. 1 6 r min = 2σ = 3.822Å
2 2. Generation of Neighbor List Fig.1 In order to decrease the calculation time, we need to generate a neighbor list that can be used for future calculation. The method is: 1) calculate the distance between every two atoms in this block. 2) examine whether this distance is smaller than cut-off distance (rcut) in Lennard-Jones potential. (We only consider the neighbors that interact with the central atom.) The corresponding subroutine named NEIGHBOR is shown below, SUBROUTINE NEIGHBOR(M,DIST,NEIGHBORS,NOFNEIGHBORS) INTEGER I, J, N INTEGER NEIGHBORS(400,400), NOFNEIGHBORS(400) N represents the total number of neighbors of each atom REAL RCUT FORTRAN CODES ARE DELETED IN THIS COPY REAL M(400,2), CONTACT DIST(400,400) ME Coordinates IF YOU NEED of each MORE atom is imported INFO. in array M RCUT = OPEN(1, FILE='NEIGHBOR_LIST') DO I=1,400 WRITE(1,*) 'NEIGHBORS OF',I,':' This loop goes through each atom 2
3 N=0 Reset N for each atom at the beginning of calculation DO J=1,400 DIST(I,J)=((M(I,1)-M(J,1))**2+(M(I,2)-M(J,2))**2)**(0.5) IF (DIST(I,J).LT. RCUT.AND. DIST(I,J) +.GT. 0.0) THEN N=N+1 Calculate the distance between each atom and atom I Exam whether J is within the interaction region and exclude I itself If so, count N NEIGHBORS(I,N)=J FORTRAN CODES Record atom ARE J DELETED as neighbor of IN I in THIS the array COPY NEIGHBORS WRITE(1,*) CONTACT J ME IF YOU NEED MORE INFO. IF NOFNEIGHBORS(I)=N Record N after examining all the possible J atoms WRITE(1,*) 'NUMBER OF NEIGHBORS:',N CLOSE(1) Several representative points are shown below to justify the generated neighbor list: NEIGHBORS OF 1 : NUMBER OF NEIGHBORS: 3 NEIGHBORS OF 2 :
4 NUMBER OF NEIGHBORS: 5 NEIGHBORS OF 22 : NUMBER OF NEIGHBORS: 8 3. Atomic Level Stress In the case of pair potentials, the atomic level stress can be written as (when particles have no velocities), i σ αβ = 1 Ω i dφ( ) k,k i d r α β ik Eq.(1) where Ωi is the average volume of one atom, α and β represent x axis and y axis respectively. Under Lennard-Jones potential, Eq.(1) can be written as (if rik < 7Å), i σ αβ = 1 Ω i 24ε k,k i r [ 2(σ ) 12 + ( σ ) 6 ] r α β ik Eq.(2) If rik > 7Å, i σ αβ = 1 [3A(r Ω ik r cut ) 2 + 2B( r cut )] r α β ik i k,k i Eq.(3) When doing the summation, we don t need to go through all the 400 atoms. Instead, we only consider the k atom in neighbor list of the i atom (based on the previous subroutine). So the subroutine to calculate atomic level stress named STRESS can be coded as, 4
5 SUBROUTINE STRESS(M,DIST,NEIGHBORS,NOFNEIGHBORS) INTEGER I, J, K INTEGER NEIGHBORS(400,400), NOFNEIGHBORS(400) REAL M(400,2), DIST(400,400) REAL FD, e, A, B, VI, RMIN REAL X, Y REAL SIG(400,4) Array SIG stores stress tensor (xx, xy, yx, yy) of all the 400 atoms OPEN(1,FILE='STRESS') RMIN = e = A = E-3 B = E-3 VI = 19.0*RMIN*19.0*RMIN/400.0 Calculate the average volume of one atom, Ωi DO I=1,400 SIG(I,1)=0 FORTRAN Reset the CODES stress tensor ARE for each DELETED i atom IN THIS COPY CONTACT ME IF YOU NEED MORE INFO. SIG(I,2)=0 SIG(I,3)=0 SIG(I,4)=0 DO J=1,NOFNEIGHBORS(I) For each I atom, go through J atom in I s neighbor list K=NEIGHBORS(I,J) Get the K atom in summation X=M(I,1)-M(K,1) Calculate rik(α) Y=M(I,2)-M(K,2) Calculate rik(β) IF (DIST(I,K).LT. 7.0) THEN If rik < 7.0Å FD=24.0*e/DIST(I,K)*(-2.0*(3.405/DIST(I,K))**12+ Calculate dφ/dr using Eq.(2) +(3.405/DIST(I,K))**6) ELSE FD=3.0*A*(DIST(I,K)-7.5)**2+2.0*B*(DIST(I,K)-7.5) Calculate dφ/dr using Eq.(3) IF SIG(I,1)=SIG(I,1)+FD*X*X/DIST(I,K)/VI SIG(I,2)=SIG(I,2)+FD*X*Y/DIST(I,K)/VI 5
6 SIG(I,4)=SIG(I,4)+FD*Y*Y/DIST(I,K)/VI SIG(I,3)=SIG(I,3)+FD*Y*X/DIST(I,K)/VI FORTRAN CODES ARE DELETED IN THIS COPY WRITE (1,*) 'ATOM',I,SIG(I,1),SIG(I,2),SIG(I,3),SIG(I,4) CONTACT ME IF YOU NEED MORE INFO. CLOSE(1) Due to the 4-fold symmetry: σxx=σyy and σxy=σyx =0 except for the 4 atoms at corner. They are shown in Fig.2. Fig.2 6
7 4. Radial Distribution Function Let ni(r) be number of atoms found at distances between r and r+δr from atom i. The average number of atoms found between r and r+δr from any atom of the system as, n(r) = 1 N n i (r) where N is the total number of atoms. Then RDF can be calculated using, N i=1 g(r) = 1 n(r) ρ 2πrΔr Eq.(4) Eq.(5) In two dimension, ρ = 1/Ωi as defined in Eq.(2). After Δr is given, we can simply go through every atom with a certain r, get ni(r) in Eq.(4) with distance between any two atoms calculated from previous subroutine, sum them together and calculate g(r) by Eq.(5). The corresponding subroutine named RDF is shown below, SUBROUTINE RDF(DIST) INTEGER I, J, K, N REAL RMIN, DELTAR, PI, RHO REAL R(1001), DIST(400,400) The array DIST is imported from previous subroutine REAL G G represents g(r) for each r OPEN(1,FILE='RDF') PI = RMIN = FORTRAN CODES ARE DELETED IN THIS COPY DELTAR = RMIN/100.0 CONTACT Δr is ME chosen IF YOU to be rmin/100 NEED MORE INFO. RHO = 400.0/(19.0*RMIN*19.0*RMIN) Calculate ρ DO I=1,1001 R(I)=(I-1.0)*0.015 r is increased from 0 to15å with each step as Å N=0 Reset Σni(r) for each r at the beginning of calculation DO J=1, 400 Go through each atom J DO K=1, 400 IF (DIST(J,K).GT. R(I)) THEN 7
8 IF (DIST(J,K).LT. R(I)+DELTAR) THEN N=N+1 Count Σni(r) if atom I is within r and r+δr of atom J IF IF FORTRAN CODES ARE DELETED IN THIS COPY CONTACT ME IF YOU NEED MORE INFO. G=N/400.0/RHO/2.0/PI/R(I)/DELTAR Calculate g(r) using Eq.(5) WRITE (1,*) R(I), G CLOSE(1) g(r) is shown in Fig.3, which illustrates 8 neighbors within the range from 0 to 15Å. 5. Main Program Fig.3 Finally the main program is coded as shown below, PROGRAM MAIN 8
9 REAL M(400,2), DIST(400,400) REAL NEIGHBORS(400,400), NOFNEIGHBORS(400) CALL SETUP(M) CALL NEIGHBOR(M,DIST,NEIGHBORS,NOFNEIGHBORS) FORTRAN CODES ARE DELETED IN THIS COPY CALL STRESS(M,DIST,NEIGHBORS,NOFNEIGHBORS) CONTACT ME IF YOU NEED MORE INFO. CALL RDF(DIST) 9
MSE 561, Atomic Modeling in Material Science Assignment 1
Depatment of Mateial Science and Engineeing, Univesity of Pennsylvania MSE 561, Atomic Modeling in Mateial Science Assignment 1 Yang Lu 1. Analytical Solution The close-packed two-dimensional stuctue is
More informationTensor Transformations and the Maximum Shear Stress. (Draft 1, 1/28/07)
Tensor Transformations and the Maximum Shear Stress (Draft 1, 1/28/07) Introduction The order of a tensor is the number of subscripts it has. For each subscript it is multiplied by a direction cosine array
More informationT S z 1,S z 2,..., S z N = S z N,S z 1,..., S z N 1
Momentum states (translationally invariant systems) A periodic chain (ring), translationally invariant eigenstates with a fixed momentum (crystal momentum ) quantum number k T n =e ik n k = m 2π, m =0,...,
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationExact diagonalization methods
Summer School on Computational Statistical Physics August 4-11, 2010, NCCU, Taipei, Taiwan Exact diagonalization methods Anders W. Sandvik, Boston University Representation of states in the computer bit
More information1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )
Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationMath 114: Make-up Final Exam. Instructions:
Math 114: Make-up Final Exam Instructions: 1. Please sign your name and indicate the name of your instructor and your teaching assistant: A. Your Name: B. Your Instructor: C. Your Teaching Assistant: 2.
More information1 Bulk Simulations in a HCP lattice
1 Bulk Simulations in a HCP lattice 1.1 Introduction This project is a continuation of two previous projects that studied the mechanism of melting in a FCC and BCC lattice. The current project studies
More informationANSWERS TO PROBLEM SET 1
CHM1485: Molecular Dynamics and Chemical Dynamics in Liquids ANSWERS TO PROBLEM SET 1 1 Properties of the Liouville Operator a. Give the equations of motion in terms of the Hamiltonian. Show that they
More informationReal Numbers. Real numbers are divided into two types, rational numbers and irrational numbers
Real Numbers Real numbers are divided into two types, rational numbers and irrational numbers I. Rational Numbers: Any number that can be expressed as the quotient of two integers. (fraction). Any number
More informationStability analysis of graphite crystal lattice with moment interactions
Proc. of XXXIV Summer School "Advanced Problems in Mechanics", St.-Petersburg, Russia. 2007. (Accepted) Stability analysis of graphite crystal lattice with moment interactions Igor E. Berinskiy berigor@mail.ru
More informationClass Notes: Solving Simultaneous Linear Equations by Gaussian Elimination. Consider a set of simultaneous linear equations:
METR/OCN 465: Computer Programming with Applications in Meteorology and Oceanography Dr Dave Dempsey Dept of Geosciences, SFSU Class Notes: Solving Simultaneous Linear Equations by Gaussian Elimination
More informationSignatures of GL n Multiplicity Spaces
Signatures of GL n Multiplicity Spaces UROP+ Final Paper, Summer 2016 Mrudul Thatte Mentor: Siddharth Venkatesh Project suggested by Pavel Etingof September 1, 2016 Abstract A stable sequence of GL n representations
More informationMSE 102, Fall 2014 Midterm #2. Write your name here [10 points]:
MSE 102, Fall 2014 Midterm #2 Write your name here [10 points]: Instructions: Answer all questions to the best of your abilities. Be sure to write legibly and state your answers clearly. The point values
More informationWorked Examples Set 2
Worked Examples Set 2 Q.1. Application of Maxwell s eqns. [Griffiths Problem 7.42] In a perfect conductor the conductivity σ is infinite, so from Ohm s law J = σe, E = 0. Any net charge must be on the
More informationOWELL WEEKLY JOURNAL
Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --
More informationRandom Variable. Pr(X = a) = Pr(s)
Random Variable Definition A random variable X on a sample space Ω is a real-valued function on Ω; that is, X : Ω R. A discrete random variable is a random variable that takes on only a finite or countably
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More information2. (25%) Suppose you have an array of 1234 records in which only a few are out of order and they are not very far from their correct positions. Which
COT 6401 The Analysis of Algorithms Midterm Test Open books and notes Name - SSN - 1. (25%) Suppose that the function F is dened for all power of 2 and is described by the following recurrence equation
More informationTopological insulator part I: Phenomena
Phys60.nb 5 Topological insulator part I: Phenomena (Part II and Part III discusses how to understand a topological insluator based band-structure theory and gauge theory) (Part IV discusses more complicated
More informationMath 175 Common Exam 2A Spring 2018
Math 175 Common Exam 2A Spring 2018 Part I: Short Form The first seven (7) pages are short answer. You don t need to show work. Partial credit will be rare and small. 1. (8 points) Suppose f(x) is a function
More informationLecture 5: Random Walks and Markov Chain
Spectral Graph Theory and Applications WS 20/202 Lecture 5: Random Walks and Markov Chain Lecturer: Thomas Sauerwald & He Sun Introduction to Markov Chains Definition 5.. A sequence of random variables
More informationWORKSHEET #13 MATH 1260 FALL 2014
WORKSHEET #3 MATH 26 FALL 24 NOT DUE. Short answer: (a) Find the equation of the tangent plane to z = x 2 + y 2 at the point,, 2. z x (, ) = 2x = 2, z y (, ) = 2y = 2. So then the tangent plane equation
More informationdt Now we will look at the E&M force on moving charges to explore the momentum conservation law in E&M.
. Momentum Conservation.. Momentum in mechanics In classical mechanics p = m v and nd Newton s law d p F = dt If m is constant with time d v F = m = m a dt Now we will look at the &M force on moving charges
More informationCSE 311: Foundations of Computing. Lecture 14: Induction
CSE 311: Foundations of Computing Lecture 14: Induction Mathematical Induction Method for proving statements about all natural numbers A new logical inference rule! It only applies over the natural numbers
More information1) Get physics under control! Magnitude of force of gravity.! Notice from symmetry is 1-D problem.!
1) Get physics under control! Magnitude of force of gravity.! Notice from symmetry is 1-D problem.! g = GM e r 2, r R e g = GM( r)r, r R e 2) Get MATLAB under control! g is a vector! in this case we will
More informationHW 3 Solution Key. Classical Mechanics I HW # 3 Solution Key
Classical Mechanics I HW # 3 Solution Key HW 3 Solution Key 1. (10 points) Suppose a system has a potential energy: U = A(x 2 R 2 ) 2 where A is a constant and R is me fixed distance from the origin (on
More informationNIELINIOWA OPTYKA MOLEKULARNA
NIELINIOWA OPTYKA MOLEKULARNA chapter 1 by Stanisław Kielich translated by:tadeusz Bancewicz http://zon8.physd.amu.edu.pl/~tbancewi Poznan,luty 2008 ELEMENTS OF THE VECTOR AND TENSOR ANALYSIS Reference
More informationComputing With Tensors: Potential Applications of Physics-Motivated Mathematics to Computer Science
Computing With Tensors: Potential Applications of Physics-Motivated Mathematics to Computer Science Martine Ceberio and Vladik Kreinovich Department of Computer Science University of Texas at El Paso El
More informationNotes on Ewald summation techniques
February 3, 011 Notes on Ewald summation techniques Adapted from similar presentation in PHY 71 he total electrostatic potential energy of interaction between point charges {q i } at the positions {r i
More informationMATH45061: SOLUTION SHEET 1 V
1 MATH4561: SOLUTION SHEET 1 V 1.) a.) The faces of the cube remain aligned with the same coordinate planes. We assign Cartesian coordinates aligned with the original cube (x, y, z), where x, y, z 1. The
More informationAVOIDING PERMUTATIONS AND THE NARAYANA NUMBERS
J Korean Math Soc 50 (2013, No 3, pp 529 541 http://dxdoiorg/104134/jkms2013503529 AVOIDING PERMUTATIONS AND THE NARAYANA NUMBERS Youngja Park and Seungkyung Park Abstract We study 132 avoiding permutations
More informationMP204 Electricity and Magnetism
MATHEMATICAL PHYSICS SEMESTER 2, REPEAT 2016 2017 MP204 Electricity and Magnetism Prof. S. J. Hands, Dr. M. Haque and Dr. J.-I. Skullerud Time allowed: 1 1 2 hours Answer ALL questions MP204, 2016 2017,
More informationThe Geometry of the Quantum Hall Effect
The Geometry of the Quantum Hall Effect Dam Thanh Son (University of Chicago) Refs: Carlos Hoyos, DTS arxiv:1109.2651 DTS, M.Wingate cond-mat/0509786 Plan Review of quantum Hall physics Summary of results
More informationCrystal Polymorphism in Hydrophobically Nanoconfined Water. O. Vilanova, G. Franzese Universitat de Barcelona
Crystal Polymorphism in Hydrophobically Nanoconfined Water O. Vilanova, G. Franzese Universitat de Barcelona Water s phase diagram Phase diagram of bulk water presents up to 15 different phases Typical
More informationThe next two questions pertain to the situation described below. Consider a parallel plate capacitor with separation d:
PHYS 102 Exams Exam 2 PRINT (A) The next two questions pertain to the situation described below. Consider a parallel plate capacitor with separation d: It is connected to a battery with constant emf V.
More informationProperties of the stress tensor
Appendix C Properties of the stress tensor Some of the basic properties of the stress tensor and traction vector are reviewed in the following. C.1 The traction vector Let us assume that the state of stress
More informationParallel Scientific Computing
IV-1 Parallel Scientific Computing Matrix-vector multiplication. Matrix-matrix multiplication. Direct method for solving a linear equation. Gaussian Elimination. Iterative method for solving a linear equation.
More informationPHY2049 Fall 2000 FINAL EXAM
PHY2049 Fall 2000 FINAL EXAM. Four charges (two of them have charge +Q, the other Q) are placed in the corners of a square with a side a. Find the magnitude of forces acting on each of the charges (Q =
More informations 1 if xπy and f(x) = f(y)
Algorithms Proessor John Rei Hash Function : A B ALG 4.2 Universal Hash Functions: CLR - Chapter 34 Auxillary Reading Selections: AHU-Data Section 4.7 BB Section 8.4.4 Handout: Carter & Wegman, "Universal
More informationFinite-temperature equation of state. T ln 2sinh h" '-
Finite-temperature equation of state * $ F(V,T) = U 0 + k B T ln 2sinh h" '- #, & )/ + % 2k B T (. 0 # Compute vibrational modes, frequencies Evaluate at a given volume V Compute F at various temperatures
More informationMultipole moments. Dipole moment. The second moment µ is more commonly called the dipole moment, of the charge. distribution and is a vector
Dipole moment Multipole moments The second moment µ is more commonly called the dipole moment, of the charge distribution and is a vector µ = µ x ˆx + µ y ŷ + µ z ẑ where the α component is given by µ
More information(a) Show that electric field and magnetic field have units (force)/area or energy/volume.
Problem. Units (a) how that electric field and magnetic field have units (force)/area or energy/volume. (b) A rule of thumb that you may need in the lab is that coaxial cable has a capcitance of 2 pf/foot.
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationCSC 125 :: Final Exam May 3 & 5, 2010
CSC 125 :: Final Exam May 3 & 5, 2010 Name KEY (1 5) Complete the truth tables below: p Q p q p q p q p q p q T T T T F T T T F F T T F F F T F T T T F F F F F F T T 6-15. Match the following logical equivalences
More informationHomework 10 Solution
CS 174: Combinatorics and Discrete Probability Fall 2012 Homewor 10 Solution Problem 1. (Exercise 10.6 from MU 8 points) The problem of counting the number of solutions to a napsac instance can be defined
More informationG13CGF NAG Fortran Library Routine Document
G13 Time Series Analysis G13CGF NAG Fortran Library Routine Document Note. Before using this routine, please read the Users Note for your implementation to check the interpretation of bold italicised terms
More information1.9 Algebraic Expressions
1.9 Algebraic Expressions Contents: Terms Algebraic Expressions Like Terms Combining Like Terms Product of Two Terms The Distributive Property Distributive Property with a Negative Multiplier Answers Focus
More informationChapter 5.8: Bessel s equation
Chapter 5.8: Bessel s equation Bessel s equation of order ν is: x 2 y + xy + (x 2 ν 2 )y = 0. It has a regular singular point at x = 0. When ν = 0,, 2,..., this equation comes up when separating variables
More informationAdvanced Analysis of Algorithms - Homework I (Solutions)
Advanced Analysis of Algorithms - Homework I (Solutions) K. Subramani LCSEE, West Virginia University, Morgantown, WV {ksmani@csee.wvu.edu} 1 Problems 1. Given an array A of n integer elements, how would
More informationLab 1: Handout GULP: an Empirical energy code
3.320/SMA 5.107/ Atomistic Modeling of Materials Spring 2003 1 Lab 1: Handout GULP: an Empirical energy code We will be using the GULP code as our energy code. GULP is a program for performing a variety
More information000=0 001=1 010=2 011=3 100=4 101=5 110=6 111=7 000=0 100=4 010=2 110=6 001=1 101=5 011=3 111=7
8 CHAPTER : Fourier Series The second idea behind the FFT concerns how to get the correct h k s into the final DFT s of size. If the subscripts to 7 of the h k s are written in binary form, as shown in
More informationAbstract. A theoretical calculation is made of the total energy of the wave vortex in a vacuum, for the general case.
Formula for the total energy in the lok. Abstract. A theoretical calculation is made of the total energy of the wave vortex in a vacuum, for the general case. Alexander Dubinyansky and Pavel Churlyaev.
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationLecture10: Plasma Physics 1. APPH E6101x Columbia University
Lecture10: Plasma Physics 1 APPH E6101x Columbia University Last Lecture - Conservation principles in magnetized plasma frozen-in and conservation of particles/flux tubes) - Alfvén waves without plasma
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationGaussian Elimination for Linear Systems
Gaussian Elimination for Linear Systems Tsung-Ming Huang Department of Mathematics National Taiwan Normal University October 3, 2011 1/56 Outline 1 Elementary matrices 2 LR-factorization 3 Gaussian elimination
More informationNumerical Implementation of Transformation Optics
ECE 5322 21 st Century Electromagnetics Instructor: Office: Phone: E Mail: Dr. Raymond C. Rumpf A 337 (915) 747 6958 rcrumpf@utep.edu Lecture #16b Numerical Implementation of Transformation Optics Lecture
More informationIntroduction: Special Relativity
Introduction: Special Relativity Observation: The speed c e.g., the speed of light is the same in all coordinate systems i.e. an object moving with c in S will be moving with c in S Therefore: If " r!
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationWithout fully opening the exam, check that you have pages 1 through 12.
Name: Section: Recitation Instructor: INSTRUCTIONS Fill in your name, etc. on this first page. Without fully opening the exam, check that you have pages 1 through 12. Show all your work on the standard
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationToday s Outline. CS 362, Lecture 13. Matrix Chain Multiplication. Paranthesizing Matrices. Matrix Multiplication. Jared Saia University of New Mexico
Today s Outline CS 362, Lecture 13 Jared Saia University of New Mexico Matrix Multiplication 1 Matrix Chain Multiplication Paranthesizing Matrices Problem: We are given a sequence of n matrices, A 1, A
More informationCourse 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra
Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra D. R. Wilkins Contents 3 Topics in Commutative Algebra 2 3.1 Rings and Fields......................... 2 3.2 Ideals...............................
More informationReconstruction of the Dirichlet to Neumann Map for a Resistor Network
Reconstruction of the Dirichlet to Neumann Map for a Resistor Network Adrian V. Mariano University of Washington Miriam A. Myjak Seattle University June 21, 23 Abstract We describe methods of reconstructing
More informationMath 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C
Math 35 Solutions for Final Exam Page Problem. ( points) (a) ompute the line integral F ds for the path c(t) = (t 2, t 3, t) with t and the vector field F (x, y, z) = xi + zj + xk. (b) ompute the line
More informationMLC Practice Final Exam
Name: Section: Recitation/Instructor: INSTRUCTIONS Fill in your name, etc. on this first page. Without fully opening the exam, check that you have pages 1 through 13. Show all your work on the standard
More informationCISC 235: Topic 1. Complexity of Iterative Algorithms
CISC 235: Topic 1 Complexity of Iterative Algorithms Outline Complexity Basics Big-Oh Notation Big-Ω and Big-θ Notation Summations Limitations of Big-Oh Analysis 2 Complexity Complexity is the study of
More informationSolutions 2017 AB Exam
1. Solve for x : x 2 = 4 x. Solutions 2017 AB Exam Texas A&M High School Math Contest October 21, 2017 ANSWER: x = 3 Solution: x 2 = 4 x x 2 = 16 8x + x 2 x 2 9x + 18 = 0 (x 6)(x 3) = 0 x = 6, 3 but x
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationMTH 234 Exam 2 November 21st, Without fully opening the exam, check that you have pages 1 through 12.
Name: Section: Recitation Instructor: INSTRUCTIONS Fill in your name, etc. on this first page. Without fully opening the exam, check that you have pages 1 through 12. Show all your work on the standard
More informationMathematical Structures Combinations and Permutations
Definitions: Suppose S is a (finite) set and n, k 0 are integers The set C(S, k) of k - combinations consists of all subsets of S that have exactly k elements The set P (S, k) of k - permutations consists
More informationarxiv: v2 [math.co] 20 Jun 2018
ON ORDERED RAMSEY NUMBERS OF BOUNDED-DEGREE GRAPHS MARTIN BALKO, VÍT JELÍNEK, AND PAVEL VALTR arxiv:1606.0568v [math.co] 0 Jun 018 Abstract. An ordered graph is a pair G = G, ) where G is a graph and is
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationCS100: DISCRETE STRUCTURES
1 CS100: DISCRETE STRUCTURES Computer Science Department Lecture 2: Functions, Sequences, and Sums Ch2.3, Ch2.4 2.3 Function introduction : 2 v Function: task, subroutine, procedure, method, mapping, v
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More informationDepartment of Mathematics and Statistics Math B: Discrete Mathematics and Its Applications Test 1: October 19, 2006
Department of Mathematics and Statistics Math 1019 3.0 B: Discrete Mathematics and Its Applications est 1: October 19, 2006 Name: (Last/amily) (irst) Student No. Instructions: 1. his test has 8 questions
More informationPhys 2025, First Test. September 20, minutes Name:
Phys 05, First Test. September 0, 011 50 minutes Name: Show all work for maximum credit. Each problem is worth 10 points. Work 10 of the 11 problems. k = 9.0 x 10 9 N m / C ε 0 = 8.85 x 10-1 C / N m e
More informationMath 171 Spring 2017 Final Exam. Problem Worth
Math 171 Spring 2017 Final Exam Problem 1 2 3 4 5 6 7 8 9 10 11 Worth 9 6 6 5 9 8 5 8 8 8 10 12 13 14 15 16 17 18 19 20 21 22 Total 8 5 5 6 6 8 6 6 6 6 6 150 Last Name: First Name: Student ID: Section:
More informationInstructions: No books. No notes. Non-graphing calculators only. You are encouraged, although not required, to show your work.
Exam 3 Math 850-007 Fall 04 Odenthal Name: Instructions: No books. No notes. Non-graphing calculators only. You are encouraged, although not required, to show your work.. Evaluate the iterated integral
More informationBinomial Coefficient Identities/Complements
Binomial Coefficient Identities/Complements CSE21 Fall 2017, Day 4 Oct 6, 2017 https://sites.google.com/a/eng.ucsd.edu/cse21-fall-2017-miles-jones/ permutation P(n,r) = n(n-1) (n-2) (n-r+1) = Terminology
More informationPage Points Score Total: 210. No more than 200 points may be earned on the exam.
Name: PID: Section: Recitation Instructor: DO NOT WRITE BELOW THIS LINE. GO ON TO THE NEXT PAGE. Page Points Score 3 18 4 18 5 18 6 18 7 18 8 18 9 18 10 21 11 21 12 21 13 21 Total: 210 No more than 200
More information1 Hooke s law, stiffness, and compliance
Non-quilibrium Continuum Physics TA session #5 TA: Yohai Bar Sinai 3.04.206 Linear elasticity I This TA session is the first of three at least, maybe more) in which we ll dive deep deep into linear elasticity
More informationUniformity of the Universe
Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of
More informationDynamic Programming: Shortest Paths and DFA to Reg Exps
CS 374: Algorithms & Models of Computation, Spring 207 Dynamic Programming: Shortest Paths and DFA to Reg Exps Lecture 8 March 28, 207 Chandra Chekuri (UIUC) CS374 Spring 207 / 56 Part I Shortest Paths
More informationElectro Magnetic Field Dr. Harishankar Ramachandran Department of Electrical Engineering Indian Institute of Technology Madras
Electro Magnetic Field Dr. Harishankar Ramachandran Department of Electrical Engineering Indian Institute of Technology Madras Lecture - 7 Gauss s Law Good morning. Today, I want to discuss two or three
More informationChapter 28 Solutions
Chapter 8 Solutions 8.1 (a) P ( V) R becomes 0.0 W (11.6 V) R so R 6.73 Ω (b) V IR so 11.6 V I (6.73 Ω) and I 1.7 A ε IR + Ir so 15.0 V 11.6 V + (1.7 A)r r 1.97 Ω Figure for Goal Solution Goal Solution
More informationFast Multipole Methods: Fundamentals & Applications. Ramani Duraiswami Nail A. Gumerov
Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov Week 1. Introduction. What are multipole methods and what is this course about. Problems from physics, mathematics,
More informationx 2 + 6x 18 x + 2 Name: Class: Date: 1. Find the coordinates of the local extreme of the function y = x 2 4 x.
1. Find the coordinates of the local extreme of the function y = x 2 4 x. 2. How many local maxima and minima does the polynomial y = 8 x 2 + 7 x + 7 have? 3. How many local maxima and minima does the
More informationPhysics 9 Spring 2012 Midterm 1 Solutions
Physics 9 Spring 22 NAME: TA: Physics 9 Spring 22 Midterm s For the midterm, you may use one sheet of notes with whatever you want to put on it, front and back. Please sit every other seat, and please
More informationSolutions to Homework 5
Solutions to Homework 5 1. Let z = f(x, y) be a twice continuously differentiable function of x and y. Let x = r cos θ and y = r sin θ be the equations which transform polar coordinates into rectangular
More informationMathathon Round 1 (2 points each)
Mathathon Round ( points each). A circle is inscribed inside a square such that the cube of the radius of the circle is numerically equal to the perimeter of the square. What is the area of the circle?
More informationVariational construction of periodic and connecting orbits in the planar Sitnikov problem. Mitsuru Shibayama(Kyoto University)
Variational construction of periodic and connecting orbits in the planar Sitnikov problem Mitsuru Shibayama(Kyoto University) 1 The three-body problem Consider the planar three-body problem which is governed
More informationAdvanced Analysis of Algorithms - Midterm (Solutions)
Advanced Analysis of Algorithms - Midterm (Solutions) K. Subramani LCSEE, West Virginia University, Morgantown, WV {ksmani@csee.wvu.edu} 1 Problems 1. Solve the following recurrence using substitution:
More informationSimpleFEM. an introduction to the Q1P0 element and its implementation. C. Thieulot. May 2011
SimpleFEM an introduction to the Q1P0 element and its implementation C. Thieulot May 2011 Warning The following presentation is about the Finite Element Method but it isn t comprehensive Warning The following
More informationMathwithsheppard.weebly.com
Unit #: Powers and Polynomials Unit Outline: Date Lesson Title Assignment Completed.1 Introduction to Algebra. Discovering the Exponent Laws Part 1. Discovering the Exponent Laws Part. Multiplying and
More information3-D Kinetics of Rigid Bodies
3-D Kinetics of Rigid Bodies Angular Momentum Generalized Newton s second law for the motion for a 3-D mass system Moment eqn for 3-D motion will be different than that obtained for plane motion Consider
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More information