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 condition reaction-diffusion system: u t = D u u + f(u, v),, x Ω, with boundary condition v t = D v v + g(u, v), All equation is in form of U t = F (U), where u can be a scalar or vector, spatial independent or dependent 1
Abstract Equation U t = F (U) Equilibrium solution: U 0 such that F (U 0 ) = 0, linearized operator: F (U 0 ) U 0 is stable if the eigenvalues of equation F (U 0 )w = λw are all with negative real parts. (Linear behavior) Since the equation U t = F (U) is approximately the linearized equation U t = F (U)(U U 0 ) near the equilibrium solution U = U 0, then U(t) C i exp(λ i t)φ i near U = U 0, where (λ i, φ i ) are the eigenvalue-eigenvector pairs of F (U 0 )w = λw. If λ 1 = max λ i, then U(t) C 1 exp(λ 1 t)φ 1 if U(0) U 0. 2
Example 1: Scalar equation u t = f(u), f(u 0 ) = 0, linearized operator: f (u 0 ) f (u 0 ) < 0 stable, otherwise unstable Example 2: Linear system ( ) u t = 2u + 3v, 2 3, linearized system: J = v t = 4u + 3v, 4 3 ( ) 2 3 Eigenvalue problem: w = λw, 4 3 eigen-pairs: ( λ 1 = ) 1, φ 1 = ((1, 1); ) λ 2 = 6, ( φ 2 ) = (3, 4) ( u(t) Solution: = c v(t) 1 e t 1 + c 1 2 e 6t 3 c 4 2 e 6t 3 4 perturbed from equilibrium. ) if 3
Example 3: Nonlinear ODE system: x t = 2x x 2 ( xy, 2 2x y x y t = 3y y 2, linearized operator: J = 2xy, 2y 3 2y 2 Equilibrium points: (0, 0), (2, 0), (3, 0) and (1, 1) For (1, 1), eigenvalues λ 1 = 1 + 2, φ 1 = (0.58, 0.82); λ 2 = 1 2, φ 2 = (0.58, 0.82) when per- ( ) x(t) So near the saddle point (1, 1), = y(t) ( ) ( ) c 2 e ( 1 2)t 0.58 1 + c 0.82 1 1 e ( 1+ 2)t turbed from equilibrium. ( 1 ) 1 ( 0.58 0.82 +c 1 e ( 1+ 2)t ) ( 0.58 0.82 ) 4
Example 4: scalar diffusion equation u t = Du xx, x (0, π), t > 0, u(t, 0) = u(t, π) = 0. Equilibrium solution: u(x) = 0. Linearized operator: F (u)[w] = Dw xx for w satisfying w(0) = w(π) = 0, Eigenvalue problem: Dw xx = λw, x (0, π), w(0) = w(π) = 0. Eigenvalue-eigenfunction pairs: λ m = Dm 2, φ m (x) = sin(mx), m = 1, 2,. Solution: u(x, t) = c m exp( Dm 2 t) sin(mx) c 1 exp( Dt) sin(x) near u = 0, and u(x) = 0 is a stable equilibrium solution 5
Example 5: scalar diffusive Malthusian equation u t = Du xx + au, x (0, π), t > 0, u(t, 0) = u(t, π) = 0. Equilibrium solution: u(x) = 0. Linearized operator: F (u)[w] = Dw xx + aw for w satisfying w(0) = w(π) = 0, Eigenvalue problem: Dw xx + aw = λw, x (0, π), w(0) = w(π) = 0. Eigenvalue-eigenfunction pairs: λ m = a Dm 2, φ m (x) = sin(mx), m = 1, 2,. Solution: u(x, t) = c m exp((a Dm 2 )t) sin(mx) c 1 exp((a D)t) sin(x) near u = 0. u(x) = 0 is a stable equilibrium solution if a < D, and it is unstable if a > D. (Thus a = D is a potential bifurcation point.) 6
Example 6: Fisher equation u t = Du xx + au(1 u), x (0, π), t > 0, u(t, 0) = u(t, π) = 0. Equilibrium solution: u(x) = 0. Linearized operator: F (u)[w] = Dw xx + aw for w satisfying w(0) = w(π) = 0, Eigenvalue problem: Dw xx + aw = λw, x (0, π), w(0) = w(π) = 0. Eigenvalue-eigenfunction pairs: λ m = a Dm 2, φ m (x) = sin(mx), m = 1, 2,. Solution: u(x, t) = c m exp((a Dm 2 )t) sin(mx) c 1 exp((a D)t) sin(x) near u = 0. u(x) = 0 is a stable equilibrium solution if a < D, and it is unstable if a > D. (Thus a = D is a bifurcation point.) 7
Fisher equation: (cont.) For a > D, Fisher equation has a positive equilibrium solution u a (x) (from local bifurcation theory, and global bifurcation of time-mapping). u a (x) is a stable equilibrium solution. (see lecture notes for a proof) 8
Diffusion is a stabilizing process? In Fisher equation: u t = Du xx + au(1 u), x (0, π), t > 0, u(t, 0) = u(t, π) = 0. Equilibrium solution u = 0 is also an equilibrium solution of u t = au(1 u), but an unstable one (since f (0) > 0). (eigenvalue λ = af (0) > 0) Equilibrium solution u(x) = 0 is stable for Fisher equation when a < D, and it is still unstable when a > D. (eigenvalues: λ m = af (0) Dm 2 ) Thus diffusion makes the eigenvalue more negative and diffusion in scalar equation has a stabilizing effect. (If stable for ODE, also stable for PDE; even unstable for ODE, still can be stable for PDE.) 9
Next model (most useful one): reaction-diffusion system u t = D u u xx + f(u, v), v t = D v v xx + g(u, v), u x (t, 0) = u x (t, 1) = 0, v x (t, 0) = v x (t, 1) = 0 If (u 0, v 0 ) is an equilibrium solution that f(u 0, v 0 ) = g(u 0, v 0 ) = 0, then (u 0, v 0 ) is also an equilibrium solution of u t = f(u, v), v t = g(u, v). Linearized operator: ODE: Jacobian J = ( du 0 PDE: diag(d u, d v ) + J = 0 d v v ( ) fu f v g u g v ) ( fu f + v g u g v ) Eigenvalue problem: D u φ + f u (u 0, v 0 )φ + f v (u 0, v 0 )ψ = λφ, D v ψ + g u (u 0, v 0 )φ + g v (u 0, v 0 )ψ = λψ, φ x (0) = φ x (1) = 0, ψ x (0) = ψ(1) = 0. 10
If all eigenvalues of J are with negative real parts (so (u 0, v 0 ) is a stable equilibrium solution for ODE), is (u 0, v 0 ) also a stable equilibrium solution for PDE? It seems so since the additional part is consist of diffusion operators only, and diffusion is supposed to stabilizing... But, as Alan Turing pointed out, (u 0, v 0 ) could be an unstable equilibrium solution for PDE even if it is stable for ODE! So diffusion has an unstable effect for such system. How is that possible? Let s calculate now... 11