In this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
|
|
- Elvin Pope
- 5 years ago
- Views:
Transcription
1 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 applications of elliptic PDEs. The steady-state of the heat equation solves Laplace s equation u t κ 2 u = 0, x, t > 0 u(x, 0) = f(x), x, u(x, t) = g(x), x, t > 0, κ 2 u = 0, x, u(x) = g(x), x. If we think of ρ(x) as a charge density, the potential, u, associated with the generated electric field satisfies Poisson s Equation 2 V = ρ(x) ǫ 0, x, V = f(x), x, where ǫ 0 is the permitivity of free space and the Electric field is given by E = V. Suppose Γ is a closed, non intersecting curve in R 3. The minimal surface with boundary Γ is the minimum of the functional J(u) := 1 + u 2 x + u2 y dxdy, 133
2 Elliptic PDEs where is the projection of the surface onto the xy plane and u is the height of the surface above the xy plane. Homework (8.1) will show the minimum satisfies Plateau s Equation (1 + u 2 y)u xx + (1 + u 2 x)u yy 2u x u y u xy = 0, u(x) = f(x) x. If u x and u y are small compared to one, Laplace s equation approximates Plateau s equation Bounded Domains Our spectral methods work for elliptic PDEs as well. Consider the PDE (approximating a soap film) on the square = {(x, y) 0 < x < π, 0 < y < π}: 2 u = 0, on u(0, y) = f(y), 0 < y < π, u(π, y) = g(y), 0 < y < π, u(x, 0) = 0, 0 < x < π, u(x, π) = 0, 0 < x < π. We find conditions on f and g so that the solution satisfies u C, () C(). Since we have homogeneous boundary conditions on the top can bottom of the square, we attempt to find a solution of the form u(x, y) = C n (x) sin(ny). i=1 The last two boundary conditions are satisfied. Putting this sum (formally) in the PDE, we get C n(x) n 2 C n sin(ny) = 0. i=1 Orthogonality implies C n n2 C n = 0 or C n (x) = c 1 e nx + c 2 e nx. Recall sinh(x) = 1 2 (enx e nx ) and cosh(x) = 1 2 (enx + e nx ). We may express the solution to C n n2 C n = 0 as C n (x) = c 1 e nx + c 2 e nx = c 1 sinh(nx) + c 2 cosh(nx) = c 1 sinh(nx) + c 2 sinh(n(π x)). The last version is the easiest to work with. Thus, so far, we have u(x, y) = a n sinh(nx) + b n sinh(n(π x) sin(ny). At x = 0 we need i=1 f(y) = b n sinh(nπ) sin(ny). i=1
3 8.2. The Fundamental Solution 135 Using orthogonality again, we find b n = Similarly, we find at x = π, a n = 2 π sinh(nπ) 2 π sinh(nπ) π 0 π 0 f(y) sin(πy) dy. g(y) sin(πy) dy. To get convergence we first apply Dirichlet s theorem. We need f, g to be continuously differentiable and their odd extension to be continuous. This requires f(0) = g(0) = f(π) = g(π) = 0. Next we note that, for large n, 1/ sinh(nπ) e nπ. When we put a n and b n back in the sum, we get terms like sinh(nx)/ sinh(nπ). For large n and 0 < x < π we see sinh(nx) sinh(nπ) en(x π). As in the parabolic case we see the Fourier coefficients converge exponentially fast, as long as x is away from the boundary, and so, u C, (). To get u(x, y) continuous on the closure of, we need both a n and b n to decay like 1/n 2 since the exponential decay goes away on the boundary. For this will will have to integrate both f and g by parts twice. This will require, in addition to the other conditions, f, g C 2 (0, π). Suppose we wanted to solve on the square = {(x, y) 0 < x < π, 0 < y < π}: 2 u = 0, on u(0, y) = f 1 (y), 0 < y < π, u(π, y) = f 2 (y), 0 < y < π, u(x, 0) = g 1 (x), 0 < x < π, u(x, π) = g 2 (x), 0 < x < π. We could do this by first solving, as above, the problems 2 u 1 = 0, 2 u 2 = 0, u 1 (0, y) = f 1 (y), 0 < y < π u 2 (0, y) = 0, 0 < y < π u 1 (π, y) = f 2 (y), 0 < y < π, u 2 (π, y) = 0, 0 < y < π, u 1 (x, 0) = 0, 0 < x < π, u 2 (x, 0) = g 1 (x), 0 < x < π, u 1 (x, π) = 0, 0 < x < π u 2 (x, π) = g 2 (x), 0 < x < π. The solution to our problem is then u = u 1 + u The Fundamental Solution We want to solve Poisson s equation 2 u = f with a variety of boundary conditions and domains - both bounded and unbounded. We follow our usual procedure for constructing a Green s function. Set Lu = 2 u. We need to compute the adjoint. For the moment let us suppose the domain
4 Elliptic PDEs is bounded subset of R N and normal (meaning the divergence theorem applies). Integrating by parts twice, we find v 2 u d = v u n dσ u v n dσ + u 2 v d, where dσ is the surface element on in R n 1 and single integral represents a multiple integral in R N. Rearranging we get Green s Second Identity: (8.1) (v 2 u u 2 v) d = v u n dσ u v n dσ. This shows L v = 2 v, and L is formally self adjoint. As before we will need to solve L v = δ(x x ). We have to know the boundary conditions for our problem before we can prescribe the boundary conditions on v. In general the solution will have the form, according to (8.1), with G(x, x ) = v(x, x ) u(x) = G(x, x )f(x ) d + u(x ) x G(x, x ) n dσ G(x, x ) x u(x ) n dσ, where x means the gradient with respect to the primed variables. Since 2 yg(x, y) = 0, x = y, we might expect the solution to be radially symmetric. This suggest trying to solve 2 u = 0 in polar coordinates. To do this, set r 2 = x x2 2 + x2 N, and u(x 1, x 2,... x N ) = U(r). Then u xi = U r xi r. Here xi r = xi r. Thus and We can easily solve u xix i = U rr x i r 2 u = = x i r + U r x2 i r r r 2, N u xix i = U rr + N 1 r i=1 1 N 1 u r N 1 r. r r 1 r N 1 r r N 1 u = 0. r This ODE is called a Cauchy-Euler equation. Written differently the equation is r 2 U rr + (N 1)rU r = 0. Its form suggests U = r α might work: we find (α(α 1) + (N 1)α)r α = 0. For N > 2, we get α = 0 and α = 2 N. Thus the general solution is c 1 + c 2 r 2 N. If N = 2, the Cauchy-Euler equation becomes 1 r (ru ) = 0. U r
5 8.3. Poisson s equation with Dirichlet and Neumann Boundaries 137 We can directly integrate this to find U = c 1 + c 2 ln r. We therefore have that a solution to 2 Φ(x) = 0 is c2 ln r N = 2 Φ(x) = c 2 r N 3, N 2 where we have dropped the trivial solution c 1. To find the constants we, formally, use Green s second identity, (8.1), and, for convenience, we set x = 0 and integrate over a ball or radius r. We find 1 = = = B r(0) B r(0) B r(0) = c 2 2π. 2 x Φ(x, x ) dx Φ(x, x ) 2 x 1 dx + x Φ(x, x ) n = B r(0) 2π 0 (1 x Φ(x, x ) Φ(x, x ) x 1) n dl c 2 r rdθ Thus c 2 = 1/2π. A similar calculation shows c 2 = 1/4π in three dimensions. In Cartesian coordinates we define the fundamental solution to be ln( x 2 +y 2 ) 2π N = 2 Φ(x) =, 1 N = 3 4π(x 2 +y 2 +z 2 ) 1/2 (check Φ solves the PDE!!). When we write Φ(x, x ) we mean Φ(x x ) as Φ is only a function of the distance between its arguments. In the next section we deal with boundaries Poisson s equation with Dirichlet and Neumann Boundaries Next we consider Poisson s equation with either Dirichlet boundary conditions. Let be an open bounded set. Poisson s equation is 2 u = f, u = g(x). An examination of Green s second identity, (8.1), we need to require G(x, x ) = 0 for x since u n is not known on the boundary. We expect the solution to be u(x) = G(x, x )f(x ) d + u(x ) x G(x, x ) n dl The only Green s function so far is the fundamental solution Φ. Let us set G(x, x ) = Φ(x, x ) + R(x, x ) and require R(x, x ) to satisfy 2 R = 0 for all x, x. Moreover, we require R(x, x ) = Φ(x, x ) for all x. To summarize
6 Elliptic PDEs Definition 8.1. The Green s function, G(x, x ), for the Dirichlet problem satisfies a) 2 x G(x, x ) = 0 for x = x. b) G(x, x ) = 0 for all x. c) The Green s function has the decomposition G(x, x ) = Φ(x, x ) + R(x, x ) where 2 x R(x, x ) = 0 for all x, x and R(x, x ) = Φ(x, x ) for all x, and Φ(x, x) is the fundamental solution. Theorem 8.2. Suppose u C 2 () solves 2 u = f, u = g(x), and G is a Green s function for the Dirichlet problem. Then u may be written u(x) = G(x, x )f(x ) dx + g(x ) x G(x, x ) n dl. Proof. We can make a rigorous proof by following Theorem 8.6 below. Next we consider 2 u = f, u n = g(x). In general this PDE will not have a solution (see homework problem (8.2)), and if a solution does exist, it is not unique since a constant may be added to it. We ignore these issues for now. Definition 8.3. The Green s function, G(x, x ), for the Neumann problem satisfies a) 2 x G(x, x ) = 0 for x = x. b) x G(x, x ) n = 0 for all x. c) The Green s function has the decomposition G(x, x ) = Φ(x, x ) + R(x, x ) where 2 x R(x, x ) = 0 for all x, x and x R(x, x ) n = x Φ(x, x ) n for all x. Theorem 8.4. Suppose u C 2 () solves 2 u = f, u n = g(x), and G is a Green s function for the Neumann problem. Then u may be written u(x) = G(x, x )f(x ) dx G(x, x )g(x ) dl. Proof. This identical to the proof above for Dirichlet problem.
7 8.3. Poisson s equation with Dirichlet and Neumann Boundaries 139 Theorem 8.5. The Green s function is symmetric. That is, G(x, x ) = G(x, x) for all x, x. Proof. This is a repeat of Theorem 6.7 in Chapter 4. Here, L = L so the proof is easier. Set u = G(x 1, x ) and v = G(x 2, x ) in Green s second identity. We find, formally, [G(x 1, x ) 2 x G(x 2, x ) G(x 2, x ) 2 x G(x 1, x )] d (8.2) = G(x 1, x ) x G(x 2, x ) n dl G(x 2, x ) x G(x 1, x ) n dl. The area integral on the left is equal to G(x 2, x 1 ) G(x 1, x 2 ) while the integrals on the right are both zero with either Dirichlet or Neumann boundary conditions. Let us prove at least one theorem rigorously. Theorem 8.6. (Representation Theorem) Suppose R 2 is a normal domain and that u C 2 (). Then for all x, u(x) can be writen as a double, single, and volume potential. That is, u(x) = u(y) y Φ(x, y) n dl y Φ(x, y) u(y) n dl y + Φ(x, y) 2 u(y) dy. Proof. The key is again Green s second identity. If u and v are smooth functions, then (u 2 v v 2 u) d = u v n dl v u n dl. We wish to apply this identity to u and Φ. However, Φ has a singularity when x = y. To avoid the singularity we remove a disk from around x. Let x be fixed and let ǫ < dist(x, ). Let B ǫ (x) be the disk of radius ǫ centered at x, and ǫ = /B ǫ (x) (all points in not in B ǫ (x)). Applying Green s second identity to u, v = Φ on ǫ. Since 2 Φ = 0 on ǫ, we get Φ(x, y) 2 u(y) dy = (8.3) ǫ u(y) y Φ(x, y) n dl y Φ(x, y) u(y) n dl y. ǫ ǫ We need to take the limit as ǫ 0. We will evaluate the limits in the order they appear in (8.3). Note that Φ(x, y) 2 u(y) dy ǫ B ǫ(x) Φ(x, y) 2 u(y) dy.
8 Elliptic PDEs Since u C 2 (), 2 u has a maximum on. That is, 2 u(y) M for some M and all y. Thus Φ(x, y) 2 ln( x y ) u(y) dy 2 u(y) dy 2π B ǫ(x) B ǫ(x) ln( x y ) M dy 2π Thus B ǫ(x) = M = M B ǫ(0) 2π ǫ lim Φ(x, y) 2 u(y) dy = ǫ 0 ǫ 0 ln(r) r dr dθ 2π 0 ln(r) 2π r dr dθ Mǫ2. Φ(x, y) 2 u(y) dy. The next two boundary integrals in (8.3) are over two surfaces: and B ǫ (x). We are only interested in the latter surface. The first boundary integral in (8.3) is u(y) 1 2π y ln x y n dl y. Here B ǫ(x) n = x y r x1 y 1 = r, x 2 y 2, r with r 2 = (x 1 y 1 ) 2 + (x 2 y 2 ) 2. Moreover, x1 y 1 y ln x y = r 2, x 2 y 2 r 2, and y ln x y n = r2 r 3. We have u(y) 1 B ǫ(x) 2π y ln x y n dl y = 1 2π 2π 0 u(x + ǫˆθ) 1 ǫ ǫdθ, where ˆθ = (cos θ, sin θ) (the integration is clockwise here). Since u is continuous, lim u(y) 1 ǫ 0 ǫ 2π y ln x y n dl y = u(y) 1 2π y ln x y n dl y = lim u(x + ǫˆθ) 1 2π 0 ǫ ǫdθ u(y) 1 2π y ln x y n dl y u(x). ǫ 0 1 2π
9 8.4. Examples 141 We leave it as an exercise to show as ǫ 0. B ǫ(x) Φ(x, y) u(y) n 0 Finally, taking the limit as ǫ 0 on both sides of (8.3) we get the representation for u stated in the theorem. Remark 8.7. If we believe 2 Φ(x, y) = δ(x y) and blindly apply Green s second identity, we get the representation for u stated in Theorem 8.6. Thus we can use delta functions as a quick way of correctly calculating improper integrals. Corollary 8.8. The Green s function in Definitions 8.1 and 8.3 satisfies u(x) = u(y) y G(x, y) n dl y G(x, y) u(y) n dl y + G(x, y) 2 u(y) dy. Proof. Since G(x, y) = Φ(x, y) + R(x, y) and R is harmonic on, R satisfies 0 = u(y) y R(x, y) n dl y R(x, y) u(y) n dl y + R(x, y) 2 u(y) dy. It therefore can be added to Φ without changing the statement in the theorem Examples Example 8.9. Let = {(x, y) x R, y > 0}, and consider 2 u = 0, u = g(x). We implement Definition 8.1. We look for a Green s function of the form G(x, x ) = 1 2π ln x x + R(x, x ) with R = Φ when x. Set x = (x, y) and x = (x, y ). We need R to be harmonic ( 2 R = 0) and R = Φ on the boundary. The only non trivial harmonic function we have is ln r. Let r 1 = x x. We will put a phoney point outside of the domain which we call x. Set r 2 = x x. Then the Green s function has the form G(x, x ) = 1 2π ln r 1 1 2π ln r 2. We need to chose x so that r 1 = r 2 when x is on the boundary (y = 0 here) and r 2 = 0 for all x. To accomplish this we set x = (x, y) (sketch a picture of the three points and imagine x on the boundary). Thus G(x, y, x, y ) = 1 4π ln (x x ) 2 + (y y ) 2 1 4π ln (x x ) 2 + (y + y ) 2 = 1 4π ln (x x ) 2 + (y y ) 2 (x x ) 2 + (y + y ) 2 It is easy to check this G satisfies the properties listed in Definition 8.1.
10 Elliptic PDEs By Theorem 8.3 or Equation (8.2) (here n = (0, 1)) u(x, y) = g(x ) x G(x, x ) n dl = = y π G y y =0 g(x )dx g(x ) (x x ) 2 + y 2 dx. To be consistent we should check our solution directly. For simplicity we take the space L () to mean the set of bounded continuous functions on. For a function g L we denote the bound on g to be g. Theorem Suppose g C(R) L (R). Then (i) u C () L (), (ii) 2 u = 0 in, (iii) lim u(x, y) = g(x 0). (x,y) (x 0,0) (iv) u C(). Proof. The integral describing u is challenging to work with since it is not only an improper integral because of the domain, but also the integrand is unbounded when y = 0. Set K(x, y, x y ) = π[(x x ) 2 + y 2 ]. Our solution is u(x, y) = K(x, y, x )g(x ) dx. For y > 0 we make the substitutions z = x x and z = y tan θ. We compute It follows K(x, y, x ) dx = u(x, y) = = = = 1. π/2 π/2 π/2 π/2 y π[(x x ) 2 + y 2 ] dx y π[z 2 + y 2 ] dz y 2 sec 2 θ πy 2 (1 + tan 2 θ) dθ dθ π K(x, y, x )g(x ) dx g K(x, y, x ) dx = g.
11 8.4. Examples 143 This shows u L (). Next we show u C(). Let x, y be arbitrary with y > 0. Set for some L, y 0 and L y. Set L = {(x, y) L x L, 0 < y 0 y L y } u n (x, y) = n n g(x )y π((x x ) 2 + y 2 ) dx. Then u n C( L ). We need to show u n converges uniformly to u on L. Then Theorem C, Theorem (1.41), applies and u C( L ). Since x, y is arbitrary, it follows u C(). Clearly u n u pointwise as n. We have, after setting Δu = u n (x, y) u(x, y) and z = x x, n y Δu g π((x x ) 2 + y 2 ) dx + = g n x = g n+l 2g y π 2g y π 2g L y π(n L). n L n L y π(z 2 + y 2 ) dz + y π(z 2 + y 2 ) dz + dz z 2 + y 2 dz z 2 n x n L n y π((x x ) 2 + y 2 ) dx y π(z 2 + y 2 ) dz y π(z 2 + y 2 ) dz This shows u n (x, y) converges uniformly to u(x, y) on L. It follows u C(). The proof u C () is similar only we repeatedly apply Theorem D, Theorem We outline one iteration. Let x, y be given and let L be defined as before. We need to check the conditions of Theorem D. Since x K(x, y, x ) and K(x, y, x ) are continuous on L, Theorem 1.8 applies and x u n (x, y) = n n x K(x, y, x )g(x ) dx is continuous on L. We have u n (x, y) converges at one point (all points) in L. To show x u n converges uniformly on L we show the sequence { x u n } is uniformly Cauchy. Note that x K(x, y, x ) = 2y π x x [(x x ) 2 + y 2 ] 2.
12 Elliptic PDEs Following the same arguments as above we have, assuming n > m x u n x u m m n 2y π 4y π g x x [(x x ) 2 + y 2 ] 2 dx + n+l m L 4L y g 2π(m L) 2 z (z 2 + y 2 ) 2 dz n m 2y π x x [(x x ) 2 + y 2 ] 2 dx for n > m. This shows the convergence in uniform and Theorem D applies. Thus lim n x u n = x lim n u n. That is, u x (x, y) = The other derivatives are found in a similar way. x K(x, y, x )g(x ) dx. To establish (iii) in the theorem fix x 0 R. Since g C(R), for all ǫ > 0 there exists a δ > 0 such that g(x) g(x 0 ) < ǫ/2 for all x such that x x 0 < δ. Thus, since K(x, y, x ) dx = 1, u(x, y) u(x 0, y) y g(x ) g(x 0 ) π (x x ) 2 + y 2 dx For the first integral we have x0+δ + x 0 δ = I + II. y g(x ) g(x 0 ) π (x x ) 2 + y 2 dx R\B δ (x 0) y g(x ) g(x 0 ) π (x x ) 2 + y 2 dx I < x0+δ x 0 δ x0+δ x 0 δ ǫ/2. K(x, y, x g(x ) g(x 0 ) dx K(x, y, x )ǫ dx To estimate the second integral we need a lower bound on x x. Let us require x to satisfy x x 0 < δ/2. The integrand in II requires requires x < x 0 δ or x > x 0 + δ, and so x x 0 δ. Thus That is, x x 0 x x + x x 0 x x + δ/2 x x x x 0. x x 1 2 x x 0.
13 8.4. Examples 145 This gives as a bound on the second integral. Indeed y g(x ) g(x 0 ) II = R\B δ (x 0) π (x x ) 2 + y 2 dx 2g y R\B δ (x 0) π (x x 0) dx 4g δπ y. Putting the two together we find u(x, y) g(x 0 ) < ǫ/2 + 4g δπ y for all x x 0 < δ/2. Thus by restricting y, for every ǫ > 0 there is a neighborhood of (x 0, 0) such that u(x, y) g(x 0 ) < ǫ, and Part (iii) in the theorem follows. Standard continuity arguments now show u C(). Example Let = {(x, y) x 2 + y 2 = a 2, a > 0}, and consider 2 u = 0, u = g(x). Again we implement Definition 8.1. We look for a Green s function of the form G(x, x ) = 1 2π ln x x + R(x, x ). Following the insight we gained in the previous example, we place a phoney point, x, somewhere in the plane so that G(x, x ) = 1 2π ln x x 1 2π ln x x = 0 when x. Moreover, the symmetry of the problem suggest putting x on the line formed by x and the origin. and Let x = r, x = ρ, x = R and φ the the angle between x and x. Then x x 2 = x x + x x 2x x = r 2 + ρ 2 2rρ cos φ, x x 2 = x x + x x 2x x = ρ 2 + R 2 2ρR cos φ. If we set x x = x x when x (when ρ = a), we would conclude R = r. This is impossible since we need x outside of the domain. Instead of looking for the ratio x x / x x = 1, we try x x / x x = γ for x. We expect γ = γ(a, r). If such a γ can be found, our Green s function is G(x, x ) = 1 2π ln x x 1 2π ln(γ x x ).
14 Elliptic PDEs Setting x x 2 = γ 2 x x 2 when ρ = a, and equating the cos φ terms, we find r = Rγ 2. Equating the remaining terms shows γ = a/r. This in turn implies γ = r a. Thus the Green s function is G(x, x ) = 1 4π ln x x 2 1 = 1 4π ln a 2 r 2 4π ln(γ2 x x 2 ) = 1 r 2 + ρ 2 2rρ cos φ ρ 2 + R 2 2ρR cos φ x x 4π ln 2 γ 2 x x 2 If x = (r cos θ, r sin θ) and x = (ρ cos θ, ρ sin θ ), then φ = θ θ. Moreover, since R = a 2 /r, G(x, x ) = 1 a 2 4π ln (r 2 + ρ 2 2rρ cos(θ θ )) r 2 ρ 2 + a 4 2ra 2 ρ cos(θ θ ) By Theorem 8.3 or Equation 8.2 the solution may be written as u(x, y) = g(x ) x G(x, x ) n dl. Here and dl = a dθ. Thus, x G(x, x ) n = G ρ (8.4) u(r, θ) = 1 2π π π = a ρ=a 2π Equation (8.4) is called Poisson s formula. 1 (r/a) 2 r 2 + a 2 2ar cos(θ θ ), a 2 r 2 r 2 + a 2 2ar cos(θ θ ) g(θ ) dθ. Theorem Equation (8.4) satisfies 2 u = 0 in and u C (). Proof. We recall Theorem 1.34 from Chapter 1. Namely, if with F and F x continuous, then f(x) = f (x) = b a b a F (x, y) dy, F (x, y) x Now we observe that, if r < a, the dominator in (8.4) is never zero and we may move as many derivatives inside the integral as we wish. That is, u C (). We leave it as an exercise to check 2 u = Maximum Principle Based on the physical problems Laplace s equation describes, we expect the maximum and minimum of a harmonic function ( 2 u = 0) to occur on the boundary. This is easy to show. dy.
15 8.5. Maximum Principle 147 Theorem Let be a open bounded set in R n. Suppose u C 2 () C() satisfies 2 u 0. Then max u = max u. Proof. Suppose for the moment 2 u < 0 in. Since u is continuous on a closed bounded set, u has a maximum. Suppose the maximum is in. By the second derivative test 2 u 0. This implies 2 u 0 contrary to our assumption. Thus the maximum occurs on the boundary. Now suppose 2 u 0 in. Set, for ǫ > 0, v = u + ǫ x2 2n. Then 2 v = 2 u ǫ < 0. By the result in the previous paragraph max where d = max u max Since, we have and the result follows. v = max v = max u + ǫ x2 2n x 2. Since this holds for all ǫ > 0, 2n max u max u. max u max u, max u + ǫd, We have not shown a solution to an elliptic PDE depends continuously on the data yet. Theorem Consider on a bounded open set. Set N = sup Proof. Set v = u + Nx2 2n. Then, in, By the previous theorem u max u max 2 u = f, u = g(x) f and M = sup g. Then, for all x, u(x) M + Nd. 2 v = 2 u N 0. v = max v = max u + Nx2 M + Nd. 2n Similarly, set v = u + Nx2 2n. Then, in, 2 v = 2 u N 0. Thus u v max v = max u v = max + Nx2 M + Nd. 2n That is, u (M + Nd). The result follows.
16 Elliptic PDEs Corollary The solution to depends continuously on f and g. 2 u = f, u = g(x) Proof. Let u i, for i = 1, 2 solve Then the difference w = u i u 2 solves 2 u i = f i, u i = g i (x). 2 w = Δf, w = Δg, where Δf = f 1 f 2 and Δg = g 1 g 2. Thus w(x) = u 1 (x) u 2 (x) sup Δg + sup Δf d, and we see small changes in f and g lead to small changes in u.
17 8.6. Homework Homework Exercise 8.1. Suppose is a normal domain. Suppose a function u is specified on the boundary,, of and we want to find the minimal surface over subject our constraint u at the boundary. That is, we want to minimize the functional (minimal surface area) J(u) := 1 + u 2 x + u2 y dxdy. Minimize J by considering an arbitrary smooth function ϕ which is zero on the boundary of. Then consider J(u + ǫϕ) := 1 + (u + ǫϕ) 2 x + (u + ǫϕ)2 y dxdy. Now compute and require J(u+ǫϕ) ǫ ǫ=0 = 0 to derive Plateau s equation (1 + u 2 y )u xx + (1 + u 2 x )u yy 2u x u y u xy = 0. Exercise 8.2. Let be a bounded, normal domain. Consider the elliptic PDE 2 u = f, x u n = g. (a) Prove no solution exists if f d = g dσ. What is the physical significance of this condition? (b) If the condition is satisfied, is the solution unique? Provide a proof either way. (c) Find a condition which will make the solution unique. Exercise 8.3. Consider the elliptic PDE on the square = {(x, y) 0 < x < π, 0 < y < π}: 2 u = 0, on (a) Find the formal solution. u x (0, y) = f(y), 0 < y < π, u x (π, y) = g(y), 0 < y < π, u(x, 0) = 0, 0 < x < π, u(x, π) = 0, 0 < x < π. (b) Find conditions on f and g so that your formal solution satisfies u C, () C().
18 Elliptic PDEs Exercise 8.4. Consider the Laplace s equation on the channel = {(x, y) 0 < x <, 0 < y < π}: 2 u = 0, on (a) Find the formal solution. u(0, y) = f(y), 0 < y < π, u(x, 0) = 0, 0 < x <, u(x, π) = 0, 0 < x <. (b) Find conditions on f so that your formal solution satisfies u C, () C(). Exercise 8.5. Let = {(x, y) R 2 : 0 < y, 0 < x}. Solve, by constructing a Green s function, 2 u = 0 0 < x, 0 < y, u(x, 0) = f(x), 0 < x, u(0, y) = 0, 0 < y, u(x, y) is bounded for large x, y and f(x) is continuous and bounded. The answer is u(x, y) = y 1 π (x x ) 2 + y 2 1 (x + x ) 2 + y 2 f(x )dx. 0 Exercise 8.6. Let = {(x, y) R 2 : 0 < y, 0 < x}. Solve, by constructing a Green s function, 2 u = 0, 0 < x, 0 < y, u n = f(x), 0 < x, u(0, y) = 0, 0 < y. (x,0) Exercise 8.7. Let be the semicircle (0 < r < a, 0 < θ < π). Find the Green s function, G(x, x ) satisfying 2 G = δ( x x ), with G(x, x ) = 0 when x is on the boundary. Exercise 8.8. Consider the PDE with Robin boundary conditions on an open set 2 u = f, u n + α(x)u = g(x) x, where α(x) > 0 on. What properties does a Green s function have for this PDE. Formulate a definition like Definition 1.4 or 1.5 in the notes. Find, formally, a representation, using the Green s function, of the solution - as in Theorem 1.7 or 1.8 in the notes. Exercise 8.9. Let be an open, normal, bounded subset of R 2. Prove that the Green s function to Poisson s equation with Dirichlet boundary conditions is unique. Hint: suppose there were two Green s functions. Set u(y) = G 1 (x, y) G 2 (x, y) for a fixed x. Show 2 u = 0 and u = 0 on. 1. Show Φ(x, y) u(y) n 0 B ǫ(x) as ǫ 0 to complete the proof of Theorem 8.3
19 8.6. Homework Make the proof the Green s function is symmetric rigorous by applying Green s second identity on ǫ and take the ǫ 0 as was done in the proof of Theorem 8.3 Exercise If is not bounded, show that the solution to 2 u = 0 may not have a maximum. Exercise Consider the Helmholtz equation on 2 u u = h u = f. (a) Suppose is normal and bounded. Use energy methods to show the solution is unique. STATE the assumptions you need on u. (b) Suppose the boundary of is fractal (Not smooth so energy methods do not apply). By using ideas from elementary calculus, develop a maximum principle to show the solution is unique. (You will have to home cook your own maximum principle.) Exercise Consider Laplace s equation on the square 0 x 1, 0 y 1 u with nonhomogeneous Neumann boundary conditions: ν = 0 on the sides {y = 0, 0 x 1}, {x = 0, 0 y 1}, {x = 1, 0 y 1}, and u ν = f(x) on the side {y = 1, 0 x 1}. (a) Show that no solution exists if the condition (b) Is the solution unique? Prove or disprove. 1 0 f(x) dx = 0 is not satisfied. Exercise Let be an open, normal, bounded subset of R n. Let u C 2 () C(). Moreover, suppose that u(x) satisfies 2 u(x) = 0, x. (a) Prove max u = max u. x x (b) Prove the solution to Poisson s equation is unique using two different methods. Exercise Suppose that u, v are C 2 () C() and solve 2 u = f 1 2 v = f 2 u = g 1 v = g 2. Moreover, let f 1 f 2 0 for all x, and g 1 g 2 0 for all x. Prove that u v 0 for all x.
MATH 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 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 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 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 informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
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 informationUsing Green s functions with inhomogeneous BCs
Using Green s functions with inhomogeneous BCs Using Green s functions with inhomogeneous BCs Surprise: Although Green s functions satisfy homogeneous boundary conditions, they can be used for problems
More informationClass Meeting # 2: The Diffusion (aka Heat) Equation
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 2: The Diffusion (aka Heat) Equation The heat equation for a function u(, x (.0.). Introduction
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 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 information4.10 Dirichlet problem in the circle and the Poisson kernel
220 CHAPTER 4. FOURIER SERIES AND PDES 4.10 Dirichlet problem in the circle and the Poisson kernel Note: 2 lectures,, 9.7 in [EP], 10.8 in [BD] 4.10.1 Laplace in polar coordinates Perhaps a more natural
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
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 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 informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH College of Informatics and Electronics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS425 SEMESTER: Autumn 25/6 MODULE TITLE: Applied Analysis DURATION OF EXAMINATION:
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 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 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 informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationLAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are
LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or
More informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationGeneral Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436
General Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436 1 Introduction You may have noticed that when we analyzed the heat equation on a bar of length 1 and I talked about the equation
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 informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationare harmonic functions so by superposition
J. Rauch Applied Complex Analysis The Dirichlet Problem Abstract. We solve, by simple formula, the Dirichlet Problem in a half space with step function boundary data. Uniqueness is proved by complex variable
More information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More informationExistence and Uniqueness
Chapter 3 Existence and Uniqueness An intellect which at a certain moment would know all forces that set nature in motion, and all positions of all items of which nature is composed, if this intellect
More informationElectrodynamics PHY712. Lecture 3 Electrostatic potentials and fields. Reference: Chap. 1 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 3 Electrostatic potentials and fields Reference: Chap. 1 in J. D. Jackson s textbook. 1. Poisson and Laplace Equations 2. Green s Theorem 3. One-dimensional examples 1 Poisson
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
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 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 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 informationChapter 3. Second Order Linear PDEs
Chapter 3. Second Order Linear PDEs 3.1 Introduction The general class of second order linear PDEs are of the form: ax, y)u xx + bx, y)u xy + cx, y)u yy + dx, y)u x + ex, y)u y + f x, y)u = gx, y). 3.1)
More informationSECTION 3.6 THE CONSTRUCTION OF THE GREEN'S FUNCTION FOR THE TWO DIMENSIONAL OPERATOR WITH DIRICHLET BOUNDARY CONDITIONS
SECTION 3.6 THE CONSTRUCTION OF THE GREEN'S FUNCTION FOR THE TWO DIMENSIONAL OPERATOR WITH DIRICHLET BOUNDARY CONDITIONS We now will construct a Green's function for the following problem: Find u such
More informationSOLUTION OF POISSON S EQUATION. Contents
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 1 2 3. The Newtonian potential from the 3rd Green
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
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 informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationContinuity. Chapter 4
Chapter 4 Continuity Throughout this chapter D is a nonempty subset of the real numbers. We recall the definition of a function. Definition 4.1. A function from D into R, denoted f : D R, is a subset of
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
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 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 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 informationMATH20411 PDEs and Vector Calculus B
MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial ifferential Equations «Viktor Grigoryan 3 Green s first identity Having studied Laplace s equation in regions with simple geometry, we now start developing some tools, which will lead
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 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 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 informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationLaplace equation. In this chapter we consider Laplace equation in d-dimensions given by. + u x2 x u xd x d. u x1 x 1
Chapter 6 Laplace equation In this chapter we consider Laplace equation in d-dimensions given by u x1 x 1 + u x2 x 2 + + u xd x d =. (6.1) We study Laplace equation in d = 2 throughout this chapter (excepting
More information= 2 x x 2 n It was introduced by P.-S.Laplace in a 1784 paper in connection with gravitational
Lecture 1, 9/8/2010 The Laplace operator in R n is = 2 x 2 + + 2 1 x 2 n It was introduced by P.-S.Laplace in a 1784 paper in connection with gravitational potentials. Recall that if we have masses M 1,...,
More informationContinuity. Chapter 4
Chapter 4 Continuity Throughout this chapter D is a nonempty subset of the real numbers. We recall the definition of a function. Definition 4.1. A function from D into R, denoted f : D R, is a subset of
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
More informationTHE METHOD OF SEPARATION OF VARIABLES
THE METHOD OF SEPARATION OF VARIABES To solve the BVPs that we have encountered so far, we will use separation of variables on the homogeneous part of the BVP. This separation of variables leads to problems
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationReview and problem list for Applied Math I
Review and problem list for Applied Math I (This is a first version of a serious review sheet; it may contain errors and it certainly omits a number of topic which were covered in the course. Let me know
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 informationHyperbolic PDEs. Chapter 6
Chapter 6 Hyperbolic PDEs In this chapter we will prove existence, uniqueness, and continuous dependence of solutions to hyperbolic PDEs in a variety of domains. To get a feel for what we might expect,
More informationMidterm 2: Sample solutions Math 118A, Fall 2013
Midterm 2: Sample solutions Math 118A, Fall 213 1. Find all separated solutions u(r,t = F(rG(t of the radially symmetric heat equation u t = k ( r u. r r r Solve for G(t explicitly. Write down an ODE for
More information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
More informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
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 informationSAMPLE FINAL EXAM SOLUTIONS
LAST (family) NAME: FIRST (given) NAME: ID # : MATHEMATICS 3FF3 McMaster University Final Examination Day Class Duration of Examination: 3 hours Dr. J.-P. Gabardo THIS EXAMINATION PAPER INCLUDES 22 PAGES
More information1. Differential Equations (ODE and PDE)
1. Differential Equations (ODE and PDE) 1.1. Ordinary Differential Equations (ODE) So far we have dealt with Ordinary Differential Equations (ODE): involve derivatives with respect to only one variable
More informationPartial Differential Equations
Partial Differential Equations Spring Exam 3 Review Solutions Exercise. We utilize the general solution to the Dirichlet problem in rectangle given in the textbook on page 68. In the notation used there
More informationElectrodynamics PHY712. Lecture 4 Electrostatic potentials and fields. Reference: Chap. 1 & 2 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 4 Electrostatic potentials and fields Reference: Chap. 1 & 2 in J. D. Jackson s textbook. 1. Complete proof of Green s Theorem 2. Proof of mean value theorem for electrostatic
More informationAPPM 2350 Final Exam points Monday December 17, 7:30am 10am, 2018
APPM 2 Final Exam 28 points Monday December 7, 7:am am, 28 ON THE FONT OF YOU BLUEBOOK write: () your name, (2) your student ID number, () lecture section/time (4) your instructor s name, and () a grading
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationPlot of temperature u versus x and t for the heat conduction problem of. ln(80/π) = 820 sec. τ = 2500 π 2. + xu t. = 0 3. u xx. + u xt 4.
10.5 Separation of Variables; Heat Conduction in a Rod 579 u 20 15 10 5 10 50 20 100 30 150 40 200 50 300 x t FIGURE 10.5.5 Example 1. Plot of temperature u versus x and t for the heat conduction problem
More informationRiemann integral and Jordan measure are generalized to unbounded functions. is a Jordan measurable set, and its volume is a Riemann integral, R n
Tel Aviv University, 214/15 Analysis-III,IV 161 1 Improper integral 1a What is the problem................ 161 1b Positive integrands................. 162 1c Newton potential.................. 166 1d Special
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 informationMathematics of Physics and Engineering II: Homework problems
Mathematics of Physics and Engineering II: Homework problems Homework. Problem. Consider four points in R 3 : P (,, ), Q(,, 2), R(,, ), S( + a,, 2a), where a is a real number. () Compute the coordinates
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 informationG: Uniform Convergence of Fourier Series
G: Uniform Convergence of Fourier Series From previous work on the prototypical problem (and other problems) u t = Du xx 0 < x < l, t > 0 u(0, t) = 0 = u(l, t) t > 0 u(x, 0) = f(x) 0 < x < l () we developed
More informationc2 2 x2. (1) t = c2 2 u, (2) 2 = 2 x x 2, (3)
ecture 13 The wave equation - final comments Sections 4.2-4.6 of text by Haberman u(x,t), In the previous lecture, we studied the so-called wave equation in one-dimension, i.e., for a function It was derived
More informationModule 7: The Laplace Equation
Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
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 informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
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 informationPrinciple of Mathematical Induction
Advanced Calculus I. Math 451, Fall 2016, Prof. Vershynin Principle of Mathematical Induction 1. Prove that 1 + 2 + + n = 1 n(n + 1) for all n N. 2 2. Prove that 1 2 + 2 2 + + n 2 = 1 n(n + 1)(2n + 1)
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationn 2 xi = x i. x i 2. r r ; i r 2 + V ( r) V ( r) = 0 r > 0. ( 1 1 ) a r n 1 ( 1 2) V( r) = b ln r + c n = 2 b r n 2 + c n 3 ( 1 3)
Sep. 7 The L aplace/ P oisson Equations: Explicit Formulas In this lecture we study the properties of the Laplace equation and the Poisson equation with Dirichlet boundary conditions through explicit representations
More informationSuggested Solution to Assignment 6
MATH 4 (6-7) partial diferential equations Suggested Solution to Assignment 6 Exercise 6.. Note that in the spherical coordinates (r, θ, φ), Thus, Let u = v/r, we get 3 = r + r r + r sin θ θ sin θ θ +
More informationMath 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
More informationHankel Transforms - Lecture 10
1 Introduction Hankel Transforms - Lecture 1 The Fourier transform was used in Cartesian coordinates. If we have problems with cylindrical geometry we will need to use cylindtical coordinates. Thus suppose
More informationBessel functions. The drum problem: Consider the wave equation (with c = 1) in a disc with homogeneous Dirichlet boundary conditions:
Bessel functions The drum problem: Consider the wave equation (with c = 1) in a disc with homogeneous Dirichlet boundary conditions: u = u t, u(r 0,θ,t) = 0, u u(r,θ,0) = f(r,θ), (r,θ,0) = g(r,θ). t (Note
More informationIntroduction and some preliminaries
1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.
More information