DYNAMICAL LOW-RANK APPROXIMATION FOR BURGERS' EQUATION WITH UNCERTAINTY

被引:6
作者
Kusch, Jonas [1 ]
Ceruti, Gianluca [2 ]
Einkemmer, Lukas [3 ]
Frank, Martin [1 ]
机构
[1] Karlsruhe Inst Technol, Karlsruhe, Germany
[2] Univ Tubingen, Tubingen, Germany
[3] Univ Innsbruck, Innsbruck, Austria
关键词
uncertainty quantification; conservation laws; hyperbolic; intrusive UQ methods; dynamical low-rank approximation; matrix projector-splitting integrator; unconventional integrator; PROJECTOR-SPLITTING INTEGRATOR; FINITE-VOLUME METHODS; CONSERVATION-LAWS; ORTHOGONAL APPROXIMATION; DIFFERENTIAL-EQUATIONS; TIME INTEGRATION; POLYNOMIAL CHAOS; FLUID; QUANTIFICATION; PROPAGATION;
D O I
10.1615/Int.J.UncertaintyQuantification.2022039345
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
Quantifying uncertainties in hyperbolic equations is a source of several challenges. First, the solution forms shocks leading to oscillatory behavior in the numerical approximation of the solution. Second, the number of unknowns required for an effective discretization of the solution grows exponentially with the dimension of the uncertainties, yielding high computational costs and large memory requirements. The number of unknowns can be reduced through generalized polynomial chaos polynomials, which allow for an efficient representation when the distribution of the uncertainties is known. These distributions are usually only available for input uncertainties such as initial conditions; therefore the efficiency of this ansatz can be lost during runtime. In this paper, we make use of the dynamical low-rank approximation (DLRA) to obtain a memorywise efficient solution approximation on a lower-dimensional manifold for Burgers' equation. We investigate and compare the use of the matrix projector-splitting integrator and the unconventional integrator for DLRA, deriving efficient time evolution equations for the spatial and uncertain basis functions, respectively. This guarantees an accurate approximation of the solution even if the underlying probability distributions change over time. The proposed methodology is analyzed for Burgers' equation equipped with uncertain initial values represented by a two-dimensional random vector. The numerical experiments show a significant reduction of the memory requirements, while important characteristics of the original system are well captured.
引用
收藏
页码:1 / 21
页数:21
相关论文
共 50 条
[31]   Efficient dynamical low-rank approximation for the Vlasov-Ampere-Fokker-Planck system [J].
Coughlin, Jack ;
Hu, Jingwei .
JOURNAL OF COMPUTATIONAL PHYSICS, 2022, 470
[32]   Dynamical low-rank integrators for second-order matrix differential equations [J].
Hochbruck, Marlis ;
Neher, Markus ;
Schrammer, Stefan .
BIT NUMERICAL MATHEMATICS, 2023, 63 (01)
[33]   A LOW-RANK PROJECTOR-SPLITTING INTEGRATOR FOR THE VLASOV-POISSON EQUATION [J].
Einkemmer, Lukas ;
Lubich, Christian .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2018, 40 (05) :B1330-B1360
[34]   A new splitting algorithm for dynamical low-rank approximation motivated by the fibre bundle structure of matrixmanifolds [J].
Billaud-Friess, Marie ;
Falco, Antonio ;
Nouy, Anthony .
BIT NUMERICAL MATHEMATICS, 2022, 62 (02) :387-408
[35]   Numerical analysis of the Burgers' equation in the presence of uncertainty [J].
Pettersson, Per ;
Iaccarino, Gianluca ;
Nordstrom, Jan .
JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 228 (22) :8394-8412
[36]   A PREDICTOR-CORRECTOR STRATEGY FOR ADAPTIVITY IN DYNAMICAL LOW-RANK APPROXIMATIONS [J].
Hauck, Cory D. ;
Schnake, Stefan .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2023, 44 (03) :971-1005
[37]   ASYMPTOTIC-PRESERVING AND ENERGY STABLE DYNAMICAL LOW-RANK APPROXIMATION FOR THERMAL RADIATIVE TRANSFER EQUATIONS [J].
Frank, Martin ;
Kusch, Jonas ;
Patwardhan, Chinmay .
MULTISCALE MODELING & SIMULATION, 2025, 23 (01) :278-312
[38]   Numerical low-rank approximation of matrix differential equations [J].
Mena, Hermann ;
Ostermann, Alexander ;
Pfurtscheller, Lena-Maria ;
Piazzola, Chiara .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2018, 340 :602-614
[39]   A new splitting algorithm for dynamical low-rank approximation motivated by the fibre bundle structure of matrix manifolds [J].
Marie Billaud-Friess ;
Antonio Falcó ;
Anthony Nouy .
BIT Numerical Mathematics, 2022, 62 :387-408
[40]   Bayesian Uncertainty Quantification for Low-Rank Matrix Completion [J].
Yuchi, Henry Shaowu ;
Mak, Simon ;
Xie, Yao .
BAYESIAN ANALYSIS, 2023, 18 (02) :491-518