Survey

* Your assessment is very important for improving the work of artificial intelligence, which forms the content of this project

Document related concepts

Transcript

II Nonlinear wave equations 2.1 Introduction • • • • Introduction Solitary waves Korteweg-deVries (KdV) equation Nonlinear Schrodinger equation Introduction Linear wave equations • Simplest (second order) linear wave equation utt – c2uxx = 0 • D’Alembert’s solution u(x,t) = f(x-ct) + g(x+ct) f, g arbitrary functions Dispersionless Dissipationless Dispersion relation w = ck Introduction Linear wave equations • Simplest Linear ut – cux = 0 or ut + cux = 0 u(x,t) = f(x+ct) or u(x,t) = f(x-ct) • Simplest Dispersive, Dissipationless ut + cux + auxxx = 0 u(x,t) = exp[i(kx – wt)] w = ck - ak3 • Simplest Nondispersive, Dissipative ut + cux - auxx = 0 u(x,t) = exp[i(kx – wt)] w = ck – iak2 Introduction Nonlinear wave equations • Simplest Nonlinear ut + (1+u)ux = 0 u(x,t) = f(x-(1+u)t) Sharpens at leading and trailing edges (shock formation) • Korteweg deVries (KdV) Equation (1895) ut + (1+u)ux + uxxx = 0 c 2 c ( x - ct - x o u(x, t) sech 2 2 Solitary wave/soliton behaviour Dispersion and tendency to shock formation in balance 2.2 Solitary waves Over one hundred and fifty years ago, while conducting experiments to determine the most efficient design for canal boats, a young Scottish engineer named John Scott Russell (1808-1882) made a remarkable scientific discovery. Here is an extract from John Scott Russell’s ‘Report on waves’ Solitary waves Russell’s report on waves “I was observing the motion of a boat which was rapidly drawn along a narrow channel by a pair of horses, when the boat suddenly stopped - not so the mass of water in the channel which it had put in motion; it accumulated round the prow of the vessel in a state of violent agitation, then suddenly leaving it behind, rolled forward with great velocity, assuming the form of a large solitary elevation, a rounded, smooth and well-defined heap of water, which continued its course along the channel apparently without change of form or diminution of speed. I followed it on horseback, and overtook it still rolling on at a rate of some eight or nine miles an hour, preserving its original figure some thirty feet long and a foot to a foot and a half in height. Its height gradually diminished, and after a chase of one or two miles I lost it in the windings of the channel. Such, in the month of August 1834, was my first chance interview with that singular and beautiful phenomenon which I have called the Wave of Translation”. 2.3 Korteweg deVries (KdV) equation • The wave of translation (or solitary wave) observed by John Scott Russell is described by a nonlinear wave equation known as the Korteweg-deVries (KdV) equation. u t - 6uu x u xxx 0 • We review various possible types of nonlinearity in wave equations before studying two specific equations – the KdV and the nonlinear Schrodinger (NLS) equations. Korteweg deVries (KdV) equation Numerical solution (strong dispersive term) 'wave.dat' 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0 -0.2 0 2 20 40 1.5 60 80 1 0.5 100 120 1400 3 2.5 4 3.5 Korteweg deVries (KdV) equation Numerical solution (weak dispersive term) 'wave.dat' 2 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0 -0.2 0 2 20 40 1.5 60 80 1 0.5 100 120 1400 3 2.5 4 3.5 Korteweg deVries (KdV) equation Effect of nonlinear term ut = -(1+u)ux 1.5 The sequence of plots at t = 0, Dt and 2Dt illustrate how a pulse forms and splits off from the leading edge of a smooth front. 1 u(Dt) -(1+u(Dt))ux(Dt) 0.5 2 4 6 1 10 -0.5 ux(Dt) -1 1.5 8 1.5 u(0) 1 u(2Dt) -(1+u(0))ux(0) 0.5 0.5 2 -0.5 -1 4 6 ux(0) 8 10 2 -0.5 -1 4 6 8 10 Korteweg deVries (KdV) equation Effect of dispersive term ut = - uxxx 1.5 1 u(0) -uxxx(0) 0.5 2 4 6 8 10 Combined effects of nonlinear and dispersive terms -0.5 -1 1.5 1 1.5 u(0) u(2Dt) u(Dt) 1 0.5 0.5 2 4 6 8 10 2 -0.5 -0.5 -1 -1 4 6 8 10 Korteweg deVries (KdV) equation Soliton simulations These simulations come from Klaus Brauer's webpage (Osnabrück) Korteweg deVries (KdV) equation Solution for PBC and sinusoidal initial conditions This animation by K. Takasaki shows the sinusoidal initial state breaking up into a soliton train. Zabusky and Kruskal (1966). Korteweg deVries (KdV) equation Analytic solution •KdV equation u t 6uux u xxx 0 •Let the solution be u = u(x,t) and consider a change of variables x = x – ct and t = t •Call the function in new variables f(x,t) •The change in u or f brought about by translations (dx, dt) or (x, t) is du u dx u dt x t df f dx f dt x t dx dx - c dt dt dt Korteweg deVries (KdV) equation Analytic solution • If we convert the change in f brought about by translations through (dx, dt) into changes in f brought about by translations through (dx, dt) df f dx - c dt f dt x t df f dx f -c f dt x t x •Since u and f represent the same function the same translation (dx, dt) must produce the same change in either. Hence u f x x u f -c f t t x Korteweg deVries (KdV) equation Analytic solution • When transforming the pde from (x, t) to (x, t) we must make the replacements x x -c t t x • In the (x, t) variables a soliton moves along the x axis as time advances •In the (x, t) variables a soliton is stationary in time provided we choose c in the transformation to be the soliton velocity Korteweg deVries (KdV) equation Analytic solution • The conventional form for the KdV equation is u t - 6uu x u xxx 0 • Travelling wave solutions have the form u(x, t) f(x ) x x - ct c is the wave velocity • Substituting for u in the KdV equation and setting the time derivative to zero we obtain - cf x - 6ff x f xxx 0 Korteweg deVries (KdV) equation Analytic solution •Integrate twice wrt x dx - cf x - 6ff x f xxx 0 - cf - 3f 2 f xx A 2 d x f cf 3f f xx dxf x A x c 2 1 2 3 - f - f f x Af B 2 2 1 2 3 c 2 f x f f Af B 2 2 Korteweg deVries (KdV) equation Analytic solution • A and B are constants of integration. In order to have a localised traveling wave packet, we impose boundary conditions: f, f x , f xx all tend to zero as |x| goes to infinity. • To ensure these conditions we set A = B = 0. Solutions also exist at zeros of the polynomial in f. • The solution with A = B = 0 obeys f x 2 f 2 2f c • Rearrange to df f 2f c dx Korteweg deVries (KdV) equation Analytic solution • Make change of variable c f - sech 2 2 df f 2f c becomes d c sech 2 tanh c - sech 2 c tanh 2 2 xo 2 x c c - 2 d c • Last term on rhs is constant of integration Korteweg deVries (KdV) equation Analytic solution • Rearrange to c x - x o 2 • Make back substitution c 2 c ( x - ct - x o f sech 2 2 2 sech e e d sech - sech tanh d 2.4 Nonlinear Schrödinger equation • The naming of the nonlinear Schrödinger (NLS) equation becomes obvious when it is compared to the time-dependent Schrödinger equation from quantum mechanics iψ P ψxx Q ψ 2 ψ 0 t 2 iψ ψxx V ψ 0 t 2m • The NLS can be derived for wave packets localised in k space for systems where the dispersion relation depends on 2 wave intensity w w (k, ) Nonlinear Schrödinger equation Derivation from dispersion relation • Consider the superposition of 2 waves of similar wavenumber and frequency ψ1 cos[(k Dk)x - (w Dw ) t] ψ 2 cos[(k Dk)x - (w Dw ) t] ψ1 ψ 2 2 cos [Dk x - Dwt] cos[kx - wt] • The result is a slow envelope wave with group velocity vg = Dw/ D k and a rapid carrier wave with velocity w/k • Simulation with Dw/Dk = 1 and w/k = 20 Nonlinear Schrödinger equation Derivation from dispersion relation • The NLS is derived from the dispersion relation for the envelope function which has a slow time variation cf the carrier waves 2 • Suppose that the dispersion relationship is w w (k, ) Make a Taylor expansion of this about ko and zero intensity ω ω - ωo k 1 ω (k k o ) 2! k 2 2 wo ... 2 2 o 2 wo (k k o ) 2 ... Nonlinear Schrödinger equation Derivation from dispersion relation • Let Ωω-ω Kk –k o ω v g k ω , 2 o 2P o o 2ω k 2 ω , 2 o o Q ω ψ2 ω ,ψ • Then the Taylor expanded dispersion relation becomes v g K PK Q 2 2 o 2 o Nonlinear Schrödinger equation Derivation from dispersion relation • Consider a wavepacket constructed from a small group of waves in slow variables X = ex, T = et e <<1 (x, t) (k, w ) exp[i(kx - wt)] dk dw - KX T (K, ) exp[i( )] dK d * exp[i(k o x - w o t)] e e - KX T (X, T) (K, ) exp[i( )] dK d * carrier wa ve e e - • The latter is the envelope function in ‘slow’ variables X,T Nonlinear Schrödinger equation Derivation from dispersion relation KX T (X, T) )] dK d (K, ) exp[i( X X - e e iK KX T )] dK d (K, ) exp[i( e - e e iK (X, T) e iK K ie X e X i ie T e T Nonlinear Schrödinger equation Derivation from dispersion relation • The dispersion relation v g K PK Q 2 2 becomes 2 2 2 ie v g - ie P - ie e Q T X X ie T - iev g X e P XX e Q 2 2 • Make further change of variables x X - v g T t eT e vg X x T t x 2 Nonlinear Schrödinger equation Derivation from dispersion relation ie T - iev g X e P XX e Q 2 2 2 becomes i P t i t Q 0 2 xx q 0 2 xx • This is the conventional form for the NLS equation. It has an envelope solution with a sech profile. (See handout) Nonlinear Schrödinger equation Application to lattice dynamics • Hooke’s Law plus additional nonlinear term 1 2 AK 4 U(r) Kr r 2 4 dU(r) F(r) -Kr AKr 3 dr • Equation of motion mu n K(u n 1 2u n u n 1 ) AK (u n 1 u n ) 3 (u n u n -1 ) 3 • Solution and dispersion relation u n R exp[i(kna - wt)] c.c. K ka 2 2 ka ω2 sin 1 6A R sin m 2 2 Nonlinear Schrödinger equation Application to lattice dynamics • We have just seen that introduction of a nonlinear term in the force law for a 1-D chain of atoms leads to a dispersion relation which depends on |R|2. At the website below, use the monatomic chain applet to see some of these localised modes. •Intrinsic localised modes in lattice dynamics of crystals K ka 2 ka sin 1 6A R sin 2 m 2 2 K ka ω2 sin m 2 ω2 1 0.8 0.6 0.4 0.2 k 0.5 1 1.5 2 2.5 3 Nonlinear Schrödinger equation Application to lattice dynamics • Click on monatomic 1-D chains and then on the link in the title to the page (works best with Internet Explorer) •You will find stationary ILM with • envelope function (c.f. solutions of NLS equation) is composed of groups of waves centred on the Brillouin zone boundary (k = p) (group velocity zero) • moving ILM composed of groups of waves centred away from the Brillouin zone boundary (group velocity nonzero) Nonlinear Schrödinger equation Application to lattice dynamics • You will also find • molecular dynamics simulations showing ILM in 3-D crystals (click on 3-D Ionic crystals) •Simulations showing ILM in 1-D chains of interacting spins Nonlinear Schrödinger equation Application to optical communications • Read the introductory articles on • Solitons in optical communications by Ablowitz et al. • Historical aspects of optical solitons by Hasegawa • Soliton propagation in optical fibres