MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20

MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20

MA 201 (2016), PDE 2 / 20

Vibrating string and the wave equation Consider a stretched string of length π with the ends fastened to the ends x = 0 and x = π. Suppose that the string is set to vibrate by displacing it from its equilibrium position and then releasing it. Let u(x,t) denote the transverse displacement at time t 0 of the point on the string at position x. In particular, u(x,0) denotes the initial shape of the string and u t (x,0) denotes the initial velocity. MA 201 (2016), PDE 3 / 20

Vibrating string and the wave equation (Contd.) Under some physical assumptions, which are not far from being realistic: We arrive at an equation known as the one-dimensional wave equation, which governs the entire process. The equation under consideration is u tt = c 2 u xx, (1) which is a linear second-order homogeneous equation. Here c represents the velocity of the wave. Note: the above equation is derived without considering any external force acting on the string. MA 201 (2016), PDE 4 / 20

Vibrating string and the wave equation (Contd.) If additional transverse forces are acting on the string (due to its weight or other impressed exterior forces), let F(x,t) denote the amount of force per unit length acting in the u-direction. Then the equation of forced vibrations of the string is u tt = c 2 u xx +F(x,t). (2) Two special cases of (2) particularly generate interest: When the external force is due to the gravitational acceleration g only (consider a string oriented horizontally), the equation becomes u tt = c 2 u xx g. (3) When the external force is due to the resistance which is proportional to the instantaneous velocity (a string vibrating in a fluid), the equation becomes u tt +2ku t = c 2 u xx, k is a positive constant. (4) MA 201 (2016), PDE 5 / 20

The D Alembert s solution of the wave equation The PDE describing the vibration of a string (one-dimensional wave equation) is probably the most easily recognized hyperbolic equation. Method of characteristics is very useful for hyperbolic equations. Two families of characteristics of hyperbolic equations, being real and distinct, are of considerable practical value. In one-dimensional progressive wave propagation, consideration of characteristics can give us a good deal of information about the propagation of the wave fronts. This solution of one-dimensional wave equation, known as D Alembert s solution, was discovered by a French mathematician named Jean Le Rond D Alembert. MA 201 (2016), PDE 6 / 20

The D Alembert solution of the wave equation (Contd.) Consider the one-dimensional wave equation: Here B 2 4AC = 4c 2 > 0. u tt = c 2 u xx. (5) The characteristics are given by Under this transformation ξ = x ct, η = x+ct. u x = u ξ ξ x +u η η x = u ξ +u η, u t = u ξ ξ t +u η η t = cu ξ +cu η, u xx = (u ξξ ξ x +u ξη η x )+(u ηξ ξ x +u ηη η x ) = u ξξ +2u ξη +u ηη, u tt = c(u ξξ ξ t +u ξη η t )+c(u ηξ ξ t +u ηη η t ), = c 2 (u ξξ 2u ξη +u ηη ). MA 201 (2016), PDE 7 / 20

The D Alembert solution of the wave equation (Contd.) Substituting these into the given equation u ξη = 0. (6) Integrating partially with respect to η: u ξ = F (ξ). Integrating partially w.r.t. ξ u = F(ξ)+G(η). The solution in physical variables: u(x,t) = F(x ct)+g(x+ct) (7) where F and G are arbitrary functions. MA 201 (2016), PDE 8 / 20

The D Alembert solution of the wave equation (Contd.) The physical interpretation of these functions is quite interesting. The functions F and G represent two progressive waves travelling in opposite directions with the speed c. To see this let us first consider the solution u = F(x ct). At t = 0, it defines the curve u = F(x), and after time t = t 1, it defines the curve u = F(x ct 1 ). But these curves are identical except that the latter is translated to the right a distance equal to ct 1. MA 201 (2016), PDE 9 / 20

Method of Characteristics (Contd.) Figure : A Progressive Wave MA 201 (2016), PDE 10 / 20

The D Alembert solution of the wave equation (Contd.) Thus the entire configuration moves along the positive direction of x-axis a distance of ct 1 in time t 1. The velocity with which the wave is propagated is, therefore, v = ct 1 t 1 = c Similarly the function G(x + ct) defines a wave progressing in the negative direction of x-axis with constant velocity c. The total solution is, therefore, the algebraic sum of these two travelling waves. Solution (7) is a very convenient representation for progressive waves which travel large distances through a uniform medium. MA 201 (2016), PDE 11 / 20

