Application of Euler and Homotopy Operators to Integrability Testing
|
|
- Hilary Fleming
- 5 years ago
- Views:
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 Willy Hereman and Douglas Poole Department of Mathematical and Computer Sciences Colorado
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations
Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado whereman@mines.edu
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations
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 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 informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions
. Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multi-dimensions Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden,
More informationSymbolic 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 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 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 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 informationThe first three (of infinitely many) conservation laws for (1) are (3) (4) D t (u) =D x (3u 2 + u 2x ); D t (u 2 )=D x (4u 3 u 2 x +2uu 2x ); D t (u 3
Invariants and Symmetries for Partial Differential Equations and Lattices Λ Ünal Göktaοs y Willy Hereman y Abstract Methods for the computation of invariants and symmetries of nonlinear evolution, wave,
More informationSymbolic Computation of Conservation Laws, 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 informationarxiv:nlin/ v1 [nlin.si] 16 Nov 2003
Centre de Recherches Mathématiques CRM Proceedings and Lecture Notes arxiv:nlin/0311029v1 [nlinsi] 16 Nov 2003 Symbolic Computation of Conserved Densities, Generalized Symmetries, and Recursion Operators
More 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 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 informationExact 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 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 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 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 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 information. 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 informationConstruction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem
Proceedings of 4th Workshop Group Analysis of Differential Equations & Integrability 2009, 1 23 Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem George W. BLUMAN,
More informationSummer 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 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 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 informationDirect construction method for conservation laws of partial differential equations Part I: Examples of conservation law classifications
Euro. Jnl of Applied Mathematics (2002), vol. 13, pp. 545 566. c 2002 Cambridge University Press DOI: 10.1017/S095679250100465X Printed in the United Kingdom 545 Direct construction method for conservation
More informationSymbolic Computation of. Conserved Densities and Fluxes for. Nonlinear Systems of Differential-Difference. Equations
Symbolic Computation of Conserved Densities and Fluxes for Nonlinear Systems of Differential-Difference Equations by Holly Eklund Copyright c 003 Holly Eklund, All rights reserved. ii A thesis submitted
More informationFinal 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 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 informationInductive 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 informationPARTIAL 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 informationProlongation structure for nonlinear integrable couplings of a KdV soliton hierarchy
Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy Yu Fa-Jun School of Mathematics and Systematic Sciences, Shenyang Normal University, Shenyang 110034, China Received
More 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 informationGradient, 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 information3 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 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 informationHomotopy Perturbation Method for Solving Systems of Nonlinear Coupled Equations
Applied Mathematical Sciences, Vol 6, 2012, no 96, 4787-4800 Homotopy Perturbation Method for Solving Systems of Nonlinear Coupled Equations A A Hemeda Department of Mathematics, Faculty of Science Tanta
More informationSymbolic 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 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 informationFinal: 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 informationSymmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota
Symmetry Methods for Differential Equations and Conservation Laws Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Santiago, November, 2010 Symmetry Groups of Differential Equations
More informationMcGill 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 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 informationDifferential 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 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 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 informationChapter 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 informationSecond 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 informationDerivatives 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 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 informationOn 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 informationGroup Actions and Cohomology in the Calculus of Variations
Group Actions and Cohomology in the Calculus of Variations JUHA POHJANPELTO Oregon State and Aalto Universities Focused Research Workshop on Exterior Differential Systems and Lie Theory Fields Institute,
More informationLeplace 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 informationLax 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 informationExtra 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 informationMultiple 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 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 informationAPPLIED 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 informationLecture 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 informationMoving 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 informationTravelling 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 informationOVER 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)
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 informationProper 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 informationFirst 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 informationFINAL 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 informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More information2. 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 informationFibonacci tan-sec method for construction solitary wave solution to differential-difference equations
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 7 (2011) No. 1, pp. 52-57 Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations
More informationA 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 informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More 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 informationConservation 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 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 informationDifferential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension
Symposium in honor of George Bluman UBC May 12 16, 2014 Differential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension Chaolu Shanghai Maritime University
More 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 informationIntegral 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 informationMultidimensional 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 informationarxiv: 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 informationChain 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 informationSuggested 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 informationNoether 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 informationChapter 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 informationMA 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 informationDirect 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 informationAn 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 informationLump-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 informationSimple 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 informationMATH 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 informationProper 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 informationLecture 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 informationarxiv: v2 [nlin.si] 3 Aug 2017
A generalised multicomponent system of Camassa-Holm-Novikov equations arxiv:1608.04604v2 [nlin.si] 3 Aug 2017 1. Introduction Diego Catalano Ferraioli 1 and Igor Leite Freire 2 Departamento de Matemática
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationMath 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 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 informationIntroduction 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 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 informationCalculus 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 informationNUMERICAL METHODS FOR SOLVING NONLINEAR EVOLUTION EQUATIONS
NUMERICAL METHODS FOR SOLVING NONLINEAR EVOLUTION EQUATIONS Thiab R. Taha Computer Science Department University of Georgia Athens, GA 30602, USA USA email:thiab@cs.uga.edu Italy September 21, 2007 1 Abstract
More informationMath 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 informationAn 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