A Low-Rank Tensor Method for PDE-Constrained Optimization with Isogeometric Analysis
From MaRDI portal
Publication:5208738
DOI10.1137/18M1227238zbMATH Open1429.65061arXiv1811.06797OpenAlexW3000212964MaRDI QIDQ5208738FDOQ5208738
Authors: Alexandra Bünger, Martin Stoll, Sergey Dolgov 
Publication date: 10 January 2020
Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)
Abstract: Isogeometric analysis (IGA) has become one of the most popular methods for the discretization of partial differential equations motivated by the use of NURBS for geometric representations in industry and science. A crucial challenge lies in the solution of the discretized equations, which we discuss in this talk with a particular focus on PDE-constrained optimization discretized using IGA. The discretization results in a system of large mass and stiffness matrices, which are typically very costly to assemble. To reduce the computation time and storage requirements, low-rank tensor methods have become a promising tool. We present a framework for the assembly of these matrices in low-rank form as the sum of a small number of Kronecker products. For assembly of the smaller matrices only univariate integration is required. The resulting low rank Kronecker product structure of the mass and stiffness matrices can be used to solve a PDE-constrained optimization problem without assembling the actual system matrices. We present a framework which preserves and exploits the low-rank Kronecker product format for both the matrices and the solution. We use the block AMEn method to efficiently solve the corresponding KKT system of the optimization problem. We show several numerical experiments with 3D geometries to demonstrate that the low-rank assembly and solution drastically reduces the memory demands and computing times, depending on the approximation ranks of the domain.
Full work available at URL: https://arxiv.org/abs/1811.06797
Recommendations
- Low rank tensor methods in Galerkin-based isogeometric analysis
- Constrained optimization with low-rank tensors and applications to parametric problems with PDEs
- Matrix generation in isogeometric analysis by low rank tensor approximation
- Tensor train based isogeometric analysis for PDE approximation on parameter dependent geometries
- Low-rank tensor methods for partial differential equations
- A Low-Rank Matrix Equation Method for Solving PDE-Constrained Optimization Problems
- Partial tensor decomposition for decoupling isogeometric Galerkin discretizations
- A low-rank in time approach to PDE-constrained optimization
- Low-Rank Tensor Approximations for Solving Multimarginal Optimal Transport Problems
Computational methods for sparse matrices (65F50) Multilinear algebra, tensor calculus (15A69) Iterative numerical methods for linear systems (65F10) Control/observation systems governed by partial differential equations (93C20)
Cites Work
- A practical guide to splines
- Tensor-train decomposition
- TT-cross approximation for multidimensional arrays
- A Multilinear Singular Value Decomposition
- Tensor Rank and the Ill-Posedness of the Best Low-Rank Approximation Problem
- A variational discretization concept in control constrained optimization: The linear-quadratic case
- Numerical solution of saddle point problems
- Optimization with PDE Constraints
- Title not available (Why is that?)
- Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement
- Efficient quadrature for NURBS-based isogeometric analysis
- Title not available (Why is that?)
- Regularization-robust preconditioners for time-dependent PDE-constrained optimization problems
- Finite elements and fast iterative solvers. With applications in incompressible fluid dynamics
- Alternating minimal energy methods for linear systems in higher dimensions
- Title not available (Why is that?)
- The density-matrix renormalization group in the age of matrix product states
- A new scheme for the tensor representation
- The alternating linear scheme for tensor optimization in the tensor train format
- Computation of extreme eigenvalues in higher dimensions using block tensor train format
- Solution of Linear Systems and Matrix Inversion in the TT-Format
- Title not available (Why is that?)
- Dynamical density-matrix renormalization group
- Title not available (Why is that?)
- A low-rank in time approach to PDE-constrained optimization
- A new design for the implementation of isogeometric analysis in Octave and Matlab: GeoPDEs 3.0
- Low-rank solvers for unsteady Stokes-Brinkman optimal control problem with random data
- Partial tensor decomposition for decoupling isogeometric Galerkin discretizations
- Low rank tensor methods in Galerkin-based isogeometric analysis
- Matrix generation in isogeometric analysis by low rank tensor approximation
- Optimal and reduced quadrature rules for tensor product and hierarchically refined splines in isogeometric analysis
- Efficient matrix computation for tensor-product isogeometric analysis: the use of sum factorization
- Low-rank space-time decoupled isogeometric analysis for parabolic problems with varying coefficients
- Canonical polyadic decomposition with a columnwise orthonormal factor matrix
- Low-rank solution to an optimization problem constrained by the Navier-Stokes equations
- Constrained optimization with low-rank tensors and applications to parametric problems with PDEs
Cited In (14)
- An incremental singular value decomposition approach for large-scale spatially parallel \& distributed but temporally serial data -- applied to technical flows
- A Low-Rank Matrix Equation Method for Solving PDE-Constrained Optimization Problems
- A numerical B-spline Galerkin method with proper generalized decomposition for reduced order modeling of partial differential equations
- Partial tensor decomposition for decoupling isogeometric Galerkin discretizations
- Easy and efficient preconditioning of the isogeometric mass matrix
- Low rank tensor methods in Galerkin-based isogeometric analysis
- A low-rank isogeometric solver based on Tucker tensors
- A multi-mode expansion method for boundary optimal control problems constrained by random Poisson equations
- New time domain decomposition methods for parabolic optimal control problems. II: Neumann-Neumann algorithms
- Uncertainty propagation of initial conditions in thermal models
- Space-time finite element discretization of parabolic optimal control problems with energy regularization
- Unstructured space-time finite element methods for optimal control of parabolic equations
- Tensor train based isogeometric analysis for PDE approximation on parameter dependent geometries
- Constrained optimization with low-rank tensors and applications to parametric problems with PDEs
Uses Software
This page was built for publication: A Low-Rank Tensor Method for PDE-Constrained Optimization with Isogeometric Analysis
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5208738)
