A GPU solver for symmetric positive-definite matrices vs. traditional codes

journal article in Web of Science

en

In Heat Transfer and Fluid Flow Laboratory in Brno, the inverse heat conduction problem (IHCP) has been extensively used to reconstruct thermal boundary conditions at hot surfaces of solid materials cooled by spraying nozzles. More than three decades of experience and cooperation with industries has proven our experimental/numerical technique to be reliable and very accurate. However, a typical calculation requires relatively long calculation time. The transient heat diffusion in a multi-material sample is the most computationally costly ingredient of the algorithm. In the present paper, the potential for speeding up our calculations is manifested by firstly benchmarking it against traditional CFD codes such as OpenFOAM (FDIC) and ANSYS Fluent (AMG). Secondly, we also unveil a unique comparison between the performance of three inhouse GPU codes each written by a different PhD student/postdoc. Chronologically listed, one student pushed his luck with a fully explicit scheme, while the other two, including us, bet on implicit methods namely the line-by-line method in OpenCL and the conjugate gradient method with the deflated truncated Neumann series preconditioner in CUDA C. (C) 2019 The Authors. Published by Elsevier Ltd.

In Heat Transfer and Fluid Flow Laboratory in Brno, the inverse heat conduction problem (IHCP) has been extensively used to reconstruct thermal boundary conditions at hot surfaces of solid materials cooled by spraying nozzles. More than three decades of experience and cooperation with industries has proven our experimental/numerical technique to be reliable and very accurate. However, a typical calculation requires relatively long calculation time. The transient heat diffusion in a multi-material sample is the most computationally costly ingredient of the algorithm. In the present paper, the potential for speeding up our calculations is manifested by firstly benchmarking it against traditional CFD codes such as OpenFOAM (FDIC) and ANSYS Fluent (AMG). Secondly, we also unveil a unique comparison between the performance of three inhouse GPU codes each written by a different PhD student/postdoc. Chronologically listed, one student pushed his luck with a fully explicit scheme, while the other two, including us, bet on implicit methods namely the line-by-line method in OpenCL and the conjugate gradient method with the deflated truncated Neumann series preconditioner in CUDA C. (C) 2019 The Authors. Published by Elsevier Ltd.

Linear solver; Inverse task; Heat transfer; GPU; CUDA; OpenFOAM

07.03.2019

PERGAMON-ELSEVIER SCIENCE LTD

OXFORD

0898-1221

78

9

2933–2943

11