Mathematical Methods - Lecture 9

Size: px
Start display at page:

Download "Mathematical Methods - Lecture 9"

Transcription

1 Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, Tel.: +33 (0) yuliya.tarabalka@inria.fr Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 1 / 33

2 Outline 1 Partial Differentiation 2 Partial Differential Equations Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 2 / 33

3 Partial Differentiation Functions of several variables f (x) is a function of one variable x We may consider functions that depend on more than one variable Example: f (x, y) = x 2 + 3xy depends on 2 variables x and y For any pair of values x, y, f (x, y) has a well-defined value Function f (x 1, x 2,..., x n ) depends on the variables x 1, x 2,..., x n Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 3 / 33

4 Partial derivatives Partial Differentiation A function f (x, y) will have a gradient in ALL directions in the xy-plane Gradient = rate of change of a function Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 4 / 33

5 Partial derivatives Partial Differentiation A function f (x, y) will have a gradient in ALL directions in the xy-plane Gradient = rate of change of a function The rates of change of f (x, y) in the positive x- and y- directions are called partial derivatives wrt x and y, respectively Partial derivatives are extremely important in a wide range of applications! Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 4 / 33

6 Partial Differentiation Partial derivatives Partial derivative of f (x, y) with respect to x is denoted by f / x: to signify that a derivative is wrt x, but to recognize that a derivative wrt y also exists f x = lim x 0 provided that the limit exists f (x + x, y) f (x, y), x Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 5 / 33

7 Partial Differentiation Partial derivatives Partial derivative of f (x, y) with respect to x is denoted by f / x: to signify that a derivative is wrt x, but to recognize that a derivative wrt y also exists f x = lim x 0 provided that the limit exists f (x + x, y) f (x, y), x The partial derivative of f with respect to y: f y = lim y 0 provided that the limit exists f (x, y + y) f (x, y), y Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 5 / 33

8 Partial Differentiation Partial derivatives - Notations Partial derivative of f (x, y) with respect to x: f x = ( f x ) = f x y Partial derivative of f with respect to y: f y = ( f y ) = f y x It is important when using partial derivatives to remember which variables are being held constant! Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 6 / 33

9 Partial derivatives Partial Differentiation The extension to the general n-variable case is straightforward: f (x 1,..., x n ) x i provided that the limit exists f (x 1,..., x i + x i,..., x n ) f (x 1,..., x i,..., x n ) = lim, x i 0 x i Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 7 / 33

10 Partial derivatives Partial Differentiation Second (and higher) partial derivatives may be defined in a similar way: x ( f x ) = 2 f x 2 = f xx x ( f y ) = 2 f x y = f xy y ( f y ) = 2 f y 2 = f yy Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 8 / 33

11 Partial Differentiation Properties of partial derivatives Provided that the second partial derivatives are continuous at the point in question: 2 f x y = 2 f y x 2 f x i x j = 2 f x j x i Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 9 / 33

12 Partial Differentiation Partial derivatives - Example Exercise: Find the first and second partial derivatives of the function f (x, y) = 2x 3 y 2 + y 3 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 10 / 33

13 Partial Differentiation Partial derivatives - Example Exercise: Find the first and second partial derivatives of the function f (x, y) = 2x 3 y 2 + y 3 The first partial derivatives: f x = 6x 2 y 2, f y = 4x 3 y + 3y 2 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 10 / 33

14 Partial Differentiation Partial derivatives - Example Exercise: Find the first and second partial derivatives of the function f (x, y) = 2x 3 y 2 + y 3 The first partial derivatives: f x = 6x 2 y 2, f y = 4x 3 y + 3y 2 The second partial derivatives: 2 f x 2 = 12xy 2, 2 f y 2 = 4x 3 + 6y, Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 10 / 33

15 Partial Differentiation Partial derivatives - Example Exercise: Find the first and second partial derivatives of the function f (x, y) = 2x 3 y 2 + y 3 The first partial derivatives: f x = 6x 2 y 2, f y = 4x 3 y + 3y 2 The second partial derivatives: 2 f x 2 = 12xy 2, 2 f y 2 = 4x 3 + 6y, 2 f x y = 2 f y x = 12x 2 y Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 10 / 33

16 The chain rule Partial Differentiation The total derivative of f (x, y) with respect to x: df dx = f dx x dx + f dy y dx The total differential of the function f (x, y): The chain rule: df = f f dx + x y dy df du = f dx x du + f dy y du Particularly useful when an equation is expressed in a parametric form Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 11 / 33

17 The chain rule Partial Differentiation Exercise: Given that x(u) = 1 + au and y(u) = bu 3, find the rate of change of f (x, y) = xe y with respect to u Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 12 / 33

