Page 73 - 6099
P. 73

0
                   the origin is zero: φ (0) = 0. If φ(x) is defined only on the half-line, its even extension is
                   defined to be
                                                        (
                                                          φ(x),       for x ≥ 0,
                                             φ even (x) =                                          (7.17)
                                                          +φ(−x),     for x ≥ 0.

                   By the same reasoning as we used above, we end up with an explicit formula for w(x, t).
                   It is                            Z
                                               1      ∞         2             2
                                  w(x, t) = √           (e −(x−y) /4kt  + e −(x+y) /4kt )φ(y)dy.   (7.18)
                                              4πkt   0
                   Notice that the only difference between (7.15) and (7.18) is a single minus sign!


                   Example 7.4 Solve (7.16) with φ(x) ≡ 1. This is the same as Example 7.3 except
                   for the single sign. So we can copy from that example:


                                              1   1                 1    1
                                  u(x, t) =    + Erf( √   x  ) +      − Erf( √   x  )  = 1
                                              2   2       4kt       2    2       4kt

                   (That was stupid: We could have guessed it!)


                   7.3      Diffusion with a source


                   In this section we solve the inhomogeneous diffusion equation on the whole line

                                     u t − ku xx = f(x, t) (−∞ < x < ∞, 0 < t < ∞),                (7.19)

                                                       u(x, 0) = φ(x)

                   with f(x, t) and φ(x) arbitrary given functions. For instance, if u(x, t) represents the
                   temperature of a rod, then φ(x) is the initial temperature distribution and f(x, t) is a
                   source (or sink) of heat provided to the rod at later times. We will show that the solution
                   of (7.19) is

                                    ∞                        t   ∞
                                  Z                        Z Z
                        u(x, t) =      S(x − y, t)φ(y)dy +         S(x − y, t − s)f(y, s)dyds.     (7.20)
                                   −∞                       0   −∞

                   Notice that there is the usual term involving the initial data φ and another term involving
                   the source f. Both terms involve the source function S. Let’s begin by explaining
                   where (7.20) comes from. Later we will actually prove the validity of the formula. Our
                   explanation is an analogy. The simplest analogy is the ODE

                                               du
                                                   + Au(t) = f(t), u(0) = φ,                       (7.21)
                                                dt

                                                                       tA
                   where A is a constant. Using the integrating factor e , the solution is
                                                             Z  t
                                              u(t) = e −tA φ +   e (s−t)A f(s)ds.                  (7.22)
                                                               0



                                                             66
   68   69   70   71   72   73   74   75   76   77   78