ISI TAŞINIMI İÇİN YAPISAL OLMAYAN AĞLARDA BİR HIZLANDIRILMIŞ SÜREKSİZ GALERKİN METDU: FORMÜLASYONU VE DOĞRULANMASI

Öz

Bu çalışmada, büyük ölçekli, birleştirilmiş sıkıştırılamaz akış ve ısı transferi problemleri için GPU ile hızlandırılmış bir yöntem sunulmuştur. Yapılandırılmamış üçgen ağlar üzerinde tanımlayıcı denklemlerini ayrıklaştırmak için yüksek dereceli, nodal süreksiz Galerkin yöntemi kullanılmıştır. Zaman ayrıklaştırması, taşınım terimlerinin açık bir şekilde ele alındığı, Stokes operatörlerinin ise örtük olarak çözüldüğü bir yarı-örtük şema kullanılarak elde edilmiştir. Basınç sistemi, tamamen GPU ile hızlandırılmış bir çoklu-ağ ön koşullandırıcının kullanıldığı, eşlenik gradyan yöntemiyle çözülmüştür. Kod, yüksek dereceli ayrıklaştırmalar için yüksek performanslı kerneller sağlayan, ölçeklenebilir kütüphane olan libParanumal üzerine yazılmıştır. Platformdan bağımsız performans taşınabilirliği OCCA (open concurrent compute abstraction) dili ile elde edilmiştir. Yapılan serbest ve karışık taşınım problemlerini içeren bir dizi sayısal test ile sunulan metodun deneysel olarak beklenen, spektral doğruluğa ulaştığını gösterilmiştir.

Anahtar Kelimeler

• süreksiz Galerkin
• GPU
• paralel
• sıkıştırılamaz akış
• ısı transferi
• yüksek seviyeli hesaplama.

AN ACCELERATED NODAL DISCONTINUOUS GALERKIN METHOD FOR THERMAL CONVECTION ON UNSTRUCTURED MESHES: FORMULATION AND VALIDATION

Öz

We present a GPU-accelerated method for large scale, coupled incompressible fluid flow and heat transfer problems. A high-order, nodal discontinuous Galerkin method is utilized to discretize governing equations on unstructured triangular meshes. A semi-implicit scheme with explicit treatment of the advective terms and implicit treatment of the split Stokes operators are used for time discretization. The pressure system is solved with a conjugate gradient method together with a fully GPU-accelerated multigrid preconditioner. The code is built on scalable libParanumal solver which is a library of high-performance kernels for high-order discretizations. Performance portability is achieved by using the open concurrent compute abstraction, OCCA. A set of numerical experiments including free and mixed convection problems indicate that our approach experimentally reaches design order of accuracy.

Anahtar Kelimeler

• discontinuous Galerkin
• GPU
• parallel
• incompressible
• heat transfer
• high-order.

Kaynakça

1. Adams, M., Brezina, M., Hu, J., and Tuminaro, R., 2003, Parallel multigrid smoothing: polynomial versus Gauss–Seidel, Journal of Computational Physics, 188(2), 593–610.
2. Arnold, D., 1982, An interior penalty finite element method with discontinuous elements. SIAM Journal on Numerical Analysis, 19(4), 742–760.
3. Chalmers, N., Karakus, A., Austin, A. P., Swirydowicz, K., and Warburton, T., 2020, libParanumal: a performance portable high-order finite element library. Release 0.4.0.
4. Chan, J., Wang, Z., Modave, A., Remacle, J. F., and Warburton, T., 2016, GPU-accelerated discontinuous Galerkin methods on hybrid meshes, Journal of Computational Physics, 318, 142–168.
5. Darekar, R. M. and Sherwin, S. J., 2001, Flow past a square-section cylinder with a wavy stagnation face, Journal of Fluid Mechanics, 426, 263–295.
6. De Vahl Davis, G., 1983, Natural convection of air in a square cavity: A bench mark numerical solution, International Journal for Numerical Methods in Fluids, 3(3), 249–264.
7. Ferrer, E. and Willden, R. H. J., 2011, A high order discontinuous Galerkin finite element solver for the incompressible Navier–Stokes equations, Computers & Fluids, 46(1), 224–230.
8. Gandham, R., Esler, K., and Zhang, Y., 2014, A GPU accelerated aggregation algebraic multigrid method, Computers & Mathematics with Applications, 68(10), 1151–1160.

