Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

Similar documents
Physics 505 Fall Homework Assignment #4 Solutions

Legendre Polynomials - Lecture 8

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Separation of Variables and a Spherical Shell with Surface Charge

$, (2.1) n="# #. (2.2)

Notes: 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.

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

Lecture 8 February 18, 2010

Physics 506 Winter 2006 Homework Assignment #6 Solutions

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions

Math 124B January 31, 2012

221B Lecture Notes Notes on Spherical Bessel Functions

Section 6: Magnetostatics

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

11.1 One-dimensional Helmholtz Equation

Higher dimensional PDEs and multidimensional eigenvalue problems

Assignment 7 Due Tuessday, March 29, 2016

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

Radiation Fields. Lecture 12

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

4 1-D Boundary Value Problems Heat Equation

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

Lecture 17 - The Secrets we have Swept Under the Rug

In Coulomb gauge, the vector potential is then given by

Cluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects.

Week 6 Lectures, Math 6451, Tanveer

David Eigen. MA112 Final Paper. May 10, 2002

14 Separation of Variables Method

6 Wave Equation on an Interval: Separation of Variables

Physics 566: Quantum Optics Quantization of the Electromagnetic Field

Gauss s law - plane symmetry

Week 5 Lectures, Math 6451, Tanveer

2M2. Fourier Series Prof Bill Lionheart

More Scattering: the Partial Wave Expansion

Math 124B January 17, 2012

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection

Lecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series

Chapter 4: Electrostatic Fields in Matter

arxiv: v1 [hep-th] 10 Dec 2018

9. EXERCISES ON THE FINITE-ELEMENT METHOD

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Course 2BA1, Section 11: Periodic Functions and Fourier Series

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

CS229 Lecture notes. Andrew Ng

Module 22: Simple Harmonic Oscillation and Torque

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers

MAT 167: Advanced Linear Algebra

C. Fourier Sine Series Overview

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

Jackson 3.3 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Wave Equation Dirichlet Boundary Conditions

1 Heat Equation Dirichlet Boundary Conditions

THE NUMERICAL EVALUATION OF THE LEVITATION FORCE IN A HYDROSTATIC BEARING WITH ALTERNATING POLES

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Math 220B - Summer 2003 Homework 1 Solutions

Chemical Kinetics Part 2. Chapter 16

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...

PhysicsAndMathsTutor.com

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

4 Separation of Variables

ANISOTROPIES OF THE MICROWAVE BACKGROUND

Dual Integral Equations and Singular Integral. Equations for Helmholtz Equation

An approximate method for solving the inverse scattering problem with fixed-energy data

Chemical Kinetics Part 2

Physics 505 Fall 2005 Homework Assignment #8 Solutions

Electromagnetism Spring 2018, NYU

MA 201: Partial Differential Equations Lecture - 10

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE

APPENDIX C FLEXING OF LENGTH BARS

Why Doesn t a Steady Current Loop Radiate?

Chapter 5. Wave equation. 5.1 Physical derivation

Term Test AER301F. Dynamics. 5 November The value of each question is indicated in the table opposite.

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Lecture Note 3: Stationary Iterative Methods

arxiv: v1 [math.ca] 6 Mar 2017

Theoretical Cosmology

APPENDIX B. Some special functions in low-frequency seismology. 2l +1 x j l (B.2) j l = j l 1 l +1 x j l (B.3) j l = l x j l j l+1

Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V

1.2 Partial Wave Analysis

ESCI 340 Physical Meteorology Radiation Lesson 5 Terrestrial Radiation and Radiation Balance Dr. DeCaria

Agenda Administrative Matters Atomic Physics Molecules

Volume 13, MAIN ARTICLES

Introduction to LMTO method

c 2007 Society for Industrial and Applied Mathematics

1D Heat Propagation Problems

1. Measurements and error calculus

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

Extraction of Astrophysical Cross Sections in the Trojan-Horse Method

b n n=1 a n cos nx (3) n=1

Elements of Kinetic Theory

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

1.1 a.) Suppose we have a conductor and place some charge density within it. For a surface S inside the conductor enclose the charge density!

MA 201: Partial Differential Equations Lecture - 11

