Green s Functions, Boundary Integral Equations and Rotational Symmetry
|
|
- Coral Harrell
- 6 years ago
- Views:
Transcription
1 Green s Functions, Boundary Integral Equations and Rotational Symmetry...or, How to Construct a Fast Solver for Stokes Equation Saibal De Advisor: Shravan Veerapaneni University of Michigan, Ann Arbor February 2, 2018 Saibal De (U of M) Axis-symmetric BIE February 2, / 42
2 Why bother?! Particulate flow omnipresent in many physical/biological processes Porous media (soil) Suspensions (proteins in cell membranes) Sedimentation Low Reynold s number = Stokes equation Saibal De (U of M) Axis-symmetric BIE February 2, / 42
3 Outline 1 Green s Functions and Boundary Integral Equations 2 Rotational Invariance and Fourier Series 3 Application to Stokes flow 4 Numerical Aspects 5 Some Results Saibal De (U of M) Axis-symmetric BIE February 2, / 42
4 Green s Functions and Boundary Integral Equations Saibal De (U of M) Axis-symmetric BIE February 2, / 42
5 Fundamental principles of electrostatics Given charge distribution, find out potential function Two main ideas: Coulomb s law Superposition u(x) = 1 Q 4π x x Saibal De (U of M) Axis-symmetric BIE February 2, / 42
6 An example: potential due to charged spherical conductor Charge Q is dumped into a spherical conductor of radius R Can only have one type of charge (positive or negative) Charge will be distributed over the surface Charge distribution will be uniform Uniform surface charge density of σ(x ) = Q 4πR 2 Potential given by u(x) = 1 σ(x ) 4π Γ x x dγ(x ) Some symmetry arguments u(x) = 1 Q 4π x for x > R Saibal De (U of M) Axis-symmetric BIE February 2, / 42
7 A More Mathematical Point of View Rewrite u(x) = Γ [ ] 1 1 4π x x σ(x ) dγ(x ) Characteristics: A linear operator Takes information about σ on surface of ball Constructs value of u in exterior of ball The kernel G(x, x ) = 1 1 4π x x is the Green s function for 3D Laplace equation Saibal De (U of M) Axis-symmetric BIE February 2, / 42
8 Green s functions Can be defined for any linear partial differential equation, Lu = f G solves: LG(x, x ) = δ(x x ) Linearity allows us to construct more solutions: u(x) = G(x, x )σ(x ) dν(x ) Saibal De (U of M) Axis-symmetric BIE February 2, / 42
9 Layer potentials Arbitrary closed, simple surface Γ Single layer potential u(x) = [Sσ](x) := G(x, x )σ(x ) dγ(x ), Γ x Γ Double layer potential u(x) = [Dµ](x) := Γ G(x, x ) ν(x ) µ(x ) dγ(x ), x Γ Saibal De (U of M) Axis-symmetric BIE February 2, / 42
10 Extending layer potentials up to boundary Theorem Γ of class C 2, σ C(Γ) = Sσ can be extended continuously everywhere, in particular to Γ Theorem Γ of class C 2, µ C(Γ) = Dµ can be extended continuously to Γ, with limit [Dµ] ± (x) = 1 µ(x) + [Dµ](x) for x Γ 2 where [Dµ] ± (x) = lim [Dµ](x ± hν(x)), x h 0+ Γ Saibal De (U of M) Axis-symmetric BIE February 2, / 42
11 Boundary integral equations Boundary value problem Lu = 0 in domain, u = f on Γ Single layer Double layer interior [Sσ](x) = f(x) + 1 µ(x) + [Dµ](x) = f(x) 2 Double layer exterior 1 µ(x) + [Dµ](x) = f(x) 2 Saibal De (U of M) Axis-symmetric BIE February 2, / 42
12 What got swept under the rug... Is integral operator u = L σ onto? Why is it enough to have sources only on the surface? Saibal De (U of M) Axis-symmetric BIE February 2, / 42
13 References BIEs and Laplace equation Kress, Linear Integral Equations, Saibal De (U of M) Axis-symmetric BIE February 2, / 42
14 Rotational Invariance and Fourier Series Saibal De (U of M) Axis-symmetric BIE February 2, / 42
15 Exploiting rotational symmetry Saibal De (U of M) Axis-symmetric BIE February 2, / 42
16 Fourier decomposition Single layer equations Γ G(x, x )σ(x ) dγ(x ) = f(x) Multiply by 1 2π e inθ and integrate w.r.t. θ: γ [ π π ] G(x, x )e in(θ θ ) dθ σ n (r, z )r dγ(r, z ) = f n (r, z) Saibal De (U of M) Axis-symmetric BIE February 2, / 42
17 Rotationally invariant kernel Need: G(x, x ) = G(θ θ, r, z, r, z ) Single layer BIE reduces to: 2π G n (r, z, r, z )σ n (r, z )r dγ(r, z ) = f n (r, z) γ Similar formulation for other layer potentials Saibal De (U of M) Axis-symmetric BIE February 2, / 42
18 Nyström discretization Integral equation 2π G n (r, z, r, z )σ n (r, z )r dγ(r, z ) = f n (r, z) γ Pick quadrature rule {(r i, z i ), w i } I i=1 on γ Discretize I 2π G n (r i, z i, r j, z j )σ n (r j, z j )r j γ (r j, z j )w j f n (r i, z i ) j=1 A linear system with unknown σ n (r j, z j ) Saibal De (U of M) Axis-symmetric BIE February 2, / 42
19 What does G n look like? Let Then Simplify G n = G n = G(x, x ) = 1 1 4π x x 1 π e int 32π 3 π r 2 + r 2 2rr cos t + (z z ) dt 2 1 π cos nt 32π 3 π r 2 + r 2 2rr cos t + (z z ) dt 2 Write in terms of special functions, G n = 2 Q rr n 1 (χ), χ = r2 + r 2 + (z z ) 2 2 2rr Saibal De (U of M) Axis-symmetric BIE February 2, / 42
20 What are all the Q stuff? Legendre differential equation: Two solutions: Maximal P λ (χ) Minimal Q λ (χ) Recurrence relations: (1 χ 2 )u 2χu + λ(λ + 1)u = 0 Q n 1/2 (χ) = 4 n 1 2n 1 χq n 3/2(χ) 2n 3 2n 1 Q n 5/2(χ) Saibal De (U of M) Axis-symmetric BIE February 2, / 42
21 References Reduction to Fourier modes Young, Hao and Martinsson, A high-order Nystrom discretization scheme for boundary integralequations defined on rotationally symmetric surfaces, Relation to Legendre Q functions Cohl and Tohline, A Compact Cylindrical Greens Function Expansion for the Solution of Potential Problems, Saibal De (U of M) Axis-symmetric BIE February 2, / 42
22 Application to Stokes flow Saibal De (U of M) Axis-symmetric BIE February 2, / 42
23 A linearization of Navier-Stokes Incompressible, isotropic with Galilean invariant stress: t u + (u )u = 1 p + µ 2 u + g ρ 0 ρ 0 u = 0 Drop the acceleration and advection term, no external force: p + µ 2 u = 0 u = 0 Saibal De (U of M) Axis-symmetric BIE February 2, / 42
24 The Stokeslet Green s function for Stokes equation: G(x, x ) = 1 [ I 8πµ x x + J(x, ] x ) x x 3 where Single layer BIE J(x, x ) = (x x )(x x ) G(x, x )σ(x ) dγ(x ) = f(x) Γ Saibal De (U of M) Axis-symmetric BIE February 2, / 42
25 Making J rotationally invariant Define cos φ sin φ 0 U(φ) = sin φ cos φ U(φ) rotates a vector about z-axis by angle φ Define J(x, x ) = U(θ)J(x, x )U(θ ) = [U(θ)(x x )][U(θ )(x x )] Saibal De (U of M) Axis-symmetric BIE February 2, / 42
26 Rotationally invariant Stokes BIE Single layer equation, Rotationally invariant kernel Γ G(x, x )σ(x ) dγ(x ) = f(x) G(x, x ) = U(θ) G(x, x )U(θ ) Transformed equation G(x, x )[U(θ )σ(x )] dγ(x ) = [U(θ)f(x)] Γ Saibal De (U of M) Axis-symmetric BIE February 2, / 42
27 References Stokes equation and related BIEs Pozrikidis, Boundary integral and singularity methods for linearized viscous flow, 1992 Saibal De (U of M) Axis-symmetric BIE February 2, / 42
28 Numerical Aspects Saibal De (U of M) Axis-symmetric BIE February 2, / 42
29 Algorithm for single body Set up quadrature on γ; generate grid points on Γ Compute kernel matrices associated to G n Apply rigid motion, apply transformation to f Compute Fourier coefficients for f Solve the linear systems G n σ n = f n Compute σ, undo transformation and rigid motion to get σ Evaluate field u = Sσ Saibal De (U of M) Axis-symmetric BIE February 2, / 42
30 Issue: Selecting appropriate quadrature The kernel G n are singular, need to adapt Nyström Solution: Use specially designed quadratures Kapur-Rokhlin Alpert Kolm-Rokhlin (Gauss) Saibal De (U of M) Axis-symmetric BIE February 2, / 42
31 Issue: Instability of Q recurrences Recurrence relations for Q n 1/2 (χ) is unstable for all χ > 1! Solution 1: use the recursion formula differently Miller s algorithm Solution 2: don t use recursion for all values of χ For χ (say), use recursion For χ > , use FFT to evaluate integral Saibal De (U of M) Axis-symmetric BIE February 2, / 42
32 Issue: Multi-body simulations Formally, N b L ij σ j = f i j=1 After discretization, L 11 L L 1Nb σ 1 f 1 L 21 L L 2Nb σ = f 2. System is too big! L Nb 1 L Nb 2... L Nb N b σ Nb f Nb Saibal De (U of M) Axis-symmetric BIE February 2, / 42
33 Issue: Multi-body simulations Solution: Split L into two parts: L = D + B Rewrite (D + B)σ = f = (I + D 1 B)σ = D 1 f Use iterative solvers: σ k+1 = D 1 f D 1 Bσ k...or, more sophisticated GMRES Saibal De (U of M) Axis-symmetric BIE February 2, / 42
34 Issue: Particle-particle interactions Evaluating Bσ k is expensive Solution: Use FMM to reduce complexity for far off particles Saibal De (U of M) Axis-symmetric BIE February 2, / 42
35 References Singular quadrature on curves Hao, Barnett, Martinsson and Young, High-order accurate Nystrom discretization of integral equations with weakly singular kernels on smooth curves in the plane, 2011 Fast Multipole Methods Saibal De (U of M) Axis-symmetric BIE February 2, / 42
36 Some Results Saibal De (U of M) Axis-symmetric BIE February 2, / 42
37 Error for single body simulations Saibal De (U of M) Axis-symmetric BIE February 2, / 42
38 Run-time for single body simulation Saibal De (U of M) Axis-symmetric BIE February 2, / 42
39 Example flow 1: interior domain Saibal De (U of M) Axis-symmetric BIE February 2, / 42
40 Example flow 2: exterior domain Saibal De (U of M) Axis-symmetric BIE February 2, / 42
41 Next Steps Move beyond fixed bodies Resolve collision between particles Saibal De (U of M) Axis-symmetric BIE February 2, / 42
42 Thanks for Listening! Questions? Saibal De (U of M) Axis-symmetric BIE February 2, / 42
lim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationCHAPTER 3 POTENTIALS 10/13/2016. Outlines. 1. Laplace s equation. 2. The Method of Images. 3. Separation of Variables. 4. Multipole Expansion
CHAPTER 3 POTENTIALS Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 Outlines 1. Laplace s equation 2. The Method of Images 3. Separation of Variables 4. Multipole Expansion
More informationHIGH-ORDER ACCURATE METHODS FOR NYSTRÖM DISCRETIZATION OF INTEGRAL EQUATIONS ON SMOOTH CURVES IN THE PLANE
HIGH-ORDER ACCURATE METHODS FOR NYSTRÖM DISCRETIZATION OF INTEGRAL EQUATIONS ON SMOOTH CURVES IN THE PLANE S. HAO, A. H. BARNETT, P. G. MARTINSSON, AND P. YOUNG Abstract. Boundary integral equations and
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationExplicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems
z Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems Johan Helsing Lund University Talk at Integral equation methods: fast algorithms and applications, Banff,
More informationThe Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations
The Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations George Biros Lexing Ying and Denis Zorin Courant Institute of Mathematical Sciences, New York University, NY 10012,
More informationRapid evaluation of electrostatic interactions in multi-phase media
Rapid evaluation of electrostatic interactions in multi-phase media P.G. Martinsson, The University of Colorado at Boulder Acknowledgements: Some of the work presented is joint work with Mark Tygert and
More informationHigh-order quadratures for boundary integral equations: a tutorial
High-order quadratures for boundary integral equations: a tutorial CBMS conference on fast direct solvers 6/23/14 Alex Barnett (Dartmouth College) Slides accompanying a partly chalk talk. Certain details,
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
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 informationFAST INTEGRAL EQUATION METHODS FOR THE MODIFIED HELMHOLTZ EQUATION
FAST INTEGRAL EQUATION METHODS FOR THE MODIFIED HELMHOLTZ EQUATION by Bryan Quaife B.Sc., University of Calgary, 2004 M.Sc., University of Calgary, 2006 a Thesis submitted in partial fulfillment of the
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationLarge scale simulations of gravity induced sedimentation of slender fibers
Large scale simulations of gravity induced sedimentation of slender fibers Katarina Gustavsson Royal Institute of Technology (KTH) Stockholm Simu mulations of rigid fibers Fibers in Waterland 800 fibers
More informationA Fast Solver for the Stokes Equations. Tejaswin Parthasarathy CS 598APK, Fall 2017
A Fast Solver for the Stokes Equations Tejaswin Parthasarathy CS 598APK, Fall 2017 Stokes Equations? http://www.earthtimes.org/ Continuity wallpapers-s.blogspot.com @ @t + @( u i) @ i =0 http://testingstufftonight.blogspot.com/
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationA HIGH-ORDER ACCURATE ACCELERATED DIRECT SOLVER FOR ACOUSTIC SCATTERING FROM SURFACES
A HIGH-ORDER ACCURATE ACCELERATED DIRECT SOLVER FOR ACOUSTIC SCATTERING FROM SURFACES JAMES BREMER,, ADRIANNA GILLMAN, AND PER-GUNNAR MARTINSSON Abstract. We describe an accelerated direct solver for the
More informationThe spectral transform method
The spectral transform method by Nils Wedi European Centre for Medium-Range Weather Forecasts wedi@ecmwf.int Advanced Numerical Methods for Earth-System Modelling Slide 1 Advanced Numerical Methods for
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 informationPartial Differential Equations. Examples of PDEs
Partial Differential Equations Almost all the elementary and numerous advanced parts of theoretical physics are formulated in terms of differential equations (DE). Newton s Laws Maxwell equations Schrodinger
More informationNonlinear Evolution of a Vortex Ring
Nonlinear Evolution of a Vortex Ring Yuji Hattori Kyushu Institute of Technology, JAPAN Yasuhide Fukumoto Kyushu University, JAPAN EUROMECH Colloquium 491 Vortex dynamics from quantum to geophysical scales
More informationChapter 5 Fast Multipole Methods
Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationarxiv: v1 [math.na] 8 Oct 2014
SPECTRALLY-ACCURATE QUADRATURES FOR EVALUATION OF LAYER POTENTIALS CLOSE TO THE BOUNDARY FOR THE 2D STOKES AND LAPLACE EQUATIONS ALEX BARNETT, BOWEI WU, AND SHRAVAN VEERAPANENI arxiv:4.287v [math.na] 8
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 informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationCHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD
CHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD For the second chapter in the thesis, we start with surveying the mathematical background which is used directly in the Boundary Element
More informationA quadrature rule for the singular kernels of Laplace and Stokes equations over a class of parameterizable surfaces
A quadrature rule for the singular kernels of Laplace and Stokes equations over a class of parameterizable surfaces Oana Marin Linné FLOW Centre, KTH-SCI (NA), SE-100 44 Stockholm, Sweden Abstract A third
More informationξ,i = x nx i x 3 + δ ni + x n x = 0. x Dξ = x i ξ,i = x nx i x i x 3 Du = λ x λ 2 xh + x λ h Dξ,
1 PDE, HW 3 solutions Problem 1. No. If a sequence of harmonic polynomials on [ 1,1] n converges uniformly to a limit f then f is harmonic. Problem 2. By definition U r U for every r >. Suppose w is a
More informationarxiv: v1 [math.na] 15 Sep 2016
An efficient preconditioner for the fast simulation of a 2D Stokes flow in porous media arxiv:169.4484v1 [math.na] 15 Sep 216 Pieter Coulier, Bryan Quaife, and Eric Darve July 18, 218 Abstract We consider
More informationChapter 22 Gauss s Law. Copyright 2009 Pearson Education, Inc.
Chapter 22 Gauss s Law Electric Flux Gauss s Law Units of Chapter 22 Applications of Gauss s Law Experimental Basis of Gauss s and Coulomb s Laws 22-1 Electric Flux Electric flux: Electric flux through
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationMULTIPOLE EXPANSIONS IN THE PLANE
MULTIPOLE EXPANSIONS IN THE PLANE TSOGTGEREL GANTUMUR Contents 1. Electrostatics and gravitation 1 2. The Laplace equation 2 3. Multipole expansions 5 1. Electrostatics and gravitation Newton s law of
More informationLecture 1: Introduction
CBMS Conference on Fast Direct Solvers Dartmouth College June June 7, 4 Lecture : Introduction Gunnar Martinsson The University of Colorado at Boulder Research support by: Many thanks to everyone who made
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationFast numerical methods for solving linear PDEs
Fast numerical methods for solving linear PDEs P.G. Martinsson, The University of Colorado at Boulder Acknowledgements: Some of the work presented is joint work with Vladimir Rokhlin and Mark Tygert at
More information1. (3) Write Gauss Law in differential form. Explain the physical meaning.
Electrodynamics I Midterm Exam - Part A - Closed Book KSU 204/0/23 Name Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try to tell about the physics involved,
More informationFFTs in Graphics and Vision. Fast Alignment of Spherical Functions
FFTs in Graphics and Vision Fast Alignment of Spherical Functions Outline Math Review Fast Rotational Alignment Review Recall 1: We can represent any rotation R in terms of the triplet of Euler angles
More informationOn the hydrodynamic diffusion of rigid particles
On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:
More informationLearning the Linear Dynamical System with ASOS ( Approximated Second-Order Statistics )
Learning the Linear Dynamical System with ASOS ( Approximated Second-Order Statistics ) James Martens University of Toronto June 24, 2010 Computer Science UNIVERSITY OF TORONTO James Martens (U of T) Learning
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationFinite and Boundary Element Methods in Acoustics
Finite and Boundary Element Methods in Acoustics W. Kreuzer, Z. Chen, H. Waubke Austrian Academy of Sciences, Acoustics Research Institute ARI meets NuHAG Kreuzer, Chen, Waubke (ARI) FEM-BEM-FMM ARI meets
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
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 informationFast direct solvers for elliptic PDEs
Fast direct solvers for elliptic PDEs Gunnar Martinsson The University of Colorado at Boulder PhD Students: Tracy Babb Adrianna Gillman (now at Dartmouth) Nathan Halko (now at Spot Influence, LLC) Sijia
More informationSUMMARY PHYSICS 707 Electrostatics. E(x) = 4πρ(x) and E(x) = 0 (1)
SUMMARY PHYSICS 707 Electrostatics The basic differential equations of electrostatics are E(x) = 4πρ(x) and E(x) = 0 (1) where E(x) is the electric field and ρ(x) is the electric charge density. The field
More informationBoundary Element Method for Stokes Flow in Incompressible Newtonian Fluids. Nurbek Akhmetuly Tazhimbetov
Boundary Element Method for Stokes Flow in Incompressible Newtonian Fluids by Nurbek Akhmetuly Tazhimbetov Submitted to the Department of Mathematics in partial fulfillment of the requirements for the
More informationSolutions to PS 2 Physics 201
Solutions to PS Physics 1 1. ke dq E = i (1) r = i = i k eλ = i k eλ = i k eλ k e λ xdx () (x x) (x x )dx (x x ) + x dx () (x x ) x ln + x x + x x (4) x + x ln + x (5) x + x To find the field for x, we
More informationFFTs in Graphics and Vision. The Laplace Operator
FFTs in Graphics and Vision The Laplace Operator 1 Outline Math Stuff Symmetric/Hermitian Matrices Lagrange Multipliers Diagonalizing Symmetric Matrices The Laplacian Operator 2 Linear Operators Definition:
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationMath 575-Lecture Fundamental solutions for point force: Stokeslet, stresslet
Math 575-Lecture 2 Fundamental solutions for point force: tokeslet, stresslet Consider putting a point singular force with strength F in a tokes flow of infinite extent. We then need to solve P µ U = δ(x)f,
More informationMath 671: Tensor Train decomposition methods
Math 671: Eduardo Corona 1 1 University of Michigan at Ann Arbor December 8, 2016 Table of Contents 1 Preliminaries and goal 2 Unfolding matrices for tensorized arrays The Tensor Train decomposition 3
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationAn explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces
An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces Helsing, Johan; Karlsson, Anders Published in: Journal of Computational Physics DOI: 10.1016/j.jcp.014.04.053
More informationCurvilinear coordinates
C Curvilinear coordinates The distance between two points Euclidean space takes the simplest form (2-4) in Cartesian coordinates. The geometry of concrete physical problems may make non-cartesian coordinates
More informationJunior-level Electrostatics Content Review
Junior-level Electrostatics Content Review Please fill out the following exam to the best of your ability. This will not count towards your final grade in the course. Do your best to get to all the questions
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 informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationHigh-Order Corrected Trapezoidal Quadrature Rules for Functions with a Logarithmic Singularity in 2-D
High-Order Corrected Trapezoidal Quadrature Rules for Functions with a Logarithmic Singularity in 2-D Juan C. Aguilar 1 Yu Chen 2 April 24, 2002 Abstract In this report we construct correction coefficients
More informationA Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles
A Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles Aleksandar Donev Courant Institute, New York University APS DFD Meeting San Francisco, CA Nov 23rd 2014 A. Donev (CIMS)
More informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
More informationLecture 3. Electric Field Flux, Gauss Law
Lecture 3. Electric Field Flux, Gauss Law Attention: the list of unregistered iclickers will be posted on our Web page after this lecture. From the concept of electric field flux to the calculation of
More informationProblem Set #4: 4.1,4.7,4.9 (Due Monday, March 25th)
Chapter 4 Multipoles, Dielectrics Problem Set #4: 4.,4.7,4.9 (Due Monday, March 5th 4. Multipole expansion Consider a localized distribution of charges described by ρ(x contained entirely in a sphere of
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationThe Helmholtz theorem at last!
Problem. The Helmholtz theorem at last! Recall in class the Helmholtz theorem that says that if if E =0 then E can be written as E = φ () if B =0 then B can be written as B = A (2) (a) Let n be a unit
More informationProblem Set #5: 5.7,5.9,5.13 (Due Monday, April 8th)
Chapter 5 Magnetostatics Problem Set #5: 5.7,5.9,5.13 (Due Monday, April 8th) 5.1 Biot-Savart Law So far we were concerned with static configuration of charges known as electrostatics. We now switch to
More informationCourse Requirements. Course Mechanics. Projects & Exams. Homework. Week 1. Introduction. Fast Multipole Methods: Fundamentals & Applications
Week 1. Introduction. Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov What are multipole methods and what is this course about. Problems from phsics, mathematics,
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationThe Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011
The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor
More informationFast direct solvers for elliptic PDEs
Fast direct solvers for elliptic PDEs Gunnar Martinsson The University of Colorado at Boulder Students: Adrianna Gillman (now at Dartmouth) Nathan Halko Sijia Hao Patrick Young (now at GeoEye Inc.) Collaborators:
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationFast Multipole Methods for Incompressible Flow Simulation
Fast Multipole Methods for Incompressible Flow Simulation Nail A. Gumerov & Ramani Duraiswami Institute for Advanced Computer Studies University of Maryland, College Park Support of NSF awards 0086075
More informationChapter 4. The dominated convergence theorem and applications
Chapter 4. The dominated convergence theorem and applications The Monotone Covergence theorem is one of a number of key theorems alllowing one to exchange limits and [Lebesgue] integrals (or derivatives
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 informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationProblem Set #3: 2.11, 2.15, 2.21, 2.26, 2.40, 2.42, 2.43, 2.46 (Due Thursday Feb. 27th)
Chapter Electrostatics Problem Set #3:.,.5,.,.6,.40,.4,.43,.46 (Due Thursday Feb. 7th). Coulomb s Law Coulomb showed experimentally that for two point charges the force is - proportional to each of the
More informationUnique equilibrium states for geodesic flows in nonpositive curvature
Unique equilibrium states for geodesic flows in nonpositive curvature Todd Fisher Department of Mathematics Brigham Young University Fractal Geometry, Hyperbolic Dynamics and Thermodynamical Formalism
More informationPhysics 11b Lecture #3. Electric Flux Gauss s Law
Physics 11b Lecture #3 lectric Flux Gauss s Law What We Did Last Time Introduced electric field by Field lines and the rules From a positive charge to a negative charge No splitting, merging, or crossing
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationAnalysis of eddy currents in a gradient coil
Analysis of eddy currents in a gradient coil J.M.B. Kroot Eindhoven University of Technology P.O.Box 53; 56 MB Eindhoven, The Netherlands Abstract To model the z-coil of an MRI-scanner, a set of circular
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part IB Thursday 7 June 2007 9 to 12 PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question in Section
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 informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationFastMMLib: a generic Fast Multipole Method library. Eric DARRIGRAND. IRMAR Université de Rennes 1
FastMMLib: a generic Fast Multipole Method library Eric DARRIGRAND joint work with Yvon LAFRANCHE IRMAR Université de Rennes 1 Outline Introduction to FMM Historically SVD principle A first example of
More informationImplicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments
Implicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments James V. Lambers September 24, 2008 Abstract This paper presents a reformulation
More informationFast multipole boundary element method for the analysis of plates with many holes
Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational
More informationDynamic Programming. Macro 3: Lecture 2. Mark Huggett 2. 2 Georgetown. September, 2016
Macro 3: Lecture 2 Mark Huggett 2 2 Georgetown September, 2016 Three Maps: Review of Lecture 1 X = R 1 + and X grid = {x 1,..., x n } X where x i+1 > x i 1. T (v)(x) = max u(f (x) y) + βv(y) s.t. y Γ 1
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationBefore you begin read these instructions carefully.
MATHEMATICAL TRIPOS Part IB Thursday, 6 June, 2013 9:00 am to 12:00 pm PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationAverage Electrostatic Potential over a Spherical Surface
Average Electrostatic Potential over a Spherical Surface EE 141 Lecture Notes Topic 8 Professor K. E. Oughstun School of Engineering College of Engineering & Mathematical Sciences University of Vermont
More information