Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps
|
|
- Andrew Craig
- 5 years ago
- Views:
Transcription
1 Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps Yuri A. Kuznetsov joint work with W. Govaerts, A. Dhooge(Gent), and E. Doedel (Montreal) LCBIF p.1/47
2 Contents References Limit cycles of ODEs and their local bifurcations Standard approach using Poincaré maps Continuation and normal form analysis using BVPs Open problems LCBIF p.2/47
3 1. References Kuznetsov, Yu. A., Govaerts, W., Doedel, E. J., and Dhooge, A. Numerical periodic normalization for codim 1 bifurcations of limit cycles. SIAM J. Numer. Analysis 43 (2005), Govaerts, W., Kuznetsov, Yu.A., and Dhooge, A. Numerical continuation of bifurcations of limit cycles in MATLAB. SIAM J. Sci. Comput. 27 (2005), Doedel, E.J., Govaerts, W., Kuznetsov, Yu.A., and Dhooge, A. Numerical continuation of branch points of equilibria and periodic orbits. Int. J. Bifurcation & Chaos 15 (2005), Doedel, E.J., Govaerts, W., and Kuznetsov, Yu.A. Computation of periodic solution bifurcations in ODEs using bordered systems. SIAM J. Numer. Analysis 41 (2003), LCBIF p.3/47
4 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Re µ LCBIF p.4/47
5 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Fold (LPC): µ 1 = 1; Re µ LCBIF p.4/47
6 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Fold (LPC): µ 1 = 1; Flip (PD): µ 1 = 1; Re µ LCBIF p.4/47
7 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Fold (LPC): µ 1 = 1; Flip (PD): µ 1 = 1; Torus (NS): µ 1,2 = e ±iθ 0. θ 0 Re µ LCBIF p.4/47
8 2.1. Generic LPC bifurcation: µ 1 = 1 C 1 C 2 C 0 W c β W c 0 W c β β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.5/47
9 2.2. Generic PD bifurcation: µ 1 = 1 W c β C 0 W0 c Wβ c C 1 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.6/47
10 2.3. Generic NS bifurcation: µ 1,2 = e ±iθ 0 C 0 T 2 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.7/47
11 3. Standard approach using Poincaré maps Poincaré map and its fixed points Continuation of PD and LP bifurcations of maps Continuation of NS bifurcation of maps Smooth normal forms of maps on center manifolds Difficulties using Poincaré maps Reference on finite-dimensional bordering techniques: Govaerts, W. Numerical Methods for Bifurcations of Dynamical Equilibria. SIAM, LCBIF p.8/47
12 3.1. Poincaré map and its fixed points Poincaré map P : R n 1 R m R n 1 A limit cycle C 0 corresponds to a fixed point: W c α y y 0 P(y, α) P(y 0, α) y 0 = 0. C 0 Σ Write P(y 0 + η, α 0 ) = y 0 + Aη B(η, η) C(η, η, η) + O( η 4 ), where A = P y (y 0, α 0 ) and, for i = 1, 2,..., n 1, B i (η, ζ) = C i (η, ζ, ω) = n 1 j,k=1 n 1 j,k,l=1 2 P i (y, α 0 ) y j y k 3 P i (y, α 0 ) y j y k y l η j ζ k, y=y0 η j ζ k ω l. y=y0 LCBIF p.9/47
13 3.2. Continuation of PD and LC bifurcations of maps Defining system: (y, α) R n 1 R 2 P(y, α) y = 0, g(y, α) = 0, where g is defined by solving P y(y, α) ± I n 1 w 1 v1 T 0 v g = 0 1 Vectors v 1, w 1 R n 1 are adapted to make the linear system nonsingular. (g y, g α ) can be computed efficiently using the adjoint linear system. LCBIF p.10/47
14 3.3. Continuation of NS bifurcation of maps Defining system: (y, α, κ) R n 1 R 2 R P(y, α) y = 0, g 11 (y, α, κ) = 0, g 22 (y, α, κ) = 0, where κ = cos θ 0 and g ij are defined by solving I n 1 2κP y (y, α) + [P y (y, α)] 2 w 1 w 2 r s v1 T 0 0 g 11 g 12 v2 T 0 0 g 21 g 22 = Vectors v 1,2, w 1,2 R n 1 are adapted to ensure unique solvability. Efficient computation of derivatives of g ij. LCBIF p.11/47
15 3.4. Smooth normal forms of maps on center manifolds Eigenvectors Normal form Critical coefficients LP Aq = q A T p = p p, q = 1 ξ β + ξ + bξ 2 +O(ξ 3 ), ξ R b = 1 2 p, B(q, q) PD Aq = q A T p = p p, q = 1 ξ (1 + β)ξ + cξ 3 +O(ξ 4 ), ξ R c = 1 6 p, C(q, q, q) + 3B(q, h 2) h 2 = (I n A) 1 B(q, q) NS Aq = e iθ 0 q A T p = e iθ 0 p e iνθ0 1 ν = 1, 2, 3, 4 ξ ξe iθ ( 1 + β + d ξ 2 ) + O( ξ 4 ), ξ C d = 1 2 e iθ 0 p, C(q, q, q) + 2B(q, h 11 ) + B( q, h 20 ) h 11 = (I n A) 1 B(q, q) p, q = 1 h 20 = (e 2iθ 0 I n A) 1 B(q, q) LCBIF p.12/47
16 LP-coefficient b Assume y 0 = 0 and write P(H) = AH B(H, H) + O( H 3 ), and locally represent the center manifold W0 c H : R R n 1, as the graph of a function u = H(ξ) = ξq h 2ξ 2 + O(ξ 3 ), ξ R, h 2 R n 1. The restriction of the Poincaré map to W0 c is ξ G(ξ) = ξ + bξ 2 + O(ξ 3 ). The invariance of the center manifold W0 c means P(H(ξ)) = H(G(ξ)). LCBIF p.13/47
17 A(ξq h 2ξ 2 ) B(ξq, ξq) + O( ξ 3 ) = (ξ + bξ 2 )q h 2ξ 2 + O( ξ 3 ) The ξ-terms give the identity: Aq = q. The ξ 2 -terms give the equation for h 2 : (A I n 1 )h 2 = B(q, q) + 2 bq. It is singular and its Fredholm solvability implies where A T p = p, p, q = 1. 1 b = p, B(q, q), 2 LCBIF p.14/47
18 3.5. Difficulties using Poincaré maps The computation of P can be unstable numerically near saddle or repelling limit cycles. Finite differences give low accuracy for B and C even if the cycle C 0 is stable. Simultaneous solving variational equations for the components of A, B, and C is possible but costly and difficult to implement. LCBIF p.15/47
19 4. Continuation and normal form analysis using BVPs Continuation of limit cycles in one parameter Simple bifurcation points Continuation of bifurcations in two parameters Periodic normalization on center manifolds LCBIF p.16/47
20 4.1. Continuation of limit cycles in one parameter Defining system (BVP with IC) [Keller & Doedel, 1981]: u(τ) T f(u(τ), α) = 0, τ [0, 1], u(0) u(1) = 0, 1 0 ũ(τ), u(τ) dτ = 0, where ũ is a reference periodic solution. Linearization w.r.t. (u, T, α): D T f x (u, α) f(u, α) T f α (u, α) δ 0 δ Int ũ 0 0 LCBIF p.17/47
21 Discretization via orthogonal collocation Mesh points: 0 = τ 0 < τ 1 < < τ N = 1. Basis points: τ i,j = τ i + j m (τ i+1 τ i ), i = 0, 1,..., N 1, j = 0, 1,..., m. Approximation: u (i) (τ) = m j=0 ui,j l i,j (τ), τ [τ i, τ i+1 ], where l i,j (τ) are the Lagrange basis polynomials and u i,m = u i+1,0. Collocation [de Boor & Swartz, 1973]: ( m j=0 ui,j l i,j i,k)) (ζ T f( m j=0 ui,j l i,j (ζ i,k ), α) = 0, u 0,0 u N 1,m = 0, m j=0 σ i,j ũ i,j, u i,j = 0, N 1 i=0 where ζ i,k, k = 1, 2,..., m, are the Gauss points. LCBIF p.18/47
22 Sparse Jacobian matrix u 0,0 u 0,1 u 1,0 u 1,1 u 2,0 u 2,1 u 3,0 T α LCBIF p.19/47
23 4.2. Simple bifurcation points Φ(τ) T f x (u(τ), α 0 )Φ(τ) = 0, Φ(0) = I n, Ψ(τ) + T fx T (u(τ), α 0 )Ψ(τ) = 0, Ψ(0) = I n. LPC: (Φ(1) I n )q 0 = 0, (Φ(1) I n )q 1 = q 0, (Ψ(1) I n )p 0 = 0, (Ψ(1) I n )p 1 = p 0. PD: (Φ(1) + I n )q 2 = 0, (Ψ(1) + I n )p 2 = 0. NS: κ = cos θ 0 (Φ(1) e iθ 0 I n )(q 3 + iq 4 ) = 0, (Ψ(1) e iθ 0 I n )(p 3 + ip 4 ) = 0. We have (I n 2κΦ(1) + Φ 2 (1))q 3,4 = 0. LCBIF p.20/47
24 4.3. Continuation of bifurcations in two parameters PD and LPC: (u, T, α) C 1 ([0, 1], R n ) R R 2 u(τ) T f(u(τ), α) = 0, τ [0, 1], u(0) u(1) = 0, 1 0 ũ(τ), u(τ) dτ = 0, G[u, T, α] = 0. NS: (u, T, α, κ) C 1 ([0, 1], R n ) R R 2 R u(τ) T f(u(τ), α) = 0, τ [0, 1], u(0) u(1) = 0, 1 0 ũ(τ), u(τ) dτ = 0, G 11 [u, T, α, κ] = 0, G 22 [u, T, α, κ] = 0. LCBIF p.21/47
25 PD-continuation There exist v 01, w 01 C 0 ([0, 1], R n ), and w 02 R n, such that N 1 : C 1 ([0, 1], R n ) R C 0 ([0, 1], R n ) R n R, D T f x (u, α) w 01 N 1 = δ 0 δ 1 w 02, Int v01 0 is one-to-one and onto near a simple PD bifurcation point. Define G by solving N 1 v 0 = 0 G. 1 LCBIF p.22/47
26 The BVP for (v, G) can be written in the classical form v(τ) T f x (u(τ), α)v(τ) + Gw 01 (τ) = 0, τ [0, 1], v(0) + v(1) + Gw 02 = 0, 1 0 v 01(τ), v(τ) dτ 1 = 0. One can take and w 02 = 0 w 01 (τ) = Ψ(τ)p 2, v 01 (τ) = Φ(τ)q 2. LCBIF p.23/47
27 LPC-continuation There exist v 01, w 01 C 0 ([0, 1], R n ), w 02 R n, and v 02, w 03 R such that N 2 : C 1 ([0, 1], R n ) R 2 C 0 ([0, 1], R n ) R n R 2, D T f x (u, α) f(u, α) w 01 δ 0 δ 1 0 w 02 N 2 =, Int f(u,α) 0 w 03 Int v01 v 02 0 is one-to-one and onto near a simple LPC bifurcation point. 0 v Define G by solving N 2 S = 0. 0 G 1 LCBIF p.24/47
28 NS-continuation There exist v 01, v 02, w 11, w 12 C 0 ([0, 2], R n ), and w 21, w 22 R n, such that N 3 : C 1 ([0, 2], R n ) R 2 C 0 ([0, 2], R n ) R n R 2, D T f x (u, α) w 11 w 12 δ 0 2κδ 1 + δ 2 w 21 w 22 N 3 =, Int v Int v is one-to-one and onto near a simple NS bifurcation point. 0 0 r s Define G jk by solving N 3 G 11 G 12 = G 21 G LCBIF p.25/47
29 Remarks on continuation of bifurcations After discretization via orthogonal collocation, all linear BVPs for G s have sparsity structure that is identical to that of the linearization of the BVP for limit cycles. For each defining system holds: Simplicity of the bifurcation + Transversality Regularity of the defining BVP. Jacobian matrix of each (discretized) defining BVP can be efficiently computed using adjoint linear BVP. Border adaptation using solutions of the adjoint linear BVPs. Actually implemented in MATCONT, also with compiled C-codes for the Jacobian matrices. LCBIF p.26/47
30 4.4. Periodic normalization on center manifolds Parameter-dependent periodic normal forms for LPC, PD, and NS Critical normal form coefficients References on periodic normal forms: Iooss, G. Global characterization of the normal form for a vector field near a closed orbit. J. Diff. Equations 76 (1988), Iooss, G. and Adelmeyer, M. Topics in Bifurcation Theory and Applications. World Scientific, LCBIF p.27/47
31 LPC bifurcation T 0 -periodic normal form on Wα: c dτ dt = 1 + ν(α) ξ + a(α)ξ 2 + O(ξ 3 ), dξ dt = β(α) + b(α)ξ 2 + O(ξ 3 ), where a, b R, sign(b) = sign( b). C 1 C 2 C 0 W c β W c 0 W c β β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.28/47
32 PD bifurcation 2T 0 -periodic normal form on Wα: c dτ dt = 1 + ν(α) + a(α)ξ 2 + O(ξ 4 ), dξ dt = β(α)ξ + c(α)ξ 3 + O(ξ 4 ), where a, c R, sign(c) = sign( c). W c β C 0 W0 c Wβ c C 1 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.29/47
33 NS bifurcation T 0 -periodic normal form on Wα: c dτ dt dξ dt = 1 + ν(α) + a(α) ξ 2 + O( ξ 4 ), ( = β(α) + iθ(α) T (α) where a R, d C, sign(re(d)) = sign(re( d)). ) ξ + d(α)ξ ξ 2 + O( ξ 4 ), C 0 T 2 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.30/47
34 Critical normal form coefficients At a codimension-one point write f(x 0 (t)+v, α 0 ) = f(x 0 (t), α 0 )+A(t)v+ 1 2 B(t; v, v)+1 6 C(t; v, v, v)+o( v 4 ), where A(t) = f x (x 0 (t), α 0 ) and the components of the multilinear functions B and C are given by and B i (t; u, v) = C i (t; u, v, w) = n j,k=1 n j,k,l=1 2 f i (x, α 0 ) x j x k 3 f i (x, α 0 ) x j x k x l for i = 1, 2,..., n. These are T 0 -periodic in t. u j v k x=x0 (t) u j v k w l, x=x0 (t) LCBIF p.31/47
35 Fold (LPC): µ 1 = 1 Critical center manifold W c 0 : τ [0, T 0], ξ R x = x 0 (τ) + ξv(τ) + H(τ, ξ), where H(T 0, ξ) = H(0, ξ), τ C 0 H(τ, ξ) = 1 2 h 2(τ)ξ 2 + O(ξ 3 ) Critical periodic normal form on W c 0 : dτ dt dξ dt = 1 ξ + aξ 2 + O(ξ 3 ), = bξ 2 + O(ξ 3 ), where a, b R, while the O(ξ 3 )-terms are T 0 -periodic in τ. ξ W c 0 LCBIF p.32/47
36 LPC: Eigenfunctions implying where ϕ satisfies v(τ) A(τ)v(τ) f(x 0 (τ), α 0 ) = 0, τ [0, T 0 ], v(0) v(t 0 ) = 0, T0 0 v(τ), f(x 0 (τ), α 0 ) dτ = 0, T0 0 ϕ (τ), f(x 0 (τ), α 0 ) dτ = 0, ϕ (τ) + A T (τ)ϕ (τ) = 0, τ [0, T 0 ], ϕ (0) ϕ (T 0 ) = 0, T0 0 ϕ (τ), v(τ) dτ 1 = 0. LCBIF p.33/47
37 LPC: Computation of b Substitute into Collect dx dt = x ξ dξ dt + x τ dτ dt ξ 0 : ẋ 0 = f(x 0, α 0 ), ξ 1 : v A(τ)v = ẋ 0, ξ 2 : ḣ 2 A(τ)h 2 = B(τ; v, v) 2af(x 0, α 0 ) + 2 v 2bv. Fredholm solvability condition b = 1 2 T0 0 ϕ (τ), B(τ; v(τ), v(τ)) + 2A(τ)v(τ) dτ. LCBIF p.34/47
38 LPC: Example in MATCONT (ABC-reaction model) u 1 = u 1 + p 1 (1 u 1 )e u 3, u 2 = u 2 + p 1 (1 u 1 p 5 u 2 )e u 3, u 3 = u 3 p 3 u 3 + p 1 p 4 (1 u 1 + p 2 p 5 u 2 )e u 3. p 3 = 1.5 p 4 = 8.0 p 5 = CPC 1.1 p b = 0 at CPC (cusp) p 1 LCBIF p.35/47
39 Flip (PD): µ 1 = 1 Critical center manifold W c 0 : τ [0, 2T 0], ξ R x = x 0 (τ) + ξw(τ) + H(τ, ξ), where H(2T 0, ξ) = H(0, ξ), τ C 0 H(τ, ξ) = 1 2 h 2(τ)ξ h 3(τ)ξ 3 + O(ξ 4 ) ξ Critical periodic normal form on W c 0 : dτ dt dξ dt = 1 + aξ 2 + O(ξ 4 ), = cξ 3 + O(ξ 4 ), where a, c R, while the O(ξ 4 )-terms are 2T 0 -periodic in τ. W c 0 LCBIF p.36/47
40 PD: Eigenfunctions w(τ) = v(τ), τ [0, T 0 ], v(τ T 0 ), τ [T 0, 2T 0 ], w (τ) = v (τ), τ [0, T 0 ], v (τ T 0 ), τ [T 0, 2T 0 ], with v(τ) A(τ)v(τ) = 0, τ [0, T 0 ], v(0) + v(t 0 ) = 0, T0 0 v(τ), v(τ) dτ 1 = 0, v (τ) + A T (τ)v (τ) = 0, τ [0, T 0 ], v (0) + v (T 0 ) = 0, T0 0 v (τ), v(τ) dτ 1/2 = 0. LCBIF p.37/47
41 PD: Quadratic terms ξ 2 : ḣ 2 A(τ)h 2 = B(τ; w, w) 2aẋ 0, τ [0, 2T 0 ]. Since Ker ( d dτ A(τ)) = span{w, ψ = ẋ 0 }, we must have where ψ satisfies 2T0 0 w (τ), B(τ; w(τ), w(τ)) 2aẋ 0 (τ) dτ = 0, 2T0 0 ψ (τ), B(τ; w(τ), w(τ)) 2aẋ 0 (τ) dτ = 0, ψ (τ) + A T (τ)ψ (τ) = 0, τ [0, T 0 ], ψ (0) ψ (T 0 ) = 0, T0 0 ψ (τ), f(x 0 (τ), α 0 ) dτ 1/2 = 0, and is extended to [T 0, 2T 0 ] by periodicity. LCBIF p.38/47
42 PD: Computation of a and h 2 a = 1 2 The first Fredholm condition holds identically for all a, while the second gives 2T0 0 ψ (τ), B(τ; w(τ), w(τ) dτ = T0 0 ψ (τ), B(τ; v(τ), v(τ) dτ. Define h 2 on [0, T 0 ] as the unique solution to ḣ 2 (τ) A(τ)h 2 (τ) B(τ; v(τ), v(τ)) + 2af(x 0 (τ), α 0 ) = 0, h 2 (0) h 2 (T 0 ) = 0, T0 0 ψ (τ), h 2 (τ) dτ = 0, and extend it by periodicity to [T 0, 2T 0 ]. LCBIF p.39/47
43 PD: Computation of c Cubic terms: ξ 3 : ḣ 3 A(τ)h 3 = C(τ; w, w, w) + 3B(τ; w, h 2 ) 6aẇ 6cw. The Fredholm solvability condition implies 6c = or 2T0 0 2T0 0 w (τ), C(τ; w(τ), w(τ), w(τ)) + 3B(τ; w(τ), h 2 (τ)) dτ w (τ), 6aA(τ)w(τ) dτ c = 1 3 T0 0 v (τ), C(τ; v(τ), v(τ), v(τ))+3b(τ; v(τ), h 2 (τ)) 6aA(τ)v(τ) dτ LCBIF p.40/47
44 Torus (NS): µ 1,2 = e ±iθ 0 with e iνθ0 1, ν = 1, 2, 3, 4 Critical center manifold W c 0 : τ [0, T 0], ξ C x = x 0 (τ) + ξv(τ) + ξ v(τ) + H(τ, ξ, ξ), H(T 0, ξ, ξ) = H(0, ξ, ξ), H(τ, ξ, ξ) = 1 2 h 20(τ)ξ 2 + h 11 (τ)ξ ξ h 02(τ) ξ h 30(τ)ξ h 21(τ)ξ 2 ξ h 12(τ)ξ ξ h 03(τ) ξ 3 + O( ξ 4 ). Critical periodic normal form on W0 c: dτ = 1 + a ξ 2 + O( ξ 4 ), dt dξ = iθ 0 ξ + dξ ξ 2 + O( ξ 4 ), dt T 0 where a R, d C, and the O( ξ 4 )-terms are T 0 -periodic in τ. LCBIF p.41/47
45 NS: Complex eigenfunctions and v(τ) A(τ)v(τ) + iθ 0 T 0 v(τ) = 0, τ [0, T 0 ], v(0) v(t 0 ) = 0, T0 0 v(τ), v(τ) dτ 1 = 0. v (τ) + A T (τ)v (τ) iθ 0 T 0 v (τ) = 0, τ [0, T 0 ], v (0) v (T 0 ) = 0, T0 0 v (τ), v(τ) dτ 1 = 0. LCBIF p.42/47
46 NS: Quadratic terms ξ 2 ξ0 : ḣ 20 A(τ)h iθ 0 T 0 h 20 = B(τ; v, v) Since e 2iθ 0 is not a multiplier of the critical cycle, the BVP ḣ 20 (τ) A(τ)h 20 (τ) + 2iθ 0 T 0 h 20 (τ) B(τ; v(τ), v(τ)) = 0, has a unique solution on [0, T 0 ]. ξ 2 : ḣ 11 A(τ)h 11 = B(τ; v, v) aẋ 0 Here Ker ( ) d dτ A(τ) h 20 (0) h 20 (T 0 ) = 0. = span(ϕ = ẋ 0 ). LCBIF p.43/47
47 NS: Computation of a and h 11 Define ϕ as the unique solution of ϕ (τ) + A T (τ)ϕ (τ) = 0, τ [0, T 0 ], ϕ (0) ϕ (T 0 ) = 0, T0 0 ϕ (τ), f(x 0 (τ), α 0 ) dτ 1 = 0. Fredholm solvability: a = T 0 0 ϕ (τ), B(τ; v(τ), v(τ) dτ. Then find h 11 on [0, T 0 ] from the BVP ḣ 11 (τ) A(τ)h 11 (τ) B(τ; v(τ), v(τ)) + af(x 0 (τ), α 0 ) = 0, h 11 (0) h 11 (T 0 ) = 0, T0 0 ϕ (τ), h 11 (τ) dτ = 0. LCBIF p.44/47
48 NS: Computation of d Cubic terms: ξ 2 ξ : ḣ21 Ah 21 + iθ 0 T 0 h 21 = 2B(τ; h 11, v) + B(τ; h 20, v) Fredholm solvability condition: + C(τ; v, v, v) 2a v 2dv. d = 1 2 T v (τ), C(τ; v(τ), v(τ), v(τ)) dτ T0 0 a v (τ), B(τ; h 11 (τ), v(τ)) + B(τ; h 20 (τ), v(τ)) dτ T0 0 v (τ), A(τ)v(τ) dτ + iaθ 0 T 0. LCBIF p.45/47
49 Remarks on numerical periodic normalization Only the derivatives of f(x, α 0 ) are used, not those of the Poincaré map P(y, α 0 ). Detection of codim 2 points is easy. After discretization via orthogonal collocation, all linear BVPs involved have the standard sparsity structure. One can re-use solutions to linear BVPs appearing in the continuation to compute the normal form coefficients. Actually implemented in MATCONT. LCBIF p.46/47
50 5. Open problems Automatic differentiation of the Poincaré map vs. BVPs. Periodic normal forms for codim 2 bifurcations of limit cycles. Branch switching at codim 2 points to limit cycle codim 1 continuation. LCBIF p.47/47
NBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011
NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets
More informationBIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs
BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange
More informationLecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)
Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic
More informationContinuation of cycle-to-cycle connections in 3D ODEs
HET p. 1/2 Continuation of cycle-to-cycle connections in 3D ODEs Yuri A. Kuznetsov joint work with E.J. Doedel, B.W. Kooi, and G.A.K. van Voorn HET p. 2/2 Contents Previous works Truncated BVP s with projection
More information( Bx + f(x, y) Cy + g(x, y)
Chapter 6 Center manifold reduction The previous chaper gave a rather detailed description of bifurcations of equilibria and fixed points in generic one-parameter families of ODEs and maps with minimal
More informationAn Introduction to Numerical Continuation Methods. with Application to some Problems from Physics. Eusebius Doedel. Cuzco, Peru, May 2013
An Introduction to Numerical Continuation Methods with Application to some Problems from Physics Eusebius Doedel Cuzco, Peru, May 2013 Persistence of Solutions Newton s method for solving a nonlinear equation
More informationNumerical Continuation of Bifurcations - An Introduction, Part I
Numerical Continuation of Bifurcations - An Introduction, Part I given at the London Dynamical Systems Group Graduate School 2005 Thomas Wagenknecht, Jan Sieber Bristol Centre for Applied Nonlinear Mathematics
More informationIMPROVED HOMOCLINIC PREDICTOR FOR BOGDANOV-TAKENS BIFURCATION
International Journal of Bifurcation and Chaos c World Scientific Publishing Company IMPROVED HOMOCLINIC PREDICTOR FOR BOGDANOV-TAKENS BIFURCATION YU.A. KUZNETSOV, H.G.E. MEIJER Department of Applied Mathematics,
More informationThe fold-flip bifurcation
The fold-flip bifurcation Yu.A. Kuznetsov, H.G.E. Meijer, and L. van Veen March 0, 003 Abstract The fold-flip bifurcation occurs if a map has a fixed point with multipliers + and simultaneously. In this
More informationAn Introduction to Numerical Continuation Methods. with Applications. Eusebius Doedel IIMAS-UNAM
An Introduction to Numerical Continuation Methods with Applications Eusebius Doedel IIMAS-UNAM July 28 - August 1, 2014 1 Persistence of Solutions Newton s method for solving a nonlinear equation [B83]
More informationcomputer session iv: Two-parameter bifurcation analysis of equilibria and limit cycles with matcont
computer session iv: Two-parameter bifurcation analysis of equilibria and limit cycles with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht May 16,
More informationStability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks
International Journal of Difference Equations ISSN 0973-5321, Volume 4, Number 1, pp. 155 163 (2009) http://campus.mst.edu/ijde Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationLINEAR SECOND-ORDER EQUATIONS
LINEAR SECON-ORER EQUATIONS Classification In two independent variables x and y, the general form is Au xx + 2Bu xy + Cu yy + u x + Eu y + Fu + G = 0. The coefficients are continuous functions of (x, y)
More informationConnecting Orbits with Bifurcating End Points
Connecting Orbits with Bifurcating End Points Thorsten Hüls Fakultät für Mathematik Universität Bielefeld Postfach 100131, 33501 Bielefeld Germany huels@mathematik.uni-bielefeld.de June 24, 2004 Abstract
More informationImproved Homoclinic Predictor for Bogdanov Takens Bifurcation
International Journal of Bifurcation and Chaos, Vol. 24, No. 4 (214) 14557 (12 pages) c World Scientific Publishing Company DOI: 1.1142/S2181274145576 Improved Homoclinic Predictor for Bogdanov Takens
More informationHomoclinic saddle to saddle-focus transitions in 4D systems
Faculty of Electrical Engineering, Mathematics & Computer Science Homoclinic saddle to saddle-focus transitions in 4D systems Manu Kalia M.Sc. Thesis July 2017 Assessment committee: Prof. Dr. S. A. van
More informationComputational Methods in Dynamical Systems and Advanced Examples
and Advanced Examples Obverse and reverse of the same coin [head and tails] Jorge Galán Vioque and Emilio Freire Macías Universidad de Sevilla July 2015 Outline Lecture 1. Simulation vs Continuation. How
More informationBifurcation Analysis of Large Equilibrium Systems in Matlab
Bifurcation Analysis of Large Equilibrium Systems in Matlab David S. Bindel 1, James W. Demmel 2, Mark J. Friedman 3, Willy J.F. Govaerts 4, and Yuri A. Kuznetsov 5 1 Department of Electrical Engineering
More informationElements of Applied Bifurcation Theory
Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Introduction to Dynamical Systems 1 1.1 Definition of a dynamical system 1 1.1.1 State space 1 1.1.2
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationDepartment of mathematics MA201 Mathematics III
Department of mathematics MA201 Mathematics III Academic Year 2015-2016 Model Solutions: Quiz-II (Set - B) 1. Obtain the bilinear transformation which maps the points z 0, 1, onto the points w i, 1, i
More informationtutorial ii: One-parameter bifurcation analysis of equilibria with matcont
tutorial ii: One-parameter bifurcation analysis of equilibria with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht February 13, 2018 1 This session
More informationStable cycles in a Cournot duopoly model of Kopel
Journal of Computational and Applied Mathematics 18 (008 47 58 www.elsevier.com/locate/cam Stable cycles in a Cournot duopoly model of Kopel W. Govaerts a,, R. Khoshsiar Ghaziani a,b a Department of Applied
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationElements of Applied Bifurcation Theory
Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Second Edition With 251 Illustrations Springer Preface to the Second Edition Preface to the First Edition vii ix 1 Introduction to Dynamical Systems
More informationTranscendental cases in stability problem. Hamiltonian systems
of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial
More informationBIFURCATIONS AND STRANGE ATTRACTORS IN A CLIMATE RELATED SYSTEM
dx dt DIFFERENTIAL EQUATIONS AND CONTROL PROCESSES N 1, 25 Electronic Journal, reg. N P23275 at 7.3.97 http://www.neva.ru/journal e-mail: diff@osipenko.stu.neva.ru Ordinary differential equations BIFURCATIONS
More informationReproducing formulas associated with symbols
Reproducing formulas associated with symbols Filippo De Mari Ernesto De Vito Università di Genova, Italy Modern Methods of Time-Frequency Analysis II Workshop on Applied Coorbit space theory September
More informationCONTINUATION OF CONNECTING ORBITS IN 3D-ODES: (II) CYCLE-TO-CYCLE CONNECTIONS arxiv: v1 [math.ds] 1 Apr 2008
CONTINUATION OF CONNECTING ORBITS IN 3D-ODES: (II) CYCLE-TO-CYCLE CONNECTIONS arxiv:0804.0179v1 [math.ds] 1 Apr 2008 E.J. DOEDEL 1, B.W. KOOI 2, YU.A. KUZNETSOV 3, and G.A.K. van VOORN 2 1 Department of
More informationHomoclinic Orbits of Planar Maps: Asymptotics and Mel nikov Functions
Department of Mathematics Mathematical Sciences Homoclinic Orbits of Planar Maps: Asymptotics and Mel nikov Functions A thesis submitted for the degree of Master of Science Author: Dirk van Kekem Project
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationDescription of the extension ddebiftool_nmfm
Description of the extension ddebiftool_nmfm Maikel Bosschaert, Bram Wage and Yuri Kuznetsov 16 September 2015 This addendum describes the normal form extension to DDE-BifTool, a bifurcation analysis toolbox
More informationHopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 65 71 Hopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation Faridon AMDJADI Department
More informationElliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 9, Number 2, pp. 227 237 (2014) http://campus.mst.edu/adsa Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationSummary of topics relevant for the final. p. 1
Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a
More informationUsing MatContM in the study of a nonlinear map in economics
Journal of Phsics: Conference Series PAPER OPEN ACCESS Using MatContM in the stud of a nonlinear map in economics To cite this article: Niels Neirnck et al 016 J. Phs.: Conf. Ser. 69 0101 Related content
More informationAvailable online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More informationGlobal bifurcation diagram for the Kerner-Konhäuser traffic flow model
International Journal of Bifurcation and Chaos c World Scientific Publishing Company Global bifurcation diagram for the Kerner-Konhäuser traffic flow model Joaquín Delgado Mathematics Department. UAM Iztapalapa.
More informationLectures on Dynamical Systems. Anatoly Neishtadt
Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances
More informationx R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ).
Continuation of equilibria Problem Parameter-dependent ODE ẋ = f(x, λ), x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x,
More informationOn the Stability of the Best Reply Map for Noncooperative Differential Games
On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University
More informationLECTURE NOTES IN PARTIAL DIFFERENTIAL EQUATIONS. Fourth Edition, February by Tadeusz STYŠ. University of Botswana
i LECTURE NOTES IN PARTIAL DIFFERENTIAL EQUATIONS Fourth Edition, February 2011 by Tadeusz STYŠ University of Botswana ii Contents 1 Solution of Partial Differential Equations 1 1.1 The General Solution
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationChanging coordinates to adapt to a map of constant rank
Introduction to Submanifolds Most manifolds of interest appear as submanifolds of others e.g. of R n. For instance S 2 is a submanifold of R 3. It can be obtained in two ways: 1 as the image of a map into
More informationBifurcation Analysis of Neuronal Bursters Models
Bifurcation Analysis of Neuronal Bursters Models B. Knowlton, W. McClure, N. Vu REU Final Presentation August 3, 2017 Overview Introduction Stability Bifurcations The Hindmarsh-Rose Model Our Contributions
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
More informationIntroduction to Bifurcation and Normal Form theories
Introduction to Bifurcation and Normal Form theories Romain Veltz / Olivier Faugeras October 9th 2013 ENS - Master MVA / Paris 6 - Master Maths-Bio (2013-2014) Outline 1 Invariant sets Limit cycles Stable/Unstable
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationDistance Between Ellipses in 2D
Distance Between Ellipses in 2D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationUNIVERSIDADE DE SÃO PAULO
UNIVERSIDADE DE SÃO PAULO Instituto de Ciências Matemáticas e de Computação ISSN 010-577 GLOBAL DYNAMICAL ASPECTS OF A GENERALIZED SPROTT E DIFFERENTIAL SYSTEM REGILENE OLIVEIRA CLAUDIA VALLS N o 41 NOTAS
More informationMath 489AB Exercises for Chapter 1 Fall Section 1.0
Math 489AB Exercises for Chapter 1 Fall 2008 Section 1.0 1.0.2 We want to maximize x T Ax subject to the condition x T x = 1. We use the method of Lagrange multipliers. Let f(x) = x T Ax and g(x) = x T
More informationSurvey of strong normal-internal k : l resonances in quasi-periodically driven oscillators for l = 1, 2, 3.
June, : WSPC - Proceedings Trim Size: in x in SPT-broer Survey of strong normal-internal k : l resonances in quasi-periodically driven oscillators for l =,,. H.W. BROER and R. VAN DIJK Institute for mathematics
More information(V.C) Complex eigenvalues and other topics
V.C Complex eigenvalues and other topics matrix Let s take a closer look at the characteristic polynomial for a 2 2 namely A = a c f A λ = detλi A = λ aλ d bc b d = λ 2 a + dλ + ad bc = λ 2 traλ + det
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationGlobal bifurcation diagram for the Kerner-Konhäuser traffic flow model
International Journal of Bifurcation and Chaos c World Scientific Publishing Company Global bifurcation diagram for the Kerner-Konhäuser traffic flow model Joaquín Delgado Mathematics Department. UAM Iztapalapa.
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationChapter 3. Local behavior of nonlinear systems. 3.1 Principle of linearized stability Linearized stability of fixed points of maps
Chapter 3 Local behavior of nonlinear systems In Chapter 2 we explained how, for linear dynamical systems, the qualitative as well as the quantitative behaviour can be described in terms of eigenvalues
More informationAnalytic invariants associated with a parabolic fixed point in C 2
Ergod. Th. & Dynam. Sys. (2008), 28, 1815 1848 doi:10.1017/s0143385707001046 c 2008 Cambridge University Press Printed in the United Kingdom Analytic invariants associated with a parabolic fixed point
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationCLASSIFICATIONS OF THE FLOWS OF LINEAR ODE
CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE PETER ROBICHEAUX Abstract. The goal of this paper is to examine characterizations of linear differential equations. We define the flow of an equation and examine
More informationHyperbolic PDEs. Chapter 6
Chapter 6 Hyperbolic PDEs In this chapter we will prove existence, uniqueness, and continuous dependence of solutions to hyperbolic PDEs in a variety of domains. To get a feel for what we might expect,
More information6 Linear Equation. 6.1 Equation with constant coefficients
6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ
More informationDynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationContinuation and bifurcation analysis of delay differential equations
Continuation and bifurcation analysis of delay differential equations Dirk Roose 1 and Robert Szalai 2 1 Department of Computer Science, Katholieke Universiteit Leuven, Belgium 2 Department of Engineering
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationfrom delay differential equations to ordinary differential equations (through partial differential equations)
from delay differential equations to ordinary differential equations (through partial differential equations) dynamical systems and applications @ BCAM - Bilbao Dimitri Breda Department of Mathematics
More informationAnalysis of the Takens-Bogdanov bifurcation on m parameterized vector fields
Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationCANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION
CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION WARREN WECKESSER Department of Mathematics Colgate University Hamilton, NY 3346 E-mail: wweckesser@mail.colgate.edu Cartwright and Littlewood discovered
More informationLecturer: Bengt E W Nilsson
2009 05 07 Lecturer: Bengt E W Nilsson From the previous lecture: Example 3 Figure 1. Some x µ s will have ND or DN boundary condition half integer mode expansions! Recall also: Half integer mode expansions
More informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
More informationLecture Notes on Numerical Analysis of Nonlinear Equations
Lecture Notes on Numerical Analysis of Nonlinear Equations Eusebius J Doedel Department of Computer Science, Concordia University, Montreal, Canada Numerical integrators can provide valuable insight into
More informationComparison of the Weyl integration formula and the equivariant localization formula for a maximal torus of Sp(1)
Comparison of the Weyl integration formula and the equivariant localization formula for a maximal torus of Sp() Jeffrey D. Carlson January 3, 4 Introduction In this note, we aim to prove the Weyl integration
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationCoherent structures near the boundary between excitable and oscillatory media
Coherent structures near the boundary between excitable and oscillatory media Jeremy Bellay University of Minnesota Department of Computer Science 200 Union St. S.E. Minneapolis, MN 55455, USA Arnd Scheel
More informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
More informationDefining Equations for Bifurcations and Singularities
Defining Equations for Bifurcations and Singularities John Guckenheimer and Yongjian Xiang Mathematics Department, Ithaca, NY 14853 For Vladimir Arnold on the occasion of his 65th birthday July 1, 2002
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationFirst order Partial Differential equations
First order Partial Differential equations 0.1 Introduction Definition 0.1.1 A Partial Deferential equation is called linear if the dependent variable and all its derivatives have degree one and not multiple
More informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
More informationAhmed Abbas Mizeal and Mudhir A. Abdul Hussain
ARCHIVUM MATHEMATICUM (BRNO) Tomus 8 (01) 7 37 TWO-MODE BIFURCATION IN SOLUTION OF A PERTURBED NONLINEAR FOURTH ORDER DIFFERENTIAL EQUATION Ahmed Abbas Mizeal and Mudhir A. Abdul Hussain Abstract. In this
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationMatcont Tutorial: ODE GUI version
Matcont Tutorial: ODE GUI version Hil Meijer Nottingham, Oct-Nov, 2011 If you want to get credit for solving a complicated mathematical problem, you will have to provide a full proof. But if you re trying
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationOPTIMAL CONTROL OF A CLASS OF PARABOLIC
OPTIMAL CONTROL OF A CLASS OF PARABOLIC TIME-VARYING PDES James Ng a, Ilyasse Aksikas a,b, Stevan Dubljevic a a Department of Chemical & Materials Engineering, University of Alberta, Canada b Department
More informationDirect and indirect methods for optimal control problems and applications in engineering
Direct and indirect methods for optimal control problems and applications in engineering Matthias Gerdts Computational Optimisation Group School of Mathematics The University of Birmingham gerdtsm@maths.bham.ac.uk
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 informationGlobal Bifurcation Diagram for the Kerner-Konhäuser Traffic Flow Model
arxiv:1311.4119v1 [math.ds] 17 Nov 2013 Global Bifurcation Diagram for the Kerner-Konhäuser Traffic Flow Model Joaquín Delgado Patricia Saavedra Keywords: Continuous traffic flow. Traveling waves. Bautin
More informationSparse representations from moments
Sparse representations from moments Bernard Mourrain Inria Méditerranée, Sophia Antipolis BernardMourrain@inriafr Sparse representation of signals Given a function or signal f (t): decompose it as f (t)
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationStability of flow past a confined cylinder
Stability of flow past a confined cylinder Andrew Cliffe and Simon Tavener University of Nottingham and Colorado State University Stability of flow past a confined cylinder p. 1/60 Flow past a cylinder
More information