Dispersive Systems. 1) Schrödinger equation 2) Cubic Schrödinger 3) KdV 4) Discreterised hyperbolic equation 5) Discrete systems.

Similar documents
Scalar Conservation Laws

Well-posedness of the generalized Proudman-Johnson equation without viscosity. Hisashi Okamoto RIMS, Kyoto University

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations-I!

, u denotes uxt (,) and u. mean first partial derivatives of u with respect to x and t, respectively. Equation (1.1) can be simply written as

Chapter 2. First Order Scalar Equations

4.2 Continuous-Time Systems and Processes Problem Definition Let the state variable representation of a linear system be

PH2130 Mathematical Methods Lab 3. z x

1 First Order Partial Differential Equations

Huazhong Tang 1 and Gerald Warnecke Introduction ANOTEON(2K + 1)-POINT CONSERVATIVE MONOTONE SCHEMES

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

4.6 One Dimensional Kinematics and Integration

EXERCISES FOR SECTION 1.5

DESIGN OF TENSION MEMBERS

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

2. Nonlinear Conservation Law Equations

first-order circuit Complete response can be regarded as the superposition of zero-input response and zero-state response.

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

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

Class Meeting # 10: Introduction to the Wave Equation

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Modelling Traffic Flow with Constant Speed using the Galerkin Finite Element Method

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

Riemann Solvers and Numerical Methods for Fluid Dynamics

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

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

Chapter 3 Boundary Value Problem

System of Linear Differential Equations

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

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

Convergence of the Neumann series in higher norms

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1

Existence Theory of Second Order Random Differential Equations

ES.1803 Topic 22 Notes Jeremy Orloff

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

Predator - Prey Model Trajectories and the nonlinear conservation law

Two Coupled Oscillators / Normal Modes

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

Linearly implicit structure-preserving schemes for Hamiltonian systems

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Chapter 1 Fundamental Concepts

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Localization and Map Making

Compatible Versus Regular Well-Posed Linear Systems

On Gronwall s Type Integral Inequalities with Singular Kernels

Example on p. 157

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

Traveling Waves. Chapter Introduction

DISCRETE GRONWALL LEMMA AND APPLICATIONS

Applied Mathematics Letters. Oscillation results for fourth-order nonlinear dynamic equations

Hamilton Jacobi equations

Homotopy Perturbation Method for Solving Some Initial Boundary Value Problems with Non Local Conditions

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).

Lecture 20: Riccati Equations and Least Squares Feedback Control

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

Asymptotic Solution of the Anti-Plane Problem for a Two-Dimensional Lattice

Math Final Exam Solutions

û s L u t 0 s a ; i.e., û s 0

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

Some Basic Information about M-S-D Systems

Conservation Laws and Hamiltonian Symmetries of Whitham-Broer-Kaup Equations

Riemann Function and Methods of Group Analysis

Chapter 8 The Complete Response of RL and RC Circuits

ME 391 Mechanical Engineering Analysis

Surfaces in the space E

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

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

Journal of Mathematical Analysis and Applications

ON THE OSCILLATION OF THIRD ORDER FUNCTIONAL DIFFERENTIAL EQUATIONS. Cairo University, Orman, Giza 12221, Egypt

Physics 1402: Lecture 22 Today s Agenda

PROJECTS WITH APPLICATIONS OF DIFFERENTIAL EQUATIONS AND MATLAB

-e x ( 0!x+1! ) -e x 0!x 2 +1!x+2! e t dt, the following expressions hold. t

Fractional Method of Characteristics for Fractional Partial Differential Equations

6. Stochastic calculus with jump processes

Chapter 5: Control Volume Approach and Continuity Principle Dr Ali Jawarneh

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

Cash Flow Valuation Mode Lin Discrete Time

On Oscillation of a Generalized Logistic Equation with Several Delays

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

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

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

The Application of Optimal Homotopy Asymptotic Method for One-Dimensional Heat and Advection- Diffusion Equations

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

SELBERG S CENTRAL LIMIT THEOREM ON THE CRITICAL LINE AND THE LERCH ZETA-FUNCTION. II

Solutions to Assignment 1

Math 10B: Mock Mid II. April 13, 2016

Predator-Prey Models with Diffusion and Time Delay

The Euler-Lagrange Approach for Steady and Unsteady Flows. M. Sommerfeld. www-mvt.iw.uni-halle.de. Title. Zentrum für Ingenieurwissenschaften

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

Some Regularity Properties of Three Dimensional Incompressible Magnetohydrodynamic Flows

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

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

A Direct Method for Solving Nonlinear PDEs and. New Exact Solutions for Some Examples

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

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

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

Transcription:

Dispersive Sysems 1) Schrödinger eqaion ) Cbic Schrödinger 3) KdV 4) Discreerised hyperbolic eqaion 5) Discree sysems KdV + + ε =, = ( ) ( )

d d + = d d =, =. ( ) = ( ) DISCONTINUITY, prescribed cri Collision of characerisic lines ( ) > ( ) 1 No coninos solion for > cri

Solions of eq.-s of compressible flow develop disconiniies (shocks). Sokes, Airy, Riemann: Conservaion laws + f + g + h = y z ( ) ( ) f = f, g = g, ec.