18 The chain rule Partial Differentiation Exercise: Given that x(u) = 1 + au and y(u) = bu 3, find the rate of change of f (x, y) = xe y with respect to u Using the chain rule: Substituting for x and y: df du = (e y )a + ( xe y )3bu 2 df du = e bu3 (a 3bu 2 3bau 3 ) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 12 / 33

19 What is a partial differential equation? Partial differential equation (PDE) is a differential equation that contains unknown multivariable functions and their partial derivatives PDEs are used to describe a wide variety of phenomena: sound, heat, electrodynamics, quantum mechanics... Examples: 1 transport: u x + u y = 0 2 shock wave: u x + uu y = 0 3 vibrating bar: u tt + u xxxx = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 13 / 33

20 Partial differential equation The most general second-order PDE in two independent variables is: F (x, y, u, u x, u y, u xx, u xy, u yy ) = 0 A solution of a PDE is a function y(x, y,...) that satisfies the equation identically, at least in some region of the x, y,... variables Most of the important PDEs of physics are second-order and linear For linear homogeneous PDEs, the superposition principle applies Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 14 / 33

21 Partial differential equation - Example Exercise: Find all u(x, y) satisfying the equation u xx = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 15 / 33

22 Partial differential equation - Example Exercise: Find all u(x, y) satisfying the equation u xx = 0 We could integrate once to get u x = constant? Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 15 / 33

23 Partial differential equation - Example Exercise: Find all u(x, y) satisfying the equation u xx = 0 We could integrate once to get u x = constant? Since there is another variable y, that s not really right. We get: u x (x, y) = f (y), where f (y) is arbitrary Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 15 / 33

24 Partial differential equation - Example Exercise: Find all u(x, y) satisfying the equation u xx = 0 We could integrate once to get u x = constant? Since there is another variable y, that s not really right. We get: where f (y) is arbitrary u x (x, y) = f (y), We integrate again to get u(x, y) = f (y)x + g(y). This is the general solution. Note that there are two arbitrary functions in the solution Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 15 / 33

25 Partial differential equation - Example 2 Exercise: Solve the PDE u xx + u = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 16 / 33

26 Partial differential equation - Example 2 Exercise: Solve the PDE u xx + u = 0 Again, it s an ODE with an extra variable y Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 16 / 33

27 Partial differential equation - Example 2 Exercise: Solve the PDE u xx + u = 0 Again, it s an ODE with an extra variable y The solution is u = f (y) cos x + g(y) sin x, where f (y) and g(y) are two arbitrary functions of y Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 16 / 33

28 Partial differential equation - Example 2 Exercise: Solve the PDE u xx + u = 0 Again, it s an ODE with an extra variable y The solution is u = f (y) cos x + g(y) sin x, where f (y) and g(y) are two arbitrary functions of y Moral: A PDE has arbitrary functions in its solution Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 16 / 33

29 Important PDEs Partial Differential Equations Transport equation: u t + cu x = 0 Wave equation: u tt = c 2 u xx Diffusion equation: u t = ku xx Laplace equation: u xx + u yy = 0. It is often written as: 2 u = 0 or u = 0, where = 2 is the Laplace operator Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 17 / 33

30 Numerical methods to solve PDEs The three most widely used numerical methods: 1 Finite difference method The simplest method to learn and use Functions are represented by their values at certain grid points, and derivatives are approximated through differences in these values 2 Finite element method 3 Finite volume method Other methods exist: method of lines, spectral methods, multigrid methods,... Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 18 / 33

31 Partial Differential Equations Transport equation: intuition and derivation Yuliya Tarabalka Differential Equations 19 / 33

32 Transport equation: general solution Solve the transport equation u t + cu x = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

33 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

34 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 I see a stationary wave Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

35 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 I see a stationary wave moving coordinate ξ = x ct Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

36 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 I see a stationary wave moving coordinate ξ = x ct u(x, t) = v(ξ, t) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

37 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 I see a stationary wave moving coordinate ξ = x ct u(x, t) = v(ξ, t) Derivatives using the chain rule: u t = v ξ ξ t + v t t t v = c ξ + v t Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

38 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 I see a stationary wave moving coordinate ξ = x ct u(x, t) = v(ξ, t) Derivatives using the chain rule: u t = v ξ ξ t + v t t t v = c ξ + v t u x = v ξ ξ x + v t t x = v ξ Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

39 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 I see a stationary wave moving coordinate ξ = x ct u(x, t) = v(ξ, t) Derivatives using the chain rule: u t = v ξ ξ t + v t t t v = c ξ + v t u x = v ξ ξ x + v t t x = v ξ u t + cu x = cv ξ + v t + cv ξ = v t = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

