High-order Multiscale Preconditioner for Elasticity of Arbitrary Structures

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Khan, Sabit Mahmood, Mehmani, Yashar
Formato: Preprint
Publicado: 2025
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866916967093370880
author Khan, Sabit Mahmood
Mehmani, Yashar
author_facet Khan, Sabit Mahmood
Mehmani, Yashar
contents We present a two-level preconditioner for solving linear systems arising from the discretization of the elliptic, linear-elastic deformation equation, in displacement unknowns, over domains that have arbitrary geometric and topological complexity and heterogeneity in material properties (including fractures). The preconditioner is an algebraic translation of the high-order pore-level multiscale method (hPLMM) proposed recently by the authors, wherein a domain is decomposed into non-overlapping subdomains, and local basis functions are numerically computed over the subdomains to construct a high-quality coarse space (or prolongation matrix). The term "high-order" stands in contrast to the recent low-order PLMM preconditioner, where BCs of local basis problems assume rigidity of all interfaces shared between subdomains. In hPLMM, interfaces are allowed to deform, through the use of suitable mortar spaces, thereby capturing local bending/twisting moments under challenging loading conditions. Benchmarked across a wide range of complex (porous) structures and material heterogeneities, we find hPLMM exhibits superior performance in Krylov solvers than PLMM, as well as state-of-the-art Schwarz and multigrid preconditioners. Applications include risk analysis of subsurface CO2/H2 storage and optimizing porous materials for batteries, prosthetics, and aircraft.
format Preprint
id arxiv_https___arxiv_org_abs_2509_19777
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle High-order Multiscale Preconditioner for Elasticity of Arbitrary Structures
Khan, Sabit Mahmood
Mehmani, Yashar
Numerical Analysis
Computational Physics
We present a two-level preconditioner for solving linear systems arising from the discretization of the elliptic, linear-elastic deformation equation, in displacement unknowns, over domains that have arbitrary geometric and topological complexity and heterogeneity in material properties (including fractures). The preconditioner is an algebraic translation of the high-order pore-level multiscale method (hPLMM) proposed recently by the authors, wherein a domain is decomposed into non-overlapping subdomains, and local basis functions are numerically computed over the subdomains to construct a high-quality coarse space (or prolongation matrix). The term "high-order" stands in contrast to the recent low-order PLMM preconditioner, where BCs of local basis problems assume rigidity of all interfaces shared between subdomains. In hPLMM, interfaces are allowed to deform, through the use of suitable mortar spaces, thereby capturing local bending/twisting moments under challenging loading conditions. Benchmarked across a wide range of complex (porous) structures and material heterogeneities, we find hPLMM exhibits superior performance in Krylov solvers than PLMM, as well as state-of-the-art Schwarz and multigrid preconditioners. Applications include risk analysis of subsurface CO2/H2 storage and optimizing porous materials for batteries, prosthetics, and aircraft.
title High-order Multiscale Preconditioner for Elasticity of Arbitrary Structures
topic Numerical Analysis
Computational Physics
url https://arxiv.org/abs/2509.19777