9722

Absolute Stability of an Integration Method

Consider the differential equation , where , and the integration scheme .
If , one recovers the explicit Euler integration scheme.
If , one recovers the implicit Euler integration scheme.
If , one recovers the Crank–Nicholson integration scheme.
Absolute stability means that the global error does not grow without bound.
Absolute stability of the numerical integration method requires the following condition:
, where the dimensionless time step is .
Consider the Euler integration scheme (i.e., ); the global error, , satisfies , where is the local truncation error at time step . This equation reveals the origin of exponential growth in the global error. Indeed, at each time step the global error is multiplied by . If is outside the disk with center at of radius 1, the global error will grow out of bounds and the numerical method will fail to give an accurate solution of the differential equation. Thus, the choice of the time step that gives good solutions is restricted to a specific domain and the method is not absolutely stable.
The present Demonstration shows the region plot in the plane for user-set values of . With , the light blue region corresponds to and the blue curve is where . Thus the light blue region represents choices of that give absolute stability.
Ifis real, one can conclude that the explicit Euler method is not absolutely stable, while both the implicit Euler and Crank–Nicholson methods are absolutely stable.

SNAPSHOTS

  • [Snapshot]
  • [Snapshot]
  • [Snapshot]

DETAILS

K. J. Beers, Numerical Methods for Chemical Engineering, New York: Cambridge University Press, 2007.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.









 
RELATED RESOURCES
Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+