A Low-Rank QTT-based Finite Element Method for Elasticity Problems

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Benvenuti, Elena, Manzini, Gianmarco, Nale, Marco, Pizzolato, Simone
Format: Preprint
Publié: 2025
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866929674998775808
author Benvenuti, Elena
Manzini, Gianmarco
Nale, Marco
Pizzolato, Simone
author_facet Benvenuti, Elena
Manzini, Gianmarco
Nale, Marco
Pizzolato, Simone
contents We present an efficient and robust numerical algorithm for solving the two-dimensional linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. Our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented, which includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method's performance in the numerical approximation of three representative test cases.
format Preprint
id arxiv_https___arxiv_org_abs_2501_07778
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A Low-Rank QTT-based Finite Element Method for Elasticity Problems
Benvenuti, Elena
Manzini, Gianmarco
Nale, Marco
Pizzolato, Simone
Numerical Analysis
We present an efficient and robust numerical algorithm for solving the two-dimensional linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. Our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented, which includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method's performance in the numerical approximation of three representative test cases.
title A Low-Rank QTT-based Finite Element Method for Elasticity Problems
topic Numerical Analysis
url https://arxiv.org/abs/2501.07778