Outline Outline 4 Definitions 4 Definitions 4 Phase Plane 4 - - PowerPoint PPT Presentation

outline outline 4 definitions 4 definitions 4 phase plane
SMART_READER_LITE
LIVE PREVIEW

Outline Outline 4 Definitions 4 Definitions 4 Phase Plane 4 - - PowerPoint PPT Presentation

Outline Outline 4 Definitions 4 Definitions 4 Phase Plane 4 Phase Plane 4 Attractor and Stability 4 Attractor and Stability 4 Bifurcation 4 Bifurcation ME 639-Turbulence G. Ahmadi ME 639-Turbulence G. Ahmadi Engineering Equation of


slide-1
SLIDE 1

1

  • G. Ahmadi

ME 639-Turbulence

  • G. Ahmadi

ME 639-Turbulence

Outline Outline 4 4Definitions Definitions 4 4Phase Plane Phase Plane 4 4Attractor and Stability Attractor and Stability 4 4Bifurcation Bifurcation

  • G. Ahmadi

ME 639-Turbulence

Simple Pendulum Simple Pendulum

mg

2

Engineering Systems Engineering Engineering Systems Systems Dynamical Dynamical Systems Systems State State State

θ θ & ,

Phase Plane Phase Plane Phase Plane

θ θ versus

  • f

Plot &

l

  • G. Ahmadi

ME 639-Turbulence

