Simulation of Contaminant Dispersion on the Cray X1: Verification and Implementation

Stabilized finite element formulation developed for simulation of dispersion of contaminants is implemented on the Cray X1. The stabilization is based on the SUPG and PSPG techniques. The governing equations are the incompressible Navier-Stokes equations coupled with the heat and mass transfer equations. The Boussinesq approximation in the momentum equation accounts for the density change due to thermal expansion. Fully implicit nonlinear systems of equations are solved iteratively using the matrix-free GMRES technique. The stabilized finite element formulation is parallelized and vectorized on the Cray X1. The three-dimensional validation problem involves transient simulation of flow past a building with source point releasing trances. Two-dimensional problems are simulated to compare the numerical results with analytical solutions.

[1]  Jalal Abedi,et al.  Simulation technique for wave generation , 2003 .

[2]  Tayfun E. Tezduyar,et al.  Parallel fluid dynamics computations in aerospace applications , 1995 .

[3]  Graham F. Carey,et al.  Penalty approximation of stokes flow , 1982 .

[4]  Tayfun E. Tezduyar,et al.  Simulation of multiple spheres falling in a liquid-filled tube , 1996 .

[5]  Francis A Kulacki,et al.  Steady and transient thermal convection in a fluid layer with uniform volumetric energy sources , 1977, Journal of Fluid Mechanics.

[6]  Thomas J. R. Hughes,et al.  A globally convergent matrix-free algorithm for implicit time-marching schemes arising in finite element analysis in fluids , 1991 .

[7]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[8]  J. Reisner,et al.  Simulations of flow around a cubical building: comparison with towing-tank data and assessment of radiatively induced thermal effects , 2001 .

[9]  S. Murakami,et al.  Three-dimensional numerical simulation of air flow around a cubic model by means of large eddy simulation , 1987 .

[10]  Marek Behr,et al.  Enhanced-Discretization Interface-Capturing Technique (EDICT) for computation of unsteady flows with interfaces , 1998 .

[11]  C. Kato,et al.  Large eddy simulation of unsteady turbulent wake of a circular cylinder using the finite element method , 1991 .

[12]  Jalal Abedi,et al.  Parallel finite element simulation of mooring forces on floating objects , 2003 .

[13]  Jalal Abedi,et al.  Stabilized finite element formulation of buoyancy driven incompressible flows , 2002 .

[14]  Tayfun E. Tezduyar,et al.  Stabilized-finite-element/interface-capturing technique for parallel computation of unsteady flows with interfaces , 2000 .

[15]  William H. Snyder,et al.  A Comparison of Numerical and Physical Modeling of Stable Atmospheric Flow and Dispersion around a Cubical Building , 1996 .

[16]  Tayfun E. Tezduyar,et al.  Flow simulation and high performance computing , 1996 .

[17]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[18]  Richard J Goldstein,et al.  Thermal convection in a horizontal fluid layer with uniform volumetric energy sources , 1972, Journal of Fluid Mechanics.

[19]  David E. Stock,et al.  The Numerical Simulation of Airflow and Dispersion in Three-Dimensional Atmospheric Recirculation Zones , 1991 .