Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
|
|
- Bernice Fowler
- 6 years ago
- Views:
Transcription
1 Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and spherical coordinate systems. In this lecture separation in cylindrical coordinates is studied, although Laplaces s equation is also separable in up to 22 other coordinate systems as previously tabulated. Thus one chooses the system in which the appropriate boundary conditions can be applied. The application of the boundary conditions (Dirichlet of Neumann on a closed surface) results in the unique solution to a specific problem. In cylindrical coordinates the divergence of the gradient of the potential yields; 2 V (ρ, φ, z) = ρ 2 V ρ 2 + V ρ + (1/ρ) 2 V φ V z 2 = Look for solutions by separation of variables; V = R(ρ)Ψ(φ)Z(z) As previously, one finds 2 separation constants, k and ν, which result in the three ode equations with 2 of these providing 2 eigenfunctions. d 2 Z dz 2 k 2 Z = d 2 Ψ dφ 2 + ν 2 Ψ = d 2 R dρ 2 + (1/ρ) dr dρ + (k2 (ν/ρ) 2 )R = Solution of the later 2 equations result in eigenfunctions. Solutions are; Z e ±kz Ψ e ±iνφ R J ν (kρ), N ν (kρ) 1
2 Figure 1: An example of the Cylindrical Bessel function J ν (x) as a function of x showing the oscillatory behavior 2 Bessel Functions In the last section, J ν (kρ), N ν (kρ) are the 2 linearly independent solutions to the separated ode radial equation. This is Bessel s equation with Bessel functions as solutions. Bessel functions oscillate, but are not periodic like harmonic functions, see Figure 1. Thus the harmonic function solution for Ψ and the Bessel function solution for R result in eigenfunctions when the boundary conditions are imposed. The Bessel functions, J ν (x), are regular at x =, while the Bessel (Neumann) functions, N ν (x), are singular at x =. From the requirement that the solution is single valued as φ 2π, ie the solution must not change when φ is replaced by φ + 2π, the values of ν are integral and this produces eigenfuntions in Ψ. The series solution for the Bessel function J ν can be found by the method of Frobenius. However, the second linearly independent equation is not easily obtained when n is an integer. The series representation of the Bessel function takes the form; J ν (x) = s= ( 1) s s!(ν + s)! (x/2)ν+2s = xν 2 n n! x ν+2 2 ν+2 (ν + 1)! + For integral values of ν one can show J ν recurrence relations; J ν 1 (x) + J ν+1 (x) = ν x J ν(x) J ν 1 (x) = x ν J ν(x) + dj ν(x) dx d dx [xν J ν (x)] = x ν J ν 1 (x) = ( 1) ν J ν. Bessel functions also satisfy the 2
3 At times an integral representation is useful. J ν (x) = (1/π) π dθ e ix cos(θ) cos(νθ) Bessel functions are orthogonal, but they are not normalized. a kdk J ν (kρ) J ν (kρ ) = δ(ρ ρ )/ρ ρdρ J ν (α νk ρ/a) J ν (α νk ρ/a) = (a 2 /2)[J ν+1 (α νk )] 2 δ kk In the above α νk are the zeros of the Bessel function of order ν where k numbers these zeros. As the Bessel functions form a complete set, any function can be expanded in a Bessel series for a finite interval in ρ or a Bessel integral if the space is infinite. F(ρ) = kdk A(k) J ν (kρ) 3 Integral Representation Let; f n = (1/π) π dω cos(xsin(ω) nω) Then look at the relations for f n+1 and f n 1. Subtracting these we obtain; f n+1 f n 1 = (1/π) cos(xsin(ω) nω) cos(ω) + sin(xsin(ω) nω) sin(ω) f n+1 f n 1 = (2/π) π π dω [cos(xsin(ω) nω) cos(ω) + sin(xsin(ω) nω) sin(ω) dω sin[x(ω) nω)] sin(ω) The term on the right is twice f n. Thus the recursion relation for the Bessel function is reproduced, so identify f n J n. The integral representation is; J n (x) = (1/π) π dω cos(xsin(ω) nω) = (1/2π) π π i(x sin(ω) nω) dω e This implies that the Bessel function, J n, is the n th Fourier coefficient of the expansion; e i2 sin(ω) = J n e inω. The following expressions for the generating function are obtained. 3
4 cos(z sin(ω)) = sin(z sin(ω)) = Observe that J n = ( 1) n J n J n cos(nω) J n sin(nω) 4 Generating function Begin with the series in the last section; cos(z sin(ω)) = sin(z sin(ω)) = Let ω = π/2 to obtain; J n cos(nω) J n sin(nω) cos(z) = J 2J 2 + 2J 4 + sin(z) = 2J 1 2J 3 + Also let e iω = t so that e iω e iω = 2i sin(ω). Then; e iz sin(ω) = e z(t 1/t)/2 = This is the generating function for J n. G(z, t) = e z(t 1/t)/2 = J n t n J n (z) t n One then obtains the value of J n (x) by determining the coefficient of the n th power in the series. Thus; d n dt n [G(z, t)] t= = n! J n (z) The generating function can be used to establish the Bessel power series, and the recursion relations. 4
5 L z V = Vo V = x V = a y Figure 2: The geometry of a cylinder with one endcap held at potential V = V (ρ) and the other sides grounded 5 Limiting values Limiting Values for the Bessel functions are needed to apply boundary conditions. J n (x) lim x (x/2)ν Γ(ν + 1) N (x) lim x (2/π)ln(x) N ν (x) lim x (1/π)Γ(ν)(x/2) ν J ν (x) lim x 2/πxcos(x νπ/2 π/4) N ν (x) lim x 2/πxsin(x νπ/2 π/4) H 1 ν (x) = J ν + in ν lim x 2/πxe i(x νπ/2 π/4) H 2 ν(x) = J ν in ν lim x 2/πz e i(x νπ/2 π/4) 6 Example Find the solution for the interior of a cylindrical shell with the top end cap held at a potential V = V (ρ) and all the other surfaces grounded, Figure 2. The solution has the form; 5
6 Figure 3: A graphical representation of the above solution V = νn A νkn J ν (k n ρ)sinh(k n z)e iνφ In this case the solution is independent of the angle φ so ν =. Note that N ν is excluded in the solution because it must be finite as ρ =. Also choose sinh(k n z) to satisfy the boundary condition at z =. The reduced solution is; V = n A n J (k n ρ)sinh(k n z) Now V = for ρ = a. This means that; J (k n a) = The values of k n a are the zeros of the Bessel function J (k n a). The first few are, α n = 2.448, 5.521, ,. Then when z = L, A n is obtained using the orthogonality of the Bessel functions. The graphic form of the solution is shown in Figure 3. A n = 2 a 2 [J 1 (k n a)] 2 sinh(k n L) a ρdρ J (k n ρ) V (ρ) 7 Example As another example, find the potential inside a cylinder when the potential is specified on the end caps and the cylindrical wall is at zero potential, Figure 4. Note that the figure shows the coordinate origin at the center of the cylinder. Thus the solution is symmetric for ±z in order to match the boundary conditions. The boundary conditions are; 6
7 V = Vo sin( φ ) z a L y x V = L V = Vo sin( φ ) Figure 4: The geometry of a problem with endcaps held at potential V = ±V sin(φ) V = V sin(φ) z = L V = V sin(φ) z = -L V = ρ = a V = on the (x, y) plane when z = The solution must have the form; V = νn A νn J ν (k n ρ) sin(νφ) sinh(k n z) Here solutions cannot contain N νn (kρ) which are infinite at the origin. To match the boundary at z = ±L a a term sin(νφ) which requires ν = 1 is needed. Then require that the Bessel function, J 1n (k n a) =. This determines the zeros of the Bessel function of order 1. Write these as α 1n so that k n = α 1n /a. The solution has the form; V = n A n J 1 (α 1n ρ/a) sinh(αz/a) sin(φ) Finally match the boundary condition at z = ±L where V = V sin(φ). Use orthogonality to obtain; (L/2)[J 2 (α 1n )] 2 A n = 1 sinh(αl/a) a ρdρ J 1 (α 1n ρ/a) V 7
8 z b c a V = f( ρ ) V = y x Figure 5: The geometry for the problem of two concentric cylinders 8 Example As another example, look for a solution of concentric cylinders with the boundary conditions; r = a, c and z = V = z = b V = f(ρ) This geometry is shown in Figure 5. Choose the solution to have the form; V = A n sinh(k n z) G (k n ρ) n Here G is a superposition of the Bessel and Neumann functions. G = [ J (k n ρ) J (k n c) N (k n ρ) N (k n c) ] This choice results in G = when ρ = c, ie on the cylindrical surface of the inner cylinder. Choose ν = because the potential is independent of φ, ie the problem is azimuthally symmetric. Now choose the values of k n to make G = when ρ = a. This selects a set of zeros, α νn, of the function, G. Because G is a solution to the cylindrical Sturm-Liouville problem, these functions form a complete orthogonal set, as do all solutions to Bessel s equation which satisfy the boundary conditions. This points out that the solutions of the radial ode are placed in forms one of which is regular at ρ = and one which is not. However, other forms are possible such as G above. As all Bessel functions, G has oscillating properties, but of course the location of the zeros differs for the various representations. Finally, use orthogonality to obtain the coefficients in the above equation. 8
9 Here; a H A n = [1/sinh(k n b)] c H = a c ρdρ G 2 (α n ρ/a) ρdρ V (ρ) G(α n ρ/a) 9 Example Consider the problem with the cylindrical wall held at potential V = f(z) and the endcaps grounded. This geometry is shown in Figure 6. The boundary conditions are; z =, b V = ρ = a V = f(z) In this case the hyperbolic function in z cannot be used to match the boundary conditions. However let k ik so that this function becomes harmonic at the expense of making the argument of the Bessel function imaginary. Note here that the problem is 2-D since the solution is independent of angle. Thus there is only one eigenfunction, and this now occurs for the z coordinate. An eigenvalue of k = nπ/b is determined for this harmonic equation. The radial ode with a Bessel function solution of imaginary argument result in an eigenfunction. sin(nπz/b) with n integral The solution has the form; V = n=1 A n sin(nπz/b) J (inπρ/b) In this case use the orthogonality of the harmonic functions rather than the Bessel functions to find the expansion coefficients. This results in; A n = 2 bj (inπa/b) b dz f(z) sin(nπz/b) Bessel functions of imaginary argument are called Modified Bessel functions. They do not oscillate as a function of the argument and cannot be eigenfunctions. 9
10 b V = z a z V = f( ) y x V = Figure 6: The geometry for the problem where a potential is placed on the cylindrical surface and the end caps grounded 1 Dual Boundary Conditions The equation for the potential is 2 V =. In cylindrical coordinates, Figure 7; (1/ρ) ρ [ρ V ρ ] + (1/ρ2 ) 2 V θ V z 2 = Assume a solution of the form; V = [ ] sin(nθ) V n (z, ρ) n= cos(nθ) In the above, one uses the harmonic functions sine, cosine or a combination of both to satisfy the angular boundary conditions. Use sin(θ) as an example. 2π V n (z, ρ) = 1/π dθ V (z, ρ, θ) sin(nθ) This results in the pde for V n. (1/ρ) ρ [ρ V n ρ ] (n2 /ρ 2 ) V n + 2 V n z 2 = Apply a Hankel transform. That is, expand the function, V (ρ, θ, z) using the completeness integral given in the section on properties of Bessel functions. ρ dρ J n (λρ) [(1/ρ) ρ [ρ V n ρ ] (n2 /ρ 2 ) V n + 2 V n z 2 ] = V n (λ, z) = ρ dρ V n J n (λρ) 1
11 z x F a y Figure 7: The geometry of a flat, conducting disk of radius a held at a potential F The solution is symmetric for z z, so choose z >. Use eigenfunctions for J n (λρ) and cos(nθ). This means the functional form for z is e λz. Choose n = to simplify the exposition of the Hankel transform given above. V n= (ρ, z) = λ dλ A(λ) e λz J (λρ) The boundary conditions are applied in two steps. Step 1 - z = and ρ < a V (ρ, ) = F(ρ) = Step 2 - z = and ρ > a V ρ z= = λ dλ A(λ) J (λρ) λ 2 dλ A(λ) J (λρ) = E ρ = In the second equation above, E ρ is the field in the ˆρ direction. The above equations then form a pair of integral equations which must be simultaneously solved. 11
12 11 Numerical Solutions Separation of variables provides an analytic solution when the boundaries of the potential surfaces are the same as those obtained by choosing each variable of the separation geometry as a constant in order to project out a surface in the coordinate system. Of course Laplace s equation also must separate into ode equations each involving only one of these variables. While analytic solutions provide insight into more realistic problems, they are not always useful in obtaining detailed information. Thus techniques to obtain accurate numerical solution of Laplace s (and Poisson s) equation are necessary. First consider a result of Gauss theorem. Integrate Laplace s equation over the volume which bounds the potential surface. dτ 2 V = V d σ = In the above, σ is the surface which encloses the volume τ. In the case of a spherical surface, d σ = R 2 dω ˆr. Substitute this in the above to write; R 2 d dr dω V = This equation means that dωv is constant and. In Cartesian coordinates divide space into a grid with cells of the dimensions (δx, δy, δz). From the above analysis, the potential at the center of a cell will approximately equal the average of the potential over the enclosing surface. In Cartesian coordinates the potential at a point is approximately the average of the sum of the potentials over its nearest neighbors. V l,m,n = 1/6[V l 1,m,n + V l+1,m,n + V l,m 1,n + V l,m+1,n + V l,m 1,n + V l,m,n 1 + V l,m,n+1 ] Begin a numerical solution by taking the exact potential values on the surface and assign initial values to the interior potential at all grid points. These initial values can be any guess. The average values at each point are then obtained, but one keeps the correct potential values on the surface. The process is then iterated to convergence. This procedure technique is called the relaxation method. It is stable with respect to iteration and converges rapidly to the correct potential inside a volume. The technique (finite element analysis) is generally applied to any process which is described by Laplace s equation. This includes a number of physical problems in addition to electrostatics. If charge is present, a solution to Poisson s equation is required. For a sphere of radius, r, the potential at the center relative to the surface is; V = ρr 2 /(6ǫ ) 12
13 Figure 8: An 3-D numerical example showing contour lines of constant potential of a geometry having grounded metal boundaries and wires held at potential This is included in the equation above when computing the average. As an example, Figure 8 shows a numerical valuation of the potential of a set of grounded metal boundaries and wires which are held at constant potential. 13
Bessel Functions - Lecture 7
1 Introduction We study the ode; Bessel Functions - Lecture 7 x 2 f ν + xf ν + (x2 ν 2 )f ν = This is a Sturm-Liouville problem where we look for solutions as the variable ν is changed. The equation has
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 information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationElectrodynamics I Midterm - Part A - Closed Book KSU 2005/10/17 Electro Dynamic
Electrodynamics I Midterm - Part A - Closed Book KSU 5//7 Name Electro Dynamic. () Write Gauss Law in differential form. E( r) =ρ( r)/ɛ, or D = ρ, E= electricfield,ρ=volume charge density, ɛ =permittivity
More informationAppendix A Vector Analysis
Appendix A Vector Analysis A.1 Orthogonal Coordinate Systems A.1.1 Cartesian (Rectangular Coordinate System The unit vectors are denoted by x, ŷ, ẑ in the Cartesian system. By convention, ( x, ŷ, ẑ triplet
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationBoundary-Value Problems: Part II
Chapter 3 Boundary-Value Problems: Part II Problem Set #3: 3., 3.3, 3.7 Due Monday March. th 3. Spherical Coordinates Spherical coordinates are used when boundary conditions have spherical symmetry. The
More informationSurvival Guide to Bessel Functions
Survival Guide to Bessel Functions December, 13 1 The Problem (Original by Mike Herman; edits and additions by Paul A.) For cylindrical boundary conditions, Laplace s equation is: [ 1 s Φ ] + 1 Φ s s s
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
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 informationMATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued).
MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued). Bessel s differential equation of order m 0: z 2 d2 f dz 2 + z df dz + (z2 m 2 )f = 0 The equation is considered on the
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 informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
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 informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationPHY752, Fall 2016, Assigned Problems
PHY752, Fall 26, Assigned Problems For clarification or to point out a typo (or worse! please send email to curtright@miami.edu [] Find the URL for the course webpage and email it to curtright@miami.edu
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More information(a) Determine the general solution for φ(ρ) near ρ = 0 for arbitary values E. (b) Show that the regular solution at ρ = 0 has the series expansion
Problem 1. Curious Wave Functions The eigenfunctions of a D7 brane in a curved geometry lead to the following eigenvalue equation of the Sturm Liouville type ρ ρ 3 ρ φ n (ρ) = E n w(ρ)φ n (ρ) w(ρ) = where
More informationChapter 4 Boundary Value Problems in Spherical and Cylindrical Coordinates
Chapter 4 Boundary Value Problems in Spherical and Cylindrical Coordinates 4. Laplace s Equation in spherical coordinates Laplace s equation in spherical coordinates, (r, θ, φ), has the form where r r
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More information1 Bessel, Neumann, and Hankel Functions: J ν (x), N ν (x), H (1)
Bessel, Neumann, and Hankel Functions: J ν x), N ν x), H ) ν x), H ) ν x) Bessel functions are solutions of the following differential equation: x y +xy +x ν )y =.) Any two of the following functions are
More informationPARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS
PARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS NAKHLE H. ASMAR University of Missouri PRENTICE HALL, Upper Saddle River, New Jersey 07458 Contents Preface vii A Preview of Applications and
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationModified Bessel functions : Iα, Kα
Modified Bessel functions : Iα, Kα The Bessel functions are valid even for complex arguments x, and an important special case is that of a purely imaginary argument. In this case, the solutions to the
More informationCHAPTER 3 POTENTIALS 10/13/2016. Outlines. 1. Laplace s equation. 2. The Method of Images. 3. Separation of Variables. 4. Multipole Expansion
CHAPTER 3 POTENTIALS Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 Outlines 1. Laplace s equation 2. The Method of Images 3. Separation of Variables 4. Multipole Expansion
More informationSteady and unsteady diffusion
Chapter 5 Steady and unsteady diffusion In this chapter, we solve the diffusion and forced convection equations, in which it is necessary to evaluate the temperature or concentration fields when the velocity
More informationFOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS
fc FOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS Second Edition J. RAY HANNA Professor Emeritus University of Wyoming Laramie, Wyoming JOHN H. ROWLAND Department of Mathematics and Department
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationPHYS 502 Lecture 8: Legendre Functions. Dr. Vasileios Lempesis
PHYS 502 Lecture 8: Legendre Functions Dr. Vasileios Lempesis Introduction Legendre functions or Legendre polynomials are the solutions of Legendre s differential equation that appear when we separate
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More information(b) Show that the charged induced on the hemisphere is: Q = E o a 2 3π (1)
Problem. Defects This problem will study defects in parallel plate capacitors. A parallel plate capacitor has area, A, and separation, D, and is maintained at the potential difference, V = E o D. There
More informationMore on Bessel functions. Infinite domain, δ-function normalization
More on Bessel functions Infinite domain, δ-function normalization Consider Bessel s equation on the domain < ρ < as R ½ Bessel s equation, (375) or (393), says ρ d ρ dj ν(kρ) + k 2 ν 2 ρ 2 J ν (kρ) =
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 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 information1 Introduction. Green s function notes 2018
Green s function notes 8 Introduction Back in the "formal" notes, we derived the potential in terms of the Green s function. Dirichlet problem: Equation (7) in "formal" notes is Φ () Z ( ) ( ) 3 Z Φ (
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General Friday, 30 May, 2014 9:00 am to 12:00 pm MATHEMATICS (2 Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationFields of a Dipole Near a Layered Substrate
Appendix C Fields of a Dipole Near a Layered Substrate µ z θ µ 1 ε 1 µ 2 ε 2 µ 3 ε 3 d z o x,y Figure C.1: An electric dipole with moment µ is located at r o = (,, z o ) near a layered substrate. The fields
More information3 Chapter. Gauss s Law
3 Chapter Gauss s Law 3.1 Electric Flux... 3-2 3.2 Gauss s Law (see also Gauss s Law Simulation in Section 3.10)... 3-4 Example 3.1: Infinitely Long Rod of Uniform Charge Density... 3-9 Example 3.2: Infinite
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 informationSOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4
SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with
More informationFinal Examination Solutions
Math. 42 Fulling) December 22 Final Examination Solutions Calculators may be used for simple arithmetic operations only! Laplacian operator in polar coordinates: Some possibly useful information 2 u =
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationl=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.
Electrodynamics I Exam - Part A - Closed Book KSU 15/11/6 Name Electrodynamic Score = 14 / 14 points Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try
More informationProblem 1. Arfken
All Arfken problems are reproduced at the end of this assignment Problem 1. Arfken 11.2.11. Problem 2. (a) Do Arfken 11.3.3 (b) The integral Integrals 4 3i 3+4i (x 2 iy 2 ) dz (1) is not defined without
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 informationSeparation of Variables in Polar and Spherical Coordinates
Separation of Variables in Polar and Spherical Coordinates Polar Coordinates Suppose we are given the potential on the inside surface of an infinitely long cylindrical cavity, and we want to find the potential
More informationReview of Electrostatics
Review of Electrostatics 1 Gradient Define the gradient operation on a field F = F(x, y, z) by; F = ˆx F x + ŷ F y + ẑ F z This operation forms a vector as may be shown by its transformation properties
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 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 information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
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 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 informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General) Friday, 31 May, 2013 9:00 am to 12:00 pm MATHEMATICS (2) Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More 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 informationReview of Electrostatics. Define the gradient operation on a field F = F(x, y, z) by;
Review of Electrostatics 1 Gradient Define the gradient operation on a field F = F(x, y, z) by; F = ˆx F x + ŷ F y + ẑ F z This operation forms a vector as may be shown by its transformation properties
More informationPhysics 505 Fall 2003
Physics 505 Fall 2003 G. Raithel December 8, 2003 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationMathematical Methods for Physics
Mathematical Methods for Physics Peter S. Riseborough June 8, 8 Contents Mathematics and Physics 5 Vector Analysis 6. Vectors................................ 6. Scalar Products............................
More informationIntegrals in cylindrical, spherical coordinates (Sect. 15.7)
Integrals in clindrical, spherical coordinates (Sect. 15.7 Integration in spherical coordinates. Review: Clindrical coordinates. Spherical coordinates in space. Triple integral in spherical coordinates.
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 informationPhysics 505 Fall 2005 Homework Assignment #8 Solutions
Physics 55 Fall 5 Homework Assignment #8 Solutions Textbook problems: Ch. 5: 5., 5.4, 5.7, 5.9 5. A circular current loop of radius a carrying a current I lies in the x-y plane with its center at the origin.
More informationPhysics 505 Fall Homework Assignment #9 Solutions
Physics 55 Fall 25 Textbook problems: Ch. 5: 5.2, 5.22, 5.26 Ch. 6: 6.1 Homework Assignment #9 olutions 5.2 a) tarting from the force equation (5.12) and the fact that a magnetization M inside a volume
More informationFinal Examination Solutions
Math. 42 Fulling) 4 December 25 Final Examination Solutions Calculators may be used for simple arithmetic operations only! Laplacian operator in polar coordinates: Some possibly useful information 2 u
More informationHEAT CONDUCTION USING GREEN S FUNCTIONS
HEAT CONDUCTION USING GREEN S FUNCTIONS Preface to the first edition Preface to the second edition Author Biographies Nomenclature TABLE OF CONTENTS FOR SECOND EDITION December 2009 Page viii x xii xiii
More informationLaplace s equation in polar coordinates. Boundary value problem for disk: u = u rr + u r r. r 2
Laplace s equation in polar coordinates Boundary value problem for disk: u = u rr + u r r + u θθ = 0, u(a, θ) = h(θ). r 2 Laplace s equation in polar coordinates Boundary value problem for disk: u = u
More informationSolutions to PS 2 Physics 201
Solutions to PS Physics 1 1. ke dq E = i (1) r = i = i k eλ = i k eλ = i k eλ k e λ xdx () (x x) (x x )dx (x x ) + x dx () (x x ) x ln + x x + x x (4) x + x ln + x (5) x + x To find the field for x, we
More informationSummary: Applications of Gauss Law
Physics 2460 Electricity and Magnetism I, Fall 2006, Lecture 15 1 Summary: Applications of Gauss Law 1. Field outside of a uniformly charged sphere of radius a: 2. An infinite, uniformly charged plane
More informationPhysics 116C Solutions to the Practice Final Exam Fall The method of Frobenius fails because it works only for Fuchsian equations, of the type
Physics 6C Solutions to the Practice Final Exam Fall 2 Consider the differential equation x y +x(x+)y y = () (a) Explain why the method of Frobenius method fails for this problem. The method of Frobenius
More informationONE AND MANY ELECTRON ATOMS Chapter 15
See Week 8 lecture notes. This is exactly the same as the Hamiltonian for nonrigid rotation. In Week 8 lecture notes it was shown that this is the operator for Lˆ 2, the square of the angular momentum.
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More 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 information1 A complete Fourier series solution
Math 128 Notes 13 In this last set of notes I will try to tie up some loose ends. 1 A complete Fourier series solution First here is an example of the full solution of a pde by Fourier series. Consider
More informationMULTIVARIABLE INTEGRATION
MULTIVARIABLE INTEGRATION (SPHERICAL POLAR COORDINATES) Question 1 a) Determine with the aid of a diagram an expression for the volume element in r, θ, ϕ. spherical polar coordinates, ( ) [You may not
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (I)
Laplace s Equation in Cylindrical Coordinates and Bessel s Equation I) 1 Solution by separation of variables Laplace s equation is a key equation in Mathematical Physics. Several phenomena involving scalar
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More information0.4. math. r r = r 2 + r 2 2r r cos(θ ) 0.5. Coordiates.
.. Electrostatics... Multipole. F = E(r) = E(r) = S n qq i r r i 3 (r r i) i= n i= q i r r i 3 (r r i) d 3 r ρ(r ) 3 (r r ) ds E(r) = 4πQ enclosed E(r) = 4πρ(r ) E(r) = Φ(r) d 3 r ρ(r ) r dl E W = q(φ(r
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationProblems for MATH-6300 Complex Analysis
Problems for MATH-63 Complex Analysis Gregor Kovačič December, 7 This list will change as the semester goes on. Please make sure you always have the newest version of it.. Prove the following Theorem For
More information# Points Score Total 100
Name: PennID: Math 241 Make-Up Final Exam January 19, 2016 Instructions: Turn off and put away your cell phone. Please write your Name and PennID on the top of this page. Please sign and date the pledge
More informationLecture 4.6: Some special orthogonal functions
Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics
More informationSyllabus: for Complex variables
EE-2020, Spring 2009 p. 1/42 Syllabus: for omplex variables 1. Midterm, (4/27). 2. Introduction to Numerical PDE (4/30): [Ref.num]. 3. omplex variables: [Textbook]h.13-h.18. omplex numbers and functions,
More informationElectric fields in matter
Electric fields in matter November 2, 25 Suppose we apply a constant electric field to a block of material. Then the charges that make up the matter are no longer in equilibrium: the electrons tend to
More informationSolutions to Problems in Jackson, Classical Electrodynamics, Third Edition. Chapter 2: Problems 11-20
Solutions to Problems in Jackson, Classical Electrodynamics, Third Edition Homer Reid December 8, 999 Chapter : Problems - Problem A line charge with linear charge density τ is placed parallel to, and
More informationxy 2 e 2z dx dy dz = 8 3 (1 e 4 ) = 2.62 mc. 12 x2 y 3 e 2z 2 m 2 m 2 m Figure P4.1: Cube of Problem 4.1.
Problem 4.1 A cube m on a side is located in the first octant in a Cartesian coordinate system, with one of its corners at the origin. Find the total charge contained in the cube if the charge density
More informationThe Divergence Theorem
Math 1a The Divergence Theorem 1. Parameterize the boundary of each of the following with positive orientation. (a) The solid x + 4y + 9z 36. (b) The solid x + y z 9. (c) The solid consisting of all points
More informationMA 351 Fall 2008 Exam #3 Review Solutions 1. (2) = λ = x 2y OR x = y = 0. = y = x 2y (2x + 2) = 2x2 + 2x 2y = 2y 2 = 2x 2 + 2x = y 2 = x 2 + x
MA 5 Fall 8 Eam # Review Solutions. Find the maimum of f, y y restricted to the curve + + y. Give both the coordinates of the point and the value of f. f, y y g, y + + y f < y, > g < +, y > solve y λ +
More informationOn an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University
On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University njrose@math.ncsu.edu 1. INTRODUCTION. The classical eigenvalue problem for the Legendre Polynomials
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationS. Ghorai 1. Lecture XV Bessel s equation, Bessel s function. e t t p 1 dt, p > 0. (1)
S Ghorai 1 1 Gamma function Gamma function is defined by Lecture XV Bessel s equation, Bessel s function Γp) = e t t p 1 dt, p > 1) The integral in 1) is convergent that can be proved easily Some special
More informationSummary for Vector Calculus and Complex Calculus (Math 321) By Lei Li
Summary for Vector alculus and omplex alculus (Math 321) By Lei Li 1 Vector alculus 1.1 Parametrization urves, surfaces, or volumes can be parametrized. Below, I ll talk about 3D case. Suppose we use e
More informationSections minutes. 5 to 10 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed.
MTH 34 Review for Exam 4 ections 16.1-16.8. 5 minutes. 5 to 1 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed. Review for Exam 4 (16.1) Line
More informationMath 337, Summer 2010 Assignment 5
Math 337, Summer Assignment 5 Dr. T Hillen, University of Alberta Exercise.. Consider Laplace s equation r r r u + u r r θ = in a semi-circular disk of radius a centered at the origin with boundary conditions
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 informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More information