Simple Pendulum Simple Pendulum Simple Pendulum Equation of Motion Equation of Equation of Motion Motion System of Differential System of Differential Equations Equations

) f(x x t , = &

⎩ ⎨ ⎧ θ ω − ω ζω − = ω ω = θ sin 2

2

  • &

&

l g

2

  • =

ω

⎩ ⎨ ⎧ θ ω − ω ζω − = ω ω = θ

2

  • 2

& &

Linear Pendulum Linear Pendulum Linear Pendulum

slide-2
SLIDE 2

2

  • G. Ahmadi

ME 639-Turbulence

Time Does Not Time Does Not Appear Explicitly Appear Explicitly Autonomous Systems Autonomous Autonomous Systems Systems Nonautonomous Systems Nonautonomous Nonautonomous Systems Systems Time Appears Time Appears Explicitly Explicitly Orbit/ Trajectory Orbit/ Orbit/ Trajectory Trajectory Curves in Phase Space Curves in Phase Space Solutions of Equation Solutions of Equation

  • f Motion
  • f Motion
  • G. Ahmadi

ME 639-Turbulence

Vectors Tangent Vectors Tangent to Trajectories in to Trajectories in Phase Space Phase Space

Flow Flow Flow Null Cline Null Cline Null Cline

Lines in Phase Space Lines in Phase Space for for

xi = &

⎩ ⎨ ⎧ ζ θ ω − = ω = ω = θ 2 / sin

  • &

& Pendulum Pendulum

  • G. Ahmadi

ME 639-Turbulence

Intersection

  • f Clines

Intersection Intersection

  • f Clines
  • f Clines

Equilibrium Points Equilibrium Points

Poincare Section Poincare Map Poincare Section Poincare Section Poincare Map Poincare Map

  • G. Ahmadi

ME 639-Turbulence

Non-Wandering Set

A set of points that orbits starting form this set come arbitrary close and arbitrary often to any point in the set

Non Non-

  • Wandering Set

Wandering Set

A set of points that orbits starting form A set of points that orbits starting form this set come arbitrary close and this set come arbitrary close and arbitrary often to any point in the set arbitrary often to any point in the set

  • Fixed (stationary) points
  • Limit cycle (periodic)
  • Quasi periodic
  • Chaotic (bounded, non-periodic)
  • Fixed (stationary) points

Fixed (stationary) points

  • Limit cycle (periodic)

Limit cycle (periodic)

  • Quasi periodic

Quasi periodic

  • Chaotic (bounded, non

Chaotic (bounded, non-

  • periodic)

periodic)

slide-3
SLIDE 3

3

  • G. Ahmadi

ME 639-Turbulence

Nonlinear Dynamical Systems

  • Qualitative Behavior
  • Non-wandering sets
  • Stability of non-wandering sets
  • Changes in the number of non-wandering sets

Nonlinear Dynamical Systems Nonlinear Dynamical Systems

  • Qualitative Behavior

Qualitative Behavior

  • Non

Non-

  • wandering sets

wandering sets

  • Stability of non

Stability of non-

  • wandering sets

wandering sets

  • Changes in the number of non

Changes in the number of non-

  • wandering sets

wandering sets

Bifurcation

Appearance and Disappearance of Non-wandering sets

Bifurcation Bifurcation

Appearance and Disappearance of Appearance and Disappearance of Non Non-

  • wandering sets

wandering sets

  • G. Ahmadi

ME 639-Turbulence

Lyapunov (Marginal) Stability

A non-wandering set (NWS) is Lyapunov stable if every orbit starting in its neighborhood stays in its neighborhood.

Lyapunov Lyapunov (Marginal) Stability (Marginal) Stability

A non A non-

  • wandering set (NWS) is

wandering set (NWS) is Lyapunov Lyapunov stable if every orbit starting in its stable if every orbit starting in its neighborhood stays in its neighborhood. neighborhood stays in its neighborhood.

Asymptotic Stability

A NWS is asymptotically stable if in addition to Lyapunov stability every orbit in its neighborhood approaches the NWS.

Asymptotic Stability Asymptotic Stability

A NWS is asymptotically stable if in A NWS is asymptotically stable if in addition to addition to Lyapunov Lyapunov stability every orbit in stability every orbit in its neighborhood approaches the NWS. its neighborhood approaches the NWS.

  • G. Ahmadi

ME 639-Turbulence

Attractors

Asymptotically stable non-wandering sets are called attractors.

Attractors Attractors

Asymptotically stable non Asymptotically stable non-

  • wandering

wandering sets are called attractors. sets are called attractors.

Basin of Attraction

The set of all initial states that approach the attractor.

Basin of Attraction Basin of Attraction

The set of all initial states that The set of all initial states that approach the attractor. approach the attractor.

  • G. Ahmadi

ME 639-Turbulence

Given an Orbit Given an Orbit Given an Orbit ) f(x x t ,

  • =

&

) t (

  • x

xo(t) is asymptotically stable if ∆x(t) decays x xo

  • (t

(t) is asymptotically ) is asymptotically stable if stable if ∆ ∆x(t) decays x(t) decays x x x f x

⋅ = ∆ d ) ( d dt ) ( d

s t

s

e x x ∆ = ∆

λ

I x f λ − d d del

s

λ

Solution Solution

slide-4
SLIDE 4

4

  • G. Ahmadi

ME 639-Turbulence

Stable Orbit if all are negative Stable Orbit if all are negative Stable Orbit if all are negative

d d tr d d del < > x f x f

s

λ

Routh-Hurwitz Theorem For two-D all the roots are negative if Routh Routh-

  • Hurwitz Theorem

Hurwitz Theorem For two For two-

  • D all the roots

D all the roots are negative if are negative if

  • G. Ahmadi

ME 639-Turbulence

Two real eigenvalues

  • f the same sign

Two real Two real eigenvalues eigenvalues

  • f the same sign
  • f the same sign

Nodes Nodes

Pair of complex conjugate eigenvalues Pair of complex Pair of complex conjugate conjugate eigenvalues eigenvalues

Spirals Spirals

Two real eigenvalues with opposite sign Two real Two real eigenvalues eigenvalues with opposite sign with opposite sign

Saddle Saddle

  • G. Ahmadi

ME 639-Turbulence

A change in the number of attractors of a A change in the number of attractors of a nonlinear dynamical systems with the change nonlinear dynamical systems with the change

  • f a system parameter is called
  • f a system parameter is called bifurcation

bifurcation. .

Bifurcation Bifurcation is associated with the is associated with the change of stability of an change of stability of an attractor attractor. . In a bifurcation point, at least one In a bifurcation point, at least one eigenvalue eigenvalue of the

  • f the Jacobian

Jacobian will attain a will attain a zero real part. zero real part.

  • G. Ahmadi

ME 639-Turbulence

In a stationary bifurcation, a single In a stationary bifurcation, a single real real eigenvalue eigenvalue crosses the boundary crosses the boundary

  • f stability.
  • f stability.
slide-5
SLIDE 5

5

  • G. Ahmadi

ME 639-Turbulence

Hopf Hopf bifurcation occurs when a bifurcation occurs when a conjugated complex pair crosses the conjugated complex pair crosses the boundary of stability. boundary of stability.

  • G. Ahmadi

ME 639-Turbulence

u ) u ( u

2

− µ = &

For For µ µ<0, u=0 is a stable equilibrium solution. <0, u=0 is a stable equilibrium solution. For For µ µ>0, u=0 is an unstable equilibrium >0, u=0 is an unstable equilibrium solution, and u= solution, and u= ±µ ±µ1/2

1/2 are stable solutions.

are stable solutions. At At µ µ=0, a =0, a suppercritical suppercritical Pitchfork Pitchfork bifurcation occurs. bifurcation occurs.

Consider a dynamical system Consider a Consider a dynamical system dynamical system

  • G. Ahmadi

ME 639-Turbulence

For For µ µ>0, u=0 is an unstable equilibrium solution. >0, u=0 is an unstable equilibrium solution.

For For µ µ<0, u=0 is a stable equilibrium solution, <0, u=0 is a stable equilibrium solution, and u= and u= ± ±( (-

  • µ

µ) )1/2

1/2 are unstable solutions.

are unstable solutions. At At µ µ=0, a =0, a subcritical subcritical Pitchfork Pitchfork bifurcation occurs. bifurcation occurs.

u ) u ( u

2

+ µ = &

Consider a dynamical system Consider a Consider a dynamical system dynamical system

  • G. Ahmadi

ME 639-Turbulence

For For µ µ<0, u=0 is a stable equilibrium solution, <0, u=0 is a stable equilibrium solution, and u= and u=µ µ is an unstable solution. is an unstable solution.

u ) u ( u − µ = &

Consider a dynamical system Consider a Consider a dynamical system dynamical system At At µ µ=0, the two solution exchange stability =0, the two solution exchange stability and a and a transcritical transcritical bifurcation occurs. bifurcation occurs. For For µ µ>0, u=0 is an unstable solution, and u= >0, u=0 is an unstable solution, and u=µ µ is is a stable equilibrium a stable equilibrium solution. solution.