Dispersion relations, stability and linearization

Save this PDF as:
Size: px
Start display at page:

Transcription

1 Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x <, t > 0}, and it satisfies a linear, constant coefficient partial differential equation such as the usual wave or diffusion equation. It happens that these type of equations have special solutions of the form or equivalently, u(x, t) = exp(ikx iωt), (1) u(x, t) = exp(σt + ikx). (2) We typically look for solutions of the first kind (1) when wave-like behavior which oscillates in time is expected, whereas (2) is used to investigate growth or decay in time. Plugging either (1) or (2) into the equation yields an algebraic relationship of the form ω = ω(k) or σ = σ(k), called the dispersion relation. It characterizes the dynamics of spatially oscillating modes of the form exp(ikx). Here are a couple standard examples. For the wave equation u tt = c 2 u xx, we plug in a wavelike solution (1) to get ω 2 exp(ikx iωt) = c 2 k 2 exp(ikx iωt), or ω(k) = ±ck. This means there are traveling wave solutions of the form u = exp(ik(x ± ct)), which we might have guessed from the d Alembert formula. For the diffusion equation u t = Du xx, we use (2) to get σ(k) = Dk 2. Since this is negative for all k, it is expected that solutions which are superpositions of (2) also decay. This is consistent with the fundamental solution representation for the diffusion equation. 1.1 Stability For dispersion relations of the form σ = σ(k) stemming from (2), the sign of the real part of σ indicates whether the solution will grow or decay in time. If the real part of σ(k) is negative for all k values, then any superposition of solutions of the form exp(σt + ikx) will also appear to decay. On the other hand, if the real part of σ(k) is positive for some values of k, then over time some components of a superposition will grow exponentially. The former case is called stable, whereas the latter is unstable. If the maximum of the real part of σ is exactly zero, the situation is called marginally stable. It is more difficult to assess the long term behavior in this case. 1.2 Phase and group velocity of waves For dispersion relations of the form ω(k), a solution of the form (1) can be written ( [ u(x, t) = exp ik x ω(k) ]) k t, (3) which we notice are waves traveling at speed ω(k)/k; this is known as the phase velocity. If the phase velocity is different for each k, a superposition of many different waves will appear to spread out or disperse. Surprisingly, if the superposition contains only wavenumbers near a central wavenumber k 0, the wave does not appear to move at the phase velocity σ(k 0 )/k 0, but at a different speed. To see how this happens, consider an initial condition which is a superposition of many different modes u(x, 0) = 1 A(k)e ikx dk,

2 where we think of A(k) as the amount of wavenumber k. If u solves an equation with a given dispersion relation ω = ω(k) then the complete solution can be written u(x, t) = A(k)e ikx iω(k)t dk. (4) In general, the integral oscillates considerably as k is varied, and so one expects significant cancellation. Suppose that A(k) is concentrated about some wavenumber k 0 ; for example, A(k) = exp( (k k 0 ) 2 /ɛ) where ɛ is small. One might think that A(k) could be replaced with a delta function δ(k k 0 ) as a good approximation, and then one recovers u exp(ik 0 x iω(k 0 )t). This turns out to be too simple, since the sideband wavenumbers near k 0 play a role as t becomes large. Instead, we Taylor expand ω(k) ω(k 0 ) + ω (k 0 )(k k 0 ), and plug into (4), giving u(x, t) e it[ω (k 0 )k 0 ω(k 0 )] A(k)e ik(x ω (k 0 )t) dk. The first factor just oscillates in time, whereas the integral is a traveling wave of the form f(x ω (k 0 )t). This means that the entire superposition appears to travel at the group velocity ω (k 0 ), which is in general different than the phase velocity ω(k 0 )/k 0. Only in the case where the dispersion relation is linear ω = ck do the phase and group velocities coincide. Take as an example the Schrödinger equation iu t + u xx = 0. Looking for solutions of the form u = exp(ikx iωt), one gets exp(ikx iωt)[i( iω) k 2 ] = 0, therefore ω = k 2. The phase velocity is ω(k)/k) = k, so a superposition of waves will disperse. The group velocity is ω (k) = 2k, which is twice the phase velocity. 2 Steady state solutions Suppose we have a (possibly nonlinear) PDE of the form u t = R(u, u x, u xx,...) (5) A steady state (or equilibrium) solution u 0 is one for which (u 0 ) t 0, so that it solves R(u 0, (u 0 ) x,...) = 0. (6) In addition, (5) and (6) might also satisfy boundary conditions. Note that a steady state solution solves an ordinary differential equation, rather then a PDE. Example. Often, the solution to (6) is just a constant u = u 0 in space as well as time. For the diffusion equation with Dirichlet-type boundary conditions u t = u xx, u(0, t) = 2 = u(1, t), (7) 2

