DOI QR코드

DOI QR Code

고차 정확도 수치기법의 GPU 계산을 통한 효율적인 압축성 유동 해석

EFFICIENT COMPUTATION OF COMPRESSIBLE FLOW BY HIGHER-ORDER METHOD ACCELERATED USING GPU

  • 장태규 (서울대학교 기계항공공학부) ;
  • 박진석 (서울대학교 항공우주신기술연구소) ;
  • 김종암 (서울대학교 기계항공공학부)
  • Chang, T.K. (Department of Aerospace Engineering, Seoul National University) ;
  • Park, J.S. (Institute of Advanced Aerospace Technology, Seoul National University) ;
  • Kim, C. (Department of Aerospace Engineering, Seoul National University)
  • 투고 : 2014.06.25
  • 심사 : 2014.09.11
  • 발행 : 2014.09.30

초록

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

키워드

참고문헌

  1. 1998, Cockburn, B. and Shu, C.-W., "The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V: Multidimensional Systems," J. of Comp. Physics, Vol.141, pp.199-224. https://doi.org/10.1006/jcph.1998.5892
  2. 2007, Huynh, H.T., "A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods," 18th AIAA Computational Fluid Dynamics Conference, AIAA2007-4079.
  3. 2009, Wang, Z.J. and Gao, H., "A Unifying Lifting Collocation Penalty Formulation Including the Discontinuous Galerkin, Spectral Volume/Difference Methods for Conservation Laws on Mixed Grids," J. of Comp. Physics, Vol.228, pp.8161-8186. https://doi.org/10.1016/j.jcp.2009.07.036
  4. 2008, Yoon, S.-H., Kim, C. and Kim, K.-H., "Multi-dimensional Limiting process for Three Dimensional Flow Physics Analyses," J. Comp. Physics, Vol.227, pp.6001-6043. https://doi.org/10.1016/j.jcp.2008.02.012
  5. 2010, Park, J.S., Yoon, S.-H. and Kim, C., "Multi-dimensional Limiting Process for Hyperbolic Conservation Laws on Unstructured Grids," J. of Comp. Physics, Vol.229, pp.788-812. https://doi.org/10.1016/j.jcp.2009.10.011
  6. 2012, Park, J.S. and Kim, C., "Multi-dimensional limiting process for finite volume methods on unstructured grids," Computers and Fluids, Vol.65, pp.8-24. https://doi.org/10.1016/j.compfluid.2012.04.015
  7. 2014, Park, J.S. and Kim, C., "Higher-order Multi-dimensional Limiting Strategy for Discontinuous Galerkin Methods in Compressible Inviscid and Viscous Flows," Computers & Fluids, Vol.96, pp.377-396. https://doi.org/10.1016/j.compfluid.2013.11.030
  8. 2014, Park, J.S., Chang, T.K., and Kim, C., "Higher-order Multi-dimensional Limiting Strategy for Correction Procedure via Reconstruction," AIAA 52nd Aerospace Sciences Meeting, AIAA2014-0772.
  9. 2012, NVIDIA, CUDA C Programming Guide, Ver.5.0.
  10. 2009, Corrigan, A., Camelli, F. and Lohner, R., "Running Unstructured Grid Based CFD Solvers on Modern Graphics Hardware," 19th AIAA Computational Fluid Dynamics Conference, AIAA2009-4001.
  11. 2009, Kloockner, A., Warburton, T. and Hesthaven, J.S., "Nodal Discontinuous Galerkin Methods on Graphics Processors," J. of Comp. Physics, Vol.228, pp.7863-7882. https://doi.org/10.1016/j.jcp.2009.06.041
  12. 2013, Zimmerman, B.J. and Wang, Z.J., "The Efficient Implementation of Correction Procedure via Reconstruction with GPU Computing," 21st AIAA Computational Fluid Dynamics Conference, AIAA2013-2692.
  13. 1999, Karniadakis, G.E. and Sherwin, S.J., "Spectral/hp element methods for CFD," Oxford University Press, USA.
  14. 2003, Kim, S.-S., Kim, C., Rho, O.-H. and Hong, S.K., "Cures for the Shock Instability: Development of a Shock-stable Roe Scheme," J. of Comp. Physics, Vol.185, pp.342-374. https://doi.org/10.1016/S0021-9991(02)00037-2
  15. 2013, Wang, Z.J. and Gao, H., "A Conservative Correction Procedure via Reconstruction Formulation with the Chain-Rule Divergence Evaluation," J. of Comp. Physics, Vol.232, pp.7-13. https://doi.org/10.1016/j.jcp.2012.08.030
  16. 2001, Spiteri, R.J. and Ruuth, S.J., "A New Class of Optimal High-Order Strong Stability-Preserving Time Discretization Methods," SIAM J. on Num. Analysis, Vol.40, No.2, pp.469-491.
  17. 2000, Chang, S.-M. and Chang, K.-S., "On the Shock-Vortex Interaction in Schardin's Problem," Shock Waves, Vol.10, pp.333-343. https://doi.org/10.1007/s001930000061