1,218 research outputs found
Numerical Solution of ODEs and the Columbus' Egg: Three Simple Ideas for Three Difficult Problems
On computers, discrete problems are solved instead of continuous ones. One
must be sure that the solutions of the former problems, obtained in real time
(i.e., when the stepsize h is not infinitesimal) are good approximations of the
solutions of the latter ones. However, since the discrete world is much richer
than the continuous one (the latter being a limit case of the former), the
classical definitions and techniques, devised to analyze the behaviors of
continuous problems, are often insufficient to handle the discrete case, and
new specific tools are needed. Often, the insistence in following a path
already traced in the continuous setting, has caused waste of time and efforts,
whereas new specific tools have solved the problems both more easily and
elegantly. In this paper we survey three of the main difficulties encountered
in the numerical solutions of ODEs, along with the novel solutions proposed.Comment: 25 pages, 4 figures (typos fixed
Implicit-Explicit multistep methods for hyperbolic systems with multiscale relaxation
We consider the development of high order space and time numerical methods
based on Implicit-Explicit (IMEX) multistep time integrators for hyperbolic
systems with relaxation. More specifically, we consider hyperbolic balance laws
in which the convection and the source term may have very different time and
space scales. As a consequence the nature of the asymptotic limit changes
completely, passing from a hyperbolic to a parabolic system. From the
computational point of view, standard numerical methods designed for the
fluid-dynamic scaling of hyperbolic systems with relaxation present several
drawbacks and typically lose efficiency in describing the parabolic limit
regime. In this work, in the context of Implicit-Explicit linear multistep
methods we construct high order space-time discretizations which are able to
handle all the different scales and to capture the correct asymptotic behavior,
independently from its nature, without time step restrictions imposed by the
fast scales. Several numerical examples confirm the theoretical analysis
STICAP: A linear circuit analysis program with stiff systems capability. Volume 1: Theory manual
STICAP (Stiff Circuit Analysis Program) is a FORTRAN 4 computer program written for the CDC-6400-6600 computer series and SCOPE 3.0 operating system. It provides the circuit analyst a tool for automatically computing the transient responses and frequency responses of large linear time invariant networks, both stiff and nonstiff (algorithms and numerical integration techniques are described). The circuit description and user's program input language is engineer-oriented, making simple the task of using the program. Engineering theories underlying STICAP are examined. A user's manual is included which explains user interaction with the program and gives results of typical circuit design applications. Also, the program structure from a systems programmer's viewpoint is depicted and flow charts and other software documentation are given
Numerical integration and other techniques for computer aided network design programming Final technical report, 1 Jan. 1970 - 1 Jan. 1971
Matrix method and stiffly stable algorithms in numerical integration for computer aided network design programmin
The numerical solution of stiff differential equations
AbstractThis paper first discusses the conditions in which a set of differential equations should give stable solutions, starting with linear systems assuming that these do not differ greatly in this respect from non-linear systems. Methods of investigating the stability of particular systems are briefly discussed. Most real biochemical systems are known from observation to be stable, but little is known of the regions over which stability persists; moreover, models of biochemical systems may not be stable, because of inaccurate choice of parameter values.The separate problem of stability and accuracy in numerical methods of approximating the solution of systems of non-linear equations is then treated. Stress is laid on the consistently unsatisfactory results given by explicit methods for systems containing “stiff” equations, and implicit multistep methods are particularly recommended for this class of problem, which is likely to include many biochemical model systems. Finally, an iteration procedure likely to give convergence both in multistep methods and in the steady-state approach is recommended, and areas in which improvement in methods is likely to occur are outlined
- …