3 it is easy to see that u(x, t) = 2 is a solution which does not depend on time or the space variable. In general, however, equilibria may depend on x; for example, for the diffusion equation with mixed boundary conditions u t = u xx, u(0, t) = 0, u x (1, t) = 1, (8) the equilibrium solution solves a two-point boundary value problem whose solution is easily obtained as u 0 = x. (u 0 ) xx = 0, u 0 (0) = 0, (u 0 ) x (1) = 1, (9) Example. Fisher s equation is a nonlinear diffusion equation u t = u xx + u(1 u), < x <. (10) We can easily find two constant solutions u(x, t) = u 0. They solve u 0 (1 u 0 ) = 0 so that u 0 = 0, 1. This is one hallmark of nonlinear equations: they often possess numerous steady state solutions. Example. A similar nonlinear diffusion equation is the Allen-Cahn equation u t = u xx + 2u(1 u 2 ), < x <. (11) Constant solutions solve u 0 (1 u 2 0 ) = 0 so that u 0 = 0, ±1. We can also find non-constant steady state solutions to (11) by imposing the conditions lim u(x) = 1, lim x u(x) = 1. x A steady solution u(x, t) = u(x) solves u xx + 2u(1 u 2 ) = 0. This is a so-called Hamiltonian equation, and can be turned into a first order equation by multiplying by u x and integrating. The result is u xx u x + 2u(1 u 2 )u x dx = 1 2 u2 x + u u4 + C = 0, where we used the facts that u xx u x = 1 2 (u2 x) x and f (u)u x = f(u) x. Imposing the conditions at x = ± gives C = 1/2. The first order equation can now be written u x = u 4 2u = 1 u 2, which can be solved by separating variables du 1 u 2 = dx, therefore 1 2 ln 1 + u 1 u = x + c so that finally u(x) = tanh(x + c). Since c is arbitrary, there are an infinite number of steady state solutions which are just translations of one another. 3

4 Example. It was observed long ago that shallow water waves can appear as isolated pulses traveling at speed V. If one runs along at the same speed, the wave appears steady. This can be modeled by the Korteweg-de Vries (KdV) equation Steady solutions u(x, t) = u(x) solve u t V u x + 6uu x + u xxx = 0. V u x + 6uu x + u xxx = 0. If we suppose lim x ± u(x) = 0, this equation can be directly integrated to give V u + 3u 2 + u xx = 0. This is also a Hamiltonian equation, so we turn it into a first order equation by the same trick as in the previous example, giving 1 2 u2 x V 2 u2 + u 3 = 0 where the constant of integration was zero using the far field conditions. This equation can now be solved by separation of variables, which gives u(x) = V 2 sech2 ( V 2 (x + c) where again c allows for translation. These are the famous KdV solitons. ), 3 Linearization It is frequently useful to approximate a nonlinear equation with a linear one, since we know a lot more about linear equations. If u 0 (x) is a steady state of an equation of the form (5), then we can look for solutions of the form u(x, t) = u 0 (x) + ɛw(x, t) (12) where ɛ presumed to be a small parameter. Plugging into (5) and keeping only the terms of order ɛ always gives us a linear, time dependent equation for w. This equation is called the linearization of (5) about u 0. The linearization can be used in a variety of ways. If it happens to admit a dispersion relation, this can tell us about stability or the dispersive nature of the wave components. Linearization also provides a basis for a wide variety of nonlinear techniques, which are beyond the scope of this discussion. It should be acknowledged that linearization has its limits: if w(x, t) becomes large, then the justification for ignoring nonlinear terms proportional to higher powers of ɛ breaks down. Consider again Fisher s equation (10), where we had previously found steady state solutions u 0 = 0, 1. Inserting the perturbed solution (12) into (10) for u 0 = 0, one gets ɛw t = ɛw xx + ɛw ɛ 2 w 2. Keeping only terms of order ɛ, we get the linearized Fisher equation w t = w xx + w. (13) 4