40 Transport equation: general solution Solve the transport equation u t + cu x = 0 You see a wave with speed c: u t + cu x = 0 I see a stationary wave: v t = 0 moving coordinate ξ = x ct u(x, t) = v(ξ, t) Derivatives using the chain rule: u t = v ξ ξ t + v t t t v = c ξ + v t u x = v ξ ξ x + v t t x = v ξ u t + cu x = cv ξ + v t + cv ξ = v t = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 20 / 33

41 Transport equation: general solution Usual coordinates (x, t): u(x, t) u t + cu x = 0 Moving coordinates (ξ, t) v(ξ, t) v t = 0 ξ = x ct Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 21 / 33

42 Transport equation: general solution Usual coordinates (x, t): u(x, t) u t + cu x = 0 Moving coordinates (ξ, t) v(ξ, t) v t = 0 ξ = x ct Procedure: solve v t = 0, use to find u(x, t) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 21 / 33

43 Transport equation: general solution Usual coordinates (x, t): u(x, t) u t + cu x = 0 Moving coordinates (ξ, t) v(ξ, t) v t = 0 ξ = x ct Procedure: solve v t = 0, use to find u(x, t) v t = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 21 / 33

44 Transport equation: general solution Usual coordinates (x, t): u(x, t) u t + cu x = 0 Moving coordinates (ξ, t) v(ξ, t) v t = 0 ξ = x ct Procedure: solve v t = 0, use to find u(x, t) v t = 0 v = F (ξ) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 21 / 33

45 Transport equation: general solution Usual coordinates (x, t): u(x, t) u t + cu x = 0 Moving coordinates (ξ, t) v(ξ, t) v t = 0 ξ = x ct Procedure: solve v t = 0, use to find u(x, t) v t = 0 v = F (ξ) v(ξ, t) = u(x, t) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 21 / 33

46 Transport equation: general solution Usual coordinates (x, t): u(x, t) u t + cu x = 0 Moving coordinates (ξ, t) v(ξ, t) v t = 0 ξ = x ct Procedure: solve v t = 0, use to find u(x, t) v t = 0 v = F (ξ) v(ξ, t) = u(x, t) u(x, t) = F (ξ) = F (x ct) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 21 / 33

47 Transport equation: general solution Usual coordinates (x, t): u(x, t) u t + cu x = 0 Moving coordinates (ξ, t) v(ξ, t) v t = 0 ξ = x ct Procedure: solve v t = 0, use to find u(x, t) v t = 0 v = F (ξ) v(ξ, t) = u(x, t) u(x, t) = F (ξ) = F (x ct) Conclusion: General solution of u t + cu x = 0 is F (x ct), where F is a once differentiable function Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 21 / 33

48 Initial and boundary conditions Because PDEs typically have so many solutions, we single out one solution by imposing auxiliary conditions These conditions are motivated by the considered problem (ex.: physics) They come in two varieties: initial conditions and boundary conditions Yuliya Tarabalka Differential Equations 22 / 33

49 Initial and boundary conditions An initial condition specifies the physical state at a particular time t 0. We can have one or several initial conditions: u(x, t 0 ) = φ(x) = φ(x, y, z) u u(x, t 0 ) = φ(x) and t (x, t 0) = ψ(x), Examples: φ(x) is the initial position, or the initial temperature,... Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 23 / 33

50 Initial and boundary conditions In each physical domain there is a domain D in which the PDE is valid. Example: For the vibrating string, D is the interval 0 < x < l, where l is the length of the string. The boundary of D consists only of the two points x = 0 and x = l To determine the solution, it is necessary to specify some boundary conditions Three most important kinds of boundary conditions: 1 Dirichlet condition: u is specified 2 Neumann condition: the normal derivative u/ n is specified Example: u/ n = g(x, t) 3 Robin condition: u/ n + au is specified, a = a(x, y, z, t) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 24 / 33

51 Initial and boundary conditions Three most important kinds of boundary conditions: 1 Dirichlet condition: u is specified 2 Neumann condition: the normal derivative u/ n is specified Example: u/ n = g(x, t) 3 Robin condition: u/ n + au is specified, a = a(x, y, z, t) In one-dimensional problems the boundary consists of two endpoints boundary conditions take the simple form (D) u(0, t) = g(t) and u(l, t) = h(t) (N) u x u (0, t) = g(t) and (l, t) = h(t) x Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 25 / 33

52 Well-posed problems Partial Differential Equations Well-posed problems consist of: a PDE in a domain a set of initial and/or boundary conditions or other auxiliary conditions that enjoy the following fundamental properties 1 Existence: There exists at least one solution u(x, t) satisfying all conditions 2 Uniqueness: There is at most one solution 3 Stability: If the data are changed a little, the solution changes only a little Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 26 / 33

53 Transport equation: method of characteristics Technique for solving PDE by reducing to ODE Yuliya Tarabalka Differential Equations 27 / 33

