t 2 B F x,t n dsdt t u x,t dxdt

Similar documents
2. Nonlinear Conservation Law Equations

An Introduction to Malliavin calculus and its applications

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Solutions to Assignment 1

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

System of Linear Differential Equations

Chapter 2. First Order Scalar Equations

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

10. State Space Methods

Chapter 6. Systems of First Order Linear Differential Equations

Math 334 Fall 2011 Homework 11 Solutions

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Class Meeting # 10: Introduction to the Wave Equation

Convergence of the Neumann series in higher norms

Hamilton Jacobi equations

EXERCISES FOR SECTION 1.5

Lecture 20: Riccati Equations and Least Squares Feedback Control

6.2 Transforms of Derivatives and Integrals.

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

1 Solutions to selected problems

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

4.6 One Dimensional Kinematics and Integration

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Undetermined coefficients for local fractional differential equations

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

Physics for Scientists & Engineers 2

Monotonic Solutions of a Class of Quadratic Singular Integral Equations of Volterra type

Two Coupled Oscillators / Normal Modes

Chapter Three Systems of Linear Differential Equations

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

Chapter 3 Boundary Value Problem

Announcements: Warm-up Exercise:

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity

CHAPTER 2 Signals And Spectra

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

15. Vector Valued Functions

DISCRETE GRONWALL LEMMA AND APPLICATIONS

The expectation value of the field operator.

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

A Sharp Existence and Uniqueness Theorem for Linear Fuchsian Partial Differential Equations

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

CONTRIBUTION TO IMPULSIVE EQUATIONS

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Math 333 Problem Set #2 Solution 14 February 2003

arxiv: v1 [math.ca] 15 Nov 2016

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

Example on p. 157

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

Differential Equations

Heat kernel and Harnack inequality on Riemannian manifolds

Positive continuous solution of a quadratic integral equation of fractional orders

KINEMATICS IN ONE DIMENSION

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

The Arcsine Distribution

ψ(t) = V x (0)V x (t)

Optimality Conditions for Unconstrained Problems

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

On Gronwall s Type Integral Inequalities with Singular Kernels

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Waves are naturally found in plasmas and have to be dealt with. This includes instabilities, fluctuations, waveinduced

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES

Unsteady Flow Problems

DIFFERENTIAL GEOMETRY HW 5

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Ordinary Differential Equations

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

6. Stochastic calculus with jump processes

CHARACTERIZATION OF REARRANGEMENT INVARIANT SPACES WITH FIXED POINTS FOR THE HARDY LITTLEWOOD MAXIMAL OPERATOR

Cash Flow Valuation Mode Lin Discrete Time

Admissible Solution for Hyperbolic Conservation Laws

Structural Dynamics and Earthquake Engineering

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0.

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation

Generalized Chebyshev polynomials

Endpoint Strichartz estimates

17. Weak and Strong Derivatives and Sobolev Spaces For this section, let Ω be an open subset of R d, p,q,r [1, ], L p (Ω) =

Math Final Exam Solutions

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

arxiv: v1 [math.pr] 19 Feb 2011

From Particles to Rigid Bodies

Empirical Process Theory

556: MATHEMATICAL STATISTICS I

1. Introduction The present paper is concerned with the scalar conservation law in one space dimension: u t + (A(u)) x = 0 (1.0.1)

Expert Advice for Amateurs

Some Basic Information about M-S-D Systems

MA 366 Review - Test # 1

Transcription:

Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion. Le he concenraion of conaminan a x, in U be denoed by ux, and le he vecor field F x, denoe he flux field for he conaminan. hen we can sae a conservaion principle for he conaminan as follows. For an arbirary ball B in U and an arbirary ime inerval 1, 0, we can asser ha ux, dx ux, 1 dx B F x, n dsd fx,dxd A B C D where he inegrals A,B represen he amoun of conaminan inside B a ime, 1, respecively, inegral C represens he amoun of conaminan flowing ou of B during he ime inerval 1,,(i.e., n here denoes he ouward uni normal o he boundary of B) and he inegral D represens he amoun of conaminan ha is creaed or desroyed inside B during he ime inerval 1, as a resul of a source densiy given by fx,. Using he divergence heorem o conver inegral C o an inegral over B, and using he he fundamenal heorem of calculus o wrie he difference A minus B as ux, dx ux, 1 dx ux,dxd we ge ux, divf f dxd 0 B Since his holds for all B in U and every ime inerval 1, 0,,we can conclude ha ux, divf x, fx, 0 x, U. Nex, we observe ha if he conaminan is being ranspored hrough he medium due o a combinaion of diffusion and convecion by a velociy field V x,, hen he flux field is relaed o he conaminan concenraion by, F x, ax, ux, ux,v x,, where ax, denoes a maerial dependen propery called "diffusiviy". he diffusiviy may be a ensor (n by n marix) in anisoropic maerials, a scalar in maerials which are isoropic bu no homogeneous, and for isoropic and homogeneous maerials, he diffusiviy is a consan. hen he conservaion saemen becomes, ux, div ax, u V udivv u fx, x, U. which means ha he concenraion saisfies an equaion of he form where ux, Lux, fx, in U Lux, divax,u V x, udivv x,u. he diffusiviy is only physically meaningful if i is nonnegaive, and we will assume furher ha i is sricly posiive a every poin in U.hen he operaor L is said o be ellipic and he operaor L is said o be parabolic. Evidenly, parabolic operaors arise in modelling he ime evoluion of physical processes ha behave like diffusion-convecion siuaion described here. 1

