Wave Equations Explicit Formulas In this lecture we derive the representation formulas for the wave equation in the whole space:
|
|
- Terence McDaniel
- 6 years ago
- Views:
Transcription
1 Nov. 07 Wave Equations Explicit Formulas In this lecture we derive the representation formulas for the wave equation in the whole space: u u t t u = 0, R n 0, ) ; u x, 0) = g x), u t x, 0) = h x). ) It turns out that the properties of the solutions depend on the dimension. More specifically, there are three cases: n =, n > even; n > odd. We will discuss in detail the three representative cases: n =,, 3 the order is actually n =, 3,, for reasons that will be clear soon).. n =. We consider the D wave equation u t t u x x = 0, R 0, ) ; u x, 0) = g x), u t x, 0) = h x). ) This equation can be solved via the following change of variables: and to make things clearer we set ũ ξ, η) = u x, t). With this change of variable we compute Therefore ξ = ; η =, 3) u t = ũ ξ ξ t + ũ η η t = ũ ξ ũ η 4) u t t = ũ ξξ ũ ξη + ũ ηη 5) u x = ũ ξ + ũ η 6) u x x = ũ ξξ + ũ ξη + ũ ηη 7) u t t u x x = 0 ũ ξη = 0 ũ ξ, η) = φ ξ) + ψ η) φ ) + ψ ). 8) Now using the initial values we have which yields This is d Alembert s formula. φ x) + ψ x) = g x) ; φ x) ψ x) = h x) 9) [ g ) + g ) ] + Theorem. Assume g C R), h C R), define u by Then i. u C R [ 0, ) ) ; ii. u t t u x x = 0 in R 0, ) ; [ g ) + g ) ] + iii. u takes the correct boundary values: x, t) x 0, 0) x, t) x 0, 0) h y) dy. 0) h y) dy. ) g x 0 ) ; ) u t x, t) = h x 0 ). 3) Proof. The proof is by direct calculation and is left as an exercise. Remark. It is easy to generalize the above theorem to the case g C k, h C k u C k. 4)
2 But in general u cannot be smoother. For example, consider the case h = g, then g ). It is clear that u cannot have better regularity than g. Remark 3. As we have shown in the exercises after the lecture on distributions, the formula remains true even for distributional solutions of the D wave equation. φ ) + ψ ) 5). Spherical means and Euler-Poisson-Darboux equation. The case n is much more complicated. The idea is the reduce the wave equation to a D equation which can be solved explicitly. The reduction is fulfilled through introducing the following auxiliary functions. Let u = u x, t). We define at each x R n, U x; r, t) u w, t) ds w, 6) B r x ) G x; r) g w) ds w, 7) B r x ) H x; r) h w) ds w. 8) B r x ) Note that when u is continuous, we can recover u from U by taking r 0. It turns out that U x; r, t) as a function of r and t satisfies a D equation which can be further transformed in to a D wave equation in the first quadrant) and enjoys the same regularity as u. Lemma 4. Euler-Poisson-Darboux equation) Fix x R n. Let u x, t) C m, m solves the wave equation. Then U x; r, t) u w, t) ds w 9) belongs to C m R + [ 0, ) ), and satisfies B r x) U t t U r r n r U r = 0 R + 0, ) ; U r, 0) = G r), U t r, 0) = H r). 0) Remark 5. Notice that r r n r r is just in R n with radial symmetry. Proof. Recall that U r x; r, t) = r n B r x) B r x) y u y, t) dy = n α n) r n B r x ) y u y, t) dy ) This shows U C, and we can define U r x; 0, t) = 0. Differentiating w. r. t r again, [ U r r x; r, t) = d dr n α n) r n y u y, t) dy B r x ) = n y u + n B r B r x ) This shows U C and also U r r x; 0, t) can be defined. We further have U t t x; r, t) = u t t = using the equation. B r x ) ] B r x ) B r x) y u. ) y u 3)
3 3. n = 3, Kirchhoff s formula. Let U, G, H be the spherical means. We set Some calculation yields Ũ = r U, G = r G, H = r H. 4) Ũ t t Ũ r r = 0 R + 0, ) ; Ũ = G, Ũ t = H. 5) Remark 6. Note that here we used the fact that n = 3. Thus we need to solve the wave equation in the first quadrant. Example 7. Consider the wave equation in the first quadrant: Let u t t u x x = 0, x > 0, ; u x, 0) = g x), u t x, 0) = h x), u = 0 for x = 0,. 6) ũ x, t) = { u x, t) x > 0 u x, t) x < 0 7) and define similarly g, h. Then it is clear that ũ solves the wave equation with initial values g, h. Thus we have ũ x, t) = [ g ) + g ) ] + h y) dy. 8) Therefore the solution to the original problem is [ g ) + g ) ] + [ g ) g t x) ] + t+ x t x h y) dy x t 0 h y) dy t x 0 Now for our purpose, we only need the case t r remember that finally we will let r 0 and recover u from U). In this case Ũ x; r, t) = [ ] G r + t) G t r) + t+ r H y) dy. 30) We have Further computation yields Ũ x; r, t) r 0 r = G t) + H t) = t which is the Kirchhoff s formula. t B t x ) B t x ) g w) ds w ) + t t r B t x) 9) h w) ds w. 3 ) [ t h w) + g w) + g w) w x) ] ds w 3) 4. n=, Method of descent and Poisson s formula. It is not possible to simplify as we did in the n = 3 case. Instead, we use the so-called method of descent, which treats the solution u x, t) of the D wave equation as a solution to the 3D equation. We set and define ḡ, h similarly. ū x, x, x 3, t) u x, x, t). 33)
4 Using the Kirchhoff s formula we have ū x, t) = B t x ) B t x ) where x = x, x 3 ) and B t x ) is the ball in R 3. From definitions of the variaous bar-ed functions, we have 4 π t where w = ) y, ± t y. B t x ) t h w ) + ḡ w ) + x ḡ w ) w x ) ds w. 34) t h y) + g y) + y g y) y x) ds w 35) Finally, let D t x) denote the ball in R centered at x with radius t, we have t h y) + g y) + y g y) y x) 4 π t ) D t x ) y x dy / t This is the Poisson s formula. = D t D t x ) t g y) + t h y) + t g y x) t y x ) / dy. 36) Remark 8. Huygens Principle) We notice that the behavior of the solutions for the D and 3D wave equations are drastically different. In D, u x, t) depends on initial data in the whole ball D t x) while in 3D it only depends on the data on the boundary of the ball B t x). Or equivalently, in 3D the effect of a vibration is only felt at the front of its propagation while in D it is felt forever after the front passed. This is the so-called Huygens principle. Remark 9. For general n, we define Ũ r, t) = ) k r r r k U x; r, t) ) 37) and define G, H accordingly. Some calculation yields the solution ) ) n 3 ) ) n 3 t n g ds + t n γ n t t t t t B t x ) B t x ) h ds ) 38) for n odd, where γ n = 3 ) t t t γ n ) n n ). Then the method of descent yields t n B t B t ) g y) dy t y x ) + / t t ) n for n even, where γ n = 4 n ) n. See L. C. Evans Partial Differential Equations pp for details. t n B t Remark 0. Nonhomogeneous problem) For the nonhomogeneous problem we use the Duhamel s principle, obtaining B t ) h y) dy t y x ). / u = f, u = 0, u t = 0, 39) 0 t u x, t; s) ds 40). If u x, t) also depends on data in the whole ball B t x) in 3D, we would not be able to clearly hear anything!
5 where u x, t; s) solves In particular, we have n = : n = 3: u t t u = 0, u x, s; s) = 0, u t x, s; s) = f, s). 4 ) 4 π 0 t x s B t x ) x + s f y, t s) dy ds. 4) f y, t y x ) y x dy. 43) Exercises. Here the integrand is called the retarded potential. Exercise. Prove the following theorem for D wave equations: Theorem. A ssume g C R), h C R), define u b y u x, t) = [ g ) + g ) ] + h y) dy. 44) T hen i. u C R [ 0, ) ) ; ii. u t t u x x = 0 in R 0, ) ; iii. u takes the correct boundary values: u x, t) = g x 0) ; 45) x, t) x 0, 0) x, t) x 0, 0) u t x, t) = h x 0). 46) Exercise. Equipartition of energy) Let u C R [ 0, ) ) solve the initial-value problem for the D wave equation u t t u x x = 0, u = g, u t = h; 47) Suppose g, h have compact support. Let k t) u t x, t) dx, p t) u x x, t) dx 48) be the kinetic and potential energy. P rove that for t large enough, k t) = p t) = constant. Hint: Use d Alembert s formula to compute the energy directly. )
Wave Equations: Explicit Formulas In this lecture we derive the representation formulas for the wave equation in the whole space:
Math 57 Fall 009 Lecture 7 Sep. 8, 009) Wave Equations: Explicit Formulas In this lecture we derive the representation formulas for the wave equation in the whole space: u u t t u = 0, R n 0, ) ; u x,
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 informationLecture 7: The wave equation: higher dimensional case
Lecture 7: The wave equation: higher dimensional case Some auxiliary facts The co area formula. Let be a bounded domain and be a function of class. Denote by : the level set of. Let be any integrable function.
More informationu( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)
M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation
More informationwhere F denoting the force acting on V through V, ν is the unit outnormal on V. Newton s law says (assume the mass is 1) that
Chapter 5 The Wave Equation In this chapter we investigate the wave equation 5.) u tt u = and the nonhomogeneous wave equation 5.) u tt u = fx, t) subject to appropriate initial and boundary conditions.
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationn 2 xi = x i. x i 2. r r ; i r 2 + V ( r) V ( r) = 0 r > 0. ( 1 1 ) a r n 1 ( 1 2) V( r) = b ln r + c n = 2 b r n 2 + c n 3 ( 1 3)
Sep. 7 The L aplace/ P oisson Equations: Explicit Formulas In this lecture we study the properties of the Laplace equation and the Poisson equation with Dirichlet boundary conditions through explicit representations
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 informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
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 informationMA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20
MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20 MA 201 (2016), PDE 2 / 20 Vibrating string and the wave equation Consider a stretched string of length
More 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 informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
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 informationMTH 847: PDE I (Fall 2017) Exam 2,
MTH 847: PDE I (Fall 2017) Exam 2, 2017.12.8 Name: Standard exam rules apply: You are not allowed to give or receive help from other students. All electronic devices must be turned off for the duration
More informationMATH 6337: Homework 8 Solutions
6.1. MATH 6337: Homework 8 Solutions (a) Let be a measurable subset of 2 such that for almost every x, {y : (x, y) } has -measure zero. Show that has measure zero and that for almost every y, {x : (x,
More informationAPPLIED PARTIAL DIFFERENTIAL EQUATIONS
APPLIED PARTIAL DIFFERENTIAL EQUATIONS AN I N T R O D U C T I O N ALAN JEFFREY University of Newcastle-upon-Tyne ACADEMIC PRESS An imprint of Elsevier Science Amsterdam Boston London New York Oxford Paris
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 informationLecture Notes Math 632, PDE Spring Semester Sigmund Selberg Visiting Assistant Professor Johns Hopkins University
Lecture Notes Math 63, PDE Spring Semester 1 Sigmund Selberg Visiting Assistant Professor Johns Hopkins University CHAPTER 1 The basics We consider the equation 1.1. The wave equation on R 1+n u =, where
More informationLast Update: April 7, 201 0
M ath E S W inter Last Update: April 7, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections.. 5 in the textbook.
More informationTravel Time Tomography and Tensor Tomography, I
Travel Time Tomography and Tensor Tomography, I Plamen Stefanov Purdue University Mini Course, MSRI 2009 Plamen Stefanov (Purdue University ) Travel Time Tomography and Tensor Tomography, I 1 / 17 Alternative
More information6.3 Fundamental solutions in d
6.3. Fundamental solutions in d 5 6.3 Fundamental solutions in d Since Lapalce equation is invariant under translations, and rotations (see Exercise 6.4), we look for solutions to Laplace equation having
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 informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationThe second-order 1D wave equation
C The second-order D wave equation C. Homogeneous wave equation with constant speed The simplest form of the second-order wave equation is given by: x 2 = Like the first-order wave equation, it responds
More informationMATH325 - QUANTUM MECHANICS - SOLUTION SHEET 11
MATH35 - QUANTUM MECHANICS - SOLUTION SHEET. The Hamiltonian for a particle of mass m moving in three dimensions under the influence of a three-dimensional harmonic oscillator potential is Ĥ = h m + mω
More informationYAO LIU. f(x t) + f(x + t)
DUNKL WAVE EQUATION & REPRESENTATION THEORY OF SL() YAO LIU The classical wave equation u tt = u in n + 1 dimensions, and its various modifications, have been studied for centuries, and one of the most
More informationThe wave front set of a distribution
The wave front set of a distribution The Fourier transform of a smooth compactly supported function u(x) decays faster than any negative power of the dual variable ξ; that is for every number N there exists
More informationFOURIER TRANSFORMS OF SURFACE MEASURE ON THE SPHERE MATH 565, FALL 2017
FOURIER TRANSFORMS OF SURFACE MEASURE ON THE SPHERE MATH 565, FALL 17 1. Fourier transform of surface measure on the sphere Recall that the distribution u on R n defined as u, ψ = ψ(x) dσ(x) S n 1 is compactly
More informationFOURIER METHODS AND DISTRIBUTIONS: SOLUTIONS
Centre for Mathematical Sciences Mathematics, Faculty of Science FOURIER METHODS AND DISTRIBUTIONS: SOLUTIONS. We make the Ansatz u(x, y) = ϕ(x)ψ(y) and look for a solution which satisfies the boundary
More informationR. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant
R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant CONTENTS I. Introductory Remarks S1. General Information about the Variety of Solutions.
More informationThis is example 3 on page 44 of BGH and example (b) on page 66 of Troutman.
Chapter 4 The brachistochrone This is example 3 on page 44 of BGH and example (b) on page 66 of Troutman. We seek the shape of a frictionless wire starting at the origin and ending at some point (, d)
More informationSummer 2017 MATH Solution to Exercise 5
Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationM4P52 Manifolds, 2016 Problem Sheet 1
Problem Sheet. Let X and Y be n-dimensional topological manifolds. Prove that the disjoint union X Y is an n-dimensional topological manifold. Is S S 2 a topological manifold? 2. Recall that that the discrete
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
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 informationCoordinate Systems and Canonical Forms
Appendix D Coordinate Systems and Canonical Forms D.1. Local Coordinates Let O be an open set in R n. We say that an n-tuple of smooth realvalued functions defined in O, (φ 1,...,φ n ), forms a local coordinate
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 informationDIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17
DIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17 6. Geodesics A parametrized line γ : [a, b] R n in R n is straight (and the parametrization is uniform) if the vector γ (t) does not depend on t. Thus,
More informationRegularity for Poisson Equation
Regularity for Poisson Equation OcMountain Daylight Time. 4, 20 Intuitively, the solution u to the Poisson equation u= f () should have better regularity than the right hand side f. In particular one expects
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES The functions that we have met so far can be described by expressing one variable explicitly in terms of another variable. y For example,, or y = x sin x,
More informationWeek 01 : Introduction. A usually formal statement of the equality or equivalence of mathematical or logical expressions
1. What are partial differential equations. An equation: Week 01 : Introduction Marriam-Webster Online: A usually formal statement of the equality or equivalence of mathematical or logical expressions
More informationTopics in Stochastic Geometry. Lecture 4 The Boolean model
Institut für Stochastik Karlsruher Institut für Technologie Topics in Stochastic Geometry Lecture 4 The Boolean model Lectures presented at the Department of Mathematical Sciences University of Bath May
More informationHomework 3 solutions Math 136 Gyu Eun Lee 2016 April 15. R = b a
Homework 3 solutions Math 136 Gyu Eun Lee 2016 April 15 A problem may have more than one valid method of solution. Here we present just one. Arbitrary functions are assumed to have whatever regularity
More informationIn this lecture we calculate moments and products of inertia of some simple geometric figures. (Refer Slide Time: 0:22)
Engineering Mechanics Professor Manoj K Harbola Department of Physics Indian Institute of Technology Kanpur Module 4 Lecture No 37 Properties of plane surfaces VIII: second moment and product of an area,
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 informationLecture 3. Econ August 12
Lecture 3 Econ 2001 2015 August 12 Lecture 3 Outline 1 Metric and Metric Spaces 2 Norm and Normed Spaces 3 Sequences and Subsequences 4 Convergence 5 Monotone and Bounded Sequences Announcements: - Friday
More informationCovariant electrodynamics
Lecture 9 Covariant electrodynamics WS2010/11: Introduction to Nuclear and Particle Physics 1 Consider Lorentz transformations pseudo-orthogonal transformations in 4-dimentional vector space (Minkowski
More informationf(x)e ikx dx. (19.1) ˆf(k)e ikx dk. (19.2)
9 Fourier transform 9 A first look at the Fourier transform In Math 66 you all studied the Laplace transform, which was used to turn an ordinary differential equation into an algebraic one There are a
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 informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
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 information12.6 Polar Coordinates
1.6 Polar Coordinates Consider Poisson s equation on the unit disk. ( 1 r u ) + 1 u r r r r = f(r, θ) θ with 0 r 1, and 0 θ π. We approximate the equation by ( 1 u i+1j u ij r i+ 1 r r i u ij u i 1j r
More informationTaylor and Laurent Series
Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x
More informationStrauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u
Strauss PDEs e: Setion 3.4 - Exerise 3 Page 1 of 13 Exerise 3 Solve u tt = u xx + os x, u(x, ) = sin x, u t (x, ) = 1 + x. Solution Solution by Operator Fatorization Bring u xx to the other side. Write
More informationBACKGROUND IN SYMPLECTIC GEOMETRY
BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationExploring Substitution
I. Introduction Exploring Substitution Math Fall 08 Lab We use the Fundamental Theorem of Calculus, Part to evaluate a definite integral. If f is continuous on [a, b] b and F is any antiderivative of f
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 information2.5 Compound Inequalities
Section.5 Compound Inequalities 89.5 Compound Inequalities S 1 Find the Intersection of Two Sets. Solve Compound Inequalities Containing and. Find the Union of Two Sets. 4 Solve Compound Inequalities Containing
More informationMath 575-Lecture 26. KdV equation. Derivation of KdV
Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not
More informationMATH 52 FINAL EXAM SOLUTIONS
MAH 5 FINAL EXAM OLUION. (a) ketch the region R of integration in the following double integral. x xe y5 dy dx R = {(x, y) x, x y }. (b) Express the region R as an x-simple region. R = {(x, y) y, x y }
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 informationWave equation on manifolds and finite speed of propagation
Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we
More informationOn the Uniqueness of Weak Solutions to the 2D Euler Equations
On the Uniqueness of Weak Solutions to the 2D Euler Equations (properties of the flow) 2009 SIAM Conference on Analysis of PDEs University of California, Riverside 9 December 2009 Bounded vorticity and
More informationANALYSIS QUALIFYING EXAM FALL 2017: SOLUTIONS. 1 cos(nx) lim. n 2 x 2. g n (x) = 1 cos(nx) n 2 x 2. x 2.
ANALYSIS QUALIFYING EXAM FALL 27: SOLUTIONS Problem. Determine, with justification, the it cos(nx) n 2 x 2 dx. Solution. For an integer n >, define g n : (, ) R by Also define g : (, ) R by g(x) = g n
More informationThe Principle of Least Action
The Principle of Least Action In their never-ending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made
More informationLecture 45: The Eigenvalue Problem of L z and L 2 in Three Dimensions, ct d: Operator Method Date Revised: 2009/02/17 Date Given: 2009/02/11
Page 757 Lecture 45: The Eigenvalue Problem of L z and L 2 in Three Dimensions, ct d: Operator Method Date Revised: 2009/02/17 Date Given: 2009/02/11 The Eigenvector-Eigenvalue Problem of L z and L 2 Section
More informationExercise Set 4. D s n ds + + V. s dv = V. After using Stokes theorem, the surface integral becomes
Exercise Set Exercise - (a) Let s consider a test volum in the pellet. The substract enters the pellet by diffusion and some is created and disappears due to the chemical reaction. The two contribute to
More informationSupplement. The Extended Complex Plane
The Extended Complex Plane 1 Supplement. The Extended Complex Plane Note. In section I.6, The Extended Plane and Its Spherical Representation, we introduced the extended complex plane, C = C { }. We defined
More informationLecture 5 - Hausdorff and Gromov-Hausdorff Distance
Lecture 5 - Hausdorff and Gromov-Hausdorff Distance August 1, 2011 1 Definition and Basic Properties Given a metric space X, the set of closed sets of X supports a metric, the Hausdorff metric. If A is
More information5/2/2012 Miscellaneous exercises Orbital motion 1
5/2/22 Miscellaneous exercises Orbital motion Orbital motion This exercise explores the motion of a particle in a spherical Coulomb potential. The equations of motion are similar to those for a planet
More informationMultiple state optimal design problems with explicit solution
Multiple state optimal design problems with explicit solution Marko Vrdoljak Department of Mathematics Faculty of Science University of Zagreb, Croatia Krešimir Burazin Department of Mathematics University
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationCONVERGENCE OF EXTERIOR SOLUTIONS TO RADIAL CAUCHY SOLUTIONS FOR 2 t U c 2 U = 0
Electronic Journal of Differential Equations, Vol. 206 (206), No. 266, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu CONVERGENCE OF EXTERIOR SOLUTIONS TO RADIAL CAUCHY
More informationRepresentation of Functions as Power Series
Representation of Functions as Power Series Philippe B. Laval KSU Today Philippe B. Laval (KSU) Functions as Power Series Today / Introduction In this section and the next, we develop several techniques
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationLecturer: Bengt E W Nilsson
9 3 19 Lecturer: Bengt E W Nilsson Last time: Relativistic physics in any dimension. Light-cone coordinates, light-cone stuff. Extra dimensions compact extra dimensions (here we talked about fundamental
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationBoundary Layers: Stratified Fluids
Boundary Layers: Stratified Fluids Lecture 3 by Jeroen Hazewinkel continued from lecture 2 Using w I = /(2σS) 2 T I, the interior of the cylinder is be described by 2σS 2 T I + 1/2 4 2 This result can
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationHomology of a Cell Complex
M:01 Fall 06 J. Simon Homology of a Cell Complex A finite cell complex X is constructed one cell at a time, working up in dimension. Each time a cell is added, we can analyze the effect on homology, by
More informationLecture: Lorentz Invariant Dynamics
Chapter 5 Lecture: Lorentz Invariant Dynamics In the preceding chapter we introduced the Minkowski metric and covariance with respect to Lorentz transformations between inertial systems. This was shown
More informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMore On Carbon Monoxide
More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations of CO 1 / 26 More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations
More informationLecture 10 - Moment of Inertia
Lecture 10 - oment of Inertia A Puzzle... Question For any object, there are typically many ways to calculate the moment of inertia I = r 2 dm, usually by doing the integration by considering different
More informationCompression on the digital unit sphere
16th Conference on Applied Mathematics, Univ. of Central Oklahoma, Electronic Journal of Differential Equations, Conf. 07, 001, pp. 1 4. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu
More informationLagrangian for Central Potentials
Physics 411 Lecture 2 Lagrangian for Central Potentials Lecture 2 Physics 411 Classical Mechanics II August 29th 2007 Here we will review the Lagrange formulation in preparation for the study of the central
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationv( x) u( y) dy for any r > 0, B r ( x) Ω, or equivalently u( w) ds for any r > 0, B r ( x) Ω, or ( not really) equivalently if v exists, v 0.
Sep. 26 The Perron Method In this lecture we show that one can show existence of solutions using maximum principle alone.. The Perron method. Recall in the last lecture we have shown the existence of solutions
More informationMATH 312 Section 4.3: Homogeneous Linear Equations with Constant Coefficients
MATH 312 Section 4.3: Homogeneous Linear Equations with Constant Coefficients Prof. Jonathan Duncan Walla Walla College Spring Quarter, 2007 Outline 1 Getting Started 2 Second Order Equations Two Real
More information9 More on the 1D Heat Equation
9 More on the D Heat Equation 9. Heat equation on the line with sources: Duhamel s principle Theorem: Consider the Cauchy problem = D 2 u + F (x, t), on x t x 2 u(x, ) = f(x) for x < () where f
More 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 information7.1 Indefinite Integrals Calculus
7.1 Indefinite Integrals Calculus Learning Objectives A student will be able to: Find antiderivatives of functions. Represent antiderivatives. Interpret the constant of integration graphically. Solve differential
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationCorrections to the First Printing: Chapter 6. This list includes corrections and clarifications through November 6, 1999.
June 2,2 1 Corrections to the First Printing: Chapter 6 This list includes corrections and clarifications through November 6, 1999. We should have mentioned that our treatment of differential forms, especially
More information