Blog

Learn something new with our collection of blog posts centered around various topics on PSpice® technology.

Lorenz Oscillator- The Lorenz System

“It may happen that small differences in the initial conditions produce very great ones in the final phenomena. A small error in the former will produce an enormous error in the future. Prediction becomes impossible”. Edward Lorenz

 

The meteorologist, Edward Lorenz, in 1961, discovered chaos theory inadvertently whilst modeling atmospheric convection and weather patterns. He published a paper in 1963 where he showed how the model produced chaotic behavior by changing the model from a linear to a nonlinear model. To save time during a simulation run on a very basic desktop computer, he started from the middle of a previous simulation session and to his surprise found the output was completely different from a previous output. By truncating five-digit simulation numbers down to three-digit numbers he discovered one of the principles of chaos theory “Sensitive to Initial Conditions (SIC)”. A small change in initial conditions i.e., reducing the number of places of decimals, produced a bigger change in the system as it evolves. What is often quoted in the literature is Lorenz’s famous quote in a presentation in 1972 to the American Association of Scientific research “Does the flap of a butterfly’s wings in Brazil set off a tornado in Texas”, which illustrates how very small change in the initial condition produces much bigger changes later. Edward Lorenz showed that for chaos to exist in autonomous ordinary differential equation systems, three variables and two quadratic non-linear elements must be present. However, since then the number of variables can be reduced to two for non-autonomous chaotic  system, and one, for discrete chaotic systems. The Lorenz weather system models clockwise and anticlockwise thermal convective air mass rotations in a closed container and the ordinary differential equations used in this model were originally created by Barry Saltzman (1962). Any n-th order nonlinear differential equation can be written as n separate coupled first-order ODEs. For example, the following coupled equations are a reduced form of the Navier-Stokes equations that model the Lorenz chaotic system that we will examine shortly:

 

                                             Lorenz_oscillator_equation1.png

 

By coupled we mean that the equations are linked to each other and the parameters. The original parameter used by Lorenz are: Prandtl number = 10, Rayleigh number, R = 28, and β = B = 8/3 and are all greater than zero. This simplified Lorenz model has since been proved to be chaotic by B Hassard and J Zhang, 1994 but also by Tucker in 1999. The parameters were given these names because Lorenz derived his equations from the Raleigh-Bénard convection flow (convection in a liquid heated in a container from below, generates Bénard cells). We rewrite these equations in a different form shortly, where we will express the variables as a function of time and introduce initial conditions. The Lorenz equations represent the movement of a horizontal layer of a fluid like air, being heated from below. The speed of the rising air mass, x, is positive for clockwise rotation and negative for anti-clockwise rotations. The temperature difference between descending and ascending air masses is y, and z is the deviation of the temperature profile from a linear scale.

 

3.1 Simulating the Lorenz system using ABM parts

 

To simulate the Lorenz system in PSpice, we rearrange the equations in integral form and introduce initial condition conditions x(t0), y(t0), z(t0) for each system variable. We do this because we will implement these equations using integrators using analog behavioral models (ABM) parts which will be replaced by electronic inverting integrators circuits later. The initial conditions will be dropped for simplicity later but will be included in the simulation model as initial conditions on the integrators. Hence, the first equation, P(y-x), becomes –P(x-y). We could use non-inverting integrators but it’s a little more complicated and left as an exercise to the interested reader. Later, we will see how to use the four-quadrant AD633 multiplier IC to realize the nonlinear product terms. This IC may also be configured to achieve integration using current outputs but that’s much later. The Lorenz equations with initial conditions, t0 are:

 

             Lorenz_oscillator_equation2.png

 

The simplest approach when investigating a chaotic system for the first time is to always use ABM parts, as shown in Figure 1 and later with actual integrated circuit models.  This approach is also useful when the circuit becomes complex and in some cases when the electronic part models cause convergence.

 

                Lorenz_oscillator_figure1.png

 

                                       Figure 1: The Lorenz ABM schematic

 