54 Transport equation: method of characteristics Technique for solving PDE by reducing to ODE x(t) = position of the moving observer Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 27 / 33

55 Transport equation: method of characteristics Technique for solving PDE by reducing to ODE x(t) = position of the moving observer How does u(x, t) change from observer s perspective? Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 27 / 33

56 Transport equation: method of characteristics Technique for solving PDE by reducing to ODE x(t) = position of the moving observer How does u(x, t) change from observer s perspective? d du u(x(t), t) = dt dt = u dx x dt + u t Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 27 / 33

57 Transport equation: method of characteristics Technique for solving PDE by reducing to ODE x(t) = position of the moving observer How does u(x, t) change from observer s perspective? d du u(x(t), t) = dt dt = u dx x dt + u t 0 = cu x + u t Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 27 / 33

58 Transport equation: method of characteristics Technique for solving PDE by reducing to ODE x(t) = position of the moving observer How does u(x, t) change from observer s perspective? d du u(x(t), t) = dt dt = u dx x dt + u t 0 = cu x + u t dx/dt = c - moving with speed c du/dt = 0 - u not changing Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 27 / 33

59 Transport equation: method of characteristics du dt u t + cu x = 0 = 0 along curves given by dx dt = c dx dt = c x = ct + x 0 Need initial condition u(x, 0) = f (x) Along x ct = x 0, u(x, t) = f (x 0 ) = f (x ct) Solution: u(x, t) = f (x ct) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 28 / 33

60 Transport equation: Exercise Solve the initial value problem (PDE + initial condition): u t + cu x = 1 u(x, 0) = sin x Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 29 / 33

61 Transport equation: Exercise Solve the initial value problem (PDE + initial condition): du dt = u dx x dt + u t u t + cu x = 1 u(x, 0) = sin x Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 29 / 33

62 Transport equation: Exercise Solve the initial value problem (PDE + initial condition): du dt = u dx x dt + u t u t + cu x = 1 u(x, 0) = sin x 1 = u x c + u t Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 29 / 33

63 Transport equation: Exercise Solve the initial value problem (PDE + initial condition): du dt = u dx x dt + u t u t + cu x = 1 u(x, 0) = sin x 1 = u x c + u t dx/dt = c du/dt = 1 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 29 / 33

64 Transport equation: Exercise Solve the initial value problem (PDE + initial condition): du dt = u dx x dt + u t u t + cu x = 1 u(x, 0) = sin x 1 = u x c + u t dx/dt = c du/dt = 1 du dt = 1 u = t + A along characteristic lines Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 29 / 33

65 Transport equation: Exercise Solve the initial value problem (PDE + initial condition): du dt = u dx x dt + u t u t + cu x = 1 u(x, 0) = sin x 1 = u x c + u t dx/dt = c du/dt = 1 du dt = 1 u = t + A along characteristic lines Using an initial condition: u(x, 0) = 0 + sin x A = sin x Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 29 / 33

66 Transport equation: Exercise Solve the initial value problem (PDE + initial condition): du dt = u dx x dt + u t u t + cu x = 1 u(x, 0) = sin x 1 = u x c + u t dx/dt = c du/dt = 1 du dt = 1 u = t + A along characteristic lines Using an initial condition: u(x, 0) = 0 + sin x A = sin x Solution: u(x, t) = t + sin(x ct) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 29 / 33

67 Solving PDE using Laplace transforms Laplace transforms are well suited to solve linear PDE with constant coefficients If we have a function of two variables w(x, t): L{w(x, t)} = W (x, s) = 0 e st w(x, t)dt Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 30 / 33

68 Solving PDE using Laplace transforms Laplace transforms are well suited to solve linear PDE with constant coefficients If we have a function of two variables w(x, t): L{w(x, t)} = W (x, s) = Transform of partial derivatives: 0 e st w(x, t)dt L{w t (x, t)} = sl{w(x, t)} w(x, 0) = sw (x, s) w(x, 0) L{w tt (x, t)} = s 2 L{w(x, t)} sw(x, 0) w t (x, 0) L{w x (x, t)} = L{w(x, t)} x L{w xx (x, t)} = 2 L{w(x, t)} x 2 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 30 / 33

69 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 31 / 33

70 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 We have: L{w x } + 2xL{w t } = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 31 / 33

71 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 We have: L{w x } + 2xL{w t } = 0 Apply TOD to obtain: W x + 2x[sW (x, s) w(x, 0)] = 0 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 31 / 33

72 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 We have: L{w x } + 2xL{w t } = 0 Apply TOD to obtain: W x + 2x[sW (x, s) w(x, 0)] = 0 Apply initial conditions: W x + 2xsW = 0 (ODE, 1st-order, linear) Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 31 / 33

