351.bwaves
Dr. Mark Kremenetsky, <mdk [at] sgi.com>
Silicon Graphics International Corp.
46600 Landing Parkway
Fremont, CA 94538, USA
Computational Fluid Dynamics
351.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 biconjugate gradient stabilized (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 test case is a 3-D laminar flow with Reynolds Number equal to 1E+05 and Prandtl Number equals to 0.72. In the reference data set, there are 53*64*2048 computational cells in the 1st, 2nd and 3rd dimensions, respectively. The problem is solved using by implicit time-marching scheme with a CFL number equals to 2.
The transient nature of the flow and iterative solver makes bwaves a difficult problem to validate. In SPEC OMP2012 this has been addressed by comparing three different outputs. These are:
Fortran 77
none
Last updated: June 20, 2012