Parallel algorithms for computational continuum dynamics (Q1085702)

This research was undertaken to develop parallel algorithms for explicit and implicit, Lagrangian and Eulerian finite difference schemes for computational continuum dynamics in one spatial dimension. The explicit conservation equations in the Lagrangian reference frame were readily reformulated for concurrent processing. An implicit solution was derived for these equations. The parallelism is achieved via a block implicit numerical scheme. A rezoning algorithm was employed with each Lagrangian integration step to transform the mesh back to the Eulerian reference frame. The algorithmic development path lead to a parallelization of the processing in blocks of the finite discretization zones. At each step of this research project, the derived numerical methods provided effective algorithms for exploiting the architectural advantages of a MIMD machine such as the HEP H1000 (heterogeneous element processor) computer. The computational timing data show significant speed-up with the number of processes.