9. Gandham, R., Medina, D., and Warburton, T., 2015, GPU accelerated discontinuous Galerkin methods for shallow water equations, Communications in Computational Physics, 18(1), 37–64.
10. Hesthaven, J. S. and Warburton, T., 2008, Nodal discontinuous Galerkin methods: algorithms, analysis, and applications, Springer.
11. Hossain, M. Z., Cantwell, C. D., and Sherwin, S. J., 2021, A spectral/hp element method for thermal convection, International Journal for Numerical Methods in Fluids, 93(7), 2380–2395.
12. Karakus, A., Chalmers, N., Hesthaven, J. S., and Warburton, T., 2019a, Discontinuous Galerkin discretizations of the Boltzmann–BGK equations for nearly incompressible flows: Semi-analytic time stepping and absorbing boundary layers, Journal of Computational Physics, 390, 175–202.
13. Karakus, A., Chalmers, N., Swirydowicz, K., and Warburton, T., 2019b, A GPU accelerated discontinuous Galerkin incompressible flow solver, Journal of Computational Physics, 390, 380–404.
14. Karakus, A., Warburton, T., Aksel, M. H., and Sert, C., 2016a, A GPU-accelerated adaptive discontinuous Galerkin method for level set equation, International Journal of Computational Fluid Dynamics, 30(1), 56–68.
15. Karakus, A., Warburton, T., Aksel, M. H., and Sert, C., 2016b, A GPU accelerated level set reinitialization for an adaptive discontinuous Galerkin method, Computers & Mathematics with Applications, 72(3), 755–767.
16. Kumar, A. and Pothérat, A., 2020, Mixed baroclinic convection in a cavity, Journal of Fluid Mechanics, 88. Publisher: Cambridge University Press.
17. Medina, D. S., St-Cyr, A., and Warburton, T., 2014, OCCA: A unified approach to multi-threading languages. arXiv:1403.0968.
18. Modave, A., St-Cyr, A., and Warburton, T., 2016, GPU performance analysis of a nodal discontinuous Galerkin method for acoustic and elastic models, Computers & Geosciences, 91, 64–76.
19. Notay, Y., 2006, Aggregation-based algebraic multilevel preconditioning, SIAM journal on matrix analysis and applications, 27(4), 998–1018.
20. Notay, Y., 2010, An aggregation-based algebraic multigrid method, Electronic transactions on numerical analysis, 37(6), 123–146.
21. Roca, X., Nguyen, N. C., and Peraire, J., 2011, GPU-accelerated sparse matrix-vector product for a hybridizable discontinuous Galerkin method. In Aerospace Sciences Meetings. American Institute of Aeronautics and Astronautics, AIAA 2011–687.
22. Saha, S., Klewicki, J. C., Ooi, A. S. H., and Blackburn, H. M., 2015, Comparison of thermal scaling properties between turbulent pipe and channel flows via DNS, International Journal of Thermal Sciences, 89, 43–57.
23. Shahbazi, K., 2005, An explicit expression for the penalty parameter of the interior penalty method, Journal of Computational Physics, 205(2), 401–407.
24. Shahbazi, K., Fischer, P. F., and Ethier, C. R., 2007, A high-order discontinuous Galerkin method for the unsteady incompressible Navier-Stokes equations, Journal of Computational Physics, 222(1), 391– 407.
25. Stokos, K., Vrahliotis, S., Pappou, T., and Tsangaris, S., 2015, Development and validation of an incompressible Navier-Stokes solver including convective heat transfer, International Journal of Numerical Methods for Heat & Fluid Flow, 25(4), 861–886.
26. Warburton, T., 2006, An explicit construction of interpolation nodes on the simplex, Journal of Engineering Mathematics, 56(3), 247–262.
27. Wheeler, M. F., 1978, An elliptic collocation-finite element method with interior penalties, SIAM Journal on Numerical Analysis, 15(1), 152–161.
28. Swirydowicz, K., Chalmers, N., Karakus, A., and Warburton, T., 2019, Acceleration of tensor-product ´ operations for high-order finite element methods, The International Journal of High Performance Computing Applications, 33(4).