TY - JOUR

T1 - Accelerated composite distribution function methods for computational fluid dynamics using GPU

AU - Smith, Matt-Hew

AU - Chen, Yen Chih

PY - 2013/1/1

Y1 - 2013/1/1

N2 - The Kinetic Theory of Gases has long been established as a useful tool for the solution of modern Computational Fluid Dynamics (CFD) problems. Together with the Finite Volume Method, such approaches have been popular in CFD for over 30 years, with techniques such as the Equilibrium Flux Method (EFM) or Kinetic Flux Vector Splitting (KFVS), Equilibrium Interface Method (EIM) together with more recent developments. One of the disadvantages to using such an approach are the expensive exponential (exp(-x2)) and error function (erf(x)) evaluations often associated with the moments taken around the distribution functions for the computation of interface fluxes. One common approach for avoiding such expenses is to employ discrete velocities in the flux calculation, taking moments around these rather than a continuous distribution function. In this talk we will discuss how we can approximate the governing particle velocity distribution function with a series of Composite Distribution Functions (CDF's) - made of more than one distribution function - to simplify the moment equations. The resulting expressions are then applied to multi-dimensional computation using Graphics Processing Units (GPU's), to which the application is well suited due to the simplicity of the flux expressions and locality of the schemes. Very high levels of speedup are demonstrated using C2075 (Fermi) and newer Kepler GPU architectures when compared to modern Xeon E5 processing cores.

AB - The Kinetic Theory of Gases has long been established as a useful tool for the solution of modern Computational Fluid Dynamics (CFD) problems. Together with the Finite Volume Method, such approaches have been popular in CFD for over 30 years, with techniques such as the Equilibrium Flux Method (EFM) or Kinetic Flux Vector Splitting (KFVS), Equilibrium Interface Method (EIM) together with more recent developments. One of the disadvantages to using such an approach are the expensive exponential (exp(-x2)) and error function (erf(x)) evaluations often associated with the moments taken around the distribution functions for the computation of interface fluxes. One common approach for avoiding such expenses is to employ discrete velocities in the flux calculation, taking moments around these rather than a continuous distribution function. In this talk we will discuss how we can approximate the governing particle velocity distribution function with a series of Composite Distribution Functions (CDF's) - made of more than one distribution function - to simplify the moment equations. The resulting expressions are then applied to multi-dimensional computation using Graphics Processing Units (GPU's), to which the application is well suited due to the simplicity of the flux expressions and locality of the schemes. Very high levels of speedup are demonstrated using C2075 (Fermi) and newer Kepler GPU architectures when compared to modern Xeon E5 processing cores.

UR - http://www.scopus.com/inward/record.url?scp=84976260249&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84976260249&partnerID=8YFLogxK

M3 - Conference article

AN - SCOPUS:84976260249

VL - 14-17-October-2013

JO - Proceedings of Science

JF - Proceedings of Science

SN - 1824-8039

M1 - 068

T2 - 1st International Workshop on Computational Science and Engineering, IWCSE 2013

Y2 - 14 October 2013 through 17 October 2013

ER -