5 This is a diffusion equation with a linear source term. Substituting (2) into it, one gets a dispersion relation σ = k Since σ > 0 when k < 1, we might expect that an arbitrary initial condition for w has components that both grow and decay exponentially. We would therefore say that u 0 = 0 is linearly unstable. If we linearize about u 0 = 1, we get a different result. Setting u = 1 + ɛw gives so that the linearization is now ɛw t = ɛw xx ɛw ɛ 2 w 2, w t = w xx w. (14) This equation has a dispersion relation σ = k 2 1, which is always negative. Thus any initial condition of the form u(x, 0) = 1 + ɛw(x, 0) in the original equation should evolve in time in a way that has w decay and u 1 as t. We call this situation linearly stable. 3.1 Examples of linearization Example. A model for flame-front propagation is the Kuramoto-Sivashinsky equation u t = u xxxx u xx u2 x. It is easy to see that u 0 = 0 is a steady state solution. Letting u = u 0 + ɛw, we get Therefore the linearization is ɛw t = ɛw xxxx ɛw xx + ɛ w2 x. w t = w xxxx w xx. Looking for a dispersion relation of the form w = exp(σt + ikx), we get σ(k) = k 4 + k 2. Since σ > 0 for k < 1, u = 0 is an unstable equilibrium. Any small initial condition containing wavenumbers in this range will therefore always grow. Of course, the solution will eventually get so large that the approximation leading to the linearization is not valid. At this point, more advanced techniques are needed to investigate the behavior. Example. A thin liquid film of height h(x, t) evolves according to the equation h t = ( h 3 [ h xx + Ah 3 ] x )x, (15) where A describes the role of intermolecular forces (e.g. van der Waals forces). Any constant solution h(x, t) = h 0 is a steady state. Setting h(x, t) = h 0 + ɛw, we can use the Taylor series approximations (h 0 + ɛw) 3 = h ɛ3h 2 0w + O(ɛ 2 ), (h 0 + ɛw) 3 = h 3 0 ɛ3h 4 0 w + O(ɛ2 ), where O(ɛ 2 ) refers to terms which are smaller than a constant times ɛ 2. Then plugging h = h 0 +ɛw into equation (15), we get ɛw t = ( (h ɛ3h 2 0w)[ ɛw xx + h 3 0 ɛ3ah 4 0 w] ) x x + O(ɛ2 ), 5

6 so that retaining the ɛ size terms, w t = h 3 0( w xxxx 3Ah 4 0 w xx). The corresponding dispersion relation is found from w = exp(σt + ikx), giving σ(k) = h 3 0( k 4 + 3Ah 4 0 k2 ), which has a band of unstable wavenumbers k < h 2 0 3A, provided A > 0. On the other hand, if A < 0, the film is stable. Example. The so-called sine-gordon equation is u tt = c 2 u xx sin(u). If we set u = ɛw(x, t), we use sin(ɛw) = ɛw + O(ɛ 2 ) which produces the linearized equation w tt = c 2 w xx w. This is a wave type equation, so it makes sense to use the form w(x, t) = exp(ikx iωt) for the dispersion relation. Plugging in we get ω 2 = c 2 k 2 1, ω(k) = ± 1 + c 2 k 2. Since the phase velocity ± 1 + c 2 k 2 /k is not constant, waves disperse. On the other hand, they do not grow or decay in time. 6

Dispersion relations, linearization and linearized dynamics in PDE models

Dispersion relations, linearization and linearized dynamics in PDE models 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient

Travelling waves. Chapter 8. 1 Introduction

Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part

Solitons : An Introduction

Solitons : An Introduction p. 1/2 Solitons : An Introduction Amit Goyal Department of Physics Panjab University Chandigarh Solitons : An Introduction p. 2/2 Contents Introduction History Applications Solitons

Math 575-Lecture 26. KdV equation. Derivation of KdV

Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not

Introduction and some preliminaries

1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.

Solutions of differential equations using transforms