We are going o consider he problem ux, Lux, fx, ux, gx, ux,0 u 0 x in U on 0, on U 0. Such a problem is called an iniial value problem for an evoluion equaion, (IVP for shor). We shall show ha a modificaion of he approach we used for he ellipic BVP s can be employed o solve hese problems. Firs i will be necessary o describe he funcion spaces in which he variable is viewed differenly from he space variables, x. H Valued Funcions of In dealing wih evoluion equaions we will be faced wih funcions of x, where x is in R n and is in R 1 or in R 1. I will urn ou o be useful o view such funcions as funcions of, aking heir values in an appropriae space of funcions of he variable x; i.e., 0, U x ux, H u, funcion of wih values in H he funcion u(x,) will be wrien as u where his is inerpreed o mean ha he value a each is u,, a funcion which is an elemen of he funcion space H. For us, mos of he funcion spaces will be Hilber spaces. H Valued Disribuions: D [0,:H] For H a Hilber space, a linear mapping, J : C c 0, H defines a disribuion in wih values in H if and only if n 0 in C c 0, implies J n 0 in H We will use he noaion D0, o denoe C c 0,. For f : 0, H, locally inegrable, define J f 0 fd D0,. hen J f is easily seen o be a disribuion in wih values in H. he noaion for his will be J f D 0, : H, and we shall use he noaion L 1 loc 0, : H o denoe locally inegrable funcions of, wih values in H. hese are he regular H-valued disribuions. For J D 0, : H, he derivaive, J D 0, : H, is defined by J for all D0,. More generally, n 1 n n for all D0,. Now define L 0, : H o be he se of all measurable funcions of, 0, wih values in H such ha

u L 0,:H 0 uh d. heorem 1- he naural inclusion, L 0, : H D 0, : H, is a coninuous injecion; i.e., f n f in L 0, : H implies J fn f in H D0, and f 0 in H D0, implies f 0 in L 0, : H If we define linear spaces wih norms, C0, : H coninuous funcions of, 0, wih values in H u C0,:H max 0 u H H 1 0, : H u L 0, : H : u L 0, : H u H 1 0,:H 0 uh u H d hen i is shown on p 86 in he ex, ha for any u in H 1 0, : H, 1. here exiss û in C0, : H such ha u û a.e., in (0,). ii) for all s, 0,, u us 0 u zdz 3. iii) u C0,:H Cu H 1 0,:H he argumens are very close o hose ha show ha funcions in H 1 a,b are coninuous. For Hilber spaces, H 1 and H wih H 1 coninuously included in H i is easy o show ha In paricular, for he siuaion D 0, : H 1 D 0, : H V H V f L 0, : V D 0, : V implies f D 0, : V D 0, : V In he evoluion equaions we are going o sudy, he ellipic operaor L is a bounded linear map from V ono is dual, V. hen a naural choice of soluion space for he evoluion equaions is W0, u L 0, : V : u L 0, : V i.e., if u is in W0, hen u and Lu boh belong o L 0, : V. heorem - W0, is a Hilber space for he inner produc f,g W0. 0 f,gv f,g V proof- Suppose f n is a Cauchy sequence in W0,. hen he definiion of he inner produc implies 3

