Algebraic Stabilization of Explicit Numerical Integration for Extremely Stiff Reaction Networks

Details Algebraic Stabilization of Explicit Numerical Integration for Extremely Stiff Reaction Networks Algebraic Stabilization of Explicit Numerical Integration for Extremely Stiff Reaction Networks

2011-12-20

arxiv.org/abs/1112.4778v2

Astronomy and Astrophysics

Astrophysics

Solar and Stellar Astrophysics

Updated references

Scientific paper

In contrast to the prevailing view in the literature, it is shown that even extremely stiff sets of ordinary differential equations may be solved efficiently by explicit methods if limiting algebraic solutions are used to stabilize the numerical integration. The stabilizing algebra differs essentially for systems well-removed from equilibrium and those near equilibrium. Explicit asymptotic and quasi-steady-state methods that are appropriate when the system is only weakly equilibrated are examined first. These methods are then extended to the case of close approach to equilibrium through a new implementation of partial equilibrium approximations. Using stringent tests with astrophysical thermonuclear networks, evidence is provided that these methods can deal with the stiffest networks, even in the approach to equilibrium, with accuracy and integration timestepping comparable to that of implicit methods. Because explicit methods can execute a timestep faster and scale more favorably with network size than implicit algorithms, our results suggest that algebraically-stabilized explicit methods might enable integration of larger reaction networks coupled to fluid dynamics than has been feasible previously for a variety of disciplines.

Affiliated with

Also associated with

No associations

Rating

Algebraic Stabilization of Explicit Numerical Integration for Extremely Stiff Reaction Networks does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Algebraic Stabilization of Explicit Numerical Integration for Extremely Stiff Reaction Networks, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Algebraic Stabilization of Explicit Numerical Integration for Extremely Stiff Reaction Networks will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-58732