SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman
|
|
- Hollie Carroll
- 5 years ago
- Views:
Transcription
1 . 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 Kruskal Fest Symposium in Applied Mathematics: Nonlinear Waves, Dynamics, Asymtotic Analysis and Physical Applications Boulder, Colorado August 3-6, 1995
2 I. INTRODUCTION Symbolic Software Painlevé test for systems of ODEs and PDEs (Macsyma & Mathematica) Conservation laws of systems of evolution equations (Mathematica) Solitons via Hirota s method (Macsyma & Mathematica) Lie symmetries for systems of ODEs and PDEs (Macsyma) Purpose of the programs Study of integrability of nonlinear PDEs Exact solutions as bench mark for numerical algorithms Classification of nonlinear PDEs Lie symmetries solutions via reductions
3 Collaborators Ünal Göktaş, Chris Elmer, Wuning Zhuang (MS students) Ameina Nuseir (Ph.D student) Mark Coffey (CU-Boulder) Tony Miller & Tracy Otto (BS students)
4 II. SYMBOLIC SOFTWARE Program 1 Macsyma Painlevé Integrability Test for Systems of ODEs and PDEs Integrability of (a systems of) ODEs or PDEs requires that the only movable singularities in its solution are poles Definition: A single equation or system has the Painlevé Property if its solution in the complex plane has no worse singularities than movable poles Aim: Verify whether or not the system of equations satisfies the necessary criteria to have the Painlevé Property For simplicity, consider the case of a single PDE The solution f expressed as a Laurent series f = g α should only have movable poles k=0 u k g k
5 Steps of the Painlevé Test Step 1: 1. Substitute the leading order term into the given equation f u 0 g α 2. Determine the integer α < 0 by balancing the most singular terms in g 3. Calculate u 0 Step 2: 1. Substitute the generic terms f u 0 g α + u r g α+r into the equation, retaining its most singular terms 2. Require that u r is arbitrary 3. Determine the corresponding values of r > 0 called resonances
6 Step 3: 1. Substitute the truncated expansion f = g α R k=0 u k g k into the complete equation (R represents the largest resonance) 2. Determine u k unambiguously at the non-resonance levels 3. Check whether or not the compatibility condition is satisfied at resonance levels An equation or system has the Painlevé Property and is conjectured to be integrable if: 1. Step 1 thru 3 can be carried out consistently with α < 0 and with positive resonances 2. The compatibility conditions are identically satisfied for all resonances The above algorithm does not detect essential singularities
7 Painlevé Integrability Test Painlevé test for 3rd order equations by Hajee (Reduce, 1982) Painlevé program (parts) by Hlavatý (Reduce, 1986) ODE Painlevé by Winternitz & Rand (Macsyma, 1986) PDE Painlevé by Hereman & Van den Bulck (Macsyma, 1987) Painlevé test by Conte & Musette (AMP, 1988) Painlevé analysis by Renner (Reduce, 1992) Painlevé test for systems of ODEs and PDEs by Hereman, Elmer and Göktaş (Macsyma, , under development) Painlevé test for single ODEs and PDEs by Hereman, Miller and Otto (Mathematica, 1995, under development) Painlevé test for systems of ODEs and PDEs by Hereman and Miller (Mathematica, , planned)
8 Painlevé Test of Systems System of ODEs due to Akhiezer u xx + 2u 3 2au + 2bv = 0 v xx + 2v 3 2av + 2bu = 0 u(x) = g α 1 (x) k=0 v(x) = g α 2 (x) Leading orders: α 1 = α 2 = 1 Resonances: r = 1, r = 4 Coefficients: k=0 u k g k (x) v k g k (x) u 0 = v 0 = i u 1 = v 1 = 0 u 2 = v 2 = i(b a)/3 u 3 = v 3 = 0 At level 4: compatibility condition is satisfied Coefficients u 4 and v 4 are arbitrary and independent The system passes the Painlevé test
9 Carleman system (simplified version) u t + u x c uv = 0 v t v x + c uv = 0 u(x, t) = g α 1 (x, t) k=0 v(x, t) = g α 2 (x, t) k=0 Leading orders: α 1 = α 2 = 1 Resonances: r = 1, r = 1 Coefficients: u k (x, t) g k (x, t) v k (x, t) g k (x, t) u 0 = 1 c (g t g x ), v 0 = 1 c (g t + g x ) At level 1: compatibility condition is satisfied However u 1 and v 1 dependent on each other u 1 = g xx g tt + c v 1 (g t g x ) c (g t + g x ) The system passes the Painlevé test
10 Coupled KdV Equations (Hirota-Satsuma system) u t 3uu x + 6vv x 1 2 u xxx = 0 v t + v xxx + 3uv x = 0 Use Kruskal simplification, g(x, t) = x h(t), in u(x, t) = g α 1 (x, t) k=0 v(x, t) = g α 2 (x, t) k=0 u k (x, t) g k (x, t) v k (x, t) g k (x, t) Leading orders: First Branch α 1 = α 2 = 2 Resonances: r = 2, 1, 3, 4, 6, 8 Coefficients: u 0 = 4, v 0 = 2 u 1 = v 1 = 0 u 2 = h t /3, v 2 = 2h t /3 u 5 = h tt + 30v 3 h t, v 5 = h tt 12v 3 h t u 7 = v 4t + 12v 3 v 4, v 7 = 4h th tt + v 3 (8h 2 t 84v 4 ) + 21v 4t At level 3: compatibility condition is satisfied Here u 3 = v 3, one function is arbitary
11 At level 4: compatibility condition is satisfied Here u 4 = 2v 4, one function is arbitrary At level 6: compatibility condition is satisfied Here u 6 = v 2 3t + 24v 6 3v 3, 12 one function is arbitrary At level 8: compatibility condition is satisfied Here u 8 = h ttt + 6v 3 h tt + (12v 3t + 198v3)h 2 2 t 3024v 8 756v 4, 756 one function is arbitrary
12 Leading orders: Second Branch α 1 = 2, α 2 = 1 Resonances: r = 1, 0, 1, 4, 5, 6 Coefficients: u 0 = 2, v 0 free u 1 = 0, v 1 free u 2 = (2h t + 3v 2 0)/6, v 2 = (4v 0 h t + 3v 3 0)/12 u 3 = v 0 v 1, v 3 = v 0t + 3v 2 0v 1 12 At level 0: compatibility condition is satisfied Coefficient v 0 is arbitary At level 1: compatibility condition is satisfied Coefficient v 1 is arbitary At level 4: compatibility condition is satisfied Here u 4 = 16v 0h 2 t + 24v0h 3 t 24v 1t + 9v u 4 v 0, 288 one function is arbitrary
13 At level 5: compatibility condition is satisfied Here u 5 = 2h tt 12v 0 v 1 h t + 3v 0 v 0t 9v 0 3 v 1 18 one function is arbitrary At level 6: compatibility condition is satisfied Here u 6 = ( 64v 0h 3 t 144v 3 0h 2 t +(96v 1t 108v u 4 v 0 )h t +24v 0 2 v 1t 16v 0tt + 288v 3 0v v u 4 v u 6 v 0 ) one function is arbitrary The system passes the Painlevé test
14 Program 2 Mathematica Conserved Densities Purpose Compute polynomial-type conservation laws of single evolution equations and systems of evolution equations For simplicity, consider a single evolution equation u t = F(u, u x, u xx,..., u nx ) Conservation law is of the form ρ t + J x = 0 both ρ(u, u x, u 2x,..., u nx ) and J(u, u x, u 2x,..., u nx ) are polynomials in their arguments Consequently P = + provided J vanishes at infinity ρ dx = constant
15 Conservation laws describe the conservation of fundamental physical quantities such as mass, linear momentum, total energy (compare with constants of motion in mechanics) For nonlinear PDEs, the existence of a sufficiently large (in principal infinite) number of conservation laws assures complete integrability Tool to test numerical integrators for PDEs Example Consider the KdV equation, u t + uu x + u 3x = 0 Conserved densities: ρ 1 = u ρ 2 = u 2 ρ 3 = u 3 3u 2 x. ρ 6 = u 6 60u 3 u 2 x 30u 4 x + 108u 2 u 2 2x u3 2x uu2 3x u2 4x. Integrable equations have many conservation laws
16 Algorithm and Implementation Consider the scaling (weights) of the KdV u 2 x 2, Compute building blocks of ρ 3 (i) Start with building block u 3 t 3 x 3 Divide by u and differentiate twice (u 2 ) 2x Produces the list of terms [u 2 x, uu 2x ] [u 2 x] Second list: remove terms that are total derivative with respect to x or total derivative up to terms earlier in the list Divide by u 2 and differentiate twice (u) 4x Produces the list: [u 4x ] [ ] [ ] is the empty list
17 Gather the terms: ρ 3 = u 3 + c[1]u 2 x where the constant c 1 must be determined (ii) Compute ρ 3 t = 3u2 u t + 2c 1 u x u xt Replace u t by (uu x + u xxx ) and u xt by (uu x + u xxx ) x (iii) Integrate the result with respect to x Carry out all integrations by parts ρ 3 t = [ 3 4 u4 + (c 1 3)uu 2 x + 3u 2 u xx c 1 u 2 xx+ 2c 1 u x u xxx ] x (c 1 + 3)u 3 x The last non-integrable term must vanish Thus, c 1 = 3 Result: (iv) Expression [...] yields ρ 3 = u 3 3u 2 x J 3 = 3 4 u4 6uu 2 x + 3u 2 u xx + 3u 2 xx 6u x u xxx
18 Computer building blocks of ρ 6 (i) Start with u 6 Divide by u and differentiate twice (u 5 ) 2x produces the list of terms [u 3 u 2 x, u 4 u 2x ] [u 3 u 2 x] Next, divide u 6 by u 2, and compute (u 4 ) 4x Corresponding list: [u 4 x, uu 2 xu 2x, u 2 u 2 2x, u 2 u x u 3x, u 3 u 4x ] [u 4 x, u 2 u 2 2x] Proceed with ( u6 u 3 ) 6x = (u 3 ) 6x, ( u6 u 4 ) 8x = (u 2 ) 8x and ( u6 u 5 ) 10x = (u) 10x Obtain the lists: [u 3 2x, u x u 2x u 3x, uu 2 3x, u 2 xu 4x, uu 2x u 4x, uu x u 5x, u 2 u 6x ] [u 3 2x, uu 2 3x] [u 2 4x, u 3x u 5x, u 2x u 6x, u x u 7x, uu 8x ] [u 2 4x] and [u 10x ] [ ] Gather the terms: ρ 6 = u 6 + c 1 u 3 u 2 x + c 2 u 4 x + c 3 u 2 u 2 2x + c 4 u 3 2x + c 5 uu 2 3x + c 6 u 2 4x
19 where the constants c i must be determined (ii) Compute t ρ 6 Replace u t, u xt,..., u nx,t by (uu x + u xxx ),... (iii) Integrate the result with respect to x Carry out all integrations by parts Require that non-integrabe part vanishes Set to zero all the coefficients of the independent combinations involving powers of u and its derivatives with respect to x Solve the linear system for unknowns c 1, c 2,..., c 6 Result: ρ 6 = u 6 60u 3 u 2 x 30u 4 x + 108u 2 u 2 2x u3 2x uu2 3x u2 4x (iv) Flux J 6 can be computed by substituting the constants into the integrable part of ρ 6
20 Conserved Densities Conserved densities by Ito & Kako (Reduce, 1985, 1994) Conserved densities in DELiA by Bocharov (Pascal, 1990) Conserved densities by Gerdt (Reduce, 1993) Conserved densities by Roelofs, Sanders and Wang (Reduce 1994, Maple 1995) Conserved densities by Hereman and Göktaş (Mathematica, ) Conservation laws by Wolf (Reduce, 1995)
21 A Class of Fifth-order Evolution Equations Special cases: u t + αu 2 u x + βu x u 2x + γuu 3x + u 5x = 0 α = 30 β = 20 γ = 10 Lax α = 5 β = 5 γ = 5 Sawada Kotera or Caudry Dodd Gibbon α = 20 β = 25 γ = 10 Kaup Kuperschmidt α = 2 β = 6 γ = 3 Ito Under what conditions for the parameters α, β, and γ does this equation admit many conservation laws? If α = 3 10 γ2 and β = 2γ then there is a sequence (without gaps!) of conserved densities (Lax case) If α = 1 5 γ2 and β = γ then there is a sequence (with gaps!) of conserved densities (SK case) If α = 1 5 γ2 and β = 5 2γ then there is a sequence (with gaps!) of conserved densities (KK case) If α = 2β2 +7βγ 4γ 2 45 or β = 2γ then there is a conserved density of degree 4 Combining both conditions gives α = 2γ2 9 and β = 2γ (Ito case)
22 Conserved Densities of Systems of Evolution Eqs. Coupled KdV Equations (Hirota-Satsuma system) u t a(u xxx + 6uu x ) 2bvv x = 0 v t + v xxx + 3uv x = 0 u 2 x 2, v 2 x 2 ρ 1 = u ρ 2 = u bv2 ρ 3 = (1 + a)(u u2 x) + b(uv 2 v 2 x) and e.g. ρ 4 = u 4 2uu 2 x u2 xx b(u2 v uvv xx uv2 x 2 3 v2 xx) b2 v 4 provided a = 1 2 There are infinitely many more conservation laws
23 The Ito system u t u xxx 6uu x 2vv x = 0 v t 2u x v 2uv x = 0 u 2 x 2, v 2 x 2 ρ 1 = c 1 u + c 2 v ρ 2 = u 2 + v 2 ρ 3 = 2u 3 + 2uv 2 u 2 x ρ 4 = 5u 4 + 6u 2 v 2 + v 4 10uu 2 x + 2v 2 u 2x + u 2 2x and there are infinitely many more conservation laws
24 The Drinfel d-sokolov system u t + 3vv x = 0 v t + vu x + 2uv x + 2v 3x = 0 u 2 x 2, 2 v free, choose v x 2 ρ 1 = u ρ 2 = v 2 ρ 3 = 2 9 u3 + uv u2 x 3 2 v2 x ρ 4 = 1 2 u4 1 6 u2 v v4 1 6 uu2 x + uv 2 x 5 12 v2 u 2x u2 2x 3 4 v2 2x and there are presumably infinitely many more conservation laws
25 The dispersiveless long-wave system u t + vu x + uv x = 0 v t + u x + vv x = 0 u free, v free, but u 2v choose u x and 2v x ρ 1 = v ρ 2 = u ρ 3 = uv ρ 4 = u 2 + uv 2 ρ 5 = 3u 2 v + uv 3 ρ 6 = 1 3 u3 + u 2 v uv4 ρ 7 = u 3 v + u 2 v uv5 ρ 8 = 1 3 u4 + 2u 3 v 2 + u 2 v uv6 always the same set irrespective the choice of weights
26 Further Examples A generalized Schamel equation n 2 u t + (n + 1)(n + 2)unu 2 x + u xxx = 0 n positive integer ρ 1 = u, ρ 2 = u 2 no further conservation laws The Boussinesq equation ρ 3 = 1 2 u2 x n2 2 u2+ 2 n u tt bu xx + 3uu xx + 3u 2 x + au 4x = 0 a and b are real constants Write as system of evolution equations u t + v x = 0 v t + bu x 3uu x au 3x = 0 u 2 x 2, b 2 x 2, v free, take v 3 x 3
27 ρ 1 = u ρ 2 = u ρ 3 = bu ρ 4 = bc 1 u + c 2 uv ρ 5 = b 2 c 1 u + c 2 (bu 2 u 3 + v 2 + au 2 x) and there are infinitely many more conservation laws
28 A modified vector derivative NLS equation B t + x (B2 B ) + αb 0 B 0 B x + e x 2 B x 2 = 0 Replace vector equation by u t + ( u(u 2 + v 2 ) + βu v x )x = 0 v t + ( v(u 2 + v 2 ) + u x )x = 0 u and v denote the components of B parallel and perpendicular to B 0 and β = αb 2 0 2u x, 2v x, The first 6 conserved densities are: ρ 1 = c 1 u + c 2 v ρ 2 = u 2 + v 2 β x ρ 3 = 1 2 (u2 + v 2 ) 2 uv x + u x v + βu 2 ρ 4 = 1 4 (u2 + v 2 ) (u2 x + v 2 x) u 3 v x + v 3 u x + β 4 (u4 v 4 )
29 ρ 5 = 1 4 (u2 + v 2 ) (u xv xx u xx v x ) (uu x + vv x ) (u2 + v 2 )(u 2 x + v 2 x) (u 2 + v 2 ) 2 (uv x u x v) + β 5 (2u2 x 4u 3 v x + 2u 6 + 3u 4 v 2 v 6 ) + β2 5 u4 ρ 6 = 7 16 (u2 + v 2 ) (u2 xx + v 2 xx) 5 2 (u2 + v 2 )(u x v xx u xx v x ) + 5(u 2 + v 2 )(uu x + vv x ) (u2 + v 2 ) 2 (u 2 x + v 2 x) (u2 + v 2 ) 3 (uv x u x v) + β 8 (5u8 + 10u 6 v 2 10u 2 v 6 5v u 2 u 2 x 12u 5 v x + 60uv 4 v x 20v 2 v 2 x) + β2 4 (u6 + v 6 )
30 Table 1 Conserved Densities for Sawada-Kotera and Lax equations Density Sawada-Kotera equation Lax equation ρ 1 u u ρ u2 ρ u3 u 2 x 1 3 u3 1 6 u2 x ρ u4 9 4 uu2 x u2 2x 1 4 u4 1 2 uu2 x u2 2x ρ u5 u 2 u 2 x uu2 2x 1 70 u2 3x ρ u u3 u 2 x 17 8 u4 x + 6u 2 u 2 2x 1 6 u6 5 3 u3 u 2 x 5 36 u4 x u2 u 2 2x +2u 3 2x 21 8 uu2 3x u2 4x u3 2x 1 14 uu2 3x u2 4x ρ u7 9u 4 u 2 x 54 5 uu4 x u3 u 2 2x 1 7 u7 5 2 u4 u 2 x 5 6 uu4 x + u 3 u 2 2x u2 xu 2 2x uu3 2x u2 u 2 3x u2 xu 2 2x uu3 2x 3 14 u2 u 2 3x u 2xu 2 3x uu2 4x 9 35 u2 5x 5 42 u 2xu 2 3x uu2 4x u2 5x ρ u8 7 2 u5 u 2 x u2 u 4 x u4 u 2 2x uu2 xu 2 2x u2 u 3 2x u4 2x u3 u 2 3x 1 4 u2 xu 2 3x 5 6 uu 2xu 2 3x u2 u 2 4x u 2xu 2 4x uu2 5x u2 6x
31 Table 2 Conserved Densities for Kaup-Kuperschmidt and Ito equations Density Kaup-Kuperschmidt equation Ito equation ρ 1 u u ρ u 2 2 ρ 3 u u2 x ---- ρ 4 u uu2 x u2 2x u uu2 x u2 2x ρ ρ 6 u u3 u 2 x u4 x u2 u 2 2x u3 2x uu2 3x u2 4x ρ 7 u u4 u 2 x uu4 x u3 u 2 2x u2 xu 2 2x uu3 2x u2 u 2 3x u 2xu 2 3x uu2 4x u2 5x ρ
32 Example 3 Macsyma/Mathematica Solitons Hirota s Method Hirota s Direct Method allows to construct soliton solutions of nonlinear evolution equations wave equations coupled systems Test conditions for existence of soliton solutions Examples: Korteweg-de Vries equation (KdV) u t + 6uu x + u 3x = 0 Kadomtsev-Petviashvili equation (KP) (u t + 6uu x + u 3x ) x + 3u 2y = 0 Sawada-Kotera equation (SK) u t + 45u 2 u x + 15u x u 2x + 15uu 3x + u 5x = 0
33 Hirota s Method Korteweg-de Vries equation u t + 6uu x + u 3x = 0 Substitute u(x, t) = 2 2 ln f(x, t) x 2 Integrate with respect to x Bilinear form ff xt f x f t + ff 4x 4f x f 3x + 3f 2 2x = 0 B(f f) def = ( D x D t + D 4 x) (f f) = 0 Introduce the bilinear operator D m x D n t (f g) = ( x x ) m ( t t ) n f(x, t) g(x, t ) x =x,t =t Use the expansion f = 1 + n=1 ɛn f n Substitute f into the bilinear equation
34 Collect powers in ɛ (book keeping parameter) Start with O(ɛ 0 ) : B(1 1) = 0 O(ɛ 1 ) : B(1 f 1 + f 1 1) = 0 O(ɛ 2 ) : B(1 f 2 + f 1 f 1 + f 2 1) = 0 O(ɛ 3 ) : B(1 f 3 + f 1 f 2 + f 2 f 1 + f 3 1) = 0 O(ɛ 4 ) : B(1 f 4 + f 1 f 3 + f 2 f 2 + f 3 f 1 + f 4 1) = 0 O(ɛ n ) : B( n j=0 f j f n j ) = 0 with f 0 = 1 f 1 = N exp(θ i) = N exp (k i x ω i t + δ i ) i=1 i=1 k i, ω i and δ i are constants Dispersion law ω i = k 3 i (i = 1, 2,..., N) If the original PDE admits a N-soliton solution then the expansion will truncate at level n = N
35 Consider the case N=3 Terms generated by B(f 1, f 1 ) determine f 2 = a 12 exp(θ 1 + θ 2 ) + a 13 exp(θ 1 + θ 3 ) + a 23 exp(θ 2 + θ 3 ) = a 12 exp [(k 1 + k 2 ) x (ω 1 + ω 2 ) t + (δ 1 + δ 2 )] + a 13 exp [(k 1 + k 3 ) x (ω 1 + ω 3 ) t + (δ 1 + δ 3 )] + a 23 exp [(k 2 + k 3 ) x (ω 2 + ω 3 ) t + (δ 2 + δ 3 )] Calculate the constants a 12, a 13 and a 23 a ij = (k i k j ) 2 (k i + k j ) 2 i, j = 1, 2, 3 Terms from B(f 1 f 2 + f 2 f 1 ) determine with f 3 = b 123 exp(θ 1 + θ 2 + θ 3 ) = b 123 exp [(k 1 +k 2 +k 3 )x (ω 1 +ω 2 +ω 3 )t+(δ 1 +δ 2 +δ 3 )] b 123 = a 12 a 13 a 23 = (k 1 k 2 ) 2 (k 1 k 3 ) 2 (k 2 k 3 ) 2 (k 1 + k 2 ) 2 (k 1 + k 3 ) 2 (k 2 + k 3 ) 2 Subsequently, f i = 0 for i > 3 Set ɛ = 1 f = 1 + exp θ 1 + exp θ 2 + exp θ 3 + a 12 exp(θ 1 + θ 2 ) + a 13 exp(θ 1 + θ 3 ) + a 23 exp(θ 2 + θ 3 ) + b 123 exp(θ 1 + θ 2 + θ 3 )
36 Return to the original u(x, t) Single soliton solution u(x, t) = 2 2 ln f(x, t) x 2 f = 1 + e θ, k, ω and δ are constants and ω = k 3 Substituting f into θ = kx ωt + δ Take k = 2K u(x, t) = 2 2 ln f(x, t) x 2 = 2( f xxf f 2 x f 2 ) u = 2K 2 sech 2 K(x 4K 2 t + δ)
37 Two-soliton solution f = 1 + e θ 1 + e θ 2 + a 12 e θ 1+θ 2 θ i = k i x ω i t + δ i with ω i = k 3 i, (i = 1, 2) and a 12 = (k 1 k 2 ) 2 (k 1 +k 2 ) 2 Select e δ i = c2 i k i e k ix ω i t+ i f = 1 4 fe 1 2 ( θ 1 + θ 2 ) θ i = k i x ω i t + i c 2 i = k 2 + k 1 k 2 k 1 k i Return to u(x, t) u(x, t) = ũ(x, t) = 2 2 ln f(x, t) x 2 = k 2 2 k k2cosech 2 2 θ2 2 + k1sech 2 2 θ1 2 (k 2 coth θ 2 2 k 1 tanh θ 1 2 ) 2
38 Program 4 Macsyma Lie-point Symmetries System of m differential equations of order k i (x, u (k) ) = 0, i = 1, 2,..., m with p independent and q dependent variables x = (x 1, x 2,..., x p ) IR p u = (u 1, u 2,..., u q ) IR q The group transformations have the form x = Λ group (x, u), ũ = Ω group (x, u) where the functions Λ group and Ω group are to be determined Look for the Lie algebra L realized by the vector field α = p i=1 ηi (x, u) x i + q l=1 ϕ l (x, u) u l
39 Procedure for finding the coefficients Construct the k th prolongation pr (k) α of the vector field α Apply it to the system of equations Request that the resulting expression vanishes on the solution set of the given system pr (k) α i j =0 i, j = 1,..., m This results in a system of linear homogeneous PDEs for η i and ϕ l, with independent variables x and u ( determining equations) Procedure thus consists of two major steps: deriving the determining equations solving the determining equations
40 Procedure for Computing the Determining Equations Use multi-index notation J = (j 1, j 2,..., j p ) IN p, to denote partial derivatives of u l u l J where J = j 1 + j j p J u l x 1 j 1 x2 j 2... xp j p, u (k) denotes a vector whose components are all the partial derivatives of order 0 up to k of all the u l Steps: (1) Construct the k th prolongation of the vector field pr (k) α = α + q l=1 u l J J ψj l (x, u (k) ), 1 J k The coefficients ψ J l of the first prolongation are: ψ J i l = D i ϕ l (x, u) p j=1 ul J j D i η j (x, u), where J i is a p tuple with 1 on the i th position and zeros elsewhere D i is the total derivative operator D i = + q x i l=1 J ul J+J i u l J, 0 J k
41 Higher order prolongations are defined recursively: ψ J+J i l = D i ψ J l p j=1 ul J+J j D i η j (x, u), J 1 (2) Apply the prolonged operator pr (k) α to each equation i (x, u (k) ) = 0 Require that pr (k) α vanishes on the solution set of the system pr (k) α i j =0 = 0 i, j = 1,..., m (3) Choose m components of the vector u (k), say v 1,..., v m, such that: (a) Each v i is equal to a derivative of a u l (l = 1,..., q) with respect to at least one variable x i (i = 1,..., p). (b) None of the v i is the derivative of another one in the set. (c) The system can be solved algebraically for the v i in terms of the remaining components of u (k), which we denoted by w: (d) The derivatives of v i, v i = S i (x, w), i = 1,..., m. v i J = D J S i (x, w),
42 where D J D j 1 1 D j D j p p, can all be expressed in terms of the components of w and their derivatives, without ever reintroducing the v i or their derivatives. For instance, for a system of evolution equations u i t(x 1,..., x p 1, t) = F i (x 1,..., x p 1, t, u (k) ), i = 1,..., m, where u (k) involves derivatives with respect to the variables x i but not t, choose v i = u i t. (4) Eliminate all v i and their derivatives from the expression prolonged vector field, so that all the remaining variables are independent (5) Obtain the determining equations for η i (x, u) and ϕ l (x, u) by equating to zero the coefficients of the remaining independent derivatives u l J.
43 Example: The Korteweg-de Vries Equation one equation (parameter a) u t + auu x + u xxx = 0 two independent variables t and x one dependent variable u vector field α = η x x + ηt t + ϕu u Format for SYMMGRP.MAX variables x[1] = x, x[2] = t, u[1] = u equation e1 : u[1, [0, 1]] + a u[1] u[1, [1, 0]] + u[1, [3, 0]] variable to be eliminated v1 : u[1, [0, 1]] coefficients of vectorfield in SYMMGRP.MAX: eta[1] = η x, eta[2] = η t and phi[1] = ϕ u
44 There are only eight determining equations phi[1] x[2] eta[2] u[1] eta[2] x[1] eta[1] u[1] = 0 = 0 = 0 2 phi[1] u[1] 2 = 0 2 phi[1] u[1] x[1] 2 eta[1] = 0 x[1] phi[1] x[1] 3 phi[1] + u[1] x[1] = 0
45 3 3 phi[1] u[1] x[1] 2 eta[1] x[2] 3 eta[1] x[1] u[1] eta[1] x[1] + phi[1] = 0 u[1] eta[2] x[2] phi[1] u[1] x[1] 2 The solution in the original variables eta[1] x[2] η x = k 1 + k 3 t k 4 x η t = k 2 3k 4 t ϕ u = k k 4 u The four infinitesimal generators are 3 eta[1] x[1] 3 u[1] eta[1] x[1] + phi[1] = 0 G 1 = x G 2 = t G 3 = t x + u G 4 = x x + 3 t t 2 u u
46 Equation is invariant under: translations G 1 and G 2 Galilean boost G 3 scaling G 4 Computation of the flows corresponding to G 1 thru G 4 shows that for any solution u = f(x, t) of the KdV equation the transformed solutions ũ = f(x ɛ, t) ũ = f(x, t ɛ) ũ = f(x ɛ, t) + ɛ will solve the KdV equation ũ = e 2ɛ f(e ɛ x, e 3ɛ t) Note that ɛ is the parameter of the transformation group
47 III. PLANS FOR THE FUTURE Extension of Symbolic Software Packages (Macsyma/Mathematica) Lie symmetries of differential-difference equations Solver for systems of linear, homogeneous PDEs (Hereman) Painlevé test for systems of PDEs (Elmer, Göktaş & Coffey) Solitons via Hirota s method for bilinear equations (Zhuang) Simplification of Hirota s method (Hereman & Nuseir) Conservation laws of PDEs with variable coefficients (Göktaş) Lax pairs, special solutions,... New Software Wavelets (prototype/educational tool) Other methods for Differential Equations
Presentation 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 informationSymbolic 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 informationNIST Electromagnetic Field Division Seminar SYMBOLIC COMPUTING IN RESEARCH. Willy Hereman
. NIST Electromagnetic Field Division Seminar SYMBOLIC COMPUTING IN RESEARCH Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, CO 80401-1887 Tuesday, March 10,
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 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 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 informationSolving 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
More informationContinuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman
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 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/0311028v1 [nlin.si] 16 Nov 2003 Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for
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 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 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 informationMultiple-Soliton Solutions for Extended Shallow Water Wave Equations
Studies in Mathematical Sciences Vol. 1, No. 1, 2010, pp. 21-29 www.cscanada.org ISSN 1923-8444 [Print] ISSN 1923-8452 [Online] www.cscanada.net Multiple-Soliton Solutions for Extended Shallow Water Wave
More informationarxiv: v1 [nlin.si] 23 Aug 2007
arxiv:0708.3247v1 [nlin.si] 23 Aug 2007 A new integrable generalization of the Korteweg de Vries equation Ayşe Karasu-Kalkanlı 1), Atalay Karasu 2), Anton Sakovich 3), Sergei Sakovich 4), Refik Turhan
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 informationExact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially Integrable Equations
Thai Journal of Mathematics Volume 5(2007) Number 2 : 273 279 www.math.science.cmu.ac.th/thaijournal Exact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 60 (00) 3088 3097 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Symmetry
More informationCoupled KdV Equations of Hirota-Satsuma Type
Journal of Nonlinear Mathematical Physics 1999, V.6, N 3, 255 262. Letter Coupled KdV Equations of Hirota-Satsuma Type S.Yu. SAKOVICH Institute of Physics, National Academy of Sciences, P.O. 72, Minsk,
More informationIntroduction to the Hirota bilinear method
Introduction to the Hirota bilinear method arxiv:solv-int/9708006v1 14 Aug 1997 J. Hietarinta Department of Physics, University of Turku FIN-20014 Turku, Finland e-mail: hietarin@utu.fi Abstract We give
More informationThe Construction of Alternative Modified KdV Equation in (2 + 1) Dimensions
Proceedings of Institute of Mathematics of NAS of Ukraine 00, Vol. 3, Part 1, 377 383 The Construction of Alternative Modified KdV Equation in + 1) Dimensions Kouichi TODA Department of Physics, Keio University,
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 informationResearch Article Application of the Cole-Hopf Transformation for Finding Exact Solutions to Several Forms of the Seventh-Order KdV Equation
Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 21, Article ID 194329, 14 pages doi:1.1155/21/194329 Research Article Application of the Cole-Hopf Transformation for Finding
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 AND SOFTWARE FOR THE PAINLEVÉ TEST AND RECURSION OPERATORS FOR NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS
SYMBOLIC ALGORITHMS AND SOFTWARE FOR THE PAINLEVÉ TEST AND RECURSION OPERATORS FOR NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS by Douglas E. Baldwin Copyright 2004 by Douglas E. Baldwin All Rights Reserved
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 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 informationPainlevé Test for the Certain (2+1)-Dimensional Nonlinear Evolution Equations. Abstract
Painlevé Test for the Certain (2+1)-Dimensional Nonlinear Evolution Equations T. Alagesan and Y. Chung Department of Information and Communications, Kwangju Institute of Science and Technology, 1 Oryong-dong,
More information-Expansion Method For Generalized Fifth Order KdV Equation with Time-Dependent Coefficients
Math. Sci. Lett. 3 No. 3 55-6 04 55 Mathematical Sciences Letters An International Journal http://dx.doi.org/0.785/msl/03039 eneralized -Expansion Method For eneralized Fifth Order KdV Equation with Time-Dependent
More informationSecond Order Lax Pairs of Nonlinear Partial Differential Equations with Schwarzian Forms
Second Order Lax Pairs of Nonlinear Partial Differential Equations with Schwarzian Forms Sen-yue Lou a b c, Xiao-yan Tang a b, Qing-Ping Liu b d, and T. Fukuyama e f a Department of Physics, Shanghai Jiao
More informationEXACT SOLITARY WAVE AND PERIODIC WAVE SOLUTIONS OF THE KAUP-KUPERSCHMIDT EQUATION
Journal of Applied Analysis and Computation Volume 5, Number 3, August 015, 485 495 Website:http://jaac-online.com/ doi:10.11948/015039 EXACT SOLITARY WAVE AND PERIODIC WAVE SOLUTIONS OF THE KAUP-KUPERSCHMIDT
More informationA supersymmetric Sawada-Kotera equation
A supersymmetric Sawada-Kotera equation arxiv:0802.4011v2 [nlin.si] 7 Dec 2008 Kai Tian and Q. P. Liu Department of Mathematics, China University of Mining and Technology, Beijing 100083, P.R. China Abstract
More informationDoes the Supersymmetric Integrability Imply the Integrability of Bosonic Sector?
Does the Supersymmetric Integrability Imply the Integrability of Bosonic Sector? Z I E M O W I T P O P O W I C Z Instytut Fizyki Teoretycznej Uniwersytet Wrocławski POLAND 15.07.2009-21.07.2009 Shihezi
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 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 informationBoussineq-Type Equations and Switching Solitons
Proceedings of Institute of Mathematics of NAS of Ukraine, Vol. 3, Part 1, 3 351 Boussineq-Type Equations and Switching Solitons Allen PARKER and John Michael DYE Department of Engineering Mathematics,
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 informationResearch Article Soliton Solutions for the Wick-Type Stochastic KP Equation
Abstract and Applied Analysis Volume 212, Article ID 327682, 9 pages doi:1.1155/212/327682 Research Article Soliton Solutions for the Wick-Type Stochastic KP Equation Y. F. Guo, 1, 2 L. M. Ling, 2 and
More informationA 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
More informationNew Exact Solutions for the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli Equation
Applied Mathematical Sciences, Vol. 6, 2012, no. 12, 579-587 New Exact Solutions for the (2+1)-dimensional Boiti-Leon-Manna-Pempinelli Equation Ying Li and Desheng Li School of Mathematics and System Science
More informationExact Solutions for a Fifth-Order Two-Mode KdV Equation with Variable Coefficients
Contemporary Engineering Sciences, Vol. 11, 2018, no. 16, 779-784 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.8262 Exact Solutions for a Fifth-Order Two-Mode KdV Equation with Variable
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 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 informationThe (2+1) and (3+1)-Dimensional CBS Equations: Multiple Soliton Solutions and Multiple Singular Soliton Solutions
The (+1) and (3+1)-Dimensional CBS Equations: Multiple Soliton Solutions and Multiple Singular Soliton Solutions Abdul-Majid Wazwaz Department of Mathematics, Saint Xavier University, Chicago, IL 60655,
More informationBäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations
arxiv:nlin/0102001v1 [nlin.si] 1 Feb 2001 Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations Ayşe Karasu (Kalkanli) and S Yu Sakovich Department
More informationGeneralized 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
More informationNew Approach of ( Ǵ/G ) Expansion Method. Applications to KdV Equation
Journal of Mathematics Research; Vol. 6, No. ; ISSN 96-9795 E-ISSN 96-989 Published by Canadian Center of Science and Education New Approach of Ǵ/G Expansion Method. Applications to KdV Equation Mohammad
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
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 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 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 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 informationNumerical solutions of the small dispersion limit of KdV, Whitham and Painlevé equations
Numerical solutions of the small dispersion limit of KdV, Whitham and Painlevé equations Tamara Grava (SISSA) joint work with Christian Klein (MPI Leipzig) Integrable Systems in Applied Mathematics Colmenarejo,
More informationResearch Article Solving Nonlinear Partial Differential Equations by the sn-ns Method
Abstract and Applied Analysis Volume 2012, Article ID 340824, 25 pages doi:10.1155/2012/340824 Research Article Solving Nonlinear Partial Differential Equations by the sn-ns Method Alvaro H. Salas FIZMAKO
More informationNew Formal Solutions of Davey Stewartson Equation via Combined tanh Function Method with Symmetry Method
Commun. Theor. Phys. Beijing China 7 007 pp. 587 593 c International Academic Publishers Vol. 7 No. April 5 007 New Formal Solutions of Davey Stewartson Equation via Combined tanh Function Method with
More informationHirota Direct Method Graham W. Griffiths 1 City University, UK. oo0oo
Hirota Direct Method Graham W. Griffiths 1 City University, UK. oo0oo Contents 1 Introduction 1 Applications of the Hirota method 3.1 Single-soliton solution to the KdV equation..........................
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 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 informationEXACT SOLUTION TO TIME FRACTIONAL FIFTH-ORDER KORTEWEG-DE VRIES EQUATION BY USING (G /G)-EXPANSION METHOD. A. Neamaty, B. Agheli, R.
Acta Universitatis Apulensis ISSN: 1582-5329 http://wwwuabro/auajournal/ No 44/2015 pp 21-37 doi: 1017114/jaua20154403 EXACT SOLUTION TO TIME FRACTIONAL FIFTH-ORDER KORTEWEG-DE VRIES EQUATION BY USING
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 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 informationBäcklund transformations for fourth-order Painlevé-type equations
Bäcklund transformations for fourth-order Painlevé-type equations A. H. Sakka 1 and S. R. Elshamy Department of Mathematics, Islamic University of Gaza P.O.Box 108, Rimae, Gaza, Palestine 1 e-mail: asakka@mail.iugaza.edu
More informationThe Solitary Wave Solutions of Zoomeron Equation
Applied Mathematical Sciences, Vol. 5, 011, no. 59, 943-949 The Solitary Wave Solutions of Zoomeron Equation Reza Abazari Deparment of Mathematics, Ardabil Branch Islamic Azad University, Ardabil, Iran
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
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 informationNew 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
More informationSymmetry classification of KdV-type nonlinear evolution equations
arxiv:nlin/0201063v2 [nlin.si] 16 Sep 2003 Symmetry classification of KdV-type nonlinear evolution equations F. Güngör Department of Mathematics, Faculty of Science and Letters, Istanbul Technical University,
More informationBiao Li a,b, Yong Chen a,b, and Hongqing Zhang a,b. 1. Introduction
Auto-Bäcklund Transformations and Exact Solutions for the Generalized Two-Dimensional Korteweg-de Vries-Burgers-type Equations and Burgers-type Equations Biao Li a,b, Yong Chen a,b, and Hongqing Zhang
More informationSOLITON SOLUTIONS OF SHALLOW WATER WAVE EQUATIONS BY MEANS OF G /G EXPANSION METHOD
Journal of Applied Analysis and Computation Website:http://jaac-online.com/ Volume 4, Number 3, August 014 pp. 1 9 SOLITON SOLUTIONS OF SHALLOW WATER WAVE EQUATIONS BY MEANS OF G /G EXPANSION METHOD Marwan
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 informationSolution of the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(2007) No.3,pp.227-234 Solution of the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition
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 informationThe General Form of Linearized Exact Solution for the KdV Equation by the Simplest Equation Method
Applied and Computational Mathematics 015; 4(5): 335-341 Published online August 16 015 (http://www.sciencepublishinggroup.com/j/acm) doi: 10.11648/j.acm.0150405.11 ISSN: 38-5605 (Print); ISSN: 38-5613
More informationResearch Article Exact Solutions of φ 4 Equation Using Lie Symmetry Approach along with the Simplest Equation and Exp-Function Methods
Abstract and Applied Analysis Volume 2012, Article ID 350287, 7 pages doi:10.1155/2012/350287 Research Article Exact Solutions of φ 4 Equation Using Lie Symmetry Approach along with the Simplest Equation
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 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 informationOn critical behaviour in systems of Hamiltonian partial differential equations
On critical behaviour in systems of Hamiltonian partial differential equations arxiv:1311.7166v1 [math-ph] 7 Nov 13 B. Dubrovin, T. Grava, C. Klein, A. Moro Abstract We study the critical behaviour of
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 informationBuilding Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients
Journal of Physics: Conference Series OPEN ACCESS Building Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients To cite this article: M Russo and S R Choudhury 2014
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 informationDispersionless integrable systems in 3D and Einstein-Weyl geometry. Eugene Ferapontov
Dispersionless integrable systems in 3D and Einstein-Weyl geometry Eugene Ferapontov Department of Mathematical Sciences, Loughborough University, UK E.V.Ferapontov@lboro.ac.uk Based on joint work with
More informationPeriodic and Soliton Solutions for a Generalized Two-Mode KdV-Burger s Type Equation
Contemporary Engineering Sciences Vol. 11 2018 no. 16 785-791 HIKARI Ltd www.m-hikari.com https://doi.org/10.12988/ces.2018.8267 Periodic and Soliton Solutions for a Generalized Two-Mode KdV-Burger s Type
More information(Received 05 August 2013, accepted 15 July 2014)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.18(2014) No.1,pp.71-77 Spectral Collocation Method for the Numerical Solution of the Gardner and Huxley Equations
More informationMATH 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
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 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 informationexp Φ ξ -Expansion Method
Journal of Applied Mathematics and Physics, 6,, 6-7 Published Online February 6 in SciRes. http://www.scirp.org/journal/jamp http://dx.doi.org/.6/jamp.6. Analytical and Traveling Wave Solutions to the
More informationIntegrability approaches to differential equations
XXIV Fall Workshop on Geometry and Physics Zaragoza, CUD, September 2015 What is integrability? On a fairly imprecise first approximation, we say integrability is the exact solvability or regular behavior
More informationA NEW VARIABLE-COEFFICIENT BERNOULLI EQUATION-BASED SUB-EQUATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS
U.P.B. Sci. Bull., Series A, Vol. 76, Iss., 014 ISSN 1-707 A NEW VARIABLE-COEFFICIENT BERNOULLI EQUATION-BASED SUB-EQUATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS Bin Zheng 1 In this paper,
More informationThe extended homogeneous balance method and exact 1-soliton solutions of the Maccari system
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol., No., 014, pp. 83-90 The extended homogeneous balance method and exact 1-soliton solutions of the Maccari system Mohammad
More informationarxiv:nlin/ v2 [nlin.si] 14 Sep 2001
Second order Lax pairs of nonlinear partial differential equations with Schwarz variants arxiv:nlin/0108045v2 [nlin.si] 14 Sep 2001 Sen-yue Lou 1,2,3, Xiao-yan Tang 2,3, Qing-Ping Liu 1,4,3 and T. Fukuyama
More informationRelation between Periodic Soliton Resonance and Instability
Proceedings of Institute of Mathematics of NAS of Ukraine 004 Vol. 50 Part 1 486 49 Relation between Periodic Soliton Resonance and Instability Masayoshi TAJIRI Graduate School of Engineering Osaka Prefecture
More informationAre 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
More informationCS205B / CME306 Homework 3. R n+1 = R n + tω R n. (b) Show that the updated rotation matrix computed from this update is not orthogonal.
CS205B / CME306 Homework 3 Rotation Matrices 1. The ODE that describes rigid body evolution is given by R = ω R. (a) Write down the forward Euler update for this equation. R n+1 = R n + tω R n (b) Show
More informationHomotopy perturbation method for the Wu-Zhang equation in fluid dynamics
Journal of Physics: Conference Series Homotopy perturbation method for the Wu-Zhang equation in fluid dynamics To cite this article: Z Y Ma 008 J. Phys.: Conf. Ser. 96 08 View the article online for updates
More informationThe tanh - coth Method for Soliton and Exact Solutions of the Sawada - Kotera Equation
Volume 117 No. 13 2017, 19-27 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu The tanh - oth Method for Soliton and Exat Solutions of the Sawada -
More informationInvariant Sets and Exact Solutions to Higher-Dimensional Wave Equations
Commun. Theor. Phys. Beijing, China) 49 2008) pp. 9 24 c Chinese Physical Society Vol. 49, No. 5, May 5, 2008 Invariant Sets and Exact Solutions to Higher-Dimensional Wave Equations QU Gai-Zhu, ZHANG Shun-Li,,2,
More informationPeriodic, hyperbolic and rational function solutions of nonlinear wave equations
Appl Math Inf Sci Lett 1, No 3, 97-101 (013 97 Applied Mathematics & Information Sciences Letters An International Journal http://dxdoiorg/101785/amisl/010307 Periodic, hyperbolic and rational function
More informationSymmetry Classification of KdV-Type Nonlinear Evolution Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 1, 125 130 Symmetry Classification of KdV-Type Nonlinear Evolution Equations Faruk GÜNGÖR, Victor LAHNO and Renat ZHDANOV Department
More informationarxiv: v1 [math-ph] 25 Jul Preliminaries
arxiv:1107.4877v1 [math-ph] 25 Jul 2011 Nonlinear self-adjointness and conservation laws Nail H. Ibragimov Department of Mathematics and Science, Blekinge Institute of Technology, 371 79 Karlskrona, Sweden
More information