Continuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman
|
|
- Loraine Long
- 5 years ago
- Views:
Transcription
1 Continuous and Discrete Homotopy Operators with Applications in Integrability Testing Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA home/whereman/ Colloquium Talk Centre de Recherches Mathématiques Montréal, Canada Thursday, December 4, 2003, 15:30 Collaborator: Michael Colagrosso Students: Adam Ringler, Ingo Kabirschke Francis Martin & Kara Namanny Talk dedicated to Ryan Sayers Research supported in part by NSF under Grant CCR
2 OUTLINE Part I: Continuous Case Integration by Parts on the Jet Space (by hand) + Mathematica Experiment Exactness or Integrability Criterion: Continuous Euler Operator Continuous Homotopy Operator Application of Continuous Homotopy Operator Demo of Mathematica software Part II: Discrete Case Inverting the Total Difference Operator (by hand) Exactness or Total Difference Criterion: Discrete Euler Operator Discrete Homotopy Operator Application of Discrete Homotopy Operator Demo of Mathematica Software Future Research 2
3 Problem Statement For continuous case: Given, for example, f = 3 u v 2 sin(u) u 3 sin(u) 6 v v cos(u) + 2 u u cos(u) + 8 v v Find F so that f = D x F or F = f dx. Result: F = 4 v 2 + u 2 cos(u) 3 v 2 cos(u) Can this be done without integration by parts? Can the problem be reduced to a single integral in one variable? For discrete case: Given, for example, f n = +1 v n vn v n+1 + vn v n+2 +1 v n Find F n so that f n = F n = F n+1 F n or F n = 1 f n. Result: F n = vn v n + +1 v n + +2 v n+1. How can this be done algorithmically? Can this be done in the same way as the continuous case? 3
4 Part I: Continuous Case Integration by Parts on the Jet Space Given f involving u(x) and v(x) and their derivatives f = 3 u v 2 sin(u) u 3 sin(u) 6 v v cos(u)+2 u u cos(u)+8 v v Find F so that f = D x F or F = f dx. Integrate by parts (compute F by hand) f F 8 v v 4 v 2 2 u u cos(u) u 2 cos(u) u 3 sin(u) 6 v v cos(u) 3 v 2 cos(u) 3 u v 2 sin(u) Integral: F = 4 v 2 + u 2 cos(u) 3 v 2 cos(u) Remark: For simplicity we denote f(u, u, u,, u (m) ) as f(u). 4
5 Exactness Criterion: Continuous Euler Operator (variational derivative) Definition (exactness): A function f(u) is exact, i.e. can be integrated fully, if there exists a function F (u), such that f(u) = D x F (u) or equivalently F (u) = D 1 x f(u) = f(u) dx. x D x is the (total) derivative with respect to x. Theorem (exactness or integrability test): A necessary and sufficient condition for a function f to be exact, i.e. the derivative of another function, is that L (0) u (f) 0 where L (0) u is the continuous Euler operator (variational derivative) defined by L (0) u = m 0 ( D x ) k k=0 u (k) = u D x u + D2 x u + + ( 1)m 0 D m 0 x where m 0 is the order (of f). u (m 0) Proof: See calculus of variations (derivation of Euler-Lagrange equations the forgotten case!). 5
6 Example: Apply the continuous Euler operator to f(u) = 3u v 2 sin(u) u 3 sin(u) 6v v cos(u)+2u u cos(u)+8v v Here u = (u, v). For component u (order 2): L (0) u (f) = u (f) D x u (f) + D2 x u (f) = 3u v 2 cos(u) u 3 cos(u) + 6v v sin(u) 2u u sin(u) D x [3v 2 sin(u) 3u 2 sin(u) + 2u cos(u)] + D 2 x[2u cos(u)] = 3u v 2 cos(u) u 3 cos(u) + 6v v sin(u) 2u u sin(u) 0 [3u v 2 cos(u) + 6v v sin(u) 3u 3 cos(u) 6u u sin(u) 2u u sin(u) + 2u cos(u)] +[ 2u cos(u) 6u u sin(u) + 2u cos(u)] For component v (order 2): L (0) v (f) = v (f) D x v (f) + D2 x v (f) = 6u v sin(u) 6v cos(u) D x [ 6v cos(u) + 8v ] + D 2 x[8v ] = 6u v sin(u) 6v cos(u) [6u v sin(u) 6v cos(u)+8v ]+8v 0 6
7 Computation of the integral F Definition (higher Euler operators): The continuous higher Euler operators are defined by L (i) u = m i k ( D x ) k i i u (k) k=i These Euler operators all terminate at some maximal order m i. Examples (for component u) : L (0) u = u D x u + D2 x u D3 x u + + ( 1)m 0 D m 0 x u (m 0) L (1) u = u 2D x u +3D2 x u 4D3 x u + (4) ( 1)m 1 m 1 D m 1 1 x L (2) u = u 3D x u +6D2 x u (4) 10D3 x u (5)+ +( 1)m 2 L (3) u = u 4D x u (4)+10D2 x u (5) 20D3 x u + (6) ( 1)m 3 Similar formulae for component L (i) v u (m 1) m 2 D m 2 2 x 2 u (m 2) m 3 D m 3 3 x 3 u (m 3) 7
8 Definition (homotopy operator): The continuous homotopy operator is defined by where H u = 1 0 N f r (u) = p r r=1 f r (u)[λu] dλ λ i=0 D i x[u r L (i+1) u r ] p r is the maximum order of u r in f N is the number of dependent variables f r (u)[λu] means that in f r (u) one replaces u λu, u λu, etc. Example: For a two-component system (N = 2) where u = (u, v): with and H u = 1 0 {f 1(u)[λu] + f 2 (u)[λu]} dλ λ f 1 (u) = p 1 i=0 f 2 (u) = p 2 i=0 D i x[ul (i+1) u ] D i x[vl (i+1) v ] Theorem (integration via homotopy operator): Given an integrable function f F = D 1 x f = f dx = H u (f) Proof: Olver s book Applications of Lie Groups to Differential Equations, p Proof is given in terms of differential forms. Work of Henri Poincaré ( ), George de Rham (1950), and Ian Anderson & Tom Duchamp (1980). Proofs based on calculus: Deconinck and Hereman. 8
9 Example: Apply the continuous homotopy operator to integrate f(u) = 3u v 2 sin(u) u 3 sin(u) 6vv cos(u) + 2u u cos(u) + 8v v For component u (order 2): i L (i+1) u (f) D i x ( ul (i+1) u (f) ) 0 u f 2D x ( u f) =3v 2 sin u 3u 2 sin u + 2u cos u 2D x (2u cos u) =3v 2 sin u + u 2 sin u 2u cos u 3uv 2 sin u + uu 2 sin u 2uu cos u 1 u f =2u cos u D x [2uu cos u] =2u 2 cos u + 2uu cos u 2uu 2 sin u Hence, f 1 (u)(f) = 3uv 2 sin(u) uu 2 sin(u) + 2u 2 cos(u) For component v (order 2): i L (i+1) v (f) D i x[vl (i+1) v (f)] 0 6v cos(u) + 8v 2D x [8v ] = 6v cos(u) 8v 6v 2 cos(u) 8vv 1 8v D x [8vv ] = 8v 2 + 8vv Hence, f 2 (u)(f) = 6v 2 cos(u) + 8v 2 The homotopy operator leads to an integral for (auxiliary) variable λ. (Use standard integration by parts to work the integral). F (u) = H u (f) = 1 0 {f 1(u)(f)[λu] + f 2 (u)(f)[λu]} dλ λ = 1 0 [3λ2 uv 2 sin(λu) λ 2 uu 2 sin(λu) + 2λu 2 cos(λu) 6λv 2 cos(λu) + 8λv 2 ] dλ = 4v 2 + u 2 cos(u) 3v 2 cos(u) 9
10 Application: Conserved densities and fluxes for PDEs with transcendental nonlinearities Definition (conservation law): D t ρ + D x J = 0 conserved density ρ and flux J. (on PDE) Example: Sine-Gordon system (type u t = F) has scaling symmetry In terms of weights: u t v t = v = u xx + α sin(u) (t, x, u, v, α) (λ 1 t, λ 1 x, λ 0 u, λv, λ 2 α) w(d x ) = 1, w(d t ) = 1, w(u) = 0, w(v) = 1, w(α) = 2 Conserved densities and fluxes ρ (1) = 2α cos(u) + v 2 + u x 2 J (1) = 2u x v ρ (2) = u x v J (2) = [ 1 2 v u x 2 α cos(u)] ρ (3) = 12 cos(u)vu x + 2v 3 u x + 2vu x 3 16v x u 2x ρ (4) = 2 cos 2 (u) 2 sin 2 (u) + v 4 + 6v 2 u 2 x + u 4 x + 4 cos(u)v cos(u)u 2 x 16v 2 x 16u 2 2x. are all scaling invariant! Remark: J (3) and J (4) are not shown (too long). 10
11 Algorithm for Conserved Densities and Fluxes Example: Density and flux of rank 2 for sine-gordon system Step 1: Construct the form of the density ρ = αh 1 (u) + h 2 (u)v 2 + h 3 (u)u x 2 + h 4 (u)u x v where h i (u) are unknown functions. Step 2: Determine the functions h i Compute E = D t ρ = ρ t + ρ (u)[f] = ρ t + m1 k=0 (on PDE) ρ D k u xu t + m2 kx k=0 ρ v kx D k xv t Since E = D t ρ = D x J, the expression E must be integrable. Require that L (0) u (E) 0 and L (0) v (E) 0. Solve the system of linear mixed system (algebraic eqs. and ODEs): h 2 (u) h 3 (u) = 0 h 2 (u) = 0 h 3 (u) = 0 h 4 (u) = 0 h 2 (u) = 0 h 4 (u) = 0 2h 2 (u) h3 (u) = 0 2h 2 (u) h3 (u) = 0 h 1 (u) + 2 sin(u)h2 (u) = 0 h 1 (u) + 2 sin(u)h2 (u) + 2 cos(u)h2 (u) = 0 11
12 Solution: h 1 (u) = 2c 1 cos(u) + c 3 h 2 (u) = h 3 (u) = c 1 h 4 (u) = c 2 (with arbitrary constants c i ). Substitute in ρ ρ = c 1 (2α cos(u) + v 2 + u x 2 ) + c 2 (u x v) + c 3 α Step 3: Compute the flux J First, compute E =D t ρ = c 1 ( 2αu t sin u + 2vv t + 2u x u xt ) + c 2 (u xt v + u x v t ) = c 1 ( 2αv sin u + 2v (u 2x + α sin u) + 2u x v x ) +c 2 (v x v + u x (u 2x + α sin u)) = c 1 (2u 2x v + 2u x v x ) + c 2 (vv x + u x u 2x + αu x sin u) Since E = D t ρ = D x J, one must integrate f = E. Apply the homotopy operator for each component of u = (u, v). For component u (order 2): i L (i+1) u (f) D i x ( ul (i+1) u (f) ) 0 2c 1 v x + c 2 (u 2x α sin u) 2c 1 uv x + c 2 (uu 2x αu sin u) 1 2c 1 v c 2 u x 2c 1 (u x v + uv x ) c 2 (u 2 x + uu 2x ) Hence, f 1 (u)(f) = 2c 1 u x v c 2 (u 2 x + αu sin u) 12
13 For component v (order 1): i L (i+1) v (f) D i ( x vl (i+1) v (f) ) 0 2c 1 u x c 2 v 2c 1 u x v c 2 v 2 Hence, f 2 (u)(f) = 2c 1 u x v c 2 v 2 The homotopy operator leads to an integral for (one) variable λ : J(u) = H u (f) = 1 (f 0 1(u)(f)[λu] + f 2 (u)(f)[λu]) dλ λ = ( 1 4c1 λu 0 x v + c 2 (λu 2 x + αu sin(λu) + λv 2 ) ) dλ = 2c 1 u x v c v u2 x α cos u Split the density and flux in independent pieces (for c 1 and c 2 ): ρ (1) = 2α cos u + v 2 + u x 2 J (1) = 2u x v ρ (2) = u x v J (2) = 1 2 v2 1 2 u2 x + α cos u Remark: Computation of J (3) and J (4) requires integration with the homotopy operator! 13
14 Computer Demos (1) Use continuous homotopy operator to integrate f = 3 u v 2 sin(u) u 3 sin(u) 6 v v cos(u)+2 u u cos(u)+8 v v (2) Compute densities of rank 8 and fluxes for 5th-order Korteweg-de Vries equation with three parameters: u t + αu 2 u x + βu x u 2x + γuu 3x + u 5x = 0 (α, β, γ are nonzero constant parameters). (3) Compute density of rank 4 and flux for sine-gordon system: u t v t = v = u xx + α sin(u) 14
15 Analogy PDEs and DDEs Continuous Case (PDEs) Semi-discrete Case (DDEs) System u t = F(u, u x, u 2x,...) =F(..., 1,, +1,...) Conservation Law D t ρ + D x J = 0 ρ n + J n+1 J n = 0 Symmetry D t G = F (u)[g] D t G = F ( )[G] = F(u + ɛg) ɛ ɛ=0 = F(u ɛ n + ɛg) ɛ=0 Recursion Operator D t R + [R, F (u)] = 0 D t R + [R, F ( )] = 0 Table 1: Conservation Laws and Symmetries KdV Equation Volterra Lattice Equation u t = 6uu x + u 3x = (+1 1 ) Densities ρ = u, ρ = u 2 ρ n =, ρ n = ( 1u 2 n+ +1 ) ρ = u u2 x ρ n = 1 3 u3 n+ +1 ( ) Symmetries G=u x, G=6uu x + u 3x G = +1 ( ) G=30u 2 u x + 20u x u 2x 1 ( ) +10uu 3x + u 5x Recursion Operator R = D 2 x + 4u + 2u x D 1 x R = (I + D)( D D 1 ) (D I) 1 1 Table 2: Prototypical Examples 15
16 Part II: Discrete Case Definitions (shift and total difference operators): D is the up-shift (forward or right-shift) operator if for F n DF n = F n+1 = F n n n+1 D 1 the down-shift (backward or left-shift) operator if D 1 F n = F n 1 = F n n n 1 = D I is the total difference operator F n = (D I)F n = F n+1 F n D (up-shift operator) corresponds the differential operator D x D x F (x) F n+1 F n x = F n x (set x = 1) For k > 0, D k = D D D (k times). Similarly, D k = D 1 D 1 D 1. Problem to be solved: Continuous case: Given f. Find F so that f = D x F or F = D 1 x f = f dx. Discrete case: Given f n. Find F n so that f n = F n = F n+1 F n or F n = 1 f n. 16
17 Inverting the Operator Given f n involving and v n and their shifts: f n = +1 v n vn+u 2 n+1 +2 v n+1 +vn+1+u 2 n+3 v n+2 +1 v n Find F n so that f n = F n = F n+1 F n or F n = 1 f n. Invert the operator (compute F n by hand) f n F n v 2 n v 2 n v 2 n+1 +1 v n +1 v n v n+1 +1 v n +1 v n ++2 v n+1 +2 v n+1 +2 v n+1 +3 v n+2 Result: F n = v 2 n + +1 v n + +1 v n + +2 v n+1. Remarks: We denote f(, +1, +2,, +p ) as f( ). Assume that all negative shifts have been removed via up-shifting Replace f n = 2 v n v n+3 by f n = D 2 f n = v n+2 v n+5. 17
18 Total Difference Criterion: Discrete Euler Operator (variational derivative) Definition (exactness): A function f n ( ) is exact, i.e. a total difference, if there exists a function F n ( ), such that f n = F n or equivalently F n = 1 f n. D is the up-shift operator. Theorem (exactness or total difference test): A necessary and sufficient condition for a function f n to be exact, i.e. a total difference, is that L (0) (f n ) 0, where L (0) is the discrete Euler operator (variational derivative) defined by L (0) = m 0 D k k=0 +k = + D 1 ( ) + D 2 ( ) + + D m 0 ( ) m0 = ( m0 D k ) k=0 L (0) = (I + D 1 + D D m 0 ) where m 0 is the highest forward shift (in f n ). 18
19 Example: Apply the discrete Euler operator to f n ( ) = +1 v n v 2 n v n+1 +v 2 n+1++3 v n+2 +1 v n Here = (, v n ). For component (highest shift 3): L (0) (f n ) = [I + D 1 + D 2 + D 3 ](f n ) =[ +1 v n ] + [ 1 v n v n v n 1 ]+[ 1 v n 1 ]+[v n 1 ] 0 For component v n (highest shift 2): L (0) v n (f n ) = [I + D 1 + D 2 ](f n ) v n = [ +1 2v n +1 ] + [ v n ] + [+1 ] 0 19
20 Computation of F n Definition (higher Euler operators): The discrete higher Euler operators are defined by L (i) = ( mi k D k ) i These Euler operators all terminate at some maximal shifts m i. Examples (for component ): L (0) = L (1) = L (2) = L (3) = k=i (I + D 1 + D 2 + D D m 0 ) (D 1 + 2D 2 + 3D 3 + 4D m 1 D m 1 ) (D 2 + 3D 3 + 6D D m 2(m 2 1) D m 2 ) (D 3 + 4D D D 6 m D m 3 ) 3 Similar formulae for L (i) v n. 20
21 Definition (homotopy operator): The discrete homotopy operator is defined by where H un = 1 0 f r,n ( ) = p r N r=1 f r,n ( )[λ ] dλ λ i=0 (D I) i [u r,n L (i+1) u r,n ] p r is the maximum shift of u r,n in f n N is the number of dependent variables f r,n ( )[λ ] means that in f r,n ( ) one replaces λ, +1 λ+1, etc. Example: For a two-component system (N = 2) where = (, v n ): with and H un = 1 0 {f 1,n( )[λ ] + f 2,n ( )[λ ]} dλ λ f 1,n ( ) = p 1 (D I) i [ L (i+1) i=0 ] f 2,n ( ) = p 2 (D I) i [v n L (i+1) i=0 v n ] Theorem (total difference via homotopy operator): Given a function f n which is a total difference, then F n = 1 f n = H un (f n ) Proof: Recent work by Mansfield and Hydon on discrete variational bi-complexes. Proof is given in terms of differential forms. Proof based on calculus: Deconinck and Hereman. 21
22 Higher Euler Operators Side by Side Continuous Case (for component u) L (0) u = u D x + D 2 x D 3 x + u x u 2x u 3x 2D x + 3D 2 x 4D 3 x + u x u 2x u 3x u 4x 3D x + 6D 2 x 10D 3 x + u 2x u 3x u 4x u 5x 4D x + 10D 2 x 20D 3 x + u 3x u 4x u 5x u 6x L (1) u = L (2) u = L (3) u = Discrete Case (for component ) L (0) = L (1) = L (2) = L (3) = (I + D 1 + D 2 + D 3 + ) (D 1 + 2D 2 + 3D 3 + 4D 4 + ) (D 2 + 3D 3 + 6D D 5 + ) (D 3 + 4D D D 6 + ) 22
23 Homotopy Operators Side by Side Continuous Case (for components u and v) with and H u = 1 0 {f 1(u)[λu] + f 2 (u)[λu]} dλ λ f 1 (u) = p 1 i=0 f 2 (u) = p 2 i=0 D i x[ul (i+1) u ] D i x[vl (i+1) v ] Discrete Case (for components and v n ) with and H un = 1 0 {f 1,n( )[λ ] + f 2,n ( )[λ ]} dλ λ f 1,n ( ) = p 1 (D I) i [ L (i+1) i=0 ] f 2,n ( ) = p 2 (D I) i [v n L (i+1) i=0 v n ] 23
24 Example: Apply the discrete homotopy operator to f n ( ) = +1 v n v 2 n v n+1 +v 2 n+1++3 v n+2 +1 v n For component (highest shift 3): i L (i+1) (f n ) (D I) i [ L (i+1) (f n )] 0 1 v n v n +2v n 1 1 v n v n +2 v n v n 1 +3v n 1 +1 v n +3+1 v n 1 v n 1 3 v n 1 2 v n 1 +2 v n+1 +1 v n +1 v n + v n 1 Hence, f 1,n ( )(f n ) = 2 +1 v n + +1 v n + +2 v n+1 For component v n (highest shift 2): i L (i+1) v n (f n ) (D I) i [v n L (i+1) v n (f n )] v n v n +2v 2 n+2+1 v n v n+1 +1 v n Hence, f 2,n ( )(f n ) = +1 v n + 2v 2 n + +1 v n + +2 v n+1 The homotopy operator leads to an integral for variable λ. (Use standard integration by parts to work the integral). F n ( ) = H un (f n ) = 1 {f 0 1,n( )(f n )[λ ] + f 2,n ( )(f n )[λ ]} dλ λ = 1 0 [2λv2 n + 3λ 2 +1 v n + 2λ+1 v n + 2λ+2 v n+1 ] dλ = vn v n + +1 v n + +2 v n+1 24
25 Application: Conserved densities and fluxes for DDEs Definition (conservation law): D t ρ n + J n = D t ρ n + J n+1 J n = 0 conserved density ρ n and flux J n. Example The Toda lattice (type = F) : has scaling symmetry In terms of weights: = v n 1 v n v n = v n ( +1 ) (t,, v n ) (λ 1 t, λ, λ 2 v n ). (on DDE) w( d dt ) = 1, w() = w(+1 ) = 1, w(v n ) = w(v n 1 ) = 2. Conserved densities and fluxes ρ (0) n = ln(v n ) J n (0) ρ (1) n = J n (1) ρ (2) n = 1 2 u2 n + v n J n (2) ρ (3) n = 1 3 u3 n + (v n 1 + v n ) J n (3) are all scaling invariant! = = v n 1 = v n 1 = 1 v n 1 + v 2 n 1 25
26 Algorithm for Conserved Densities and Fluxes Example: Density of rank 3 for Toda system Step 1: Construct the form of the density. ρ n = c 1 u 3 n + c 2 v n 1 + c 3 v n where c i are unknown constants. Step 2: Determine the constants c i. Compute E n = D t ρ n = ρ n t + ρ n( )[F] (on DDE) = (3c 1 c 2 )u 2 nv n 1 + (c 3 3c 1 )u 2 nv n + (c 3 c 2 )v n 1 v n +c 2 1 v n 1 + c 2 v 2 n 1 c 3 +1 v n c 3 v 2 n Compute Ẽn = DE n to remove negative shift n 1. Since Ẽn = J n, the expression Ẽn must be a total difference. Require L (0) (Ẽn) = (I + D 1 + D 2 )(Ẽn) = = 2(3c 1 c 2 ) v n 1 + 2(c 3 3c 1 ) v n +(c 2 c 3 ) 1 v n 1 + (c 2 c 3 )+1 v n 0 (D + I + D 1 )(E n ) and L (0) v n (Ẽn) = (I + D 1 )(Ẽn) = (D + I)(E n ) v n v n = (3c 1 c 2 )u 2 n+1 + (c 3 c 2 )v n+1 + (c 2 c 3 ) +1 +2(c 2 c 3 )v n + (c 3 3c 1 )u 2 n + (c 3 c 2 )v n
27 Solve the linear system S = {3c 1 c 2 = 0, c 3 3c 1 = 0, c 2 c 3 = 0}. Solution: 3c 1 = c 2 = c 3 Choose c 1 = 1 3, and c 2 = c 3 = 1. Substitute in ρ n ρn = 1 3 u3 n + (v n 1 + v n ) Step 3: Compute the flux J n. Start from Ẽn = +1 v n v 2 n v n+1 + v 2 n+1 Apply the discrete homotopy operator to f n = Ẽn. For component (highest shift 2): i L (i+1) (f n ) (D I) i ( L (i+1) ( Ẽn)) 0 1 v n v n 1 v n v n 1 1 v n 1 +1 v n 1 v n 1 Hence, f 1,n ( )(f n ) = 2 +1 v n For component v n (highest shift 1): i L (i+1) v n (f n ) (D I) i (v n L (i+1) v n ( Ẽn)) v n v n +1 +2vn 2 Hence, f 2,n ( )(f n ) = +1 v n + 2v 2 n J n = H un (f n ) = 1 (f 0 1,n( )(f n )[λ ] + f 2,n (f n )( )[λ ]) dλ λ = 1 0 (3λ2 +1 v n + 2λvn) 2 dλ = +1 v n + vn. 2 Final Result: J n = D 1 Jn = 1 v n 1 + v 2 n 1 27
28 Computer Demos (1) Use discrete homotopy operator to compute F n = 1 f n for f n = +1 v n vn+u 2 n+1 +2 v n+1 +vn+1+u 2 n+3 v n+2 +1 v n (2) Compute density of rank 4 and flux for Toda system: = v n 1 v n v n = v n ( +1 ) (3) Compute density of rank 2 for Ablowitz-Ladik system: i = κu n ( ) (u n is the complex conjugate of ). This is an integrable discretization of the NLS equation: iu t + u xx + κu 2 u = 0 Take equation and its complex conjugate. Treat and v n = u n as dependent variables. Absorb i in t: = v n ( ) v n = (v n+1 2v n + v n 1 ) v n (v n+1 + v n 1 ). 28
29 Future Research Generalize continuous homotopy operator in multi-dimensions (x, y, z,...). Problem (in three dimensions): Given E = J = J (1) x + J (2) y + J (3) z. Find J = (J (1), J (2), J (3) ). Application: Compute densities and fluxes of multi-dimensional systems of PDEs (in t, x, y, z). Generalize discrete homotopy operator in multi-dimensions (n, m,...). 29
30 Higher Euler Operators Side by Side Continuous Case (for component u) L (0) u = u D x + D 2 x D 3 x + u x u 2x u 3x 2D x + 3D 2 x 4D 3 x + u x u 2x u 3x u 4x 3D x + 6D 2 x 10D 3 x + u 2x u 3x u 4x u 5x 4D x + 10D 2 x 20D 3 x + u 3x u 4x u 5x u 6x L (1) u = L (2) u = L (3) u = Discrete Case (for component ) L (0) = L (1) = L (2) = L (3) = (I + D 1 + D 2 + D 3 + ) (D 1 + 2D 2 + 3D 3 + 4D 4 + ) (D 2 + 3D 3 + 6D D 5 + ) (D 3 + 4D D D 6 + ) 30
31 Homotopy Operators Side by Side Continuous Case (for components u and v) with and H u = 1 0 {f 1(u)[λu] + f 2 (u)[λu]} dλ λ f 1 (u) = p 1 i=0 f 2 (u) = p 2 i=0 D i x[ul (i+1) u ] D i x[vl (i+1) v ] Discrete Case (for components and v n ) with and H un = 1 0 {f 1,n( )[λ ] + f 2,n ( )[λ ]} dλ λ f 1,n ( ) = p 1 (D I) i [ L (i+1) i=0 ] f 2,n ( ) = p 2 (D I) i [v n L (i+1) i=0 v n ] 31
Symbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN
. Symbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN Joint work with ÜNAL GÖKTAŞ and Grant Erdmann Graduate Seminar Department
More informationApplication of Euler and Homotopy Operators to Integrability Testing
Application of Euler and Homotopy Operators to Integrability Testing Willy Hereman Department of Applied Mathematics and Statistics Colorado School of Mines Golden, Colorado, U.S.A. http://www.mines.edu/fs
More informationarxiv:nlin/ v1 [nlin.si] 16 Nov 2003
Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/0311029v1 [nlinsi] 16 Nov 2003 Symbolic Computation of Conserved Densities, Generalized Symmetries, and Recursion Operators
More informationSolving 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/
More informationThe first three (of infinitely many) conservation laws for (1) are (3) (4) D t (u) =D x (3u 2 + u 2x ); D t (u 2 )=D x (4u 3 u 2 x +2uu 2x ); D t (u 3
Invariants and Symmetries for Partial Differential Equations and Lattices Λ Ünal Göktaοs y Willy Hereman y Abstract Methods for the computation of invariants and symmetries of nonlinear evolution, wave,
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions
Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions Willy Hereman and Douglas Poole Department of Mathematical and Computer Sciences Colorado
More informationSymbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations
Symbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations Ünal Göktaş and Willy Hereman Department of Computer Engineering Turgut Özal University, Keçiören, Ankara, Turkey
More informationAlgorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of. Nonlinear Evolution and Differential-Difference.
Algorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of Nonlinear Evolution and Differential-Difference Equations by Ünal GÖKTAŞ A thesis submitted to the Faculty and the
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations
Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations Willy Hereman Department of Applied Mathematics and Statistics Colorado School of Mines Golden, Colorado whereman@mines.edu
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations
Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado whereman@mines.edu
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions
. Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden,
More informationSYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman
. SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado
More informationSymbolic Computation of Conservation Laws, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations
Symbolic Computation of Conservation Laws, Generalized Symmetries, and Recursion Operators for Nonlinear Differential-Difference Equations Ünal Götaş and Willy Hereman Abstract Algorithms for the symbolic
More informationSymbolic Computation of. Conserved Densities and Fluxes for. Nonlinear Systems of Differential-Difference. Equations
Symbolic Computation of Conserved Densities and Fluxes for Nonlinear Systems of Differential-Difference Equations by Holly Eklund Copyright c 003 Holly Eklund, All rights reserved. ii A thesis submitted
More informationSymbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions
1 Symbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions Douglas Poole and Willy Hereman Department of Mathematical and Computer Sciences, Colorado
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationarxiv:nlin/ v1 [nlin.si] 16 Nov 2003
Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/0311028v1 [nlin.si] 16 Nov 2003 Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for
More informationHamiltonian 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
More informationSymbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions
1 Symbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions Douglas Poole and Willy Hereman Department of Mathematical and Computer Sciences, Colorado
More informationHamiltonian 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
More information2. 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
More informationConstruction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem
Proceedings of 4th Workshop Group Analysis of Differential Equations & Integrability 2009, 1 23 Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem George W. BLUMAN,
More informationRESEARCH ARTICLE. A symbolic algorithm for computing recursion operators of nonlinear PDEs
International Journal of Computer Mathematics Vol. 00, No. 00, Month 200, 1 27 RESEARCH ARTICLE A symbolic algorithm for computing recursion operators of nonlinear PDEs D.E. Baldwin and W. Hereman Department
More informationSymbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs.
Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs. Douglas Baldwin, Willy Hereman,
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationNUMERICAL METHODS FOR SOLVING NONLINEAR EVOLUTION EQUATIONS
NUMERICAL METHODS FOR SOLVING NONLINEAR EVOLUTION EQUATIONS Thiab R. Taha Computer Science Department University of Georgia Athens, GA 30602, USA USA email:thiab@cs.uga.edu Italy September 21, 2007 1 Abstract
More informationA Very General Hamilton Jacobi Theorem
University of Salerno 9th AIMS ICDSDEA Orlando (Florida) USA, July 1 5, 2012 Introduction: Generalized Hamilton-Jacobi Theorem Hamilton-Jacobi (HJ) Theorem: Consider a Hamiltonian system with Hamiltonian
More informationA Note on Nonclassical Symmetries of a Class of Nonlinear Partial Differential Equations and Compatibility
Commun. Theor. Phys. (Beijing, China) 52 (2009) pp. 398 402 c Chinese Physical Society and IOP Publishing Ltd Vol. 52, No. 3, September 15, 2009 A Note on Nonclassical Symmetries of a Class of Nonlinear
More informationB.7 Lie Groups and Differential Equations
96 B.7. LIE GROUPS AND DIFFERENTIAL EQUATIONS B.7 Lie Groups and Differential Equations Peter J. Olver in Minneapolis, MN (U.S.A.) mailto:olver@ima.umn.edu The applications of Lie groups to solve differential
More informationIntroduction 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
More informationA Recursion Formula for the Construction of Local Conservation Laws of Differential Equations
A Recursion Formula for the Construction of Local Conservation Laws of Differential Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada December
More informationMACSYMA 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
More informationA Short historical review Our goal The hierarchy and Lax... The Hamiltonian... The Dubrovin-type... Algebro-geometric... Home Page.
Page 1 of 46 Department of Mathematics,Shanghai The Hamiltonian Structure and Algebro-geometric Solution of a 1 + 1-Dimensional Coupled Equations Xia Tiecheng and Pan Hongfei Page 2 of 46 Section One A
More informationAbout Integrable Non-Abelian Laurent ODEs
About Integrable Non-Abelian Laurent ODEs T. Wolf, Brock University September 12, 2013 Outline Non-commutative ODEs First Integrals and Lax Pairs Symmetries Pre-Hamiltonian Operators Recursion Operators
More informationFibonacci tan-sec method for construction solitary wave solution to differential-difference equations
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 7 (2011) No. 1, pp. 52-57 Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations
More informationLOWELL JOURNAL NEWS FROM HAWAII.
U # 2 3 Z7 B 2 898 X G V ; U» VU UU U z z 5U G 75 G } 2 B G G G G X G U 8 Y 22 U B B B G G G Y «j / U q 89 G 9 j YB 8 U z U - V - 2 z j _ ( UU 9 -! V G - U GB 8 j - G () V B V Y 6 G B j U U q - U z «B
More informationPotential Symmetries and Differential Forms. for Wave Dissipation Equation
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 42, 2061-2066 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2013.36163 Potential Symmetries and Differential Forms for Wave Dissipation
More informationThe 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
More informationPresentation for the Visiting Committee & Graduate Seminar SYMBOLIC COMPUTING IN RESEARCH. Willy Hereman
Presentation for the Visiting Committee & Graduate Seminar SYMBOLIC COMPUTING IN RESEARCH Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, CO 80401-1887 Monday,
More informationPainlevé analysis and some solutions of variable coefficient Benny equation
PRAMANA c Indian Academy of Sciences Vol. 85, No. 6 journal of December 015 physics pp. 1111 11 Painlevé analysis and some solutions of variable coefficient Benny equation RAJEEV KUMAR 1,, R K GUPTA and
More informationENERGY PRESERVING INTEGRATION OF BI-HAMILTONIAN PARTIAL DIFFERENTIAL EQUATIONS
TWMS J. App. Eng. Math. V., N.,, pp. 75-86 ENERGY PRESERVING INTEGRATION OF BI-HAMILTONIAN PARTIAL DIFFERENTIAL EQUATIONS B. KARASÖZEN, G. ŞİMŞEK Abstract. The energy preserving average vector field (AVF
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 informationSYMBOLIC COMPUTATION OF LAX PAIRS OF NONLINEAR PARTIAL DIFFERENCE EQUATIONS
SYMBOLIC COMPUTATION OF LAX PAIRS OF NONLINEAR PARTIAL DIFFERENCE EQUATIONS by Terry J. Bridgman c Copyright by Terry J. Bridgman, 2018 All Rights Reserved A thesis submitted to the Faculty and the Board
More informationConditional symmetries of the equations of mathematical physics
W.I. Fushchych, Scientific Works 2003, Vol. 5, 9 16. Conditional symmetries of the equations of mathematical physics W.I. FUSHCHYCH We briefly present the results of research in conditional symmetries
More informationHomotopy Perturbation Method for Solving Systems of Nonlinear Coupled Equations
Applied Mathematical Sciences, Vol 6, 2012, no 96, 4787-4800 Homotopy Perturbation Method for Solving Systems of Nonlinear Coupled Equations A A Hemeda Department of Mathematics, Faculty of Science Tanta
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 informationv(x, 0) = g(x) where g(x) = f(x) U(x). The solution is where b n = 2 g(x) sin(nπx) dx. (c) As t, we have v(x, t) 0 and u(x, t) U(x).
Problem set 4: Solutions Math 27B, Winter216 1. The following nonhomogeneous IBVP describes heat flow in a rod whose ends are held at temperatures u, u 1 : u t = u xx < x < 1, t > u(, t) = u, u(1, t) =
More informationDiagonalization 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
More informationÉlie Cartan s Theory of Moving Frames
Élie Cartan s Theory of Moving Frames Orn Arnaldsson Department of Mathematics University of Minnesota, Minneapolis Special Topics Seminar, Spring 2014 The Story of Symmetry Felix Klein and Sophus Lie
More informationProlongation structure for nonlinear integrable couplings of a KdV soliton hierarchy
Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy Yu Fa-Jun School of Mathematics and Systematic Sciences, Shenyang Normal University, Shenyang 110034, China Received
More informationMath RE - Calculus I Trigonometry Limits & Derivatives Page 1 of 8. x = 1 cos x. cos x 1 = lim
Math 0-0-RE - Calculus I Trigonometry Limits & Derivatives Page of 8 Trigonometric Limits It has been shown in class that: lim 0 sin lim 0 sin lim 0 cos cos 0 lim 0 cos lim 0 + cos + To evaluate trigonometric
More informationHamiltonian 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
More informationResearch Article Equivalent Lagrangians: Generalization, Transformation Maps, and Applications
Journal of Applied Mathematics Volume 01, Article ID 86048, 19 pages doi:10.1155/01/86048 Research Article Equivalent Lagrangians: Generalization, Transformation Maps, and Applications N. Wilson and A.
More informationCanonically conjugate variables for the periodic Camassa-Holm equation
arxiv:math-ph/21148v2 15 Mar 24 Canonically conjugate variables for the periodic Camassa-Holm equation Alexei V. Penskoi Abstract The Camassa-Holm shallow water equation is known to be Hamiltonian with
More informationEnergy Preserving Numerical Integration Methods
Energy Preserving Numerical Integration Methods Department of Mathematics and Statistics, La Trobe University Supported by Australian Research Council Professorial Fellowship Collaborators: Celledoni (Trondheim)
More informationGroup Actions and Cohomology in the Calculus of Variations
Group Actions and Cohomology in the Calculus of Variations JUHA POHJANPELTO Oregon State and Aalto Universities Focused Research Workshop on Exterior Differential Systems and Lie Theory Fields Institute,
More informationRecurrence in the KdV Equation? A. David Trubatch United States Military Academy/Montclair State University
Nonlinear Physics: Theory and Experiment IV Gallipoli, Lecce, Italy 27 June, 27 Recurrence in the KdV Equation? A. David Trubatch United States Military Academy/Montclair State University david.trubatch@usma.edu
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationHamiltonian 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 Lecture 2 Recall: the main goal is to compare
More informationArtificial boundary conditions for dispersive equations. Christophe Besse
Artificial boundary conditions for dispersive equations by Christophe Besse Institut Mathématique de Toulouse, Université Toulouse 3, CNRS Groupe de travail MathOcéan Bordeaux INSTITUT de MATHEMATIQUES
More informationGroup analysis, nonlinear self-adjointness, conservation laws, and soliton solutions for the mkdv systems
ISSN 139-5113 Nonlinear Analysis: Modelling Control, 017, Vol., No. 3, 334 346 https://doi.org/10.15388/na.017.3.4 Group analysis, nonlinear self-adjointness, conservation laws, soliton solutions for the
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 informationSymmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota
Symmetry Methods for Differential Equations and Conservation Laws Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Santiago, November, 2010 Symmetry Groups of Differential Equations
More informationON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS
Proceedings of the International Conference on Difference Equations, Special Functions and Orthogonal Polynomials, World Scientific (2007 ON THE SYMMETRIES OF INTEGRABLE PARTIAL DIFFERENCE EQUATIONS ANASTASIOS
More informationInvariant Variational Problems & Invariant Curve Flows
Invariant Variational Problems & Invariant Curve Flows Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Oxford, December, 2008 Basic Notation x = (x 1,..., x p ) independent variables
More informationMATH20411 PDEs and Vector Calculus B
MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411
More information2.4 Eigenvalue problems
2.4 Eigenvalue problems Associated with the boundary problem (2.1) (Poisson eq.), we call λ an eigenvalue if Lu = λu (2.36) for a nonzero function u C 2 0 ((0, 1)). Recall Lu = u. Then u is called an eigenfunction.
More informationChapter 2 Boundary and Initial Data
Chapter 2 Boundary and Initial Data Abstract This chapter introduces the notions of boundary and initial value problems. Some operator notation is developed in order to represent boundary and initial value
More informationCLASSIFICATION 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
More informationThe Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations
Symmetry in Nonlinear Mathematical Physics 1997, V.1, 185 192. The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations Norbert EULER, Ove LINDBLOM, Marianna EULER
More informationSoliton and Periodic Solutions to the Generalized Hirota-Satsuma Coupled System Using Trigonometric and Hyperbolic Function Methods.
ISSN 1749-889 (print), 1749-897 (online) International Journal of Nonlinear Science Vol.14(01) No.,pp.150-159 Soliton and Periodic Solutions to the Generalized Hirota-Satsuma Coupled System Using Trigonometric
More informationMA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation
MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation Lecture 12 Lecture 12 MA 201, PDE (2016) 1 / 24 Formal Solution of
More informationComplete integrability of a modified vector derivative nonlinear Schrödinger equation
Complete integrability of a modified vector derivative nonlinear Schrödinger equation arxiv:solv-int/9407002v1 14 Jul 1994 Ralph Willox 1, Willy Hereman 2 and Frank Verheest 3 1 Dienst Theoretische Natuurkunde,
More informationBreaking soliton equations and negative-order breaking soliton equations of typical and higher orders
Pramana J. Phys. (2016) 87: 68 DOI 10.1007/s12043-016-1273-z c Indian Academy of Sciences Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders ABDUL-MAJID
More informationarxiv: v2 [nlin.si] 3 Aug 2017
A generalised multicomponent system of Camassa-Holm-Novikov equations arxiv:1608.04604v2 [nlin.si] 3 Aug 2017 1. Introduction Diego Catalano Ferraioli 1 and Igor Leite Freire 2 Departamento de Matemática
More informationPartial Differential Equations, Winter 2015
Partial Differential Equations, Winter 215 Homework #2 Due: Thursday, February 12th, 215 1. (Chapter 2.1) Solve u xx + u xt 2u tt =, u(x, ) = φ(x), u t (x, ) = ψ(x). Hint: Factor the operator as we did
More informationDifferential Equation Types. Moritz Diehl
Differential Equation Types Moritz Diehl Overview Ordinary Differential Equations (ODE) Differential Algebraic Equations (DAE) Partial Differential Equations (PDE) Delay Differential Equations (DDE) Ordinary
More informationDifferential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension
Symposium in honor of George Bluman UBC May 12 16, 2014 Differential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension Chaolu Shanghai Maritime University
More informationIntroduction to the Calculus of Variations
236861 Numerical Geometry of Images Tutorial 1 Introduction to the Calculus of Variations Alex Bronstein c 2005 1 Calculus Calculus of variations 1. Function Functional f : R n R Example: f(x, y) =x 2
More informationExam in Numerical Methods (MA2501)
Norwegian University of Science and Technology Department of Mathematical Sciences Page 1 of 7 MA251 Numeriske Metoder Olivier Verdier (contact: 48 95 2 66) Exam in Numerical Methods (MA251) 211-5-25,
More informationBoundary value problems for integrable equations compatible with the symmetry algebra
Boundary value problems for integrable equations compatible with the symmetry algebra Burak Gürel, Metin Gürses, and Ismagil Habibullin Citation: J. Math. Phys. 36, 6809 (1995); doi: 10.1063/1.531189 View
More informationConstructive Algebra for Differential Invariants
Constructive Algebra for Differential Invariants Evelyne Hubert Rutgers, November 2008 Constructive Algebra for Differential Invariants Differential invariants arise in equivalence problems and are used
More informationDirect construction method for conservation laws of partial differential equations Part I: Examples of conservation law classifications
Euro. Jnl of Applied Mathematics (2002), vol. 13, pp. 545 566. c 2002 Cambridge University Press DOI: 10.1017/S095679250100465X Printed in the United Kingdom 545 Direct construction method for conservation
More informationAyuntamiento de Madrid
9 v vx-xvv \ ü - v q v ó - ) q ó v Ó ü " v" > - v x -- ü ) Ü v " ñ v é - - v j? j 7 Á v ü - - v - ü
More informationThe (G'/G) - Expansion Method for Finding Traveling Wave Solutions of Some Nonlinear Pdes in Mathematical Physics
Vol.3, Issue., Jan-Feb. 3 pp-369-376 ISSN: 49-6645 The ('/) - Expansion Method for Finding Traveling Wave Solutions of Some Nonlinear Pdes in Mathematical Physics J.F.Alzaidy Mathematics Department, Faculty
More informationNew methods of reduction for ordinary differential equations
IMA Journal of Applied Mathematics (2001) 66, 111 125 New methods of reduction for ordinary differential equations C. MURIEL AND J. L. ROMERO Departamento de Matemáticas, Universidad de Cádiz, PO Box 40,
More informationOn universality of critical behaviour in Hamiltonian PDEs
Riemann - Hilbert Problems, Integrability and Asymptotics Trieste, September 23, 2005 On universality of critical behaviour in Hamiltonian PDEs Boris DUBROVIN SISSA (Trieste) 1 Main subject: Hamiltonian
More informationGroup classification of nonlinear wave equations
JOURNAL OF MATHEMATICAL PHYSICS 46, 053301 2005 Group classification of nonlinear wave equations V. Lahno a State Pedagogical University, 36000 Poltava, Ukraine R. Zhdanov b Institute of Mathematics of
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More informationContinuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China
Continuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China the 3th GCOE International Symposium, Tohoku University, 17-19
More informationSymmetries and Group Invariant Reductions of Integrable Partial Difference Equations
Proceedings of 0th International Conference in MOdern GRoup ANalysis 2005, 222 230 Symmetries and Group Invariant Reductions of Integrable Partial Difference Equations A. TONGAS, D. TSOUBELIS and V. PAPAGEORGIOU
More informationSymmetry Reductions of (2+1) dimensional Equal Width. Wave Equation
Authors: Symmetry Reductions of (2+1) dimensional Equal Width 1. Dr. S. Padmasekaran Wave Equation Asst. Professor, Department of Mathematics Periyar University, Salem 2. M.G. RANI Periyar University,
More informationSoliton Solutions of the Time Fractional Generalized Hirota-satsuma Coupled KdV System
Appl. Math. Inf. Sci. 9, No., 17-153 (015) 17 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.1/amis/090 Soliton Solutions of the Time Fractional Generalized Hirota-satsuma
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationIntroduction to invariant and equivariant problems
Introduction to invariant and equivariant problems The curve completion problem Consider the curve completion problem, which is a subproblem of the much more complex inpainting problem. Suppose we are
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationLinearization of Mirror Systems
Journal of Nonlinear Mathematical Physics 00, Volume 9, Supplement 1, 34 4 Proceedings: Hong Kong Linearization of Mirror Systems Tat Leung YEE Department of Mathematics, The Hong Kong University of Science
More informationDispersion relations, stability and linearization
Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential
More informationRNDr. Petr Tomiczek CSc.
HABILITAČNÍ PRÁCE RNDr. Petr Tomiczek CSc. Plzeň 26 Nonlinear differential equation of second order RNDr. Petr Tomiczek CSc. Department of Mathematics, University of West Bohemia 5 Contents 1 Introduction
More informationDmitry Pelinovsky Department of Mathematics, McMaster University, Canada
Spectrum of the linearized NLS problem Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaboration: Marina Chugunova (McMaster, Canada) Scipio Cuccagna (Modena and Reggio Emilia,
More information