A framework for exploring numerical solutions of advection–reaction–diffusion equations using a GPU-based approach
Autor: | Allen Sanderson, Miriah Meyer, Chris R. Johnson, Robert M. Kirby |
---|---|
Rok vydání: | 2008 |
Předmět: |
Mathematical optimization
Partial differential equation Computation Numerical analysis General Engineering Graphics processing unit Theoretical Computer Science Euler equations Computational science symbols.namesake Computational Theory and Mathematics Modeling and Simulation symbols Euler's formula Crank–Nicolson method Computer Vision and Pattern Recognition Central processing unit Software ComputingMethodologies_COMPUTERGRAPHICS Mathematics |
Zdroj: | Computing and Visualization in Science. 12:155-170 |
ISSN: | 1433-0369 1432-9360 |
Popis: | In this paper we describe a general purpose, graphics processing unit (GP-GPU)-based approach for solving partial differential equations (PDEs) within advection–reaction–diffusion models. The GP-GPU-based approach provides a platform for solving PDEs in parallel and can thus significantly reduce solution times over traditional CPU implementations. This allows for a more efficient exploration of various advection–reaction–diffusion models, as well as, the parameters that govern them. Although the GPU does impose limitations on the size and accuracy of computations, the PDEs describing the advection–reaction–diffusion models of interest to us fit comfortably within these constraints. Furthermore, the GPU technology continues to rapidly increase in speed, memory, and precision, thus applying these techniques to larger systems should be possible in the future. We chose to solve the PDEs using two numerical approaches: for the diffusion, a first-order explicit forward Euler solution and a semi-implicit second order Crank–Nicholson solution; and, for the advection and reaction, a first-order explicit solution. The goal of this work is to provide motivation and guidance to the application scientist interested in exploring the use of the GP-GPU computational framework in the course of their research. In this paper, we present a rigorous comparison of our GPU-based advection–reaction–diffusion code model with a CPU-based analog, finding that the GPU model out-performs the CPU implementation in one-to-one comparisons. |
Databáze: | OpenAIRE |
Externí odkaz: |