Transcription:

Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that the answer otained in this way agrees with the direct soution from the differentia equation. Reca that Proem 3. asks for the potentia etween two concentric spheres of radii a and with > a, where the upper hemisphere of the inner sphere and the ower hemisphere of the outer sphere are maintained at potentia V, and where the other hemispheres are at zero potentia. Since this proem invoves the potentia etween two spheres, we use the Dirichet Green s function G x, x = +,m a + r< a+ r + < r> + r > + Y m ΩY m Ω Because there are no charges etween the spheres, the Green s function soution for the potentia ony invoves the surface integra Φ x = S Φ x G n da Here, we note a sutety in that the oundary surface is actuay disconnected, and incudes oth the inner sphere of radius a and the outer sphere of radius. This means that the potentia may e written as a sum of two contriutions Φ x = r =a Φ x G n a dω r = We now compute the norma derivatives of the Green s function and G r n = G r< =a r = =r =a = a a,m a + r + a Φ x G n dω + r Y m ΩY m Ω G r n = G r> = r = =r = = r a a +,m a + Y m ΩY m Ω r

Inserting these expressions into yieds Φ x =,m +,m V a Ω Y m Ω dω V Ω Y m Ω dω a a + a + r + r r a a + + r a Y m Ω Y m Ω This is the genera expression for the soution to the oundary vaue proem where V a Ω is the potentia on the inner sphere and V Ω is the potentia on the outer sphere. For the upper/ower hemispheres proem, we note that azimutha symmetry aows us to restrict the m vaues to m = ony. In this case, the spherica harmonic expansion reduces to a Legendre poynomia expansion Φ x = + V a ζp ζdζ V ζp ζdζ + a + a + r + r a + r a a + r a P cos θ + P cos θ where ζ = cos θ. Since V a = V for ζ > and V = V for ζ <, this aove expression reduces to Φ x = + V N a + a + r a + r + + V N r a a a + r = + V N a + a + r a + r r a + P cos θ + P cos θ a r + P cos θ where N = { = P ζdζ = j+ Γj πj! = j odd

If desired, the potentia may e rearranged to read Φ x = + V N a + + + a a + r + + + a + r = V + V j+ 4j Γj a j a j j= 4 πj! a 4j + r P cos θ a j r j + P j cos θ which agrees with the soution to Proem 3. that we have found earier. 3.7 The Dirichet Green function for the unounded space etween the panes at z = and z = L aows discussion of a point charge or a distriution of charge etween parae conducting panes hed at zero potentia. a Using cyindrica coordinates show that one form of the Green function is G x, x = 4 L n= m= e imφ φ sin nπz sin L nπz L nπ nπ I m L ρ < K m L ρ > In cyindrica coordinates, the poar direction φ is periodic with period π. This suggests that the Green s function coud e expanded as a Fourier series in e imφ. Simiary, the oundary conditions G = at z = and z = L motivates the use of a Fourier sine series sinnπz/l in the z coordinate. More precisey, a compete Fourier expansion in φ and z woud give G x, x = gρ, ρ e imφ e im φ nπz n πz sin sin L L m,n,m,n However, it turns out that m and m and n and n do not need to e chosen to e independent. This can e seen from the Green s function equation given here as a differentia equation in x x G x, x = δ 3 x x In cyindrica coordinates, this reads ρ ρ ρ ρ + ρ φ + z Gρ, φ, z; ρ, φ, z = ρ δρ ρ δφ φ δz z 3

Using the competeness reations and suggests that we take n= m= e imφ φ = πδφ φ 4 nπz nπz sin sin L L = L δz z G x, x = n= m= gρ, ρ e imφ φ sin nπz nπz sin L L 5 Sustituting this decomposition into 3 gives d ρ dρ ρ d dρ m ρ Making the sustitution nπ L gρ, ρ = 4 Lρ δρ, ρ x = nπρ L converts the homogeneous part of this to a modified Besse equation d dx + x d dx + m x gx, x = 4 Lx δx, x At this stage, the soution ecomes standard. Noting that the modified Besse function I m x ows up as x and the function K m x ows up as x, we are eft with { gx, x AIm x x < x = BK m x x > x where the coefficients A and B are determined y the matching conditions g < = g >, d dx g < = d dx g > + 4 Lx at x = x. This system may e soved to yied A = 4 K m x Lx I mx K m x I m x K mx B = 4 I m x Lx I mx K m x I m x K mx

Noting that the modified Besse functions satisfy the Wronskian formua I ν xk νx I νxk ν x = x finay gives { Im xk m x x < x gx, x = 4 L I m x K m x x > x = 4 L I mx < K m x > where x < = minx, x, x > = maxx, x Converting x ack to ρ and sustituting into 5 then gives the desired Dirichet Green s function G x, x = 4 L n= m= e imφ φ sin nπz nπz sin L L I m nπρ< L nπρ> K m L Show that an aternative form of the Green function is G x, x = m= dk e imφ φ J m kρj m kρ sinhkz < sinhkl z > sinhkl This aternative form of the Green s function is derived y expanding in φ and ρ instead of φ and z. For the ρ expansion, we use the integra reation kj ν kρj ν kρ dk = ρ δρ ρ aong with the competeness reation 4 to motivate the decomposition G x, x = m= kdk g k z, z e imφ φ J m kρj m kρ 6 Since the Besse function J m kρ satisfies the Besse equation d dρ + ρ d dρ + k m ρ J m kρ = the sustitution of 6 into the Greens function equation 3 gives d dz k g k z, z = δz z

Since g k z, z vanishes at z = and z = L, this is a standard one-dimensiona Green s function proem. Writing { g k z, z A sinhkz z < z = B sinhkl z z > z we find that the matching and jump conditions ecome A sinhkz = B sinhkl z, A coshkz = B coshkl z + k This may e soved to give so that A = k g k z, z = Sustituting this into 6 then yieds G x, x = m= k sinhkl z, B = sinhkz sinhkl k sinhkl k sinhkl sinhkz < sinhkl z > dk e imφ φ J m kρj m kρ sinhkz < sinhkl z > sinhkl 3.6 Consider the Green function appropriate for Neumann oundary conditions for the voume V etween the concentric spherica surfaces defined y r = a and r =, a <. To e ae to use.46 for the potentia, impose the simpe constraint.45. Use an expansion in spherica harmonics of the form G x, x = g r, r P cos γ = where g r, r = r </r + > + f r, r. a Show that for >, the radia Green function has the symmetric form g r, r = r < r> + + + + a + rr + + a + r r rr + a+ + + r + r + There are severa approaches to this proem. However, we first consider the Neumann oundary condition.45 G x x n = ndy S

For this proem with two oundaries, the surface area S must e the area of oth oundaries ie it is the tota area surrounding the voume. Hence S = a +, and in particuar this is uniform constant in the anges. As a resut, this wi ony contriute to the = term in the expansion of the Green s function. More precisey, we coud write G x x n = ndy g r, r n P cos γ ndy = a + P cos γ Since the Legendre poynomias are orthogona, this impies that g r, r n = ndy a + δ, Noting that the outward norma is either in the ˆr or the ˆr direction for the sphere at a or, respectivey, we end up with two oundary condition equations g r, r r = a a + δ g r, r, r = a + δ, 7 Now that we have written down the oundary conditions for g r, r, we proceed to otain its expicit form. The suggestion of the proem is to write Since g r, r = r < r + > x x = + f r, r r< r> + P cos γ we see that the first term in g r, r is designed to give the singuar source deta function. The remaining term F x, x = f r, r P cos γ then soves the homogeneous equation x F x, x =. But we know how to sove Lapace s equation in spherica coordinates, and the resut is that the radia function must e of the form f r, r = A r + B r + Note that we are taking the Green s function equation to act on the x variae, where x may e thought of as a parameter constant giving the ocation of the deta function source. We thus have g r, r = r < r + > + A r + B r + 8

A that remains is to use the oundary conditions 7 to sove for A and B. For the inside sphere at a, we have whie for the outside sphere we have a r + + A a + B a + = δ, a + 9 + r + + A + B + = δ, a + For we rewrite these equations as a + + A a + = + /r + + B + r which may e soved to give A = B + + a + r = + a + Inserting this into 8 yieds g r, r = r < r + > + r + a + = r < r> + + + + a + rr + + = + a + = + + a + r< + + + a + /r + + a + + r a/r + + + / a + + / + a/r + a + + + r + a + + + r < r > + + + a + r + < + r a + r r + a + r rr + a+ + r + + a + r < r > + + + r < r + > r > + + + a + r > + r< + r r + which is vaid for. Note that in the ast few ines we have een ae to rewrite the Green s function in terms of a product of ur < and vr > where u and v satisfies Neumann oundary conditions at r = a and r =, respectivey. r + >

This is reated to another possie method of soving this proem. Using the Legendre identity = + P cos γ = δφ φ δcos θ cos θ the Green s function equation may e reduced to the one-dimensiona proem d d dr r + g dr r, r = + δr r Using the genera method for the Sturm-Liouvie proem, the Green s function is given y g r, r = + A u r < v r > where ur and vr sove the homogeneous equation and the constant A is fixed y the Wronskian, W u, v = A /r. For the oundary conditions 7 are homogeneous u r r =a = v r r = = It is easy to see that these are satisfied y ur = r + Computing the Wronskian gives which aows us to identify a + + r + vr = r + + + r + u v u v = + a+ + + r A = + + + a + This gives the resut of the ast ine of. Show that for = g r, r = a r > a + r + fr where fr is aritrary. Show expicity in.46 that answers for the potentia Φ x are independent of fr. The = case invoves a non-homogeneous oundary condition. Hence the resut of wi not work. Of course, we can sti work out the one-dimensiona deta

function proem with matching and jump conditions at r = r. However it is more direct to return to 9 and and to simpy sove those conditions for =. Both 9 and resut in a B = a + whie eaving A competey undetermined. Finay, since r is thought of as a parameter, this indicates that A = fr can e an aritrary function of r. The = Green s function is given y 8 g r, r = r > a a + r + fr Incidentay, we note that without the inhomogeneous Neumann oundary condition term /S there wi e no soution to the system 9 and for = uness is taken to. This demonstrates the inconsistency of simpy setting G/ n = for the Neumann Green s function. Note that, y setting fr = a /a + r we otain a symmetrica Green s function g r, r = a r > a + r + r On the other hand, the choice of fr is unphysica. This arises ecause, for the Neumann Green s function, the fr contriution to the potentia is given y Φ x = ρ x fr d 3 x + Φ x ɛ V S n fr da = fr ρ x d 3 x ɛ E x dâ ɛ V S = fr q enc ɛ E x dâ = ɛ y Gauss aw. It is important not to mix up r and r in this derivation. 3.7 Appy the Neumann Green function of Proem 3.6 to the situation in which the norma eectric fied is E r = E cos θ at the outer surface r = and is E r = on the inner surface r = a. a Show that the eectrostatic potentia inside the voume V is r cos θ Φ x = E p 3 + a3 r 3 where p = a/. Find the components of the eectric fied cos θ E r r, θ = E p 3 a3 sin θ r 3, E θ r, θ = E p 3 S + a3 r 3

Since there is no charge etween the spheres, the soution to e oundary vaue proem is given y Φ x = Φ x S n G x, x da = E r Ω G x, x dω r = = E G x, x cos θ dω By writing = E r = = P cos γ = r = + g r, r P cos γ cos θ dω m Y m ΩY m Ω and noting that cos θ = /3Y Ω, we end up with the expansion Φ x = E 3,m g r, Y m Ω + = E g r, Y Ω 3 3 = E cos θ g r, 3 Y m Ω Y Ω dω where we have used orthogonaity of the spherica harmonics. Inserting = into then gives Φ x = E cos θ r 3 3 a 3 + a3 3 r = E r cos θ a/ 3 + a3 r 3 3 This is the potentia for a constant eectric fied comined with an eectric dipoe. Defining p = a/, the components of the eectric fied are E r = Φ r = E cos θ p 3 a3 r 3, E θ = Φ r θ = E sin θ p 3 + a3 r 3 Note that the oundary conditions E r r=a = and E r r= = E cos θ are oviousy satisfied. On the other hand, the parae component of the fied, E θ, is non-vanishing on oth surfaces except at the poes. Physicay, this indicates that these surfaces are not conductors. Cacuate the Cartesian or cyindrica components of the fied, E z ad E ρ, and make a sketch or computer pot of the ines of eectric force for a typica case of p =.5.

Rewriting 3 as Φ x = E p 3 z + a3 z r 3 = E p 3 z + a 3 z ρ + z 3/ we otain E z = Φ z = E p 3 + a3 3ẑ r 3 E ρ = Φ ρ = E p 3 3a3 ẑ ˆρ r 3 where ˆρ = ρ/r, ẑ = z/r, and r = ρ + z. As indicated aove, this corresponds to a constant eectric fied comined with an eectric dipoe. For E >, a sketch of the eectric fied ines ooks ike - - - - This sketch indicates that the radia component of the eectric fied vanishes at the surface of the inner sphere.