73 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 We have: L{w x } + 2xL{w t } = 0 Apply TOD to obtain: W x + 2x[sW (x, s) w(x, 0)] = 0 Apply initial conditions: W x + 2xsW = 0 (ODE, 1st-order, linear) General solution of the separable equation: W (x, s) = A(s)e sx2 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 31 / 33

74 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 We have: L{w x } + 2xL{w t } = 0 Apply TOD to obtain: W x + 2x[sW (x, s) w(x, 0)] = 0 Apply initial conditions: W x + 2xsW = 0 (ODE, 1st-order, linear) General solution of the separable equation: W (x, s) = A(s)e sx2 Obtain A(s): W (0, s) = L{w(0, t)} = L{t} = 1/s 2 A(s) = 1/s 2 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 31 / 33

75 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 We have: L{w x } + 2xL{w t } = 0 Apply TOD to obtain: W x + 2x[sW (x, s) w(x, 0)] = 0 Apply initial conditions: W x + 2xsW = 0 (ODE, 1st-order, linear) General solution of the separable equation: W (x, s) = A(s)e sx2 Obtain A(s): W (0, s) = L{w(0, t)} = L{t} = 1/s 2 A(s) = 1/s 2 W (x, s) = 1 s 2 e sx2 Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 31 / 33

76 Solving PDE using Laplace transforms: Example Exercise: Solve, via Laplace transforms: w x + 2xw t = 0, w(x, 0) = 0, w(0, t) = t, t 0 W (x, s) = 1 s 2 e sx2 Second shifting theorem: L 1 {e cs F (s)} = u c (t)f (t c), where the Heaviside (unit step) function: 0, t < c; u c (t) = 1, t c F (s) = 1 s 2 f (t) = t w(x, t) = L 1 { 1 s 2 e sx2 } = u x 2(t)[t x 2 ] Yuliya Tarabalka (yuliya.tarabalka@inria.fr) Differential Equations 32 / 33

77 Partial differential equations Other analytical methods exist: Separation of variables Change of variables Fundamental solution You can read more here: W. A. Strauss, Partial Differential Equations: An Introduction, John Wiley & Sons Ltd, 2008 Yuliya Tarabalka Differential Equations 33 / 33

Mathematical Methods - Lecture 7

Mathematical Methods - Lecture 7 Mathematical Methods - Lecture 7 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr

More information

First order wave equations. Transport equation is conservation law with J = cu, u t + cu x = 0, < x <.

First order wave equations. Transport equation is conservation law with J = cu, u t + cu x = 0, < x <. First order wave equations Transport equation is conservation law with J = cu, u t + cu x = 0, < x

More information

Chapter 3 Second Order Linear Equations

Chapter 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 information

Introduction to Partial Differential Equations

Introduction 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

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE

CLASSIFICATION 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 information

Final: Solutions Math 118A, Fall 2013

Final: 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 information

Math 212-Lecture 8. The chain rule with one independent variable

Math 212-Lecture 8. The chain rule with one independent variable Math 212-Lecture 8 137: The multivariable chain rule The chain rule with one independent variable w = f(x, y) If the particle is moving along a curve x = x(t), y = y(t), then the values that the particle

More information

First order Partial Differential equations

First order Partial Differential equations First order Partial Differential equations 0.1 Introduction Definition 0.1.1 A Partial Deferential equation is called linear if the dependent variable and all its derivatives have degree one and not multiple

More information

Chapter One: Where PDEs come from

Chapter One: Where PDEs come from Chapter One: Where PDEs come from After thinking about the meaning of a partial differential equation, we will flex our mathematical muscles by solving a few of them. Then we will see how naturally they

More information

MATH 425, FINAL EXAM SOLUTIONS

MATH 425, FINAL EXAM SOLUTIONS MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u

More information

Introduction of Partial Differential Equations and Boundary Value Problems

Introduction 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 information

Analysis III (BAUG) Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017

