Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] acfk_v1_0.gz(63 Kbytes)
Manuscript Title: BOREAS: a program for 1-D ideal fluid dynamics with shocks.
Authors: D. Odstrcil
Program title: BOREAS
Catalogue identifier: ACFK_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 31(1984)13
Programming language: Fortran.
Computer: SIEMENS 4004/150.
Operating system: BS 1000 RESP, VM/370.
RAM: 200K words
Word size: 32
Keywords: Fluid dynamics, Gas, Fct algorithm, Shock waves, One-dimensional, Explicit Finite-difference.
Classification: 12.

Subprograms used:
Cat Id Title Reference
ABUF_v1_0 OLYMPUS (ICL4/70) CPC 7(1974)245
ABUF_v2_0 OLYMPUS (IBM 370/165) CPC 9(1975)51
ABUF_v3_0 OLYMPUS (CDC 6500) CPC 10(1975)167
ABUH_v1_0 TIMER CPC 8(1974)118
ABUI_v1_0 LEDGER CPC 8(1974)123

Nature of problem:
The purpose of BOREAS is to calculate ideal fluid dynamics with shocks. A one-dimensional hydrodynamic model based on the mass, momentum and total energy conservation laws is used. Optionally a gravitational force field may be taken into account and the equation of the total energy conservation may be replaced by the pressure equation.

Solution method:
The hdrodynamic equations are solved by the FCT algorithm, which is a two-step explicit finite-difference scheme. This algorithm has a second order amplitude and fourth order phase accuracy. Solutions are possible in either Cartesian, cylindrical or spherical geometry. The mesh may be arbitrarily spaced and may be Eulerian or Lagrangian.

Version 1 of BOREAS provides a simple model for the fluid dynamics with shocks. A number of physical phenomena and situations have been omitted. They may be included in later versions or by use of the EXPERT facility. The accuracy of the solution depends on the mesh size as well as on the choice of timestep.

Unusual features:
BOREAS is written in OLYMPUS FORTRAN, i.e. Standard FORTRAN (ASA X3.9- 1966) except for the use of the NAMELIST facility and AND and ERR parameters in the READ statement. COMMON blocks are given only once and before compilation C/INSERT statements must be replaced by the COMMON blocks to which they refer. Program BOREAS uses the OLYMPUS system environment which provides an underlying program structure, supplies system-dependent features and contains some utility subprograms and standard COMMON blocks. Some FORMAT statements assume four characters per word. The main program may require a program card on some computers.

Running time:
Execution times depend mainly on the mesh size which determines the number of meshpoints to be used and maximum value of the timestep. The standard Test 1 described in the long write-up is carried out on a mesh with 51 cells for 100 timesteps and requires 39 s on the IBM 370/148 at VLS1PK Bratislava.