Solutions of differential equations using transforms Process: Take transform of equation and boundary/initial conditions in one variable. Derivatives are turned into multiplication operators. Solve (hopefully

MACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS

. MACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS Willy Hereman Mathematics Department and Center for the Mathematical Sciences University of Wisconsin at

Wave Phenomena Physics 15c. Lecture 11 Dispersion

Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed

Making Waves in Vector Calculus

Making Waves in Vector Calculus J. B. Thoo Yuba College 2014 MAA MathFest, Portland, OR This presentation was produced

A Search for the Simplest Chaotic Partial Differential Equation

A Search for the Simplest Chaotic Partial Differential Equation C. Brummitt University of Wisconsin-Madison, Department of Physics cbrummitt@wisc.edu J. C. Sprott University of Wisconsin-Madison, Department

Rogue periodic waves for mkdv and NLS equations

Rogue periodic waves for mkdv and NLS equations Jinbing Chen and Dmitry Pelinovsky Department of Mathematics, McMaster University, Hamilton, Ontario, Canada http://dmpeli.math.mcmaster.ca AMS Sectional

Lecture17: Generalized Solitary Waves

Lecture17: Generalized Solitary Waves Lecturer: Roger Grimshaw. Write-up: Andrew Stewart and Yiping Ma June 24, 2009 We have seen that solitary waves, either with a pulse -like profile or as the envelope

What we do understand by the notion of soliton or rather solitary wave is a wave prole which is very stable in the following sense

Introduction to Waves and Solitons What we do understand by the notion of soliton or rather solitary wave is a wave prole which is very stable in the following sense i) It is localized, which means it

The Orchestra of Partial Differential Equations. Adam Larios

The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline

Superposition of electromagnetic waves

Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many

Introduction LECTURE 1

LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in

Strauss PDEs 2e: Section Exercise 4 Page 1 of 6

Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where

Dispersive nonlinear partial differential equations

Dispersive nonlinear partial differential equations Elena Kartashova 16.02.2007 Elena Kartashova (RISC) Dispersive nonlinear PDEs 16.02.2007 1 / 25 Mathematical Classification of PDEs based on the FORM

ECE 487 Lecture 6 : Time-Dependent Quantum Mechanics I Class Outline:

ECE 487 Lecture 6 : Time-Dependent Quantum Mechanics I Class Outline: Time-Dependent Schrödinger Equation Solutions to thetime-dependent Schrödinger Equation Expansion of Energy Eigenstates Things you

2. As we shall see, we choose to write in terms of σ x because ( X ) 2 = σ 2 x.

Section 5.1 Simple One-Dimensional Problems: The Free Particle Page 9 The Free Particle Gaussian Wave Packets The Gaussian wave packet initial state is one of the few states for which both the { x } and

