A Validation Method for Solution of Nonlinear Differential Equations: Construction of Exact Solutions Neighboring Approximate Solutions

Title & Authors
A Validation Method for Solution of Nonlinear Differential Equations: Construction of Exact Solutions Neighboring Approximate Solutions
Lee, Sang-Chul;

Abstract
An inverse method is introduced to construct benchmark problems for the numerical solution of initial value problems. Benchmark problems constructed through this method have a known exact solution, even though analytical solutions are generally not obtainable. The solution is constructed such that it lies near a given approximate numerical solution, and therefore the special case solution can be generated in a versatile and physically meaningful fashion and can serve as a benchmark problem to validate approximate solution methods. A smooth interpolation of the approximate solution is forced to exactly satisfy the differential equation by analytically deriving a small forcing function to absorb all of the errors in the interpolated approximate solution. A multi-variable orthogonal function expansion method and computer symbol manipulation are successfully used for this process. Using this special case exact solution, it is possible to directly investigate the relationship between global errors of a candidate numerical solution process and the associated tuning parameters for a given code and a given problem. Under the assumption that the original differential equation is well-posed with respect to the small perturbations, we thereby obtain valuable information about the optimal choice of the tuning parameters and the achievable accuracy of the numerical solution. Illustrative examples show the utility of this method not only for the ordinary differential equations (ODEs) but for the partial differential equations (PDEs).
Keywords
inverse method;benchmark problem;numerical solution;tuning parameters;
Language
English
Cited by
References
1.
Suk, J, and Kim, Y., "On the Modeling of Dynamic Systems," KSAS International Journal, Vol. 2, No.1, 2001, pp. 78-92.

2.
Miele, A., and Iyer, R R, "General Technique for Solving Nonlinear, Two-Point Boundary-Value Problems via the Method of Particular Solutions," Journal of Optimization Theory and Applications, Vol. 5, No. 5, 1970, pp. 382-403.

3.
Silver, W. M., "On the Equivalence of Lagrangian and Newton-Euler Dynamics for Manipulators," The International Journal of Robotics Research, Vol. 1, No. 2, 1982, pp. 60-70.

4.
Luh, J. Y. S., Walker, M. W., and Paul, R P. c, "On-Line Computational Scheme for Mechanical Manipulators," ASME Journal of Dynamic Systems, Measurement, and Control, Vol. 102, 1980, pp. 69-76.

5.
Gear, C. W., Numerical Initial Value Problems in Ordinary Differential Equations, Prentice-Hall, Englewood Cliffs, NJ, 1971.

6.
Shampine, L. F., "Limiting Precision in Differential Equation Solvers, II: Sources of Trouble and Starting a Code," Math Comp., Vol. 32, No. 144, 1978, pp. 1115-1122.

7.
Enright, W. H., "Analysis of Error Control Strategies for Continuous Runge-Kutta Methods," SIAM J. Numer. Anal., Vol. 26, No.3, 1989, pp. 588-599.

8.
Hull, T. E., Enright, W. H., Fellen, B. M., and Sedgwick, A. E., "Comparing Numerical Methods for Ordinary Differential Equations," SIAM J. Numer. Anal., Vol. 9, No.4, 1972, pp. 603-637.

9.
Krogh, F. T., "On Testing a Subroutine for the Numerical Integration of Ordinary Differential Equations," Journal of the Association for Computing Machinery, Vol. 20, No.4, 1973, pp, 545-562.

10.
IMSL MATH, LIBRARY User's' Manual Version 1.1, IMSL Inc., 1989.

11.
Junkins, J. L., An Introduction to Optimal Estimation of Dynamical Systems, Sijhoff & Noordhoff, Alphenaan den Rijn, The Netherlands, 1978.

12.
Abramowitz, M., and Stegun, I. A., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, National Bureau of Standards, Applied Mathematics Series 55, U.S. Department of Commerce, 1972.

13.
MA CSYMA Reference Manual Version 13, Symbolics Inc., 1988.

14.
Junkins, J. L., and Lee, S., "Benchmark Problems for the Solution of Ordinary Differential Equations," Shock and Vibration Computer Programs, Vol SVM-13, 1995, pp. 497-523.

15.
Battin, R. H., An Introduction to the Mathematics and Methods of Astrodynamics, AIAA Education Series, New York, New York, 1987.

16.
Lee, S., and Junkin s, J. L., "Explicit Generalization of Lagrange' s Equations for Hybrid Coordinate Dynamical Systems," AIAA Journal of Guidance, Control, and Dynamics, Vol. 15, No.6, 1992, pp. 1443-1452.

17.
Junkins, J. L., and Kim, Y., An Introduction to Dynamics and Control of Flexible Structures, AIAA, Washington, DC, 1993.

18.
Shampine, L. F., and Gordon, M. K., Computer solution of Ordinary Differential Equations, W. H. Freeman, San Francisco, 1975.

19.
Junkins, J. L., Miller, G. W., and Jancaiti s, J. R, "A Weighting Function Approach to Modeling of Irregular Surfaces," Journal of Geophysical Research, Vol. 78, No. 11, 1973, pp. 1794-1803.

20.
Enright, W. H., Hull, T. E., and Lindberg, B., "Comparing Numerical Methods for Stiff Systems of ODEs," BIT, Vol. 15, 1975, pp. 10-48.