f07 practice midterm solutions

Upload: kapil-adhikari

Post on 10-Apr-2018

230 views

Category:

Documents


1 download

TRANSCRIPT

  • 8/8/2019 F07 Practice Midterm Solutions

    1/8

    Physics 505 Fall 2007

    Practice Midterm Solutions

    This midterm will be a two hour open book, open notes exam. Do all three problems.

    1. A rectangular box has sides of lengths a, b and c

    c

    x

    y

    z

    a

    b

    a) For the Dirichlet problem in the interior of the box, the Greens function may be

    expanded as

    G(x,y,z ; x , y , z ) =

    m =1

    n =1

    gmn (z, z )sinmx

    asin

    mxa

    sinny

    bsin

    nyb

    Write down the appropriate differential equation that gmn (z, z ) must satisfy.

    Note that sin kx satises the completeness relation

    m =1

    sinmx

    asin

    mxa

    =a2

    (x x )

    Hence the Greens function equation

    2x G( x, x ) = 4

    3 ( x x )has an expansion

    m,n

    2x gmn (z, z )sin

    mxa

    sinmx

    asin

    nyb

    sinny

    b

    = 4(z z )4ab

    m,n

    sinmx

    asin

    mxa

    sinny

    bsin

    nyb

    Working out the x and y derivatives on the left-hand side yields

    m,n

    d2

    dz 2 ma

    2

    nb

    2

    gmn (z, z )sinmx

    asin

    mxa

    sinny

    bsin

    nyb

    = 16ab

    (z z )m,n

    sinmx

    asin

    mxa

    sinny

    bsin

    nyb

  • 8/8/2019 F07 Practice Midterm Solutions

    2/8

    However, since sin kx forms an orthogonal basis, each term in this sum mustvanish by itself. This results in the differential equation

    d2

    dz 2 2mn gmn (z, z ) =

    16ab

    (z z ) (1)where mn = (m/a )

    2 + ( n/b )2 is given in part c). Note that the Fourier sineexpansion automatically satises Dirichlet boundary conditions for x and y. Theremaining boundary condition is that gmn (z, z ) vanishes whenever z or z is equalto 0 or c.

    b) Solve the Greens function equation for gmn (z, z ) subject to Dirichlet boundaryconditions and write down the result for G(x,y,z ; x , y , z ).

    We may solve the Greens function equation (1) by rst noting that the homoge-neous equation is of the form

    g (z )

    2

    mng(z ) = 0

    This is a second-order linear equation with constant coefficients admitting thefamiliar solution

    g(z ) = Ae mn z + Be mn z

    However, we want g(z ) = 0 when z = 0 or z = c. This motivates us to writeout the solutions

    u(z ) = sinh mn z 0 < z < zv(z ) = sinh[ mn (c z )] z < z < c

    The Greens function solution is then given by

    gmn (z, z ) =Au(z ) z < zBv (z ) z > z

    The matching conditions

    g< = g> ,d

    dzg< =

    ddz

    g> +16ab

    then give the system

    A sinh mn z = B sinh[ mn (c z)],A cosh mn z = B cosh[ mn (c z)] +

    16ab mn

    which may be solved to yield

    A =16 sinh[ mn (c z)]

    ab mn (cosh mn z sinh[ mn (c z)] + sinh mn z cosh[ mn (c z)])=

    16 sinh[ mn (c z)]ab mn sinh mn c

    =16

    ab mn sinh mn cv(z)

  • 8/8/2019 F07 Practice Midterm Solutions

    3/8

    and

    B =16 sinh mn z

    ab mn (cosh mn z sinh[ mn (c z)] + sinh mn z cosh[ mn (c z)])=

    16 sinh mn zab mn sinh mn c =

    16ab mn sinh mn cu(z)

    As a result, the full Greens function solution is then given by

    gmn (z, z ) =16

    ab mn sinh mn cu(z< )v(z> )

    =16

    ab mn sinh mn csinh mn z< sinh mn [(c z> )]

    Hence

    G( x, x ) =16ab

    m,n

    1 mn sinh mn c

    sinmx

    asin

    mxa

    sinny

    bsin

    nyb

    sinh mn z< sinh mn [(c z> )](2)

    c) Consider the boundary value problem where the potential on top of the box is(x,y,c ) = V (x, y ) while the potential on the other ve sides vanish. Using theGreens function obtained above, show that the potential may be written as

    (x,y,z ) =

    m =1

    n =1

    Amn sinmx

    asin

    nyb

    sinh mn z

    where mn = (m/a )2 + ( n/b )2 andAmn =

    4ab sinh mn c

    a

    0dx

    b

    0dy V (x, y )sin

    mxa

    sinny

    b

    Since we only have to worry about the potential on the top of the box (and sincewe assume there is no charge inside the box), the Greens function solution maybe written

    ( x ) = 1

    4 z = c (x ) G ( x, x )n da=

    14 z = c V (x , y ) G ( x, x )n dz dy

    (3)

  • 8/8/2019 F07 Practice Midterm Solutions

    4/8

    Noting that the outward-pointing normal n on the top of the box is in the + zdirection, we compute the normal derivative of (2)

    G ( x, x )n z = c

    =G ( x, x )

    z z = c

    = 16ab

    m,n

    1 mn sinh mn c

    sin mxa

    sin mxa

    sin nyb

    sin nyb

    sinh mn z mn cosh mn [(c z )] z = c= 4

    m,n

    4ab sinh mn c

    sinmx

    asin

    nyb

    sinh mn z

    sinmx

    asin

    nyb

    Inserting this into (3) then straightforwardly gives the desired result. Note thatthe primes may be dropped from the double integral Amn once it has been isolatedfrom the rest of the expression.

    2. The potential on the surface of a sphere of radius a is specied by

    V 0

    0

    0

    V (, ) =0, 0 < V 0 ,

    0, <

    There are no other charges in this problem.

    a) Show that the potential outside the sphere may be expressed as

    (r,, ) =l =0 ,2 ,4 ,6 ,...

    V 0 [P l +1 (cos ) P l 1 (cos )]ar

    l +1P l (cos )

    where we take P 1 (x) = 0. Note that Legendre polynomials satisfy the relation(2l + 1) P l (x) = P l +1 (x) P l 1 (x).There are several ways of solving this problem. Perhaps the most straightforwardis to realize from azimuthal symmetry that the potential necessarily admits aLegendre expansion

    ( x ) =l

    lar

    l+1P l (cos )

  • 8/8/2019 F07 Practice Midterm Solutions

    5/8

  • 8/8/2019 F07 Practice Midterm Solutions

    6/8

    The quadrupole is given by l = 2. Hence the quadrupole moment is related tothe 2 term in the expansion. Since we do not care about normalization (we onlycare to maximize the moment) it is sufficient to write

    q2 ,0 2 P 3 (cos ) P 1 (cos )We extremize this quantity by taking a derivative with respect to and settingthe result to zero. Noting that taking a derivative simply undoes the integration,we end up solving

    dq2 ,0d

    P 2 (cos )sin = 12 (3cos2 1)sin = 0

    This is solved by = 0 and = cos 1 (1/ 3). Clearly the quadrupole momentvanishes when = 0 (as the entire sphere is at a constant zero potential). Hencethe maximum is when = cos 1 (1/ 3), provided we restrict 0 / 2.Note that > / 2 does not really make sense, except in a formal manner whereV 0 V 0 (corresponding to interchanging the integration limits). Technically,the problem should have asked to maximize the magnitude of the quadrupolemoment instead of to make q2 ,0 as positive as possible (which would depend onthe sign of V 0 ).

    3. A line charge on the z axis extends from z = a to z = + a and has linear chargedensity varying as

    x

    0

    0

    | z |

    z

    a

    a

    y (z) = 0 z , 0 < z a

    0 |z| , a z < 0

    where is a positive constant. The total charge on the 0 < z a segment is Q (andthe charge on the a z < 0 segment is Q).a) Calculate all of the multipole moments of the charge distribution. Make sure to

    indicate which moments are non-vanishing.

    Noting that a uniformly charged line charge on the + z axis has charge density

    = 0

    2r 2(cos 1)

    we see that a varying line charge yields

    =(r )2r 2

    (cos 1)

  • 8/8/2019 F07 Practice Midterm Solutions

    7/8

    This may be checked by observing

    dq = d3 x = r 2 drdd (cos ) = (r )drd2 cos =1

    (Note that there is no distinction between r and z for cos = 1.) Hence, for thepositive and negative line charge, we have

    = 0 r

    2r 2[(cos 1) (cos + 1)]

    We may normalize 0 by integrating from 0 to a to obtain the total charge Q

    Q = a

    0(z) dz =

    a

    0 0 z dz =

    0 a +1

    + 1

    Hence

    = Q + 12ar 2

    ra

    [(cos 1) (cos + 1)]

    The multipole moments are then

    qlm = r l Y lm () r 2 dr d = Q

    + 1a

    a

    0r l

    ra

    dr Y lm ()[(cos 1) (cos + 1)] d2

    By azimuthal symmetry, only the m = 0 moments are non-vanishing

    ql, 0 = Q + 1

    aa l +1

    + l + 1 1

    1 2l + 14 P l (cos )[(cos 1) (cos + 1)] d(cos )= Qa l

    + 1 + l + 1 2l + 14 [P l (1) P l (1)]

    = Qa l + 1

    + l + 1 2l + 14 [1()l ]Hence all moments vanish unless l is odd and m is zero. Then

    ql, 0 = Qa l + 1

    + l + 1 2l + 1 l oddb) Write down the multipole expansion for the potential in explicit form up to the

    rst two non-vanishing terms.

  • 8/8/2019 F07 Practice Midterm Solutions

    8/8

    The multipole expansion gives

    =1

    4 0l,m

    42l + 1

    qlmY lm (, )

    r l+1

    = 14 0l odd 42l + 1 ql, 0 P

    l (cos )r l+1

    =2Q

    4 0 rl odd

    + 1 + l + 1

    ar

    lP l (cos )

    =2Q

    4 0( + 1) a

    + 2P l (cos )

    r 2+

    ( + 1) a3

    + 4P 3 (cos )

    r 4+

    =2Q

    4 0( + 1) a

    + 2cos

    r 2+

    ( + 1) a3

    2( + 4)5cos3 3cos

    r 4+

    (4)

    c) What is the dipole moment p in terms of Q, a and ?

    Note that the dipole term in (4) has the form

    =1

    4 02Q( + 1) a

    + 2zr 3

    Comparing this with the dipole expression

    =1

    4 0 p xr 3

    gives

    p =2Q( + 1) a

    + 2z

    Alternatively, one could compute directly

    p = x d3 x = a

    axQ( + 1) |z|

    a +1sgn(z) dz

    x = y =0

    The z component is the only non-vanishing component

    pz = 2 Q( + 1) a

    0|z| +1a +1

    dz =2Q( + 1) a

    + 2