1 Assignment 1: Nonlinear dynamics (due September

Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.

Pattern formation in Nikolaevskiy s equation

Stephen Cox School of Mathematical Sciences, University of Nottingham Differential Equations and Applications Seminar 2007 with Paul Matthews, Nottingham Outline What is Nikolaevskiy s equation? Outline

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian

Lecture 18: Bistable Fronts PHYS 221A, Spring 2017

Lecture 18: Bistable Fronts PHYS 221A, Spring 2017 Lectures: P. H. Diamond Notes: Xiang Fan June 15, 2017 1 Introduction In the previous lectures, we learned about Turing Patterns. Turing Instability is

13.1 Ion Acoustic Soliton and Shock Wave

13 Nonlinear Waves In linear theory, the wave amplitude is assumed to be sufficiently small to ignore contributions of terms of second order and higher (ie, nonlinear terms) in wave amplitude In such a

Diagonalization of the Coupled-Mode System.

Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute

PDEs, 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

Final Exam May 4, 2016

1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.

Dispersion in Shallow Water

Seattle University in collaboration with Harvey Segur University of Colorado at Boulder David George U.S. Geological Survey Diane Henderson Penn State University Outline I. Experiments II. St. Venant equations

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem

Stability and Shoaling in the Serre Equations. John D. Carter. March 23, Joint work with Rodrigo Cienfuegos.

March 23, 2009 Joint work with Rodrigo Cienfuegos. Outline The Serre equations I. Derivation II. Properties III. Solutions IV. Solution stability V. Wave shoaling Derivation of the Serre Equations Derivation

Solving Nonlinear Wave Equations and Lattices with Mathematica. Willy Hereman

Solving Nonlinear Wave Equations and Lattices with Mathematica Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs home/whereman/

E = 1 c. where ρ is the charge density. The last equality means we can solve for φ in the usual way using Coulomb s Law:

Photons In this section, we want to quantize the vector potential. We will work in the Coulomb gauge, where A = 0. In gaussian units, the fields are defined as follows, E = 1 c A t φ, B = A, where A and

Physics 486 Discussion 5 Piecewise Potentials

Physics 486 Discussion 5 Piecewise Potentials Problem 1 : Infinite Potential Well Checkpoints 1 Consider the infinite well potential V(x) = 0 for 0 < x < 1 elsewhere. (a) First, think classically. Potential

Periodic Solutions of the Serre Equations. John D. Carter. October 24, Joint work with Rodrigo Cienfuegos.

October 24, 2009 Joint work with Rodrigo Cienfuegos. Outline I. Physical system and governing equations II. The Serre equations A. Derivation B. Justification C. Properties D. Solutions E. Stability Physical

Srednicki Chapter 9. QFT Problems & Solutions. A. George. August 21, Srednicki 9.1. State and justify the symmetry factors in figure 9.

Srednicki Chapter 9 QFT Problems & Solutions A. George August 2, 22 Srednicki 9.. State and justify the symmetry factors in figure 9.3 Swapping the sources is the same thing as swapping the ends of the

ENGI 9420 Lecture Notes 8 - PDEs Page 8.01

ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives

Generalized bilinear differential equations

Generalized bilinear differential equations Wen-Xiu Ma Department of Mathematics and Statistics, University of South Florida, Tampa, FL 33620-5700, USA Abstract We introduce a kind of bilinear differential

Partial Differential Equations

Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but

MATH 425, HOMEWORK 3 SOLUTIONS

MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the

Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche

Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second

Stability and instability of nonlinear waves:

Stability and instability of nonlinear waves: Introduction 1. Nonlinear Waves 2. Stability problems 3. Stability of pulses and fronts 4. Stability of periodic waves 1 Nonlinear waves particular solutions

Introduction to Partial Differential Equations

Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define

13 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,

Introduction of Partial Differential Equations and Boundary Value Problems

Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions

Separation of Variables

Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us

Partial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators

PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane

Summary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:

The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free theory:

LINEAR DISPERSIVE WAVES

LINEAR DISPERSIVE WAVES Nathaniel Egwu December 16, 2009 Outline 1 INTRODUCTION Dispersion Relations Definition of Dispersive Waves 2 SOLUTION AND ASYMPTOTIC ANALYSIS The Beam Equation General Solution

d 3 k In the same non-relativistic normalization x k = exp(ikk),

PHY 396 K. Solutions for homework set #3. Problem 1a: The Hamiltonian 7.1 of a free relativistic particle and hence the evolution operator exp itĥ are functions of the momentum operator ˆp, so they diagonalize

1 Infinite-Dimensional Vector Spaces

Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators

= k, (2) p = h λ. x o = f1/2 o a. +vt (4)

Traveling Functions, Traveling Waves, and the Uncertainty Principle R.M. Suter Department of Physics, Carnegie Mellon University Experimental observations have indicated that all quanta have a wave-like

A Propagating Wave Packet The Group Velocity

Lecture 7 A Propagating Wave Packet The Group Velocity Phys 375 Overview and Motivation: Last time we looked at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds

Long-time solutions of the Ostrovsky equation

Long-time solutions of the Ostrovsky equation Roger Grimshaw Centre for Nonlinear Mathematics and Applications, Department of Mathematical Sciences, Loughborough University, U.K. Karl Helfrich Woods Hole

New explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation

Physics Letters A 07 (00) 107 11 www.elsevier.com/locate/pla New explicit solitary wave solutions for ( + 1)-dimensional Boussinesq equation and ( + 1)-dimensional KP equation Yong Chen, Zhenya Yan, Honging

Hamiltonian partial differential equations and Painlevé transcendents

Winter School on PDEs St Etienne de Tinée February 2-6, 2015 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN SISSA, Trieste Cauchy problem for evolutionary PDEs with

Notes on Fourier Series and Integrals Fourier Series

Notes on Fourier Series and Integrals Fourier Series et f(x) be a piecewise linear function on [, ] (This means that f(x) may possess a finite number of finite discontinuities on the interval). Then f(x)

A review of stability and dynamical behaviors of differential equations:

A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that

Wave interactions. J. Vanneste. School of Mathematics University of Edinburgh, UK. vanneste < > - + Wave interactions p.

Wave interactions J. Vanneste School of Mathematics University of Edinburgh, UK www.maths.ed.ac.uk/ vanneste Wave interactions p.1/50 Linear waves in conservative systems Small amplitude motion: linearize

Are Solitary Waves Color Blind to Noise?

Are Solitary Waves Color Blind to Noise? Dr. Russell Herman Department of Mathematics & Statistics, UNCW March 29, 2008 Outline of Talk 1 Solitary Waves and Solitons 2 White Noise and Colored Noise? 3

Slow Modulation & Large-Time Dynamics Near Periodic Waves

Slow Modulation & Large-Time Dynamics Near Periodic Waves Miguel Rodrigues IRMAR Université Rennes 1 France SIAG-APDE Prize Lecture Jointly with Mathew Johnson (Kansas), Pascal Noble (INSA Toulouse), Kevin

A new integrable system: The interacting soliton of the BO

Phys. Lett., A 204, p.336-342, 1995 A new integrable system: The interacting soliton of the BO Benno Fuchssteiner and Thorsten Schulze Automath Institute University of Paderborn Paderborn & Germany Abstract

6 Non-homogeneous Heat Problems

6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems

Making Waves in Multivariable Calculus

Making Waves in Multivariable Calculus J. B. Thoo Yuba College 2014 CMC3 Fall Conference, Monterey, Ca This presentation

The Euler Equation of Gas-Dynamics

The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u

Chapter 1. Harmonic Oscillator. 1.1 Energy Analysis

Chapter 1 Harmonic Oscillator Figure 1.1 illustrates the prototypical harmonic oscillator, the mass-spring system. A mass is attached to one end of a spring. The other end of the spring is attached to

The Whitham Equation. John D. Carter April 2, Based upon work supported by the NSF under grant DMS

April 2, 2015 Based upon work supported by the NSF under grant DMS-1107476. Collaborators Harvey Segur, University of Colorado at Boulder Diane Henderson, Penn State University David George, USGS Vancouver

Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory

Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Gennady El 1, Roger Grimshaw 1 and Noel Smyth 2 1 Loughborough University, UK, 2 University of Edinburgh, UK

Hamiltonian partial differential equations and Painlevé transcendents

The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN

Optical Solitons. Lisa Larrimore Physics 116

Lisa Larrimore Physics 116 Optical Solitons An optical soliton is a pulse that travels without distortion due to dispersion or other effects. They are a nonlinear phenomenon caused by self-phase modulation

Problem Set 1 October 30, 2017

1. e π can be calculated from e x = x n. But that s not a good approximation method. The n! reason is that π is not small compared to 1. If you want O(0.1) accuracy, then the last term you need to include

Quantization of scalar fields

Quantization of scalar fields March 8, 06 We have introduced several distinct types of fields, with actions that give their field equations. These include scalar fields, S α ϕ α ϕ m ϕ d 4 x and complex

Harmonic oscillator. U(x) = 1 2 bx2

Harmonic oscillator The harmonic oscillator is a familiar problem from classical mechanics. The situation is described by a force which depends linearly on distance as happens with the restoring force

Fourier Series and Integrals

Fourier Series and Integrals Fourier Series et f(x) beapiece-wiselinearfunctionon[, ] (Thismeansthatf(x) maypossessa finite number of finite discontinuities on the interval). Then f(x) canbeexpandedina

Model Equation, Stability and Dynamics for Wavepacket Solitary Waves

p. 1/1 Model Equation, Stability and Dynamics for Wavepacket Solitary Waves Paul Milewski Mathematics, UW-Madison Collaborator: Ben Akers, PhD student p. 2/1 Summary Localized solitary waves exist in the

2. Examples of Integrable Equations

Integrable equations A.V.Mikhailov and V.V.Sokolov 1. Introduction 2. Examples of Integrable Equations 3. Examples of Lax pairs 4. Structure of Lax pairs 5. Local Symmetries, conservation laws and the

MATH 425, FINAL EXAM SOLUTIONS

MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u

INTRODUCTION TO PDEs

INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial

STUDENT ARTICLE 4 Soliton Solutions of Integrable Systems and Hirota s Method

STUDENT ARTICLE 4 Soliton Solutions of Integrable Systems and Hirota s Method Abstract Justin M. Curry Massachusetts Institute of Technology 8 Cambridge, MA 139 jmcurry@mit.edu In this paper we investigate

FAST COMMUNICATION THREE-DIMENSIONAL LOCALIZED SOLITARY GRAVITY-CAPILLARY WAVES

COMM. MATH. SCI. Vol. 3, No., pp. 89 99 c 5 International Press FAST COMMUNICATION THREE-DIMENSIONAL LOCALIZED SOLITARY GRAVITY-CAPILLARY WAVES PAUL A. MILEWSKI Abstract. In a weakly nonlinear model equation

Introduction to the Wave Equation(s)

Chapter 1 Introction to the Wave Equation(s) 1.1 First Order Linear Wave Equation First, x represents space and t represents time. Consider on an infinite domain ( 1), the linear first order wave equation

THE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))

