Application of Euler and Homotopy Operators to Integrability Testing

Size: px
Start display at page:

Download "Application of Euler and Homotopy Operators to Integrability Testing"

Transcription

1 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. home/whereman/ The Ninth IMACS International Conference on Nonlinear Evolution Equations and Wave Phenomena: Computation and Theory Athens, Georgia Thursday, April 2, 2015, 4:40p.m.

2 Acknowledgements Douglas Poole (Chadron State College, Nebraska) Ünal Göktaş (Turgut Özal University, Turkey) Mark Hickman (Univ. of Canterbury, New Zealand) Bernard Deconinck (Univ. of Washington, Seattle) Several undergraduate and graduate students Research supported in part by NSF under Grant No. CCF This presentation was made in TeXpower

3 Outline Motivation of the Research Part I: Continuous Case Calculus-based formulas for the continuous homotopy operator (in multi-dimensions) See: Anco & Bluman (2002); Hereman et al. (2007); Olver (1986) Symbolic integration by parts and inversion of the total divergence operator Application: symbolic computation of conservation laws of nonlinear PDEs in multiple space dimensions

4 Part II: Discrete Case Simple formula for the homotopy operator See: Hereman et al. (2004); Hydon & Mansfield (2004); Mansfield & Hydon (2002) Symbolic summation by parts and inversion of the forward difference operator Analogy: continuous and discrete formulas Application: symbolic computation of conservation laws of nonlinear DDEs Conclusions and Future Work

5 Motivation of the Research Conservation Laws for Nonlinear PDEs System of evolution equations of order M u t = F(u (M) (x)) with u = (u, v, w,...) and x = (x, y, z). Conservation law in (1+1)-dimensions D t ρ + D x J = 0 where = means evaluated on the PDE. Conserved density ρ and flux J.

6 Conservation law in (2+1)-dimensions D t ρ + J = D t ρ + D x J 1 + D y J 2 = 0 Conserved density ρ and flux J = (J 1, J 2 ). Conservation law in (3+1)-dimensions D t ρ + J = D t ρ + D x J 1 + D y J 2 + D z J 3 = 0 Conserved density ρ and flux J = (J 1, J 2, J 3 ).

7 Notation Computations on the Jet Space Independent variables: t (time), x = (x, y, z) Dependent variables u = (u (1), u (2),..., u (j),..., u (N) ) In examples: u = (u, v, w, θ, h,...) Partial derivatives u kx = k u x k, u kx ly = k+l u x k y l, etc. Examples: u xxxxx = u 5x = 5 u x 5 u xx yyyy = u 2x 4y = Differential functions 6 u x 2 y 4 Example: f = uvv x + x 2 u 3 xv x + u x v xx

8 Total derivatives: D t, D x, D y,... Example: Let f = uvv x + x 2 u 3 xv x + u x v xx Then D x f = f x + u f x u + u xx +v x f v + v xx f u x f + v xxx v x f v xx = 2xu 3 xv x + u x (vv x ) + u xx (3x 2 u 2 xv x + v xx ) +v x (uv x ) + v xx (uv + x 2 u 3 x) + v xxx (u x ) = 2xu 3 xv x + vu x v x + 3x 2 u 2 xv x u xx + u xx v xx +uv 2 x + uvv xx + x 2 u 3 xv xx + u x v xxx