The ABM INTEG part has its gain set to {–Invert_tau}, where the time constant defined in a PARAM part. The gain of an inverting integrator is set negative and equal to the inverse of the time constant because in the next blog we will use inverting electronic part models. Making the gain as the inverse of the time constant enables us to set a frequency close to the main operating frequency of the system, and in the section on scaling later in this chapter, we will see the need for scaling and how to scale both amplitude and time variables. Figure 2 shows how the three variables vary over time

 

                   Lorenz_oscillator_figure2.png

                        Figure 2: The Lorenz x, y and z output

 

Figure 3 displays the spectra for the variables and is obtained by pressing the FFT icon when in the Probe environment (appears after simulation). Note we may express the spectral components as dB(V(x)*V(x)), to ensure the smaller spectral component are emphasized, otherwise the main spectral component of 1.32 Hz will tend to mask the other components. Select the Add Trace icon (a small graph-like icon), and then select dB( ) from the list of functions in the right panel. Select between the brackets and type v(x)*v(x). However, I just checked before going to press and the dB() function doesn’t now seem to allow the product of variables in the brackets so instead use 10*log10(v(x)*v(x)).  Increasing the simulation run time will increase the frequency resolution but obviously increases the simulation time. The z output is quite periodic and is therefore said to be chaotically weak.

 

                  Lorenz_oscillator_figure3.png

                                                            Figure 3: Power Spectra

 

 

After Probe simulation and to automate the process of placing multiple windows and additional plots, we use the log command which records all keys selected. Select log in Probe File menu. When this is selected, it will record  log file of all keys pressed when doing tasks such as separating the windows. This is useful where you intend to repeat the simulation for parameter variations, for example. If you change the integrator gain from 100 to 1, you will change the system operating frequency hence you will need to change the simulation settings. I have included a log file I created called Lorenz.cmd which will plot all the variables and is called from the File menu in Probe. Select Run Commands and then select Lorenz.cmd. This will plot Figure 1 and 2 but also it will plot the famous Lorenz strange attractor by changing time to the x variable.

Fri, 2017-06-30 04:06
anonymous
Offline
Last seen: 3 months 2 days ago
Joined: 2017-06-30 04:02

Nice Post

Thanks for this amazing info. I appreciate your hard work here. As a student I buy custom essays sometimes struggle with deadlines, especially when I have so many responsibilities. Luckily, there is this paper service may decide any paper tasks as soon as possible.