THE WAVE EQUATION The aim is to derive a mathematical model that describes small vibrations of a tightly stretched flexible string for the one-dimensional case, or of a tightly stretched membrane for the

Lecture 21: The one dimensional Wave Equation: D Alembert s Solution

Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without explicit written permission from the copyright owner. 1 Lecture 21: The one dimensional

Finite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations

Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO

MAGIC058 & MATH64062: Partial Differential Equations 1

MAGIC58 & MATH6462: Partial Differential Equations Section 6 Fourier transforms 6. The Fourier integral formula We have seen from section 4 that, if a function f(x) satisfies the Dirichlet conditions,

Solving the Generalized Kaup Kupershmidt Equation

Adv Studies Theor Phys, Vol 6, 2012, no 18, 879-885 Solving the Generalized Kaup Kupershmidt Equation Alvaro H Salas Department of Mathematics Universidad de Caldas, Manizales, Colombia Universidad Nacional

Salmon: 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

Vector Solitons and Their Internal Oscillations in Birefringent Nonlinear Optical Fibers

Vector Solitons and Their Internal Oscillations in Birefringent Nonlinear Optical Fibers By Jianke Yang In this article, the vector solitons in birefringent nonlinear optical fibers are studied first.

Final: 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

Spectral stability of periodic waves in dispersive models