Analysis III (BAUG) Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017 Analysis III (BAUG Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017 Question 1 et a 0,..., a n be constants. Consider the function. Show that a 0 = 1 0 φ(xdx. φ(x = a 0 + Since the integral

More information

INTRODUCTION TO PDEs

INTRODUCTION 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 information

An Introduction to Partial Differential Equations

An Introduction to Partial Differential Equations An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)

More information

Partial Differential Equations, Winter 2015

Partial Differential Equations, Winter 2015 Partial Differential Equations, Winter 215 Homework #2 Due: Thursday, February 12th, 215 1. (Chapter 2.1) Solve u xx + u xt 2u tt =, u(x, ) = φ(x), u t (x, ) = ψ(x). Hint: Factor the operator as we did

More information

Partial Differential Equations

Partial Differential Equations Partial Differential Equations Lecture Notes Dr. Q. M. Zaigham Zia Assistant Professor Department of Mathematics COMSATS Institute of Information Technology Islamabad, Pakistan ii Contents 1 Lecture 01

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 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 information

Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine

Mathé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 information

Math 5587 Lecture 2. Jeff Calder. August 31, Initial/boundary conditions and well-posedness

Math 5587 Lecture 2. Jeff Calder. August 31, Initial/boundary conditions and well-posedness Math 5587 Lecture 2 Jeff Calder August 31, 2016 1 Initial/boundary conditions and well-posedness 1.1 ODE vs PDE Recall that the general solutions of ODEs involve a number of arbitrary constants. Example

More information

Math 4263 Homework Set 1

Math 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 information

The first order quasi-linear PDEs

The first order quasi-linear PDEs Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.

More information

MATH 220: Problem Set 3 Solutions

MATH 220: Problem Set 3 Solutions MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )

More information

Lecture Introduction

Lecture 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 information

Partial Differential Equations

Partial 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 information

LECTURE NOTES FOR MATH 124A PARTIAL DIFFERENTIAL EQUATIONS

LECTURE NOTES FOR MATH 124A PARTIAL DIFFERENTIAL EQUATIONS LECTURE NOTES FOR MATH 124A PARTIAL DIFFERENTIAL EQUATIONS S. SETO 1. Motivation for PDEs 1.1. What are PDEs? An algebraic equation is an equation which only involves algebraic operations, e.g. x 2 1 =.

More information

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: 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 information

Differential equations, comprehensive exam topics and sample questions

Differential 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 information

Introduction and some preliminaries

Introduction and some preliminaries 1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.

More information

Wave Equation With Homogeneous Boundary Conditions

Wave 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 information

Math 311, Partial Differential Equations, Winter 2015, Midterm

Math 311, Partial Differential Equations, Winter 2015, Midterm Score: Name: Math 3, Partial Differential Equations, Winter 205, Midterm Instructions. Write all solutions in the space provided, and use the back pages if you have to. 2. The test is out of 60. There

More information

An Introduction to Numerical Methods for Differential Equations. Janet Peterson

An Introduction to Numerical Methods for Differential Equations. Janet Peterson An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences

More information

PDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation

PDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation PDEs, Homework #3 Solutions. Use Hölder s inequality to show that the solution of the heat equation u t = ku xx, u(x, = φ(x (HE goes to zero as t, if φ is continuous and bounded with φ L p for some p.

More information

HOMEWORK 4 1. P45. # 1.

HOMEWORK 4 1. P45. # 1. HOMEWORK 4 SHUANGLIN SHAO P45 # Proof By the maximum principle, u(x, t x kt attains the maximum at the bottom or on the two sides When t, x kt x attains the maximum at x, ie, x When x, x kt kt attains

More information

MA 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 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 information

Strauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u

Strauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u Strauss PDEs e: Setion 3.4 - Exerise 3 Page 1 of 13 Exerise 3 Solve u tt = u xx + os x, u(x, ) = sin x, u t (x, ) = 1 + x. Solution Solution by Operator Fatorization Bring u xx to the other side. Write

More information

Solving First Order PDEs

Solving First Order PDEs Solving Ryan C. Trinity University Partial Differential Equations January 21, 2014 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.

More information

MATH 220: MIDTERM OCTOBER 29, 2015

MATH 220: MIDTERM OCTOBER 29, 2015 MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and

More information

Lecture Notes in Mathematics. A First Course in Quasi-Linear Partial Differential Equations for Physical Sciences and Engineering Solution Manual

Lecture Notes in Mathematics. A First Course in Quasi-Linear Partial Differential Equations for Physical Sciences and Engineering Solution Manual Lecture Notes in Mathematics A First Course in Quasi-Linear Partial Differential Equations for Physical Sciences and Engineering Solution Manual Marcel B. Finan Arkansas Tech University c All Rights Reserved

More information

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012 MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to

More information

MA 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 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 information

THE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))

THE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t))) THE WAVE EQUATION The aim is to derive a mathematical model that describes small vibrations of a tightly stretched flexible string for the one-dimensional case, or of a tightly stretched membrane for the

More information

Solving First Order PDEs

Solving First Order PDEs Solving Ryan C. Trinity University Partial Differential Equations Lecture 2 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.

More information

Homework for Math , Fall 2016

Homework 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 information

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the 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: 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 information

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference

More information

Diffusion on the half-line. The Dirichlet problem

Diffusion on the half-line. The Dirichlet problem Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection

More information

Classification of partial differential equations and their solution characteristics

Classification 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 information

9 More on the 1D Heat Equation

9 More on the 1D Heat Equation 9 More on the D Heat Equation 9. Heat equation on the line with sources: Duhamel s principle Theorem: Consider the Cauchy problem = D 2 u + F (x, t), on x t x 2 u(x, ) = f(x) for x < () where f

