Notes for Lecture 10
|
|
- Colin Ferguson
- 5 years ago
- Views:
Transcription
1 February 2, 26 Notes for Lecture Introduction to grid-based methods for solving Poisson and Laplace Equations Finite difference methods The basis for most grid-based finite difference methods is the Taylor s expansion: Φ(r + u) Φ(r) + u Φ(r) + 2! (u )2 Φ(r) + 3! (u )3 Φ(r) +! (u ) Φ(r) +. () We will work out some explicit formulae for a 2-dimensional regular grid with h denoting the step length. For the 2-dimensional Poisson equation we have 2 ρ(x, y) Φ(x, y). (2) ε We note that a sum of surrounding edge values gives: S A Φ(x + h, y) + Φ(x h, y) + Φ(x, y + h) + Φ(x, y h) (3) Φ(x, y) + h 2 2 Φ(x, y) + h 2 x + Φ(x, y) + (h 6...). y Similarly, a sum of surrounding corner values gives: S B Φ(x + h, y + h) + Φ(x h, y + h) + Φ(x + h, y h) + Φ(x h, y h) () ( Φ(x, y) + 2h 2 2 Φ(x, y) + h 6 x + y ) Φ(x, y) + (h 6...). x We note that we can combine these two results into the relation S A + S B 5Φ(x, y) + 3h2 2 2 Φ(x, y) + h Φ(x, y) + (h 6...). (5) This result can be written in the form; Φ(x, y) 5 S A 2 S B 3h2 ρ(x, y) + h 2 ρ(x, y). (6) ε ε
2 Φ(x,a)V 2 3 Φ(,y) 5 6 Φ(a,y) Φ(x,) Figure : 3 3 grid for solution of the Poisson equation within a 2-dimensional square. In general, the right hand side of this equation is known, and most of the left hand side of the equation, except for the boundary values are unknown. It can be used to develop a set of linear equations for the values of Φ(x, y) on the grid points. For example, consider a solution to the Laplace equation in the square region x a, y a which Φ(x, ) Φ(, y) Φ(a, y) and Φ(x, a) V. We will first analyze this system with a mesh of 9 points. In this case, φ 5 Φ( a 2, a 2 ) is unknown, while φ φ 2 φ 3 and φ φ 6 φ 7 φ 8 φ 9. For this example, Eq. 6 states φ 5 5 (φ 2 + φ + φ 6 + φ 8 ) + 2 (φ + φ 3 + φ 7 + φ 9 ) 3 V. (7) This results is within 2% of the exact answer of Φ( a 2, a 2 ).25V. If analyze this same system with the next more accurate grid, using the symmetry of the system Φ(x, y) Φ(a x, y), we have now 6 unknown values {φ 5, φ 6, φ 8, φ 9, φ, φ 2 } and boundary values φ φ 2 φ 3 and φ φ 7 φ φ 3 φ φ 5. This results in the following relations between the grid points: φ 5 5 (φ 2 + φ + φ 6 + φ 8 ) 2 (φ + φ 3 + φ 7 + φ 9 ), (8) φ 6 5 (φ 3 + φ 5 + φ 5 + φ 9 ) 2 (φ 2 + φ 2 + φ 8 + φ 8 ), (9) φ 8 5 (φ 5 + φ 7 φ 9 + φ ) 2 (φ + φ 6 + φ + φ 2 ), () φ 9 5 (φ 6 + φ 8 + φ 8 + φ 2 ) 2 (φ 5 + φ 5 + φ + φ ), ()
3 Φ(x,a)V Φ(,y) Φ(a,y) Φ(x,) Figure 2: 5 5 grid for solution of the Poisson equation within a 2-dimensional square. φ 5 (φ 8 + φ + φ 2 + φ ) 2 (φ 7 + φ 9 + φ 3 + φ 5 ), (2) φ 2 5 (φ 9 + φ + φ + φ 5 ) 2 (φ 8 + φ 8 + φ + φ ). (3) These equations can be cast into the form of a matrix problem which can be easily solved using Maple: /5 /5 /2 φ 5 3/ 2/5 / /5 φ 6 3/ /5 /2 /5 /5 /2 φ 8 V / /5 2/5 / /5 φ 9. () /5 /2 /5 φ / /5 2/5 φ 2 The solution to these equations and the exact results are found to be: φ 5 φ 6 φ 8 φ 9 φ φ V ; (exact) V. (5)
4 We see that the accuracy has improved considerably with the new mesh. Introduction to Finite element method The finite element approach is based on an expansion of the unknown electrostatic potential in terms of known grid-based functions of fixed shape. In two dimensions, using the indices {i, j} to reference the grid, we can denote the shape functions as {φ ij (x, y)}. The finite element expansion of the potential in two dimensions can take the form: πε Φ(x, y) ij ψ ij φ ij (x, y), (6) where ψ ij represents the amplitude associated with the shape function φ ij (x, y). The amplitude values can be determined for a given solution of the Poisson equation: by solving a linear algebra problem of the form 2 (πε Φ(x, y)) πρ(x, y), (7) M kl,ij ψ ij G kl, (8) ij where M kl,ij φ kl (x, y) φ ij (x, y) and G kl φ kl (x, y) πρ(x, y). (9) In obtaining this result, we have assumed that the boundary values vanish. In order for this result to be useful, we need to be able evaluate the integrals for M kl,ij and for G kl. In the latter case, we need to know the form of the charge density. The form of M kl,ij only depends upon the form of the shape functions. If we take these functions to be: where X i (x) φ ij (x, y) X i (x)y j (y), (2) { x x i h for xi h x x i + h otherwise and Y j (y) has a similar expression in the variable y. Then M kl,ij [ dxk (x) dx i (x) Y l(y)y j (y) + X k (x)x i (x) dy l(y) There are four types of non-trivial contributions to these values: xi +h xi +h, (2) ] dy j (y). (22) (X i (x)) 2 h ( u ) 2 du 2h 3, (23) (X i (x)x i+ (x)) h ( u)udu h 6, (2)
5 and xi +h xi +h ( dxi (x) ( dxi (x) ) 2 h du 2 h, (25) ) dx i+ (x) du h h. (26) These basic ingredients lead to the following distinct values for the matrix: 8 for k i and l j 3 M kl,ij for k i ± and/or l j ± 3 otherwise. (27) For problems in which the boundary values are, Eq. 8 then can be used to find all of the interior amplitudes ψ ij. In order to use this technique to solve the boundary value problem discussed above, we have to make one modification. The boundary value of Φ(x, a) V is not consistent with the derivation of Eq. (9), however, since we are only interested in the region y a, we can extend our numerical analysis to the region y a + h and require Φ(x, a + h) in addition to Φ(x, a) V. Using the same indexing as in Fig. 2, this means that ψ ψ 2 ψ 3 V. The finite element approach for this problem thus can be put into the matrix form for analysis by Maple: 8/3 /3 /3 /3 ψ 5 2/3 8/3 2/3 /3 ψ 6 /3 /3 8/3 /3 /3 /3 ψ 8 V 2/3 /3 2/3 8/3 2/3 /3 ψ 9. (28) /3 /3 8/3 /3 ψ 2/3 /3 2/3 8/3 ψ 2 The solution to these equations and the exact results are found to be: ψ 5 ψ 6 ψ 8 ψ 9 ψ ψ V ; (exact) V. (29) We see that the results are similar to those obtained using the finite difference approach.
Electrodynamics PHY712. Lecture 4 Electrostatic potentials and fields. Reference: Chap. 1 & 2 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 4 Electrostatic potentials and fields Reference: Chap. 1 & 2 in J. D. Jackson s textbook. 1. Complete proof of Green s Theorem 2. Proof of mean value theorem for electrostatic
More informationElectrodynamics PHY712. Lecture 3 Electrostatic potentials and fields. Reference: Chap. 1 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 3 Electrostatic potentials and fields Reference: Chap. 1 in J. D. Jackson s textbook. 1. Poisson and Laplace Equations 2. Green s Theorem 3. One-dimensional examples 1 Poisson
More informationSpezielle Funktionen 3 Orthogonale Polynome
Spezielle Funktionen 3 Orthogonale Polynome M. Gutting 28. Mai 2015 3.5 Application of the Legendre Polynomials in Electrostatics Let x R 3 and ρ(x) be a charge distribution with total charge R 3 ρ(x)dx.
More informationChapter 1. Introduction to Electrostatics
Chapter. Introduction to Electrostatics. Electric charge, Coulomb s Law, and Electric field Electric charge Fundamental and characteristic property of the elementary particles There are two and only two
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationSummary of time independent electrodynamics
hapter 10 Summary of time independent electrodynamics 10.1 Electrostatics Physical law oulomb s law charges as origin of electric field Superposition principle ector of the electric field E(x) in vacuum
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationAlgorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations
dmitry.lyakhov@kaust.edu.sa Symbolic Computations May 4, 2018 1 / 25 Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations Dmitry A. Lyakhov 1 1 Computational Sciences
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationNotes 19 Gradient and Laplacian
ECE 3318 Applied Electricity and Magnetism Spring 218 Prof. David R. Jackson Dept. of ECE Notes 19 Gradient and Laplacian 1 Gradient Φ ( x, y, z) =scalar function Φ Φ Φ grad Φ xˆ + yˆ + zˆ x y z We can
More informationMAT 419 Lecture Notes Transcribed by Eowyn Cenek 6/1/2012
(Homework 1: Chapter 1: Exercises 1-7, 9, 11, 19, due Monday June 11th See also the course website for lectures, assignments, etc) Note: today s lecture is primarily about definitions Lots of definitions
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationModern Algebra Math 542 Spring 2007 R. Pollack Solutions for HW #5. 1. Which of the following are examples of ring homomorphisms? Explain!
Modern Algebra Math 542 Spring 2007 R. Pollack Solutions for HW #5 1. Which of the following are examples of ring homomorphisms? Explain! (a) φ : R R defined by φ(x) = 2x. This is not a ring homomorphism.
More informationSelected electrostatic problems (Lecture 14)
Selected electrostatic problems (Lecture 14) February 1, 2016 236/441 Lecture outline In this lecture, we demonstrate how to solve several electrostatic problems related to calculation of the fields generated
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationCOMPLETION OF A METRIC SPACE
COMPLETION OF A METRIC SPACE HOW ANY INCOMPLETE METRIC SPACE CAN BE COMPLETED REBECCA AND TRACE Given any incomplete metric space (X,d), ( X, d X) a completion, with (X,d) ( X, d X) where X complete, and
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week 9 November 13 November Deadline to hand in the homeworks: your exercise class on week 16 November 20 November Exercises (1) Show that if T B(X, Y ) and S B(Y, Z)
More informationLecture 2 Some Sources of Lie Algebras
18.745 Introduction to Lie Algebras September 14, 2010 Lecture 2 Some Sources of Lie Algebras Prof. Victor Kac Scribe: Michael Donovan From Associative Algebras We saw in the previous lecture that we can
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
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 informationLecture-05 Perturbation Theory and Feynman Diagrams
Lecture-5 Perturbation Theory and Feynman Diagrams U. Robkob, Physics-MUSC SCPY639/428 September 3, 218 From the previous lecture We end up at an expression of the 2-to-2 particle scattering S-matrix S
More informationIntroduction to Electrostatics
Chapter 1 Introduction to Electrostatics Problem Set #1: 1.5, 1.7, 1.12 (Due Monday Feb. 11th) 1.1 Electric field Coulomb showed experimentally that for two point charges the force is -proportionaltoeachofthecharges,
More informationLecture 17. Higher boundary regularity. April 15 th, We extend our results to include the boundary. Let u C 2 (Ω) C 0 ( Ω) be a solution of
Lecture 7 April 5 th, 004 Higher boundary regularity We extend our results to include the boundary. Higher a priori regularity upto the boundary. Let u C () C 0 ( ) be a solution of Lu = f on, u = ϕ on.
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationLent 2019 VECTOR CALCULUS EXAMPLE SHEET 1 G. Taylor
Lent 29 ECTOR CALCULU EXAMPLE HEET G. Taylor. (a) The curve defined parametrically by x(t) = (a cos 3 t, a sin 3 t) with t 2π is called an astroid. ketch it, and find its length. (b) The curve defined
More informationl=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.
Electrodynamics I Exam - Part A - Closed Book KSU 15/11/6 Name Electrodynamic Score = 14 / 14 points Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try
More informationUniversity of Illinois at Chicago Department of Physics. Electricity & Magnetism Qualifying Examination
University of Illinois at Chicago Department of Physics Electricity & Magnetism Qualifying Examination January 7, 28 9. am 12: pm Full credit can be achieved from completely correct answers to 4 questions.
More informationn 2 xi = x i. x i 2. r r ; i r 2 + V ( r) V ( r) = 0 r > 0. ( 1 1 ) a r n 1 ( 1 2) V( r) = b ln r + c n = 2 b r n 2 + c n 3 ( 1 3)
Sep. 7 The L aplace/ P oisson Equations: Explicit Formulas In this lecture we study the properties of the Laplace equation and the Poisson equation with Dirichlet boundary conditions through explicit representations
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationSolutions to Problems in Jackson, Classical Electrodynamics, Third Edition. Chapter 2: Problems 11-20
Solutions to Problems in Jackson, Classical Electrodynamics, Third Edition Homer Reid December 8, 999 Chapter : Problems - Problem A line charge with linear charge density τ is placed parallel to, and
More information1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian. which is given in spherical coordinates by (2)
1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian operator (1) 2 = 2 x + 2 2 y + 2 2 z 2 which is given in spherical coordinates by (2) 2 = 1 ( r 2 ) + 1 r 2 r r r 2 sin θ θ and in
More informationLECTURE-13 : GENERALIZED CAUCHY S THEOREM
LECTURE-3 : GENERALIZED CAUCHY S THEOREM VED V. DATAR The aim of this lecture to prove a general form of Cauchy s theorem applicable to multiply connected domains. We end with computations of some real
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationGradient operator. In our calculation of dφ along the vector ds, we see that it can be described as the scalar product
Gradient operator In our calculation of dφ along the vector ds, we see that it can be described as the scalar product ( φ dφ = x î + φ y ĵ + φ ) z ˆk ( ) u x dsî + u y dsĵ + u z dsˆk We take dφ = φ ds
More informationElectric and magnetic multipoles
Electric and magnetic multipoles Trond Saue Trond Saue (LCPQ, Toulouse) Electric and magnetic multipoles Virginia Tech 2017 1 / 22 Multipole expansions In multipolar gauge the expectation value of the
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationVector and Tensor Calculus
Appendices 58 A Vector and Tensor Calculus In relativistic theory one often encounters vector and tensor expressions in both three- and four-dimensional form. The most important of these expressions are
More informationNPTEL web course on Complex Analysis. A. Swaminathan I.I.T. Roorkee, India. and. V.K. Katiyar I.I.T. Roorkee, India
NPTEL web course on Complex Analysis A. Swaminathan I.I.T. Roorkee, India and V.K. Katiyar I.I.T. Roorkee, India A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 1 / 20 Complex Analysis Module: 2:
More informationS12.1 SOLUTIONS TO PROBLEMS 12 (ODD NUMBERS)
OLUTION TO PROBLEM 2 (ODD NUMBER) 2. The electric field is E = φ = 2xi + 2y j and at (2, ) E = 4i + 2j. Thus E = 2 5 and its direction is 2i + j. At ( 3, 2), φ = 6i + 4 j. Thus the direction of most rapid
More informationEnergy method for wave equations
Energy method for wave equations Willie Wong Based on commit 5dfb7e5 of 2017-11-06 13:29 Abstract We give an elementary discussion of the energy method (and particularly the vector field method) in the
More informationContinuum Electrostatics for Ionic Solutions with Nonuniform Ionic Sizes
Continuum Electrostatics for Ionic Solutions with Nonuniform Ionic Sizes Bo Li Department of Mathematics and Center for Theoretical Biological Physics (CTBP) University of California, San Diego, USA Supported
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationLagrangian submanifolds and generating functions
Chapter 4 Lagrangian submanifolds and generating functions Motivated by theorem 3.9 we will now study properties of the manifold Λ φ X (R n \{0}) for a clean phase function φ. As shown in section 3.3 Λ
More informationRegularity for Poisson Equation
Regularity for Poisson Equation OcMountain Daylight Time. 4, 20 Intuitively, the solution u to the Poisson equation u= f () should have better regularity than the right hand side f. In particular one expects
More informationParticle in Cell method
Particle in Cell method Birdsall and Langdon: Plasma Physics via Computer Simulation Dawson: Particle simulation of plasmas Hockney and Eastwood: Computer Simulations using Particles we start with an electrostatic
More informationDulac-Cherkas function and its application to estimation of limit cycles number for some classes of planar autonomous systems
Dulac-Cherkas function and its application to estimation of limit cycles number for some classes of planar autonomous systems Castro Urdiales / September 12-16, 2011 Alexander Grin Yanka Kupala State University
More informationPart I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
More informationLecture 7 From Dirac equation to Feynman diagramms. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2
Lecture 7 From Dirac equation to Feynman diagramms SS2011: Introduction to Nuclear and Particle Physics, Part 2 2 1 Dirac equation* The Dirac equation - the wave-equation for free relativistic fermions
More informationu( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)
M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 11
WiSe 22..23 Prof. Dr. A-S. Smith Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg Problem. Theoretische Physik 2:
More informationPHY752, Fall 2016, Assigned Problems
PHY752, Fall 26, Assigned Problems For clarification or to point out a typo (or worse! please send email to curtright@miami.edu [] Find the URL for the course webpage and email it to curtright@miami.edu
More informationLecture 7. Capacitors and Electric Field Energy. Last lecture review: Electrostatic potential
Lecture 7. Capacitors and Electric Field Energy Last lecture review: Electrostatic potential V r = U r q Q Iclicker question The figure shows cross sections through two equipotential surfaces. In both
More informationCS Lecture 19. Exponential Families & Expectation Propagation
CS 6347 Lecture 19 Exponential Families & Expectation Propagation Discrete State Spaces We have been focusing on the case of MRFs over discrete state spaces Probability distributions over discrete spaces
More informationMath Exam III - Spring
Math 3 - Exam III - Spring 8 This exam contains 5 multiple choice questions and hand graded questions. The multiple choice questions are worth 5 points each and the hand graded questions are worth a total
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationJackson 2.3 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson.3 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell PROBLEM: A straight-line charge with constant linear charge λ is located perpendicular to the x-y plane in
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationNumerical Solution of Laplace Equation By Gilberto E. Urroz, October 2004
Numerical Solution of Laplace Equation By Gilberto E. Urroz, October 004 Laplace equation governs a variety of equilibrium physical phenomena such as temperature distribution in solids, electrostatics,
More informationSolutions of the Harmonic Oscillator Equation in a B-Polynomial Basis
Adv. Studies Theor. Phys., Vol. 3, 2009, no. 12, 451-460 Solutions of the Harmonic Oscillator Equation in a B-Polynomial Basis Muhammad I. Bhatti Department of Physics and Geology University of Texas-Pan
More informationMultivariate Distribution Models
Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is
More informationIntroduction to Electrostatics
Introduction to Electrostatics Charles Augustin de Coulomb (1736-1806) December 23, 2000 Contents 1 Coulomb s Law 2 2 Electric Field 4 3 Gauss s Law 7 4 Differential Form of Gauss s Law 9 5 An Equation
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More information(d). Why does this imply that there is no bounded extension operator E : W 1,1 (U) W 1,1 (R n )? Proof. 2 k 1. a k 1 a k
Exercise For k 0,,... let k be the rectangle in the plane (2 k, 0 + ((0, (0, and for k, 2,... let [3, 2 k ] (0, ε k. Thus is a passage connecting the room k to the room k. Let ( k0 k (. We assume ε k
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationModule 3: BASICS OF CFD. Part A: Finite Difference Methods
Module 3: BASICS OF CFD Part A: Finite Difference Methods THE CFD APPROACH Assembling the governing equations Identifying flow domain and boundary conditions Geometrical discretization of flow domain Discretization
More informationNumerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions
Numerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions Z. Jomaa C. Macaskill August 8, 28 Abstract We describe a 2-D finite difference algorithm for inverting
More informationChapter 4. Multiple Integrals. 1. Integrals on Rectangles
Chapter 4. Multiple Integrals Let be a rectangle in I 2 given by 1. Integrals on ectangles = [a 1, b 1 [a 2, b 2 := {(x, y) I 2 : a 1 x b 1, a 2 y b 2 }. Let P 1 be a partition of [a 1, b 1 : P 1 : a 1
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationLecture : PSPACE IP
IITM-CS6845: Theory Toolkit February 16, 2012 Lecture 22-23 : PSPACE IP Lecturer: Jayalal Sarma.M.N. Scribe: Sivaramakrishnan.N.R. Theme: Between P and PSPACE 1 Interactive Protocol for #SAT In the previous
More information18.440: Lecture 19 Normal random variables
18.440 Lecture 19 18.440: Lecture 19 Normal random variables Scott Sheffield MIT Outline Tossing coins Normal random variables Special case of central limit theorem Outline Tossing coins Normal random
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Midterm Exam, Fall x i x dx i x i x x i x dx
Physics 74 Graduate Quantum Mechanics Solutions to Midterm Exam, Fall 4. [ points] Consider the wave function x Nexp x ix (a) [6] What is the correct normaliation N? The normaliation condition is. exp,
More information1 Review of Electrostatics
1 Review of Electrostatics The first 11 sections you can omit 1.1 and 1.13) of this chapter in Jackson should only be revision; little is done there that you shouldn t have seen before presumably...).
More informationIntroduction to Probability and Stocastic Processes - Part I
Introduction to Probability and Stocastic Processes - Part I Lecture 2 Henrik Vie Christensen vie@control.auc.dk Department of Control Engineering Institute of Electronic Systems Aalborg University Denmark
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
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 informationMathematical Notes for E&M Gradient, Divergence, and Curl
Mathematical Notes for E&M Gradient, Divergence, and Curl In these notes I explain the differential operators gradient, divergence, and curl (also known as rotor), the relations between them, the integral
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationEE243 Advanced Electromagnetic Theory Lec #3: Electrostatics (Apps., Form),
EE4 Advanced Electromagnetic Theory Lec #: Electrostatics Apps., Form, Electrostatic Boundary Conditions Energy, Force and Capacitance Electrostatic Boundary Conditions on Φ Image Solutions Eample Green
More informationReflection of Plane Electromagnetic Wave from Conducting Plane
Reflection of Plane Electromagnetic Wave from Conducting Plane Zafar Turakulov August 19, 2014 Abstract The phenomenon of reflection from conducting surface is considered in terms of exact solutions of
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationCompatible Discretization Schemes
Compatible Discretization Schemes Marc Gerritsma 1 1 Collaborators Artur Palha, Guido Oud, Jasper Kreeft & Mick Bouman Marc Gerritsma Burgers Course CFD I, 25-29 January 2010 1 / 81 What if... Why mimetic
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationChapter 12: Bivariate & Conditional Distributions
Chapter 12: Bivariate & Conditional Distributions James B. Ramsey March 2007 James B. Ramsey () Chapter 12 26/07 1 / 26 Introduction Key relationships between joint, conditional, and marginal distributions.
More informationAdditional Mathematical Tools: Detail
Additional Mathematical Tools: Detail September 9, 25 The material here is not required, but gives more detail on the additional mathmatical tools: coordinate systems, rotations, the Dirac delta function
More informationSolutions Definition 2: a solution
Solutions As was stated before, one of the goals in this course is to solve, or find solutions of differential equations. In the next definition we consider the concept of a solution of an ordinary differential
More informationlaplace s method for ordinary differential equations
Physics 24 Spring 217 laplace s method for ordinary differential equations lecture notes, spring semester 217 http://www.phys.uconn.edu/ rozman/ourses/p24_17s/ Last modified: May 19, 217 It is possible
More information