in dispersive models Collaborators: Th. Gallay, E. Lombardi T. Kapitula, A. Scheel in dispersive models One-dimensional nonlinear waves Standing and travelling waves u(x ct) with c = 0 or c 0 periodic

MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20

MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20 MA 201 (2016), PDE 2 / 20 Vibrating string and the wave equation Consider a stretched string of length

Yair Zarmi Physics Department & Jacob Blaustein Institutes for Desert Research Ben-Gurion University of the Negev Midreshet Ben-Gurion, Israel

PERTURBED NONLINEAR EVOLUTION EQUATIONS AND ASYMPTOTIC INTEGRABILITY Yair Zarmi Physics Department & Jacob Blaustein Institutes for Desert Research Ben-Gurion University of the Negev Midreshet Ben-Gurion,

Numerical methods for complex systems. S. V. Gurevich

Numerical methods for complex systems S. V. Gurevich December 16, 2008 2 Part I Partial Differential Equations 3 Chapter 1 Intorduction 1.1 Definition, Notation and Classification A differential equation

Lecture 4: Oscillators to Waves

Matthew Schwartz Lecture 4: Oscillators to Waves Review two masses Last time we studied how two coupled masses on springs move If we take κ=k for simplicity, the two normal modes correspond to ( ) ( )

B = 0. E = 1 c. E = 4πρ

Photons In this section, we will treat the electromagnetic field quantum mechanically. We start by recording the Maxwell equations. As usual, we expect these equations to hold both classically and quantum

Strauss PDEs 2e: Section Exercise 1 Page 1 of 6

Strauss PDEs 2e: Section 3 - Exercise Page of 6 Exercise Carefully derive the equation of a string in a medium in which the resistance is proportional to the velocity Solution There are two ways (among