9 Example 1: The Zakharov-Kuznetsov Equation u t + αuu x + β(u xx + u yy ) x = 0 models ion-sound solitons in a low pressure uniform magnetized plasma. Conservation laws: ) ( ) ) D t (u + D α2 x u2 + βu xx + D y (βu xy = 0 D t (u 2) ( ) + D 2α x 3 u3 β(u 2 x u2 y ) + 2βu(u xx + u yy ) ( ) + D y 2βu x u y = 0

10 More conservation laws (ZK equation): ) ( D t (u 3 3β α (u2 x + u 2 y) + D x 3u 2 ( α 4 u2 + βu xx ) 6βu(u 2 x + u 2 y) + 3β2 α (u2 xx u 2 yy) 6β2 + D y ( 3βu 2 u xy + 6β2 α u xy(u xx + u yy ) ) α (u x(u xxx + u xyy ) + u y (u xxy + u yyy )) ) = 0 ( ) D t tu 2 2 α xu + D x (t( 2α 3 u3 β(u 2 x u2 y ) + 2βu(u xx + u yy )) ) ( ) x(u 2 + 2β α u xx) + 2β α u x + D y 2β(tu x u y + 1 α xu xy) = 0

11 Example 2: Shallow Water Wave Equations P. Dellar, Phys. Fluids 15 (2003) u t + (u )u + 2 Ω u + (θh) 1 2 h θ = 0 θ t + u ( θ) = 0 h t + (uh) = 0 In components: u t + uu x + vu y 2 Ωv hθ x + θh x = 0 v t + uv x + vv y + 2 Ωu hθ y + θh y = 0 θ t + uθ x + vθ y = 0 h t + hu x + uh x + hv y + vh y = 0 where u(x, y, t), θ(x, y, t) and h(x, y, t)

12

13 First few densities-flux pairs of SWW system: ρ (1) = h ρ (2) = h θ J (1) = h J (2) = h θ ρ (3) = h θ 2 J (3) = h θ 2 ρ (4) = h (u 2 + v 2 + hθ) ρ (5) = θ (2Ω + v x u y ) J (5) = 1 2 θ J (4) = h u v u v u v u (u2 + v 2 + 2hθ) v (v 2 + u 2 + 2hθ) 4Ωu 2uu y + 2uv x hθ y 4Ωv + 2vv x 2vu y + hθ x

14 Generalizations: D t ( f(θ)h ) + Dx ( f(θ)hu ) + Dy ( f(θ)hv ) = 0 ( D t g(θ)(2ω + vx u x ) ) ( ) + D 1 x 2 g(θ)(4ωu 2uu y + 2uv x hθ y ) ( ) + D 1 y 2 g(θ)(4ωv 2u yv + 2vv x + hθ x ) = 0 for any functions f(θ) and g(θ)

15 Conservation Laws for Nonlinear Differential-Difference Equations (DDEs) System of DDEs u n =F(, u n 1, u n, u n+1, ) Conservation law in (1 + 1) dimensions D t ρ n + J n = D t ρ n + J n+1 J n = 0 conserved density ρ n and flux J n Example: Toda lattice u n = v n 1 v n v n = v n (u n u n+1 )

16 First few densities-flux pairs for Toda lattice: ρ (0) n =ln(v n ) J n (0) =u n ρ (1) n =u n J n (1) =v n 1 ρ (2) n = 1 2 u2 n + v n J n (2) =u n v n 1 ρ (3) n = 1 3 u3 n + u n (v n 1 + v n ) J n (3) =u n 1 u n v n 1 + vn 1 2

17 PART I: CONTINUOUS CASE Problem Statement Continuous case in 1D: Example: For u(x) and v(x) f =8v x v xx u 3 x sin u +2u x u xx cos u 6vv x cos u +3u x v 2 sin u Question: Can the expression be integrated? If yes, find F = f dx (so, f = D x F ) Result (by hand): F = 4 vx 2 + u2 x cos u 3 v2 cos u

18 Continuous case in 2D or 3D: Example: For u(x, y) and v(x, y) f = u x v y u 2x v y u y v x + u xy v x Question: Is there an F so that f = Div F? If yes, find F. Result (by hand): F = (uv y u x v y, uv x + u x v x ) Can this be done without integration by parts? Can the computation be reduced to a single integral in one variable?

19 Tools from the Calculus of Variations Definition: A differential function f is a exact iff f = Div F. Special case (1D): f = D x F. Question: How can one test that f = Div F? Theorem (exactness test): f = Div F iff L u (j) (x) f 0, j = 1, 2,..., N. N is the number of dependent variables. The Euler operator annihilates divergences

20 Euler operator in 1D (variable u(x)): L u(x) = M ( D x ) k k=0 = u D x u kx u x + D 2 x u xx D 3 x Euler operator in 2D (variable u(x, y)): L u(x,y) = M x k=0 M y ( D x ) k ( D y ) l l=0 = u D x + D 2 x u x D y u xx +D x D y u y +D 2 y u xy u kx ly u yy D 3 x u xxx + u xxx

21 Application: Testing Exactness Continous Case Example: f =8v x v xx u 3 x sin u +2u x u xx cos u 6vv x cos u +3u x v 2 sin u where u(x) and v(x) f is exact After integration by parts (by hand): F = f dx = 4 vx 2 + u 2 x cos u 3 v 2 cos u

22 Exactness test with Euler operator: f =8v x v xx u 3 x sin u +2u x u xx cos u 6vv x cos u +3u x v 2 sin u L u(x) f = f u D x L v(x) f = f v D x f u x + D 2 x f v x + D 2 x f u xx 0 f v xx 0

23 Question: How can one compute F = Div 1 f? Theorem (integration by parts): In 1D: If f is exact then F = D 1 x f = f dx = H u(x) f In 2D: If f is a divergence then F = Div 1 f = (H (x) u(x,y) f, H(y) u(x,y) f) The homotopy operator inverts total derivatives and divergences!

24 Peter Olver s Book

25 Homotopy Concept in Olver s Book

26 Homotopy Formula in Olver s Book

27 Zoom into Homotopy Formula in Olver s Book

28

29 Homotopy Operator in 1D (variable x): H u(x) f = 1 0 N (I u (j)f)[λu] dλ λ j=1 with integrand I u (j)f = M x (j) k=1 k 1 i=0 ix ( D x) k (i+1) u (j) f u (j) kx (I u (j)f)[λu] means that in I u (j)f one replaces u λu, u x λu x, etc. More general: u λ(u u 0 ) + u 0 u x λ(u x u x 0) + u x 0 etc.

30 Homotopy Operator in 2D (variables x and y): H (x) u(x,y) f = 1 0 N j=1 (I (x) u (j) f)[λu] dλ λ H (y) u(x,y) f = 1 0 N j=1 (I (y) u (j) f)[λu] dλ λ where for dependent variable u(x, y) M x M y ( k 1 l i+j I u (x) f = i u ix jy k=1 l=0 i=0 j=0 ( D x ) k i 1 ( D y ) l j ) )( k+l i j 1 k i 1 ( k+l ) k f u kx ly )

31 Application of Homotopy Operator in 1D Example: f =8v x v xx u 3 x sin u +2u xu xx cos u 6vv x cos u +3u x v 2 sin u Goal: Find F = 4 v 2 x + u 2 x cos u 3 v 2 cos u Easy to verify: f = D x F Compute I u f = u f u x + (u x I ud x ) f u xx = uu 2 x sin u + 3uv 2 sin u + 2u 2 x cos u

32 Similarly, I v f = v f v x + (v x I vd x ) f v xx Finally, = 6v 2 cos u + 8v 2 x 1 F = H u(x) f = (I u f + I v f) [λu] dλ 0 λ 1 ( = 3λ 2 uv 2 sin(λu) λ 2 uu 2 x sin(λu) 0 +2λu 2 x cos(λu) 6λv2 cos(λu) + 8λv 2 x ) dλ = 4v 2 x + u2 x cos u 3v2 cos u

33 Application of Homotopy Operator in 2D Example: f = u x v y u xx v y u y v x + u xy v x By hand: F = (uv y u x v y, uv x + u x v x ) Easy to verify: f = Div F Compute I u (x) f = u f + (u x I ud x ) f u x ( ) + 1 f 2 u yi 1 2 ud y u xy u xx = uv y u yv x u x v y uv xy

34 Similarly, I v (x) f = v f = u y v + u xy v v x Hence, F 1 = H (x) u(x,y) f = 1 = 1 0 λ ( I (x) u ) f + I v (x) f [λu] dλ 0 λ ( uv y u yv x u x v y uv xy u y v + u xy v = 1 2 uv y u yv x 1 2 u xv y uv xy 1 2 u yv u xyv ) dλ

35 Analogously, F 2 = H (y) u(x,y) f = 1 = 1 0 ( λ ( ) I u (y) f + I v (y) f [λu] dλ 0 λ ) ) ( uv x 12 uv xx + 12 u xv x + λ (u x v u xx v) dλ = 1 2 uv x 1 4 uv xx u xv x u xv 1 2 u xxv So, F= 1 4 2uv y + u y v x 2u x v y + uv xy 2u y v + 2u xy v 2uv x uv xx + u x v x + 2u x v 2u xx v

36 Let K= F F then K= 1 4 2uv y u y v x 2u x v y uv xy + 2u y v 2u xy v 2uv x + uv xx + 3u x v x 2u x v + 2u xx v then Div K = 0 Also, K = (D y φ, D x φ) with φ = 1 4 (2uv uv x 2u x v) (curl in 2D) After removing the curl term K: F = F + K = (uv y u x v y, uv x + u x v x ) Avoid curl terms algorithmically!

37 Why does this work? Sketch of Derivation and Proof (in 1D with variable x, and for one component u) Definition: Degree operator M Mf = M i=0 u ix f u ix = u f u +u x f u x +u 2x f u 2x + +u Mx f u Mx f is of order M in x Example: f = u p u q xu r 3x (p, q, r non-negative integers) g = Mf = 3 i=0 u ix f u ix = (p + q + r) u p u q xu r 3x Application of M computes the total degree

38 Theorem (inverse operator) M 1 g(u) = 1 0 g[λu] dλ λ Proof: d M dλ g[λu] = i=0 g[λu] dλu ix λu ix dλ = 1 λ M i=0 Integrate both sides with respect to λ u ix g[λu] u ix = 1 λ Mg[λu] 1 0 d dλ g[λu] dλ = g[λu] λ=1 λ=0 = 1 0 = g(u) g(0) Mg[λu] dλ λ = M 1 0 g[λu] dλ λ Assuming g(0) = 0, M 1 g(u) = 1 0 g[λu] dλ λ

39 Example: If g(u) = (p + q + r) u p u q xu r 3x, then g[λu] = (p + q + r)λ p+q+r u p u q xu r 3x Hence, M 1 g = 1 0 (p + q + r) λ p+q+r 1 u p u q x ur 3x dλ = u p u q xu r 3x λ p+q+r λ=1 = λ=0 up u q xu r 3x

40 Theorem: If f is an exact differential function, then F = Dx 1 f = f dx = H u(x) f Proof: Multiply L u(x) f = M ( D x ) k k=0 f u kx by u to restore the degree. Split off u f u. Split off u x f u x. Lastly, split off u Mx Integrate by parts. Repeat the process. f u Mx.

41 ul u(x) f = u = u f u D x = u f u + u x =... M ( D x ) k k=0 u f u kx M ( D x ) k 1 k=1 f u x D x M +u x ( D x ) k 2 k=2 u f u kx f u kx M ( D x ) k 1 k=1 M + u x ( D x ) k 1 f k=1 u kx M + u 2x ( D x ) k 2 k=2 f u kx f u kx

42 f u x u Mx = u f u + u x M D x u ( D x ) k 1 = k= u (M 1)x M M i=0 = Mf D x = 0 u ix f u ix D x M 1 i=0 k=m M 1 u ix f u Mx f u kx + u x ( D x ) k M i=0 M k=i+1 u ix M k=i+1 M ( D x ) k 2 k=2 f u kx ( D x ) k (i+1) ( D x ) k (i+1) f u kx f u kx f u kx

43

44 Mf = D x M 1 i=0 u ix M k=i+1 ( D x ) k (i+1) f u kx Apply M 1 and use M 1 D x = D x M 1. M 1 M f = D x M 1 ( D x ) k (i+1) u ix i=0 k=i+1 f u kx Apply Dx 1 and use the formula for M 1. 1 M 1 M F = Dx 1 f = u ix ( D x ) k (i+1) = 1 0 = H u(x) f k=1 0 i=0 k=i+1 M k 1 u ix ( D x ) k (i+1) i=0 f u kx f u kx [λu] dλ λ [λu] dλ λ

45

46 Application 1: Zakharov-Kuznetsov Equation Computation of Conservation Laws u t + αuu x + β(u xx + u yy ) x = 0 Step 1: Compute the dilation invariance ZK equation is invariant under scaling symmetry (t, x, y, u) ( t λ 3, x λ, y λ, λ2 u)= ( t, x, ỹ, ũ) λ is an arbitrary parameter. Hence, the weights of the variables are W (u) = 2, W (D t ) = 3, W (D x ) = 1, W (D y ) = 1.

47 A conservation law is invariant under the scaling symmetry of the PDE. W (u) = 2, W (D t ) = 3, W (D x ) = 1, W (D y ) = 1. For example, ) ( D t (u 3 3β α (u2 x +u 2 y) + D x 3u 2 ( α 4 u2 +βu xx ) 6βu(u 2 x + u 2 y) ) + 3β2 α (u2 xx u 2 xy) 6β2 α (u x(u xxx +u xyy ) + u y (u xxy +u yyy )) ( ) + D y 3βu 2 u xy + 6β2 α u xy(u xx + u yy ) = 0 Rank (ρ) = 6, Rank (J) = 8. Rank (conservation law) = 9.

48 Compute the density of selected rank, say, 6. Step 2: Construct the candidate density For example, construct a density of rank 6. Make a list of all terms with rank 6: {u 3, u 2 x, uu xx, u 2 y, uu yy, u x u y, uu xy, u 4x, u 3xy, u 2x2y, u x3y, u 4y } Remove divergences and divergence-equivalent terms. Candidate density of rank 6: ρ = c 1 u 3 + c 2 u 2 x + c 3 u 2 y + c 4 u x u y

49 Step 3: Compute the undetermined coefficients Compute D t ρ = ρ t + ρ (u)[u t ] = ρ t + = M x k=0 M y l=0 ρ u kx ly D k x Dl y u t ( 3c 1 u 2 I + 2c 2 u x D x + 2c 3 u y D y + c 4 (u y D x + u x D y ) ) Substitute u t = (αuu x + β(u xx + u yy ) x. ) u t

50 E = D t ρ = 3c 1 u 2 (αuu x + β(u xx + u xy ) x ) + 2c 2 u x (αuu x + β(u xx + u yy ) x ) x + 2c 3 u y (αuu x + β(u xx + u yy ) x ) y + c 4 (u y (αuu x + β(u xx + u yy ) x ) x + u x (αuu x + β(u xx + u yy ) x ) y ) Apply the Euler operator (variational derivative) M x M y L u(x,y) E = ( D x ) k ( D y ) l E k=0 l=0 u kx ly ( = 2 (3c 1 β + c 3 α)u x u yy + 2(3c 1 β + c 3 α)u y u xy ) +2c 4 αu x u xy +c 4 αu y u xx +3(3c 1 β +c 2 α)u x u xx 0

51 Solve a parameterized linear system for the c i : 3c 1 β + c 3 α = 0, c 4 α = 0, 3c 1 β + c 2 α = 0 Solution: c 1 = 1, c 2 = 3β α, c 3 = 3β α, c 4 = 0 Substitute the solution into the candidate density ρ = c 1 u 3 + c 2 u 2 x + c 3 u 2 y + c 4 u x u y Final density of rank 6: ρ = u 3 3β α (u2 x + u 2 y)

52 Step 4: Compute the flux Use the homotopy operator to invert Div: ( ) J=Div 1 E = H (x) u(x,y) E, H(y) u(x,y) E where with I (x) u E = H (x) M x k=1 u(x,y) E = 1 M y l=0 ( k 1 i=0 0 l j=0 (I u (x) E)[λu] dλ λ ( i+j i u ix jy ( D x ) k i 1 ( D y ) l j ) )( k+l i j 1 k i 1 ( k+l ) k E u kx ly Similar formulas for H (y) u(x,y) E and I(y) u E. )

53 Let A = αuu x + β(u xxx + u xyy ) so that E = 3u 2 A 6β α u xa x 6β α u ya y Then, ( J = H (x) u(x,y) ) E, H(y) u(x,y) E = ( 3α 4 u4 + βu 2 (3u xx + 2u yy ) 2βu(3u 2 x + u 2 y) + 3β2 4α u(u 2x2y + u 4y ) β2 α u x( 7 2 u xyy + 6u xxx ) β2 α u y(4u xxy u yyy) + β2 α (3u2 xx u2 xy u2 yy) + 5β2 4α u xxu yy, βu 2 u xy 4βuu x u y 3β2 4α u(u x3y + u 3xy ) β2 4α u x(13u xxy + 3u yyy ) 5β2 4α u y(u xxx + 3u xyy ) + 9β2 4α u xy(u xx + u yy ) )

54 Application 2: Shallow Water Wave Equations Computation of Conservation Laws Quick Recapitulation Conservation law in (2+1) dimensions D t ρ + J = D t ρ + D x J 1 + D y J 2 = 0 conserved density ρ and flux J = (J 1, J 2 ) Example: Shallow water wave (SWW) equations u t + uu x + vu y 2 Ωv hθ x + θh x = 0 v t + uv x + vv y + 2 Ωu hθ y + θh y = 0 θ t + uθ x + vθ y = 0 h t + hu x + uh x + hv y + vh y = 0

55 Typical density-flux pair: ρ (5) = θ (v x u y + 2Ω) J (5) = 1 2 θ 4Ωu 2uu y + 2uv x hθ y 4Ωv + 2vv x 2vu y + hθ x

56 Step 1: Construct the form of the density The SWW equations are invariant under the scaling symmetries (x, y, t, u, v, θ, h, Ω) (λ 1 x, λ 1 y, λ 2 t, λu, λv, λθ, λh, λ 2 Ω) and (x, y, t, u, v, θ, h, Ω) (λ 1 x, λ 1 y, λ 2 t, λu, λv, λ 2 θ, λ 0 h, λ 2 Ω) Construct a candidate density, for example, ρ = c 1 Ωθ + c 2 u y θ + c 3 v y θ + c 4 u x θ + c 5 v x θ which is scaling invariant under both symmetries.

57 Step 2: Determine the constants c i Compute E = D t ρ and remove time derivatives E = ( ρ u x u tx + ρ u y u ty + ρ v x v tx + ρ v y v ty + ρ θ θ t) = c 4 θ(uu x + vu y 2Ωv hθ x + θh x ) x + c 2 θ(uu x + vu y 2Ωv hθ x + θh x ) y + c 5 θ(uv x + vv y + 2Ωu hθ y + θh y ) x + c 3 θ(uv x + vv y + 2Ωu hθ y + θh y ) y + (c 1 Ω + c 2 u y + c 3 v y + c 4 u x + c 5 v x )(uθ x + vθ y ) Require that L u(x,y) E = L v(x,y) E = L θ(x,y) E = L h(x,y) E 0.

58 Solution: c 1 = 2, c 2 = 1, c 3 = c 4 = 0, c 5 = 1 gives ρ = ρ (5) = θ (2Ω u y + v x ) Step 3: Compute the flux J E = θ(u x v x + uv xx + v x v y + vv xy + 2Ωu x θ xh y u x u y uu xy u y v y u yy v +2Ωv y 1 2 θ yh x ) + 2Ωuθ x + 2Ωvθ y uu y θ x u y vθ y + uv x θ x + vv x θ y Apply the 2D homotopy operator: J = (J 1, J 2 ) = Div 1 E = (H (x) u(x,y) E, H(y) u(x,y) E)

59 Compute I u (x) E = u E + u x ( ) 1 E 2 u yi 1 2 ud y u xy Similarly, compute = uv x θ + 2Ωuθ u2 θ y uu y θ I (x) v E = vv y θ v2 θ y + uv x θ I (x) θ E = 1 2 θ2 h y + 2Ωuθ uu y θ + uv x θ I (x) h E = 1 2 θθ yh

60 Next, J 1 = H (x) = = u(x,y) E 1 ( I u (x) ) E + I(x) v E + I (x) θ E + I (x) h E [λu] dλ λ ( ( 4λΩuθ + λ 2 3uv x θ u2 θ y 2uu y θ + vv y θ v2 θ y θ2 h y 1 2 θθ yh) ) dλ = 2Ωuθ 2 3 uu yθ+ uv x θ+ 1 3 vv yθ+ 1 6 u2 θ y v2 θ y 1 6 hθθ y h yθ 2

61 Analogously, J 2 = H (y) u(x,y) E Hence, = 2Ωvθ vv xθ vu y θ 1 3 uu xθ 1 6 u2 θ x 1 6 v2 θ x hθθ x 1 6 h xθ 2 J = J (5) = 1 12Ωuθ 4uu yθ+6uv x θ+2vv y θ+(u 2 +v 2 )θ y hθθ y +h y θ Ωvθ+4vv x θ 6vu y θ 2uu x θ (u 2 +v 2 )θ x +hθθ x h x θ 2

62 There are curl terms in J Indeed, subtract K where Div K = 0 Here K= 1 6 (2uu yθ+2vv y θ+u 2 θ y +v 2 θ y +2hθθ y +h y θ 2 ) 2vv x θ+2uu x θ+u 2 θ x +v 2 θ x +2hθθ x +h x θ 2 Note that K = (D y φ, D x φ) with φ = (hθ 2 + u 2 θ + v 2 θ) (i.e., curl in 2D).

63 After removing the curl term K, J = J (5) = 1 2 θ 4Ωu 2uu y + 2uv x hθ y 4Ωv + 2vv x 2vu y + hθ x

64 PART II: DISCRETE CASE Problem Statement Discrete case in 1D: Example: f n = u n u n+1 v n vn+u 2 n+1 u n+2 v n+1 +vn+1+u 2 n+3 v n+2 u n+1 v n Question: Can the expression be summed by parts? If yes, find F n = 1 f n (so, f n = F n = F n+1 F n ) Result (by hand): F n =vn 2 + u n u n+1 v n + u n+1 v n + u n+2 v n+1 How can this be done algorithmically? Can this be done as in the continuous case?

65 Tools from the Discrete Calculus of Variations Definitions: D is the up-shift (forward or right-shift) operator DF n = F n+1 = F n n n+1 D 1 the down-shift (backward or left-shift) operator D 1 F n = F n 1 = F n n n 1 = D I is the forward difference operator F n = (D I)F n = F n+1 F n Problem to be solved: Given f n. Find F n = 1 f n (so f n = F n = F n+1 F n )

66 Analogy Continuous & Discrete Cases Euler Operators Continuous Case Discrete Case L u(x) = M ( D x ) k k=0 u k x L un = M k=0 = u n D k M k=0 u n+k D k

67 Analogy Continuous & Discrete Cases Homotopy Operators & Integrands Continuous Case Discrete Case 1 H u(x) f= 0 N (I u (j)f)[λu] dλ λ H u n f n = j=1 1 0 N j=1 (I u (j) n f n )[λu n ] dλ λ I u (j)f = k 1 i=0 u (j) M (j) k=1 ix ( D x) k (i+1) f u (j) kx I (j) u f n = n k i=1 M (j) k=1 D i u (j) n+k f n u (j) n+k

68 Euler Operators Side by Side Continuous Case (for component u) L u = u D x u x + D 2 x u 2x D 3 x Discrete Case (for component u n ) L un = u n + D 1 u n+1 + D 2 u n+2 + D 3 = u n (I + D 1 + D 2 + D 3 + ) u 3x + u n+3 +

69 Homotopy Operators Side by Side Continuous Case (for components u and v) with I u f = H u(x) f = M (1) k=1 1 0 i=0 (I u f + I v f) [λu] dλ λ k 1 u ix ( D x ) k (i+1) f u kx and I v f = M (2) k=1 k 1 i=0 v ix ( D x ) k (i+1) f v kx

70 Discrete Case (for components u n and v n ) with H un f n = I un f n = 1 0 M (1) k=1 (I un f n + I vn f n ) [λu n ] dλ λ k i=1 D i u n+k f n u n+k and I vn f n = M (2) k=1 k i=1 D i v n+k f n v n+k

71 Analogy of Definitions & Theorems Continuous Case (PDE) Semi-discrete Case (DDE) u t =F(u, u x, u 2x, ) u n =F(, u n 1, u n, u n+1, ) D t ρ + D x J = 0 D t ρ n + J n = 0 Definition: f n is exact iff f n = F n = F n+1 F n Theorem (exactness test): f n = F n iff L un f n 0 Theorem (summation with homotopy operator): If f n is exact then F n = 1 f n = H un (f n )

72 Example: Testing Exactness Discrete Case f n = u n u n+1 v n v 2 n+u n+1 u n+2 v n+1 +v 2 n+1+u n+3 v n+2 u n+1 v n f n is exact After summation by parts (done by hand): F n = v 2 n + u n u n+1 v n + u n+1 v n + u n+2 v n+1 Easy to verify: f n = F n = F n+1 F n

73 Exactness test with Euler operator: For component u n (highest shift 3): L un f n = ( I + D 1 + D 2 + D 3) f n u n = u n+1 v n u n 1 v n 1 + u n+1 v n v n 1 Similarly, 0 + u n 1 v n 1 + v n 1 L vn f n = v n ( I + D 1) f n = u n u n+1 + 2v n u n u n+1 2v n 0

74 Application of Discrete Homotopy Operator Example: f n = u n u n+1 v n v 2 n+u n+1 u n+2 v n+1 +v 2 n+1+u n+3 v n+2 u n+1 v n Here, M (1) = 3 and M (2) = 2. Compute I un f n = (D 1 ) u n+1 f n u n+1 + (D 1 + D 2 ) u n+2 f n u n+2 + (D 1 + D 2 + D 3 ) u n+3 f n u n+3 = 2u n u n+1 v n + u n+1 v n + u n+2 v n+1

75 I vn f n = (D 1 ) v n+1 f n v n+1 Finally, + (D 1 + D 2 ) v n+2 f n v n+2 = u n u n+1 v n + 2v 2 n + u n+1 v n + u n+2 v n+1 F n = = (I un f n + I vn f n ) [λu n ] dλ λ ) (2λvn 2 + 3λ 2 u n u n+1 v n + 2λu n+1 v n + 2λu n+2 v n+1 dλ = v 2 n + u n u n+1 v n + u n+1 v n + u n+2 v n+1

76 Application: Computation of Conservation Laws System of DDEs u n =F(, u n 1, u n, u n+1, ) Conservation law D t ρ n + J n = 0 conserved density ρ n and flux J n

77 Example: Toda lattice u n = v n 1 v n v n = v n (u n u n+1 ) Typical density-flux pair: ρ (3) n = 1 3 u3 n + u n (v n 1 + v n ) J (3) n = u n 1 u n v n 1 + v 2 n 1

78 Computation Conservation Laws for Toda Lattice Step 1: Construct the form of the density The Toda lattice is invariant under scaling symmetry (t, u n, v n ) (λ 1 t, λu n, λ 2 v n ) Construct a candidate density, for example, ρ n = c 1 u 3 n + c 2 u n v n 1 + c 3 u n v n which is scaling invariant under the symmetry

79 Step 2: Determine the constants c i Compute E n = D t ρ n and remove time derivatives E n = (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 u n 1 u n v n 1 + c 2 v 2 n 1 c 3 u n u n+1 v n c 3 v 2 n Compute Ẽn = DE n to remove negative shift n 1 Require that L un Ẽ n = L vn Ẽ n 0 Solution: c 1 = 1 3, c 2 = c 3 = 1 gives ρ n = ρ (3) n = 1 3 u3 n + u n (v n 1 + v n )

80 Step 3: Compute the flux J n Ẽ n = DE n = u n u n+1 v n + vn 2 u n+1 u n+2 v n+1 vn+1 2 Apply the homotopy operator J n = DJ n = 1 (Ẽn) = H un (Ẽn) Compute I un Ẽ n = (D 1 Ẽn ) u n+1 + (D 1 + D 2 ) u n+2 u n+1 = 2u n u n+1 v n Ẽn u n+2 Likewise, I vn Ẽ n = (D 1 ) v n+1 Ẽn v n+1 = (u n u n+1 v n + 2v 2 n)

81 Next, compute J n = = ) (I un Ẽ n + I vn Ẽ n [λu n ] dλ λ (3λ 2 u n u n+1 v n + 2λv 2 n) dλ = u n u n+1 v n + vn 2 Finally, backward shift J n = D 1 ( J n ) given J n = J (3) n = u n 1 u n v n 1 + v 2 n 1

82 Conclusions and Future Work The power of Euler and homotopy operators: Testing exactness Integration by parts: D 1 x and Div 1 Integration of non-exact expressions Example: f = u x v + uv x + u 2 u xx fdx = uv + u 2 u xx dx Use other homotopy formulas (moving terms amongst the components of the flux; prevent curl terms)

83 Homotopy operator approach pays off for computing irrational fluxes Example: short pulse equation (nonlinear optics) u xt = u + (u 3 ) xx = u + 6uu 2 x + 3u 2 u xx with non-polynomial conservation law ( ) D t 1 + 6u 2 x D x (3u 2 ) 1 + 6u 2 x = 0 Continue the implementation in Mathematica Software: whereman

84 Thank You

85 Publications 1. D. Poole and W. Hereman, Symbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions, J. Symb. Comp. 46(12), (2011). 2. W. Hereman, P. J. Adams, H. L. Eklund, M. S. Hickman, and B. M. Herbst, Direct Methods and Symbolic Software for Conservation Laws of Nonlinear Equations, In: Advances of Nonlinear Waves and Symbolic Computation, Ed.: Z. Yan, Nova Science Publishers, New York (2009), Chapter 2, pp

86 3. W. Hereman, M. Colagrosso, R. Sayers, A. Ringler, B. Deconinck, M. Nivala, and M. S. Hickman, Continuous and Discrete Homotopy Operators and the Computation of Conservation Laws. In: Differential Equations with Symbolic Computation, Eds.: D. Wang and Z. Zheng, Birkhäuser Verlag, Basel (2005), Chapter 15, pp W. Hereman, B. Deconinck, and L. D. Poole, Continuous and discrete homotopy operators: A theoretical approach made concrete, Math. Comput. Simul. 74(4-5), (2007).

87 5. W. Hereman, Symbolic computation of conservation laws of nonlinear partial differential equations in multi-dimensions, Int. J. Quan. Chem. 106(1), (2006).

Symbolic 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 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 information

Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations

Symbolic 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 information

Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations

Symbolic 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 information

Continuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman

Continuous 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 information

Symbolic 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 . 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 information

Symbolic Computation of Recursion Operators of Nonlinear Differential-Difference Equations

Symbolic 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 information

Symbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions

Symbolic 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 information

Symbolic computation of conservation laws for nonlinear partial differential equations in multiple space dimensions

Symbolic 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 information

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

The 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

The 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 information

Symbolic 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 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 information

arxiv:nlin/ v1 [nlin.si] 16 Nov 2003

arxiv: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 information

Alexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011

Alexei 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 information

Solving Nonlinear Wave Equations and Lattices with Mathematica. Willy Hereman

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

More information

Exact solutions of the two-dimensional Boussinesq and dispersive water waves equations

Exact solutions of the two-dimensional Boussinesq and dispersive water waves equations Advances in Fluid Mechanics VIII 293 Exact solutions of the two-dimensional Boussinesq and dispersive water waves equations F. P. Barrera 1, T. Brugarino 2 & F. Montano 1 1 Dip. di Ingegneria dei Trasporti,

More information

Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.

Conservation 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 information

SYMBOLIC 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 . 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 information

arxiv:nlin/ v1 [nlin.si] 16 Nov 2003

arxiv: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 information

Algorithmic 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. 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 information

. p.1/31. Department Mathematics and Statistics Arizona State University. Don Jones

. p.1/31. Department Mathematics and Statistics Arizona State University. Don Jones Modified-Truncation Finite Difference Schemes. p.1/31 Department Mathematics and Statistics Arizona State University Don Jones dajones@math.asu.edu . p.2/31 Heat Equation t u = λu xx + f . p.2/31 Heat

More information

Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem

Construction 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 information

Summer 2017 MATH Solution to Exercise 5

Summer 2017 MATH Solution to Exercise 5 Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a

More information

Homotopy perturbation method for the Wu-Zhang equation in fluid dynamics

Homotopy 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 information

A 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 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 information

Direct construction method for conservation laws of partial differential equations Part I: Examples of conservation law classifications

Direct 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 information

Symbolic 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 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 information

Final Exam May 4, 2016

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

More information

Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders

Breaking 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 information

Inductive and Recursive Moving Frames for Lie Pseudo-Groups

Inductive and Recursive Moving Frames for Lie Pseudo-Groups Inductive and Recursive Moving Frames for Lie Pseudo-Groups Francis Valiquette (Joint work with Peter) Symmetries of Differential Equations: Frames, Invariants and Applications Department of Mathematics

More information

PARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007

PARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007 PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.

More information

Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy

Prolongation 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 information

RESEARCH ARTICLE. A symbolic algorithm for computing recursion operators of nonlinear PDEs

RESEARCH 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 information

Gradient, Divergence and Curl in Curvilinear Coordinates

Gradient, Divergence and Curl in Curvilinear Coordinates Gradient, Divergence and Curl in Curvilinear Coordinates Although cartesian orthogonal coordinates are very intuitive and easy to use, it is often found more convenient to work with other coordinate systems.

More information

3 Applications of partial differentiation

3 Applications of partial differentiation Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives

More information

Hamiltonian partial differential equations and Painlevé transcendents

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

More information

Homotopy Perturbation Method for Solving Systems of Nonlinear Coupled Equations

Homotopy 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 information

Symbolic Computation and New Soliton-Like Solutions of the 1+2D Calogero-Bogoyavlenskii-Schif Equation

Symbolic Computation and New Soliton-Like Solutions of the 1+2D Calogero-Bogoyavlenskii-Schif Equation MM Research Preprints, 85 93 MMRC, AMSS, Academia Sinica, Beijing No., December 003 85 Symbolic Computation and New Soliton-Like Solutions of the 1+D Calogero-Bogoyavlenskii-Schif Equation Zhenya Yan Key

More information

Symbolic algorithms for the Painlevé test, special solutions, and recursion operators for nonlinear PDEs.

Symbolic 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 information

Final: Solutions Math 118A, Fall 2013

Final: Solutions Math 118A, Fall 2013 Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or

More information

Symmetry 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 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 information

McGill University April 20, Advanced Calculus for Engineers

McGill University April 20, Advanced Calculus for Engineers McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student

More information

Group analysis, nonlinear self-adjointness, conservation laws, and soliton solutions for the mkdv systems

Group 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 information

Differential equations, comprehensive exam topics and sample questions

Differential equations, comprehensive exam topics and sample questions Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions

More information

The Solitary Wave Solutions of Zoomeron Equation

The 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 information

A First Course of Partial Differential Equations in Physical Sciences and Engineering

A First Course of Partial Differential Equations in Physical Sciences and Engineering A First Course of Partial Differential Equations in Physical Sciences and Engineering Marcel B. Finan Arkansas Tech University c All Rights Reserved First Draft August 29 2 Preface Partial differential

More information

Chapter 3 Second Order Linear Equations

Chapter 3 Second Order Linear Equations Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,

More information

Second Order ODEs. Second Order ODEs. In general second order ODEs contain terms involving y, dy But here only consider equations of the form

Second Order ODEs. Second Order ODEs. In general second order ODEs contain terms involving y, dy But here only consider equations of the form Second Order ODEs Second Order ODEs In general second order ODEs contain terms involving y, dy But here only consider equations of the form A d2 y dx 2 + B dy dx + Cy = 0 dx, d2 y dx 2 and F(x). where

More information

Derivatives and Integrals

Derivatives and Integrals Derivatives and Integrals Definition 1: Derivative Formulas d dx (c) = 0 d dx (f ± g) = f ± g d dx (kx) = k d dx (xn ) = nx n 1 (f g) = f g + fg ( ) f = f g fg g g 2 (f(g(x))) = f (g(x)) g (x) d dx (ax

More information

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE

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

More information

On the classification of certain curves up to projective tranformations

On the classification of certain curves up to projective tranformations On the classification of certain curves up to projective tranformations Mehdi Nadjafikhah Abstract The purpose of this paper is to classify the curves in the form y 3 = c 3 x 3 +c 2 x 2 + c 1x + c 0, with

More information

Group Actions and Cohomology in the Calculus of Variations

Group 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 information

Leplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math EECE

Leplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math EECE Leplace s Analyzing the Analyticity of Analytic Analysis Engineering Math EECE 3640 1 The Laplace equations are built on the Cauchy- Riemann equations. They are used in many branches of physics such as

More information

Lax Representations for Matrix Short Pulse Equations

Lax Representations for Matrix Short Pulse Equations Lax Representations for Matrix Short Pulse Equations Z. Popowicz arxiv:1705.04030v1 [nlin.si] 11 May 017 May 1, 017 Institute of Theoretical Physics, University of Wrocław, Wrocław pl. M. Borna 9, 50-05

More information

Extra Problems and Examples

Extra Problems and Examples Extra Problems and Examples Steven Bellenot October 11, 2007 1 Separation of Variables Find the solution u(x, y) to the following equations by separating variables. 1. u x + u y = 0 2. u x u y = 0 answer:

More information

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015 Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction

More information

The (G'/G) - Expansion Method for Finding Traveling Wave Solutions of Some Nonlinear Pdes in Mathematical Physics

The (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 information

APPLIED MATHEMATICS. Part 1: Ordinary Differential Equations. Wu-ting Tsai

APPLIED MATHEMATICS. Part 1: Ordinary Differential Equations. Wu-ting Tsai APPLIED MATHEMATICS Part 1: Ordinary Differential Equations Contents 1 First Order Differential Equations 3 1.1 Basic Concepts and Ideas................... 4 1.2 Separable Differential Equations................

More information

Lecture Notes in Mathematics. A First Course in Quasi-Linear Partial Differential Equations for Physical Sciences and Engineering Solution Manual

Lecture Notes in Mathematics. A First Course in Quasi-Linear Partial Differential Equations for Physical Sciences and Engineering Solution Manual Lecture Notes in Mathematics A First Course in Quasi-Linear Partial Differential Equations for Physical Sciences and Engineering Solution Manual Marcel B. Finan Arkansas Tech University c All Rights Reserved

More information

Moving Frame Derivation of the Fundamental Equi-Affine Differential Invariants for Level Set Functions

Moving Frame Derivation of the Fundamental Equi-Affine Differential Invariants for Level Set Functions Moving Frame Derivation of the Fundamental Equi-Affine Differential Invariants for Level Set Functions Peter. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@umn.edu http://www.math.umn.edu/

More information

Travelling Wave Solutions and Conservation Laws of Fisher-Kolmogorov Equation

Travelling Wave Solutions and Conservation Laws of Fisher-Kolmogorov Equation Gen. Math. Notes, Vol. 18, No. 2, October, 2013, pp. 16-25 ISSN 2219-7184; Copyright c ICSRS Publication, 2013 www.i-csrs.org Available free online at http://www.geman.in Travelling Wave Solutions and

More information

OVER 150 J-HOP To Q JCATS S(MJ) BY NOON TUESDAY CO-ED S LEAGUE HOLDS MEETING

OVER 150 J-HOP To Q JCATS S(MJ) BY NOON TUESDAY CO-ED S LEAGUE HOLDS MEETING F F V à y y > * y y y! F! * * F k y è 2 3 U D Y F B U Y 3 * 93 G U P B PU FF; JH H D V 50 JHP QJ (J) BY UDY P G y D; P U F P 0000 GUU VD PU F B U q F' yyy " D Py D x " By ; B x; x " P 93 ; Py y y F H j

More information

(Received 05 August 2013, accepted 15 July 2014)

(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 information

Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein

Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems Institute for Mathematics and Scientific Computing, Austria DISC Summerschool 5 Outline of the talk POD and singular value decomposition

More information

First Integrals/Invariants & Symmetries for Autonomous Difference Equations

First Integrals/Invariants & Symmetries for Autonomous Difference Equations Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 3, 1253 1260 First Integrals/Invariants & Symmetries for Autonomous Difference Equations Leon ARRIOLA Department of Mathematical

More information

FINAL EXAM, MATH 353 SUMMER I 2015

FINAL EXAM, MATH 353 SUMMER I 2015 FINAL EXAM, MATH 353 SUMMER I 25 9:am-2:pm, Thursday, June 25 I have neither given nor received any unauthorized help on this exam and I have conducted myself within the guidelines of the Duke Community

More information

Introduction of Partial Differential Equations and Boundary Value Problems

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

More information

2. Second-order Linear Ordinary Differential Equations

2. Second-order Linear Ordinary Differential Equations Advanced Engineering Mathematics 2. Second-order Linear ODEs 1 2. Second-order Linear Ordinary Differential Equations 2.1 Homogeneous linear ODEs 2.2 Homogeneous linear ODEs with constant coefficients

More information

Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations

Fibonacci 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 information

A Generalized Extended F -Expansion Method and Its Application in (2+1)-Dimensional Dispersive Long Wave Equation

A Generalized Extended F -Expansion Method and Its Application in (2+1)-Dimensional Dispersive Long Wave Equation Commun. Theor. Phys. (Beijing, China) 6 (006) pp. 580 586 c International Academic Publishers Vol. 6, No., October 15, 006 A Generalized Extended F -Expansion Method and Its Application in (+1)-Dimensional

More information

The Orchestra of Partial Differential Equations. Adam Larios

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

More information

Hamiltonian partial differential equations and Painlevé transcendents

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

More information

Conservation Laws of Surfactant Transport Equations

Conservation Laws of Surfactant Transport Equations Conservation Laws of Surfactant Transport Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Winter 2011 CMS Meeting Dec. 10, 2011 A. Cheviakov

More information

2. Examples of Integrable Equations

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

More information

Differential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension

Differential 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 information

Invariant Sets and Exact Solutions to Higher-Dimensional Wave Equations

Invariant 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 information

Integral Theorems. September 14, We begin by recalling the Fundamental Theorem of Calculus, that the integral is the inverse of the derivative,

Integral Theorems. September 14, We begin by recalling the Fundamental Theorem of Calculus, that the integral is the inverse of the derivative, Integral Theorems eptember 14, 215 1 Integral of the gradient We begin by recalling the Fundamental Theorem of Calculus, that the integral is the inverse of the derivative, F (b F (a f (x provided f (x

More information

Multidimensional partial differential equation systems: Nonlocal symmetries, nonlocal conservation laws, exact solutions

Multidimensional partial differential equation systems: Nonlocal symmetries, nonlocal conservation laws, exact solutions JOURNAL OF MATHEMATICAL PHYSICS 51, 103522 2010 Multidimensional partial differential equation systems: Nonlocal symmetries, nonlocal conservation laws, exact solutions Alexei F. Cheviakov 1,a and George

More information

arxiv: v1 [math.ap] 25 Aug 2009

arxiv: v1 [math.ap] 25 Aug 2009 Lie group analysis of Poisson s equation and optimal system of subalgebras for Lie algebra of 3 dimensional rigid motions arxiv:0908.3619v1 [math.ap] 25 Aug 2009 Abstract M. Nadjafikhah a, a School of

More information

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics 3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)

More information

Suggested Solution to Assignment 7

Suggested Solution to Assignment 7 MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the

More information

Noether Symmetries and Conservation Laws For Non-Critical Kohn-Laplace Equations on Three-Dimensional Heisenberg Group

Noether Symmetries and Conservation Laws For Non-Critical Kohn-Laplace Equations on Three-Dimensional Heisenberg Group arxiv:0706.1745v2 [math.ap] 13 Feb 2008 Noether Symmetries and Conservation Laws For Non-Critical Kohn-Laplace Equations on Three-Dimensional Heisenberg Group Igor Leite Freire Institute of Mathematics,

More information

Chapter 4. MAC Scheme. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes. u = 0

Chapter 4. MAC Scheme. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes. u = 0 Chapter 4. MAC Scheme 4.1. MAC Scheme and the Staggered Grid. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes equation in the velocity-pressure formulation: (4.1.1)

More information

MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014

MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014 MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014 Dr. E. Jacobs The main texts for this course are Calculus by James Stewart and Fundamentals of Differential Equations by Nagle, Saff

More information

Direct and inverse problems on conservation laws

Direct and inverse problems on conservation laws Direct and inverse problems on conservation laws Roman O. Popovych Wolfgang Pauli Institute, Vienna, Austria & Institute of Mathematics of NAS of Ukraine, Kyiv, Ukraine joint work with Alexander Bihlo,

More information

An improved collocation method based on deviation of the error for solving BBMB equation

An improved collocation method based on deviation of the error for solving BBMB equation Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 2, 2018, pp. 238-247 An improved collocation method based on deviation of the error for solving BBMB equation Reza

More information

Lump-type solutions to nonlinear differential equations derived from generalized bilinear equations

Lump-type solutions to nonlinear differential equations derived from generalized bilinear equations International Journal of Modern Physics B Vol. 30 2016) 1640018 16 pages) c World Scientific Publishing Company DOI: 10.1142/S021797921640018X Lump-type solutions to nonlinear differential equations derived

More information

Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations

Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute

More information

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012 MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to

More information

Proper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein

Proper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein Proper Orthogonal Decomposition for PDE Constrained Optimization Institute of Mathematics and Statistics, University of Constance Joined work with F. Diwoky, M. Hinze, D. Hömberg, M. Kahlbacher, E. Kammann,

More information

Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225

Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225 Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225 Dr. Asmaa Al Themairi Assistant Professor a a Department of Mathematical sciences, University of Princess Nourah bint Abdulrahman,

More information

arxiv: v2 [nlin.si] 3 Aug 2017

arxiv: 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 information

Problem 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 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 information

Math 46, Applied Math (Spring 2008): Final

Math 46, Applied Math (Spring 2008): Final Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the

More information

Symmetry Classification of KdV-Type Nonlinear Evolution Equations

Symmetry 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 information

Introduction and some preliminaries

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

More information

A symmetry-based method for constructing nonlocally related partial differential equation systems

A 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 information

Calculus of Variations and Computer Vision

Calculus of Variations and Computer Vision Calculus of Variations and Computer Vision Sharat Chandran Page 1 of 23 Computer Science & Engineering Department, Indian Institute of Technology, Bombay. http://www.cse.iitb.ernet.in/ sharat January 8,

More information

NUMERICAL METHODS FOR SOLVING NONLINEAR EVOLUTION EQUATIONS

NUMERICAL 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 information

Math 234 Final Exam (with answers) Spring 2017

Math 234 Final Exam (with answers) Spring 2017 Math 234 Final Exam (with answers) pring 217 1. onsider the points A = (1, 2, 3), B = (1, 2, 2), and = (2, 1, 4). (a) [6 points] Find the area of the triangle formed by A, B, and. olution: One way to solve

More information

An Introduction to Partial Differential Equations

An Introduction to Partial Differential Equations An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)

More information