Dynamics of Adaptive Time-Stepping ODE solvers
|
|
- Janis Douglas
- 6 years ago
- Views:
Transcription
1 Dynamics of Adaptive Time-Stepping ODE solvers A.R. Humphries Joint work with: NICK CHRISTODOULOU PAUL FENTON RATNAM VIGNESWARAN ARH acknowledges support of the Engineering and Physical Sciences Research Council (UK), Leverhulme Trust (UK), and McGill. Kansas Dec 2002 p.1/32
2 Abstract For efficiency, variable time-stepping methods are often used to numerically integrate dynamical systems. The flow on chaotic attractors is often organised by the unstable manifolds of the fixed points, and it is thus necessary to obtain good numerical approximations in the neighbourhood of fixed points to reproduce the dynamics. However the standard adaptive algorithm typically fails to do this. Implicit methods designed for stiff problems are also unsuitable; they typically destroy the structure of the unstable manifold unless very small step-sizes are used. We will present examples to illustrate these poor dynamical behaviours, together with theoretical results on the approximation of stable/unstable manifolds, and suggest a phase space/stability based improvement to the standard algorithm. Kansas Dec 2002 p.2/32
3 1. Approximation of a Dynamical System with a Fixed Step-Size Runge-Kutta Method Contents 2. Approximation of Local Unstable Manifolds with a Fixed Step-Size Runge-Kutta Method 3. Approximation of a Dynamical System with a Variable Step-Size Runge-Kutta Embedded Pair 4. Approximation of Local Unstable Manifolds with a Variable Step-Size Runge-Kutta Pair 5. Phase Space Stability Error Control Kansas Dec 2002 p.3/32
4 Approximation of a dynamical system with a fixed step-size Runge-Kutta method The dynamical system u(t) = f(u(t)), u(0) = U IR d, has solution operator S( ) and so u(t) = S(t)U for all t. Kansas Dec 2002 p.4/32
5 Approximation of a dynamical system with a fixed step-size Runge-Kutta method The dynamical system u(t) = f(u(t)), u(0) = U IR d, has solution operator S( ) and so u(t) = S(t)U for all t. c A b T The order p Runge-Kutta method has evolution map S h. Example Forward Euler is defined by u n+1 = u n + hf(u n ) S h (u n ), n 0, u 0 = U. Kansas Dec 2002 p.4/32
6 Approximation of a dynamical system with a fixed step-size Runge-Kutta method The dynamical system u(t) = f(u(t)), u(0) = U IR d, has solution operator S( ) and so u(t) = S(t)U for all t. c A b T The order p Runge-Kutta method has evolution map S h. S h (u) advances the numerical solution with step-size h u n+1 = S h (u n ), n 0, u 0 = U. Each u n is an approximation of S(nh)U = u(nh). Kansas Dec 2002 p.4/32
7 Organisation of the flow 50 Phase Space Dormand Prince 5(4) applied to Lorenz Flow in forward time organised by fixed points and unstable manifolds. u So consider approximation of unstable manifolds by numerical methods u 2 u 1 Kansas Dec 2002 p.5/32
8 Organisation of the flow 50 Phase Space Dormand Prince 5(4) applied to Lorenz Flow in forward time organised by fixed points and unstable manifolds. u So consider approximation of unstable manifolds by numerical methods u u 2 Lorenz vector field is f(x, y, z) = (σ(y x), rx y xz, xy bz). Relationship of flow to fixed points obvious from figure. Kansas Dec 2002 p.5/32
9 Approximation of Local Unstable Manifolds with a Fixed Step-Size Runge-Kutta Method Stable Manifold W s (u^) The unstable manifold of equilibrium point û is W s,δ (u^) δ W u (u^) W u (û) = {u IR d : S( t)u û 0 as t }. W u,δ (u^) Unstable Manifold Let δ > 0. The local unstable manifold of û is W u,δ (û) = {u W u (û) : S( t)u û δ t 0}. Kansas Dec 2002 p.6/32
10 Approximation of Local Unstable Manifolds with a Fixed Step-Size Runge-Kutta Method Stable Manifold W s (u^) The unstable manifold of equilibrium point û is W s,δ (u^) δ W u (u^) W u (û) = {u IR d : S( t)u û 0 as t }. W u,δ (u^) Unstable Manifold Let δ > 0. The local unstable manifold of û is W u,δ (û) = {u W u (û) : S( t)u û δ t 0}. Generate {u n } n 0 using a fixed step-size h. unstable h-manifold of û is The Stable h Manifold W s h (u^) W u h (û) = {u IRd u 0 = u, u k û 0 as k }. W s,δ h (u^) δ W u h (u^) Let δ > 0. The local unstable h-manifold of û is W u,δ h (u^) Unstable h Manifold W u,δ h (û) = {u W h u (û) u 0 = u, u k û δ k}. Kansas Dec 2002 p.6/32
11 Local Unstable Manifold Theorem Theorem Apply a fixed step-size Runge-Kutta method of order p to u = f(u). Let û be a hyperbolic equilibrium and f C p+1 (IR d ). Then there exists C, H, > 0 such that δ (0, ), h (0, H) the following holds: for each u W u,δ (û), there exists a u h W u,δ h (û) such that u u h Ch p u û 2 ; and for each u h W u,δ h (û), there exists u W u,δ (û) such that u u h Ch p u h û 2. Kansas Dec 2002 p.7/32
12 Proof Let Df(û) be the Jacobian of f evaluated at û. Shift the coordinates v = u û. Linearise the solution operator and the Runge-Kutta evolution map about 0 S(h)v = exp(hdf(û))v+g h (v), Ŝ h (v) = R(hDf(û))v+N h (v), where R is matrix generalisation of linear stability function. Show that both W u,δ (û) and W u,δ h (û) are indeed manifolds. Show that both W u,δ (û) and W u,δ h (û) are representable as graphs. Show that the graphs are close. Kansas Dec 2002 p.8/32
13 Remarks A parallel result holds for local stable manifolds. Kansas Dec 2002 p.9/32
14 A parallel result holds for local stable manifolds. Remarks The result is a generalisation of Beyn s result: there exists C > 0 such that u u h Ch p. Kansas Dec 2002 p.9/32
15 A parallel result holds for local stable manifolds. Remarks The result is a generalisation of Beyn s result: there exists C > 0 such that u u h Ch p. Theorem 1 implies that W u,δ (û) and W u,δ h (û) are tangential at the fixed point, and so is a form of local (un)stable manifold theorem. This follows from û 2 on RHS of equations; so distance between numerical and exact manifolds depends on the square of the distance from the fixed point. Kansas Dec 2002 p.9/32
16 Approximation of a Dynamical System with a Variable Step-Size Runge-Kutta Embedded Pair Consider embedded Runge-Kutta pair with p p = 1. c A u n+1 = S hn (u n ) b T order p ũ n+1 = S hn (u n ) bt order p Kansas Dec 2002 p.10/32
17 Approximation of a Dynamical System with a Variable Step-Size Runge-Kutta Embedded Pair Consider embedded Runge-Kutta pair with p p = 1. c A u n+1 = S hn (u n ) b T order p ũ n+1 = S hn (u n ) bt order p S h (u) advances the numerical solution u n+1 = S hn (u n ), n 0, u 0 = U. Kansas Dec 2002 p.10/32
18 Approximation of a Dynamical System with a Variable Step-Size Runge-Kutta Embedded Pair Consider embedded Runge-Kutta pair with p p = 1. c A u n+1 = S hn (u n ) b T order p ũ n+1 = S hn (u n ) bt order p S h (u) advances the numerical solution u n+1 = S hn (u n ), n 0, u 0 = U. Note Does not define dynamical system on IR d as h n varies with n. Kansas Dec 2002 p.10/32
19 Local error approximation With user-defined tolerance, 0 < τ 1, step h n chosen by E(u n, h n ) τ, where E(u n, h n ) = 1 h ρ (u n+1 ũ n+1 ). n with ρ = 0 error per step (EPS) or ρ = 1 error per unit step (EPUS). Kansas Dec 2002 p.11/32
20 Local error approximation With user-defined tolerance, 0 < τ 1, step h n chosen by E(u n, h n ) τ, where E(u n, h n ) = 1 h ρ (u n+1 ũ n+1 ). n with ρ = 0 error per step (EPS) or ρ = 1 error per unit step (EPUS). Algorithm attempts to ensure E(u n, h n ) γτ, γ (0, 1) safety factor Leads to trouble near fixed points since f(u n ) = 0 implies E(u n, h n ) = 0. Kansas Dec 2002 p.11/32
21 Let h k n be a candidate for h n. Step-changing Algorithm If error control condition is satisfied, set h n = h k n, update solution and set [ ( ) ] h 0 γτ (1/p) n+1 = min h max, h n. E(u n, h n ) If error control condition not satisfied set h k+1 n = ( γτ E(u n, h k n) ) (1/p) h k n. p = min(p, p) + 1 ρ γ (0, 1) safety factor. Kansas Dec 2002 p.12/32
22 Adaptive Time-Stepping Algorithm as Dynamical Let τ > 0 be a (user-defined) error tolerance. An acceptable step-size h for u IR d satisfies S h (u) S h (u) τ. The maximum step-size (independent of τ) is h max. System Kansas Dec 2002 p.13/32
23 Adaptive Time-Stepping Algorithm as Dynamical Let τ > 0 be a (user-defined) error tolerance. An acceptable step-size h for u IR d satisfies S h (u) S h (u) τ. System The maximum step-size (independent of τ) is h max. Construct the sequence {(u n, h n )} n 0 using the algorithmic map S τ : IR d (0, h max ] IR d (0, h max ] (u n+1, h n+1 ) = S τ (u n, h n ). S τ finds an acceptable step-size h n and advances the solution. Kansas Dec 2002 p.13/32
24 Variable Time-Stepping Algorithm as Dynamical System S τ : IR d (0, h max ] IR d (0, h max ] (u n+1, h n+1 ) = S τ (u n, h n ). Kansas Dec 2002 p.14/32
25 Variable Time-Stepping Algorithm as Dynamical System S τ : IR d (0, h max ] IR d (0, h max ] (u n+1, h n+1 ) = S τ (u n, h n ). S τ is discontinuous. Kansas Dec 2002 p.14/32
26 Variable Time-Stepping Algorithm as Dynamical S τ : IR d (0, h max ] IR d (0, h max ] S τ is discontinuous. (u n+1, h n+1 ) = S τ (u n, h n ). System (A. Stuart & H. Lamba) There exists γ (0, 1) and C > 0 such that ( ) τ γ h C. f(u) Kansas Dec 2002 p.14/32
27 Variable Time-Stepping Algorithm as Dynamical S τ : IR d (0, h max ] IR d (0, h max ] S τ is discontinuous. (u n+1, h n+1 ) = S τ (u n, h n ). System (A. Stuart & H. Lamba) There exists γ (0, 1) and C > 0 such that ( ) τ γ h C. f(u) (G. Hall & D.J. Higham) f(u) 0 stability restricts the step-size. Kansas Dec 2002 p.14/32
28 Stable Fixed Point Example Consider the method RK1(2) applied to the linear system u = [ ], u = [ u 1 u 2 ], u(0) = [ 1, 10 4] T y n y(t) RK1(2) method with Standard (0, 0) stable fixed point. y y 1 For this method, the numerical solution gives persistent spurious oscillations and the y 1 component has O(τ) oscillation about the fixed point. Kansas Dec 2002 p.15/32
29 u = Methods RK2(3) and RK4(5) applied to Saddle ( ) u, u = ( u 1 u 2 ) Point Example, u(0) = ( 0.99, 10 10) T. Kansas Dec 2002 p.16/32
30 u = Methods RK2(3) and RK4(5) applied to Saddle ( ) u, u = ( u 1 u 2 ) Point Example, u(0) = ( 0.99, 10 10) T RK2(3) with Standard y n y(t) RK2(3) numerical solution does not pass close to fixed point or the local unstable manifold. y y 1 Kansas Dec 2002 p.16/32
31 u = Methods RK2(3) and RK4(5) applied to Saddle ( ) u, u = ( u 1 u 2 ) Point Example, u(0) = ( 0.99, 10 10) T RK2(3) with Standard y n y(t) RK2(3) numerical solution does not pass close to fixed point or the local unstable manifold. y RK4(5) with Standard y n y(t) y 1 RK4(5) has spurious oscillations about the unstable manifold. Numerical solution can ultimately end up either side of 0.02 the unstable manifold. y y 1 Kansas Dec 2002 p.16/32
32 Stable Saddle Point Example! Consider 2-stage RK1(2) method with stability domain 0 1/2 1/ Im(z) 1 0 R(z) < Re(z) Kansas Dec 2002 p.17/32
33 Stable Saddle Point Example! Consider 2-stage RK1(2) method with stability domain 0 1/2 1/ Im(z) 1 0 R(z) < 1 1 Apply this method to u = ( Re(z) ) u, u(0) = U IR 2. Kansas Dec 2002 p.17/32
34 Stable Saddle Point Example! 2 x 10 4 Phase Space using RK1(2). u(0) = (0.1,1e 05) 2 x 10 4 Phase Space using RK1(2). u(0) = (0.1,1e 05) 2.5 Step size Approximate Exact Acceptable Rejection u 2 (t) 0 u 2 (t) 0 h Approximate Exact u (t) u (t) 1 x t In an O(τ)-neighbourhood of the origin, the step-size oscillates about 1. Numerical solution near stable manifold becomes trapped near fixed point. Spurious stable invariant object in numerical flow. Kansas Dec 2002 p.18/32
35 Stable Saddle Point Example! Stability function is R(z) = 1 + z z 2. With z = hλ and here λ = ±1. Consider fixed step-size. For h (0, 1), the numerical manifolds are W s h (0) = {(x, y) IR2 y = 0} and W u h (0) = {(x, y) IR2 x = 0}. For h (1, 2), the numerical manifolds are W s h (0) = {(x, y) IR2 x = 0} and W u h (0) = {(x, y) IR2 y = 0}. Kansas Dec 2002 p.19/32
36 Stable Saddle Point Example! Stability function is R(z) = 1 + z z 2. With z = hλ and here λ = ±1. Consider fixed step-size. For h (0, 1), the numerical manifolds are W s h (0) = {(x, y) IR2 y = 0} and W u h (0) = {(x, y) IR2 x = 0}. For h (1, 2), the numerical manifolds are W s h (0) = {(x, y) IR2 x = 0} and W u h (0) = {(x, y) IR2 y = 0}. When h crosses 1, the manifolds are reversed. In adaptive algorithm this creates a chaotic attractor which persists for all τ > 0. Important to keep the step-size below linear (un)stability limits. Kansas Dec 2002 p.19/32
37 Approximation of Local Unstable Manifolds with a Variable Step-Size Runge-Kutta Pair Let σ > 0. The local unstable set of û, Wτ u,σ (û), is the set of (u, h) such that there exists a backward orbit under S τ {(u n, h n )} n=0 IR d (0, h max ], such that (u, h) = (u 0, h 0 ), u n û σ n 0, and u n û as n. Let σ δ > 0. W u,σ τ,δ (û) is the set of (u, h) such that there exists a finite backward orbit under S τ {(u n, h n )} N n=0 IR d (0, h max ], such that (u, h) = (u 0, h 0 ), u n û σ u N W u,δ h max (û). n = 0,..., N, and Kansas Dec 2002 p.20/32
38 Approximation of Local Unstable Manifolds with a Variable Step-Size Runge-Kutta Pair σ u* σ ~ W u (u*) W u τ (u*) W u,σ~ h max (u*) Kansas Dec 2002 p.20/32
39 An adaptive "Stable Manifold" Theorem Lemma There exists δ = O(τ) and h max sufficiently small and independent of τ such that W u,σ τ,δ u,σ (û) = Wτ (û). Theorem Apply an RKp( p) method with p p = 1 to u = f(u). Let û be a hyperbolic equilibrium and f C max{p, p} (IR d ). Then σ +, H + > 0 such that for h max (0, H + ) & σ (0, σ + ) d H (W u,σ (û), P u W u,σ τ (û)) 0 as τ 0 where P u : (u, h) IR d+1 u IR d is the projection operator. That is, the local unstable manifolds of the dynamical system and the unstable set of the RKp( p) are close for small τ. Kansas Dec 2002 p.21/32
40 Phase Space Stability Error Control Standard algorithm performs well during finite-time integration with fixed initial condition. Kansas Dec 2002 p.22/32
41 Phase Space Stability Error Control Standard algorithm performs well during finite-time integration with fixed initial condition. However unless h max is less than the linear stability limit the algorithm admits spurious fixed points; Kansas Dec 2002 p.22/32
42 Phase Space Stability Error Control Standard algorithm performs well during finite-time integration with fixed initial condition. However unless h max is less than the linear stability limit the algorithm admits spurious fixed points; performs badly around a stable fixed point; Kansas Dec 2002 p.22/32
43 Phase Space Stability Error Control Standard algorithm performs well during finite-time integration with fixed initial condition. However unless h max is less than the linear stability limit the algorithm admits spurious fixed points; performs badly around a stable fixed point; performs badly near saddle points. Kansas Dec 2002 p.22/32
44 Phase Space Stability Error Control Standard algorithm performs well during finite-time integration with fixed initial condition. However unless h max is less than the linear stability limit the algorithm admits spurious fixed points; performs badly around a stable fixed point; performs badly near saddle points. Don t want to restrict h max on whole phase space because of poor behaviour near fixed points Kansas Dec 2002 p.22/32
45 Phase Space Stability Error Control Standard algorithm performs well during finite-time integration with fixed initial condition. However unless h max is less than the linear stability limit the algorithm admits spurious fixed points; performs badly around a stable fixed point; performs badly near saddle points. Don t want to restrict h max on whole phase space because of poor behaviour near fixed points Don t want to introduce expensive algorithm to compute linear stability limit near fixed points. Kansas Dec 2002 p.22/32
46 Phase Space Stability Error Control Standard algorithm performs well during finite-time integration with fixed initial condition. However unless h max is less than the linear stability limit the algorithm admits spurious fixed points; performs badly around a stable fixed point; performs badly near saddle points. Introduce new phase space based error control to automatically control the step-size relative to the stability limit. Kansas Dec 2002 p.22/32
47 Phase Space (PS θ ) Error Control We demand at each step the phase space (PS θ ) error control u n+1 u n h n [(1 θ)f(u n ) + θf(u n+1 )] ϕh n (1 θ)f(u n ) + θf(u n+1 ), ϕ (0, 1). Kansas Dec 2002 p.23/32
48 Phase Space (PS θ ) Error Control We demand at each step the phase space (PS θ ) error control u n+1 u n h n [(1 θ)f(u n ) + θf(u n+1 )] ϕh n (1 θ)f(u n ) + θf(u n+1 ), ϕ (0, 1). h n (1 θ)f(u n ) + θf(u n+1 ) is approximation to arc length evolved over step. Kansas Dec 2002 p.23/32
49 Phase Space (PS θ ) Error Control We demand at each step the phase space (PS θ ) error control u n+1 u n h n [(1 θ)f(u n ) + θf(u n+1 )] ϕh n (1 θ)f(u n ) + θf(u n+1 ), ϕ (0, 1). h n (1 θ)f(u n ) + θf(u n+1 ) is approximation to arc length evolved over step. So PS θ error control bounds an approximation to local error by a fraction ϕ of an approximation to solution arc length in phase space. So is a phase space error control. Kansas Dec 2002 p.23/32
50 Phase Space (PS θ ) Error Control We demand at each step the phase space (PS θ ) error control u n+1 u n h n [(1 θ)f(u n ) + θf(u n+1 )] ϕh n (1 θ)f(u n ) + θf(u n+1 ), ϕ (0, 1). h n (1 θ)f(u n ) + θf(u n+1 ) is approximation to arc length evolved over step. So PS θ error control bounds an approximation to local error by a fraction ϕ of an approximation to solution arc length in phase space. So is a phase space error control. Will show it also acts as a stability control. Kansas Dec 2002 p.23/32
51 Phase Space (PS θ ) Error Control We demand at each step the phase space (PS θ ) error control u n+1 u n h n [(1 θ)f(u n ) + θf(u n+1 )] ϕh n (1 θ)f(u n ) + θf(u n+1 ), ϕ (0, 1). h n (1 θ)f(u n ) + θf(u n+1 ) is approximation to arc length evolved over step. So PS θ error control bounds an approximation to local error by a fraction ϕ of an approximation to solution arc length in phase space. So is a phase space error control. Will show it also acts as a stability control. Will combine this error control with standard error control; and demand both are satisfied at every step. Kansas Dec 2002 p.23/32
52 Key Features of this Error Control Negligible additional computation is needed; Kansas Dec 2002 p.24/32
53 Key Features of this Error Control Negligible additional computation is needed; Away from fixed points the standard error control is sufficient to ensure that the PS θ condition is satisfied. Kansas Dec 2002 p.24/32
54 Key Features of this Error Control Negligible additional computation is needed; Away from fixed points the standard error control is sufficient to ensure that the PS θ condition is satisfied. Prevents spurious fixed points; Kansas Dec 2002 p.24/32
55 Key Features of this Error Control Negligible additional computation is needed; Away from fixed points the standard error control is sufficient to ensure that the PS θ condition is satisfied. Prevents spurious fixed points; Forces convergence to stable fixed points; Kansas Dec 2002 p.24/32
56 Key Features of this Error Control Negligible additional computation is needed; Away from fixed points the standard error control is sufficient to ensure that the PS θ condition is satisfied. Prevents spurious fixed points; Forces convergence to stable fixed points; Gives stable step-size sequence with suitable step-size selection mechanism Kansas Dec 2002 p.24/32
57 Key Features of this Error Control Negligible additional computation is needed; Away from fixed points the standard error control is sufficient to ensure that the PS θ condition is satisfied. Prevents spurious fixed points; Forces convergence to stable fixed points; Gives stable step-size sequence with suitable step-size selection mechanism Good behaviour near saddle points Kansas Dec 2002 p.24/32
58 Step-size selection R(u n, h n ) = u n+1 u n h n [(1 θ)f(u n ) + θf(u n+1 )] h n (1 θ)f(u n ) + θf(u n+1 ) Next step: h θ n+1 = ( χϕ Order p 2 and θ 1/2 q = 1; Order p 3 and θ = 1/2 q = 2; χ (0, 1) is safety factor. ϕ. R(u n,h n )) 1/ qhn, where by Taylor series Kansas Dec 2002 p.25/32
59 Step-size selection R(u n, h n ) = u n+1 u n h n [(1 θ)f(u n ) + θf(u n+1 )] h n (1 θ)f(u n ) + θf(u n+1 ) Next step: h θ n+1 = ( χϕ Order p 2 and θ 1/2 q = 1; Order p 3 and θ = 1/2 q = 2; ϕ. R(u n,h n )) 1/ qhn, where by Taylor series χ (0, 1) is safety factor. New step-size selected as ] h n+1 = min [h s n+1, h θ n+1, αh n, where h s n+1 given by standard time-stepping strategy. α > 1 is a maximum step-size ratio, α = 5. Kansas Dec 2002 p.25/32
60 Linear system Theorem Consider forward Euler method under PS θ error control in with ϕ θ/(1 θ) applied to u t = Λu, Λ = Diag[λ 1, λ 2,, λ d ], u(0) = u 0 IR d where λ 1 < λ 2 < < λ d < 0. Then u n 0 as n with: Kansas Dec 2002 p.26/32
61 Linear system Theorem Consider forward Euler method under PS θ error control in with ϕ θ/(1 θ) applied to u t = Λu, Λ = Diag[λ 1, λ 2,, λ d ], u(0) = u 0 IR d where λ 1 < λ 2 < < λ d < 0. Then u n 0 as n with: 1. u n d 0 monoty as n ; Kansas Dec 2002 p.26/32
62 Linear system Theorem Consider forward Euler method under PS θ error control in with ϕ θ/(1 θ) applied to u t = Λu, Λ = Diag[λ 1, λ 2,, λ d ], u(0) = u 0 IR d where λ 1 < λ 2 < < λ d < 0. Then u n 0 as n with: 1. u n d 0 monoty as n ; 2. If λ i θ(1+ϕ) ϕ λ d, then u n i 0 & un i u n d 0 mono ty as n ; Kansas Dec 2002 p.26/32
63 Linear system Theorem Consider forward Euler method under PS θ error control in with ϕ θ/(1 θ) applied to u t = Λu, Λ = Diag[λ 1, λ 2,, λ d ], u(0) = u 0 IR d where λ 1 < λ 2 < < λ d < 0. Then u n 0 as n with: 1. u n d 0 monoty as n ; 2. If λ i θ(1+ϕ) ϕ 3. If θ(1+ϕ) ϕ λ d, then u n i 0 & un i u 0 mono ty as n ; n d [ ] λ d > λ i 2θ(1+ϕ) ϕ 1 λ d, then u n i 0 & un i u n 0; d Kansas Dec 2002 p.26/32
64 Linear system Theorem Consider forward Euler method under PS θ error control in with ϕ θ/(1 θ) applied to u t = Λu, Λ = Diag[λ 1, λ 2,, λ d ], u(0) = u 0 IR d where λ 1 < λ 2 < < λ d < 0. Then u n 0 as n with: 1. u n d 0 monoty as n ; 2. If λ i θ(1+ϕ) ϕ 3. If θ(1+ϕ) ϕ λ d, then u n i 0 & un i u 0 mono ty as n ; n d [ ] λ d > λ i 2θ(1+ϕ) ϕ 1 λ d, then u n i 0 & un i u n 0; d 4. Otherwise u n i 0 as n with lim sup n u n i u n d < ϕ θ ϕ 1+ϕ ; Kansas Dec 2002 p.26/32
65 Linear system Theorem Consider forward Euler method under PS θ error control in with ϕ θ/(1 θ) applied to u t = Λu, Λ = Diag[λ 1, λ 2,, λ d ], u(0) = u 0 IR d where λ 1 < λ 2 < < λ d < 0. Then u n 0 as n with: 1. u n d 0 monoty as n ; 2. If λ i θ(1+ϕ) ϕ 3. If θ(1+ϕ) ϕ λ d, then u n i 0 & un i u 0 mono ty as n ; n d [ ] λ d > λ i 2θ(1+ϕ) ϕ 1 λ d, then u n i 0 & un i u n 0; d 4. Otherwise u n u i 0 as n with lim sup n i n u n < ϕ d θ ; ϕ 1+ϕ 5. Let θ n be angle between u n and [0,..., 0, 1] IR d. If u 0 d 0 lim inf cos θ n 1 1 (d 1)ϕ2 n 2 θ 2 + O(ϕ3 ). Kansas Dec 2002 p.26/32
66 Remarks Can extend to arbitrary norms. Kansas Dec 2002 p.27/32
67 Can extend to arbitrary norms. Remarks Can extend to u = Au (i.e. non-diagonal). and nonlinear hyperbolic equilibria. Kansas Dec 2002 p.27/32
68 Can extend to arbitrary norms. Remarks Can extend to u = Au (i.e. non-diagonal). and nonlinear hyperbolic equilibria. Extends to non-stiff saddle points. Kansas Dec 2002 p.27/32
69 Can extend to arbitrary norms. Remarks Can extend to u = Au (i.e. non-diagonal). and nonlinear hyperbolic equilibria. Extends to non-stiff saddle points. Can extend to arbitrary methods. Kansas Dec 2002 p.27/32
70 Can extend to arbitrary norms. Remarks Can extend to u = Au (i.e. non-diagonal). and nonlinear hyperbolic equilibria. Extends to non-stiff saddle points. Can extend to arbitrary methods. Bound in 4 independent of stiffness/eigenvalues and can be made arbitray small by reducing ϕ. Kansas Dec 2002 p.27/32
71 Can extend to arbitrary norms. Remarks Can extend to u = Au (i.e. non-diagonal). and nonlinear hyperbolic equilibria. Extends to non-stiff saddle points. Can extend to arbitrary methods. Bound in 4 independent of stiffness/eigenvalues and can be made arbitray small by reducing ϕ. The exact solution is tangent to [0, 0,..., 0, 1] at fixed point, so 5 gives bound on angle between exact and numerical solutions at fixed point. Reducing ϕ makes angle arbitrarily small (independent of the stiffness/eigenvalues). Kansas Dec 2002 p.27/32
72 Proof Forward Euler method gives u n+1 = R(h n A)u n = Diag[1 + h n λ 1,, 1 + h n λ d ]u n. Kansas Dec 2002 p.28/32
73 Proof Forward Euler method gives u n+1 = R(h n A)u n = Diag[1 + h n λ 1,, 1 + h n λ d ]u n. With -norm, PS θ error control becomes θh 2 nλ 2 1 un 1 θh 2 nλ 2 2 un 2. θh 2 nλ 2 d un d ϕh n λ 1 (1 + θλ 1 h n )u n 1 λ 2 (1 + θλ 2 h n )u n 2. λ d (1 + θλ d h n )u n d. Note Unlike standard control not trivially true at near point. Kansas Dec 2002 p.28/32
74 Proof Forward Euler method gives u n+1 = R(h n A)u n = Diag[1 + h n λ 1,, 1 + h n λ d ]u n. With -norm, PS θ error control becomes θh 2 nλ 2 1 un 1 λ 1 (1 + θλ 1 h n )u n 1 θh 2 nλ 2 2 un 2 λ ϕh n 2 (1 + θλ 2 h n )u n 2.. θh 2 nλ 2 d un d λ d (1 + θλ d h n )u n d hence for some i {1, 2,..., d} h n θλ 2 i u n i ϕλ i 1 + θλ i h n u n i.. Kansas Dec 2002 p.28/32
75 Proof True if and only if h n θλ 2 i u n i ϕλ i 1 + θλ i h n u n i. h n ϕ λ i θ(1 + ϕ) ϕ λ d θ(1 + ϕ). Kansas Dec 2002 p.28/32
76 Proof True if and only if h n θλ 2 i u n i ϕλ i 1 + θλ i h n u n i. h n ϕ λ i θ(1 + ϕ) ϕ λ d θ(1 + ϕ). So the d th condition always holds and monotonic convergence in this component follows. Kansas Dec 2002 p.28/32
77 Proof True if and only if h n θλ 2 i u n i ϕλ i 1 + θλ i h n u n i. h n ϕ λ i θ(1 + ϕ) ϕ λ d θ(1 + ϕ). So the d th condition always holds and monotonic convergence in this component follows. Prove 2 and 3 by showing that d th condition implies (monotonic) convergence for these components. Kansas Dec 2002 p.28/32
78 Proof True if and only if h n θλ 2 i u n i ϕλ i 1 + θλ i h n u n i. h n ϕ λ i θ(1 + ϕ) ϕ λ d θ(1 + ϕ). So the d th condition always holds and monotonic convergence in this component follows. Prove 2 and 3 by showing that d th condition implies (monotonic) convergence for these components. Prove 4 by showing that failure of the i-th condition bounds yi n/yn d and hard work. 5 follows on noting that in limit n all components bounded in terms of d th component. Kansas Dec 2002 p.28/32
79 Stable Fixed Point Example (Revisited) 10 4 RK1(2) method with Standard 10 6 y n y(t) 10 8 y y 1 Recall solution with standard algorithm Kansas Dec 2002 p.29/32
80 Stable Fixed Point Example (Revisited) 10 4 RK1(2) with PS θ & Standard 10 5 y n y(t) 10 6 y y 1 x 10 4 With PS θ spurious oscillation is removed Kansas Dec 2002 p.29/32
81 Stable Fixed Point Example (Revisited) 10 4 RK1(2) with PS θ & Standard 10 0 RK1(2) with PS θ & Standard 10 5 y n y(t) Linear stability limit(λ =5) 10 6 y Step size PS θ Standard Rejected h n for Standard y 1 x t With PS θ spurious oscillation is removed Step-size is kept below stability limit. Kansas Dec 2002 p.29/32
82 Stable Fixed Point Example (Revisited) 10 4 RK1(2) with PS θ & Standard 10 0 RK1(2) with PS θ & Standard RK1(2) with PS θ & Standard,φ =10 1, τ = y n y(t) Linear stability limit(λ =5) y Step size 10 1 h x 10 3 PS θ Standard Rejected h n for Standard y 1 x t With PS θ spurious oscillation is removed Step-size is kept below stability limit. Step-sizes bounded near fixed point. PS θ only determines step-size near fixed point. Y X Kansas Dec 2002 p.29/32
83 Saddle Point Example (Revisited) 10 2 RK2(3) with Standard y n y(t) y y 1 Recall solution with RK2(3) standard algorithm Kansas Dec 2002 p.30/32
84 Saddle Point Example (Revisited) 10 2 RK2(3) with PS θ y n y(t) y y 1 With PS θ spurious oscillation is removed Kansas Dec 2002 p.30/32
85 Saddle Point Example (Revisited) 10 2 RK2(3) with PS θ 10 1 RK2(3) with PS θ & Standard 10 0 y n y(t) PS θ Standard Rejected h n for Standard Linear stability limit(λ =1) y Step size y t With PS θ spurious oscillation is removed Step-size is kept below stability limit. Kansas Dec 2002 p.30/32
86 Saddle Point Example (Revisited) 10 2 RK2(3) with PS θ 10 1 RK2(3) with PS θ & Standard RK2(3) with PS θ & Standard,φ =10 1, τ = y n y(t) PS θ Standard Rejected h n for Standard Linear stability limit(λ =1) y Step size h y t With PS θ spurious oscillation is removed Step-size is kept below stability limit. Step-sizes bounded near fixed point. PS θ only determines step-size near fixed point Y X Kansas Dec 2002 p.30/32
87 Large Average Step-Sizes RK1(2) Method applied to RK1(2) method with PS θ u = u. h stability limit for λ = 1 1 maximum step size (:) stability limit for λ = 10 2 average step size ( ) Step-size oscillates stability limit for λ 3 = 100 minimum step size ( ) χφ Kansas Dec 2002 p.31/32
88 Large Average Step-Sizes RK1(2) Method applied to u = Step-size oscillates. u. All steps below λ = 1 stability limit; monotonic convergence of this component. h stability limit for λ 1 = 1 stability limit for λ 2 = 10 RK1(2) method with PS θ maximum step size (:) stability limit for λ 3 = 100 average step size ( ) minimum step size ( ) χφ Kansas Dec 2002 p.31/32
89 RK1(2) Method applied to u = Step-size oscillates. u. All steps below λ = 1 stability limit; monotonic convergence of this component. Large Average Step-Sizes h stability limit for λ 1 = 1 stability limit for λ 2 = 10 RK1(2) method with PS θ maximum step size (:) stability limit for λ 3 = 100 average step size ( ) minimum step size ( ) χφ Average step-size below λ = 10 stability limit, but for ϕ > 0.1 some step-sizes above this limit. Kansas Dec 2002 p.31/32
90 RK1(2) Method applied to u = Step-size oscillates. u. All steps below λ = 1 stability limit; monotonic convergence of this component. Large Average Step-Sizes h stability limit for λ 1 = 1 stability limit for λ 2 = 10 RK1(2) method with PS θ maximum step size (:) stability limit for λ 3 = 100 average step size ( ) minimum step size ( ) χφ Average step-size below λ = 10 stability limit, but for ϕ > 0.1 some step-sizes above this limit. Except for ϕ tiny, average and maximum step-sizes above λ = 100 stability limit. But convergence to fixed point enforced. Kansas Dec 2002 p.31/32
91 Standard algorithm behaves poorly near saddle points. Stiff methods do not resolve problem for saddles. Conclusions Kansas Dec 2002 p.32/32
92 Standard algorithm behaves poorly near saddle points. Stiff methods do not resolve problem for saddles. Conclusions PS θ proved to give correct behaviour near stable fixed points gives correct behaviour near non-stiff saddles Kansas Dec 2002 p.32/32
93 Standard algorithm behaves poorly near saddle points. Stiff methods do not resolve problem for saddles. Conclusions PS θ proved to give correct behaviour near stable fixed points Ongoing gives correct behaviour near non-stiff saddles PS θ currently being implemented in standard ODE solver Kansas Dec 2002 p.32/32
CS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationMath 660 Lecture 4: FDM for evolutionary equations: ODE solvers
Math 660 Lecture 4: FDM for evolutionary equations: ODE solvers Consider the ODE u (t) = f(t, u(t)), u(0) = u 0, where u could be a vector valued function. Any ODE can be reduced to a first order system,
More informationMTH 452/552 Homework 3
MTH 452/552 Homework 3 Do either 1 or 2. 1. (40 points) [Use of ode113 and ode45] This problem can be solved by a modifying the m-files odesample.m and odesampletest.m available from the author s webpage.
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More information= f(u), u(0) = U (1.1)
THE EFFECTIVE STABILITY OF ADAPTIVE TIMESTEPPING ODE SOLVERS HARBIR LAMBA Abstract. We consider the behaviour of certain adaptive timestepping methods, based upon embedded explicit Runge-Kutta pairs, when
More informationWhat we ll do: Lecture 21. Ordinary Differential Equations (ODEs) Differential Equations. Ordinary Differential Equations
What we ll do: Lecture Ordinary Differential Equations J. Chaudhry Department of Mathematics and Statistics University of New Mexico Review ODEs Single Step Methods Euler s method (st order accurate) Runge-Kutta
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More information8.1 Introduction. Consider the initial value problem (IVP):
8.1 Introduction Consider the initial value problem (IVP): y dy dt = f(t, y), y(t 0)=y 0, t 0 t T. Geometrically: solutions are a one parameter family of curves y = y(t) in(t, y)-plane. Assume solution
More informationSelected HW Solutions
Selected HW Solutions HW1 1 & See web page notes Derivative Approximations. For example: df f i+1 f i 1 = dx h i 1 f i + hf i + h h f i + h3 6 f i + f i + h 6 f i + 3 a realmax 17 1.7014 10 38 b realmin
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationModule 4: Numerical Methods for ODE. Michael Bader. Winter 2007/2008
Outlines Module 4: for ODE Part I: Basic Part II: Advanced Lehrstuhl Informatik V Winter 2007/2008 Part I: Basic 1 Direction Fields 2 Euler s Method Outlines Part I: Basic Part II: Advanced 3 Discretized
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
More informationChapter 5 Exercises. (a) Determine the best possible Lipschitz constant for this function over 2 u <. u (t) = log(u(t)), u(0) = 2.
Chapter 5 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rjl/fdmbook Exercise 5. (Uniqueness for
More informationSolving Models With Off-The-Shelf Software. Example Of Potential Pitfalls Associated With The Use And Abuse Of Default Parameter Settings
An Example Of Potential Pitfalls Associated With The Use And Abuse Of Default Parameter Settings Ric D. Herbert 1 Peter J. Stemp 2 1 Faculty of Science and Information Technology, The University of Newcastle,
More information2015 Holl ISU MSM Ames, Iowa. A Few Good ODEs: An Introduction to Modeling and Computation
2015 Holl Mini-Conference @ ISU MSM Ames, Iowa A Few Good ODEs: An Introduction to Modeling and Computation James A. Rossmanith Department of Mathematics Iowa State University June 20 th, 2015 J.A. Rossmanith
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationModeling and nonlinear parameter estimation with Kronecker product representation for coupled oscillators and spatiotemporal systems
Physica D 227 (27) 78 99 wwwelseviercom/locate/physd Modeling and nonlinear parameter estimation with Kronecker product representation for coupled oscillators and spatiotemporal systems Chen Yao, Erik
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
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 informationSolving Ordinary Differential Equations
Solving Ordinary Differential Equations Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Initial Value Problem Euler s Method Runge-Kutta Methods Multistep Methods
More informationOrdinary Differential Equations
Ordinary Differential Equations Introduction: first order ODE We are given a function f(t,y) which describes a direction field in the (t,y) plane an initial point (t 0,y 0 ) We want to find a function
More informationEquilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains
Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains J. Földes Department of Mathematics, Univerité Libre de Bruxelles 1050 Brussels, Belgium P. Poláčik School
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationPeriodic Sinks and Observable Chaos
Periodic Sinks and Observable Chaos Systems of Study: Let M = S 1 R. T a,b,l : M M is a three-parameter family of maps defined by where θ S 1, r R. θ 1 = a+θ +Lsin2πθ +r r 1 = br +blsin2πθ Outline of Contents:
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 5/Part A - 23 November,
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationStability regions of Runge-Kutta methods. Stephan Houben Eindhoven University of Technology
Stability regions of Runge-Kutta methods Stephan Houben Eindhoven University of Technology February 19, 2002 1 Overview of the talk 1. Quick review of some concepts 2. Stability regions 3. L-stability
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationAdaptive Integration of Stiff ODE
INFORMATICA, 2000, Vol. 11, No. 4, 371 380 371 2000 Institute of Mathematics and Informatics, Vilnius Adaptive Integration of Stiff ODE Raimondas ČIEGIS, Olga SUBOČ, Vilnius Gediminas Technical University
More informationDepartment of Economics. Issn Discussion paper 30/09
Department of Economics Issn 1441-5429 Discussion paper 30/09 Solving Macroeconomic Models with "Off-the-Shelf" Software: An Example of Potential Pitfalls Ric D. Herbert a and Peter J. Stemp b,* Abstract:
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationODEs. PHY 688: Numerical Methods for (Astro)Physics
ODEs ODEs ODEs arise in many physics problems Classifications: As with the other topics, there are a large number of different methods Initial value problems Boundary value problems Eigenvalue problems
More informationLorenz like flows. Maria José Pacifico. IM-UFRJ Rio de Janeiro - Brasil. Lorenz like flows p. 1
Lorenz like flows Maria José Pacifico pacifico@im.ufrj.br IM-UFRJ Rio de Janeiro - Brasil Lorenz like flows p. 1 Main goals The main goal is to explain the results (Galatolo-P) Theorem A. (decay of correlation
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 informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationInvestigation on the Most Efficient Ways to Solve the Implicit Equations for Gauss Methods in the Constant Stepsize Setting
Applied Mathematical Sciences, Vol. 12, 2018, no. 2, 93-103 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.711340 Investigation on the Most Efficient Ways to Solve the Implicit Equations
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationThe Existence of Chaos in the Lorenz System
The Existence of Chaos in the Lorenz System Sheldon E. Newhouse Mathematics Department Michigan State University E. Lansing, MI 48864 joint with M. Berz, K. Makino, A. Wittig Physics, MSU Y. Zou, Math,
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
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 informationODEs. PHY 688: Numerical Methods for (Astro)Physics
ODEs ODEs ODEs arise in many physics problems Classifications: As with the other topics, there are a large number of different methods Initial value problems Boundary value problems Eigenvalue problems
More informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationComputers, Lies and the Fishing Season
1/47 Computers, Lies and the Fishing Season Liz Arnold May 21, 23 Introduction Computers, lies and the fishing season takes a look at computer software programs. As mathematicians, we depend on computers
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
More informationComplex system approach to geospace and climate studies. Tatjana Živković
Complex system approach to geospace and climate studies Tatjana Živković 30.11.2011 Outline of a talk Importance of complex system approach Phase space reconstruction Recurrence plot analysis Test for
More informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More informationParallel Methods for ODEs
Parallel Methods for ODEs Levels of parallelism There are a number of levels of parallelism that are possible within a program to numerically solve ODEs. An obvious place to start is with manual code restructuring
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 informationON THE CONTINUITY OF GLOBAL ATTRACTORS
ON THE CONTINUITY OF GLOBAL ATTRACTORS LUAN T. HOANG, ERIC J. OLSON, AND JAMES C. ROBINSON Abstract. Let Λ be a complete metric space, and let {S λ ( ) : λ Λ} be a parametrised family of semigroups with
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationHandout 2: Invariant Sets and Stability
Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationStudy guide: Generalizations of exponential decay models
Study guide: Generalizations of exponential decay models Hans Petter Langtangen 1,2 Center for Biomedical Computing, Simula Research Laboratory 1 Department of Informatics, University of Oslo 2 Sep 13,
More informationStudy guide: Generalizations of exponential decay models. Extension to a variable coecient; Crank-Nicolson
Extension to a variable coecient; Forward and Backward Euler Study guide: Generalizations of exponential decay models Hans Petter Langtangen 1,2 Center for Biomedical Computing, Simula Research Laboratory
More informationNumerical Methods for Differential Equations
CHAPTER 5 Numerical Methods for Differential Equations In this chapter we will discuss a few of the many numerical methods which can be used to solve initial value problems and one-dimensional boundary
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More informationOrdinary Differential Equations: Initial Value problems (IVP)
Chapter Ordinary Differential Equations: Initial Value problems (IVP) Many engineering applications can be modeled as differential equations (DE) In this book, our emphasis is about how to use computer
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 informationIntroduction Knot Theory Nonlinear Dynamics Topology in Chaos Open Questions Summary. Topology in Chaos
Introduction Knot Theory Nonlinear Dynamics Open Questions Summary A tangled tale about knot, link, template, and strange attractor Centre for Chaos & Complex Networks City University of Hong Kong Email:
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
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 information