Page 101 - 6099
P. 101

The boundary and initial conditions become
                                   v(0, t) = 0,  v x (h, t) = 0,  v(x, 0) = f(x) − µ(x, 0).

                   Thus we have a heat equation with the source s(x, t)−µ t (x, t). We could apply separation
                   of variables to find a solution of the form
                                                        ∞       r
                                                       X           2      (2n − 1)πx
                                    u(x, t) = µ(x, t) +    u n (t)   sin                .
                                                                   h          2h
                                                       n=1

                   9.5.2     Direct eigenfunction expansion.

                   Alternatively we could seek a direct eigenfunction expansion of u(x, t).

                                                   ∞       r
                                                  X           2      (2n − 1)πx
                                         u(x, t) =    u n (t)   sin                .
                                                              h           2h
                                                   n=1
                   Note that the eigenfunctions satisfy the homogeneous boundary conditions while u(x, t)
                   does not. If we choose any fixed time t = t 0 and form the periodic extension of the
                   function u(x, t 0 ) to define it for x outside the range (0, h), then this function will have
                   jump discontinuities. This means that our eigenfunction expansion will not converge
                   uniformly. We are not allowed to differentiate the series with respect to x. We can’t just
                   plug the series into the partial differential equation to determine the coefficients. Instead,
                   we will multiply Equation 9.11, by an eigenfunction and integrate from x = 0 to x = h.
                   To avoid differentiating the series with respect to x, we will use integration by parts to
                                                                                                    2
                   move derivatives from u(x, t) to the eigenfunction. (We will denote λ n =  (2n−1)π  .)
                                                                                               2h
                             r    Z  h                            r    Z  h
                                2         p                          2         p
                                      sin( λ n x)(u t − κu xx ) dx =       sin( λ n x)s(x, t) dx
                                h  0                                 h  0
                                 r                         r          Z
                                    2  h       p      i h     2  p       h       p
                           0
                         u (t) −      κ u x sin( λ n x)  +      κ λ n     u x cos( λ n x) dx = s n (t)
                           n
                                    h                  0      h        0
                                         r                     r
                                           2                     2  p    h      p      i h
                                  0                n
                                 u (t) −     κ(−1) u x (h, t) +    κ λ n u cos( λ n x)    +
                                  n
                                           h                     h                       0
                                             r        Z  h
                                                2              p
                                           +     κλ n    u sin( λ n x) dx = s n (t)
                                                h      0
                                       r                 r
                                         2                  2  p
                                0                n
                               u (t) −     κ(−1) β(t) −      κ λ n u(0, t) + κλ n u n (t) = s n (t)
                                n
                                         h                  h
                                                     r
                                                        2   p
                                  0
                                                                             n
                                 u (t) + κλ n u n (t) =  κ     λ n α(t) + (−1) β(t) + s n (t)
                                  n
                                                        h
                   Now we have an ordinary differential equation for each of the u n (t). We obtain initial
                   conditions for them using the initial condition for u(x, t).
                                                    ∞       r
                                                   X           2     p
                                         u(x, 0) =     u n (0)   sin( λ n x) = f(x)
                                                               h
                                                   n=1
                                                             94
   96   97   98   99   100   101   102   103   104   105   106