More information

13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)

13 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 information

Chapter 2 Boundary and Initial Data

Chapter 2 Boundary and Initial Data Chapter 2 Boundary and Initial Data Abstract This chapter introduces the notions of boundary and initial value problems. Some operator notation is developed in order to represent boundary and initial value

More information

Homework 3 solutions Math 136 Gyu Eun Lee 2016 April 15. R = b a

Homework 3 solutions Math 136 Gyu Eun Lee 2016 April 15. R = b a Homework 3 solutions Math 136 Gyu Eun Lee 2016 April 15 A problem may have more than one valid method of solution. Here we present just one. Arbitrary functions are assumed to have whatever regularity

More information

Diff. Eq. App.( ) Midterm 1 Solutions

Diff. Eq. App.( ) Midterm 1 Solutions Diff. Eq. App.(110.302) Midterm 1 Solutions Johns Hopkins University February 28, 2011 Problem 1.[3 15 = 45 points] Solve the following differential equations. (Hint: Identify the types of the equations

More information

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt

More information

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known

More information

Lecture No 1 Introduction to Diffusion equations The heat equat

Lecture No 1 Introduction to Diffusion equations The heat equat Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and

More information

Chapter 3. Second Order Linear PDEs

Chapter 3. Second Order Linear PDEs Chapter 3. Second Order Linear PDEs 3.1 Introduction The general class of second order linear PDEs are of the form: ax, y)u xx + bx, y)u xy + cx, y)u yy + dx, y)u x + ex, y)u y + f x, y)u = gx, y). 3.1)

More information

Functions of Several Variables

Functions of Several Variables Functions of Several Variables Partial Derivatives Philippe B Laval KSU March 21, 2012 Philippe B Laval (KSU) Functions of Several Variables March 21, 2012 1 / 19 Introduction In this section we extend

More information

Lecture 19: Heat conduction with distributed sources/sinks

Lecture 19: Heat conduction with distributed sources/sinks Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without explicit written permission from the copyright owner. 1 ecture 19: Heat conduction

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 6. Solve the completely inhomogeneous diffusion problem on the half-line

Strauss PDEs 2e: Section Exercise 2 Page 1 of 6. Solve the completely inhomogeneous diffusion problem on the half-line 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

More information

Tutorial 2. Introduction to numerical schemes

Tutorial 2. Introduction to numerical schemes 236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2

More information

Section 12.6: Non-homogeneous Problems

Section 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 information

Boundary-value Problems in Rectangular Coordinates

Boundary-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 information

6 Non-homogeneous Heat Problems

6 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 information

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C. Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference

More information

Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225

Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225 Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225 Dr. Asmaa Al Themairi Assistant Professor a a Department of Mathematical sciences, University of Princess Nourah bint Abdulrahman,

More information

Applications of the Maximum Principle

Applications 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 information

Strauss PDEs 2e: Section Exercise 1 Page 1 of 6

Strauss PDEs 2e: Section Exercise 1 Page 1 of 6 Strauss PDEs e: Setion.1 - Exerise 1 Page 1 of 6 Exerise 1 Solve u tt = u xx, u(x, 0) = e x, u t (x, 0) = sin x. Solution Solution by Operator Fatorization By fatoring the wave equation and making a substitution,

More information

HOMEWORK 5. Proof. This is the diffusion equation (1) with the function φ(x) = e x. By the solution formula (6), 1. e (x y)2.

HOMEWORK 5. Proof. This is the diffusion equation (1) with the function φ(x) = e x. By the solution formula (6), 1. e (x y)2. HOMEWORK 5 SHUANGLIN SHAO. Section 3.. #. Proof. This is the diffusion equation with the function φx e x. By the solution formula 6, vx, t e x y e x+y φydy e x y e x+y e x y y dy e y dy e x+y y dy To compute

More information

Hyperbolic PDEs. Chapter 6

Hyperbolic 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 information

Propagation of discontinuities in solutions of First Order Partial Differential Equations

Propagation of discontinuities in solutions of First Order Partial Differential Equations Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in

More information

MATH3203 Lecture 1 Mathematical Modelling and ODEs

MATH3203 Lecture 1 Mathematical Modelling and ODEs MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1

More information

Lecture Notes on Numerical Schemes for Flow and Transport Problems

Lecture Notes on Numerical Schemes for Flow and Transport Problems Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute

More information

Final Exam May 4, 2016

Final Exam May 4, 2016 1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.

More information

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables. Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These

More information

MATH 425, HOMEWORK 5, SOLUTIONS

