SOLUTION OF POISSON S EQUATION. Contents
|
|
- Sharon Young
- 5 years ago
- Views:
Transcription
1 SOLUTION OF POISSON S EQUATION CRISTIAN E. GUTIÉRREZ OCTOBER 5, 2013 Contents 1. Differentiation under the integral sign 1 2. The Newtonian potential is C The Newtonian potential from the 3rd Green formula 4 4. Second derivatives 5 5. Solution in all space 9 6. The Newtonian potential is not necessarily 2nd order differentiable Helmholtz decomposition 12 References Differentiation under the integral sign A pretty general result that is often useful is the following. Theorem 1. Let kx, y) be a function defined for all x O R m open and for a.e. y such that kx, ) L 1 ) for all x O. Then a) The function ux) = kx, y) dy is well defined for all x O. b) Fix 1 i m. Assume that xi kx, y) exists for all x O and for almost all y, and there exists g L 1 ) such that xi kx, y) gy), x O, and for almost all y. Assume further that xi k, y) CO) for almost all y. 1
2 2 C. E. GUTIÉRREZ OCTOBER 5, 2013 Then the derivative xi u exists for all x O and xi ux) = xi kx, y) dy, for all x O. c) If the gradient D x k, y) CO) and there exists g L 1 ) such that D x kx, y) gy), x O, and for almost all y, then D x ux) = D x kx, y) dy, for all x O. d) If xj xi k, y) CO) for a.e. y, and there exists g L 1 ) such that then xj xi kx, y) gy), x O, and for almost all y, xj xi ux) = xj xi kx, y) dy. Proof. It follows immediately writing down the differential increment, using the mean value theorem, the continuity assumption and the Lebesgue dominated convergence theorem. 2. The Newtonian potential is C 1 Theorem 2. If ρ L 1 ) L ), then wx) = ρy)γx y) dy is well defined at each x and w C 1 ). Proof. Let ηt) be a smooth function such that ηt) = 0 for t 1, 0 η 1 everywhere, and ηt) = 1 for t 2. Let n 3, Γx) = 1 for n = 2, Γx) = log x ) x n 2 and consider the function Γ ɛ x) = Γx) η x /ɛ), where ɛ > 0. Step 1 If ρ L 1 ), then the function w ɛ x) = ρy)γ ɛ x y) dy is well defined for all x. This follows inmediately since Γ ɛ is bounded.
3 SOLUTION OF POISSON S EQUATION October 5, Step 2 If ρ L 1 ), then w ɛ C 1 ). First prove using the previous theorem that w ɛ has derivative at each point taking kx, y) = ρy)γ ɛ x y). Notice that the function Fx, y) = D x Γx y)η x y /ɛ )) Cɛ) for ɛ x y 2ɛ and Fx, y) = 0 otherwise. Therefore D x kx, y) Cɛ)ρy) for all x, y. Second we show that xi wx) is continuous in x for each ɛ fixed, i = 1,, n. Write j w ɛ x) j w ɛ z) = ρy) xj Γ ɛ x y) xj Γ ɛ z y) ) dy R n = ρy) D ) ) xj Γ ɛ tx + 1 t)z y) x z) dy. Differentiating we obtain as before that D ) xj Γ ɛ x) Cɛ) for all x, and therefore j w ɛ x) j w ɛ z) Cɛ) ρy) dy x z and the continuity of j w ɛ follows. Step 3 If ρ L 1 ) L ), then w ɛ x) wx) = ρy)γx y) dy uniformly for x, where wx) is well defined for all x. Clearly, wx) is well defined since Γ is integrable near the origin, Γ is bounded away from the origin, ρ is bounded and integrable. Write wx) w ɛ x) = ρy) Γx y) 1 η x y /ɛ) ) dy. So wx) w ɛ x) ρy) Γx y) dy C ρ ɛ 2. Step 4 If ρ L 1 ) L ), then Dw ɛ v uniformly in, where vx) = ρy)d x Γx y)) dy = v 1 x),, v n x)),
4 4 C. E. GUTIÉRREZ OCTOBER 5, 2013 where the last integral is well defined. It is clear that v is well defined. Write j w ɛ x) v j x) = ρy) xj Γ ɛ x y) xj Γx y) ) dy R n = ρy) ) xj Γx y)η x y /ɛ) + Γx y) xj η x y /ɛ) xj Γx y) ) dy R n = ρy) xj Γx y) η x y /ɛ) 1 ) ) + Γx y) ) xj η x y /ɛ) dy. So j w ɛ x) v j x) ρy) C n ρ xj Γx y) + Γx y) 1 ) ɛ η x y /ɛ) dy ) 1 x y dy = C ρ n 1 x y n 2 ɛ. ɛ Step 5 As a consequence of steps 2-4, the functions w and v are continous when ρ L 1 ) L ). Step 6 If ρ L 1 ) L ), then w C 1 ), and Dwx) = vx). By the fundamental theorem of calculus, w ɛ x 1, x 2,, x n ) = x1 and from the uniform convergence y 1 x1 w ɛ t, x 2,, x n ) dt + w ɛ y 1, x 2,, x n ), wx 1, x 2,, x n ) = x1 y 1 v 1 t, x 2,, x n ) dt + wy 1, x 2,, x n ). 3. The Newtonian potential from the 3rd Green formula 1 Let n 3 and Γx) = ω n 1 2 n) x 2 n, where ω n 1 is the surface area of the unit sphere in. The third Green formula reads Γx y) uy) = ux) Γx y) u ) ηx) η x) dσx) + Γx y) ux) dx, where u C 2 ) C 1 ), and u L 1 ) and the domain is sufficiently regular. Suppose u C 2 ) solves the Poisson equation u = ρ
5 SOLUTION OF POISSON S EQUATION October 5, in all space, with ρ L 1 ) L ). Applying the third Green formula in the ball B R y) yields Γx y) uy) = ux) Γx y) u ) B R y) ηx) η x) dσx) + Γx y) ρx) dx. B R y) Suppose in addition that ux) = O1/ x ), Dux) = O1/ x 2 ), when x. Then, when R, the surface integral tends to zero and we obtain uy) = Γx y) ρx) dx, for each y. 4. Second derivatives We assume f : R is locally Hölder continuous, that is, there exists 0 < α 1 such that for each K compact there is a constant C K > 0 such that f x) f y) C K x y α, x, y K. We also assume f is bounded in. The goal is to prove the representation formula 1) D ij wx) = D ij Γ)x y) f y) f x)) dy f x) 0 D i Γx y)ν j y) dσy), 0 for the second derivatives of the Newtonian potential, where 0 is any domain for which the divergence theorem holds, with 0 ; and f x) = f x) for x and f x) = 0 in 0 \. Assume n 3, let Γx) = C n x 2 n be the fundamental solution, η : R R is smooth, ηt) = 0 for t 1, 0 η 1, ηt) = 1 for t 2 and 0 η 2. Define for x ux) = D ij Γ)x y) f y) f x)) dy f x) D i Γx y)ν j y) dσy). 0 0 Since D ij Γx) C n x n, and f is bounded and locally Hölder continuous, then the function ux) is well defined for all x. Let ɛ > 0 and define v ɛ x) = D i Γx y)ηx y)/ɛ) f y) dy,
6 6 C. E. GUTIÉRREZ OCTOBER 5, 2013 and let vx) = D i wx), where wx) = Γx y) f y) dy. If f is integrable in, then, as in Theorem 2, v ɛ C 1 ), and if in addition f is bounded, w C 1 ). We have D j v ɛ x) = D xj Di Γx y)ηx y)/ɛ) ) f y) dy = D xj Di Γx y)ηx y)/ɛ) ) f y) dy 0 = D xj Di Γx y)ηx y)/ɛ) ) f y) f x)) dy 0 + f x) D xj Di Γx y)ηx y)/ɛ) ) dy 0 = D xj Di Γx y)ηx y)/ɛ) ) f y) f x)) dy 0 f x) D yj Di Γx y)ηx y)/ɛ) ) dy 0 = D xj Di Γx y)ηx y)/ɛ) ) f y) f x)) dy 0 f x) D i Γx y)ηx y)/ɛ)ν j y) dσy) 0 from the divergence theorem. Since x, if we take ɛ distx, 0 )/2, then x y 2ɛ for y 0 and so ηx y)/ɛ) = 1. So D j v ɛ x) = D xj Di Γx y)ηx y)/ɛ) ) f y) f x)) dy 0 f x) D i Γx y)ν j y) dσy). 0
7 SOLUTION OF POISSON S EQUATION October 5, Then subtracting we get ux) D j v ɛ x) = Dij Γx y) D xj Di Γx y)ηx y)/ɛ) )) f y) f x)) dy 0 [ ) = D xj 1 ηx y)/ɛ) Di Γx y) ] f y) f x)) dy 0 [ ) = D xj 1 ηx y)/ɛ) Di Γx y) ] f y) f x)) dy [ = D ij Γx y) 1 ηx y)/ɛ) ) D i Γx y) 1 ɛ η x y /ɛ) x ] j y j f x y y) f x)) dy. Let K be compact, and K = {y : disty, K) distk, )/2}. We have K is compact, K, and if ɛ < distk, )/2, then B ɛ x) K for all x K. Therefore estimating the last integral we obtain for x K [ DijΓx ux) D j v ɛ x) y) Di + Γx y) ] C f y) f x) dy ɛ C K [ C x y + C C n x y n 1 ɛ ] x y α dy C K ɛ α. Therefore D j v ɛ u uniformly on compact subsets of as ɛ 0, so u is continuous in. In addition, v ɛ v uniformly in this follows directly by subtracting). This implies that w C 2 ) and ux) = D ij wx) for x by the fundamental theorem of calculus, because we write v ɛ x 1,, x j 1, y j, x j+1,, x n ) = yj x j D j v ɛ x 1,, x j 1, t, x j+1,, x n ) dt + v ɛ x), and pass to the limit as ɛ 0. So we have proved the following representation formula for the second derivatives of the Newtonian potential w for x : 2) D ij wx) = D ij Γ)x y) f y) f x)) dy f x) D i Γx y)ν j y) dσy), 0 0 under the assumption that f is bounded in and locally Hölder continuous for some 0 < α 1. We next prove that w = f in. Let x, take R sufficiently large such that B R x), and apply 2) with 0 = B R x), then D ii wx) = D ij Γ)x y) f y) f x)) dy f x) x y <R D i Γx y)ν j y) dσy), x y =R
8 8 C. E. GUTIÉRREZ OCTOBER 5, 2013 adding over 1 i n and using that Γx) = 0 for x 0, we obtain wx) = f x) DΓx y) νy) dσy). x y =R We have DΓx) = C n 2 n) x y x and νy) = x n x y. Then wx) = f x)c n2 n)ω n, 1, we are 2 n)ω n where ω n is the surface area of the unit sphere in. Since C n = done. We also have that wx) = 0 for x, which follows differentiating twice under the integral sign the justification follows from Theorem 1). So we have proved the following theorem. Theorem 3. Let be a bounded open set, ρ be a bounded function in. Then 1) w C 1 ); 2) w C 2 \ ); 3) w = 0 in \ ; 4) if in addition ρ is Hölder continuous, then w C 2 ) and w = ρ in. 5) the second derivatives of w might not exist in. Actually, the last item was not proved and we leave it as an exercise. For that consider in R 3 the ball B0, R). Assume that the ball is homogeneous with density λ. The Newtonian potential due to this ball is given by 1 Ux) = G λ B0,R) y x dy. a) Show that Gλ 4 1 Ux) = 3 πr3 x, for x > R Gλ 2 3 π3r2 x 2 ), for 0 x R. HINT: to simplify the calculation notice that UOx) = Ux) for any rotation O around the origin, that is, U is a radial function. b) Prove that the first order derivatives of U exist and are continuous everywhere. c) Prove that all the second derivatives of U exist and are continuous in R 3 \ {x : x = R}. Prove that they do not exist on x = R.
9 SOLUTION OF POISSON S EQUATION October 5, d) Prove that Gλ4π, for 0 x < R Ux) = div DUx) = 0, for x > R. 5. Solution in all space Let wx) = Γx y)ρy) dy. We shall prove that under some conditions on ρ, we have that w is C 2 everywhere and wx) = ρx) for all x. Fix R > 0 and let φ C ) be such that φy) = 1 for y 2R, φy) = 0 for y 3R and 0 φ 1. Write wx) = Γx y)ρy)φy) dy + Γx y)ρy)1 φy)) dy We have := w 1 x) + w 2 x). w 1 x) = Γx y)ρy)φy) dy B 3R 0) and from what we proved before we have that if ρ is locally Hölder continuous in, then w 1 C 2 B 3R 0)) and w 1 x) = ρx)φx) for all x B 3R 0). Therefore w 1 x) = ρx) for x B 2R 0). We shall prove that w 2 x) = 0 for all x B R 0). If ρ L 1 ) L ), then from Theorem 2 we have that w 2 C 1 ) and 3) Dw 2 x) = D x Γx y)ρy)1 φy)) dy. Let x 0 B R 0), then B R/2 x 0 ) B 3R/2 0). Notice that the integrand in w 2 vanishes for y 2R, and therefore the domain of integration reduces to the set y 2R. We want to differentiate under the integral sign in 3), and for that we need to verify the hypotheses in Theorem 1. Indeed, we have that 2 C x j x j Γx y) x y, n for x y, and since x y R/2 for x B R x 0 ) and y 2R xj x j Γx y) 2n C there. So if we let gy) = 2n C ρy)1 φy)), then we obtain that Rn hx, y) := xj x j Γx y)ρy)1 φy)) gy)
10 10 C. E. GUTIÉRREZ OCTOBER 5, 2013 for all x B R/2 x 0 ) and for all y. In addition, h, y) CB R/2 x 0 )) for all y. Therefore, we can differentiate under the integral sign and we obtain xj x j w 2 x) = xj x j Γx y)ρy)1 φy)) dy, for all x B R/2 x 0 ) and so w 2 x) = Γx y)ρy)1 φy)) dy = 0, for all x B R/2 x 0 ), for each x 0 B R 0). Since R is arbitrary, we therefore have proved the following theorem: Theorem 4. If ρ defined in is locally Hölder continuous of order α, for some 0 < α 1, and ρ L 1 ) L ), then the Newtonian potential wx) = Γx y)ρy) dy is C 2 in all space, and satisfies the Poisson equation wx) = ρx) in. 6. The Newtonian potential is not necessarily 2nd order differentiable We show that there exists f continuous in B 2 0) such that wx) = Γx y) f y) dy B 2 0) is not second order differentiable at the origin. The example is from [GT83, Exercise 4.9]. Assume n 3. Let c k be a sequence such that c k 0 and k and c k = +. Let η C 0 B 20)) with ηx) = 1 for x < 1, and let P be a harmonic polynomial of degree two in with D α P 0 for some α = 2, e.g., P = x 1 x 2. Let f x) = c k ηp)2 k x). Given x 0 in B 2 0) there exists a unique integer N 0 such that 2 N x < 2 1 N. Notice that ηp)x) = 0 for x > 2 and for x < 1. We have f x) = = c N ηp)2 N x). {k:2 k x <2 1 k }
11 SOLUTION OF POISSON S EQUATION October 5, This function is clearly continuous in B 2 0). We write wx) = = Γx y) f y) dy 2 k y <2 1 k c k Γx y) ηp)2 2 k y) dy = k y <2 1 k c k I k x). Changing variables we have I k x) = 2 2k 1 z <2 Γ2 k x z) ηp)z) dz = 2 2k Γ2 k x z) ηp)z) dz and from the 3rd Green formula since ηp has compact support we get as in Section 3 that Γ2 k x z) ηp)z) dz = ηp)2 k x). Therefore wx) = c k 2 2k ηp)2 k x). Given x < 2 there exists a unique N 0 integer such that 2 N x < 2 1 N. So 2 k x 2 k N 2 for k N + 1, and so ηp)2 k x) = 0 for k N + 1. On the other hand, 2 k x < 2 k 2 1 N 1 if k N 1 and hence η2 k x) = 1 for k N 1. Therefore Thus N 1 wx) = c k 2 2k P2 k x) + c N 2 2N ηp)2 N x), for 2 N x < 2 1 N. N 1 D α wx) = c k D α P2 k x) + c N D α ηp)2 N x), α = 2, for 2 N x < 2 1 N. If P = x 1 x 2, then N 1 D 12 wx) = c k + c N D 12 ηp)2 N x), for 2 N x < 2 1 N. If x 0, then N and therefore D 12 wx) +.
12 12 C. E. GUTIÉRREZ OCTOBER 5, Helmholtz decomposition Let F be a vector field in R 3 we shall prove that under appropriate hypotheses F = R + S with div S = 0 and curl R = 0. This means that F = Du + curl H for some function u and some field H. Recall that if G is a vector field with div G = 0 in R 3, then there exists a field H such that G = curl H. Also, if G is a field with curl G = 0, then there exists a function u such that G = Du. In addition, if G = curl H 1 = curl H 2, then H 1 H 2 = Dφ for some function φ. Assuming the desired decomposition, we have div F = div S + div R = div R. Since curl R = 0, there exists u with R = Du, and so div F = div Du = u. Therefore, if div F satisfies the hypotheses of Theorem 4 and let u be the Newtonian potential of div F, then u = div F. So if we let R = Du, and S = F R, then curl R = 0, and div S = div F div R = div F div Du = 0. So we have: Theorem 5 Helmholtz decomposition). Let F : R 3 R 3 be a field such that div F is locally Hölder continuous in R 3, uniformly bounded and integrable over all space. If wx) = Γx y) div Fy) dy, R 3 then Fx) = Dwx) + curl Hx) for some field H. The field H is determined up to the gradient of an arbitrary function φ. If H is a field we define the Laplacian of F by H = Ddiv H) curl curl H). If H = a 1 x)i + a 2 x)j + a 3 x)k, then is easy to see calculating the components of H that H = a 1 x)i + a 2 x)j + a 3 x)k. If we assume div H = 0, then H = curl curl H) = curl F by the Helmholtz decomposition. But the equation H = curl F is a Poisson equation now between
13 SOLUTION OF POISSON S EQUATION October 5, vectors, that is, three scalar Poisson equations a i = curl F) i, i = 1, 2, 3. Solving each one as before with the Newtonian potential we get that Hx) = Γx y) curl Fy) dy, R 3 that is, H is the Newtonian potential now a vector) of curl F. This means that if for a given field F we know div F and curl F, then F is automatically determined by the formula ) ) Fx) = D Γx y) div Fy) dy + curl Γx y) curl Fy) dy. R 3 R 3 References [GT83] D. Gilbarg and N. S. Trudinger. Elliptic Partial Differential Equations of Second Order. Springer-Verlag, New York, 2nd edition, Department of Mathematics, Temple University, Philadelphia, PA address: gutierre@temple.edu
THE GREEN FUNCTION. Contents
THE GREEN FUNCTION CRISTIAN E. GUTIÉRREZ NOVEMBER 5, 203 Contents. Third Green s formula 2. The Green function 2.. Estimates of the Green function near the pole 2 2.2. Symmetry of the Green function 3
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationNOTES ON SCHAUDER ESTIMATES. r 2 x y 2
NOTES ON SCHAUDER ESTIMATES CRISTIAN E GUTIÉRREZ JULY 26, 2005 Lemma 1 If u f in B r y), then ux) u + r2 x y 2 B r y) B r y) f, x B r y) Proof Let gx) = ux) Br y) u r2 x y 2 Br y) f We have g = u + Br
More informationSOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD. 1. Dirichlet integral
SOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD CRISTIAN E. GUTIÉRREZ FEBRUARY 3, 29. Dirichlet integral Let f C( ) with open and bounded. Let H = {u C ( ) : u = f on } and D(u) = Du(x) 2 dx.
More informationELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS)
ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS) MATANIA BEN-ARTZI. BOOKS [CH] R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol. Interscience Publ. 962. II, [E] L.
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationCHAPTER 2. Laplace s equation
18 CHAPTER 2 Laplace s equation There can be but one option as to the beauty and utility of this analysis by Laplace; but the manner in which it has hitherto been presented has seemed repulsive to the
More informationElliptic PDEs of 2nd Order, Gilbarg and Trudinger
Elliptic PDEs of 2nd Order, Gilbarg and Trudinger Chapter 2 Laplace Equation Yung-Hsiang Huang 207.07.07. Mimic the proof for Theorem 3.. 2. Proof. I think we should assume u C 2 (Ω Γ). Let W be an open
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
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 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 informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationElementary Theory and Methods for Elliptic Partial Differential Equations. John Villavert
Elementary Theory and Methods for Elliptic Partial Differential Equations John Villavert Contents 1 Introduction and Basic Theory 4 1.1 Harmonic Functions............................... 5 1.1.1 Mean Value
More informationMaxwell s equations for electrostatics
Maxwell s equations for electrostatics October 6, 5 The differential form of Gauss s law Starting from the integral form of Gauss s law, we treat the charge as a continuous distribution, ρ x. Then, letting
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 information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More information4 Divergence theorem and its consequences
Tel Aviv University, 205/6 Analysis-IV 65 4 Divergence theorem and its consequences 4a Divergence and flux................. 65 4b Piecewise smooth case............... 67 4c Divergence of gradient: Laplacian........
More information2. Function spaces and approximation
2.1 2. Function spaces and approximation 2.1. The space of test functions. Notation and prerequisites are collected in Appendix A. Let Ω be an open subset of R n. The space C0 (Ω), consisting of the C
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationNonlinear aspects of Calderón-Zygmund theory
Ancona, June 7 2011 Overture: The standard CZ theory Consider the model case u = f in R n Overture: The standard CZ theory Consider the model case u = f in R n Then f L q implies D 2 u L q 1 < q < with
More informationTHE DIRICHLET PROBLEM FOR THE LAPLACE OPERATOR
THE DIRICHLET PROBLEM FOR THE LAPLACE OPERATOR Stefano Meda Università di Milano-Bicocca c Stefano Meda 2013 ii A Francesco Contents I The Dirichlet problem via Perron s method 1 1 The classical Dirichlet
More informationMa 1c Practical - Solutions to Homework Set 7
Ma 1c Practical - olutions to omework et 7 All exercises are from the Vector Calculus text, Marsden and Tromba (Fifth Edition) Exercise 7.4.. Find the area of the portion of the unit sphere that is cut
More informationPractice Problems for Exam 3 (Solutions) 1. Let F(x, y) = xyi+(y 3x)j, and let C be the curve r(t) = ti+(3t t 2 )j for 0 t 2. Compute F dr.
1. Let F(x, y) xyi+(y 3x)j, and let be the curve r(t) ti+(3t t 2 )j for t 2. ompute F dr. Solution. F dr b a 2 2 F(r(t)) r (t) dt t(3t t 2 ), 3t t 2 3t 1, 3 2t dt t 3 dt 1 2 4 t4 4. 2. Evaluate the line
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationRegularity of the p-poisson equation in the plane
Regularity of the p-poisson equation in the plane Erik Lindgren Peter Lindqvist Department of Mathematical Sciences Norwegian University of Science and Technology NO-7491 Trondheim, Norway Abstract We
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationAsymptotic behavior of infinity harmonic functions near an isolated singularity
Asymptotic behavior of infinity harmonic functions near an isolated singularity Ovidiu Savin, Changyou Wang, Yifeng Yu Abstract In this paper, we prove if n 2 x 0 is an isolated singularity of a nonegative
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationPartial Differential Equations. William G. Faris
Partial Differential Equations William G. Faris May 17, 1999 Contents 1 The Laplace equation 5 1.1 Gradient and divergence....................... 5 1.2 Equilibrium conservation laws....................
More informationAsymptotic Behavior of Infinity Harmonic Functions Near an Isolated Singularity
Savin, O., and C. Wang. (2008) Asymptotic Behavior of Infinity Harmonic Functions, International Mathematics Research Notices, Vol. 2008, Article ID rnm163, 23 pages. doi:10.1093/imrn/rnm163 Asymptotic
More informationMINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA
MINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA SPENCER HUGHES In these notes we prove that for any given smooth function on the boundary of
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
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 informationENGI Gradient, Divergence, Curl Page 5.01
ENGI 94 5. - Gradient, Divergence, Curl Page 5. 5. The Gradient Operator A brief review is provided here for the gradient operator in both Cartesian and orthogonal non-cartesian coordinate systems. Sections
More informationThe Mollifier Theorem
The Mollifier Theorem Definition of the Mollifier The function Tx Kexp 1 if x 1 2 1 x 0 if x 1, x R n where the constant K is chosen such that R n Txdx 1, is a test function on Rn. Note that Tx vanishes,
More informationHarmonic Functions and Brownian motion
Harmonic Functions and Brownian motion Steven P. Lalley April 25, 211 1 Dynkin s Formula Denote by W t = (W 1 t, W 2 t,..., W d t ) a standard d dimensional Wiener process on (Ω, F, P ), and let F = (F
More informationg(x) = P (y) Proof. This is true for n = 0. Assume by the inductive hypothesis that g (n) (0) = 0 for some n. Compute g (n) (h) g (n) (0)
Mollifiers and Smooth Functions We say a function f from C is C (or simply smooth) if all its derivatives to every order exist at every point of. For f : C, we say f is C if all partial derivatives to
More informationClass Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem
MATH 8.52 COUSE NOTES - CLASS MEETING # 9 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem. epresentation Formula
More informationEverywhere differentiability of infinity harmonic functions
Everywhere differentiability of infinity harmonic functions Lawrence C. Evans and Charles K. Smart Department of Mathematics University of California, Berkeley Abstract We show that an infinity harmonic
More informationQuantitative Homogenization of Elliptic Operators with Periodic Coefficients
Quantitative Homogenization of Elliptic Operators with Periodic Coefficients Zhongwei Shen Abstract. These lecture notes introduce the quantitative homogenization theory for elliptic partial differential
More informationA PROPERTY OF SOBOLEV SPACES ON COMPLETE RIEMANNIAN MANIFOLDS
Electronic Journal of Differential Equations, Vol. 2005(2005), No.??, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A PROPERTY
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationNotes on Partial Differential Equations. John K. Hunter. Department of Mathematics, University of California at Davis
Notes on Partial Differential Equations John K. Hunter Department of Mathematics, University of California at Davis Contents Chapter 1. Preliminaries 1 1.1. Euclidean space 1 1.2. Spaces of continuous
More informationApplications of the Maximum Principle
Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is
More informationJim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt
Jim Lambers MAT 28 ummer emester 212-1 Practice Final Exam olution 1. Evaluate the line integral xy dx + e y dy + xz dz, where is given by r(t) t 4, t 2, t, t 1. olution From r (t) 4t, 2t, t 2, we obtain
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 informationTHE DIRICHLET PROBLEM WITH BM O BOUNDARY DATA AND ALMOST-REAL COEFFICIENTS
THE DIRICHLET PROBLEM WITH BM O BOUNDARY DATA AND ALMOST-REAL COEFFICIENTS ARIEL BARTON Abstract. It is known that a function, harmonic in a Lipschitz domain, is the Poisson extension of a BMO function
More information18.02 Multivariable Calculus Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.02 Multivariable Calculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. V7. Laplace's Equation and
More informationCorrections to the First Printing: Chapter 6. This list includes corrections and clarifications through November 6, 1999.
June 2,2 1 Corrections to the First Printing: Chapter 6 This list includes corrections and clarifications through November 6, 1999. We should have mentioned that our treatment of differential forms, especially
More informationANALYSIS QUALIFYING EXAM FALL 2017: SOLUTIONS. 1 cos(nx) lim. n 2 x 2. g n (x) = 1 cos(nx) n 2 x 2. x 2.
ANALYSIS QUALIFYING EXAM FALL 27: SOLUTIONS Problem. Determine, with justification, the it cos(nx) n 2 x 2 dx. Solution. For an integer n >, define g n : (, ) R by Also define g : (, ) R by g(x) = g n
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More informationwhere F denoting the force acting on V through V, ν is the unit outnormal on V. Newton s law says (assume the mass is 1) that
Chapter 5 The Wave Equation In this chapter we investigate the wave equation 5.) u tt u = and the nonhomogeneous wave equation 5.) u tt u = fx, t) subject to appropriate initial and boundary conditions.
More information< 1. max x B(0,1) f. ν ds(y) Use Poisson s formula for the ball to prove. (r x ) x y n ds(y) (x B0 (0, r)). 1 nα(n)r n 1
7 On the othe hand, u x solves { u n in U u on U, so which implies that x G(x, y) x n ng(x, y) < n (x B(, )). Theefoe u(x) fo all x B(, ). G ν ds(y) + max g G x ν ds(y) + C( max g + max f ) f(y)g(x, y)
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More information14 Divergence, flux, Laplacian
Tel Aviv University, 204/5 Analysis-III,IV 240 4 Divergence, flux, Laplacian 4a What is the problem................ 240 4b Integral of derivative (again)........... 242 4c Divergence and flux.................
More informationMollifiers and its applications in L p (Ω) space
Mollifiers and its alications in L () sace MA Shiqi Deartment of Mathematics, Hong Kong Batist University November 19, 2016 Abstract This note gives definition of mollifier and mollification. We illustrate
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationMath 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C
Math 35 Solutions for Final Exam Page Problem. ( points) (a) ompute the line integral F ds for the path c(t) = (t 2, t 3, t) with t and the vector field F (x, y, z) = xi + zj + xk. (b) ompute the line
More informationHARMONIC FUNCTIONS TSOGTGEREL GANTUMUR
HARMONIC FUNCTIONS TSOGTGEREL GANTUMUR Abstract. In these notes, we explore the fundamental properties of harmonic functions, by using relatively elementary methods. Contents. Introduction 2. Green s identities
More informationSecond Order ODEs. Second Order ODEs. In general second order ODEs contain terms involving y, dy But here only consider equations of the form
Second Order ODEs Second Order ODEs In general second order ODEs contain terms involving y, dy But here only consider equations of the form A d2 y dx 2 + B dy dx + Cy = 0 dx, d2 y dx 2 and F(x). where
More informationMethods on Nonlinear Elliptic Equations
Wenxiong Chen and Congming Li Methods on Nonlinear Elliptic Equations SPIN AIMS internal project number, if known Monograph September 22, 2008 AIMS Preface In this book we intend to present basic materials
More informationPartial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces
Partial Differential Equations, nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Shih-Hsin Chen, Yung-Hsiang Huang 7.8.3 Abstract In these exercises always denote an open set of with smooth boundary. As
More informationComplex Analysis, Stein and Shakarchi Meromorphic Functions and the Logarithm
Complex Analysis, Stein and Shakarchi Chapter 3 Meromorphic Functions and the Logarithm Yung-Hsiang Huang 217.11.5 Exercises 1. From the identity sin πz = eiπz e iπz 2i, it s easy to show its zeros are
More informationTHE HARDY LITTLEWOOD MAXIMAL FUNCTION OF A SOBOLEV FUNCTION. Juha Kinnunen. 1 f(y) dy, B(x, r) B(x,r)
Appeared in Israel J. Math. 00 (997), 7 24 THE HARDY LITTLEWOOD MAXIMAL FUNCTION OF A SOBOLEV FUNCTION Juha Kinnunen Abstract. We prove that the Hardy Littlewood maximal operator is bounded in the Sobolev
More informationLecture notes: Introduction to Partial Differential Equations
Lecture notes: Introduction to Partial Differential Equations Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 Classification of Partial Differential
More informationMULTIPLE SOLUTIONS FOR CRITICAL ELLIPTIC PROBLEMS WITH FRACTIONAL LAPLACIAN
Electronic Journal of Differential Equations, Vol. 016 (016), No. 97, pp. 1 11. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu MULTIPLE SOLUTIONS
More informationSuggested Solution to Assignment 7
MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationPERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 28, 2003, 207 222 PERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS Fumi-Yuki Maeda and Takayori Ono Hiroshima Institute of Technology, Miyake,
More informationAdvanced Calculus: MATH 410 Professor David Levermore 28 November 2006
Advanced Calculus: MATH 410 Professor David Levermore 28 November 2006 1. Uniform Continuity Uniform continuity is a very useful concept. Here we introduce it in the context of real-valued functions with
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 informationSobolev Spaces. Chapter Hölder spaces
Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect
More informationELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia
ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia Abstract This paper is concerned with the study of scattering of
More informationRecall that if X is a compact metric space, C(X), the space of continuous (real-valued) functions on X, is a Banach space with the norm
Chapter 13 Radon Measures Recall that if X is a compact metric space, C(X), the space of continuous (real-valued) functions on X, is a Banach space with the norm (13.1) f = sup x X f(x). We want to identify
More informationMATH 228: Calculus III (FALL 2016) Sample Problems for FINAL EXAM SOLUTIONS
MATH 228: Calculus III (FALL 216) Sample Problems for FINAL EXAM SOLUTIONS MATH 228 Page 2 Problem 1. (2pts) Evaluate the line integral C xy dx + (x + y) dy along the parabola y x2 from ( 1, 1) to (2,
More informationMath 220A Homework 4 Solutions
Math 220A Homework 4 Solutions Jim Agler 26. (# pg. 73 Conway). Prove the assertion made in Proposition 2. (pg. 68) that g is continuous. Solution. We wish to show that if g : [a, b] [c, d] C is a continuous
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationComparison Results for Semilinear Elliptic Equations in Equimeasurable Domains
Comparison Results for Semilinear Elliptic Equations in Equimeasurable Domains François HAMEL Aix-Marseille University & Institut Universitaire de France In collaboration with Emmanuel RUSS Workshop on
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationRegularity of Weak Solution to Parabolic Fractional p-laplacian
Regularity of Weak Solution to Parabolic Fractional p-laplacian Lan Tang at BCAM Seminar July 18th, 2012 Table of contents 1 1. Introduction 1.1. Background 1.2. Some Classical Results for Local Case 2
More information4.1. Dirichlet, Poisson and Neumann boundary value problems
Chapter 4 Laplace s equation 4.. Dirichlet, Poisson and Neumann boundary value problems The most commonly occurring form of problem that is associated with Laplace s equation is a boundary value problem,
More informationTHE LAPLACIAN AND THE DIRICHLET PROBLEM
THE LAPLACIAN AND THE DIRICHLET PROBLEM Stefano Meda Università di Milano-Bicocca c Stefano Meda 2016 ii A Francesco Contents 0.1 Background and preliminary results................ 1 0.1.1 Convolution..........................
More informationMATH 425, HOMEWORK 3 SOLUTIONS
MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the
More informationVISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS
VISCOSITY SOLUTIONS OF ELLIPTIC EQUATIONS LUIS SILVESTRE These are the notes from the summer course given in the Second Chicago Summer School In Analysis, in June 2015. We introduce the notion of viscosity
More informationMath Exam IV - Fall 2011
Math 233 - Exam IV - Fall 2011 December 15, 2011 - Renato Feres NAME: STUDENT ID NUMBER: General instructions: This exam has 16 questions, each worth the same amount. Check that no pages are missing and
More informationVector Calculus, Maths II
Section A Vector Calculus, Maths II REVISION (VECTORS) 1. Position vector of a point P(x, y, z) is given as + y and its magnitude by 2. The scalar components of a vector are its direction ratios, and represent
More informationGRAND SOBOLEV SPACES AND THEIR APPLICATIONS TO VARIATIONAL PROBLEMS
LE MATEMATICHE Vol. LI (1996) Fasc. II, pp. 335347 GRAND SOBOLEV SPACES AND THEIR APPLICATIONS TO VARIATIONAL PROBLEMS CARLO SBORDONE Dedicated to Professor Francesco Guglielmino on his 7th birthday W
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More information