3 the Laplace Equation

Total Page:16

File Type:pdf, Size:1020Kb

3 the Laplace Equation 3 The Laplace Equation Introduction – Well Posed Problems. Uniqueness – Harmonic Functions – Fundamental Solution and Newtonian Potential – The Green Function – Uniqueness in Unbounded Domains – Surface Potentials 3.1 Introduction The Laplace equation Δu = 0 occurs frequently in applied sciences, in particular in the study of the steady state phenomena. Its solutions are called harmonic functions. For instance, the equilibrium position of a perfectly elastic membrane is a harmonic function as it is the velocity potential of a homogeneous fluid. Also, the steady state temperature of a homogeneous and isotropic body is a harmonic function and in this case Laplace equation constitutes the stationary counterpart (time independent) of the diffusion equation. Slightly more generally, Poisson’s equation Δu = f plays an important role in the theory of conservative fields (electrical, magnetic, gravitational,...) where the vector field is derived from the gradient of a potential. For example, let E be a force field due to a distribution of electric charges in a domain Ω ⊂ R3. Then, in standard units, div E =4πρ,whereρ represents the density of the charge distribution. When a potential u exists such that ∇u = −E, then Δu =div∇u = −4πρ, which is Poisson’s equation. If the electric field is created by charges located outside Ω,thenρ =0inΩ and u is harmonic therein. Analogously, the potential of a gravitational field due to a mass distribution is a harmonic function in a region free from mass. In dimension two, the theories of harmonic and holomorphic functions are strictly connected1. Indeed, the real and the imaginary part of a holomorphic 1 A complex function f = f (z)isholomorphic in an open subset Ω of the complex plane if for every z0 ∈ Ω, the limit f (z) − f (z0) lim = f (z0) z→z0 z − z0 Salsa S. Partial Differential Equations in Action: From Modelling to Theory c Springer-Verlag 2008, Milan 3.2 Well Posed Problems. Uniqueness 103 function are harmonic. For instance, since the functions zm = rm (cos mθ + i sin mθ) ,m∈ N, (r, θ polar coordinates) are holomorphic in the whole plane C, the functions u (r, θ)=rm cos mθ and v (r, θ)=rm sin mθ m ∈ N, are harmonic in R2 (called elementary harmonics). In Cartesian coordinates, they are harmonic polynomials; for m =1, 2, 3 we find x, y, xy, x2 − y2,x3 − 3xy2, 3x2y − y3. Other examples are u (x, y)=eαx cos αy, v (x, y)=eαx sin αy (α ∈ R), the real and imaginary parts of f (z)=eiαz, both harmonic in R2,and u (r, θ)=logr, v (r, θ)=θ, ≡ R2\ the real and imaginary parts of f (z)=log0 z log r + iθ,harmonicin (0, 0) and R2\{θ =0} , respectively. In this chapter we present the formulation of the most important well posed problems and the classical properties of harmonic functions, focusing mainly on dimensions two and three. As in Chapter 2, we emphasize some probabilistic as- pects, exploiting the connection among random walks, Brownian motion and the Laplace operator. A central notion is the concept of fundamental solution, that we develop in conjunction with the very basic elements of the so called potential theory. 3.2 Well Posed Problems. Uniqueness Consider the Poisson equation Δu = f in Ω (3.1) where Ω ⊂ Rn is a bounded domain. The well posed problems associated with equation (3.1) are the stationary counterparts of the corresponding problems for the diffusion equation. Clearly here there is no initial condition. On the boundary ∂Ω we may assign: • Dirichlet data u = g, (3.2) exists and it is finite. 104 3 The Laplace Equation • Neumann data ∂ν u = h, (3.3) where ν is the outward normal unit vector to ∂Ω, • a Robin (radiation) condition ∂ν u + αu = h (α>0), (3.4) • a mixed condition; for instance, u = g on ΓD (3.5) ∂ν u = h on ΓN , where ΓD ∪ ΓN = ∂Ω, ΓD ∩ ΓN = ∅,andΓN is a relatively open subset of ∂Ω. When g = h = 0 we say that the above boundary conditions are homogeneous. We give some interpretations. If u is the position of a perfectly flexible mem- brane and f is an external distributed load (vertical force per unit surface), then (3.1) models a steady state. The Dirichlet condition corresponds to fixing the position of the membrane at its boundary. Robin condition describes an elastic attachment at the boundary while a homogeneous Neumann condition corresponds to a free vertical motion of the boundary. If u is the steady state concentration of a substance, the Dirichlet condition prescribes the level of u at the boundary, while the Neumann condition assigns the flux of u through the boundary. Using Green’s identity (1.13) we can prove the following uniqueness result. Theorem 3.1. Let Ω ⊂ Rn be a smooth, bounded domain. Then there exists at most one solution u ∈ C2 (Ω) ∩ C1 Ω of (3.1), satisfying on ∂Ω one of the conditions (3.2), (3.4) or (3.5). In the case of the Neumann condition, that is when ∂ν u = h on ∂Ω, two solutions differ by a constant. Proof. Let u and v be solutions of the same problem, sharing the same boundary data, and let w = u − v.Thenw is harmonic and satisfies homogeneous boundary conditions (one among (3.2)-(3.5)). Substituting u = v = w into (1.13) we find 2 |∇w| dx = w∂ν wdσ. Ω ∂Ω If Dirichlet or mixed conditions hold, we have w∂νwdσ=0. ∂Ω 3.3 Harmonic Functions 105 When a Robin condition holds 2 w∂νwdσ= − αw dσ ≤ 0. ∂Ω ∂Ω In any case we obtain that |∇w|2 dx ≤0. (3.6) Ω From (3.6) we infer ∇w = 0 and therefore w = u − v = constant. This concludes the proof in the case of Neumann condition. In the other cases, the constant must be zero (why?), hence u = v. Remark 3.1. Consider the Neumann problem Δu = f in Ω, ∂ν u = h on ∂Ω. Integrating the equation on Ω and using Gauss’ formula we find fdx = hdσ. (3.7) Ω ∂Ω The relation (3.7) appears as a compatibility condition on the data f and h,that has necessarily to be satisfied in order for the Neumann problem to admit a so- lution. Thus, when having to solve a Neumann problem, the first thing to do is to check the validity of (3.7). If it does not hold, the problem does not have any solution. We will examine later the physical meaning of (3.7). 3.3 Harmonic Functions 3.3.1 Discrete harmonic functions In Chapter 2 we have examined the connection between Brownian motion and diffusion equation. We go back now to the multidimensional symmetric random walk considered in Section 2.6, analyzing its relation with the Laplace operator Δ. For simplicity we will work in dimension n = 2 but both arguments and conclusions may be easily extended to any dimension n>2. We fix a time step τ>0, a space 2 step h>0 and denote by hZ the lattice of points x =(x1,x2) whose coordinates are integer multiples of h.Letp (x,t)=p (x1,x2,t) be the transition probability function, giving the probability to find our random particle at x at time t.From the total probability formula we found a difference equation for p,thatwerewrite in dimension two: 1 p (x,t+ τ)= {p (x+he ,t)+p (x−he ,t)+p (x+he ,t)+p (x−he ,t)} . 4 1 1 2 2 (3.8) We can write this formula in a more significant way by introducing the mean value operator Mh, whose action on a function u = u(x) is defined by the following 106 3 The Laplace Equation formula: 1 M f (x)= {u (x+he )+u (x−he )+u (x+he )+u (x−he )} h 4 1 1 2 2 1 = u (y) . 4 |x−y|=h 2 Note that Mhu (x) gives the average of u over the points of the lattice hZ at distance h from x. We say that these points constitute the discrete neighborhood of x of radius h. It is clear that (3.8) can be written in the form p (x,t+ τ)=Mhp (x,t) . (3.9) In (3.9) the probability p at time t + τ is determined by the action of Mh at the previous time, and then it is natural to interpret the mean value operator as the generator of the random walk. Now we come to the Laplacian. If u is twice continuously differentiable, it is not difficult to show that2 M u (x) − u (x) 1 lim h → Δu (x) . (3.10) h→0 h2 4 The formula (3.10) induces to define, for any fixed h>0, a discrete Laplace operator through the formula ∗ − Δh = Mh I ∗ where I denotes the identity operator (i.e. Iu = u). The operator Δh acts on functions u defined in the whole lattice hZ2 and, coherently, we say that u is ∗ d-harmonic (d for discrete)ifΔhu =0. Thus, the value of a d-harmonic function at any point x is given by the average of the values at the points in the discrete neighborhood of x of radius h. We can proceed further and define a discrete Dirichlet problem. Let A be a subset of hZ2. We say that A is connected if, given any couple of points x0, x1 in A,itis possible to connect them by a walk3 on hZ2 entirely contained in A. Moreover, we say that x ∈A is an interior point of A if its h−neighborhood is contained in A.ThepointsofA that are not interior are called boundary points (Fig.
Recommended publications
  • Best Subordinant for Differential Superordinations of Harmonic Complex-Valued Functions
    mathematics Article Best Subordinant for Differential Superordinations of Harmonic Complex-Valued Functions Georgia Irina Oros Department of Mathematics and Computer Sciences, Faculty of Informatics and Sciences, University of Oradea, 410087 Oradea, Romania; [email protected] or [email protected] Received: 17 September 2020; Accepted: 11 November 2020; Published: 16 November 2020 Abstract: The theory of differential subordinations has been extended from the analytic functions to the harmonic complex-valued functions in 2015. In a recent paper published in 2019, the authors have considered the dual problem of the differential subordination for the harmonic complex-valued functions and have defined the differential superordination for harmonic complex-valued functions. Finding the best subordinant of a differential superordination is among the main purposes in this research subject. In this article, conditions for a harmonic complex-valued function p to be the best subordinant of a differential superordination for harmonic complex-valued functions are given. Examples are also provided to show how the theoretical findings can be used and also to prove the connection with the results obtained in 2015. Keywords: differential subordination; differential superordination; harmonic function; analytic function; subordinant; best subordinant MSC: 30C80; 30C45 1. Introduction and Preliminaries Since Miller and Mocanu [1] (see also [2]) introduced the theory of differential subordination, this theory has inspired many researchers to produce a number of analogous notions, which are extended even to non-analytic functions, such as strong differential subordination and superordination, differential subordination for non-analytic functions, fuzzy differential subordination and fuzzy differential superordination. The notion of differential subordination was adapted to fit the harmonic complex-valued functions in the paper published by S.
    [Show full text]
  • Harmonic Functions
    Lecture 1 Harmonic Functions 1.1 The Definition Definition 1.1. Let Ω denote an open set in R3. A real valued function u(x, y, z) on Ω with continuous second partials is said to be harmonic if and only if the Laplacian ∆u = 0 identically on Ω. Note that the Laplacian ∆u is defined by ∂2u ∂2u ∂2u ∆u = + + . ∂x2 ∂y2 ∂z2 We can make a similar definition for an open set Ω in R2.Inthatcase, u is harmonic if and only if ∂2u ∂2u ∆u = + =0 ∂x2 ∂y2 on Ω. Some basic examples of harmonic functions are 2 2 2 3 u = x + y 2z , Ω=R , − 1 3 u = , Ω=R (0, 0, 0), r − where r = x2 + y2 + z2. Moreover, by a theorem on complex variables, the real part of an analytic function on an open set Ω in 2 is always harmonic. p R Thus a function such as u = rn cos nθ is a harmonic function on R2 since u is the real part of zn. 1 2 1.2 The Maximum Principle The basic result about harmonic functions is called the maximum principle. What the maximum principle says is this: if u is a harmonic function on Ω, and B is a closed and bounded region contained in Ω, then the max (and min) of u on B is always assumed on the boundary of B. Recall that since u is necessarily continuous on Ω, an absolute max and min on B are assumed. The max and min can also be assumed inside B, but a harmonic function cannot have any local extrema inside B.
    [Show full text]
  • 23. Harmonic Functions Recall Laplace's Equation
    23. Harmonic functions Recall Laplace's equation ∆u = uxx = 0 ∆u = uxx + uyy = 0 ∆u = uxx + uyy + uzz = 0: Solutions to Laplace's equation are called harmonic functions. The inhomogeneous version of Laplace's equation ∆u = f is called the Poisson equation. Harmonic functions are Laplace's equation turn up in many different places in mathematics and physics. Harmonic functions in one variable are easy to describe. The general solution of uxx = 0 is u(x) = ax + b, for constants a and b. Maximum principle Let D be a connected and bounded open set in R2. Let u(x; y) be a harmonic function on D that has a continuous extension to the boundary @D of D. Then the maximum (and minimum) of u are attained on the bound- ary and if they are attained anywhere else than u is constant. Euivalently, there are two points (xm; ym) and (xM ; yM ) on the bound- ary such that u(xm; ym) ≤ u(x; y) ≤ u(xM ; yM ) for every point of D and if we have equality then u is constant. The idea of the proof is as follows. At a maximum point of u in D we must have uxx ≤ 0 and uyy ≤ 0. Most of the time one of these inequalities is strict and so 0 = uxx + uyy < 0; which is not possible. The only reason this is not a full proof is that sometimes both uxx = uyy = 0. As before, to fix this, simply perturb away from zero. Let > 0 and let v(x; y) = u(x; y) + (x2 + y2): Then ∆v = ∆u + ∆(x2 + y2) = 4 > 0: 1 Thus v has no maximum in D.
    [Show full text]
  • Harmonic Forms, Minimal Surfaces and Norms on Cohomology of Hyperbolic 3-Manifolds
    Harmonic Forms, Minimal Surfaces and Norms on Cohomology of Hyperbolic 3-Manifolds Xiaolong Hans Han Abstract We bound the L2-norm of an L2 harmonic 1-form in an orientable cusped hy- perbolic 3-manifold M by its topological complexity, measured by the Thurston norm, up to a constant depending on M. It generalizes two inequalities of Brock and Dunfield. We also study the sharpness of the inequalities in the closed and cusped cases, using the interaction of minimal surfaces and harmonic forms. We unify various results by defining two functionals on orientable closed and cusped hyperbolic 3-manifolds, and formulate several questions and conjectures. Contents 1 Introduction2 1.1 Motivation and previous results . .2 1.2 A glimpse at the non-compact case . .3 1.3 Main theorem . .4 1.4 Topology and the Thurston norm of harmonic forms . .5 1.5 Minimal surfaces and the least area norm . .5 1.6 Sharpness and the interaction between harmonic forms and minimal surfaces6 1.7 Acknowledgments . .6 2 L2 Harmonic Forms and Compactly Supported Cohomology7 2.1 Basic definitions of L2 harmonic forms . .7 2.2 Hodge theory on cusped hyperbolic 3-manifolds . .9 2.3 Topology of L2-harmonic 1-forms and the Thurston Norm . 11 2.4 How much does the hyperbolic metric come into play? . 14 arXiv:2011.14457v2 [math.GT] 23 Jun 2021 3 Minimal Surface and Least Area Norm 14 3.1 Truncation of M and the definition of Mτ ................. 16 3.2 The least area norm and the L1-norm . 17 4 L1-norm, Main theorem and The Proof 18 4.1 Bounding L1-norm by L2-norm .
    [Show full text]
  • Complex Analysis
    Complex Analysis Andrew Kobin Fall 2010 Contents Contents Contents 0 Introduction 1 1 The Complex Plane 2 1.1 A Formal View of Complex Numbers . .2 1.2 Properties of Complex Numbers . .4 1.3 Subsets of the Complex Plane . .5 2 Complex-Valued Functions 7 2.1 Functions and Limits . .7 2.2 Infinite Series . 10 2.3 Exponential and Logarithmic Functions . 11 2.4 Trigonometric Functions . 14 3 Calculus in the Complex Plane 16 3.1 Line Integrals . 16 3.2 Differentiability . 19 3.3 Power Series . 23 3.4 Cauchy's Theorem . 25 3.5 Cauchy's Integral Formula . 27 3.6 Analytic Functions . 30 3.7 Harmonic Functions . 33 3.8 The Maximum Principle . 36 4 Meromorphic Functions and Singularities 37 4.1 Laurent Series . 37 4.2 Isolated Singularities . 40 4.3 The Residue Theorem . 42 4.4 Some Fourier Analysis . 45 4.5 The Argument Principle . 46 5 Complex Mappings 47 5.1 M¨obiusTransformations . 47 5.2 Conformal Mappings . 47 5.3 The Riemann Mapping Theorem . 47 6 Riemann Surfaces 48 6.1 Holomorphic and Meromorphic Maps . 48 6.2 Covering Spaces . 52 7 Elliptic Functions 55 7.1 Elliptic Functions . 55 7.2 Elliptic Curves . 61 7.3 The Classical Jacobian . 67 7.4 Jacobians of Higher Genus Curves . 72 i 0 Introduction 0 Introduction These notes come from a semester course on complex analysis taught by Dr. Richard Carmichael at Wake Forest University during the fall of 2010. The main topics covered include Complex numbers and their properties Complex-valued functions Line integrals Derivatives and power series Cauchy's Integral Formula Singularities and the Residue Theorem The primary reference for the course and throughout these notes is Fisher's Complex Vari- ables, 2nd edition.
    [Show full text]
  • THE IMPACT of RIEMANN's MAPPING THEOREM in the World
    THE IMPACT OF RIEMANN'S MAPPING THEOREM GRANT OWEN In the world of mathematics, scholars and academics have long sought to understand the work of Bernhard Riemann. Born in a humble Ger- man home, Riemann became one of the great mathematical minds of the 19th century. Evidence of his genius is reflected in the greater mathematical community by their naming 72 different mathematical terms after him. His contributions range from mathematical topics such as trigonometric series, birational geometry of algebraic curves, and differential equations to fields in physics and philosophy [3]. One of his contributions to mathematics, the Riemann Mapping Theorem, is among his most famous and widely studied theorems. This theorem played a role in the advancement of several other topics, including Rie- mann surfaces, topology, and geometry. As a result of its widespread application, it is worth studying not only the theorem itself, but how Riemann derived it and its impact on the work of mathematicians since its publication in 1851 [3]. Before we begin to discover how he derived his famous mapping the- orem, it is important to understand how Riemann's upbringing and education prepared him to make such a contribution in the world of mathematics. Prior to enrolling in university, Riemann was educated at home by his father and a tutor before enrolling in high school. While in school, Riemann did well in all subjects, which strengthened his knowl- edge of philosophy later in life, but was exceptional in mathematics. He enrolled at the University of G¨ottingen,where he learned from some of the best mathematicians in the world at that time.
    [Show full text]
  • 21 Laplace's Equation and Harmonic Func- Tions
    21 Laplace's Equation and Harmonic Func- tions 21.1 Introductory Remarks on the Laplacian operator Given a domain Ω ⊂ Rd, then r2u = div(grad u) = 0 in Ω (1) is Laplace's equation defined in Ω. If d = 2, in cartesian coordinates @2u @2u r2u = + @x2 @y2 or, in polar coordinates 1 @ @u 1 @2u r2u = r + : r @r @r r2 @θ2 If d = 3, in cartesian coordinates @2u @2u @2u r2u = + + @x2 @y2 @z2 while in cylindrical coordinates 1 @ @u 1 @2u @2u r2u = r + + r @r @r r2 @θ2 @z2 and in spherical coordinates 1 @ @u 1 @2u @u 1 @2u r2u = ρ2 + + cot φ + : ρ2 @ρ @ρ ρ2 @φ2 @φ sin2 φ @θ2 We deal with problems involving most of these cases within these Notes. Remark: There is not a universally accepted angle notion for the Laplacian in spherical coordinates. See Figure 1 for what θ; φ mean here. (In the rest of these Notes we will use the notation r for the radial distance from the origin, no matter what the dimension.) 1 Figure 1: Coordinate angle definitions that will be used for spherical coordi- nates in these Notes. Remark: In the math literature the Laplacian is more commonly written with the symbol ∆; that is, Laplace's equation becomes ∆u = 0. For these Notes we write the equation as is done in equation (1) above. The non-homogeneous version of Laplace's equation, namely r2u = f(x) (2) is called Poisson's equation. Another important equation that comes up in studying electromagnetic waves is Helmholtz's equation: r2u + k2u = 0 k2 is a real, positive parameter (3) Again, Poisson's equation is a non-homogeneous Laplace's equation; Helm- holtz's equation is not.
    [Show full text]
  • Chapter 2 Complex Analysis
    Chapter 2 Complex Analysis In this part of the course we will study some basic complex analysis. This is an extremely useful and beautiful part of mathematics and forms the basis of many techniques employed in many branches of mathematics and physics. We will extend the notions of derivatives and integrals, familiar from calculus, to the case of complex functions of a complex variable. In so doing we will come across analytic functions, which form the centerpiece of this part of the course. In fact, to a large extent complex analysis is the study of analytic functions. After a brief review of complex numbers as points in the complex plane, we will ¯rst discuss analyticity and give plenty of examples of analytic functions. We will then discuss complex integration, culminating with the generalised Cauchy Integral Formula, and some of its applications. We then go on to discuss the power series representations of analytic functions and the residue calculus, which will allow us to compute many real integrals and in¯nite sums very easily via complex integration. 2.1 Analytic functions In this section we will study complex functions of a complex variable. We will see that di®erentiability of such a function is a non-trivial property, giving rise to the concept of an analytic function. We will then study many examples of analytic functions. In fact, the construction of analytic functions will form a basic leitmotif for this part of the course. 2.1.1 The complex plane We already discussed complex numbers briefly in Section 1.3.5.
    [Show full text]
  • Harmonic and Real Analysis Herbert Koch Universit¨Atbonn Wintersemester 2014-2015
    Notes for Harmonic and real analysis Herbert Koch Universit¨atBonn Wintersemester 2014-2015 Recommended literature: [10, 7, 14, 13, 15] 1 Contents Chapter 1. Introduction 5 1. An example 5 Chapter 2. Fourier series 9 1. Definitions 9 2. The convolution and Young's inequality 10 3. Lp convergence of partial sums 15 4. Regularity and Fourier series 16 5. Complex Interpolation 17 6. The Hardy-Littlewood maximal function and real interpolation 21 7. Higher dimensions 30 Chapter 3. Harmonic functions and the Poisson kernel 31 1. Basic properties and the Poisson kernel 31 2. Boundary behaviour of harmonic functions 36 3. Almost everywhere convergence 37 4. Subharmonic functions 38 5. The theorems of Riesz 40 6. Conjugate functions 41 7. The Hilbert transform 43 n Chapter 4. The Fourier transform on R 47 1. Definition and first properties 47 2. The Fourier transform of Schwartz functions 48 3. The Fourier transform on tempered distributions 51 4. Oscillatory integrals and stationary phase 58 5. Quantization 64 6. Cotlar's Lemma and the Theorem of Calder´on-Vaillancourt 70 Chapter 5. Singular integrals of Calder´on-Zygmund type 75 1. The setting 75 2. The Calder´on-Zygmund theorem 77 3. Examples 81 Chapter 6. Hardy and BMO 87 1. More general maximal functions and the Hardy space Hp 87 2. The atomic decomposition 91 3. Duality and BMO 95 4. Relation to harmonic functions 98 5. Div-curl type results 100 6. Application to elliptic PDEs 101 7. Pointwise estimates and perturbations of elliptic equations 107 3 4 CONTENTS Chapter 7.
    [Show full text]
  • Harmonic Functions
    6 HARMONIC FUNCTIONS This chapter is devoted to the Laplace equation. We introduce two of its important properties, the maximum principle and the rotational invariance. Then we solve the equation in series form in rectangles, circles, and related shapes. The case of a circle leads to the beautiful Poisson formula. 6.1 LAPLACE’S EQUATION If a diffusion or wave process is stationary (independent of time), then ut 0 ≡ and utt 0. Therefore, both the diffusion and the wave equations reduce to the Laplace≡ equation: uxx 0 in one dimension = u u uxx u yy 0 in two dimensions ∇·∇ = = + = u u uxx u yy uzz 0 in three dimensions ∇·∇ = = + + = A solution of the Laplace equation is called a harmonic function. In one dimension, we have simply uxx 0, so the only harmonic functions in one dimension are u(x) A Bx. But=this is so simple that it hardly gives us a clue to what happens =in higher+ dimensions. The inhomogeneous version of Laplace’s equation u f (1) = with f a given function, is called Poisson’s equation. Besides stationary diffusions and waves, some other instances of Laplace’s and Poisson’s equations include the following. 152 6.1 LAPLACE’S EQUATION 153 1. Electrostatics. From Maxwell’s equations, one has curl E 0 and div E 4πρ, where ρ is the charge density. The first equation implies= E grad= φ for a scalar function φ (called the electric potential). Therefore,=− φ div(grad φ) div E 4πρ , = =− =− which is Poisson’s equation (with f 4πρ ). =− 2. Steady fluid flow.
    [Show full text]
  • 3 Laplace's Equation
    3 Laplace’s Equation We now turn to studying Laplace’s equation ∆u = 0 and its inhomogeneous version, Poisson’s equation, ¡∆u = f: We say a function u satisfying Laplace’s equation is a harmonic function. 3.1 The Fundamental Solution Consider Laplace’s equation in Rn, ∆u = 0 x 2 Rn: Clearly, there are a lot of functions u which satisfy this equation. In particular, any constant function is harmonic. In addition, any function of the form u(x) = a1x1 +:::+anxn for constants ai is also a solution. Of course, we can list a number of others. Here, however, we are interested in finding a particular solution of Laplace’s equation which will allow us to solve Poisson’s equation. Given the symmetric nature of Laplace’s equation, we look for a radial solution. That is, we look for a harmonic function u on Rn such that u(x) = v(jxj). In addition, to being a natural choice due to the symmetry of Laplace’s equation, radial solutions are natural to look for because they reduce a PDE to an ODE, which is generally easier to solve. Therefore, we look for a radial solution. If u(x) = v(jxj), then x u = i v0(jxj) jxj 6= 0; xi jxj which implies 1 x2 x2 u = v0(jxj) ¡ i v0(jxj) + i v00(jxj) jxj 6= 0: xixi jxj jxj3 jxj2 Therefore, n ¡ 1 ∆u = v0(jxj) + v00(jxj): jxj Letting r = jxj, we see that u(x) = v(jxj) is a radial solution of Laplace’s equation implies v satisfies n ¡ 1 v0(r) + v00(r) = 0: r Therefore, 1 ¡ n v00 = v0 r v00 1 ¡ n =) = v0 r =) ln v0 = (1 ¡ n) ln r + C C =) v0(r) = ; rn¡1 1 which implies ½ c1 ln r + c2 n = 2 v(r) = c1 (2¡n)rn¡2 + c2 n ¸ 3: From these calculations, we see that for any constants c1; c2, the function ½ c1 ln jxj + c2 n = 2 u(x) ´ c1 (3.1) (2¡n)jxjn¡2 + c2 n ¸ 3: for x 2 Rn, jxj 6= 0 is a solution of Laplace’s equation in Rn ¡ f0g.
    [Show full text]
  • Applications of the Cauchy Theory
    Chapter 4 Applications Of The Cauchy Theory This chapter contains several applications of the material developed in Chapter 3. In the first section, we will describe the possible behavior of an analytic function near a singularity of that function. 4.1 Singularities We will say that f has an isolated singularity at z0 if f is analytic on D(z0,r) \{z0} for some r. What, if anything, can be said about the behavior of f near z0? The basic tool needed to answer this question is the Laurent series, an expansion of f(z)in powers of z − z0 in which negative as well as positive powers of z − z0 may appear. In fact, the number of negative powers in this expansion is the key to determining how f behaves near z0. From now on, the punctured disk D(z0,r) \{z0} will be denoted by D (z0,r). We will need a consequence of Cauchy’s integral formula. 4.1.1 Theorem Let f be analytic on an open set Ω containing the annulus {z : r1 ≤|z − z0|≤r2}, 0 <r1 <r2 < ∞, and let γ1 and γ2 denote the positively oriented inner and outer boundaries of the annulus. Then for r1 < |z − z0| <r2, we have 1 f(w) − 1 f(w) f(z)= − dw − dw. 2πi γ2 w z 2πi γ1 w z Proof. Apply Cauchy’s integral formula [part (ii)of (3.3.1)]to the cycle γ2 − γ1. ♣ 1 2 CHAPTER 4. APPLICATIONS OF THE CAUCHY THEORY 4.1.2 Definition For 0 ≤ s1 <s2 ≤ +∞ and z0 ∈ C, we will denote the open annulus {z : s1 < |z−z0| <s2} by A(z0,s1,s2).
    [Show full text]