Inegrae over G: + = d dv F nds d G G F = f, g, h fl ( ), Weak solions, disribions + f = d d b a d + f f = a b a. y (), dy b d = s

y b y b d d + d = [ ] d + d + s l r = a y a y d y a b [ ] fd fd+ s = l r y f + f f + f + s[ ], l a b r l r So by he conservaion law, s = [ f ] [ ]

Rankine-Hgonio condiion Eample: + = f ( ) ( 1 ) 1 + = = [ ] [ ] e r e+ r 1 s = f = = e r Enropy condiion: e > r

1 + = ( ) + = log a + = RH condiion s = log e + r e e r log r

Mliply eqaion + = ( ) 3 3 R-H condiion s + = + + = = 3 3 3 3 3 e r e e r r e r e + r Which conservaion law??

Goodman-La; in + = discreise space:, k kδ ( ) ( ) Replace by qoien: symmeric difference k + k = Δ d =. d k+ 1 k 1 Dispersive Conservaion form

f k f f k+ 1 1 k+ + = Δ = k+ k k+ 1 1 3 (,) = + cos Nmerical eperimens reveal oscillaions and weak convergence for >, cri

1. Eplain oscillaory behavior. Prove weak convergence 3. Deermine weak limi 4. Wha eqaion does he weak limi saisfy?

f k k k + + k+ 1 k = k Δ f f k+ 1 1 k + = Δ 1 = 1 k+ 1 Divide by : k log + k+ 1 k 1 = k Δ k

g g k 1 1 log + k + = k Δ g k + + = k k 1 1 +

( Δ) If wold end srongly o, boh and ( ) + f = ( ) g( ) log + = wold be saisfied. B hey are incompaible for disconinos solions.

The solion of he difference eqaion desperaely ries o saisfy wo incompaible conservaion laws. Hence he oscillaions.

Limis of Semigrops. Denoe by U(, ε ) he solion operaor for he KdV eqaion + + (1) ε =, ha is, he operaor U, :,;, ; ( ε ) ( ε) ( ε) ha relaes iniial vales of solions of (1) o heir vales a. Clearly, he operaors

U(, ε ) form a semigrop (even a grop) in he parameer : (, ) (, ) = ( +, ) U s U U s ε ε ε () La & Levermore have shown ha he limi ε of U(, ε ) eis in he weak b for > cri no in he srong opology: U (, ) ( ) ε U (3). Noe ha since eqaion (1), wrien in conservaion form

1 ( ) + + ε =, is nonlinear, a weak limi of is solions ha is no a srong limi is no a solion of he limiing eqaion 1 ( ) + =. (4) So U( ) in () is no he solion of (4) wih iniial vale. Qesion. Do he limi operaors U( ) defined in (3) form a semigrop:

( ) ( ) ( ). U U s U + s (5) L & L have derived an eplici formla for U ( ), which shows ha U ( ) is no coninos in he weak opology. Tha makes i very dbios ha he limi as ε ends zero of () is (5). The qesion can be easily decided by a single calclaion.

There is an analoge of his qesion in he heory of rblence. Here operaors o be considered are he solions operaors U(, R ) for he Navier-Sokes eqaion. I srmise ha as he Reynolds nmber R, he operaors U(, R ) end weakly b no srongly o a limi U( ).

Since he limiing Eler eqaions are nonlinear conservaion laws, U ( ) is no he solion operaor for he incompressible Eler eqaion b a descripion of rblen flow. In analogy wih he zero dispersion limi of he KdV eqaion I wold srmise ha he operaors U ( ), he limis of U(, R ), do no form a semigrop, as conjecred by Heinz Kreiss.

Compacness,, + = can be wrien as conservaion law ( 1 ) + =. Iniial vale (,) ( ) zero oside a finie inerval. =, bonded and Viscosiy mehod + = ε ( ), ε >, =, ; ε,

( ) ( ) ( ) ( ),; ε =. lim, ; ε =,, ε bondedly, in 1 L norm. How does depend on iniial? and se Define ( ) ( ) U = y dy, = sp U 1. For any o solions, υ wih iniial vales, υ,

(, ) (, ) υ 1 υ 1 Proof. Inegrae viscos eqaion: and + 1 =, U U εu + 1 =. V V εv Sbrac hese eqaions, and denoe U V as D: D + 1 U + V D= ε D ( ) By maimm principle,

as claimed. ( n) D(, ) sp D, Le be any seqence of iniial dae ha end weakly, ha is in he sense of he nom, o ; we saw ha hen 1 for every ends weakly o (, ) ( n) eqaion and ( n ) (, ) >. saisfy he basic differenial

( n) ( ) ( n) ( n) + = 1 1 ( ) + = ends o in he sense of disribion. ) Therefore ( n ends o. disribion. So ( n ) ends o, in he sense of he sense of disribion. Since boh n and end o zero as, he + cons in

n consan is zero. B hen ms end o srongly, ha is in he L 1 norm. Smmary If he iniial vales of ( n) end o he iniial vales of weakly, n (, ) ends o (, ) srongly, for every >. The mapping from iniial vales (,) o heir vales (, ), >, is a compac mapping.