Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
|
|
- Joan Morgan
- 6 years ago
- Views:
Transcription
1 Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44
2 1 The Vibrating String Equation 2 Second order PDEs 3 The D Alembert solution 4 The Klein-Gordon and the telegrapher s equations 5 Solutions on the real line (Fourier Transforms) 6 Finite domains 7 Non homogeneous problems and resonances Lucio Demeio - DIISM wave equation 2 / 44
3 The vibrating string equation Lucio Demeio - DIISM wave equation 3 / 44
4 The vibrating string equation m= ρ x T (x+ x) β α T ( x) u( x+ x,t) u( x,t) x x+ x Lucio Demeio - DIISM wave equation 4 / 44
5 The vibrating string equation m= ρ x T (x+ x) β α T ( x) u( x+ x,t) u( x,t) x x+ x Newton s law for m along x and y gives 0 = T (x + x) cos β T (x) cos α No horizontal displacement T (x + x) cos β = T (x) cos α T 0 m 2 u = T (x + x) sin β T (x) sin α (1) t2 Lucio Demeio - DIISM wave equation 4 / 44
6 The vibrating string equation With m = ρ x and by substituting from (1): 2 u t 2 = T 0 tan β tan α ρ x = v 2 u x u (x + x, t) x (x, t) v 2 2 u x x 2 2 u t 2 v2 2 u x 2 = 0 (2) which is called wave equation (in 1 D). Lucio Demeio - DIISM wave equation 5 / 44
7 The vibrating string equation With m = ρ x and by substituting from (1): 2 u t 2 = T 0 tan β tan α ρ x = v 2 u x u (x + x, t) x (x, t) v 2 2 u x x 2 2 u t 2 v2 2 u x 2 = 0 (2) which is called wave equation (in 1 D). Generalizing to 2 or 3 dimensions: 2 u t 2 v2 u = 0 The wave equation describes virtually all wave phenomena: sound waves, light waves, waves in fluids, gases and plasmas, waves in solids, etc. Lucio Demeio - DIISM wave equation 5 / 44
8 Quasilinear second order PDEs Lucio Demeio - DIISM wave equation 6 / 44
9 Quasilinear second order PDEs General form in 2 variables The wave equation is a particular case of a second-order PDE; the most general form is a 2 u x + 2 b 2 u 2 x y + c 2 u y = d (3) 2 where a = a(x, y, u, u/ x, u/ y), and so for b, c and d (quasilinearity). Lucio Demeio - DIISM wave equation 7 / 44
10 Quasilinear second order PDEs General form in 2 variables The wave equation is a particular case of a second-order PDE; the most general form is a 2 u x + 2 b 2 u 2 x y + c 2 u y = d (3) 2 where a = a(x, y, u, u/ x, u/ y), and so for b, c and d (quasilinearity). Boundary Value and Initial Value Problems A differential equation possesses a family of solutions, and auxiliary conditions need to be specified in order to assure uniqueness of the solution. There are two sorts of auxiliary conditions: Lucio Demeio - DIISM wave equation 7 / 44
11 Quasilinear second order PDEs General form in 2 variables The wave equation is a particular case of a second-order PDE; the most general form is a 2 u x + 2 b 2 u 2 x y + c 2 u y = d (3) 2 where a = a(x, y, u, u/ x, u/ y), and so for b, c and d (quasilinearity). Boundary Value and Initial Value Problems A differential equation possesses a family of solutions, and auxiliary conditions need to be specified in order to assure uniqueness of the solution. There are two sorts of auxiliary conditions: Initial Value (Cauchy) Problem: when the auxiliary conditions are imposed on the function and its derivatives along a line belonging to the domain of the independent variables; Lucio Demeio - DIISM wave equation 7 / 44
12 Quasilinear second order PDEs General form in 2 variables The wave equation is a particular case of a second-order PDE; the most general form is a 2 u x + 2 b 2 u 2 x y + c 2 u y = d (3) 2 where a = a(x, y, u, u/ x, u/ y), and so for b, c and d (quasilinearity). Boundary Value and Initial Value Problems A differential equation possesses a family of solutions, and auxiliary conditions need to be specified in order to assure uniqueness of the solution. There are two sorts of auxiliary conditions: Initial Value (Cauchy) Problem: when the auxiliary conditions are imposed on the function and its derivatives along a line belonging to the domain of the independent variables; Boundary Value Problem: the solution is sought for in a domain Ω R n and the auxiliary conditions are imposed on the function (or its derivatives) on the boundary Ω of the domain. Lucio Demeio - DIISM wave equation 7 / 44
13 Quasilinear second order PDEs Initial Value and Boundary Value Problems Cauchy Problem: y O Ω γ x u(x, y) = h(x, y) u (x, y) = φ(x, y) x u (x, y) = ψ(x, y) y On γ y Boundary Value Problem: Ω O Ω x u(x, y) = f(x, y) or u (x, y) = ψ(x, y) n On Ω Lucio Demeio - DIISM wave equation 8 / 44
14 Quasilinear second order PDEs Cauchy problem Existence and uniqueness of the solution starting from assigned values of u, u/ x and u/ y on a smooth curve in the (x, y) plane. Let { x = f(s) y = g(s) be a parametric representation of a curve γ in the (x, y) plane. Lucio Demeio - DIISM wave equation 9 / 44
15 Quasilinear second order PDEs Cauchy problem Existence and uniqueness of the solution starting from assigned values of u, u/ x and u/ y on a smooth curve in the (x, y) plane. Let { x = f(s) y = g(s) be a parametric representation of a curve γ in the (x, y) plane. Cauchy problem We must prescribe u(x(s), y(s)) = h(s) (4) u (x(s), y(s)) = φ(s) x (5) u (x(s), y(s)) = ψ(s) y (6) and find the conditions for existence and uniqueness of the solution. Lucio Demeio - DIISM wave equation 9 / 44
16 Quasilinear second order PDEs Let us differentiate (5) and (6) w.r. to s; by adding these to the differential equation (3) we have the following system: 2 u x 2 f (s) + 2 u x y g (s) = φ (s) (7) 2 u x y f (s) + 2 u y 2 g (s) = ψ (s) (8) a 2 u x b 2 u x y + c 2 u = d, (9) y2 where the partial derivatives are to be considered as unknowns. Lucio Demeio - DIISM wave equation 10 / 44
17 Quasilinear second order PDEs Let us differentiate (5) and (6) w.r. to s; by adding these to the differential equation (3) we have the following system: 2 u x 2 f (s) + 2 u x y g (s) = φ (s) (7) 2 u x y f (s) + 2 u y 2 g (s) = ψ (s) (8) a 2 u x b 2 u x y + c 2 u = d, (9) y2 where the partial derivatives are to be considered as unknowns. Unique solution if the determinant f (s) g (s) 0 = 0 f (s) g (s) a 2 b c = a [g (s)] 2 2 b f (s) g (s) + c [f (s)] 2 0 If = 0, γ is called characteristic curve. Lucio Demeio - DIISM wave equation 10 / 44
18 Classification Lucio Demeio - DIISM wave equation 11 / 44
19 Classification The condition a [g (s)] 2 2 b f (s) g (s) + c [f (s)] 2 = 0 is equivalent to a [ g ] 2 (s) f 2 b g (s) (s) f (s) + c = 0 Lucio Demeio - DIISM wave equation 12 / 44
20 Classification The condition a [g (s)] 2 2 b f (s) g (s) + c [f (s)] 2 = 0 is equivalent to a [ g ] 2 (s) f 2 b g (s) (s) f (s) + c = 0 If b 2 a c > 0 we have g (s) f (s) = b ± b2 a c a (10) and the equation is called hyperbolic. Lucio Demeio - DIISM wave equation 12 / 44
21 Classification The condition a [g (s)] 2 2 b f (s) g (s) + c [f (s)] 2 = 0 is equivalent to a [ g ] 2 (s) f 2 b g (s) (s) f (s) + c = 0 If b 2 a c > 0 we have g (s) f (s) = b ± b2 a c a (10) and the equation is called hyperbolic. If b 2 a c < 0 we have no solutions for g (s)/f (s) and the equation is called elliptic Lucio Demeio - DIISM wave equation 12 / 44
22 Classification The condition a [g (s)] 2 2 b f (s) g (s) + c [f (s)] 2 = 0 is equivalent to a [ g ] 2 (s) f 2 b g (s) (s) f (s) + c = 0 If b 2 a c > 0 we have g (s) f (s) = b ± b2 a c a (10) and the equation is called hyperbolic. If b 2 a c < 0 we have no solutions for g (s)/f (s) and the equation is called elliptic If b 2 a c = 0 we have and the equation is called parabolic. g (s) f (s) = b a Lucio Demeio - DIISM wave equation 12 / 44
23 Classification Characteristics Thus, hyperbolic equations possess two families of (real) characteristic curves in the (x, y) plane; Lucio Demeio - DIISM wave equation 13 / 44
24 Classification Characteristics Thus, hyperbolic equations possess two families of (real) characteristic curves in the (x, y) plane; parabolic equations have one family; Lucio Demeio - DIISM wave equation 13 / 44
25 Classification Characteristics Thus, hyperbolic equations possess two families of (real) characteristic curves in the (x, y) plane; parabolic equations have one family; elliptic equations have none (or have complex characteristics); Lucio Demeio - DIISM wave equation 13 / 44
26 Classification Characteristics Thus, hyperbolic equations possess two families of (real) characteristic curves in the (x, y) plane; parabolic equations have one family; elliptic equations have none (or have complex characteristics); If a, b and c depend only on x and y (linear equation), equation (10) can be written as dy dx = b ± b2 a c λ ±(x, y) (11) a and it becomes a set of two differential equations for the characteristic curves. Lucio Demeio - DIISM wave equation 13 / 44
27 Classification Characteristics Thus, hyperbolic equations possess two families of (real) characteristic curves in the (x, y) plane; parabolic equations have one family; elliptic equations have none (or have complex characteristics); If a, b and c depend only on x and y (linear equation), equation (10) can be written as dy dx = b ± b2 a c λ ±(x, y) (11) a and it becomes a set of two differential equations for the characteristic curves. Role of the characteristic curves In time dependent problems, one can say, qualitatively, that the characteristics are the lines along which solutions are transported ; Lucio Demeio - DIISM wave equation 13 / 44
28 Classification Characteristics Thus, hyperbolic equations possess two families of (real) characteristic curves in the (x, y) plane; parabolic equations have one family; elliptic equations have none (or have complex characteristics); If a, b and c depend only on x and y (linear equation), equation (10) can be written as dy dx = b ± b2 a c λ ±(x, y) (11) a and it becomes a set of two differential equations for the characteristic curves. Role of the characteristic curves In time dependent problems, one can say, qualitatively, that the characteristics are the lines along which solutions are transported ; The presence of two families of characteristics in hyperbolic equations corresponds to counter-propagating wave forms; Lucio Demeio - DIISM wave equation 13 / 44
29 Classification Role of the characteristic curves the presence of one family of characteristics in parabolic equations corresponds to relaxation towards a statistical equilibrium (e.g., in gases); Lucio Demeio - DIISM wave equation 14 / 44
30 Classification Role of the characteristic curves the presence of one family of characteristics in parabolic equations corresponds to relaxation towards a statistical equilibrium (e.g., in gases); elliptic equations, which do not possess characteristics, describe static problems (e.g., static electric fields due to a charge distribution, steady-state temperature distributions, linear elasticity problems and more). Lucio Demeio - DIISM wave equation 14 / 44
31 Classification Role of the characteristic curves the presence of one family of characteristics in parabolic equations corresponds to relaxation towards a statistical equilibrium (e.g., in gases); elliptic equations, which do not possess characteristics, describe static problems (e.g., static electric fields due to a charge distribution, steady-state temperature distributions, linear elasticity problems and more). Prototypical equations Wave equation: Heat equation: 2 u t 2 v2 2 u x 2 = 0 b 2 a c = v 2 > 0, hyperbolic equation u t κ 2 u x 2 = 0 b 2 a c = 0, parabolic equation Laplace equation: 2 u x u y 2 = 0 b 2 a c = 1 < 0, elliptic equation Lucio Demeio - DIISM wave equation 14 / 44
32 Classification Canonical form Lu = a 2 u x b 2 u x y + c 2 u y 2 Lucio Demeio - DIISM wave equation 15 / 44
33 Classification Canonical form Lu = a 2 u x b 2 u x y + c 2 u y 2 Polynomial form (a, b, c constants for simplicity) p L (λ) = aλ b λ + c = a(λ λ 1 ) (λ λ 2 ) and, by analogy, ( ) ( ) Lu = a x λ 1 y x λ 2 u y Lucio Demeio - DIISM wave equation 15 / 44
34 Classification Reduction to canonical form (after linear transformation to (ξ, η)): (λ 1 λ 2 R) : Lu 2 u ξ η or also Lu 2 u ξ 2 2 u η 2 (λ 1 = λ 2 R) : Lu 2 u η 2 (λ 1 = λ 2 C) : Lu 2 u ξ u η 2 hyperbolic case parabolic case elliptic case Lucio Demeio - DIISM wave equation 16 / 44
35 The D Alembert solution Lucio Demeio - DIISM wave equation 17 / 44
36 Back to the wave equation Characteristics From equation (11) with a = 1, b = 0, c = v 2 dy dx = ±v x + v t = ξ x v t = η Lucio Demeio - DIISM wave equation 18 / 44
37 Back to the wave equation Characteristics From equation (11) with a = 1, b = 0, c = v 2 dy dx = ±v x + v t = ξ x v t = η x+v t= ξ x-v t= η Lucio Demeio - DIISM wave equation 18 / 44
38 Back to the wave equation Characteristics From equation (11) with a = 1, b = 0, c = v 2 dy dx = ±v x + v t = ξ x v t = η x+v t= ξ x-v t= η Canonical Form The transformation (x, t) (ξ, η) leads to: 2 u t 2 v2 2 u x 2 = 0 = 2 U ξ η = 0 Lucio Demeio - DIISM wave equation 18 / 44
39 Back to the wave equation Characteristics From equation (11) with a = 1, b = 0, c = v 2 dy dx = ±v x + v t = ξ x v t = η x+v t= ξ x-v t= η Canonical Form The transformation (x, t) (ξ, η) leads to: 2 u t 2 v2 2 u x 2 = 0 = 2 U ξ η = 0 General Solution U(ξ, η) = F (ξ) + G(η) (12) u(x, t) = F (x + v t) + G(x v t) (13) Lucio Demeio - DIISM wave equation 18 / 44
40 D Alembert s form of the solution Initial value problem We prescribe the initial data u(x, 0) = h(x) (14) u (x, 0) = ψ(x) t (15) Lucio Demeio - DIISM wave equation 19 / 44
41 D Alembert s form of the solution Initial value problem We prescribe the initial data u(x, 0) = h(x) (14) u (x, 0) = ψ(x) t (15) Equations (12) and (13) give: F (x) + G(x) = h(x) (16) v [F (x) G (x)] = ψ(x) (17) Lucio Demeio - DIISM wave equation 19 / 44
42 D Alembert s form of the solution Initial value problem We prescribe the initial data u(x, 0) = h(x) (14) u (x, 0) = ψ(x) t (15) Equations (12) and (13) give: F (x) + G(x) = h(x) (16) v [F (x) G (x)] = ψ(x) (17) and, by differentiating (16): F (x) = v h (x) + ψ(x) 2 v G (x) = v h (x) ψ(x) 2 v (18) (19) Lucio Demeio - DIISM wave equation 19 / 44
43 D Alembert s form of the solution After integration F (x) = h(x) 2 G(x) = h(x) v 1 2 v x 0 x 0 ψ(λ)dλ (20) ψ(λ)dλ (21) Lucio Demeio - DIISM wave equation 20 / 44
44 D Alembert s form of the solution After integration F (x) = h(x) 2 G(x) = h(x) v 1 2 v x 0 x 0 ψ(λ)dλ (20) ψ(λ)dλ (21) We finally obtain the d Alembert solution of the Initial Value Problem: u(x, t) = F (x + v t) + G(x v t) = = 1 1 [h(x + v t) + h(x v t)] v x+v t x v t ψ(λ)dλ (22) Lucio Demeio - DIISM wave equation 20 / 44
45 D Alembert s form of the solution After integration F (x) = h(x) 2 G(x) = h(x) v 1 2 v x 0 x 0 ψ(λ)dλ (20) ψ(λ)dλ (21) We finally obtain the d Alembert solution of the Initial Value Problem: u(x, t) = F (x + v t) + G(x v t) = = 1 1 [h(x + v t) + h(x v t)] v x+v t x v t ψ(λ)dλ (22) Works well on infinite domains; on a finite domain, we must use it piecewise and take reflections into account. Lucio Demeio - DIISM wave equation 20 / 44
46 Dependence domain and Influence cone Domain of dependence From the d Alembert form (22), we see that u(x, t) depends upon the values of the functions h and ψ in the interval [x v t, x + v t] on the x axis. This is called domain of dependence. Lucio Demeio - DIISM wave equation 21 / 44
47 Dependence domain and Influence cone Domain of dependence From the d Alembert form (22), we see that u(x, t) depends upon the values of the functions h and ψ in the interval [x v t, x + v t] on the x axis. This is called domain of dependence. The meaning is that the solution at position x and time t depends upon the initial data within [x v t, x + v t] and not outside this interval. Lucio Demeio - DIISM wave equation 21 / 44
48 Dependence domain and Influence cone Domain of dependence From the d Alembert form (22), we see that u(x, t) depends upon the values of the functions h and ψ in the interval [x v t, x + v t] on the x axis. This is called domain of dependence. The meaning is that the solution at position x and time t depends upon the initial data within [x v t, x + v t] and not outside this interval. Cone of influence We also see that the initial condition at x = x 0 influences the solution at later times within the cone (a triangle in a plane) delimited by the surfaces x 0 v t = 0 and x 0 + v t = 0. This is called cone of influence. Lucio Demeio - DIISM wave equation 21 / 44
49 Dependence domain and Influence cone Domain of dependence From the d Alembert form (22), we see that u(x, t) depends upon the values of the functions h and ψ in the interval [x v t, x + v t] on the x axis. This is called domain of dependence. The meaning is that the solution at position x and time t depends upon the initial data within [x v t, x + v t] and not outside this interval. Cone of influence We also see that the initial condition at x = x 0 influences the solution at later times within the cone (a triangle in a plane) delimited by the surfaces x 0 v t = 0 and x 0 + v t = 0. This is called cone of influence. These concepts are consistent with the experience of a finite propagation speed. Lucio Demeio - DIISM wave equation 21 / 44
50 Dependence domain and Influence cone t ( x,t) x- vt x+vt Domain of dependence x Lucio Demeio - DIISM wave equation 22 / 44
51 Dependence domain and Influence cone t ( x,t) x - vt 0 0 x +vt x- vt x+vt Domain of dependence x 0 Cone of influence x Lucio Demeio - DIISM wave equation 22 / 44
52 The Klein-Gordon and the telegrapher s equations Lucio Demeio - DIISM wave equation 23 / 44
53 The Klein-Gordon and the telegrapher s equations The Klein-Gordon equation If the vibrating string is subjected to a uniformly distributed elastic force, we obtain the Klein-Gordon equation: 2 u t 2 v2 2 u + γu = 0 (23) x2 where γ is a constant (e.g., proportional to the elastic constant). u( x,t) x Lucio Demeio - DIISM wave equation 24 / 44
54 The Klein-Gordon and the telegrapher s equations The telegrapher s equation Propagation of waves rarely happens without dissipation; when taken into account, we obtain the telegrapher s equation (or transmission line equation): where λ is a damping constant. 2 u t 2 v2 2 u x λ u + γu = 0 (24) t u( x,t) u( x,t) x x Lucio Demeio - DIISM wave equation 25 / 44
55 Solutions on the real line (Fourier Transforms) Lucio Demeio - DIISM wave equation 26 / 44
56 Solutions by Fourier Transforms The wave equation Assume that u(x, t) can be written as a Fourier integral (which means that u(x, t) 0 fast enough as x ± ) u(x, t) = û(k, t) e i k x dk û(k, t) = 1 2π and substitute in the wave equation. One obtains 2 û t 2 + ω(k)2 û = 0 with ω(k) = k v. The general solution is given by û(k, t) = A(k) e iωt + B(k)e iωt u(x, t) e i k x dx Inversion of the Fourier transform gives back the d Alembert solution. Lucio Demeio - DIISM wave equation 27 / 44
57 Solutions by Fourier Transforms The wave equation Assume that u(x, t) can be written as a Fourier integral (which means that u(x, t) 0 fast enough as x ± ) u(x, t) = û(k, t) e i k x dk û(k, t) = 1 2π and substitute in the wave equation. One obtains 2 û t 2 + ω(k)2 û = 0 with ω(k) = k v. The general solution is given by û(k, t) = A(k) e iωt + B(k)e iωt u(x, t) e i k x dx Inversion of the Fourier transform gives back the d Alembert solution. However, the Fourier form says something very important: the solution of the wave equation can be written as a superposition of plave waves with constant group velocity v g = dω/dk = v, so there is no dispersion. Lucio Demeio - DIISM wave equation 27 / 44
58 Solutions by Fourier Transforms The Klein-Gordon Equation Again, let u(x, t) = û(k, t) e i k x dk û(k, t) = 1 2π and substitute in the Klein-Gordon equation. One obtains 2 û t 2 + ω(k)2 û = 0 with ω(k) = k 2 v 2 + γ. The general solution is again given by û(k, t) = A(k) e iωt + B(k)e iωt u(x, t) e i k x dx Lucio Demeio - DIISM wave equation 28 / 44
59 Solutions by Fourier Transforms The Klein-Gordon Equation Again, let u(x, t) = û(k, t) e i k x dk û(k, t) = 1 2π and substitute in the Klein-Gordon equation. One obtains 2 û t 2 + ω(k)2 û = 0 with ω(k) = k 2 v 2 + γ. The general solution is again given by û(k, t) = A(k) e iωt + B(k)e iωt u(x, t) e i k x dx Now the group velocity v g = dω/dk is not constant and there is dispersion: waves with different k s propagate at different speeds and the shape (signal) acquires a distorsion in time. Lucio Demeio - DIISM wave equation 28 / 44
60 Solutions by Fourier Transforms The telegrapher s equation Again, with u(x, t) = and γ = 0, one obtains û(k, t) e i k x dk û(k, t) = 1 2π 2 û t 2 + 2λ û t + ω(k)2 û = 0 with ω(k) = k v. The general solution is again given by u(x, t) e i k x dx û(k, t) = e λt [ A 1 (k) e Ωt + A 2 (k) e Ωt] for k λ v = e λt [ B 1 (k) e iνt + B 2 (k)e iνt] for k > λ v where Ω = λ 2 k 2 v 2 and ν = k 2 v 2 λ 2. Lucio Demeio - DIISM wave equation 29 / 44
61 The solution shows both dispersion and diffusion. Lucio Demeio - DIISM wave equation 29 / 44 Solutions by Fourier Transforms The telegrapher s equation Again, with u(x, t) = and γ = 0, one obtains û(k, t) e i k x dk û(k, t) = 1 2π 2 û t 2 + 2λ û t + ω(k)2 û = 0 with ω(k) = k v. The general solution is again given by u(x, t) e i k x dx û(k, t) = e λt [ A 1 (k) e Ωt + A 2 (k) e Ωt] for k λ v = e λt [ B 1 (k) e iνt + B 2 (k)e iνt] for k > λ v where Ω = λ 2 k 2 v 2 and ν = k 2 v 2 λ 2.
62 An example Initial condition We compare solutions of the wave equation, the Klein-Gordon equation and the telegrapher s equation with initial condition u(x, 0) = e x2 /2 u (x, 0) = 0 t (black line: wave eq., red line: Klein-Gordon, blue line: telegr.) Lucio Demeio - DIISM wave equation 30 / 44
63 An example Initial condition We compare solutions of the wave equation, the Klein-Gordon equation and the telegrapher s equation with initial condition u(x, 0) = e x2 /2 u (x, 0) = 0 t (black line: wave eq., red line: Klein-Gordon, blue line: telegr.) Lucio Demeio - DIISM wave equation 30 / 44
64 The wave equation on a finite domain Lucio Demeio - DIISM wave equation 31 / 44
65 The wave equation on a finite domain Vector space with scalar product Consider now the wave equation in the interval 0 x l. On a finite domain, we must impose two boundary conditions (the equation is of 2 nd order in x). Let u(0, t) = u(l, t) = 0 (25) be the homogeneous Dirichlet boundary conditions and consider the linear operator L u = 2 u x 2 defined on all twice differentiable functions on [0, l] with boundary conditions (25). Equipped with the scalar product (u, v) = l 0 u(x) v(x) dx it becomes a vector space with scalar product (pre-hilbert space) and it is easily seen that L is a self-adjoint operator w.r. to this scalar product. Lucio Demeio - DIISM wave equation 32 / 44
66 The wave equation on a finite domain Eigenfunctions Then, L has a set of real eigenvalues, kn 2 = (nπ/l) 2, for n = 1, 2,... and real orthogonal eigenfunctions φ n (x), given by 2 φ n (x) = l sin k n x (26) which form an orthonormal basis for the vector space, that is (φ n, φ m ) = l 0 φ n (x) φ m (x) dx = δ mn (27) Lucio Demeio - DIISM wave equation 33 / 44
67 The wave equation on a finite domain Eigenfunctions Then, L has a set of real eigenvalues, kn 2 = (nπ/l) 2, for n = 1, 2,... and real orthogonal eigenfunctions φ n (x), given by 2 φ n (x) = l sin k n x (26) which form an orthonormal basis for the vector space, that is (φ n, φ m ) = l 0 φ n (x) φ m (x) dx = δ mn (27) Eigenfunction expansion Then, any function of this vector space can be expressed as 2 u(x, t) = c n (t) φ n (x) = c n (t) l sin k n x (28) n=1 Lucio Demeio - DIISM wave equation 33 / 44 n=1
68 The wave equation on a finite domain Solution of the wave equation By substituting (28) into the wave equation (2) we obtain { cn(t) φ n(x) v 2 c n(t) φ n(x) } = 0 n=1 { cn(t) φ n(x) + v 2 c n(t) L φ } n(x) = 0 n=1 { cn(t) φ n(x) + v 2 c n(t) kn 2 φ } n(x) = 0 n=1 { cn(t) + ωn 2 c } n(t) φ n(x) = 0 with ω n = v k n n=1 By taking the scalar product with φ m we obtain c m(t) + ωm 2 c m(t) = 0 whose general solution is c m(t) = A m cos ω mt + B m sin ω mt Lucio Demeio - DIISM wave equation 34 / 44
69 The wave equation on a finite domain Solution of the wave equation The general solution of the wave equation (2) can then be written as a superposition of all eigenfunctions as u(x, t) = {A n cos ω nt + B n sin ω nt} φ n(x) (29) n=1 Each of the terms in the sum is called mode of vibration or free vibrations. The general solution is thus a linear superposition of the vibration modes. Lucio Demeio - DIISM wave equation 35 / 44
70 The wave equation on a finite domain Solution of the wave equation The general solution of the wave equation (2) can then be written as a superposition of all eigenfunctions as u(x, t) = {A n cos ω nt + B n sin ω nt} φ n(x) (29) n=1 Each of the terms in the sum is called mode of vibration or free vibrations. The general solution is thus a linear superposition of the vibration modes. Initial conditions The wave equation must be accompanied by the initial conditions (14)-(15): u(x, 0) = h(x) u (x, 0) = ψ(x) t from which the coefficents A n and B n can be determined: A n = l 0 h(x) φ n(x) dx B n = 1 ω n l 0 ψ(x) φ n(x) dx Lucio Demeio - DIISM wave equation 35 / 44
71 The wave equation on a finite domain The Klein-Gordon equation and the telegrapher s equation By following similar steps, we may write the general solutions of the Klein-Gordon equation (23) and the telegrapher s equation (24): u(x, t) = {A n cos ω n t + B n sin ω n t} φ n (x) (K.-G.) (30) n=1 u(x, t) = e λt n=1 {A n cos ν n t + B n sin ν n t} φ n (x) (telegr.) (31) where the frequencies now are ω n = v 2 k 2 n + γ and ν n = v 2 k 2 n λ 2. We again observe dispersion in the Klein-Gordon equation and both dispersion and diffusion in the telegrapher s equation. Lucio Demeio - DIISM wave equation 36 / 44
72 The wave equation on a finite domain The Klein-Gordon equation and the telegrapher s equation By following similar steps, we may write the general solutions of the Klein-Gordon equation (23) and the telegrapher s equation (24): u(x, t) = {A n cos ω n t + B n sin ω n t} φ n (x) (K.-G.) (30) n=1 u(x, t) = e λt n=1 {A n cos ν n t + B n sin ν n t} φ n (x) (telegr.) (31) where the frequencies now are ω n = v 2 k 2 n + γ and ν n = v 2 k 2 n λ 2. We again observe dispersion in the Klein-Gordon equation and both dispersion and diffusion in the telegrapher s equation. Numerical simulations FiniteDomain.W, FiniteDomain.KG, FiniteDomain.tele, Lucio Demeio - DIISM wave equation 36 / 44
73 Non homogeneous problems Lucio Demeio - DIISM wave equation 37 / 44
74 Non homogeneous problems Generalities Non homogeneous terms arise whenever nonhomogeneous boundary conditions or external forces are present. As customary in these cases, the solution is a sum of the general solution of the homogeneous problem and a particular solution of the complete equation. Lucio Demeio - DIISM wave equation 38 / 44
75 Non homogeneous problems Generalities Non homogeneous terms arise whenever nonhomogeneous boundary conditions or external forces are present. As customary in these cases, the solution is a sum of the general solution of the homogeneous problem and a particular solution of the complete equation. Nonhomogeneous boundary conditions 2 u t 2 u 2 v2 x = 0 2 u(0, t) = 0 u(l, t) = A General solution: u(x, t) = u p(x) + ũ(x, t) with u p(x) and ũ such that v 2 u p = 0 u p(0) = 0 stationary solution u p(l) = A which gives u p(x) = A x/l. 2 ũ t 2 ũ 2 v2 x = 0 2 ũ(0, t) = ũ(l, t) = 0 Lucio Demeio - DIISM wave equation 38 / 44
76 Non homogeneous problems The solution then is: u(x, t) = A x l + {A n cos ω n t + B n sin ω n t} φ n (x) n=1 where now the coefficients A n and B n are given by A n = l 0 [ h(x) A x ] φ n (x) dx B n = 1 l ω n l 0 ψ(x) φ n (x) dx Lucio Demeio - DIISM wave equation 39 / 44
77 Non homogeneous problems The solution then is: u(x, t) = A x l + {A n cos ω n t + B n sin ω n t} φ n (x) n=1 where now the coefficients A n and B n are given by A n = l 0 [ h(x) A x ] φ n (x) dx B n = 1 l ω n l 0 ψ(x) φ n (x) dx Numerical simulations FiniteDomain.WNH Lucio Demeio - DIISM wave equation 39 / 44
78 Non homogeneous problems External loads, e.g. gravity 2 u t 2 v2 2 u x 2 + k = 0 u(0, t) = 0 u(l, t) = A General solution: u(x, t) = u p (x) + ũ(x, t) with u p (x) and ũ such that v 2 u p + k = 0 stationary solution u p (0) = 0 u p (l) = 0 which gives u p (x) = k x (x l) 2 v2 2 ũ t 2 v2 2 ũ x 2 = 0 ũ(0, t) = ũ(l, t) = 0 Lucio Demeio - DIISM wave equation 40 / 44
79 Non homogeneous problems The solution then is: u(x, t) = k x 2 v 2 (x l) + {A n cos ω n t + B n sin ω n t} φ n (x) n=1 where now the coefficients A n and B n are given by A n = l 0 [ h(x) k x ] (x l) 2 v2 φ n (x) dx B n = 1 ω n l 0 ψ(x) φ n (x) dx Lucio Demeio - DIISM wave equation 41 / 44
80 Non homogeneous problems The solution then is: u(x, t) = k x 2 v 2 (x l) + {A n cos ω n t + B n sin ω n t} φ n (x) n=1 where now the coefficients A n and B n are given by A n = l 0 [ h(x) k x ] (x l) 2 v2 φ n (x) dx B n = 1 ω n l 0 ψ(x) φ n (x) dx Numerical simulations FiniteDomain.WNH Lucio Demeio - DIISM wave equation 41 / 44
81 Non homogeneous problems Exteral excitations and resonances Resonance phenomena occur when harmonic external sources (forcing terms) or harmonic boundary conditions are present: Sources: 2 u t 2 u 2 v2 = f(x, t) x2 u(0, t) = 0 u(l, t) = 0 Non hom. b.c.: 2 u t 2 u 2 v2 x = 0 2 u(0, t) = 0 u(l, t) = f(t) Consider the case of harmonic boundary conditions with f(t) = A 0 sin µt. The unknown function can be written as u(x, t) = A 0 x l sin µt + ũ(x, t) Then, ũ obeys the equation (looks like a problem with source) 2 ũ t 2 ũ 2 v2 x = x 2 A0 l µ2 sin µt (32) ũ(0, t) = ũ(l, t) = 0 Lucio Demeio - DIISM wave equation 42 / 44
82 Non homogeneous problems External excitation We expand ũ in the usual eigenfunctions: ũ(x, t) = After substituting into equation (32): n=1 c n(t) φ n(x). n=1 { cn(t) + ω 2 n c n(t) } φ n(x) = A 0 x l µ2 sin µt with ω n = v k n and, by taking scalar products with the φ m s, c m(t) + ω 2 m c m(t) = A 0 γ m l µ 2 sin µt with γ m = (x, φ m(x)) whose general solution is c m(t) = A m cos ω mt + B m sin ω mt + A0 l γ m µ 2 sin µt ωm 2 µ 2 Lucio Demeio - DIISM wave equation 43 / 44
83 Non homogeneous problems Solution of the wave equation The general solution of the wave equation (2) can then be written as a superposition of all eigenfunctions as u(x, t) = A 0 x l sin µt + n=1 { A n cos ω nt + B n sin ω nt + A0 l } γ n µ 2 sin µt φ ωn 2 µ 2 n(x) With the initial conditions u(x, 0) = h(x) we have for the coefficients A n and B n: u (x, 0) = ψ(x) t A n = l 0 B n = 1 ω n h(x) φ n(x) dx { l [ 0 ψ(x) A 0 µ x l Numerical simulations: FiniteDomain.WNH ] A0 µ γn φ n(x) dx l (1 + µ2 ω 2 n µ 2 )} Lucio Demeio - DIISM wave equation 44 / 44
Separation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationENGI 4430 PDEs - d Alembert Solutions Page 11.01
ENGI 4430 PDEs - d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationBoundary-value Problems in Rectangular Coordinates
Boundary-value Problems in Rectangular Coordinates 2009 Outline Separation of Variables: Heat Equation on a Slab Separation of Variables: Vibrating String Separation of Variables: Laplace Equation Review
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationThere are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.
Preliminary Examination (Solutions): Partial Differential Equations, 1 AM - 1 PM, Jan. 18, 16, oom Discovery Learning Center (DLC) Bechtel Collaboratory. Student ID: There are five problems. Solve four
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More information2 u 1-D: 3-D: x + 2 u
c 2013 C.S. Casari - Politecnico di Milano - Introduction to Nanoscience 2013-14 Onde 1 1 Waves 1.1 wave propagation 1.1.1 field Field: a physical quantity (measurable, at least in principle) function
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationBOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES
1 BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES 1.1 Separable Partial Differential Equations 1. Classical PDEs and Boundary-Value Problems 1.3 Heat Equation 1.4 Wave Equation 1.5 Laplace s Equation
More informationPDE and Boundary-Value Problems Winter Term 2014/2015
PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 12 Saarland University 15. Dezember 2014 c Daria Apushkinskaya (UdS) PDE and BVP lecture 12 15. Dezember 2014 1 / 24 Purpose of Lesson To introduce
More informationName: Math Homework Set # 5. March 12, 2010
Name: Math 4567. Homework Set # 5 March 12, 2010 Chapter 3 (page 79, problems 1,2), (page 82, problems 1,2), (page 86, problems 2,3), Chapter 4 (page 93, problems 2,3), (page 98, problems 1,2), (page 102,
More informationVibrating Strings and Heat Flow
Vibrating Strings and Heat Flow Consider an infinite vibrating string Assume that the -ais is the equilibrium position of the string and that the tension in the string at rest in equilibrium is τ Let u(,
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationLecture6. Partial Differential Equations
EP219 ecture notes - prepared by- Assoc. Prof. Dr. Eser OĞAR 2012-Spring ecture6. Partial Differential Equations 6.1 Review of Differential Equation We have studied the theoretical aspects of the solution
More informationMA 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
More information1 Basic Second-Order PDEs
Partial Differential Equations A. Visintin a.a. 2011-12 These pages are in progress. They contain: an abstract of the classes; notes on some (few) specific issues. These notes are far from providing a
More informationHyperbolic PDEs. Chapter 6
Chapter 6 Hyperbolic PDEs In this chapter we will prove existence, uniqueness, and continuous dependence of solutions to hyperbolic PDEs in a variety of domains. To get a feel for what we might expect,
More informationObservation problems related to string vibrations. Outline of Ph.D. Thesis
Observation problems related to string vibrations Outline of Ph.D. Thesis András Lajos Szijártó Supervisor: Dr. Ferenc Móricz Emeritus Professor Advisor: Dr. Jenő Hegedűs Associate Professor Doctoral School
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationMathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 2 The wave equation Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.0 28/09/2018 1 Learning objectives of this lecture Understand the fundamental properties of the wave equation
More informationWave Equation Modelling Solutions
Wave Equation Modelling Solutions SEECS-NUST December 19, 2017 Wave Phenomenon Waves propagate in a pond when we gently touch water in it. Wave Phenomenon Our ear drums are very sensitive to small vibrations
More informationBASIC WAVE CONCEPTS. Reading: Main 9.0, 9.1, 9.3 GEM 9.1.1, Giancoli?
1 BASIC WAVE CONCEPTS Reading: Main 9.0, 9.1, 9.3 GEM 9.1.1, 9.1.2 Giancoli? REVIEW SINGLE OSCILLATOR: The oscillation functions you re used to describe how one quantity (position, charge, electric field,
More information1. Partial differential equations. Chapter 12: Partial Differential Equations. Examples. 2. The one-dimensional wave equation
1. Partial differential equations Definitions Examples A partial differential equation PDE is an equation giving a relation between a function of two or more variables u and its partial derivatives. The
More informationMath 124A October 11, 2011
Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although
More informationSolving the Heat Equation (Sect. 10.5).
Solving the Heat Equation Sect. 1.5. Review: The Stationary Heat Equation. The Heat Equation. The Initial-Boundary Value Problem. The separation of variables method. An example of separation of variables.
More informationCHAPTER 4. Introduction to the. Heat Conduction Model
A SERIES OF CLASS NOTES FOR 005-006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 4 A COLLECTION OF HANDOUTS ON PARTIAL DIFFERENTIAL EQUATIONS
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationSECTION (See Exercise 1 for verification when both boundary conditions are Robin.) The formal solution of problem 6.53 is
6.6 Properties of Parabolic Partial Differential Equations SECTION 6.6 265 We now return to a difficulty posed in Chapter 4. In what sense are the series obtained in Chapters 4 and 6 solutions of their
More informationPartial Differential Equations
Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but
More informationMA Chapter 10 practice
MA 33 Chapter 1 practice NAME INSTRUCTOR 1. Instructor s names: Chen. Course number: MA33. 3. TEST/QUIZ NUMBER is: 1 if this sheet is yellow if this sheet is blue 3 if this sheet is white 4. Sign the scantron
More informationSeparation of variables in two dimensions. Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 )
Separation of variables in two dimensions Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 ) Separation of variables in two dimensions Overview of method: Consider linear, homogeneous
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
More informationWave Equation With Homogeneous Boundary Conditions
Wave Equation With Homogeneous Boundary Conditions MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 018 Objectives In this lesson we will learn: how to solve the
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More informationA Guided Tour of the Wave Equation
A Guided Tour of the Wave Equation Background: In order to solve this problem we need to review some facts about ordinary differential equations: Some Common ODEs and their solutions: f (x) = 0 f(x) =
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationApplications of the Maximum Principle
Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More information(The) Three Linear Partial Differential Equations
(The) Three Linear Partial Differential Equations 1 Introduction A partial differential equation (PDE) is an equation of a function of 2 or more variables, involving 2 or more partial derivatives in different
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More informationTopic: Introduction to Green s functions
Lecture notes on Variational and Approximate Methods in Applied Mathematics - A Peirce UBC 1 Topic: Introduction to Green s functions (Compiled 2 September 212) In this lecture we provide a brief introduction
More informationMETHODS OF ENGINEERING MATHEMATICS
METHODS OF ENGINEERING MATHEMATICS Edward J. Hang Kyung K. Choi Department of Mechanical Engineering College of Engineering The University of Iowa Iowa City, Iowa 52242 METHODS OF ENGINEERING MATHEMATICS
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationOn the Standard Linear Viscoelastic model
On the Standard Linear Viscoelastic model M. Pellicer (Universitat de Girona) Work in collaboration with: J. Solà-Morales (Universitat Politècnica de Catalunya) (bounded problem) B. Said-Houari (Alhosn
More informationMATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:
MATH 251 Final Examination December 19, 2012 FORM A Name: Student Number: Section: This exam has 17 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
More informationMA 201: Method of Separation of Variables Finite Vibrating String Problem Lecture - 11 MA201(2016): PDE
MA 201: Method of Separation of Variables Finite Vibrating String Problem ecture - 11 IBVP for Vibrating string with no external forces We consider the problem in a computational domain (x,t) [0,] [0,
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
More informationMA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation
MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation Lecture 12 Lecture 12 MA 201, PDE (2016) 1 / 24 Formal Solution of
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationSuperposition of electromagnetic waves
Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationSturm-Liouville Theory
More on Ryan C. Trinity University Partial Differential Equations April 19, 2012 Recall: A Sturm-Liouville (S-L) problem consists of A Sturm-Liouville equation on an interval: (p(x)y ) + (q(x) + λr(x))y
More informationBoundary value problems for partial differential equations
Boundary value problems for partial differential equations Henrik Schlichtkrull March 11, 213 1 Boundary value problem 2 1 Introduction This note contains a brief introduction to linear partial differential
More informationDr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationAPPM GRADUATE PRELIMINARY EXAMINATION PARTIAL DIFFERENTIAL EQUATIONS SOLUTIONS
Thursday August 24, 217, 1AM 1PM There are five problems. Solve any four of the five problems. Each problem is worth 25 points. On the front of your bluebook please write: (1) your name and (2) a grading
More informationIn this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.
Eigen Function Expansion and Applications. In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation. a/ The theory. b/ Example: Solving the Euler equation in two ways.
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationFundamentals of Fluid Dynamics: Waves in Fluids
Fundamentals of Fluid Dynamics: Waves in Fluids Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/ tzielins/ Institute
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More informationTP computing lab - Integrating the 1D stationary Schrödinger eq
TP computing lab - Integrating the 1D stationary Schrödinger equation September 21, 2010 The stationary 1D Schrödinger equation The time-independent (stationary) Schrödinger equation is given by Eψ(x)
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationAdditional Homework Problems
Additional Homework Problems These problems supplement the ones assigned from the text. Use complete sentences whenever appropriate. Use mathematical terms appropriately. 1. What is the order of a differential
More information