410.bwaves
Dr. Mark Kremenetsky, <mdk [at] sgi.com>
Silicon Graphics, Inc
1500 Crittenden Lane
Mountain View, CA 94043, USA
Computational Fluid Dynamics
410.bwaves numerically simulates blast waves in three dimensional transonic transient laminar viscous flow.
The initial configuration of the blast waves problem consists of a high pressure and density region at the center of a cubic cell of a periodic lattice, with low pressure and density elsewhere. Periodic boundary conditions are applied to the array of cubic cells forming an infinite network. Initially, the high pressure volume begins to expand in the radial direction as classical shock waves. At the same time, the expansion waves move to fill the void at the center of the cubic cell. When the expanding flow reaches the boundaries, it collides with its periodic images from other cells, thus creating a complex structure of interfering nonlinear waves. These processes create a nonlinear damped periodic system with energy being dissipated in time. Finally, the system will come to an equilibrium and steady state.
The algorithm implemented is an unfactored solver for the implicit solution of the compressible Navier-Stokes equations using the Bi-CGstab algorithm, which solves systems of non-symmetric linear equations iteratively.
The input file describes the grid size, flow parameters, initial boundary condition and number of time steps. The three data sets, test, train and ref, differ only in grid size and number of time steps.
The transient nature of the flow and iterative solver makes bwaves a difficult problem to validate. In SPEC CPU2006 this has been addressed by comparing three different outputs. These are:
Fortran 77
none
Last updated: 27 June 2006