High-order Tensor-Train Finite Volume Method for Shallow Water Equations

In this paper, we introduce a high-order tensor-train (TT) finite volume method for the Shallow Water Equations (SWEs). We present the implementation of the \(3^{rd}\) order Upwind and the \(5^{th}\) order Upwind and WENO reconstruction schemes in the TT format. It is shown in detail that the linear...

Full description

Saved in:
Bibliographic Details
Published inarXiv.org
Main Authors Danis, Mustafa Engin, Truong, Duc P, DeSantis, Derek, Petersen, Mark, Rasmussen, Kim O, Alexandrov, Boian S
Format Paper
LanguageEnglish
Published Ithaca Cornell University Library, arXiv.org 07.08.2024
Subjects
Acceleration
Algorithms
Essentially non-oscillatory schemes
Finite volume method
Fluxes
Interpolation
Mathematical analysis
Reconstruction
Shallow water equations
Solvers
Taylor series
Tensors
Thickness
Upwind schemes (mathematics)
Online AccessGet full text

Cover

More Information
Summary:In this paper, we introduce a high-order tensor-train (TT) finite volume method for the Shallow Water Equations (SWEs). We present the implementation of the \(3^{rd}\) order Upwind and the \(5^{th}\) order Upwind and WENO reconstruction schemes in the TT format. It is shown in detail that the linear upwind schemes can be implemented by directly manipulating the TT cores while the WENO scheme requires the use of TT cross interpolation for the nonlinear reconstruction. In the development of numerical fluxes, we directly compute the flux for the linear SWEs without using TT rounding or cross interpolation. For the nonlinear SWEs where the TT reciprocal of the shallow water layer thickness is needed for fluxes, we develop an approximation algorithm using Taylor series to compute the TT reciprocal. The performance of the TT finite volume solver with linear and nonlinear reconstruction options is investigated under a physically relevant set of validation problems. In all test cases, the TT finite volume method maintains the formal high-order accuracy of the corresponding traditional finite volume method. In terms of speed, the TT solver achieves up to 124x acceleration of the traditional full-tensor scheme.
ISSN:2331-8422