An overlapping domain decomposition method for parametric Stokes and Stokes-Darcy problems via proper generalized decomposition

A strategy to construct physics-based local surrogate models for parametric Stokes flows and coupled Stokes-Darcy systems is presented. The methodology relies on the proper generalized decomposition (PGD) method to reduce the dimensionality of the parametric flow fields and on an overlapping domain...

Descripción completa

Detalles Bibliográficos
Autores: Discacciati, Marco|||0000-0001-8343-8953, Evans, Ben J., Giacomini, Matteo|||0000-0001-6094-5944
Tipo de recurso: artículo
Fecha de publicación:2026
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:dnet:upcommonspor::fc13d82fb915f6e84b33f3323d0e2c77
Acceso en línea:https://hdl.handle.net/2117/461076
https://dx.doi.org/10.1016/j.cma.2026.118878
Access Level:acceso abierto
Palabra clave:Reduced order models
Proper generalized decomposition
Overlapping domain decomposition
Non-intrusiveness
Stokes-Darcy
Parametric flows
Àrees temàtiques de la UPC::Matemàtiques i estadística::Anàlisi numèrica::Mètodes en elements finits
Descripción
Sumario:A strategy to construct physics-based local surrogate models for parametric Stokes flows and coupled Stokes-Darcy systems is presented. The methodology relies on the proper generalized decomposition (PGD) method to reduce the dimensionality of the parametric flow fields and on an overlapping domain decomposition (DD) paradigm to reduce the number of globally coupled degrees of freedom in space. The DD-PGD approach provides a non-intrusive framework in which end-users only need access to the matrices arising from the (finite element) discretization of the full-order problems in the subdomains. The traces of the finite element functions used for the discretization within the subdomains are employed to impose arbitrary Dirichlet boundary conditions at the interface, without introducing auxiliary basis functions. The methodology is seamless to the choice of the discretization schemes in space, being compatible with both LBB-compliant finite element pairs and stabilized formulations, and the DD-PGD paradigm is transparent to the employed overlapping DD approach. The local surrogate models are glued together in the online phase by solving a parametric interface system to impose continuity of the subdomain solutions at the interfaces, without introducing Lagrange multipliers to enforce the continuity in the entire overlap and without solving any additional physical problem in the reduced space. Numerical results are presented for parametric single-physics (Stokes-Stokes) and multi-physics (Stokes-Darcy) systems, showcasing the accuracy, robustness, and computational efficiency of DD-PGD.