Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

FORTRAN oriented simulation system for the general solution of partial differential equations

Technical Report ·
OSTI ID:4936628
The traditional continuous simulation system is merely a program package which handles the integration of simultaneous ordinary differential equations, using stardard integration algorithms. Such a program may be used to integrate partial differential equations only if the user first discretizes the equations in space, thus transforming them into an ordinary differential equation system. This is frequently an onerous task. A program designed to handle the solution of partial differential equations in a generalized manner is described. The user supplies the partial differential equations in a FORTRAN sub-routine together with all the necessary boundary conditions, and the program completes the solution automatically, using the method of discrete ordinates. In this method, each partial differential equation is converted into a set of coupled ordinary differential equations in time which are then integrated by a standard algorithm. The user may specify the number of spatial divisions, and whether three-point or five-point difference approximations to the spatial derivatives are to be used. The paper concludes with a parametric investigation of the accuracy obtained from different methods of spatial division by comparing computed results to analytical solutions. (auth)
Research Organization:
Atomic Energy of Canada Ltd., Chalk River, Ontario. Chalk River Nuclear Labs.
Sponsoring Organization:
Sponsor not identified
NSA Number:
NSA-29-015087
OSTI ID:
4936628
Report Number(s):
AECL--4603; CONF-730704--3
Country of Publication:
Canada
Language:
English