Some Aspects of Solutions of Partial Differential Equations
|
|
- Ginger Little
- 5 years ago
- Views:
Transcription
1 Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum , Kerala Sakthivel@iist.ac.in Periyar University, Salem February 22, 2013 K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 1 / 24
2 Plan of the Talk 1 Formulation of Diffusion Model 2 Solution by Fourier Method 3 Nonlinear PDEs - Conservation Equations - Method of Characteristics 4 Weak Solutions for PDEs 5 Weak Solutions by Variational Methods K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 2 / 24
3 Formulation of Diffusion Model Consider A motionless liquid filling a straight tube or pipe A chemical substance, say die, which is diffusing through the liquid Let u(z, t) be the concentration of the die at position z of the pipe at time t. On the region x to x + x, the mass of the dye is x+ x x+ x M(t) = u(z, t)dz, and so dm u = (z, t)dz. x dt x t Fick s Law: Flux goes from region of higher concentration to the region of lower concentration. The rate of motion is propositional to the concentration gradient. dm dt = x+ x x u (z, t)dz = Net Change of Concentration t = k u x u (x + x, t) k (x, t). x K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 3 / 24
4 Formulation of Diffusion Model... Now, x+ x x u t (z, t)dz = k[ u x where k proportionality constant. u (x + x, t) (x, t)], x Applying Mean Value Theorem (!!!) on x < ζ < x + x, we get u [ u t (ζ, t) = k u ] (x + x, t) x x (x, t) / x. Taking x 0, we have u t (x, t) = k 2 u (x, t). x 2 K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 4 / 24
5 Formulation of Diffusion Model... Now, x+ x x u t (z, t)dz = k[ u x where k proportionality constant. u (x + x, t) (x, t)], x Applying Mean Value Theorem (!!!) on x < ζ < x + x, we get u [ u t (ζ, t) = k u ] (x + x, t) x x (x, t) / x. Taking x 0, we have u t (x, t) = k 2 u (x, t). x 2 In the multidimensional case, we get u t = k u K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 4 / 24
6 Formulation of Diffusion Model... Now, x+ x x u t (z, t)dz = k[ u x where k proportionality constant. u (x + x, t) (x, t)], x Applying Mean Value Theorem (!!!) on x < ζ < x + x, we get u [ u t (ζ, t) = k u ] (x + x, t) x x (x, t) / x. Taking x 0, we have u t (x, t) = k 2 u (x, t). x 2 In the multidimensional case, we get u t = k u + f (x, t), where f is a source" or sink" of the dye. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 4 / 24
7 Fourier Transform and its Properties Let f L 1 (R). The Fourier transform of f is defined as f (ω) := F [f (x)] = 1 2π e iωx f (x)dx. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 5 / 24
8 Fourier Transform and its Properties Let f L 1 (R). The Fourier transform of f is defined as Inverse Fourier Transform: f (ω) := F [f (x)] = 1 2π f (x) := F 1 [ f (ω)] = 1 2π e iωx f (x)dx. e iωx f (ω)dω. Some Properties: If f is a continuous piecewise differentiable function with f, f x, f xx L 1 (R) and f (x) = 0, then lim x F [f x (x)] = iωf [f (x)] and F [f xx (x)] = ω 2 F [f (x)]!!!. Fourier transform is linear, that is, f, g, L 1 (R) and a, b R, then F [af + bg] = af [f ] + bf [g]. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 5 / 24
9 Fourier Transform and its Properties... Note that F [f (x)g(x)] F [f (x)]f [g(x)], whereas F [(f g)(x)] = F [f (x)]f [g(x)], where (f g) is the convolution" of f and g defined as (f g)(x) := 1 2π f (x y)g(y)dy = (g f )(x). K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 6 / 24
10 Fourier Transform and its Properties... Note that F [f (x)g(x)] F [f (x)]f [g(x)], whereas F [(f g)(x)] = F [f (x)]f [g(x)], where (f g) is the convolution" of f and g defined as (f g)(x) := 1 2π 1-D Heat Conduction Model (in an infinite rod): f (x y)g(y)dy = (g f )(x). u t (x, t) = ku xx (x, t), < x <, 0 < t < u(x, 0) = φ(x), < x <. Applying Fourier transform F [u t (x, t)] = kf [u xx (x, t)], F [u(x, 0)] = F [φ(x)] K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 6 / 24
11 1-D Heat Conduction Model Recall that F [u t (x, t)] = = t 1 2π e [ 1 2π iξx u (x, t)dx t ] e iξx u(x, t)dx Taking φ(ξ) = F [φ(x)], we get the following ODE": dû dt (t) = kξ2 û(t) with data û(0) = φ(ξ). = û (ξ, t). t K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 7 / 24
12 1-D Heat Conduction Model Recall that F [u t (x, t)] = = t 1 2π e [ 1 2π iξx u (x, t)dx t ] e iξx u(x, t)dx Taking φ(ξ) = F [φ(x)], we get the following ODE": dû dt (t) = kξ2 û(t) with data û(0) = φ(ξ). = û (ξ, t). t Solution to the transformed ODE: F [u(x, t)] = û(t) = φ(ξ)e ktξ2. Note that: (f g)(x) = F 1[ F [f (x)]f [g(x)] ] and F [e a2 x 2 ] = 1 2a e ξ2 /4a 2. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 7 / 24
13 Solution of 1-D Heat Conduction Model Taking inverse Fourier transform u(x, t) = F 1 [ φ(ξ)e ktξ2 ] = 1 F 1[ ] F [φ(x)]f [e x 2 /4kt ] 2kt = φ(x) [ 1 ] e x 2 /4kt 2kt = 1 4πkt φ(y)e (x y)2 /4kt dy. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 8 / 24
14 Solution of 1-D Heat Conduction Model Taking inverse Fourier transform u(x, t) = F 1 [ φ(ξ)e ktξ2 ] = 1 F 1[ ] F [φ(x)]f [e x 2 /4kt ] 2kt = φ(x) = 1 4πkt [ 1 ] e x 2 /4kt 2kt φ(y)e (x y)2 /4kt dy. Fundamental Solution (or Heat Kernel) of 1-D Heat Equation: Multidimensional case: Φ(x, t) = 1 4πt e x 2 /4t, k = 1, t > 0. 1 Ψ(x, t) = /4t (4πt) d/2 e x 2, x R d, t > 0, where x = (x 1, x 2,, x d ) and x = x x x 2 d. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 8 / 24
15 Properties of Fundamental Solution For t > 0, the function Ψ(x, t) > 0 is an infinitely differentiable function of x and t. (Recall Smoothing Property) Ψ t (x, t) = Ψ(x, t), x R d, t > 0 and Ψ(x, t)dx = 1, t > 0. R d For any continuous and bounded function g(x) : lim Ψ(x, t)g(x)dx = g(0). t 0 R d K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 9 / 24
16 Properties of Fundamental Solution For t > 0, the function Ψ(x, t) > 0 is an infinitely differentiable function of x and t. (Recall Smoothing Property) Ψ t (x, t) = Ψ(x, t), x R d, t > 0 and Ψ(x, t)dx = 1, t > 0. R d For any continuous and bounded function g(x) : lim Ψ(x, t)g(x)dx = g(0). t 0 R d Since Ψ(x, t) is the probability density of a Gaussian random variable, Ψ(x, t)dx = 1, for all t > 0!!!. R d Moreover, by Lebesgue dominated convergence theorem Ψ(x, t)g(x)dx R d 1 = (4πt) d/2 e x 2 /4t g(x)dx = 1 R d π d/2 e y 2 g(y 4t)dx R d 1 = g(0)dx = g(0) as t 0. (π) d/2 e y 2 R d K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 9 / 24
17 Conservation Equations - Traffic Flow Model Consider a stretch of highway on which cars are moving from left to right. Assume that no exit or entrance of ramps. u(x, t) density of cars at x f (x, t) flux of cars at x (cars per minute passing the point x) Change in the number of cars in [a, b] = d dt b Using flux, change in the number of cars in [a, b] = = f (a, t) f (b, t) = b a a u(x, t)dx f (x, t)dx, x by fundamental theorem of calculus. The last two integrals yield b a u (x, t)dx = t The interval length [a, b] is arbitrary, b a f (x, t)dx. x u t (x, t) + f x (x, t) = 0. (Conservation Equation) K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 10 / 24
18 Traffic Flow Model The amount of cars passing a given point is generally a function of density u, that is, f (u) : for example, f (u) = u 2 (quadratic flow rate). From conservation equation and chain rule: u t + 2uu x = 0. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 11 / 24
19 Traffic Flow Model The amount of cars passing a given point is generally a function of density u, that is, f (u) : for example, f (u) = u 2 (quadratic flow rate). From conservation equation and chain rule: Consider the model u t + 2uu x = 0. u t + 2uu x = 0, < x <, 0 < t < with the initial density of cars!!! 1 x 0 u(x, 0) = u 0 (x) = 1 x 0 < x < 1 0 x 1 K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 11 / 24
20 Solution by Method of Characteristics Let x(t) be a curve in x t plane and u(x(t), t) is the function of x and t. How does u change from the perspective of x(t)? By chain rule du dx = u x dt dt + u t; but taking traffic flow model into account So dx dt u t + 2uu x = 0 { dx dt = 2u du dt = 0 = 2u x = 2u(x, t)t + x 0 Characteristic equation Along such characteristic u is constant. From the initial density of cars, the characteristic curve starting from (x 0, 0) is x = 2u(x 0, 0)t + x 0 K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 12 / 24
21 Solution by Method of Characteristics... Case 1: For x 0 0, u(x 0, 0) = 1, the characteristics are x = 2t + x 0 or t = 1 2 (x x 0). Case 2: For 0 < x 0 < 1 u(x 0, 0) = 1 x 0, the characteristics are x = 2(1 x 0 )t + x 0 or t = 1 x x x 0 Case 3: For 1 x 0 <, u(x 0, 0) = 0, the characteristics are vertical lines given by x = x 0. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 13 / 24
22 Solution by Method of Characteristics... Case 1: For x 0 0, u(x 0, 0) = 1, the characteristics are x = 2t + x 0 or t = 1 2 (x x 0). Case 2: For 0 < x 0 < 1 u(x 0, 0) = 1 x 0, the characteristics are x = 2(1 x 0 )t + x 0 or t = 1 x x x 0 Case 3: For 1 x 0 <, u(x 0, 0) = 0, the characteristics are vertical lines given by x = x 0. Observations: Note that characteristics run together starting at t = 1/2. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 13 / 24
23 Solution by Method of Characteristics... Case 1: For x 0 0, u(x 0, 0) = 1, the characteristics are x = 2t + x 0 or t = 1 2 (x x 0). Case 2: For 0 < x 0 < 1 u(x 0, 0) = 1 x 0, the characteristics are x = 2(1 x 0 )t + x 0 or t = 1 x x x 0 Case 3: For 1 x 0 <, u(x 0, 0) = 0, the characteristics are vertical lines given by x = x 0. Observations: Note that characteristics run together starting at t = 1/2. So except for short time t = 1/2, the solution is not continuous. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 13 / 24
24 Solution by Method of Characteristics... Case 1: For x 0 0, u(x 0, 0) = 1, the characteristics are x = 2t + x 0 or t = 1 2 (x x 0). Case 2: For 0 < x 0 < 1 u(x 0, 0) = 1 x 0, the characteristics are x = 2(1 x 0 )t + x 0 or t = 1 x x x 0 Case 3: For 1 x 0 <, u(x 0, 0) = 0, the characteristics are vertical lines given by x = x 0. Observations: Note that characteristics run together starting at t = 1/2. So except for short time t = 1/2, the solution is not continuous. When characteristics run together, we have a phenomenon of Shock Waves" (discontinuous solutions) K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 13 / 24
25 Solution by Method of Characteristics... Case 1: For x 0 0, u(x 0, 0) = 1, the characteristics are x = 2t + x 0 or t = 1 2 (x x 0). Case 2: For 0 < x 0 < 1 u(x 0, 0) = 1 x 0, the characteristics are x = 2(1 x 0 )t + x 0 or t = 1 x x x 0 Case 3: For 1 x 0 <, u(x 0, 0) = 0, the characteristics are vertical lines given by x = x 0. Observations: Note that characteristics run together starting at t = 1/2. So except for short time t = 1/2, the solution is not continuous. When characteristics run together, we have a phenomenon of Shock Waves" (discontinuous solutions) Shock waves occurs due to the flux that grows very large as a function of density u. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 13 / 24
26 Notion of Weak Solution for PDEs Test Functions : C 0 () Let R d be a bounded domain. The set of all functions beloging to C () and compactly supported in. Support of a function ϕ(x) : supp(ϕ) = {x : ϕ(x) 0}. Example (Compact Support) Suppose d = 1. { ϕ(x) = e 1 x 2 1 if x < 1. 0 elsewhere. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 14 / 24
27 Notion of Weak Solution for PDEs Test Functions : C 0 () Let R d be a bounded domain. The set of all functions beloging to C () and compactly supported in. Support of a function ϕ(x) : supp(ϕ) = {x : ϕ(x) 0}. Example (Compact Support) Suppose d = 1. { ϕ(x) = e 1 x 2 1 if x < 1. 0 elsewhere. Example (Weak Differentiability): Suppose d = 1 and = ( 1, 1). The function u(x) = x, x is not differentiable in the classical sense. But it has a weak derivative!!! { 1 if x < 0 Du = 1 if x > 0. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 14 / 24
28 Weak Differentiability... For any ϕ C 0 (), integrating by parts, we get: 1 1 = u(x)dϕ(x)dx 0 = = u(x)dϕ(x)dx u(x)dϕ(x)dx Du(x)ϕ(x)dx + uϕ 0 1 Du(x)ϕ(x)dx + uϕ 1 Du(x)ϕ(x)dx [u(0)]ϕ(0). Note that ϕ(1) = ϕ( 1) = 0 and [u(0)] = u(0 + ) u(0 ) = 0 since u is continuous at x = K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 15 / 24
29 Sobolev Space H 1 () H 1 () := {v L 2 (), v L 2 ()}. The space H 1 () is a separable Hilbert space equipped with inner product (u, v) H 1 () = uvdx + u vdx. Norm on H 1 (): u 2 H 1 () = Space H 1 0 () := {u H1 () : u = 0}. u 2 dx + u 2 dx. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 16 / 24
30 Weak Solution for Poisson Equation Let R d be a bounded domain with smooth boundary. } u(x) = f (x) in u(x) = 0 (1) where f (x) is the external force; x = (x 1, x 2,..., x d ) and is the Laplace operator = 2 x x xd 2. Weak and Classical solutions : Assume u C 2 () and ϕ C0 (). Multiplying Poisson equation by ϕ and integrating on, u(x)ϕ(x)dx = f (x)ϕ(x)dx. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 17 / 24
31 Weak and Classical solutions... Green s identity: u u(x)ϕ(x)dx = ν (x)ϕ(x)ds + u ϕdx where ν is the unit normal vector outward to. Further, u ϕdx = u ϕ ν ds u ϕ(x)dx We get the following identities : u ϕdx = f (x)ϕ(x)dx (2) u ϕ(x)dx = f (x)ϕ(x)dx (3) Thus, if u C 2 () is a solution of Poisson equation then for any ϕ C0 () the above identities hold. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 18 / 24
32 Weak and Classical solutions... Conversely for any ϕ C0 (), u C2 () satisfies (3), we also get ( u f )ϕ(x)dx = 0. Since ϕ is arbitrary, u f = 0 in, i.e, u is a solution of Poisson equation. Note : The integral (2) makes sense if u H 1 () whereas for (3) we only need u L 2 (). Weak Solution: A function u H 1 () is said to be a weak solution of (1), if the following identity holds: u ϕdx = f ϕdx for any ϕ C0 () or H1 0 () since C 0 () is dense in H1 0 (). K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 19 / 24
33 Optimization Result from Linear Algebra Consider a system Ax = b, where A is n n matrix, symmetric and positive definite and b R n. Recall that for x, y R n x y =< x, y >= x T y = n x i y i. i=1 A vector x R n is a solution of Ax = b iff it is global minimizer of φ(x) = 1 2 < Ax, x > < b, x >. Assume that φ has a global minimizer!!. For any y R n, ɛ R the variation Φ(ɛ) := φ(x + ɛy) achieves its minimum at ɛ = 0. Φ (0) = d dɛ φ(x + ɛy) ɛ=0 = 0. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 20 / 24
34 Optimization Result from Linear Algebra... Now Φ (ɛ) = d [ ] 1 < A(x + ɛy), x + ɛy > < b, x + ɛy > dɛ 2 = < A(x + ɛy), y > < b, y > But Φ (0) = 0 implying that < Ax, y > < b, y >=< Ax b, y >= 0. Since y is arbitrary, x is a solution of Ax = b. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 21 / 24
35 Variational Approach for Poisson Equation Consider a functional J[v] = 1 v 2 dx fvdx. 2 If u H0 1() is an extremal of the functional J[v] in H1 0 (), then u is a weak solution of the Dirichlet problem (1). Suppose u H0 1 () is an extremum of J[v](need a proof!!), then for any ϕ H0 1 () and ɛ R, F(ɛ) := J[u + ɛϕ] = 1 (u + ɛϕ) 2 dx f (u + ɛϕ)dx. 2 But F achieves its minimum at ɛ = 0, F (0) = dj dɛ (u + ɛϕ) ɛ=0 = 0. Note that F (ɛ) = ( u + ɛ ϕ) ϕdx f ϕdx So F (0) = 0 u ϕdx = f ϕdx. Therefore, K.Sakthivel u(iist) is a weak solution Solutions of Partial of the Differential Poisson Equationsequation. Periyar University, Salem 22 / 24
36 Uniqueness of Weak Solutions Suppose u 1 and u 2 are two weak solutions of the Poisson equation. For any ϕ H0 1 (), we have u 1 ϕdx = f ϕdx and u 2 ϕdx = f ϕdx Taking u := u 1 u 2, we arrive at u ϕdx = 0, Choosing ϕ = u, u 2 dx = 0. Applying the Poincare inequality: u 2 dx C() so that ϕ H 1 0 (). u 2 dx, u 2 dx = 0 u = 0 a.e. in. K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 23 / 24
37 Thank You! K.Sakthivel (IIST) Solutions of Partial Differential Equations Periyar University, Salem 24 / 24
Modeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationThe Heat Equation John K. Hunter February 15, The heat equation on a circle
The Heat Equation John K. Hunter February 15, 007 The heat equation on a circle We consider the diffusion of heat in an insulated circular ring. We let t [0, ) denote time and x T a spatial coordinate
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationMATH 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 informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationFOURIER TRANSFORMS. 1. Fourier series 1.1. The trigonometric system. The sequence of functions
FOURIER TRANSFORMS. Fourier series.. The trigonometric system. The sequence of functions, cos x, sin x,..., cos nx, sin nx,... is called the trigonometric system. These functions have period π. The trigonometric
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationPartial Differential Equations and Diffusion Processes
Partial Differential Equations and Diffusion Processes James Nolen 1 Department of Mathematics, Stanford University 1 Email: nolen@math.stanford.edu. Reproduction or distribution of these notes without
More informationHOMEWORK 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 informationMATH 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 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 informationDiffusion of a density in a static fluid
Diffusion of a density in a static fluid u(x, y, z, t), density (M/L 3 ) of a substance (dye). Diffusion: motion of particles from places where the density is higher to places where it is lower, due to
More informationFourier Transform & Sobolev Spaces
Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev
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 informationLecture 4: Fourier Transforms.
1 Definition. Lecture 4: Fourier Transforms. We now come to Fourier transforms, which we give in the form of a definition. First we define the spaces L 1 () and L 2 (). Definition 1.1 The space L 1 ()
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationMATH 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 informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationConservation and dissipation principles for PDEs
Conservation and dissipation principles for PDEs Modeling through conservation laws The notion of conservation - of number, energy, mass, momentum - is a fundamental principle that can be used to derive
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationSolution Sheet 3. Solution Consider. with the metric. We also define a subset. and thus for any x, y X 0
Solution Sheet Throughout this sheet denotes a domain of R n with sufficiently smooth boundary. 1. Let 1 p
More informationPartial 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 informationISSN: [Engida* et al., 6(4): April, 2017] Impact Factor: 4.116
IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY ELLIPTICITY, HYPOELLIPTICITY AND PARTIAL HYPOELLIPTICITY OF DIFFERENTIAL OPERATORS Temesgen Engida*, Dr.Vishwajeet S. Goswami
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 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 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 informationUNIVERSITY OF MANITOBA
DATE: May 8, 2015 Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationDiffusion - 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 informationSolutions of differential equations using transforms
Solutions of differential equations using transforms Process: Take transform of equation and boundary/initial conditions in one variable. Derivatives are turned into multiplication operators. Solve (hopefully
More informationLecture 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 informationScalar Conservation Laws and First Order Equations Introduction. Consider equations of the form. (1) u t + q(u) x =0, x R, t > 0.
Scalar Conservation Laws and First Order Equations Introduction. Consider equations of the form (1) u t + q(u) x =, x R, t >. In general, u = u(x, t) represents the density or the concentration of a physical
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 informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationIntroduction to Microlocal Analysis
Introduction to Microlocal Analysis First lecture: Basics Dorothea Bahns (Göttingen) Third Summer School on Dynamical Approaches in Spectral Geometry Microlocal Methods in Global Analysis University of
More informationPartial Differential Equations 2 Variational Methods
Partial Differential Equations 2 Variational Methods Martin Brokate Contents 1 Variational Methods: Some Basics 1 2 Sobolev Spaces: Definition 9 3 Elliptic Boundary Value Problems 2 4 Boundary Conditions,
More informationThe 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 informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationMATH FALL 2014
MATH 126 - FALL 2014 JASON MURPHY Abstract. These notes are meant to supplement the lectures for Math 126 (Introduction to PDE) in the Fall of 2014 at the University of California, Berkeley. Contents 1.
More informationThe Fourier Transform Method
The Fourier Transform Method R. C. Trinity University Partial Differential Equations April 22, 2014 Recall The Fourier transform The Fourier transform of a piecewise smooth f L 1 (R) is ˆf(ω) = F(f)(ω)
More informationGENERATORS WITH INTERIOR DEGENERACY ON SPACES OF L 2 TYPE
Electronic Journal of Differential Equations, Vol. 22 (22), No. 89, pp. 3. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GENERATORS WITH INTERIOR
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationMath 124A October 11, 2011
Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although
More 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 informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More 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 informationBoundary conditions. Diffusion 2: Boundary conditions, long time behavior
Boundary conditions In a domain Ω one has to add boundary conditions to the heat (or diffusion) equation: 1. u(x, t) = φ for x Ω. Temperature given at the boundary. Also density given at the boundary.
More informationFinal 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 information4 Riesz Kernels. Since the functions k i (ξ) = ξ i. are bounded functions it is clear that R
4 Riesz Kernels. A natural generalization of the Hilbert transform to higher dimension is mutiplication of the Fourier Transform by homogeneous functions of degree 0, the simplest ones being R i f(ξ) =
More information1 Fourier Integrals of finite measures.
18.103 Fall 2013 1 Fourier Integrals of finite measures. Denote the space of finite, positive, measures on by M + () = {µ : µ is a positive measure on ; µ() < } Proposition 1 For µ M + (), we define the
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationOn some weighted fractional porous media equations
On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationDiffusion Processes. Lectures INF2320 p. 1/72
Diffusion Processes Lectures INF2320 p. 1/72 Lectures INF2320 p. 2/72 Diffusion processes Examples of diffusion processes Heat conduction Heat moves from hot to cold places Diffusive (molecular) transport
More information1. 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 informationWeak form of Boundary Value Problems. Simulation Methods in Acoustics
Weak form of Boundary Value Problems Simulation Methods in Acoustics Note on finite dimensional description of functions Approximation: N f (x) ˆf (x) = q j φ j (x) j=1 Residual function: r(x) = f (x)
More informationMATH 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 informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationNotions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy
Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationTraffic flow problems. u t + [uv(u)] x = 0. u 0 x > 1
The flow of cars is modelled by the PDE Traffic flow problems u t + [uvu)] x = 1. If vu) = 1 u and x < a) ux, ) = u x 2 x 1 u x > 1 b) ux, ) = u e x, where < u < 1, determine when and where a shock first
More informationPDEs, 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 informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
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 informationTutorial 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 informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
More information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
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 informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationDouble Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids
Double Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids joint work with Hengguang Li Yu Qiao School of Mathematics and Information Science Shaanxi Normal University Xi an,
More informationNOTES ON NEWLANDER-NIRENBERG THEOREM XU WANG
NOTES ON NEWLANDER-NIRENBERG THEOREM XU WANG Abstract. In this short note we shall recall the classical Newler-Nirenberg theorem its vector bundle version. We shall also recall an L 2 -Hörmer-proof given
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationANALYTIC SEMIGROUPS AND APPLICATIONS. 1. Introduction
ANALYTIC SEMIGROUPS AND APPLICATIONS KELLER VANDEBOGERT. Introduction Consider a Banach space X and let f : D X and u : G X, where D and G are real intervals. A is a bounded or unbounded linear operator
More informationMethods on Nonlinear Elliptic Equations
Wenxiong Chen and Congming Li Methods on Nonlinear Elliptic Equations SPIN AIMS internal project number, if known Monograph September 22, 2008 AIMS Preface In this book we intend to present basic materials
More information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More information