Scalar Conservation Laws and First Order Equations Introduction. Consider equations of the form. (1) u t + q(u) x =0, x R, t > 0.
|
|
- Percival Parrish
- 5 years ago
- Views:
Transcription
1 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 quantity Q and q(u) its flux function. Equation (1) constitutes a link between density and flux and expresses a (scalar) conservation law for the following reason. If we consider a control interval [x 1,x 2 ], the integral x2 x 1 u(x, t)dt gives the amount of Q between x 1 and x 2 at time t. A conservation law states that, without sources or sinks, the rate of changes of Q in the interior of [x 1,x 2 ] is determined by the net flux through the end points of the interval. If the flux is modelled by a function q = q(u), the law translates into the equation (2) d x2 u(x, t)dx = q(u(x 2,t)) + q(u(x 1,t)), dt x 1 where we assume that q>(q<) for a flux along the positive (negative) direction of the x axes. Ifu are q are smooth functions, equation (2) can be written in the form x2 x 1 [u t (x, t)+q(u(x, t)) x ]dx =, which implies (1), due to the arbitrariness of the interval [x 1,x 2 ]. In the next section we go back to the model of pollution in a channel, neglecting the diffusion and choosing for q a linear function of u, namely, q(u) = vu, where v is constant. The result is a pure transport model, in which the vector vi is the advection speed. In the sequel, to introduce and motivate some important concepts and results, we shall use a nonlinear model from traffic dynamics, with speed v depending on u. Typeset by AMS-TEX 1
2 2 Linear Transport Equation 1. Pollution In a Channel Go back to the simple model for the evolution of a pollutant in a narrow channel. When diffusion and transport are both relevant we have the equation c t = Dc xx vc x where c is the concentration of the pollutant and vi is the stream speed (c >, constant). We discuss here the case of the pure transport equation (4) c t + vc x =, i.e. D =. Introducing the vector v = vi + j equation (4) can be written in the form vc x + c t = c v =, pointing out the orthogonality of c and v. But c is orthogonal to the level lines of c, along which c is constant. Therefore the level lines of c are the straight lines parallel to v, of equation x = v t + x. These straight lines are called characteristics. Let us compute c along the characteristic x = vt + x, letting w(t) =c(x + vt, t). Since ẇ(t) =vc x (x + vt, t)+c t (x + vt, t), equation (4) becomes an o.d.e. ẇ(t) =, which implies that c is constant along the characteristic. We want to determine the evolution of the concentration c, by knowing its initial profile c(x, ) = g(x). The method to compute the solution at a point (x, t), t>, is very simple. Let x = vt + x be the equation of the characteristic passing through (x, t). Go back in time along this characteristic from (x, t) until the point (x, ), of intersection with the x-axis. Since c is constant along the characteristic and c(x, ) = g(x ), we have c(x, t) =g(x )=g(x vt). Thus, if g C 1 (R), the solution of the initial value problem (4), (5) is given by (6) c(x, t) =g(x vt). The solution (6) represents a travelling wave, moving with speed v in the positive x-direction.
3 3 2. Distributed Source The equation (7) c t + vc x = f(x, t), (8) c(x, ) = g(x), describes the effect of an external distributed source along the channel. The function f represents the intensity of the source, measured in concentration per unit time. Again, to compute the value of the solution u at a point (x, t) is not difficult. Let x = x + vt be the characteristic passing through (x, t) and compute u along this characteristic, setting From (7), w satisfies the o.d.e. Thus w(t) =c(x + vt, t). ẇ(t) =vc x (x + vt, t)+c t (x + vt, t) =f(x + vt, t) w(t) =g(x )+ w() = g(x ). f(x + vs, s)ds. Letting t = t and recalling that x = x vt, we get c(x, t) =w(t) =g(x vt)+ f(x v(t s),s)ds. Proposition 1. Let g C 1 (R) and f, f x C(R R + ). The solution of the initial value problem is given by the formula c t + vc x = f(x, t), x R, t >, c(x, ) = g(x) x R. c(x, t) =g(x vt)+ f(x v(t s),s)ds.
4 4 3. Decay And Localized Source Suppose that, due to biological decomposition, the pollutant decays at the rate r(x, t) = γc(x, t), γ >. Without external sources and diffusion, the mathematical model is c t + vc x = γc, Setting c(x, ) = g(x). (9) u(x, t) =c(x, t)e γ v x we have u x = (c x + γv c ) e γ v x and u t = c t e γ v x and therefore the equation for u is From Proposition 1, we obtain and from (9) u t + vu x = u(x, ) = g(x)e γ v x. u(x, t) =g(x vt)e γ v (x vt) c(x, t) =g(x vt)e γt, which is a damped travelling wave. We now examine the effect of a source of pollutant placed at a certain point of the channel, e.g. at x =. Typically, one can think of waste material from industrial machineries. Before the machine start working, for instance before time t =, we assume that the channel is clean. We want to determine the pollutant concentration, supposing that at x = it is kept at a constant level β>, for t>. To model this source we introduce the Heaviside function { 1, t H(t) =, t <, with the boundary condition (1) c(,t)=βh(t)
5 5 and the initial condition (11) c(x, )=, for x>. As before, let u(x, t) =c(x, t)e γ v x, which is a solution of Then u t + vu x =. u(x, ) =c(x, )e γ v x =, x > u(,t)=c(,t)=βh(t). Since u is constant along the characteristics, it must be of the form (12) u(x, t) =u (x vt) where u is to be determined from the boundary condition (1) and the initial condition (11). To compute u for x<vt, observe that a characteristic leaving the t-axis from a point (x, t) carries the data βh(t). Therefore, we must have Letting s = vt, we obtain u ( vt) =βh(t). u (s) =βh( s v ) and from (12) we obtain u(x, t) =βh(t x v ). This formula gives the solution also in the sector x>vt, t> since the characteristics leaving the x-axis carry zero data and hence we deduce u = c =. This means that the pollutant has not yet reached the point x at time t, ifx>vt. Finally, recalling (9), we find c(x, t) =βh(t x v )e γ v x. Observe that in (, ) there is a jump discontinuity which is transported along the characteristic x = vt.
6 6 4. Inflow And Outflow Characteristics. A Stability Estimate. The domain in the localized source problem is the quadrant x>, t>. To uniquely determine the solution we have used the initial data on the x-axis, x> and the boundary data on the t-axis, t>. The solution to this problem is uniquely determined. This is due to the fact that, since v>, when time increases, all the characteristics carry the information (the data) towards the interior of the quadrant x>, t>. In other words, the characteristics are inflow characteristics. More generally, consider the equation u t + au x = f(x, t) in the domain x>, t>, where a is constant (a ). The characteristics are the lines x at = constant. (i) If a>, we are in the case of the pollutant model: all the characteristics are inflow and the data must be assigned on both semi-axes. (ii) If a<, the characteristics leaving the x-axis are inflow, while those leaving the t-axis are outflow. In these case the initial data alone are sufficient to uniquely determine the solution, while no data has to be assigned on the semi-axis x =, t>. Coherently, a problem in the half-strip <x<r, t>, besides the initial data, requires a data assigning on the inflow boundary, namely, { u(,t)=h (t) if a>, u(r, t) =h R (t) if a<. The solution to the resulting initial-boundary-value problem is uniquely determined at every point in the strip by its values along the characteristics. Moreover a stability estimate can be proved as follows. Consider, for instance, the case a> and the problem u t + au x = <x<r, t>, (13) u(,t)=h(t) t>, u(x, ) = g(x) <x<r. Multiply the differential equation by u and write =uu t + auu x = 1 d 2 dt u2 + a d 2 dx u2 =. Integrating in x over (,R) we obtain d dt u 2 (x, t)dx + a[u 2 (R, t) u 2 (,t)] =.
7 7 Now use the data u(,t)=h(t) and the positivity of a to obtain d dt u 2 (x, t)dx ah 2 (t). Integrating in t we have, using the initial condition u(x, ) = g(x), (14) u 2 (x, t)dx g 2 (x, t)dx + a h 2 (s)ds. Now, let u 1 and u 2 be solutions of problem (13) with initial data g 1, g 2 and boundary data h 1, h 2 on x =. Then, by linearity, w = u 1 u 2 is a solution of the problem (13), with initial data h 1 h 2 on x =. Applying the inequality (14) to w, we have [u 1 (x, t) u 2 (x, t)] 2 dx [g 1 (x) g 2 (x)] 2 dx + a [h 1 (s) h 2 (s)] 2 ds. Thus, a least-square approximation of the data controls a least-square approximation of the corresponding solutions. In this sense, the solution of the problem depends continuously on the initial data and on the boundary data on x =. Remark. The values of u on x = R do not appear on (14).
Numerical 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 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 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 informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
More informationProblem Set 1. This week. Please read all of Chapter 1 in the Strauss text.
Math 425, Spring 2015 Jerry L. Kazdan Problem Set 1 Due: Thurs. Jan. 22 in class. [Late papers will be accepted until 1:00 PM Friday.] This is rust remover. It is essentially Homework Set 0 with a few
More information1.5 First Order PDEs and Method of Characteristics
1.5. FIRST ORDER PDES AND METHOD OF CHARACTERISTICS 35 1.5 First Order PDEs and Method of Characteristics We finish this introductory chapter by discussing the solutions of some first order PDEs, more
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then 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 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 informationModeling 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 informationDiffusion equation in one spatial variable Cauchy problem. u(x, 0) = φ(x)
Diffusion equation in one spatial variable Cauchy problem. u t (x, t) k u xx (x, t) = f(x, t), x R, t > u(x, ) = φ(x) 1 Some more mathematics { if x < Θ(x) = 1 if x > is the Heaviside step function. It
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 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 informationA symmetry-based method for constructing nonlocally related partial differential equation systems
A symmetry-based method for constructing nonlocally related partial differential equation systems George W. Bluman and Zhengzheng Yang Citation: Journal of Mathematical Physics 54, 093504 (2013); doi:
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 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 informationUNIT 3 INTEGRATION 3.0 INTRODUCTION 3.1 OBJECTIVES. Structure
Calculus UNIT 3 INTEGRATION Structure 3.0 Introduction 3.1 Objectives 3.2 Basic Integration Rules 3.3 Integration by Substitution 3.4 Integration of Rational Functions 3.5 Integration by Parts 3.6 Answers
More informationAP Calculus AB. Review for Test: Applications of Integration
Name Review for Test: Applications of Integration AP Calculus AB Test Topics: Mean Value Theorem for Integrals (section 4.4) Average Value of a Function (manipulation of MVT for Integrals) (section 4.4)
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationLecture Notes on Hyperbolic Conservation Laws
Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More informationExponential Energy Decay for the Kadomtsev-Petviashvili (KP-II) equation
São Paulo Journal of Mathematical Sciences 5, (11), 135 148 Exponential Energy Decay for the Kadomtsev-Petviashvili (KP-II) equation Diogo A. Gomes Department of Mathematics, CAMGSD, IST 149 1 Av. Rovisco
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 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 informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationHyperbolic PDEs. Chapter 6
Chapter 6 Hyperbolic PDEs In this chapter we will prove existence, uniqueness, and continuous dependence of solutions to hyperbolic PDEs in a variety of domains. To get a feel for what we might expect,
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS Basic concepts: Find y(x) where x is the independent and y the dependent varible, based on an equation involving x, y(x), y 0 (x),...e.g.: y 00 (x) = 1+y(x) y0 (x) 1+x or,
More informationLecture 11: Non-linear Diffusion
Lecture 11: Non-linear Diffusion Scribe: Lou Odette - American International Group (AIG) October 17, 006 1 Non-linear Drift In the continuum limit the PDF ρ(x, t) for the x at time t of a single random
More informationMath 220A - Fall 2002 Homework 5 Solutions
Math 0A - Fall 00 Homework 5 Solutions. Consider the initial-value problem for the hyperbolic equation u tt + u xt 0u xx 0 < x 0 u t (x, 0) ψ(x). Use energy methods to show that the domain of dependence
More informationMATH 220: PROBLEM SET 1, SOLUTIONS DUE FRIDAY, OCTOBER 2, 2015
MATH 220: PROBLEM SET 1, SOLUTIONS DUE FRIDAY, OCTOBER 2, 2015 Problem 1 Classify the following PDEs by degree of non-linearity (linear, semilinear, quasilinear, fully nonlinear: (1 (cos x u x + u y =
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
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 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 information1. The accumulated net change function or area-so-far function
Name: Section: Names of collaborators: Main Points: 1. The accumulated net change function ( area-so-far function) 2. Connection to antiderivative functions: the Fundamental Theorem of Calculus 3. Evaluating
More informationFirst order Partial Differential equations
First order Partial Differential equations 0.1 Introduction Definition 0.1.1 A Partial Deferential equation is called linear if the dependent variable and all its derivatives have degree one and not multiple
More informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationStereographic projection and inverse geometry
Stereographic projection and inverse geometry The conformal property of stereographic projections can be established fairly efficiently using the concepts and methods of inverse geometry. This topic is
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 informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
More informationSeparation of variables
Separation of variables Idea: Transform a PDE of 2 variables into a pair of ODEs Example : Find the general solution of u x u y = 0 Step. Assume that u(x,y) = G(x)H(y), i.e., u can be written as the product
More informationVortex stretching in incompressible and compressible fluids
Vortex stretching in incompressible and compressible fluids Esteban G. Tabak, Fluid Dynamics II, Spring 00 1 Introduction The primitive form of the incompressible Euler equations is given by ( ) du P =
More informationS. Ghorai 1. Lecture IV Linear equations, Bernoulli equations, Orthogonal trajectories, Oblique trajectories. (p(x)y r(x))dx + dy = 0.
S. Ghorai 1 Lecture IV Linear equations, Bernoulli equations, Orthogonal trajectories, Oblique trajectories 1 Linear equations A first order linear equations is of the form This can be written as Here
More information14 Divergence, flux, Laplacian
Tel Aviv University, 204/5 Analysis-III,IV 240 4 Divergence, flux, Laplacian 4a What is the problem................ 240 4b Integral of derivative (again)........... 242 4c Divergence and flux.................
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationMath 5520 Homework 2 Solutions
Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas
More informationSolving First Order PDEs
Solving Ryan C. Trinity University Partial Differential Equations Lecture 2 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.
More informationsurface area per unit time is w(x, t). Derive the partial differential
1.2 Conduction of Heat in One-Dimension 11 1.2.4. Derive the diffusion equation for a chemical pollutant. (a) Consider the total amount of the chemical in a thin region between x and x + Ax. (b) Consider
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationQuasi-linear first order equations. Consider the nonlinear transport equation
Quasi-linear first order equations Consider the nonlinear transport equation u t + c(u)u x = 0, u(x, 0) = f (x) < x < Quasi-linear first order equations Consider the nonlinear transport equation u t +
More informationModeling Rarefaction and Shock waves
April 30, 2013 Inroduction Inroduction Rarefaction and shock waves are combinations of two wave fronts created from the initial disturbance of the medium. Inroduction Rarefaction and shock waves are combinations
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationMath 207 Honors Calculus III Final Exam Solutions
Math 207 Honors Calculus III Final Exam Solutions PART I. Problem 1. A particle moves in the 3-dimensional space so that its velocity v(t) and acceleration a(t) satisfy v(0) = 3j and v(t) a(t) = t 3 for
More informationMATH 6B Spring 2017 Vector Calculus II Study Guide Final Exam Chapters 8, 9, and Sections 11.1, 11.2, 11.7, 12.2, 12.3.
MATH 6B pring 2017 Vector Calculus II tudy Guide Final Exam Chapters 8, 9, and ections 11.1, 11.2, 11.7, 12.2, 12.3. Before starting with the summary of the main concepts covered in the quarter, here there
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
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 informationSolving First Order PDEs
Solving Ryan C. Trinity University Partial Differential Equations January 21, 2014 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.
More 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 information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
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 informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationMATH 425, HOMEWORK 5, SOLUTIONS
MATH 425, HOMEWORK 5, SOLUTIONS Exercise (Uniqueness for the heat equation on R) Suppose that the functions u, u 2 : R x R t R solve: t u k 2 xu = 0, x R, t > 0 u (x, 0) = φ(x), x R and t u 2 k 2 xu 2
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationAPPLICATIONS OF DIFFERENTIABILITY IN R n.
APPLICATIONS OF DIFFERENTIABILITY IN R n. MATANIA BEN-ARTZI April 2015 Functions here are defined on a subset T R n and take values in R m, where m can be smaller, equal or greater than n. The (open) ball
More informationMATH 320, WEEK 4: Exact Differential Equations, Applications
MATH 320, WEEK 4: Exact Differential Equations, Applications 1 Exact Differential Equations We saw that the trick for first-order differential equations was to recognize the general property that the product
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationCross-diffusion models in Ecology
Cross-diffusion models in Ecology Gonzalo Galiano Dpt. of Mathematics -University of Oviedo (University of Oviedo) Review on cross-diffusion 1 / 42 Outline 1 Introduction: The SKT and BT models The BT
More informationTHE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK. Summer Examination 2009.
OLLSCOIL NA héireann, CORCAIGH THE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK Summer Examination 2009 First Engineering MA008 Calculus and Linear Algebra
More informationSolution of Math132 Exam1
Solution of Math3 Exam. ( %) (a) A bacteria has population at time t = an its rate of growth is (t 3 + t + ) bacteria per hour after t hours. What is the population after 4 hours? Suppose p(t) is the bacteria
More informationMATH1013 Calculus I. Introduction to Functions 1
MATH1013 Calculus I Introduction to Functions 1 Edmund Y. M. Chiang Department of Mathematics Hong Kong University of Science & Technology May 9, 2013 Integration I (Chapter 4) 2013 1 Based on Briggs,
More informationAlgebra I EOC Review (Part 2)
1. Let x = total miles the car can travel Answer: x 22 = 18 or x 18 = 22 2. A = 1 2 ah 1 2 bh A = 1 h(a b) 2 2A = h(a b) 2A = h a b Note that when solving for a variable that appears more than once, consider
More informationPropagation of discontinuities in solutions of First Order Partial Differential Equations
Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in
More informationSalmon: Lectures on partial differential equations
4 Burger s equation In Lecture 2 we remarked that if the coefficients in u x, y,! "! "x + v x,y,! "! "y = 0 depend not only on x,y but also on!, then the characteristics may cross and the solutions become
More informationChapter 3. Second Order Linear PDEs
Chapter 3. Second Order Linear PDEs 3.1 Introduction The general class of second order linear PDEs are of the form: ax, y)u xx + bx, y)u xy + cx, y)u yy + dx, y)u x + ex, y)u y + f x, y)u = gx, y). 3.1)
More informatione (x y)2 /4kt φ(y) dy, for t > 0. (4)
Math 24A October 26, 2 Viktor Grigoryan Heat equation: interpretation of the solution Last time we considered the IVP for the heat equation on the whole line { ut ku xx = ( < x
More informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
More information3.1 Day 1: The Derivative of a Function
A P Calculus 3.1 Day 1: The Derivative of a Function I CAN DEFINE A DERIVATIVE AND UNDERSTAND ITS NOTATION. Last chapter we learned to find the slope of a tangent line to a point on a graph by using a
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationAPPM GRADUATE PRELIMINARY EXAMINATION PARTIAL DIFFERENTIAL EQUATIONS SOLUTIONS
Thursday August 24, 217, 1AM 1PM There are five problems. Solve any four of the five problems. Each problem is worth 25 points. On the front of your bluebook please write: (1) your name and (2) a grading
More information2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.
Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following
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 informationMAS3706 Topology. Revision Lectures, May I do not answer enquiries as to what material will be in the exam.
MAS3706 Topology Revision Lectures, May 208 Z.A.Lykova It is essential that you read and try to understand the lecture notes from the beginning to the end. Many questions from the exam paper will be similar
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More information10. Smooth Varieties. 82 Andreas Gathmann
82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationDEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF MASSACHUSETTS. MATH 233 SOME SOLUTIONS TO EXAM 2 Fall 2018
DEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF MASSACHUSETTS MATH 233 SOME SOLUTIONS TO EXAM 2 Fall 208 Version A refers to the regular exam and Version B to the make-up. Version A. A particle
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationAPPLICATIONS OF DIFFERENTIATION
4 APPLICATIONS OF DIFFERENTIATION APPLICATIONS OF DIFFERENTIATION 4.9 Antiderivatives In this section, we will learn about: Antiderivatives and how they are useful in solving certain scientific problems.
More informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
More informationScience One Math. January
Science One Math January 10 2018 (last time) The Fundamental Theorem of Calculus (FTC) Let f be continuous on an interval I containing a. 1. Define F(x) = f t dt with F (x) = f(x). on I. Then F is differentiable
More information