Page 158 - Essentials of applied mathematics for scientists and engineers
P. 158

book   Mobk070    March 22, 2007  11:7








                     148  ESSENTIALS OF APPLIED MATHEMATICS FOR SCIENTISTS AND ENGINEERS
                       with boundary conditions


                                                            r (λ n , 0) = 0
                                                           (λ n , 1) = 0


                            The solution is the Bessel function J 0 (λ n r)with λ n given by J 0 (λ n ) = 0. Thus the
                       transform of u(t, r)isasfollows:
                                                                    1


                                           H λ u(t,r) = U(t,λ n ) =   rJ 0 (λ n r)u(t,r)dr
                                                                  r=0
                       This is called a Hankel transform. The appropriate differential equation for U(t,λ n )is

                                                          dU n     2
                                                              =−λ U n
                                                                   n
                                                          dt
                       so that

                                                                       2
                                                        U n (t,λ n ) = Be  −λ t
                                                                       n
                       Applying the initial condition, we find

                                                            1

                                                      B =    rf (r)J 0 (λ n r)dr
                                                          r=0
                       and from Eq. (9.19)


                                                     ∞     1
                                                            rf (r)J 0 (λ n r)dr    2
                                                         r=0                     −λ t
                                           u(t,r) =                       J 0 (λ n r)e  n
                                                            #      # 2
                                                             J 0 (λ n r) #
                                                    n=0
                                                            #
                       Example 9.4 (The sine transform with a source). Next consider a one-dimensional transient
                       diffusion with a source term q(x):
                                                   u t = u xx + q(x)

                                                   y(0, x) = y(t, 0) = t(t,π) = 0

                       First we determine that the sine transform is appropriate. The operator   is such that


                                                           =   xx = λ
   153   154   155   156   157   158   159   160   161   162   163