Programs in Physics & Physical Chemistry
|[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.|
|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|
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.
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.
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.
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.
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|