Page 43 - 6637
P. 43

Solved Problems


                     and C characteristics through (x , t ) and let them intersect the x-axis at the
                                                            0
                            −
                                                                0
                     points A,B. Use the properties of these curves derived in part (a) to determine
                     u (x , t ) in terms of initial data at points A and B. Using a similar technique
                       t
                          0
                              0
                     to obtain u (x , τ) with 0 < τ < t, determine u(x , t ) by integration with
                                   t
                                      0
                                                                                   0
                                                                                      0
                     respect to τ, and compare this with the solution derived in class:
                                                                                    Z  x+ct
                                         1                                       1
                              u(x, t) =    (u(x + ct, 0) + u(x − ct, 0)) +                 u (τ, 0)dτ.
                                                                                             t
                                         2                                      2c   x−ct
                     Observe that this “method of characteristics” again shows that u(x , t ) de-
                                                                                                      0
                                                                                                         0
                     pends only on that part of the initial data between points A and B.
               Solution:
                 (a) We write the value of u along the line x−ct = k as a function of t: u(k+ct, t).
                     We differentiate u − cu with respect to t to see how the quantity varies.
                                          t
                                                  x
                               d                                                         2
                                  (u (k + ct, t) − cu (k + ct, t)) = cu + u − c u           xx  − cu xt
                                                                            xt
                                                       x
                                                                                   tt
                                     t
                               dt
                                                                                 2
                                                                      = u − c u     xx
                                                                           tt
                                                                      = 0
                     Thus u −cu is constant along the line x−ct = k. Now we examine u +cu                      x
                                    x
                              t
                                                                                                        t
                     along the line x + ct = k.
                              d
                                                                                          2
                                 (u (k − ct, t) + cu (k − ct, t)) = −cu + u − c u            xx  + cu xt
                                                                                    tt
                                   t
                                                                             xt
                                                      x
                             dt
                                                                                2
                                                                     = u − c u     xx
                                                                          tt
                                                                     = 0
                     u + cu is constant along the line x + ct = k.
                       t
                              x
                (b) From part (a) we know
                                u (x , t ) − cu (x , t ) = u (x − ct , 0) − cu (x − ct , 0)
                                                                                     x
                                                     0
                                                                                         0
                                                                   0
                                                                          0
                                                                t
                                                        0
                                                                                                0
                                                 x
                                     0
                                  t
                                         0
                                u (x , t ) + cu (x , t ) = u (x + ct , 0) + cu (x + ct , 0).
                                                                                     x
                                        0
                                                                t
                                                                          0
                                     0
                                 t
                                                                   0
                                                                                         0
                                                 x
                                                                                                0
                                                     0
                                                        0
                     We add these equations to find u (x , t ).
                                                                0
                                                            t
                                                                   0
                                   1
                     u (x , t )= (u (x − ct , 0)−cu (x −ct , 0)u (x + ct , 0)+cu (x +ct , 0))
                       t
                           0
                                                                                      0
                                                                                                          0
                                                                           t
                                                                                                x
                                                                               0
                                                                                                    0
                              0
                                       t
                                          0
                                                 0
                                                               0
                                                                     0
                                                           x
                                   2
                     Since t was arbitrary, we have
                             0
                                   1
                      u (x , τ)= (u (x − cτ, 0)−cu (x − cτ, 0)u (x + cτ, 0)+cu (x + cτ, 0))
                                                           x
                       t
                           0
                                                                                               x
                                                                                                   0
                                       t
                                                                            t
                                           0
                                                               0
                                                                               0
                                   2
                     for 0 < τ < t . We integrate with respect to τ to determine u(x , t ).
                                                                                                 0
                                                                                                    0
                                     0
                                                Z
                                                   t 0  1
                      u(x , t ) = u(x , 0) +           (u (x − cτ, 0) − cu (x − cτ, 0)u (x + cτ, 0)+
                                                                                                  0
                             0
                                                                              x
                                        0
                                                                                  0
                                                          t
                                                                                               t
                                                             0
                          0
                                                 0   2
                                                              39
   38   39   40   41   42   43   44   45   46   47   48