Runge-Kutta formulas are given which are suited to the tasks arising in simulation. They are methods permitting interpolation which use overlap into the succeeding step to reduce the cost of a step and its error estimate.
Solving stiff differential equations for simulation
โ Scribed by T.D. Bui
- Publisher
- Elsevier Science
- Year
- 1981
- Tongue
- English
- Weight
- 540 KB
- Volume
- 23
- Category
- Article
- ISSN
- 0378-4754
No coin nor oath required. For personal study only.
โฆ Synopsis
Computer simulation
of dynamic systems very often leads to the solution of a set of stiff ordinary differential equations.
The solution of this set of equations involves the eigenvalues of its Jacobian matrix.
The greater the spread in eigenvalues, the more time consuming the solutions become when existing numerical methods are employed.
Extremely stiff differential equations can become a very serious problem for some systems, rendering accurate numerical solutions completely uneconomic.
In this paper, we propose new techniques for solving extremely stiff systems of differential equations.
These algorithms are based on a class of implicit Runge-Kutta procedure with complete error estimate.
The new techniques are applied to solving mathematical models of the relaxation problem behind blast waves.
๐ SIMILAR VOLUMES
Based on the idea of quasi-interpolation and radial basis functions approximation, a numerical method is developed to quasi-interpolate the forcing term of di erential equations by using radial basis functions. A highly accurate approximation for the solution can then be obtained by solving the corr
When the s-stage fully implicit Runge}Kutta (RK) method is used to solve a system of n ordinary di!erential equations (ODE) the resulting algebraic system has a dimension ns. Its solution by Gauss elimination is expensive and requires 2sn/3 operations. In this paper we present an e$cient algorithm,
## Communicated by E. Y. Rodin AImraet--A study of Rosenbrock-Wanner (ROW) methods showed that they are not AN-stable. A second-order improved ROW-method, which is AN-stable and with local error estimate is presented with some numerical results.
We introduce a new method for solving very stiff sets of ordinary differential equations. The basic idea is to replace the original nonlinear equations with a set of equally stiff equations that are piecewise linear, and therefore can be solved exactly. We demonstrate the value of the method on smal