Can a Positive Real Part of an Eigenvalue Indicate a Stiff ODE?

  • Thread starter Master J
  • Start date
  • Tags
    Ode
In summary, a stiff ODE is characterized by a large difference between the maximum and minimum eigenvalues of its Jacobian. The real part of the eigenvalue controls the error in a numerical approximation, with a negative real part leading to error decay and a positive real part leading to error growth. Therefore, a positive real part can also indicate a stiff ODE, even if the error does not decay but instead grows due to certain components dominating after a certain time. This behavior is still considered indicative of a stiff ODE.
  • #1
Master J
226
0
A stiff ODE is defined as one for which the magnitude of the maximum eigenvalue of its Jacobian is much greater than that of the mininmum.

It is the real part of the eigenvalue which controls the error in an approximation when a numerical scheme is used to solve the ODE. If it is negative, the error decays away and the approximation approaches the true value for higher iterations.

My question is, is a positive real part also indicitive of a stiff ODE (in the definiton above)? If the error doesn't decay but grows instead, and some components dominate after a certain time, is this still a stiff ODE?
 
Physics news on Phys.org
  • #2


Yes, a positive real part of the eigenvalue can also indicate a stiff ODE. In this case, the error in the numerical approximation will not decay but will instead grow. This can happen when certain components of the ODE dominate after a certain time, causing the error to increase rather than decrease. This behavior is still characteristic of a stiff ODE, as it is defined by the large difference between the maximum and minimum eigenvalues of the Jacobian.
 

FAQ: Can a Positive Real Part of an Eigenvalue Indicate a Stiff ODE?

What is a stiff ODE?

A stiff ODE is an ordinary differential equation that contains two types of solutions: fast and slow. Fast solutions change rapidly while slow solutions change slowly. This results in a large difference in the time scales of the solutions, making the equation "stiff."

How do you know if an ODE is stiff?

An ODE is considered stiff if it contains solutions that are very different in magnitude and require different numerical methods to accurately solve. This can be determined by looking at the stiffness ratio, which is the ratio of the largest to the smallest eigenvalues of the Jacobian matrix of the ODE.

Why is it important to understand stiff ODEs?

Stiff ODEs are commonly found in many scientific and engineering applications, such as chemical reactions, electrical circuits, and population dynamics. Understanding stiff ODEs is important in order to accurately model and simulate these systems, as well as to choose the appropriate numerical methods for solving them.

What are some methods for solving stiff ODEs?

Some methods commonly used for solving stiff ODEs include implicit methods, such as the backward Euler method and the trapezoidal method, as well as multi-step methods, such as the Adams-Bashforth and Adams-Moulton methods. These methods are designed to handle the stiffness of the equations and provide accurate solutions.

How can I improve the accuracy of my solutions for stiff ODEs?

One way to improve the accuracy of solutions for stiff ODEs is to use adaptive step size control. This means adjusting the step size of the numerical method based on the stiffness of the equation and the error in the solution. Another method is to use higher order methods, which can provide more accurate results but may also require more computational resources.

Similar threads

Replies
5
Views
1K
Replies
8
Views
3K
Replies
13
Views
2K
Replies
4
Views
4K
Replies
5
Views
4K
Replies
2
Views
8K
Back
Top