Page 70 - 6099
P. 70

√                 √
                   by substituting q = x/ 4kt, dq = (dx)/ 4kt. Now look more carefully at the sketch of
                   (x, t) for a very small t. If we cut out the tall spike, the rest of (x, t) is very small. Thus

                                                      max S(x, t) → 0
                                                      |x|>δ

                       Notice that the value of the solution u(x, t) is a kind of weighted average of the initial
                   values around the point x. Indeed, letting z = x − y in the integral (7.6), we can write
                                          Z  ∞
                                                                    X
                                 u(x, t) =     S(x − y, t)φ(y)dy =      S(x − y i , t)φ(y i )∆y i
                                           −∞                        i
                   approximately. This is the average of the solutions S(x − y i , t) with the weights φ(y i ).
                   For very small t, the source function is a spike so that the formula exaggerates the values
                   of φ near x. For any t > 0 the solution is a spread-out version of the initial values at
                   t = 0.
                       Here’s the physical interpretation. Consider diffusion. S(x−y, t) represents the result
                   of a unit mass (say, 1 gram) of substance located at time zero exactly at the position
                   y which is diffusing (spreading out) as time advances. For any initial distribution of
                   concentration, the amount of substance initially in the interval ∆y spreads out in time
                   and contributes approximately the term S(x − y i , t)φ(y i )∆y i . All these contributions are
                   added up to get the whole distribution of matter. Now consider heat flow. S(x − y, t)
                   represents the result of a ”hot spot” at y at time 0. The hot spot is cooling off and
                   spreading its heat along the rod.
                       Another physical interpretation is brownian motion, where particles move randomly in
                   space. For simplicity, we assume that the motion is one-dimensional; that is, the particles
                   move along a tube. Then the probability that a particle which begins at position x ends
                                                                R  b
                   up in the interval (a, b) at time t is precisely  S(x − y, t)dy for some constant k. In
                                                                 a
                   other words, if we let u(x, t) be the probability density (probability per unit length) and
                   if the initial probability density is φ(x), then the probability at all later times is given by
                   formula (7.6). That is, u(x, t) satisfies the diffusion equation.
                       It is usually impossible to evaluate the integral (7.8) completely in terms of elementary
                   functions. Answers to particular problems, that is, to particular initial data φ(x), are
                   sometimes expressible in terms of the error function of statistics,
                                                                Z  x
                                                             2       −p 2
                                                 Erf(x) = √         e   dp.                          (7.9)
                                                              π  0
                   Notice that Erf(0) = 0 and lim Erf(x) = 1.
                                               x→+∞

                   Example 7.1 From (7.5) we can write Q(x, t) in terms of Erf as

                                                           1    1
                                                 Q(x, t) =   + Erf( √   x  ).
                                                           2    2       4kt
                                                                                                      −x
                   Example 7.2 Solve the diffusion equation with the initial condition u(x, O) = e .
                   To do so, we simply plug this into the general formula (7.8):
                                                            Z  ∞
                                                       1         −(x−y) /4kt −y
                                                                       2
                                          u(x, t) = √           e          e dy.
                                                      4πkt   −∞

                                                             63
   65   66   67   68   69   70   71   72   73   74   75