f n F 0 in L 0, : V : and f n G 0 in L 0, : V Now L 0, : V D 0, : V implies Bu hence However, f n F 0 in D 0, : V : and f n F 0 in D 0, : V D 0, : V D 0, : V, f n F 0 in D 0, : V. f n G 0 in L 0, : V D 0, : V, and i follows ha F 0 G 0 and F 0 W0,. Since we are sudying equaions of he form ulu for L an ellipic operaor of order, we will be paricularly ineresed in he case of W0, when V H 0 1 U (and H H 0 U, V H 1 U. I is no hard o show, using mollifiers, ha in his case C 0, : H 0 1 U is dense in W0,. hen we have he following resul, heorem 3- he naural injecion C 0, : H 0 1 U C0, : H 0 U can be exended as a coninuous injecion of Moreover, if u W0,,hen W0, C0, : H 0 U i) u, H is absoluely coninuous on 0, ii) d d u, H u,u V V and d d u,v H u,v V V v,u V V iii) here is a consan C 0 such ha u C0,:H Cu W0, proof-heurisically, his resul is nohing more han he observaion ha for u W0,, u, H u,,u, H hence d d u, H u,,u, H u,u V V where we have noed ha u L 0, : V and ha he dualiy pairing for V V is jus he Hinner produc. his can all be made rigorous as follows. Le u W0, u L 0, : H 0 1 U : u L 0, : H 1 U Exend u o, as ũ u 0 u 0 Clearly his defines a coninuous injecion of W0, ino W,. Now le a denoe a 4

smooh cu-off funcion; i.e., a C R a 1 for 0, a 0 for, and 0 a 1 for 0. Now if u C 0, : H 0 1 U, hen d/daũ, H d/daũ,,aũ, H Bu aũ,aũ V V d/daũ,s H ds aũ, H 0 from which i follows ha aũ, H aũ,s,aũ,s V V ds aũ,s aũ,s H 1 U H 1 0 U ds aũ,s H 1 U C a u,s 0 H 1 U C a u W0, aũ,s H0 1 Uds u,s H0 1 Uds where C a denoes a consan depending on he cu-off funcion, a. Now i follows ha max 0 u, H C a u W0, and, since C 0, : H 0 1 U is dense in W0,, his holds no jus for u in C 0, : H 0 1 U, bu for all u in W0,. he meaning of his resul is ha for every u in W0, here exiss a û C0, : H 0 U such ha u û a.e.,0. In paricular, for u in W0,, a soluion of an IVP, i means ha as a coninuous funcion of wih values in H 0 U, as ends o zero, u ends o u0 in H 0 U. his gives meaning o he iniial condiion, u0 u 0 H 0 U. Wih regard o boundary condiions associaed wih an IVP, any u such ha, u W0, u L 0, : H 0 1 U : u L 0, : H 1 U auomaically saisfies he homogeneous boundary condiion, 0 u u 0. An inhomogeneous boundary condiion of he form, 0 u u g, for a funcion, g, L 0, : H 1/ can be formulaed as follows. Le Gx, be such ha G L 0, : H 1 U, G L 0, : H 1 U G, G,0 in H 0 U as 0, 0 G g, a.e., 0. hen u L 0, : H 1 U saisfies he boundary condiion 0 u g, if 5

u, G, H 0 1 U,a.e., 0 ; i.e., ug u L 0, : H 0 1 U : u L 0, : H 1 U. I is worh noing ha here is no condiion of compaibiliy required beween u 0 and G,0. While he condiions u, G, H 0 1 U,a.e., 0 ; u, u 0, in H 0 U as 0 seem o require ha u 0 G, H 0 1 U, i is clear from he fac ha u 0 H 0 U ha his condiion has no meaning. In fac, if h is any funcion such ha and h u L 0, : H 0 1 U : u L 0, : H 1 U h,0 G,0, hen replacing G in he IVP wih Gh would have no effec on he boundary condiion; i.e., he value of G a 0 is no relevan o he formulaion of he IVP. For example, consider he IVP in 1-d, ux, xx ux, 0, 0 x 1, 0, ux,0 1 0 x 1, u0, u1, 0, 0. Alhough he incompaibiliy beween he iniial and boundary condiions prevens he exisence of a soluion which is coninuous on he closure of U, i does no preclude he exisence of a soluion in W0,. his soluion is hen coninuous as a funcion of wih values in H 0 U, which means ha as ends o zero, ux, ends o u 0 x 1, in he norm of H 0 U. A he same ime, for each posiive, ux, is valued in L 0, : H 0 1 U, which means ha 0 u, 0; i.e., he boundary condiion is saisfied. ha his implies no conradicion a 0,0 and 1,0 is a resul of he fac ha he funcions involved are no poinwise defined bu are defined only up o ses of measure zero. 6