Thu, 2017-08-17 05:24 (Reply to #1)
toby
toby's picture
Offline
Last seen: 3 months 4 weeks ago
Joined: 2016-09-09 06:22

Chaos

Thank you Katejohn for your nice comments-much appreciated, and sorry for the delay. I will be posting more articles on chaos very soon, so stand by for practical chaos implementation. Thanks Toby

Thu, 2017-08-31 06:22
RAHU
Offline
Last seen: 3 months 3 days ago
Joined: 2017-08-31 06:14

rahu

As soon as the connection have been developed Robux generators could influence the resources for roblox.
 

Thu, 2017-09-14 07:15
rojana
Offline
Last seen: 1 month 1 week ago
Joined: 2017-08-11 05:43

quite eventful

Year 2013 has been quite eventful. SRSG also witnessed events and evolved in terms of business verticals, technologies, people and processes.
Apple service center Delhi
Apple iphone 7 reseller in Kolkata
ipad service center in Guwahati
Mac pro service center in Kolkata
FM Channel License

Thu, 2017-09-21 05:55
mb22
Offline
Last seen: 2 months 2 weeks ago
Joined: 2017-09-21 05:17

VIGYAN DHARA is the first

VIGYAN DHARA is the first Indian IIT-JEE/ NEET/AIIMS coaching Institute which is purely based on Historical Indian Education System that is the Gurukul Tradition of Education where students leave their homes during their studies.
PMT Coaching center in hisar
IIT coaching institute in hisar

Thu, 2017-09-21 05:55
mb22
Offline
Last seen: 2 months 2 weeks ago
Joined: 2017-09-21 05:17

The Fairdeal Group of

The Fairdeal Group of Companies is a fast-growing group with the ambition of helping people across Kenya build their dream homes and offices.
Furniture manufacture in Kenya
Dealers in metal beds Nairobi
Dealers in Kids Furniture in Kenya
Dealers in Coffee Tables in Kenya

Tue, 2017-09-26 07:54
mb22
Offline
Last seen: 2 months 2 weeks ago
Joined: 2017-09-21 05:17

Solutions

Solutions Middle East, works to Envision, Execute and Empower our clients by offerings design, implementation and on-going support services.
Citrix Mobility Solutions
FOIP Connection
SMS Gateway providers
Riverbed Application
Printing software solutions
Staff Outsourcing Companies

Tue, 2017-09-26 07:54
mb22
Offline
Last seen: 2 months 2 weeks ago
Joined: 2017-09-21 05:17

school

DPS are committed to provide our students with a sound knowledge base with emphasis on building strong fundamentals in the principles of education, moral values and strong character and concrete the same by the way of continuous training with psychological methods.
Best CBSE school in Rajasthan

Thu, 2017-09-28 03:12
sheeraz Khatri
Offline
Last seen: 2 days 6 hours ago
Joined: 2017-09-28 03:11

Hero Instinct

This is a fantastic website and I can not recommend you guys enough. 

Hero Instinct

Thu, 2017-09-28 13:04
sheeraz Khatri
Offline
Last seen: 2 days 6 hours ago
Joined: 2017-09-28 03:11

Make Him Desire You

These thoughts just blew my mind. I am glad you have posted this. 

Make Him Desire You

Fri, 2017-09-29 17:17
sheeraz Khatri
Offline
Last seen: 2 days 6 hours ago
Joined: 2017-09-28 03:11

Health Benefits of Coffee

Thanks a lot for one’s intriguing write-up. It’s actually exceptional. Searching ahead for this sort of revisions. 

Health Benefits of Coffee

Sun, 2017-10-29 05:03
sheeraz Khatri
Offline
Last seen: 2 days 6 hours ago
Joined: 2017-09-28 03:11

Brazilian hair wigs on sale

A great website with interesting and unique material what else would you need. 

Brazilian hair wigs on sale

Mon, 2017-11-06 05:47
rojana
Offline
Last seen: 1 month 1 week ago
Joined: 2017-08-11 05:43

Moonlight

Moonlight Films and Theatre Studio operates under the Moonlight Films and Theatre Society , a non-profit association whose goals are to spread knowledge of and interest in theatrical and media arts through education.
Short filmmaking courses in Delhi

Wed, 2017-11-08 12:18
sheeraz Khatri
Offline
Last seen: 2 days 6 hours ago
Joined: 2017-09-28 03:11

คาสิโนออนไลน์

It proved to be Very helpful to me and I am sure to all the commentators here!  

คาสิโนออนไลน์

Fri, 2017-11-10 14:11
sheeraz Khatri
Offline
Last seen: 2 days 6 hours ago
Joined: 2017-09-28 03:11

myfinancelab

Here at this site really the fastidious material collection so that everybody can enjoy a lot. 

myfinancelab

Tue, 2017-12-05 10:50
connexion Republic
Offline
Last seen: 2 days 8 hours ago
Joined: 2017-12-05 10:47

connexion Republic

Wow what a Great Information about World Day its exceptionally pleasant educational post. a debt of gratitude is in order for the post. 

connexion Republic

Thu, 2017-12-07 03:11
sheeraz Khatri
Offline
Last seen: 2 days 6 hours ago
Joined: 2017-09-28 03:11

prix carte grise

Thanks for providing recent updates regarding the concern, I look forward to read more.  

prix carte grise

Tue, 2017-12-12 08:51
connexion Republic
Offline
Last seen: 2 days 8 hours ago
Joined: 2017-12-05 10:47

Dofollow

If you want a great experience, this is the service to buy. A+ in communication, quality of service, and delivery. I took my time to analyse the links, they are quite rare and low OBL. So if you want an actual page rank, to boost your page and domain authority try and get this link. I will definitly buy this service again.  

Dofollow

Download PSpice Lite and try it for free! Get PSpice Lite
Cadence