MATH 425, HOMEWORK 5, SOLUTIONS MATH 425, HOMEWORK 5, SOLUTIONS Exercise (Uniqueness for the heat equation on R) Suppose that the functions u, u 2 : R x R t R solve: t u k 2 xu = 0, x R, t > 0 u (x, 0) = φ(x), x R and t u 2 k 2 xu 2

More information

Suggested Solution to Assignment 7

Suggested Solution to Assignment 7 MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the

More information

Lecture 6: Introduction to Partial Differential Equations

Lecture 6: Introduction to Partial Differential Equations Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without explicit written permission from the copyright owner. 1 Lecture 6: Introduction

More information

Lecture Notes on Numerical Schemes for Flow and Transport Problems

Lecture Notes on Numerical Schemes for Flow and Transport Problems Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute

More information

M.Sc. in Meteorology. Numerical Weather Prediction

M.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 information

Heat Equation, Wave Equation, Properties, External Forcing

Heat Equation, Wave Equation, Properties, External Forcing MATH348-Advanced Engineering Mathematics Homework Solutions: PDE Part II Heat Equation, Wave Equation, Properties, External Forcing Text: Chapter 1.3-1.5 ecture Notes : 14 and 15 ecture Slides: 6 Quote

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,

More information

MATH 173: PRACTICE MIDTERM SOLUTIONS

MATH 173: PRACTICE MIDTERM SOLUTIONS MATH 73: PACTICE MIDTEM SOLUTIONS This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve all of them. Write your solutions to problems and in blue book #, and your

More information

Diffusion - The Heat Equation

Diffusion - The Heat Equation Chapter 6 Diffusion - The Heat Equation 6.1 Goal Understand how to model a simple diffusion process and apply it to derive the heat equation in one dimension. We begin with the fundamental conservation

More information

Strauss PDEs 2e: Section Exercise 4 Page 1 of 6

Strauss PDEs 2e: Section Exercise 4 Page 1 of 6 Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where

More information

Before we look at numerical methods, it is important to understand the types of equations we will be dealing with.

Before we look at numerical methods, it is important to understand the types of equations we will be dealing with. Chapter 1. Partial Differential Equations (PDEs) Required Readings: Chapter of Tannehill et al (text book) Chapter 1 of Lapidus and Pinder (Numerical Solution of Partial Differential Equations in Science

More information

A First Course of Partial Differential Equations in Physical Sciences and Engineering

A First Course of Partial Differential Equations in Physical Sciences and Engineering A First Course of Partial Differential Equations in Physical Sciences and Engineering Marcel B. Finan Arkansas Tech University c All Rights Reserved First Draft August 29 2 Preface Partial differential

More information

1. Differential Equations (ODE and PDE)

1. Differential Equations (ODE and PDE) 1. Differential Equations (ODE and PDE) 1.1. Ordinary Differential Equations (ODE) So far we have dealt with Ordinary Differential Equations (ODE): involve derivatives with respect to only one variable

More information

MATH20411 PDEs and Vector Calculus B

MATH20411 PDEs and Vector Calculus B MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411

More information

Partial Differential Equations

Partial 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 information

Some Aspects of Solutions of Partial Differential Equations

Some Aspects of Solutions of Partial Differential Equations Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in

More information

Math 220A - Fall 2002 Homework 5 Solutions

Math 220A - Fall 2002 Homework 5 Solutions Math 0A - Fall 00 Homework 5 Solutions. Consider the initial-value problem for the hyperbolic equation u tt + u xt 0u xx 0 < x 0 u t (x, 0) ψ(x). Use energy methods to show that the domain of dependence

More information

Week 4 Lectures, Math 6451, Tanveer

Week 4 Lectures, Math 6451, Tanveer 1 Diffusion in n ecall that for scalar x, Week 4 Lectures, Math 6451, Tanveer S(x,t) = 1 exp [ x2 4πκt is a special solution to 1-D heat equation with properties S(x,t)dx = 1 for t >, and yet lim t +S(x,t)

More information

x ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4.

x ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4. . The solution is ( 2 e x+ct + e x ct) + 2c x+ct x ct sin(s)dx ( e x+ct + e x ct) + ( cos(x + ct) + cos(x ct)) 2 2c 2. To solve the PDE u xx 3u xt 4u tt =, you can first fact the differential operat to

More information

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics 3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)

More information

ENGI 4430 PDEs - d Alembert Solutions Page 11.01

ENGI 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 information

MATH-UA 263 Partial Differential Equations Recitation Summary

MATH-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 information

Summer 2017 MATH Solution to Exercise 5

Summer 2017 MATH Solution to Exercise 5 Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a

More information

There are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.

There 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 information

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt

More information

Partial Differential Equations Summary

Partial Differential Equations Summary Partial Differential Equations Summary 1. The heat equation Many physical processes are governed by partial differential equations. temperature of a rod. In this chapter, we will examine exactly that.

More information