Math 6395: Hyperbolic Conservation Laws and Numerical Methods. Hyperbolic Conservation Laws are a type of PDEs arise from physics: fluid/gas dynamics.
|
|
- Tamsin Daniels
- 5 years ago
- Views:
Transcription
1 Mth 6395: Hyperbolic Conservtion Lws nd Numericl Methods Course structure Hyperbolic Conservtion Lws re type of PDEs rise from physics: fluid/gs dynmics. Distintive solution structures Other types of PDEs we met before, e.g. the time independent elliptic equtions, or time dependent prbolic equtions. Solutions re smooth. The solutions of hyperbolic problems hve distinctive feture: formtion of discontinuities (shocks). Such feture together with others fctors ttrct mny ttentions mong the scientific nd engineering community, leding to mny interesting reserch work in the PDE level, nd in the design of numericl schemes in the pst yers or so. Structure of the course: mthemticl theory numericl methods methods for liner equtions: locl trunction error, ccurcy, stbility, convergence methods for nonliner sclr equtions low order schemes with good properties corresponding to those in the PDE theory extension to high order schemes in the discontinuous Glerkin frmework. 1 Introduction Hyperbolic conservtion lws: u t + x f( u) = 0, (x, t) R d R + (1) u( x, t = 0) = u 0 ( x) (2) u: conserved quntities, R d R + R m, with d being the dimension of the problem, nd m being the number of components of u. f(u): flux function, R m R m For exmple The cse of m = 1 nd d = 1 corresponds to 1-D sclr conservtion lws. In trffic flow problems, ρ t + f(ρ) x = 0 is of the form of (1). ρ is the density of the trffic, nd f(ρ) = ρv(ρ) is the flux pss loction (mss multiple by the velocity of the trffic). m = 2 + d. In gs/fluid dynmics, the compressible Euler equtions in one dimension reds u t + x f( u) = 0, (3) where u = (ρ, ρv, E) T re the conserved quntities (mss, momentum nd energy) nd f( u) = (ρv, ρv 2 + p, (E + p)v) T re the corresponding flux functions. Here ρ is the density, v is the velocity, ρv is the momentum, E is the energy nd p is the pressure given s function of other stte vribles, known s eqution of stte. Hyperbolic equtions re type of time dependent liner/nonliner PDEs rise from science nd engineering. Most of the time, it is not possible to write down the exct solutions of the equtions. People hevily rely on numerics to ccess the solution structures in future times. 1
2 Definition 1.1. (Hyperbolic for 1-D problem) The 1-D system (1) is hyperbolic if the m m Jcobin mtrix f 1 f u 1 1 u m J f ( u) = f m u 1 of the flux function hs the following property: For ech vlue of u the eigenvlues of J f ( u) re rel, nd the mtrix is digonlizble, i.e. there is complete set of m linerly independent eigenvectors. f m u m Definition 1.2. (Hyperbolic for 2-D problem) The 2-D system in the form of u t + f( u) x + g( u) y = 0 (4) with u: R 2 R + R m, f, g : R m R m is hyperbolic if ny liner combintion of m m mtrix αj f ( u) + βj g ( u) of the flux Jcobins hs rel eigenvlues nd the mtrix is digonlizble, i.e. there is complete set of m linerly independent eigenvectors. Remrk 1.3. If the Jcobin mtrices hve distinct rel eigenvlues, it follows tht it s digonlizble. In this cse the system is clled strictly hyperbolic. If the eigenvlues re rel but not distinct, the Jcobin mtrix my not be digonlizble. Such system is clled wekly hyperbolic. If the Jcobin mtrix is digonlizble but hs complex eigenvlues, the system is not hyperbolic. Typicl fetures of solutions for nonliner hyperbolic systems. Development of discontinuities or shocks in future solutions, even with smooth initil dt. Becuse of this, solutions in clssicl sense could fil to exist. If the solutions re defined in the wek distribution sense, there could be more thn one wek solutions Entropy solution is the uniquely exist physiclly relevnt solution mong wek solutions. In the PDE level, questions to sk/nswer re: How to define entropies solutions? (vnishing viscosity method, entropy inequlities) Wht re pproprite spces (norms) for entropy solutions? (BV norm, L 1 nd L ) Wht re criteri for selecting entropy solutions mong wek solutions? (entropy conditions) Becuse of the bove mentioned difficulties t the PDE level, numericl chllenges in designing schemes for hyperbolic type PDEs re how to design schemes with good ccurcy nd stbility? ( good : high resolution/high order ccurcy, stbility in the sense of totl vritionl diminishing or mximum principle preserving) how to cpture shocks without numericl rtifcts such s oscilltions? how to know (or prove) if the numericl scheme pproximtes wek solutions, or the unique entropy solution mong wek solutions? 2 The Derivtion of Conservtion Lws We cn model the trffic flow s 1-D sclr conservtion lw. Denote density of the trffic s ρ(x, t), then the evolution of the trffic flow density cn be described by the following t ρ(x, t)dx = f(ρ(, t)) f(ρ(b, t)), (5) for ny b. Eqution (5): the rte of chnge of the density over n intervl [, b] is due to the flux flowing into the region from the left boundry x =, nd the flux flowing out of the region from the right boundry x = b. The flux pst given point is give by f(ρ(x, t)) = ρ(x, t)v(x, t), 2
3 where v(x, t) is the velocity of the trffic flow, which could be dependent on ρ. For exmple, in simple liner model, v(ρ) = v mx (1 ρ/ρ mx ). At zero density (empty rod), the velocity is v mx ; but decreses to zero s ρ pproches the mximum cpcity of the rod ρ mx (trffic jm). If we ssume f(ρ(x, t)) x exists, then Hence t ρ(x, t)dx = f(ρ(x, t)) x dx (6) ( t ρ(x, t) + f(ρ(x, t)) x)dx = 0 (7) Since the bove eqution is true for ny b, then in terms of the differentil eqution, we hve Remrk 2.1. Note tht ρ t + f(ρ) x = 0. (8) Eqution(5) smoothf(ρ) = Eqution(8); Eqution(8) Eqution(5) The integrl eqution does not require smoothness conditions on the solution. Remrk 2.2. System of conservtion lws such s Euler eqution in gs/fluid dynmics cn be modeled in similr (but more involved) wy concerning the conservtion of mss, momentum nd energy. See Chp. 5 of the book. 3 Sclr Conservtion Lws 3.1 Liner cse u t + u x = 0, u(x, t = 0) = u 0 (x), < x <, (9) chrcteristics in x t plne in deriving the exct solution. Along chrcteristics dx dt =, the solution stys constnt du dt = u t + u x x t = u t + u x = 0. (10) Hence the exct solution is u(x, t) = u 0 (x t). finite speed of propgtion (vs. infinite speed of propgtion for prbolic equtions) If u 0 (x) is smooth function, then u(x, t) is eqully smooth in spce nd in time. If u 0 (x) is not smooth (the concept of clssicl solutions for differentil equtions fils), the solution of u 0 (x t) stisfies the integrl form of the eqution, nd is solution in the wek/distribution sense. 3.2 Nonliner cse e.g. Burgers eqution f(u) = u 2 /2. u t + f(u) x = 0 (11) 3
4 3.2.1 Chrcteristics nd Shock Formtion Along chrcteristics dx dt = f (u) = u, the solution stys constnt du dt = u t + u x x t = u t + f (u)u x = 0. (12) Hence the exct solution u(x, t) = u 0 (x ), where x is the root of the chrcteristic eqution x x t 0 = f (u 0 (x )). (13) Exmple 3.1. Consider the Burgers eqution with initil condition u 0 (x) = sin(x). Plot for crossing of chrcteristics. Shock formtion time. Let d. dt = t + x x t = t + u x be the differentition long chrcteristics. Let v = u x, du dt = 0. dv dt = v t + f (u) v x = 2 u t x + u 2 u x 2. Tking xderivtive on the Burgers eqution we hve, 2 u t x = 2 f(u) x 2, Thus dv dt = f(u) 2 x 2 + u 2 u x 2 = (u x) 2 = v 2. 1 v(t ) 1 v(0) d( 1 v ) = dt v(t ) = = T v(0) T v(0) + 1 (14) v(t ) goes to infinity when T = 1 v(0). Shock forms when T = min x:u 0 (x)<0 { 1 u 0 (x)} Wek Solutions Let φ(x, t) C0(R 1 R + ) be test function, where C0 1 is the spce of function tht re continuously differentible with compct support. If we multiply the eqution (11) by the test function φ nd integrte over spce nd time, we obtin Performing integrtion by prts yields 0 0 (φ t u + φ x f(u))dxdt + φ(u t + f(u) x )dxdt = 0. (15) φ(x, 0)u(x, 0)dx = 0. (16) Definition 3.2. (Wek solution) The function u(x, t) is clled wek solution of the conservtion lws if eqution (16) holds for ll test function φ C 1 0(R R + ). Remrk 3.3. Mthemticlly the definition of wek solution in the integrl form bove is equivlent to the solution for the integrl form of the eqution over ny choices of sptil nd time intervl. The wek form bove could be esier to work with sometimes. 4
5 Specil cses. If u(x, t) is continuously differentible, then u(x, t) is clssicl solution nd wek solution. In the x t plne, if u(x, t) C 1 (R R + ) nd stisfies the strong/differentil form of eqution (11) except on curve x(t), then 0 = u(x, t)dx + f(u(b, t)) f(u(, t)) t ( = x(t) ) b u(x, t)dx + u(x, t)dx + f(u(b, t)) f(u(, t)) t x(t) + x(t) = u t (x, t)dx + u(x(t), t)x (t) + u t (x, t)dx u(x(t) +, t)x (t) + f(u(b, t)) f(u(, t)) x(t) + x(t) = u t (x, t)dx + f(u(x(t), t)) f(u(, t)) + u t (x, t)dx + f(u(b, t)) f(u(x(t) +, t)) x(t) + (f(u(x(t), t)) f(u(x(t) +, t))) + u(x(t), t)x (t) u(x(t) +, t)x (t) = (f(u(x(t), t)) f(u(x(t) +, t))) + u(x(t), t)x (t) u(x(t) +, t)x (t) (17) Let [ f(u) ]. = (f(u(x(t), t)) + f(u(x(t) +, t))) nd [ u ]. = (u(x(t), t) + u(x(t) +, t)) be the jump of f nd u cross discontinuity. From eqution (17), we hve s =. x (t) = [ f(u) ]. (18) [ u ] This is the so clled Rnkine-Hugoniot jump condition. s. = x (t) is the shock speed. Proposition 3.4. A function u(x, t) is piecewise smooth nd stisfies the PDE strong whenever u C 1. If the function stisfies the Rnkine-Hugoniot jump condition long the discontinuity curve, then u(x, t) is wek solution of eqution (11). Remrk 3.5. For curves on which the solution is continuous but not differentible, eqution (17) is stisfied. The solution is wek solution Riemnn Problem nd Not-uniqueness of Wek Solutions Exmple 3.6. Consider Burgers eqution with initil condition consist of two constnt sttes u l nd u r, u 0 (x) = { ul = 1, x < 0 u r = 1, x > 0 (19) Drwing chrcteristics: colliding into ech other (formtion of shocks). By Prop. 3.4, u(x, t) = u 0 (x) is the wek solution of the Burgers eqution. Exmple 3.7. Consider Burgers eqution with initil condition consist of two constnt sttes u l nd u r, u 0 (x) = { ul = 1, x < 0 u r = 1, x > 0 Drwing chrcteristics: spreding out. A wek solution is the rrefction wve 1, x < t x u(x, t) = t, t < x < t 1, x > t (20) (21) Yet by Prop. 3.4, u(x, t) = u 0 (x) is lso wek solution of the Burgers eqution. In fct there re infinitely mny wek solutions for this problem. 5
6 Remrk 3.8. (from wiki) A Riemnn problem, nmed fter Bernhrd Riemnn, consists of conservtion lw together with piecewise constnt dt hving single discontinuity. The Riemnn problem is very useful for the understnding of hyperbolic prtil differentil equtions like the Euler equtions becuse ll properties, such s shocks nd rrefction wves, pper s chrcteristics in the solution. It lso gives n exct solution to some complex nonliner equtions, such s the Euler equtions. In numericl nlysis, Riemnn problems pper in nturl wy in finite volume methods for the solution of eqution of conservtion lws due to the discreteness of the grid. For tht it is widely used in computtionl fluid dynmics nd in MHD simultions. In these fields Riemnn problems re clculted using Riemnn solvers Entropy Solution We hve seen tht the clssicl/strong solution my fil to exist for hyperbolic equtions; but wek solutions my not be unique. Additionl criteri re needed to pick out the unique physiclly relevnt solution, i.e. the entropy solution, mong wek solutions. In physicl models, the blnce lws only come with some physicl viscosity. For exmple, in the trffic flow modeling, the viscosity tkes the form of slow response of drivers nd utomobiles; in the fluid dynmics, the viscosity corresponds to the informl notion of thickness. For exmple, honey hs higher viscosity thn wter. Conservtion lws with viscous terms provide more physiclly relevnt models. Definition 3.9. (Vnishing viscosity method for the entropy solution) Consider the viscous eqution u ε t + f(u ε ) x = εu ε xx (22) An entropy solution of (11) is the limit (.e.) of u ε of eqution (22) when ε 0. There re severl issues to resolve to mke sense of the bove definition. For exmple, the solution for the viscous eqution (22) hs to exist nd is unique for ech ε; there hs to exist limit of u ε s ε 0; such limit should be one of the wek solutions (in fct the physiclly relevnt entropy solution) for eqution (11). These questions hve been nswered in the PDE theory, but they re out of the scope of this course. Proposition (Equivlent definition for entropy solution) A wek solution of (11) is n entropy solution if for ll convex entropy function U(u) with U (u) 0 nd the ssocited entropy flux function F (u) with F (u) = U (u)f (u), we hve U(u) t + F (u) x 0, (23) in the distribution sense. Tht is for ll φ 0, φ C 1 0(R R + ), we hve 0 Proof: see Section of the textbook. (φ t U(u) + φ x F (u))dxdt φ(x, 0)U(u(x, 0))dx 0. (24) Remrk If solution is continuously differentible (clssicl solution), then it is n entropy solution. The entropy inequlity becomes n equlity. If solution u(x, t) C 1 (R R + ) nd stisfies the strong/differentil form of eqution (11) except on curve x(t), then by using similr procedure in deriving the Rnkine-Hugoniot condition, we hve the following inequlity for the entropy solution s[ U ] + [ F (U) ] 0, for ny entropy-entropy flux pirs (25) The bove definitions for entropy solutions re mthemticlly rigorous; however, it will be difficult to directly pply these definitions to select the entropy solution mong wek solutions (e.g. for the exmple of Riemnn problem bove). Below, we look for simpler nd more prcticl entropy conditions tht one cn directly pply in exmples. Proposition (Oleinik entropy condition) A discontinuity propgting with speed s = [ f(u) ] [ u ] given by the Rnkine-Hugoniot jump condition stisfies the Oleinik entropy condition if for ll u between u l nd u r, f(u) f(u l ) u u l s f(u) f(u r) u u r (26) where u l nd u r re left nd right stte long the discontinuity respectively. 6
7 Proposition (Lx entropy condition) A discontinuity propgting with speed s given by the Rnkine- Hugoniot jump condition stisfies the Lx entropy condition if f (u l ) > s > f (u r ), (27) where u l nd u r re left nd right stte long the discontinuity respectively. Remrk One cn see tht Oleinik entropy condition implies Lx entropy condition; the converse does not hold. Lx entropy condition is necessry but not sufficient condition to single out the entropy condition. However, if f(u) is strictly convex with f (u) > 0 (or strictly concve with f (u) < 0), then the Lx entropy condition is equivlent to the Oleinik entropy condition, nd is sufficient to single out the entropy condition. In fct, under the convexity (or concvity) ssumption on f(u), the Lx entropy condition is reduced to f (u l ) > f (u r ), i.e., the chrcteristics propgting into the shock, rther thn diverging from the shock. Entropy solutions of the sclr hyperbolic eqution (11) re proved to exist nd re unique in the clss of BV L 1 L. Here BV stnds for bounded vrition. A function is of bounded vrition on given intervl [, b], if V b (f). = sup i f(x i ) f(x i+1 ) <, {x i } i is ny prtition of[, b]. If f(u) is continuously differentible, then V b (f). = f (x) dx. Proposition (BV) If u 0 (x) is function of loclly bounded vrition on (, ), then for ech t > 0, u(, t) is lso function of loclly bounded vrition on (, ), nd where s = mx x f (u). V R Ru(, t) T V R+st Rst u 0( ), Proposition (L 1 contrction property) If u(x, t) nd v(x, t) re solutions of the sclr hyperbolic eqution (11) with initil dt u 0 (x) nd v 0 (x) respectively, then u(, t) v(, t) L 1 u 0 ( ) v 0 ( ) L 1. Specificlly, consider v 0, then u(, t) L 1 u 0 ( ) L 1. Proposition (L mximum principle) If u(x, t) is solution of the sclr hyperbolic eqution (11) with initil condition u 0 (x), then mx x u(x, t) mx u 0(x), min x x u(x, t) min u 0 (x), x 7
Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004
Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More information1 2-D Second Order Equations: Separation of Variables
Chpter 12 PDEs in Rectngles 1 2-D Second Order Equtions: Seprtion of Vribles 1. A second order liner prtil differentil eqution in two vribles x nd y is A 2 u x + B 2 u 2 x y + C 2 u y + D u 2 x + E u +
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More information1.9 C 2 inner variations
46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More information1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation
1 1.1. Liner Constnt Coefficient Equtions Section Objective(s): Overview of Differentil Equtions. Liner Differentil Equtions. Solving Liner Differentil Equtions. The Initil Vlue Problem. 1.1.1. Overview
More informationKinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2)
Introduction Kinemtic Wves These re wves which result from the conservtion eqution E t + I = 0 (1) where E represents sclr density field nd I, its outer flux. The one-dimensionl form of (1) is E t + I
More informationConservation Law. Chapter Goal. 5.2 Theory
Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationAdvanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015
Advnced Clculus: MATH 410 Uniform Convergence of Functions Professor Dvid Levermore 11 December 2015 12. Sequences of Functions We now explore two notions of wht it mens for sequence of functions {f n
More informationPDE Notes. Paul Carnig. January ODE s vs PDE s 1
PDE Notes Pul Crnig Jnury 2014 Contents 1 ODE s vs PDE s 1 2 Section 1.2 Het diffusion Eqution 1 2.1 Fourier s w of Het Conduction............................. 2 2.2 Energy Conservtion.....................................
More informationHomework 4. (1) If f R[a, b], show that f 3 R[a, b]. If f + (x) = max{f(x), 0}, is f + R[a, b]? Justify your answer.
Homework 4 (1) If f R[, b], show tht f 3 R[, b]. If f + (x) = mx{f(x), 0}, is f + R[, b]? Justify your nswer. (2) Let f be continuous function on [, b] tht is strictly positive except finitely mny points
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationMapping the delta function and other Radon measures
Mpping the delt function nd other Rdon mesures Notes for Mth583A, Fll 2008 November 25, 2008 Rdon mesures Consider continuous function f on the rel line with sclr vlues. It is sid to hve bounded support
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationMA Handout 2: Notation and Background Concepts from Analysis
MA350059 Hndout 2: Nottion nd Bckground Concepts from Anlysis This hndout summrises some nottion we will use nd lso gives recp of some concepts from other units (MA20023: PDEs nd CM, MA20218: Anlysis 2A,
More informationMath 5440 Problem Set 3 Solutions
Mth 544 Mth 544 Problem Set 3 Solutions Aron Fogelson Fll, 25 1: Logn, 1.5 # 2) Repet the derivtion for the eqution of motion of vibrting string when, in ddition, the verticl motion is retrded by dmping
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationMath 5440 Problem Set 3 Solutions
Mth 544 Mth 544 Problem Set 3 Solutions Aron Fogelson Fll, 213 1: (Logn, 1.5 # 2) Repet the derivtion for the eqution of motion of vibrting string when, in ddition, the verticl motion is retrded by dmping
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationLecture 1: Introduction to integration theory and bounded variation
Lecture 1: Introduction to integrtion theory nd bounded vrition Wht is this course bout? Integrtion theory. The first question you might hve is why there is nything you need to lern bout integrtion. You
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More informationIntroduction to the Calculus of Variations
Introduction to the Clculus of Vritions Jim Fischer Mrch 20, 1999 Abstrct This is self-contined pper which introduces fundmentl problem in the clculus of vritions, the problem of finding extreme vlues
More information1 The fundamental theorems of calculus.
The fundmentl theorems of clculus. The fundmentl theorems of clculus. Evluting definite integrls. The indefinite integrl- new nme for nti-derivtive. Differentiting integrls. Tody we provide the connection
More informationLecture 3 ( ) (translated and slightly adapted from lecture notes by Martin Klazar)
Lecture 3 (5.3.2018) (trnslted nd slightly dpted from lecture notes by Mrtin Klzr) Riemnn integrl Now we define precisely the concept of the re, in prticulr, the re of figure U(, b, f) under the grph of
More informationDefinite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationThe final exam will take place on Friday May 11th from 8am 11am in Evans room 60.
Mth 104: finl informtion The finl exm will tke plce on Fridy My 11th from 8m 11m in Evns room 60. The exm will cover ll prts of the course with equl weighting. It will cover Chpters 1 5, 7 15, 17 21, 23
More informationSTURM-LIOUVILLE BOUNDARY VALUE PROBLEMS
STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS Throughout, we let [, b] be bounded intervl in R. C 2 ([, b]) denotes the spce of functions with derivtives of second order continuous up to the endpoints. Cc 2
More informationMA 124 January 18, Derivatives are. Integrals are.
MA 124 Jnury 18, 2018 Prof PB s one-minute introduction to clculus Derivtives re. Integrls re. In Clculus 1, we lern limits, derivtives, some pplictions of derivtives, indefinite integrls, definite integrls,
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationCalculus of Variations
Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function
More informationReview on Integration (Secs ) Review: Sec Origins of Calculus. Riemann Sums. New functions from old ones.
Mth 20B Integrl Clculus Lecture Review on Integrtion (Secs. 5. - 5.3) Remrks on the course. Slide Review: Sec. 5.-5.3 Origins of Clculus. Riemnn Sums. New functions from old ones. A mthemticl description
More informationWeek 10: Line Integrals
Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationMath 1431 Section M TH 4:00 PM 6:00 PM Susan Wheeler Office Hours: Wed 6:00 7:00 PM Online ***NOTE LABS ARE MON AND WED
Mth 43 Section 4839 M TH 4: PM 6: PM Susn Wheeler swheeler@mth.uh.edu Office Hours: Wed 6: 7: PM Online ***NOTE LABS ARE MON AND WED t :3 PM to 3: pm ONLINE Approimting the re under curve given the type
More informationMath 61CM - Solutions to homework 9
Mth 61CM - Solutions to homework 9 Cédric De Groote November 30 th, 2018 Problem 1: Recll tht the left limit of function f t point c is defined s follows: lim f(x) = l x c if for ny > 0 there exists δ
More informationf(x)dx . Show that there 1, 0 < x 1 does not exist a differentiable function g : [ 1, 1] R such that g (x) = f(x) for all
3 Definite Integrl 3.1 Introduction In school one comes cross the definition of the integrl of rel vlued function defined on closed nd bounded intervl [, b] between the limits nd b, i.e., f(x)dx s the
More informationConsequently, the temperature must be the same at each point in the cross section at x. Let:
HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the
More informationVariational Techniques for Sturm-Liouville Eigenvalue Problems
Vritionl Techniques for Sturm-Liouville Eigenvlue Problems Vlerie Cormni Deprtment of Mthemtics nd Sttistics University of Nebrsk, Lincoln Lincoln, NE 68588 Emil: vcormni@mth.unl.edu Rolf Ryhm Deprtment
More informationMath 360: A primitive integral and elementary functions
Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:
More informationThe Riemann Integral
Deprtment of Mthemtics King Sud University 2017-2018 Tble of contents 1 Anti-derivtive Function nd Indefinite Integrls 2 3 4 5 Indefinite Integrls & Anti-derivtive Function Definition Let f : I R be function
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More informationPart II. Analysis of PDE
Prt II. Anlysis of PDE Lecture notes for MA342H P. Krgeorgis pete@mths.tcd.ie 1/38 Second-order liner equtions Consider the liner opertor L(u) which is defined by L(u) = 1 u xx + 2 u yy +b 1 u x +b 2 u
More information1 The Riemann Integral
The Riemnn Integrl. An exmple leding to the notion of integrl (res) We know how to find (i.e. define) the re of rectngle (bse height), tringle ( (sum of res of tringles). But how do we find/define n re
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationDefinition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim
Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More informationNotes on length and conformal metrics
Notes on length nd conforml metrics We recll how to mesure the Eucliden distnce of n rc in the plne. Let α : [, b] R 2 be smooth (C ) rc. Tht is α(t) (x(t), y(t)) where x(t) nd y(t) re smooth rel vlued
More informationA. Limits - L Hopital s Rule ( ) How to find it: Try and find limits by traditional methods (plugging in). If you get 0 0 or!!, apply C.! 1 6 C.
A. Limits - L Hopitl s Rule Wht you re finding: L Hopitl s Rule is used to find limits of the form f ( x) lim where lim f x x! c g x ( ) = or lim f ( x) = limg( x) = ". ( ) x! c limg( x) = 0 x! c x! c
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationc n φ n (x), 0 < x < L, (1) n=1
SECTION : Fourier Series. MATH4. In section 4, we will study method clled Seprtion of Vribles for finding exct solutions to certin clss of prtil differentil equtions (PDEs. To do this, it will be necessry
More informationa < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1
Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More information1 The fundamental theorems of calculus.
The fundmentl theorems of clculus. The fundmentl theorems of clculus. Evluting definite integrls. The indefinite integrl- new nme for nti-derivtive. Differentiting integrls. Theorem Suppose f is continuous
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationGeneralizations of the Basic Functional
3 Generliztions of the Bsic Functionl 3 1 Chpter 3: GENERALIZATIONS OF THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 3.1 Functionls with Higher Order Derivtives.......... 3 3 3.2 Severl Dependent Vribles...............
More informationMath 31S. Rumbos Fall Solutions to Assignment #16
Mth 31S. Rumbos Fll 2016 1 Solutions to Assignment #16 1. Logistic Growth 1. Suppose tht the growth of certin niml popultion is governed by the differentil eqution 1000 dn N dt = 100 N, (1) where N(t)
More informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More information1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.
1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More information3 Conservation Laws, Constitutive Relations, and Some Classical PDEs
3 Conservtion Lws, Constitutive Reltions, nd Some Clssicl PDEs As topic between the introduction of PDEs nd strting to consider wys to solve them, this section introduces conservtion of mss nd its differentil
More informationW. We shall do so one by one, starting with I 1, and we shall do it greedily, trying
Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)
More informationg i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f
1. Appliction of functionl nlysis to PEs 1.1. Introduction. In this section we give little introduction to prtil differentil equtions. In prticulr we consider the problem u(x) = f(x) x, u(x) = x (1) where
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationLECTURE. INTEGRATION AND ANTIDERIVATIVE.
ANALYSIS FOR HIGH SCHOOL TEACHERS LECTURE. INTEGRATION AND ANTIDERIVATIVE. ROTHSCHILD CAESARIA COURSE, 2015/6 1. Integrtion Historiclly, it ws the problem of computing res nd volumes, tht triggered development
More informationx = b a N. (13-1) The set of points used to subdivide the range [a, b] (see Fig. 13.1) is
Jnury 28, 2002 13. The Integrl The concept of integrtion, nd the motivtion for developing this concept, were described in the previous chpter. Now we must define the integrl, crefully nd completely. According
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More informationSturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1
Ch.4. INTEGRAL EQUATIONS AND GREEN S FUNCTIONS Ronld B Guenther nd John W Lee, Prtil Differentil Equtions of Mthemticl Physics nd Integrl Equtions. Hildebrnd, Methods of Applied Mthemtics, second edition
More informationCalculus of Variations: The Direct Approach
Clculus of Vritions: The Direct Approch Lecture by Andrejs Treibergs, Notes by Bryn Wilson June 7, 2010 The originl lecture slides re vilble online t: http://www.mth.uth.edu/~treiberg/directmethodslides.pdf
More informationIntegrals along Curves.
Integrls long Curves. 1. Pth integrls. Let : [, b] R n be continuous function nd let be the imge ([, b]) of. We refer to both nd s curve. If we need to distinguish between the two we cll the function the
More informationSection 4.8. D v(t j 1 ) t. (4.8.1) j=1
Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationMAT612-REAL ANALYSIS RIEMANN STIELTJES INTEGRAL
MAT612-REAL ANALYSIS RIEMANN STIELTJES INTEGRAL DR. RITU AGARWAL MALVIYA NATIONAL INSTITUTE OF TECHNOLOGY, JAIPUR, INDIA-302017 Tble of Contents Contents Tble of Contents 1 1. Introduction 1 2. Prtition
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationdifferent methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).
Mth 1A with Professor Stnkov Worksheet, Discussion #41; Wednesdy, 12/6/217 GSI nme: Roy Zho Problems 1. Write the integrl 3 dx s limit of Riemnn sums. Write it using 2 intervls using the 1 x different
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationThe Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.
Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F
More informationChapters 4 & 5 Integrals & Applications
Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions
More informationAMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives
AMATH 731: Applied Functionl Anlysis Fll 214 Additionl notes on Fréchet derivtives (To ccompny Section 3.1 of the AMATH 731 Course Notes) Let X,Y be normed liner spces. The Fréchet derivtive of n opertor
More informationThe Riemann-Lebesgue Lemma
Physics 215 Winter 218 The Riemnn-Lebesgue Lemm The Riemnn Lebesgue Lemm is one of the most importnt results of Fourier nlysis nd symptotic nlysis. It hs mny physics pplictions, especilly in studies of
More information(4.1) D r v(t) ω(t, v(t))
1.4. Differentil inequlities. Let D r denote the right hnd derivtive of function. If ω(t, u) is sclr function of the sclrs t, u in some open connected set Ω, we sy tht function v(t), t < b, is solution
More information1. On some properties of definite integrals. We prove
This short collection of notes is intended to complement the textbook Anlisi Mtemtic 2 by Crl Mdern, published by Città Studi Editore, [M]. We refer to [M] for nottion nd the logicl stremline of the rguments.
More information4. Calculus of Variations
4. Clculus of Vritions Introduction - Typicl Problems The clculus of vritions generlises the theory of mxim nd minim. Exmple (): Shortest distnce between two points. On given surfce (e.g. plne), nd the
More informationMAC-solutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More informationZ b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...
Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.
More informationHeat flux and total heat
Het flux nd totl het John McCun Mrch 14, 2017 1 Introduction Yesterdy (if I remember correctly) Ms. Prsd sked me question bout the condition of insulted boundry for the 1D het eqution, nd (bsed on glnce
More informationFirst midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009
Mth 3B Review Steve 18 Mrch 2009 About the finl Fridy Mrch 20, 3pm-6pm, Lkretz 110 No notes, no book, no clcultor Ten questions Five review questions (Chpters 6,7,8) Five new questions (Chpters 9,10) No
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More information