Strauss PDEs 2e: Section 3.3 - Exercise 2 Page of 6 Exercise 2 Solve the completely inhomogeneous diffusion problem on the half-line v t kv xx = f(x, t) for < x <, < t < v(, t) = h(t) v(x, ) = φ(x), by carrying out the subtraction method begun in the text. Solution Following the subtraction method in the text, make the substitution, V (x, t) = v(x, t) h(t), in order to make the Dirichlet boundary condition homogeneous. This will allow us to use the method of reflection to solve the problem as in the previous exercise. Find the derivatives of v in terms of this new variable. V t = v t h (t) v t = V t h (t) V x = v x V xx = v xx The PDE V satisfies is then [V t h (t) kv xx = f(x, t), or V t kv xx = f(x, t) h (t), < x <, t >. The initial and boundary conditions associated with this PDE are V (x, ) = v(x, ) h() = φ(x) h() V (, t) = v(, t) h(t) = h(t) h(t) =. Solution by the Method of Reflection The strategy here is to solve this same problem over the whole line ( < x < ) using the odd extensions of the functions, f(x, t) h (t) and φ(x) h(), and then to restrict the domain to x > to get the solution to V. Once we have V (x, t), we can add h(t) to it to get v(x, t), the desired solution. The point of using the odd extensions of the functions is so that the Dirichlet boundary condition V (, t) = is satisfied. Let W (x, t) be the solution to the same problem over the whole line, W t kw xx = F odd (x, t), < x <, t > W (x, ) = Φ odd (x), where F odd (x, t) = f(x, t) h (t) x > f( x, t) h (t) x < and Φ odd (x) = φ(x) h() x > φ( x) h() x <.
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page 2 of 6 Use the fact that the diffusion equation is linear to split the problem over the whole line into two. Let W (x, t) = w(x, t) u(x, t), where w and u satisfy the following problems. w t kw xx =, < x <, t > u t ku xx = F odd (x, t), < x <, t > w(x, ) = Φ odd (x) u(x, ) = The solution for w is given in section 2.4 on page 49. w(x, t) = Φ odd (r) dr According to Duhamel s principle, the solution to the inhomogeneous diffusion equation is u(x, t) = t U(x, t s; s) ds, where U = U(x, t; s) is the solution to the associated homogeneous equation with a particular choice for the initial condition. U t ku xx =, < x <, t > U(x, ; s) = F odd (x, s) Use the solution on page 49 again to solve for U. U(x, t; s) = F odd (r, s) dr The solution to the inhomogeneous diffusion equation is then Consequently, W (x, t) = u(x, t) = = t t U(x, t s; s) ds Φ odd (r) dr t F odd (r, s) dr ds. The solution to V is then just the restriction of W to x >. V (x, t) = t Φ odd (r) dr F odd (r, s) dr ds. F odd (r, s) dr ds, x >
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page 3 of 6 Our task now is to write the solution in terms of the given functions, φ(x) and f(x, t). Split up each of the integrals into one over the negative values of x and one over the positive values of x and substitute the appropriate functions for Φ odd and F odd in these intervals. V (x, t) = t φ( r) h() dr f( r, s) h (s) dr [φ(r) h() dr [f(r, s) h (s) dr ds, x > Substitute q = r in the integrals from to and substitute q = r in the integrals from to. V (x, t) = t (x q)2 φ(q) h()( dq) (x q)2 f(q, s) h (s)( dq) Use the minus signs in front of dq to make the integrals go from to. V (x, t) = t (x q)2 φ(q) h() dq Combine the integrals like so. V (x, t) = t (x q)2 [φ(q) h() dq [f(q, s) h (s) dq ds, x > f(q, s) h (s) dq [φ(q) h() dq [f(q, s) h (s) dq ds, x > (x q)2 (x q)2 h() dq (x q)2 f(q, s) dq (x q)2 h (s) dq ds, x >
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page 4 of 6 Distribute the integral in ds. V (x, t) = t t (x q)2 (x q)2 h() dq Switch the order of integration in the last integral. (x q)2 V (x, t) = (x q)2 t t (x q)2 h() dq (x q)2 (x q)2 h (s) dq ds, x > (x q)2 Integrate the last integral by parts in order to remove the derivative from h. (x q)2 V (x, t) = (x q)2 h() dq t (x q)2 (x q)2 h(s) t s h (s) ds dq, x > t dq (x q)2 h(s) ds dq, x >
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page 5 of 6 The terms with h() cancel. V (x, t) = (x q)2 (x q)2 h() dq t t (x q)2 (x q)2 h(s) dq s=t (x q)2 h() dq t s (x q)2 h(s) ds dq, The integrand evaluated at s = t is zero because the onents of e are. What remains then is (x q)2 V (x, t) = = t t t s s Therefore, since v(x, t) = V (x, t) h(t), v(x, t) = h(t) t t s (x q)2 (x q)2 (x q)2 (x q)2 (x q)2 (x q)2 (x q)2 (x q)2 x > h(s) ds dq, x > h(s) dq ds, x >. h(s) dq ds, x >.
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page 6 of 6 The solution can be written compactly as v(x, t) = h(t) t [G(x q, t) G(x q, t) [G(x q, t s) G(x q, t s) t [G(x q, t s) G(x q, t s)h(s) dq ds, x >, s where G = G(x, t) is the Green s function for the one-dimensional diffusion equation. G(x, t) = ) ( x2