The D Alembert solution of the wave equation (Contd.) It is to be noted that the use of string problem to demonstrate the solution of the wave problem is a matter of convenience. However, any variables satisfying the wave equation possess the mathematical properties developed for the string. For instance, u(x,t) = sin(x±ct) represent sinusoidal waves traveling with speed c in the positive and negative directions respectively without change of shape. MA 201 (2016), PDE 12 / 20

The D Alembert solution of the wave equation (Contd.) MA 201 (2016), PDE 13 / 20

The D Alembert solution of the wave equation (Contd.) Let us consider the following two initial conditions for a uniform medium < x <. Displacement: u(x, 0) = φ(x), (8) Velocity: u t (x,0) = ψ(x), (9) that is we consider the vibration of a thin string of infinite length with some initial displacement and initial velocity. From the solution (7) we find that F(x)+G(x) = φ(x), (10) cf (x)+c G (x) = ψ(x), (11) for all values of x. MA 201 (2016), PDE 14 / 20

The D Alembert solution of the wave equation (Contd.) Integrating the second equation with respect to x F(x)+G(x) = 1 c x where A is an integration constant and τ is a dummy variable. From (10) and (12) F(x) = 1 2 G(x) = 1 2 Substituting these expressions into (7) [ φ(x) 1 x c [ φ(x)+ 1 x c x 0 ψ(τ)dτ +A. (12) x 0 ψ(τ)dτ x 0 ψ(τ)dτ u(x,t) = 1 2 [φ(x ct)+φ(x+ct)]+ 1 2c ] A/2, (13) ] +A/2, (14) x+ct x ct This is D Alembert s solution of one-dimensional wave equation. ψ(τ)dτ. (15) MA 201 (2016), PDE 15 / 20

The D Alembert solution of the wave equation (Contd.) Thus for a given initial displacement and velocity in the vertical direction, the wave equation is completely solved and this solution is usually called the progressive wave solution. It is easy to verify by direct substitution that u(x,t), represented by (15), is the unique solution of the wave equation (5) provided φ is twice continuously differentiable and ψ is continuously differentiable. It is important to note that the solution u(x,t) depends only on the initial values of φ at points x ct and x+ct and values of ψ between these two points. In other words, the solution does not depend at all on initial values outside this interval, x ct x x+ct. This interval is called the domain of dependence of the variables (x,t). MA 201 (2016), PDE 16 / 20

Special cases of D Alembert solution CASE I (Initial velocity zero). Suppose the string has IC u(x,0) = φ(x), < x <, u t (x,0) = 0, < x <. Then D Alembert s solution is u(x,t) = 1 2 [φ(x ct)+φ(x+ct)]. MA 201 (2016), PDE 17 / 20

Special cases of D Alembert solution (Contd.) CASE 2. (Initial displacement zero) Suppose the string has the following IC: u(x,0) = 0, < x < u t (x,0) = ψ(x), < x <. In this case, the solution is u(x,t) = 1 2c x+ct x ct ψ(τ)dτ. The solution u at (x, t) may be interpreted as integrating the initial velocity between x ct and x+ct on the initial line t = 0. MA 201 (2016), PDE 18 / 20

Special cases of D Alembert solution (Contd.) Example (Zero initial velocity) Solve the IVP: PDE: u tt = c 2 u xx, < x <, t > 0 IC: u(x, 0) = sin(x), u t (x,0) = 0. Solution: Using D Alembert s formula with φ(x) = sin(x) and ψ(x) = 0, we obtain u(x,t) = 1 2 [sin(x ct)+sin(x+ct)]. MA 201 (2016), PDE 19 / 20

Special cases of D Alembert solution (Contd.) Example (Zero initial displacement) Consider the IVP: PDE: u tt = c 2 u xx, < x <, t > 0 I.C. u(x,0) = 0, u t (x,0) = sin(x). Solution: Here the string is initially straight (u(x,0) = 0), but has a variable velocity at t = 0 (u t (x,0) = sin(x)). Thus, using D Alembert s formula with φ(x) = 0 and ψ(x) = sin(x), we obtain u(x,t) = 1 2c x+ct x ct sin(τ)dτ = 1 2c [cos(x+ct) cos(x ct)]. MA 201 (2016), PDE 20 / 20