Heat Equation, Wave Equation, Properties, External Forcing
|
|
- Melvin Norris
- 5 years ago
- Views:
Transcription
1 MATH348-Advanced Engineering Mathematics Homework Solutions: PDE Part II Heat Equation, Wave Equation, Properties, External Forcing Text: Chapter ecture Notes : 14 and 15 ecture Slides: 6 Quote of Homework Solutions PDE Part II Go down in your own way and everyday is the right day and as you rise above the fearlines in his frown you look down and hear the sound of the faces in the crowd. Pink Floyd : Fearless - Meddle (1971 Consider the one-dimensional heat equation, 1. Heat Equation on a closed and bounded spatial domain of R 1+1 (1 ( u t = u c, x x (,, t (,, c = K σρ. Equations (1-( model the time-evolution of the temperature, u = u(x, t, of a heat conducting medium in one-dimension. The object, of length, is assumed to have a homogenous thermal conductivity K, specific heat σ, and linear density ρ. That is, K, σ, ρ R +. If we consider an object of finite-length, positioned on say (,, then we must also specify the boundary conditions 1, (3 u x(, t =, u x(, t =,. astly, for the problem to admit a unique solution we must know the initial configuration of the temperature, (4 u(x, = f(x Separation of Variables : General Solution. Assume that the solution to (1-( is such that u(x, t = F (xg(t and use separation of variables to find the general solution to (1-(, which satisfies (3-(4. (1 Assume that, u(x, t = F (xg(t then u xx = F (xg(t and u t = F (xg (t and the 1-D heat equation becomes, (5 G (t c G(t = F (x F (x = λ, where we have introduced the separation constant λ. 3 From this equation we have the two ODE s, (6 (7 G (t + λc G(t =, F (x + λf (x =. Each of these ODE s can be solved through elementary methods to get, 4 (8 (9 (1 (11 λ R : G(t = Ae λct, A R, λ R + : F (x = c 1 cos( λx + c ( λx, λ R : F (x = c 3 cosh( λ x + c 4 h( λ x, λ = : F (x = c 5 + c 6x. Each of the functions F (x must also satisfy the boundary conditions, u x(, t = and u x(, t and so we won t need all of them. Notice that the boundary conditions imply that, (1 (13 u x(, t = F (G(t =, u x(, t = F (G(t =, 1 Here the boundary conditions correspond to perfect insulation of both endpoints An insulated bar is discussed in examples 4 and 5 on page This occurs in conjunction with the following argument. Since (5 must be true for all (x, t then both sides must be equal to a function that has neither t s nor x s. Hence they must be equal to a constant function. To see that this is true put an x or t on the side that has λ and test points. 4 These elementary methods are those you learned in ODE s and can be found in the ODE review of the lecture slides on separation of variables. 1
2 which gives F ( = and F ( =. 5 So, we now have to determine, which of the previous functions, F (x, satisfy these boundary conditions. To this end we have the following arguments, λ R + : F ( = c 1 λ ( λ + c λ cos( λ = c1 λ + c λ 1 = c =, λ R + : F ( = c 1 λ ( λ + c λ cos( λ = c1 λ ( λ + λ cos( λ = = c 1 λ ( λ = c1 = or ( λ =. If we consider the case that c 1 = then we have F (x = for λ R + but we should try to keep as many solutions as possible and we ignore this case. Thus assume that c 1 we have that ( λ =, which is true for λ = nπ/ and we have the following eigenvalue/eigenfunction pairs indexed by n, (14 F n(x = c n cos( λx, λ n = n π We now consider the λ R case to find that,, n = 1,, 3,.... λ R : F (x = c 3 λ h( λ + c4 λ cosh( λ = c3 + c 4 1 = = c 4 =, λ R : F (x = c 3 λ h( λ + c4 λ cosh( λ = c3 λ h( λ + λ cosh( λ = e λ e λ = c 3 λ = = c 3 =, which means that for λ R we only have the trivial solution F (x =. astly, we consider the case λ = to get, (15 λ = : F ( = F ( = c 6 = = c 5 R, which gives the last eigenpair, 6 (16 F = c R λ =. Noting that there are infinitely many λ s implies now that there are infinitely many temporal solutions (8 and we have, (17 (18 G n(t = A ne λnct, λ = n π, n = 1,, 3,.... For the case where λ = we have the ODE G (t =, whose solution isg (t = A R. Thus we have infinitely many functions, that solve the PDE, of the form: u n(x, t = F n(xg n(t, n =, 1,, 3,... Hence ce the PDE is linear superposition implies that we have the general solution, (19 ( (1 ( u(x, t = u n(x, t = u (x, t + u n(x, t n= = F (xg (x + = c A + = a + F n(xg n(t A nc n cos( λ nxe λnc t a n cos( λ nxe λnct, which is the general solution of the heat equation with the given boundary conditions. 1.. Qualitative Dynamics. Describe how the long term behavior of the general solution to (1-(4 changes as the thermal conductivity, K, is increased while all other parameters are held constant. Also, describe how the solution changes when the linear density, ρ, is increased while all other parameters are held constant. If κ(thermal conductivity is increased, the temporal solution decays faster and the system reaches equilibrium sooner. If ρ(density is increased, the temporal solution decays slower and the system takes longer to reach equilibrium. 5 We assume that G(t = because if it did then we would have u(x, t = F (xg(t = F (x =, which is called the trivial solution and is ignored ce it is already in thermal equilibrium. We care about dynamics! 6 Here we have used the subscripts to denote that these are all associated with the λ = case. We have also trivially changed c 5 to c.
3 Fourier Series : Solution to the IVP. Define, (3 f(x = (4 (5 k x, < x, k ( x, < x < and for the following questions we consider the solution, u, to the heat equation given by, (1-(, which satisfies the initial condition given by (53. 7 For = 1 and k = 1, find the particular solution to (1-( with boundary conditions (3-(4 for when the initial temperature profile of the medium is given by (53. Show that lim t u(x, t = f avg =.5. 8 To find the unknown constants present in the general solution we must apply an initial condition, u(x, = f(x. Doing so gives, u(x, = f(x = a + = a + a n cos( λ nxe λnc a n cos( λ nx, which is a Fourier coe half-range expansion of the initial condition. Thus the unknown constants are Fourier coefficients and, (6 (7 (8 (9 a = 1 a n = f(xdx, f(x cos( λ nxdx. If we note that these integrals have been done in a previous homework on Fourier series, then we known these Fourier coefficients as, a = k, 4k [ ( nπ a n = cos n π x ( 1 n 1. Moreover, if we take = k = 1 we see that lim t u(x, t = a =.5, which implies that under these insulating boundary conditions the equilibrium state for the medium is a constant function u =.5 and that this is nothing more than the average of the initial configuration. Consider the one-dimensional wave equation,. Wave Equation on a closed and bounded spatial domain of R 1+1 (3 (31 u t = c u x, x (,, t (,, c = T ρ. Equations (3-(31 model the time-evolution of the displacement,u = u(x, t, from rest, of an elastic medium in one-dimension. The object, of length, is assumed to have a homogeneous lateral tension T, and linear density ρ. That is, T, ρ R +. Assume, as well, the boundary conditions 9, (3 u x(, t =, u x(, t =, and initial conditions, (33 (34 u(x, = f(x, u t(x, = g(x. 7 When solving the following problems it would be a good idea to go back through your notes and the homework looking for similar calculations. 8 It is interesting here to note that though the initial condition f doesn t appear to satisfy the boundary conditions its periodic Fourier extension does. That is, if you draw the even periodic extension of the initial condition then you would see that the slope is not well defined at the end points. Remembering that the Fourier series averages the right and left hand limits of the periodic extension of the function f at the endpoints shows that the boundary conditions are, in fact, satisfied, ce the derivative of an average is the average of derivatives. 9 These boundary conditions imply that the object must have zero slope at each endpoint.
4 4.1. Separation of Variables : General Solution. Assume that the solution to (3-(31 is such that u(x, t = F (xg(t and use separation of variables to find the general solution to (3-(31, which satisfies (3-(34. The only difference between this problem and the heat problem above are the time-dynamics specified by the PDE. This gives a secondorder constant linear ODE in time and from this ODE we have oscillations of Fourier modes instead of exponential decay (35 Assume that, u(x, t = F (xg(t then u xx = F (xg(t and u tt = F (xg (t and the 1-D wave equation becomes, G (t c G(t = F (x F (x = λ, where we have introduced the separation constant λ. 1 From this equation we have the two ODE s, (36 (37 G (t + λc G(t =, F (x + λf (x =. Each of these ODE s can be solved through elementary methods to get, 13 (38 λ R : G(t = A 1 cos ( c ( λt + A 1 c λt, A 1, A 1 R, (39 (4 (41 λ R + : F (x = c 1 cos( λx + c ( λx, λ R : F (x = c 3 cosh( λ x + c 4 h( λ x, λ = : F (x = c 5 + c 6x. Each of the functions F (x must also satisfy the boundary conditions, u x(, t = and u x(, t and so we won t need all of them. Notice that the boundary conditions imply that, (4 (43 u x(, t = F (G(t =, u x(, t = F (G(t =, which gives F ( = and F ( =. 14 conditions. To this end we have the following arguments, So, we now have to determine, which of the previous functions, F (x, satisfy these boundary λ R + : F ( = c 1 λ ( λ + c λ cos( λ = c1 λ + c λ 1 = c =, λ R + : F ( = c 1 λ ( λ + c λ cos( λ = c1 λ ( λ + λ cos( λ = = c 1 λ ( λ = c1 = or ( λ =. If we consider the case that c 1 = then we have F (x = for λ R + but we should try to keep as many solutions as possible and we ignore this case. Thus assume that c 1 we have that ( λ =, which is true for λ = nπ/ and we have the following eigenvalue/eigenfunction pairs indexed by n, (44 We now consider the λ R case to find that, F n(x = c n cos( λx, λ n = n π, n = 1,, 3,.... λ R : F (x = c 3 λ h( λ + c4 λ cosh( λ = c3 + c 4 1 = = c 4 =, λ R : F (x = c 3 λ h( λ + c4 λ cosh( λ = c3 λ h( λ + λ cosh( λ = e λ e λ = c 3 λ = = c 3 =, which means that for λ R we only have the trivial solution F (x =. astly, we consider the case λ = to get, (45 λ = : F ( = F ( = c 6 = = c 5 R, 1 It is important to notice that the solution to the spatial portion of the problem is the same as the heat problem above. 11 Remember that in this case we have a nontrivial spatial solution for zero eigenvalue. From this you should find the associated temporal function should find that G (t = C 1 + C t. 1 This occurs in conjunction with the following argument. Since (35 must be true for all (x, t then both sides must be equal to a function that has neither t s nor x s. Hence they must be equal to a constant function. To see that this is true put an x or t on the side that has λ and test points. 13 These elementary methods are those you learned in ODE s and can be found in the ODE review of the lecture slides on separation of variables. 14 We assume that G(t = because if it did then we would have u(x, t = F (xg(t = F (x =, which is called the trivial solution and is ignored ce it is already in thermal equilibrium. We care about dynamics!
5 5 which gives the last eigenpair, 15 (46 F = c R λ =. (47 Noting that there are infinitely many λ s implies now that there are infinitely many temporal solutions (38 and we have, ( G n(t = G(t = A n cos c ( λ nt + A n c λ nt, λ n = n π, n = 1,, 3,.... For the case where λ = we have the ODE G (t =, whose solution isg (t = A + A t. Thus we have infinitely many functions, that solve the PDE, of the form: (48 u n(x, t = F n(xg n(t, n =, 1,, 3,... Hence ce the PDE is linear superposition implies that we have the general solution, (49 u(x, t = u n(x, t = u (x, t + u n(x, t (5 (51 (5 n= = F (xg (x + F n(xg n(t = c (A + A t + = a + a t + c n cos( [ ( λ nx A n cos c ( λ nt + A n c λ nt cos( [ ( λ nx a n cos c ( λ nt + a n c λ nt, which is the general solution of the one-dimensional wave equation with the given the boundary conditions... Qualitative Dynamics. Describe how the the general solution to (3-(31 changes as the tension, T, is increased while all other parameters are held constant. Also, describe how the solution changes when the linear density, ρ, is increased while all other parameters are held constant. If T (Tension is increased, the temporal solution oscillate faster. If ρ(density is increased, the temporal oscillates slower..3. Fourier Series : Solution to the IVP. Define, (53 f(x = k x, < x, k ( x, < x <. et = 1 and k = 1 and find the particular solution, which satisfies the initial displacement, f(x, given by (53 and has zero initial velocity for all points on the object. To find the unknown constants present in the general solution we must apply an initial condition, u(x, = f(x. Doing so gives, u(x, = f(x = a + a + cos( [ ( λ nx a n cos c ( λ n + a n c (54 λ n (55 = a + a n cos( λ nx, which is a Fourier coe half-range expansion of the initial condition. Thus the unknown constants are Fourier coefficients and, (56 (57 a = 1 a n = f(xdx, f(x cos( λ nxdx. If we note that these integrals have been done in a previous homework on Fourier series, then we known these Fourier coefficients as, (58 (59 a = k, 4k [ ( nπ a n = cos n π x ( 1 n Here we have used the subscripts to denote that these are all associated with the λ = case. We have also trivially changed c 5 to c.
6 6 Now we apply the initial velocity, u t(x, = g(x = to get that, u t(x, = g(x = a + cos( [ λ nx a ( nc λ n c λ n + a nc ( λ n cos c (6 λ n (61 = a + a nc λ n cos( λ nx which is a Fourier coe half-range expansion of the initial condition. Thus the unknown constants are Fourier coefficients and, (6 (63 a = 1 a n = but ce g(x = we have that a = a n = for all n. 16 c λ n g(xdx, g(x cos( λ nxdx, 3. D alembert Solution to the Wave Equation in R 1+1 Show that by direct substitution that the function u(x, t given by, (64 u(x, t = 1 [u(x ct + u(x + ct + 1 c x+ct x ct v (ydy, is a solution to the one-dimensional wave equation where u and v are the ideally elastic objects initial displacement and velocity, respectively. 17 (65 (66 We must take derivatives of u(x, t. We note the following relations, which follow from the chain-rule. [u(x ct = u [x ct = u c = t [x ct t t [u(x ct = c u (x ct [u(x ct = u x [x ct [x ct x = u = [u(x ct = u x (x ct The derivatives on the second term involving the initial-velocity are more difficult. First, note the following simplification, (67 x+ct x ct v (ydy = x+ct v (ydy x ct Also note the following formula that follows from the fundamental theorem of calculus. (68 Thus, (69 (7 and (71 (7 (73 With this in hand we now have, d dx f(x g(tdt = df dx g(f(x v (ydy. [ x±ct v (ydy = [x ± ct v(x ± ct t t [ x±ct = ±cv (x ± ct = v (ydy = c v (x ± ct, t [ x±ct v (ydy = [x ± ct v(x ± ct x x u tt = c = v (x ± ct = x [ x±ct v (ydy = v (x ± ct. [ u (x ct + u (x + ct + c [ v (x + ct v (x ct (74 c u xx = c ( 1 [ u (x ct + u (x + ct + 1 [ v (x + ct v (x ct, c which implies that u tt c u xx =. 16 It is interesting to note that if a then the displacement grows with time. This is a consequence of the boundary conditions, which require the string to be flat at the endpoints but do not require they stay fixed. 17 d x This is called the D Alembert solution to the wave equation. To do this you may want to recall the fundamental theorem of calculus, f(tdt = dx a a f(x and properties of integrals, f(xdx = f(xdx + f(xdx. a a
7 7 4. Time Dependent Boundary Conditions It makes sense to consider time-dependent interface conditions. That is, (1 and (4 subject to (75 u(, t = g(t, u(, t = h(t, t (, Show that this PDE transforms into: (76 (77 with boundary conditions and initial conditions, (78 (79 where F (x = f(x S(x, and S(x, t = w t = w c St(x, t x, x (,, t (,, c = κ ρσ. h(t + g(t x + g(t. 18 w(, t = w(, t =, w(x, = F (x, We need to find the correct function w(x, t such that (1 and (4 subject to the given boundary conditions transforms to (76-(79. Specifically, we need to find, (8 w(x, t = u(x, t T (x, such that w(, t = w(, t =. Hoping for the best we assume that T is a linear function of x, T (x = ax + b, where a and b are unknown. The left boundary condition on w implies that, (81 (8 (83 w(, t = u(, t T ( = g(t b =, while the right boundary condition on w implies that, (84 (85 (86 w(, t = u(, t T ( = h(t a b =. From this we have that, (87 (88 This gives the form for w as, (89 (9 b = g(t, a = h(t b. ( h(t b w(x, t = u(x, t ( h(t g(t = u(x, t (91 = u(x, t S(x, t, ( h(t g(t where S(x, t = x + g(t. The initial condition on w is now, (9 Finally, we have the transformations, x + g(t w(x, = u(x, S(x,. x + g(t (93 (94 (95 u t = w t + S t, u xx = w xx + S xx = w xx, 18 A similar transformation can be found for the wave equation with inhomogeneous boundary conditions. The moral here is that time-dependent boundary conditions can be transformed into externally driven (AKA Forced or inhomogeneous PDE with standard boundary conditions.
8 8 where S xx(x, t = (96 ( h(t g(t x + g(t =. Thus, the PDE on u becomes, xx w t = w c x S t, which implies that the PDE on u with time-dependent boundary conditions transforms to a PDE on w with the inhomogeneous term S t. 5. Inhomogeneous Wave Equation on a closed and bounded spatial domain of R 1+1 Consider the non-homogeneous one-dimensional wave equation, (97 (98 u t = u c + F (x, t, x x (,, t (,, c = T ρ. with boundary conditions and initial conditions, (99 (1 u(, t = u(, t =, u(x, = u t(x, =. etting F (x, t = A (ωt gives the following Fourier Series Representation of the forcing function F, (11 F (x, t = where ( nπx f n(t, (1 f n(t = A nπ (1 ( 1n (ωt Educated Fourier Series Guesg. Based on the boundary conditions we assume a Fourier e series solution. However, the time-dependence is unclear. So, assume that, (13 u(x, t = ( nπ x G n(t, where G n(t represents the unknown dynamics of the n-th Fourier mode. Ug this assumption and (11-(1, show by direct substitution that (97 yields the ODE, (14 Gn + (15 (16 We have that, u tt c u xx F (x, t = = Gn = A nπ (1 ( 1n (ωt. ( nπ x Gn(t + [ ( nπ Since x for all x we require that, G n(t + Gn(t f n(t ( nπ x G n(t ( nπ x =. f n(t ( nπx (17 (18 G n(t + Gn(t f n(t = G n + Gn = A nπ (1 ( 1n (ωt 5.. Solving for the Dynamics. The solution to (14 is given by, (19 G n(t = G h n(t + G p n(t, ( where G h cnπ ( n(t = B n cos t + Bn cnπ t is the homogeneous solution and G p n(t is the particular solution to (14.
9 Particular Solution - I. If ω cnπ/ then what would the choice for G p n(t be, assuming you were solving for G p n(t ug the method of undetermined coefficients? do not solve for these coefficients If the forcing function to (14 is given by f n(t = A nπ (1 ( 1n (ωt = α n (ωt then we should choose the particular solution G n(t = β n cos(ωt + βn (ωt for ω cnπ/. This would then give the general solution to the inhomogeneous PDE, ( nπ [ ( u(x, t = x B n cos t + Bn cnπ (11 t + β n cos(ωt + βn (ωt, where β n and β n are undetermined coefficients found by standard ODE techniques Particular Solution - II. If ω = cnπ/ then what would the choice for G p n(t be, assuming you were solving for G p n(t ug the method of undetermined coefficients? do not solve for these coefficients If the forcing function to (14 is given by f n(t = A nπ (1 ( 1n (ωt = α n (ωt and ω = cnπ/ then we should choose the particular solution G n(t = β nt cos(ωt + βnt (ωt for ω cnπ/. This would then give the general solution to the inhomogeneous PDE, ( nπ [ ( u(x, t = x B n cos t + Bn cnπ (111 t + β nt cos(ωt + βnt (ωt, where β n and β n are undetermined coefficients found by standard ODE techniques Physical Conclusions. For the Particular Solution - II, what is lim t u(x, t and what does this limit imply physically? There is a t-prefactor on this particular solution, which causes the amplitudes of oscillations to grow as a function of time. Formally, this limit becomes unbounded as t, which is the hallmark of resonance. The moral here is that systems that are prone to oscillate can be made to resonate under the right external force. This is seen in ODE s via mass-spring systems but is just as true for any ideal oscillatory systems.
Math 201 Assignment #11
Math 21 Assignment #11 Problem 1 (1.5 2) Find a formal solution to the given initial-boundary value problem. = 2 u x, < x < π, t > 2 u(, t) = u(π, t) =, t > u(x, ) = x 2, < x < π Problem 2 (1.5 5) Find
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 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 informationMath Assignment 14
Math 2280 - Assignment 14 Dylan Zwick Spring 2014 Section 9.5-1, 3, 5, 7, 9 Section 9.6-1, 3, 5, 7, 14 Section 9.7-1, 2, 3, 4 1 Section 9.5 - Heat Conduction and Separation of Variables 9.5.1 - Solve the
More informationPartial 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 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 informationMethod of Separation of Variables
MODUE 5: HEAT EQUATION 11 ecture 3 Method of Separation of Variables Separation of variables is one of the oldest technique for solving initial-boundary value problems (IBVP) and applies to problems, where
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 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 informationSeparation of Variables. A. Three Famous PDE s
Separation of Variables c 14, Philip D. Loewen A. Three Famous PDE s 1. Wave Equation. Displacement u depends on position and time: u = u(x, t. Concavity drives acceleration: u tt = c u xx.. Heat Equation.
More informationAnalysis 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 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 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 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 informationClass Meeting # 2: The Diffusion (aka Heat) Equation
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 2: The Diffusion (aka Heat) Equation The heat equation for a function u(, x (.0.). Introduction
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 informationReview For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.
Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =
More informationDiffusion 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 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 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 informationMATH 251 Final Examination August 14, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination August 14, 2015 FORM A Name: Student Number: Section: This exam has 11 questions for a total of 150 points. Show all your work! In order to obtain full credit for partial credit
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 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 informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More informationMATH 412 Fourier Series and PDE- Spring 2010 SOLUTIONS to HOMEWORK 5
MATH 4 Fourier Series PDE- Spring SOLUTIONS to HOMEWORK 5 Problem (a: Solve the following Sturm-Liouville problem { (xu + λ x u = < x < e u( = u (e = (b: Show directly that the eigenfunctions are orthogonal
More informationHomework 7 Solutions
Homework 7 Solutions # (Section.4: The following functions are defined on an interval of length. Sketch the even and odd etensions of each function over the interval [, ]. (a f( =, f ( Even etension of
More informationIntroduction to the Wave Equation
Introduction to the Ryan C. Trinity University Partial Differential Equations ecture 4 Modeling the Motion of an Ideal Elastic String Idealizing Assumptions: The only force acting on the string is (constant
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 informationStrauss 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 informationMA22S2 Lecture Notes on Fourier Series and Partial Differential Equations.
MAS Lecture Notes on Fourier Series and Partial Differential Equations Joe Ó hógáin E-mail: johog@maths.tcd.ie Main Text: Kreyszig; Advanced Engineering Mathematics Other Texts: Nagle and Saff, Zill and
More informationMATH20411 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 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 informationCHAPTER 10 NOTES DAVID SEAL
CHAPTER 1 NOTES DAVID SEA 1. Two Point Boundary Value Problems All of the problems listed in 14 2 ask you to find eigenfunctions for the problem (1 y + λy = with some prescribed data on the boundary. To
More informationMATH 251 Final Examination December 16, 2014 FORM A. Name: Student Number: Section:
MATH 2 Final Examination December 6, 204 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 0 points. In order to obtain full credit for partial credit problems, all work must
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 informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 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 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 informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More informationWeek 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 informationLecture 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 informationTHE 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 informationHomework 6 Math 309 Spring 2016
Homework 6 Math 309 Spring 2016 Due May 18th Name: Solution: KEY: Do not distribute! Directions: No late homework will be accepted. The homework can be turned in during class or in the math lounge in Pedelford
More informationLecture 24. Scott Pauls 5/21/07
Lecture 24 Department of Mathematics Dartmouth College 5/21/07 Material from last class The heat equation α 2 u xx = u t with conditions u(x, 0) = f (x), u(0, t) = u(l, t) = 0. 1. Separate variables to
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 informationLECTURE 19: SEPARATION OF VARIABLES, HEAT CONDUCTION IN A ROD
ECTURE 19: SEPARATION OF VARIABES, HEAT CONDUCTION IN A ROD The idea of separation of variables is simple: in order to solve a partial differential equation in u(x, t), we ask, is it possible to find a
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 informationHomework 7 Math 309 Spring 2016
Homework 7 Math 309 Spring 2016 Due May 27th Name: Solution: KEY: Do not distribute! Directions: No late homework will be accepted. The homework can be turned in during class or in the math lounge in Pedelford
More informationMATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination May 4, 2015 FORM A Name: Student Number: Section: This exam has 16 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all work must
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
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 informationMATH 6B Spring 2017 Vector Calculus II Study Guide Final Exam Chapters 8, 9, and Sections 11.1, 11.2, 11.7, 12.2, 12.3.
MATH 6B pring 2017 Vector Calculus II tudy Guide Final Exam Chapters 8, 9, and ections 11.1, 11.2, 11.7, 12.2, 12.3. Before starting with the summary of the main concepts covered in the quarter, here there
More informationLecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems
Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationMath 2a Prac Lectures on Differential Equations
Math 2a Prac Lectures on Differential Equations Prof. Dinakar Ramakrishnan 272 Sloan, 253-37 Caltech Office Hours: Fridays 4 5 PM Based on notes taken in class by Stephanie Laga, with a few added comments
More informationMATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:
MATH 5 Final Examination December 6, 5 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 5 points. In order to obtain full credit for partial credit problems, all work must
More informationu tt = a 2 u xx u tt = a 2 (u xx + u yy )
10.7 The wave equation 10.7 The wave equation O. Costin: 10.7 1 This equation describes the propagation of waves through a medium: in one dimension, such as a vibrating string u tt = a 2 u xx 1 This equation
More informationMath 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 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 informationMath 2930 Worksheet Final Exam Review
Math 293 Worksheet Final Exam Review Week 14 November 3th, 217 Question 1. (* Solve the initial value problem y y = 2xe x, y( = 1 Question 2. (* Consider the differential equation: y = y y 3. (a Find the
More informationTHE UNIVERSITY OF WESTERN ONTARIO. Applied Mathematics 375a Instructor: Matt Davison. Final Examination December 14, :00 12:00 a.m.
THE UNIVERSITY OF WESTERN ONTARIO London Ontario Applied Mathematics 375a Instructor: Matt Davison Final Examination December 4, 22 9: 2: a.m. 3 HOURS Name: Stu. #: Notes: ) There are 8 question worth
More informationMath 2930 Worksheet Wave Equation
Math 930 Worksheet Wave Equation Week 13 November 16th, 017 Question 1. Consider the wave equation a u xx = u tt in an infinite one-dimensional medium subject to the initial conditions u(x, 0) = 0 u t
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 informationA 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 information0 3 x < x < 5. By continuing in this fashion, and drawing a graph, it can be seen that T = 2.
04 Section 10. y (π) = c = 0, and thus λ = 0 is an eigenvalue, with y 0 (x) = 1 as the eigenfunction. For λ > 0 we again have y(x) = c 1 sin λ x + c cos λ x, so y (0) = λ c 1 = 0 and y () = -c λ sin λ
More informationMATH FALL 2014 HOMEWORK 10 SOLUTIONS
Problem 1. MATH 241-2 FA 214 HOMEWORK 1 SOUTIONS Note that u E (x) 1 ( x) is an equilibrium distribution for the homogeneous pde that satisfies the given boundary conditions. We therefore want to find
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 informationThe One-Dimensional Heat Equation
The One-Dimensional Heat Equation R. C. Trinity University Partial Differential Equations February 24, 2015 Introduction The heat equation Goal: Model heat (thermal energy) flow in a one-dimensional object
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 informationFINAL EXAM, MATH 353 SUMMER I 2015
FINAL EXAM, MATH 353 SUMMER I 25 9:am-2:pm, Thursday, June 25 I have neither given nor received any unauthorized help on this exam and I have conducted myself within the guidelines of the Duke Community
More informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
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 informationMATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section:
MATH 5 Final Examination May 3, 07 FORM A Name: Student Number: Section: This exam has 6 questions for a total of 50 points. In order to obtain full credit for partial credit problems, all work must be
More informationLECTURE 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 informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationMATH 425, HOMEWORK 3 SOLUTIONS
MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the
More informationChapter 12 Partial Differential Equations
Chapter 12 Partial Differential Equations Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2 12.1 Basic Concepts of PDEs Partial Differential Equation A
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 informationMath 220a - Fall 2002 Homework 6 Solutions
Math a - Fall Homework 6 Solutions. Use the method of reflection to solve the initial-boundary value problem on the interval < x < l, u tt c u xx = < x < l u(x, = < x < l u t (x, = x < x < l u(, t = =
More informationPeriodic functions: simple harmonic oscillator
Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos
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 informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationMAE 200B Homework #3 Solutions University of California, Irvine Winter 2005
Problem 1 (Haberman 5.3.2): Consider this equation: MAE 200B Homework #3 Solutions University of California, Irvine Winter 2005 a) ρ 2 u t = T 2 u u 2 0 + αu + β x2 t The term αu describes a force that
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation. is called the reduced equation of (N).
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
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 informationIn what follows, we examine the two-dimensional wave equation, since it leads to some interesting and quite visualizable solutions.
ecture 22 igher-dimensional PDEs Relevant section of text: Chapter 7 We now examine some PDEs in higher dimensions, i.e., R 2 and R 3. In general, the heat and wave equations in higher dimensions are given
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 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 informationA Motivation for Fourier Analysis in Physics
A Motivation for Fourier Analysis in Physics PHYS 500 - Southern Illinois University November 8, 2016 PHYS 500 - Southern Illinois University A Motivation for Fourier Analysis in Physics November 8, 2016
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 information1. Solve the boundary-value problems or else show that no solutions exist. y (x) = c 1 e 2x + c 2 e 3x. (3)
Diff. Eqns. Problem Set 6 Solutions. Solve the boundary-value problems or else show that no solutions exist. a y + y 6y, y, y 4 b y + 9y x + e x, y, yπ a Assuming y e rx is a solution, we get the characteristic
More informationWelcome to Math 257/316 - Partial Differential Equations
Welcome to Math 257/316 - Partial Differential Equations Instructor: Mona Rahmani email: mrahmani@math.ubc.ca Office: Mathematics Building 110 Office hours: Mondays 2-3 pm, Wednesdays and Fridays 1-2 pm.
More informationMATH 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 informationMath 216 Second Midterm 16 November, 2017
Math 216 Second Midterm 16 November, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material
More information1. (10 points) Find the general solution to the following second-order differential equation:
Math 307A, Winter 014 Midterm Solutions Page 1 of 8 1. (10 points) Find the general solution to the following second-order differential equation: 4y 1y + 9y = 9t. To find the general solution to this nonhomogeneous
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 informationMidterm 2: Sample solutions Math 118A, Fall 2013
Midterm 2: Sample solutions Math 118A, Fall 213 1. Find all separated solutions u(r,t = F(rG(t of the radially symmetric heat equation u t = k ( r u. r r r Solve for G(t explicitly. Write down an ODE for
More informationkg meter ii) Note the dimensions of ρ τ are kg 2 velocity 2 meter = 1 sec 2 We will interpret this velocity in upcoming slides.
II. Generalizing the 1-dimensional wave equation First generalize the notation. i) "q" has meant transverse deflection of the string. Replace q Ψ, where Ψ may indicate other properties of the medium that
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
More information