Numerically Tracking Contact Discontinuities with an Introduction for GPU Programming

[摘要] We review some of the classic numerical techniques used to analyze contact discontinuities and compare their effectiveness. Several finite difference methods (the Lax-Wendroff method, a Multidimensional Positive Definite Advection Transport Algorithm (MPDATA) method and a Monotone Upstream Scheme for Conservation Laws (MUSCL) scheme with an Artificial Compression Method (ACM)) as well as the finite element Streamlined Upwind Petrov-Galerkin (SUPG) method were considered. These methods were applied to solve the 2D advection equation. Based on our results we concluded that the MUSCL scheme produces the sharpest interfaces but can inappropriately steepen the solution. The SUPG method seems to represent a good balance between stability and interface sharpness without any inappropriate steepening. However, for solutions with discontinuities, the MUSCL scheme is superior. In addition, a preliminary implementation in a GPU program is discussed.

[发布日期] 2012-08-17

[学科分类] 数学（综合）

[关键词] ADVECTION;ALGORITHMS;COMPRESSION;CONSERVATION LAWS;FINITE DIFFERENCE METHOD;IMPLEMENTATION;PROGRAMMING;STABILITY;TRANSPORT