<p>Structured finite element methods (FEMs) based on low-rank approximation in the form of the so-called <i>quantized tensor train (QTT)</i> decomposition (QTT-FEM) have been proposed and extensively studied in the case of elliptic equations. In this work, we design a QTT-FE method for time-domain acoustic wave equations, combining stable low-rank approximation in space with a suitable conservative discretization in time. For the acoustic wave equation with a homogeneous source term in a single space dimension as a model problem, we consider its reformulation as a first-order system in time. In space, we employ a low-rank QTT-FEM discretization based on continuous piecewise linear finite elements corresponding to uniformly refined nested meshes. Time integration is performed using symplectic high-order Gauss–Legendre Runge–Kutta methods. In our numerical experiments, we investigate the energy conservation and exponential convergence of the proposed method.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Symplectic QTT-FEM solution of the one-dimensional acoustic wave equation in the time domain

  • Sara Fraschini,
  • Vladimir Kazeev,
  • Ilaria Perugia

摘要

Structured finite element methods (FEMs) based on low-rank approximation in the form of the so-called quantized tensor train (QTT) decomposition (QTT-FEM) have been proposed and extensively studied in the case of elliptic equations. In this work, we design a QTT-FE method for time-domain acoustic wave equations, combining stable low-rank approximation in space with a suitable conservative discretization in time. For the acoustic wave equation with a homogeneous source term in a single space dimension as a model problem, we consider its reformulation as a first-order system in time. In space, we employ a low-rank QTT-FEM discretization based on continuous piecewise linear finite elements corresponding to uniformly refined nested meshes. Time integration is performed using symplectic high-order Gauss–Legendre Runge–Kutta methods. In our numerical experiments, we investigate the energy conservation and exponential convergence of the proposed method.