Wavelet‐adaptive solvers on multi‐core architectures for the simulation of complex systems

We build wavelet‐based adaptive numerical methods for the simulation of advection‐dominated flows that develop multiple spatial scales, with an emphasis on fluid mechanics problems. Wavelet‐based adaptivity is inherently sequential and in this work we demonstrate that these numerical methods can be implemented in software that is capable of harnessing the capabilities of multi‐core architectures while maintaining their computational efficiency. Recent designs in frameworks for multi‐core software development allow us to rethink parallelism as task‐based, where parallel tasks are specified and automatically mapped onto physical threads. This way of exposing parallelism enables the parallelization of algorithms that were considered inherently sequential, such as wavelet‐based adaptive simulations. In this paper we present a framework that combines wavelet‐based adaptivity with the task‐based parallelism. We demonstrate the promising performance obtained by simulating various physical systems on different multi‐core architectures using up to 16 cores. Copyright © 2010 John Wiley & Sons, Ltd.

[1]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[2]  Alessandro Curioni,et al.  Billion vortex particle direct numerical simulations of aircraft wakes , 2008 .

[3]  Margaret Martonosi,et al.  Characterizing and improving the performance of Intel Threading Building Blocks , 2008, 2008 IEEE International Symposium on Workload Characterization.

[4]  Thomas Ertl,et al.  Hardware Accelerated Wavelet Transformations , 2000, VisSym.

[5]  I. Daubechies,et al.  Biorthogonal bases of compactly supported wavelets , 1992 .

[6]  Oleg V. Vasilyev,et al.  A Brinkman penalization method for compressible flows in complex geometries , 2007, J. Comput. Phys..

[7]  Olivier Roussel,et al.  A conservative fully adaptive multiresolution algorithm for parabolic PDEs , 2003 .

[8]  Oleg V. Vasilyev,et al.  An Adaptive Wavelet Collocation Method for Fluid-Structure Interaction at High Reynolds Numbers , 2005, SIAM J. Sci. Comput..

[9]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[10]  John E. Pearson,et al.  Self-replicating spots in reaction-diffusion systems , 1997 .

[11]  Y. Kaneda,et al.  Study of High-Reynolds Number Isotropic Turbulence by Direct Numerical Simulation , 2009 .

[12]  Francisco Tirado,et al.  Parallel Implementation of the 2D Discrete Wavelet Transform on Graphics Processing Units: Filter Bank versus Lifting , 2008, IEEE Transactions on Parallel and Distributed Systems.

[13]  David L. Donoho,et al.  Interpolating Wavelet Transforms , 1992 .

[14]  Linda Yang,et al.  Coarse-Grained Parallel Algorithms for Multi-Dimensional Wavelet Transforms , 2004, The Journal of Supercomputing.

[15]  David A. Bader,et al.  Computing discrete transforms on the Cell Broadband Engine , 2009, Parallel Comput..

[16]  Kai Schneider,et al.  An adaptive multiresolution scheme with local time stepping for evolutionary PDEs , 2008, J. Comput. Phys..

[17]  A. Harten.  Adaptive Multiresolution Schemes for Shock Computations , 1994 .

[18]  Oleg V. Vasilyev,et al.  Simultaneous space-time adaptive wavelet solution of nonlinear parabolic differential equations , 2006, J. Comput. Phys..

[19]  Kai Schneider,et al.  Space--time adaptive multiresolution methods for hyperbolic conservation laws: Applications to compressible Euler equations , 2009 .

[20]  Michael Voss,et al.  Optimization via Reflection on Work Stealing in TBB , 2008, 2008 IEEE International Symposium on Parallel and Distributed Processing.

[21]  Jacques Liandrat,et al.  Resolution of the 1D regularized Burgers equation using a spatial wavelet approximation , 1990 .

[22]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[23]  B. Wendroff,et al.  Approximate Riemann Solvers, Godunov Schemes and Contact Discontinuities , 2001 .

[24]  Javier Jiménez,et al.  Computing high-Reynolds-number turbulence: will simulations ever replace experiments?* , 2003 .

[25]  P. Koumoutsakos.  MULTISCALE FLOW SIMULATIONS USING PARTICLES , 2005 .

[26]  Michael Bergdorf,et al.  A Lagrangian Particle-Wavelet Method , 2006, Multiscale Model. Simul..