The Trapping Region of the Lorenz equations

Click For Summary
SUMMARY

The discussion centers on the trapping region of the Lorenz equations, specifically defined by the ellipsoid ρx² + σy² + σ(z - 2ρ)² < R. Participants explore the use of Lyapunov functions to identify trapping regions and the discrepancies encountered when applying alternative methods, such as defining radial distance r = √(x² + y² + z²). It is established that different Lyapunov functions can yield varying trapping regions, and while a minimal region can be identified, it does not guarantee uniqueness. The conversation highlights the complexity of determining minimal trapping regions and the role of radial velocity in these systems.

PREREQUISITES
  • Understanding of nonlinear systems of differential equations
  • Familiarity with Lyapunov functions and their applications
  • Knowledge of the Lorenz equations and their dynamics
  • Basic concepts of polar coordinates in mathematical analysis
NEXT STEPS
  • Research the properties and applications of Lyapunov functions in dynamical systems
  • Explore the derivation and implications of the Lorenz equations
  • Study the concept of minimal trapping regions in nonlinear dynamics
  • Learn about the conversion of Cartesian coordinates to polar coordinates in differential equations
USEFUL FOR

Mathematicians, physicists, and engineers working with dynamical systems, particularly those focusing on chaos theory and stability analysis in nonlinear differential equations.

Oliver321
Messages
59
Reaction score
5
I was dealing with nonlinear systems of differential equations like the Lorenz equations (https://en.wikipedia.org/wiki/Lorenz_system). Now there is a trapping region of this system defined by the ellipsoid ρx^2+σy^2+σ(z-2ρ)^2<R.
I wondered how this region is found and I found out that a Lyapunov function is used.
However, I tried to do it another way. I defined r=√(x^2+y^2+z^2). Now I thought: when d/dt r <0 the radial velocity points inward and so there is a trapping region. I calculated this but the result was completely different.
Why doesn’t it work with my way?

Thanks for every awnser!
 
Physics news on Phys.org
In the trapping region, isn't dr/dt often bigger than zero? The particles move around and make loops around the attractor, but dr/dt isn't constantly shrinking while doing this.
 
  • Like
Likes   Reactions: Oliver321
Office_Shredder said:
In the trapping region, isn't dr/dt often bigger than zero? The particles move around and make loops around the attractor, but dr/dt isn't constantly shrinking while doing this.

Yes that’s right, but on the edge of the trapping region all arrows should point inward (or at least be tangential to the boundary), so dr/dt should be less or equal to zero?
 
If you use different Liapunov functions, you will get different trapping regions.
 
  • Like
Likes   Reactions: Oliver321
pasmith said:
If you use different Liapunov functions, you will get different trapping regions.
But why? With both methods I can get a minimal region, which is a trapping region. Because it is minimal, it has to be the same?
 
Oliver321 said:
But why? With both methods I can get a minimal region, which is a trapping region. Because it is minimal, it has to be the same?

How do you know the region is minimal?

All you can conlcude from looking at a particular Liapunov function is that the attractor is in some region, because once trajextories enter it they cannot leave. But there's always the possibility that a different Liapunov function would give you a stirctly smaller region, and hence a closer bound on the attractor.
 
  • Like
Likes   Reactions: Oliver321
pasmith said:
How do you know the region is minimal?

All you can conlcude from looking at a particular Liapunov function is that the attractor is in some region, because once trajextories enter it they cannot leave. But there's always the possibility that a different Liapunov function would give you a stirctly smaller region, and hence a closer bound on the attractor.

Thanks for the answer!
So I thought that there is only one Liapunov function for a system. Like if i have a conservative system I can use a energy function to do the same things (and the Energy function is unique?). A liapunov function is a generalised energy function (at least so I was told). But if there are a few, this would be plausible.
Nevertheless the method with converting to polar coordinates should give a minimal region?
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 12 ·
Replies
12
Views
5K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
954
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K