An Introduction to Numerical Continuation Methods. with Application to some Problems from Physics. Eusebius Doedel. Cuzco, Peru, May 2013
|
|
- Rafe Garrett
- 5 years ago
- Views:
Transcription
1 An Introduction to Numerical Continuation Methods with Application to some Problems from Physics Eusebius Doedel Cuzco, Peru, May 2013
2 Persistence of Solutions Newton s method for solving a nonlinear equation G(u) = 0, G( ), u R n, may not converge if the initial guess is not close to a solution. However, one can put an artificial homotopy parameter in the equation. Actually, most equations already have parameters. We will discuss persistence of solutions to such equations. 1
3 The Implicit Function Theorem Let G : R n R R n satisfy (i) G(u 0, λ 0 ) = 0, u 0 R n, λ 0 R. (ii) G u (u 0, λ 0 ) is nonsingular (i.e., u 0 is an isolated solution), (iii) G and G u are smooth near u 0. Then there exists a unique, smooth solution family u(λ) such that G(u(λ), λ) = 0, for all λ near λ 0, u(λ 0 ) = u 0. NOTE: The IFT also holds in more general spaces 2
4 EXAMPLE : A Simple Homotopy. (Course demo Homotopy ) Let g(u, λ) = (u 2 1) (u 2 4) + λ u 2 e cu, where c is fixed, c = 0.1. When λ = 0 the equation has four solutions, namely, g(u, 0) = 0, u = ± 1, and u = ± 2. We have g u (u, λ) λ=0 d du (u, λ) λ=0 = 4u 3 10u. 3
5 Since we have g u (u, 0) = 4u 3 10u, g u ( 1, 0) = 6, g u ( 1, 0) = 6, g u ( 2, 0) = 12, g u ( 2, 0) = 12, which are all nonzero. Thus each of the four solutions when λ = 0 is isolated. Hence each of these solutions persists as λ becomes nonzero, ( at least for small values of λ ). 4
6 2 1 u λ Four solution families of g(u, λ) = 0. Note the fold. 5
7 NOTE: Each of the four solutions at λ = 0 is isolated. Thus each of these solutions persists as λ becomes nonzero. Only two of the four homotopies reach λ = 1. The two other homotopies meet at a fold. IFT condition (ii) is not satisfied at the fold. (Why not?) 6
8 Next we give a useful Corollary to the IFT : In our equation let G(u, λ) = 0, u, G(, ) R n, λ R, Then the equation can be written x ( ) u λ G(x) = 0, G : R n+1 R n.. DEFINITION. A solution x 0 of G(x) = 0 is regular if the matrix G 0 x G x (x 0 ), (with n rows and n + 1 columns) has maximal rank, i.e., if Rank(G 0 x) = n. 7
9 In the parameter formulation, we have Rank(G 0 x) = Rank(G 0 u G 0 λ) = n Above, and G(u, λ) = 0, (i) G 0 u is nonsingular, or (ii) N (G 0 u) denotes the null space of G 0 u, R(G 0 u) denotes the range of G 0 u, i.e., the linear space spanned by the n columns of G 0 u. dim N (G 0 u) = 1, and G 0 λ R(G0 u). 8
10 COROLLARY (to the IFT) : Let x 0 ( u 0, λ 0 ) be a regular solution of G(x) = 0. Then, near x 0, there exists a unique one-dimensional solution family x(s) with x(0) = x 0. PROOF. Since Rank( G 0 x ) = Rank( G 0 u G 0 λ ) = n, then either G 0 u is nonsingular and by the IFT we have u = u(λ) near x 0, or else we can interchange colums in the Jacobian G 0 x to see that the solution can locally be parametrized by one of the components of u. Thus a unique solution family passes through a regular solution. QED 9
11 NOTE: Such a solution family is often called a solution branch. Case (i) is where the IFT applies directly. Case (ii) is that of a simple fold. Thus even near a simple fold there is a unique solution family. However, near such a fold, the family can not be parametrized by λ. 10
12 EXAMPLES of IFT Application We give examples where the IFT shows that a given solution persists, (at least locally ), when a problem parameter is changed. We also consider cases where the conditions of the IFT are not satisfied. 11
13 EXAMPLE : The A B C Reaction. (Course demo ABC-reaction ) u 1 = u 1 + D(1 u 1 )e u 3, u 2 = u 2 + D(1 u 1 )e u 3 Dσu 2 e u 3, where u 3 = u 3 βu 3 + DB(1 u 1 )e u 3 + DBασu 2 e u 3, 1 u 1 is the concentration of A, u 2 is the concentration of B, u 3 is the temperature, α = 1, σ = 0.04, B = 8, D is the Damkohler number, β = 1.15 is the heat transfer coefficient. EXERCISE: Use the IFT to check that the zero stationary solution at D = 0 persists (locally). 12
14 7 6 5 u D Stationary solution family of the A B C reaction for β = Solid/dashed curves denote stable/unstable solutions. The red square denotes a Hopf bifurcation. 13
15 The Hopf Bifurcation Theorem THEOREM. Suppose that along a stationary solution family (u(λ), λ), of a complex conjugate pair of eigenvalues u = f(u, λ), α(λ) ± i β(λ), of f u (u(λ), λ) crosses the imaginary axis transversally, i.e., for some λ 0, α(λ 0 ) = 0, β(λ 0 ) 0, and α(λ 0 ) 0. Also assume that there are no other eigenvalues on the imaginary axis. Then there is a Hopf bifurcation, i.e., a family of periodic solutions bifurcates from the stationary solution at (u 0, λ 0 ). NOTE: The assumptions imply that f 0 u is nonsingular, so the stationary solution family can indeed be parametrized locally using λ. 14
16 7 6 5 u D Stationary (blue) and periodic (red) solution families of the A B C reaction. The periodic family contains stable solutions and terminates in a homoclinic orbit. 15
17 10 9 u u u 1 The periodic family orbit family (red) approaching a homoclinic orbit (black). The red dot is the Hopf point; the blue dot is the saddle point on the homoclinic. 16
18 EXAMPLE : A Predator-Prey Model. (Course demo Predator-prey ) u 1 = 3u 1 (1 u 1 ) u 1 u 2 λ(1 e 5u 1 ), u 2 = u 2 + 3u 1 u 2. Here u 1 may be thought of as fish and u 2 as sharks, while the term λ (1 e 5u 1 ), represents fishing, with fishing-quota λ. When λ = 0 the stationary solutions are 3u 1 (1 u 1 ) u 1 u 2 = 0 u 2 + 3u 1 u 2 = 0 (u 1, u 2 ) = (0, 0), (1, 0), ( 1 3, 2). 17
19 The Jacobian matrix is G u (u 1, u 2 ; λ) = ( 3 6u1 u 2 5λe 5u 1 u 1 3u u 1 ) G u (0, 0 ; 0) = G u (1, 0 ; 0) = ( ( ) ; real eigenvalues 3, -1 (unstable) ) ; real eigenvalues -3, 2 (unstable) G u ( 1 3, 2 ; 0) = ( ) ; complex eigenvalues 1 2 ± i (stable) All three Jacobians at λ = 0 are nonsingular. Thus, by the IFT, all three stationary points persist for (small) λ 0. 18
20 In this problem we can explicitly find all solutions: Branch I : (u 1, u 2 ) = (0, 0) Branch II : u 2 = 0, λ = 3u 1(1 u 1 ) 1 e 5u 1 ( Note that lim u 1 0 λ = lim 3(1 2u 1 ) u 1 0 5e 5u 1 = 3 5 ) Branch III : u 1 = 1 3, u 2 λ(1 e 5/3 ) = 0 u 2 = 2 3λ(1 e 5/3 ) These solution families intersect at two branch points, one of which is (u 1, u 2, λ) = (0, 0, 3/5). 19
21 fish sharks quota Stationary solution families of the predator-prey model. Solid/dashed curves denote stable/unstable solutions. Note the bifurcations (open squares), and the Hopf bifurcation (red square). 20
22 fish quota Stationary solution families of the predator-prey model, showing fish versus quota. Solid/dashed curves denote stable/unstable solutions. 21
23 Stability of branch I ( 3 5λ 0 G u (0, 0 ; λ) = 0 1 ) ; eigenvalues 3 5λ, 1. Hence the trivial solution is : and unstable if λ < 3/5, stable if λ > 3/5. Stability of branch II : This family has no stable positive solutions. 22
24 Stability of branch III : At λ H 0.67 the complex eigenvalues cross the imaginary axis: This crossing is a Hopf bifurcation, Beyond λ H there are stable periodic solutions. Their period T increases as λ increases. The period becomes infinite at λ = λ This final orbit is a heteroclinic cycle. 23
25 0.8 max fish quota Stationary (blue) and periodic (red) solution families of the predator-prey model. 24
26 sharks fish Some periodic solutions of the predator-prey model. The largest orbits are very close to a heteroclinic cycle. 25
27 From the Figure we can deduce the solution behavior for (slowly) increasing λ : Branch III is followed until λ H Periodic solutions of increasing period until λ = λ Collapse to trivial solution (Branch I). 26
28 EXAMPLE : The Gelfand-Bratu Problem. (Course demo Bratu ) u (x) + λ e u(x) = 0, x [0, 1], u(0) = u(1) = 0, defines the stationary states of a solid fuel ignition model. If λ = 0 then u(x) 0 is a solution. This problem can be thought of as an operator equation G(u; λ) = 0. We can use (a generalized) IFT to prove that there is a continuation u = u(λ), for λ small. 27
29 The linearization of G(u; λ) acting on v, i.e., G u (u; λ)v, corresponds to v (x) + λ e u(x) v = 0, v(0) = v(1) = 0, which for the solution u(x) 0 at λ = 0 becomes v (x) = 0, v(0) = v(1) = 0. Since this equation only has the zero solution v(x) 0, the IFT applies. Thus (locally) a unique solution family passes through u(x) 0, λ = 0. 28
30 u(x) dx λ Bifurcation diagram of the Gelfand-Bratu equation. 29
31 u(x) x Some solutions of the Gelfand-Bratu equation. 30
32 EXAMPLE : A Boundary Value Problem with Bifurcations. (Course demo BVP ) u + λ u(1 u) = 0, has u(x) 0 as a solution for all λ. u(0) = u(1) = 0, QUESTION: Are there more solutions? Again, this problem corresponds to an operator equation G(u; λ) = 0. Its linearization acting on v corresponds to v + λ (1 2u)v = 0, v(0) = v(1) = 0. 31
33 In particular, the linearization about the solution family u 0 is v + λ v = 0, v(0) = v(1) = 0, which for most values of λ only has the zero solution v(x) 0. However, when λ = λ k k 2 π 2, then there are nonzero solutions, namely, v(x) = sin(kπx), Thus the IFT does not apply at λ k = k 2 π 2. (We will see that these solutions are branch points.) 32
34 u at x = λ (scaled) Solution families to the nonlinear eigenvalue problem. 33
35 u x Some solutions to the nonlinear eigenvalue problem. 34
36 Parameter Continuation Here the continuation parameter is taken to be λ. Suppose we have a solution (u 0, λ 0 ) of G(u, λ) = 0, as well as the derivative u 0. Here u du dλ. We want to compute the solution u 1 at λ 1 λ 0 + λ. 35
37 "u" u 0 u (0) u 1 u (= du d λ at λ 0 ) λ λ λ 0 1 λ Graphical interpretation of parameter-continuation. 36
38 To solve the equation G(u 1, λ 1 ) = 0, for u 1 (with λ = λ 1 fixed) we use Newton s method G u (u (ν) 1, λ 1 ) u (ν) 1 = G(u (ν) 1, λ 1 ), u (ν+1) 1 = u (ν) 1 + u (ν) 1. ν = 0, 1, 2,. As initial approximation use u (0) 1 = u 0 + λ u 0. If G u (u 1, λ 1 ) is nonsingular, and λ sufficiently small, then the Newton convergence theory guarantees that this iteration will converge. 37
39 After convergence, the new derivative u 1 can be computed by solving G u (u 1, λ 1 ) u 1 = G λ (u 1, λ 1 ). This equation follows from differentiating with respect to λ at λ = λ 1. G(u(λ), λ) = 0, NOTE: u 1 can be computed without another LU-factorization of G u (u 1, λ 1 ). Thus the extra work to find u 1 is negligible. 38
40 EXAMPLE : The Gelfand-Bratu Problem. (Course demo Bratu ) u (x) + λ e u(x) = 0 for x [0, 1], u(0) = 0, u(1) = 0. If λ = 0 then u(x) 0 is an isolated solution. Discretize by introducing a mesh, 0 = x 0 < x 1 < < x N = 1, x j x j 1 = h, (1 j N), h = 1/N. The discrete equations are : u j+1 2u j + u j 1 h 2 + λ e u j = 0, j = 1,, N 1, with u 0 = u N = 0. 39
41 Let u u 1 u 2 u N 1. Then we can write the above as G( u, λ ) = 0, where G : R N 1 R R N 1. 40
42 Parameter-continuation : Suppose we have λ 0, u 0, and u 0. Set Newton s method : λ 1 = λ 0 + λ. G u (u (ν) 1, λ 1 ) u (ν) 1 = G(u (ν) 1, λ 1 ), u (ν+1) 1 = u (ν) 1 + u (ν) 1, ν = 0, 1, 2,, with After convergence find u 1 from u (0) 1 = u 0 + λ u 0. G u (u 1, λ 1 ) u 1 = G λ (u 1, λ 1 ). Repeat the above procedure to find u 2, u 3,. 41
43 Here G u ( u, λ ) = 2 h 2 + λe u 1 1 h λe u h 2 h h h 2 2 h 2 + λe u N 1. Thus we must solve a tridiagonal system for each Newton iteration. NOTE: The solution family has a fold where parameter-continuation fails! A better continuation method is pseudo-arclength continuation. There are also better discretizations, namely collocation, as used in AUTO. 42
44 Keller s Pseudo-Arclength Continuation This method allows continuation of a solution family past a fold. Suppose we have a solution (u 0, λ 0 ) of G( u, λ ) = 0, as well as the direction vector ( u 0, λ 0 ) of the solution branch. Pseudo-arclength continuation solves the following equations for (u 1, λ 1 ) : G(u 1, λ 1 ) = 0, (u 1 u 0 ) u 0 + (λ 1 λ 0 ) λ 0 s = 0. 43
45 "u" u 0 u 1 u 0 s 01 λ ( u, ) 0 λ 0 λ 0 λ 1 λ Graphical interpretation of pseudo-arclength continuation. 44
46 Solve the equations G(u 1, λ 1 ) = 0, (u 1 u 0 ) u 0 + (λ 1 λ 0 ) λ 0 s = 0. for (u 1, λ 1 ) by Newton s method: (G1 u) (ν) u 0 (G 1 λ )(ν) λ 0 ( (ν) ) u 1 λ (ν) 1 = G(u (ν) 1, λ (ν) 1 ) (u (ν) 1 u 0 ) u 0 + (λ (ν) 1 λ 0 ) λ 0 s. Compute the next direction vector by solving G1 u u 0 G 1 λ λ 0 ( u1 λ 1 ) =
47 NOTE: We can compute ( u 1, λ 1 ) with only one extra backsubstitution. The orientation of the branch is preserved if s is sufficiently small. Rescale the direction vector so that indeed u λ 2 1 = 1. 46
48 FACT: The Jacobian is nonsingular at a regular solution. ( ) u PROOF: Let x R λ n+1. Then pseudo-arclength continuation can be written more simply as G(x 1 ) = 0, (x 1 x 0 ) ẋ 0 s = 0, ( ẋ 0 = 1 ). x 1 x 0 s x 0 Pseudo-arclength continuation. 47
49 The pseudo-arclength equations are G(x 1 ) = 0, (x 1 x 0 ) ẋ 0 s = 0, ( ẋ 0 = 1 ). The matrix in Newton s method at s = 0 is ( ) G 0 x. ẋ 0 At a regular solution N (G 0 x) = Span{ẋ 0 }. We must show that ( ) G 0 x ẋ 0 is nonsingular at a regular solution. 48
50 If on the contrary ( ) G 0 x ẋ 0 is singular then G 0 x z = 0, ẋ 0 z = 0, for some vector z 0. Thus z = c ẋ 0, for some constant c. But then 0 = ẋ 0 z = c ẋ 0 ẋ0 = c ẋ 0 2 = c, so that z = 0, which is a contradiction. QED 49
51 EXAMPLE : The Gelfand-Bratu Problem. Use pseudo-arclength continuation for the discretized Gelfand-Bratu problem. Then the matrix ( ) Gx ẋ = ( ) Gu G λ u λ, in Newton s method is a bordered tridiagonal matrix :. 50
52 Following Folds At a fold the the behavior of a system can change drastically. How does the fold location change when a second parameter varies? Thus we want the compute a locus of folds in 2 parameters. We also want to compute loci of Hopf bifurcations in 2 parameters. (We do not discuss the theory here ) 51
53 Continuation of Folds Treat µ as one of the unknowns, and compute a solution family X(s) ( u(s), φ(s), λ(s), µ(s) ), to G(u, λ, µ) = 0, F(X) G u (u, λ, µ) φ = 0, φ φ 0 1 = 0, and the added continuation equation (u u 0 ) u 0 + (φ φ 0 ) φ0 + (λ λ 0 ) λ 0 + (µ µ 0 ) µ 0 s = 0. As before, ( u 0, φ0, λ0, µ 0 ), is the direction of the branch at the current solution point ( u 0, φ 0, λ 0, µ 0 ). 52
54 EXAMPLE : The A B C Reaction. (Course demo ABC-reaction ) The equations are u 1 = u 1 + D(1 u 1 )e u 3, u 2 = u 2 + D(1 u 1 )e u 3 Dσu 2 e u 3, u 3 = u 3 βu 3 + DB(1 u 1 )e u 3 + DBασu 2 e u 3, where 1 u 1 is the concentration of A, u 2 is the concentration of B, u 3 is the temperature, α = 1, σ = 0.04, B = 8, D is the Damkohler number, β is the heat transfer coefficient. 53
55 7 6 5 u D A stationary solution family for β =
56 β D The locus of folds for the A B C reaction. 55
57 Following Hopf bifurcations The extended system is f(u, λ, µ) = 0, F(u, φ, β, λ; µ) f u (u, λ, µ) φ i β φ = 0, φ φ 0 1 = 0, where F : R n C n R 2 R R n C n C, and to which we want to compute a solution family ( u, φ, β, λ, µ ), with u R n, φ C n, β, λ, µ R. Above φ 0 belongs to a reference solution ( u 0, φ 0, β 0, λ 0, µ 0 ), which typically is the latest computed solution point of a family. 56
58 β D The locus of Hopf bifurcations for A B C reaction. 57
59 7 6 5 u D Stationary solution families for β = 1.20, 1.25, 1.30,
60 Boundary Value Problems. Consider the first order system of ordinary differential equations u (t) f( u(t), µ, λ ) = 0, t [0, 1], where u( ), f( ) R n, λ R, µ R nµ, subject to boundary conditions b( u(0), u(1), µ, λ ) = 0, b( ) R n b, and integral constraints 1 0 q( u(s), µ, λ ) ds = 0, q( ) R nq. 59
61 This boundary value problem (BVP) is of the form where F( X ) = 0, X = ( u, µ, λ ), to which we add the continuation equation < X X 0, Ẋ 0 > s = 0, where X 0 represents the preceding solution on the branch. In detail, the continuation equation is 1 0 (u(t) u 0 (t)) u 0 (t) dt + (µ µ 0 ) µ 0 + (λ λ 0 ) λ 0 s = 0. 60
62 We want to solve BVP for u( ) and µ. We can think of λ as the continuation parameter. (In pseudo-arclength continuation, we don t distinguish µ and λ.) In order for problem to be formally well-posed we must have n µ = n b + n q n 0. A simple case is n q = 0, n b = n, for which n µ = 0. 61
63 Computing Periodic Solutions Periodic solutions can be computed very effectively using a BVP approach. This method also determines the period very accurately. Moreover, the technique can compute unstable periodic orbits. 62
64 Consider u (t) = f( u(t), λ ), u( ), f( ) R n, λ R. Fix the interval of periodicity by the transformation t t T. Then the equation becomes u (t) = T f( u(t), λ ), u( ), f( ) R n, T, λ R. and we seek solutions of period 1, i.e., u(0) = u(1). Note that the period T is one of the unknowns. 63
65 The above equations do not uniquely specify u and T : Assume that we have computed ( u k 1 ( ), T k 1, λ k 1 ), and we want to compute the next solution ( u k ( ), T k, λ k ). Specifically, u k (t) can be translated freely in time : If u k (t) is a periodic solution, then so is u k (t + σ), for any σ. Thus, a phase condition is needed. 64
66 An example is the Poincaré orthogonality condition (u k (0) u k 1 (0)) u k 1(0) = 0. (Below we derive a numerically more suitable phase condition.) u k-1 (0) u k-1 (0) u (0) k Graphical interpretation of the Poincaré phase condition. 65
67 An Integral Phase Condition If ũ k (t) is a solution then so is for any σ. ũ k (t + σ), We want the solution that minimizes D(σ) 1 0 ũ k (t + σ) u k 1 (t) 2 2 dt. The optimal solution ũ k (t + ˆσ), must satisfy the necessary condition D (ˆσ) = 0. 66
68 Differentiation gives the necessary condition Writing 1 0 ( ũ k (t + ˆσ) u k 1 (t) ) ũ k(t + ˆσ) dt = 0. u k (t) ũ k (t + ˆσ), gives 1 0 ( u k (t) u k 1 (t) ) u k(t) dt = 0. Integration by parts, using periodicity, gives 1 u 0 k (t) u k 1 (t) dt = 0. This is the integral phase condition. 67
69 Continuation of Periodic Solutions Pseudo-arclength continuation is also used to follow periodic solutions. It allows computation past folds along a family of periodic solutions. It also allows calculation of a vertical family of periodic solutions. For periodic solutions the continuation equation is 1 0 (u k (t) u k 1 (t)) u k 1 (t) dt + (T k T k 1 ) T k 1 + (λ k λ k 1 ) λ k 1 = s. 68
70 SUMMARY: We have the following equations for continuing periodic solutions : u k(t) = T f( u k (t), λ k ), u k (0) = u k (1), 1 0 u k (t) u k 1(t) dt = 0, withh continuation equation 1 0 (u k (t) u k 1 (t)) u k 1 (t) dt + (T k T k 1 ) T k 1 + (λ k λ k 1 ) λ k 1 = s, where u( ), f( ) R n, λ, T R. 69
71 9 8 7 max u D Stationary and periodic solution families of the A B C reaction: β =
72 Following folds along families of periodic orbits Recall that periodic orbits families can be computed using the equations u (t) T f( u(t), λ ) = 0, u(0) u(1) = 0, 1 0 u(t) u 0(t) dt = 0, where u 0 is a reference orbit, typically the latest computed orbit. The above boundary value problem is of the form where F( X, λ ) = 0, X = ( u, T ). 71
73 At a fold with respect to λ we have F X ( X, λ ) Φ = 0, where < Φ, Φ > = 1, X = ( u, T ), Φ = ( v, S ), or, in detail, v (t) T f u ( u(t), λ ) v S f( u(t), λ ) = 0, v(0) v(1) = 0, 1 0 v(t) u 0(t) dt = 0, 1 0 v(t) v(t) dt + S 2 = 1. 72
74 The complete extended system to follow a fold is F( X, λ, µ ) = 0, F X ( X, λ, µ ) Φ = 0, < Φ, Φ > 1 = 0, with two free problem parameters λ and µ. To the above we add the continuation equation < X X 0, Ẋ 0 > + < Φ Φ 0, Φ0 > + (λ λ 0 ) λ 0 + (µ µ 0 ) µ 0 s = 0. 73
75 In detail: u (t) T f( u(t), λ, µ ) = 0, 1 0 u(0) u(1) = 0, u(t) u 0(t) dt = 0, v (t) T f u ( u(t), λ, µ ) v S f( u(t), λ, µ ) = 0, with normalization 1 and continuation equation v(0) v(1) = 0, 1 (u(t) u 0 (t)) u 0 (t) dt + 0 v(t) u 0(t) dt = 0, v(t) v(t) dt + S 2 1 = 0, 1 0 (v(t) v 0 (t)) v 0 (t) dt + + (T 0 T ) T 0 + (S 0 S)Ṡ0 + (λ λ 0 ) λ 0 + (µ µ 0 ) µ 0 s = 0. 74
76 9 8 7 max u D Stationary and periodic solution families of the A B C reaction: β =
77 β D Loci of folds along periodic solution families for the A B C reaction. 76
78 9 8 7 max u D Stationary and periodic solution families of the A B C reaction: β =
79 9 8 7 max u D Stationary and periodic solution families of the A B C reaction: β =
80 9 8 7 max u D Stationary and periodic solution families of the A B C reaction: β =
81 9 8 7 max u D Stationary and periodic solution families of the A B C reaction: β =
82 9 8 7 u u u Periodic solutions along the isola for β = (Stable solutions are blue, unstable solutions are red.) 81
83 Discretization: Orthogonal Collocation Introduce a mesh { 0 = t 0 < t 1 < < t N = 1 }, where h j t j t j 1, (1 j N), Define the space of (vector) piecewise polynomials P m h as P m h { p h C[0, 1] : p h [tj 1,t j ] P m }, where P m is the space of (vector) polynomials of degree m. 82
84 The collocation method consists of finding p h P m h, µ R nµ, such that the following collocation equations are satisfied: p h(z j,i ) = f( p h (z j,i ), µ, λ ), j = 1,, N, i = 1,, m, and such that p h satisfies the boundary and integral conditions. The collocation points z j,i in each subinterval [ t j 1, t j ], are the (scaled) roots of the mth-degree orthogonal polynomial (Gauss points). 83
85 0 1 t t t t N t j-1 t j z j,1 z j,2 z j,3 t j-1 t j t j-2/3 t j-1/3 lj,3(t) lj,1(t) The mesh {0 = t 0 < t 1 < < t N = 1}, with collocation points and extended-mesh points shown for m = 3. Also shown are two of the four local Lagrange basis polynomials. 84
86 Since each local polynomial is determined by (m + 1) n, coefficients, the total number of degrees of freedom (considering λ as fixed) is (m + 1) n N + n µ. This is matched by the total number of equations : collocation : m n N, continuity : (N 1) n, constraints : n b + n q ( = n + n µ ). 85
87 Assume that the solution u(t) of the BVP is sufficiently smooth. Then the order of accuracy of the orthogonal collocation method is m, i.e., p h u = O(h m ). At the main meshpoints t j we have superconvergence : max j p h (t j ) u(t j ) = O(h 2m ). The scalar variables µ are also superconvergent. 86
88 Implementation For each subinterval [ t j 1, t j ], introduce the Lagrange basis polynomials defined by where { l j,i (t) }, j = 1,, N, i = 0, 1,, m, l j,i (t) = t j i m m k=0,k i t j i m t t j k m t j i m h j. t j k m The local polynomials can then be written m p j (t) = l j,i (t) u j i. m With the above choice of basis i=0 u j u(t j ) and u j i m where u(t) is the solution of the continuous problem. 87, u(t j i ), m
89 The collocation equations are p j(z j,i ) = f( p j (z j,i ), µ, λ ), i = 1,, m, j = 1,, N. The discrete boundary conditions are b i ( u 0, u N, µ, λ ) = 0, i = 1,, n b. The integral constraints can be discretized as N m j=1 i=0 ω j,i q k ( u j i m, µ, λ) = 0, k = 1,, n q, where the ω j,i are the Lagrange quadrature weights. 88
90 The continuation equation is 1 0 (u(t) u 0 (t)) u 0 (t) dt + (µ µ 0 ) µ 0 + (λ λ 0 ) λ 0 s = 0, where ( u 0, µ 0, λ 0 ), is the previous solution on the solution branch, and ( u 0, µ 0, λ0 ), is the normalized direction of the branch at the previous solution. The discretized continuation equation is N j=1 m i=0 ω j,i [ u j i m (u 0 ) j i m ] ( u 0 ) j i m + (µ µ 0 ) µ 0 + (λ λ 0 ) λ 0 s = 0. 89
91 Numerical Linear Algebra The complete discretization consists of m n N + n b + n q + 1, nonlinear equations, in the unknowns {u j i } R mnn+n, µ R nµ, λ R. m These equations can be solved by a Newton-Chord iteration. 90
92 We illustrate the numerical linear algebra for the case n = 2 ODEs, N = 4 mesh intervals, m = 3 collocation points, n b = 2 boundary conditions, n q = 1 integral constraint, and the continuation equation. The operations are also done on the right hand side, which is not shown. Entries marked have been eliminated by Gauss elimination. Entries marked denote fill-in due to pivoting. Most of the operations can be done in parallel. 91
93 u 0 u 1 3 u 2 3 u 1 u 2 u 3 u N µ λ The structure of the Jacobian 92
94 u 0 u 1 3 u 2 3 u 1 u 2 u 3 u N µ λ The system after condensation of parameters. 93
95 u 0 u 1 3 u 2 3 u 1 u 2 u 3 u N µ λ The preceding matrix, showing the decoupled sub-system. 94
96 u 0 u 1 3 u 2 3 u 1 u 2 u 3 u N µ λ Stage 1 of the nested dissection to solve the decoupled system. 95
97 u 0 u 1 3 u 2 3 u 1 u 2 u 3 u N µ λ Stage 2 of the nested dissection to solve the decoupled system. 96
98 u 0 u 1 3 u 2 3 u 1 u 2 u 3 u N µ λ The preceding matrix showing the final decoupled + sub-system. 97
99 u 0 u 1 3 u 2 3 u 1 u 2 u 3 u N µ λ A A B B + + A A B B The Floquet Multipliers are the eigenvalues of B 1 A. 98
100 Accuracy Test The Table shows the location of the fold in the Gelfand-Bratu problem, for 4 Gauss collocation points per mesh interval, and N mesh intervals u(x) dx N Fold location λ 99
101 Periodic Solutions of Conservative Systems EXAMPLE : A Model Conservative System. (Course demo Vertical-HB ) u 1 = u 2, u 2 = u 1 (1 u 1 ). PROBLEM: This equation has a family of periodic solutions, but no parameter! This system has a constant of motion, namely the Hamiltonian H(u 1, u 2 ) = 1 2 u u u
102 REMEDY: Introduce an unfolding term with unfolding parameter λ : u 1 = λ u 1 u 2, u 2 = u 1 (1 u 1 ). Then there is a vertical Hopf bifurcation from the trivial solution at λ =
103 norm λ Figure 1: Bifurcation diagram of the vertical Hopf bifurcation problem. 102
104 NOTE: The family of periodic solutions is vertical. The parameter λ is solved for in each continuation step. Upon solving, λ is found to be zero, up to numerical precision. One can use standard BVP continuation and bifurcation software. 103
105 u u 1 Figure 2: A phase plot of some periodic solutions. 104
106 EXAMPLE : The Circular Restricted 3-Body Problem (CR3BP). (Course demo CR3BP/Orbits ) x = 2y + x y = 2x + y z (1 µ) z = r1 3 (1 µ) (x + µ) r 3 1 (1 µ) y r 3 1 µ z r 3 2, µ y r 3 2 µ (x 1 + µ) r 3 2,, where r 1 = (x + µ) 2 + y 2 + z 2, r 2 = (x 1 + µ) 2 + y 2 + z 2. and ( x, y, z), denotes the position of the zero-mass body. For the Earth-Moon system µ
107 The CR3BP has one integral of motion, namely, the Jacobi-constant : J = x 2 + y 2 + z 2 2 U(x, y, z) µ 1 µ 2, where where U = 1 2 (x2 + y 2 ) + 1 µ r 1 + µ r 2, r 1 = (x + µ) 2 + y 2 + z 2, r 2 = (x 1 + µ) 2 + y 2 + z
108 BOUNDARY VALUE FORMULATION: x = T v x, y = T v y, z = T v z, v x = T [ 2v y + x (1 µ)(x + µ)r1 3 µ(x 1 + µ)r2 3 + λ v x ], v y = T [ 2v x + y (1 µ)yr1 3 µyr2 3 + λ v y ], v z = T [ (1 µ)zr1 3 µzr2 3 + λ v z ], with periodicity boundary conditions x(1) = x(0), y(1) = y(0), z(1) = z(0), v x (1) = v x (0), v y (1) = v y (0), v z (1) = v z (0), + phase constraint + pseudo-arclength equation. 107
109 NOTE: One can use standard BVP continuation and bifurcation software. The unfolding term λ v regularizes the continuation. λ will be zero, once solved for. Other unfolding terms are possible. 108
110 Families of Periodic Solutions of the Earth-Moon system. 109
111 The planar Lyapunov family L1. 110
112 The Halo family H1. 111
113 The Halo family H1. 112
114 The Vertical family V1. 113
115 The Axial family A1. 114
116 Stable and Unstable Manifolds EXAMPLE : Unstable Manifolds in the CR3BP. (Course demo CR3BP/Manifolds ) Small Halo orbits have one real Floquet multiplier outside the unit circle. Such Halo orbits are unstable. They have a 2D unstable manifold. 115
117 The unstable manifold can be computed by continuation. First compute a starting orbit in the manifold. Then continue the orbit keeping, for example, x(1) fixed. 116
118 Continuation, keeping the endpoint x(1) fixed. 117
119 The initial orbit can be taken to be much longer Continuation with x(1) fixed can lead to a Halo-to-torus connection! 118
120 119
121 The Halo-to-torus connection can be continued as a solution to F( X k ) = 0, where < X k X k 1, Ẋ k 1 > s = 0. X = ( Halo orbit, Floquet function, connecting orbit). 120
122 In detail, the continuation system is du dτ T uf(u(τ), µ, l) = 0, 1 0 u(1) u(0) = 0, u(τ), u 0 (τ) dτ = 0, dv dτ T ud u f(u(τ), µ, l)v(τ) + λ u v(τ) = 0, v(1) sv(0) = 0 (s = ±1), v(0), v(0) 1 = 0, dw dτ T wf(w(τ), µ, 0) = 0, w(0) (u(0) + εv(0)) = 0, w(1) x x Σ =
123 The system has 18 ODEs, 20 boundary conditions, 1 integral constraint. We need = 4 free parameters. Parameters: An orbit in the unstable manifold: T w, l, T u, x Σ Compute the unstable manifold: T w, l, T u, ε Follow a connecting orbit: λ u, l, T u, ε 122
124 First Example 123
125 124
126 125
127 126
128 127
129 128
130 129
131 130
132 131
133 132
134 133
135 134
136 135
137 136
138 Second Example 137
139 138
140 139
141 140
142 141
143 142
144 143
145 144
146 145
147 146
148 147
149 148
150 149
151 150
152 151
153 152
154 153
155 154
156 155
157 156
158 157
159 158
160 159
161 160
162 161
163 162
164 163
165 164
166 165
167 166
168 167
169 168
170 169
171 170
172 171
173 172
174 173
An 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 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 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 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 informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationAn introduction to numerical continuation with AUTO
An introduction to numerical continuation with AUTO Jennifer L. Creaser EPSRC Centre for Predictive Modelling in Healthcare University of Exeter j.creaser@exeter.ac.uk Internal 20 October 2017 AUTO Software
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 informationNumerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps
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
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 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 informationThe Computation of Periodic Solutions of the 3-Body Problem Using the Numerical Continuation Software AUTO
The Computation of Periodic Solutions of the 3-Body Problem Using the Numerical Continuation Software AUTO Donald Dichmann Eusebius Doedel Randy Paffenroth Astrodynamics Consultant Computer Science Applied
More informationNumerical techniques: Deterministic Dynamical Systems
Numerical techniques: Deterministic Dynamical Systems Henk Dijkstra Institute for Marine and Atmospheric research Utrecht, Department of Physics and Astronomy, Utrecht, The Netherlands Transition behavior
More information5.2.2 Planar Andronov-Hopf bifurcation
138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω
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 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 informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationNBA 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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
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 informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More information1 2 predators competing for 1 prey
1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2
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 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 informationComputing Periodic Orbits and their Bifurcations with Automatic Differentiation
Computing Periodic Orbits and their Bifurcations with Automatic Differentiation John Guckenheimer and Brian Meloon Mathematics Department, Ithaca, NY 14853 September 29, 1999 1 Introduction This paper
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationContinuous Threshold Policy Harvesting in Predator-Prey Models
Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July
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 informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
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 informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
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 informationIMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET
IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationIMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET
IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each
More informationMatrix & Linear Algebra
Matrix & Linear Algebra Jamie Monogan University of Georgia For more information: http://monogan.myweb.uga.edu/teaching/mm/ Jamie Monogan (UGA) Matrix & Linear Algebra 1 / 84 Vectors Vectors Vector: A
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 13
STAT 309: MATHEMATICAL COMPUTATIONS I FALL 208 LECTURE 3 need for pivoting we saw that under proper circumstances, we can write A LU where 0 0 0 u u 2 u n l 2 0 0 0 u 22 u 2n L l 3 l 32, U 0 0 0 l n l
More informationSolution to Homework 1
Solution to Homework Sec 2 (a) Yes It is condition (VS 3) (b) No If x, y are both zero vectors Then by condition (VS 3) x = x + y = y (c) No Let e be the zero vector We have e = 2e (d) No It will be false
More informationMath 4A Notes. Written by Victoria Kala Last updated June 11, 2017
Math 4A Notes Written by Victoria Kala vtkala@math.ucsb.edu Last updated June 11, 2017 Systems of Linear Equations A linear equation is an equation that can be written in the form a 1 x 1 + a 2 x 2 +...
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationAlgebra II. Paulius Drungilas and Jonas Jankauskas
Algebra II Paulius Drungilas and Jonas Jankauskas Contents 1. Quadratic forms 3 What is quadratic form? 3 Change of variables. 3 Equivalence of quadratic forms. 4 Canonical form. 4 Normal form. 7 Positive
More informationMath Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT
Math Camp II Basic Linear Algebra Yiqing Xu MIT Aug 26, 2014 1 Solving Systems of Linear Equations 2 Vectors and Vector Spaces 3 Matrices 4 Least Squares Systems of Linear Equations Definition A linear
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 informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
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 informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationNumerical Methods for Two Point Boundary Value Problems
Numerical Methods for Two Point Boundary Value Problems Graeme Fairweather and Ian Gladwell 1 Finite Difference Methods 1.1 Introduction Consider the second order linear two point boundary value problem
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 informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
More informationInvariant Manifolds of Dynamical Systems and an application to Space Exploration
Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated
More informationNonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University
Nonlinear stability of time-periodic viscous shocks Margaret Beck Brown University Motivation Time-periodic patterns in reaction-diffusion systems: t x Experiment: chemical reaction chlorite-iodite-malonic-acid
More informationTangent spaces, normals and extrema
Chapter 3 Tangent spaces, normals and extrema If S is a surface in 3-space, with a point a S where S looks smooth, i.e., without any fold or cusp or self-crossing, we can intuitively define the tangent
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationIntroduction to Algebraic and Geometric Topology Week 14
Introduction to Algebraic and Geometric Topology Week 14 Domingo Toledo University of Utah Fall 2016 Computations in coordinates I Recall smooth surface S = {f (x, y, z) =0} R 3, I rf 6= 0 on S, I Chart
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 1: Course Overview & Matrix-Vector Multiplication Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 20 Outline 1 Course
More informationStudy Guide for Linear Algebra Exam 2
Study Guide for Linear Algebra Exam 2 Term Vector Space Definition A Vector Space is a nonempty set V of objects, on which are defined two operations, called addition and multiplication by scalars (real
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More information(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,
2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a
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 information27. Topological classification of complex linear foliations
27. Topological classification of complex linear foliations 545 H. Find the expression of the corresponding element [Γ ε ] H 1 (L ε, Z) through [Γ 1 ε], [Γ 2 ε], [δ ε ]. Problem 26.24. Prove that for any
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 informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationMAT 610: Numerical Linear Algebra. James V. Lambers
MAT 610: Numerical Linear Algebra James V Lambers January 16, 2017 2 Contents 1 Matrix Multiplication Problems 7 11 Introduction 7 111 Systems of Linear Equations 7 112 The Eigenvalue Problem 8 12 Basic
More information10.34 Numerical Methods Applied to Chemical Engineering Fall Quiz #1 Review
10.34 Numerical Methods Applied to Chemical Engineering Fall 2015 Quiz #1 Review Study guide based on notes developed by J.A. Paulson, modified by K. Severson Linear Algebra We ve covered three major topics
More informationLiouville theorems for superlinear parabolic problems
Liouville theorems for superlinear parabolic problems Pavol Quittner Comenius University, Bratislava Workshop in Nonlinear PDEs Brussels, September 7-11, 2015 Tintin & Prof. Calculus (Tournesol) c Hergé
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationChapter 2 Hopf Bifurcation and Normal Form Computation
Chapter 2 Hopf Bifurcation and Normal Form Computation In this chapter, we discuss the computation of normal forms. First we present a general approach which combines center manifold theory with computation
More informationDYNAMICAL SYSTEMS PROBLEMS. asgor/ (1) Which of the following maps are topologically transitive (minimal,
DYNAMICAL SYSTEMS PROBLEMS http://www.math.uci.edu/ asgor/ (1) Which of the following maps are topologically transitive (minimal, topologically mixing)? identity map on a circle; irrational rotation of
More informationMath 321: Linear Algebra
Math 32: Linear Algebra T. Kapitula Department of Mathematics and Statistics University of New Mexico September 8, 24 Textbook: Linear Algebra,by J. Hefferon E-mail: kapitula@math.unm.edu Prof. Kapitula,
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More information(a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =
. (5 points) (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? dim N(A), since rank(a) 3. (b) If we also know that Ax = has no solution, what do we know about the rank of A? C(A)
More informationMath 471 (Numerical methods) Chapter 3 (second half). System of equations
Math 47 (Numerical methods) Chapter 3 (second half). System of equations Overlap 3.5 3.8 of Bradie 3.5 LU factorization w/o pivoting. Motivation: ( ) A I Gaussian Elimination (U L ) where U is upper triangular
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
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 informationNOTES ON LINEAR ODES
NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary
More informationComplex Behavior in Coupled Nonlinear Waveguides. Roy Goodman, New Jersey Institute of Technology
Complex Behavior in Coupled Nonlinear Waveguides Roy Goodman, New Jersey Institute of Technology Nonlinear Schrödinger/Gross-Pitaevskii Equation i t = r + V (r) ± Two contexts for today: Propagation of
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 information2.4 Eigenvalue problems
2.4 Eigenvalue problems Associated with the boundary problem (2.1) (Poisson eq.), we call λ an eigenvalue if Lu = λu (2.36) for a nonzero function u C 2 0 ((0, 1)). Recall Lu = u. Then u is called an eigenfunction.
More informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More informationSolution of Linear Equations
Solution of Linear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 7, 07 We have discussed general methods for solving arbitrary equations, and looked at the special class of polynomial equations A subclass
More informationSpike-adding canard explosion of bursting oscillations
Spike-adding canard explosion of bursting oscillations Paul Carter Mathematical Institute Leiden University Abstract This paper examines a spike-adding bifurcation phenomenon whereby small amplitude canard
More informationThe Higgins-Selkov oscillator
The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More informationSolutions Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du
Solutions hapter 9 dode 9 asic Solution Techniques 9 hoose one or more of the following differential equations, and then: (a) Solve the equation directly (b) Write down its phase plane equivalent, and
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More information