Coupled 3D-1D problems arise in many practical applications, in an attempt to reduce the computational burden in simulations where cylindrical inclusions with a small section are embedded in a much larger domain. Nonetheless the resolution of such problems can be non trivial, both from a mathematical and a geometrical standpoint. Indeed 3D-1D coupling requires to operate in non standard function spaces, and, also, simulation geometries can be complex for the presence of multiple intersecting domains. Recently, a PDE-constrained optimization based formulation has been proposed for such problems, proving a well posed mathematical formulation and allowing for the use of non conforming meshes for the discrete problem. Here an unconstrained optimization formulation of the problem is derived and an efficient gradient based solver is proposed for such formulation, along with a suitable preconditioner to speed up the iterative solver. Some numerical tests on quite complex configurations are discussed to show the viability of the method.

A gradient based resolution strategy for a PDE-constrained optimization approach for 3D-1D coupled problems / Berrone, Stefano; Grappein, Denise; Scialò, Stefano; Vicini, Fabio. - In: GEM. - ISSN 1869-2680. - 13:(2022), pp. 1-25. [10.1007/s13137-021-00192-0]

A gradient based resolution strategy for a PDE-constrained optimization approach for 3D-1D coupled problems

Berrone, Stefano;Grappein Denise;Scialò Stefano;Vicini Fabio
2022

Abstract

Coupled 3D-1D problems arise in many practical applications, in an attempt to reduce the computational burden in simulations where cylindrical inclusions with a small section are embedded in a much larger domain. Nonetheless the resolution of such problems can be non trivial, both from a mathematical and a geometrical standpoint. Indeed 3D-1D coupling requires to operate in non standard function spaces, and, also, simulation geometries can be complex for the presence of multiple intersecting domains. Recently, a PDE-constrained optimization based formulation has been proposed for such problems, proving a well posed mathematical formulation and allowing for the use of non conforming meshes for the discrete problem. Here an unconstrained optimization formulation of the problem is derived and an efficient gradient based solver is proposed for such formulation, along with a suitable preconditioner to speed up the iterative solver. Some numerical tests on quite complex configurations are discussed to show the viability of the method.
File in questo prodotto:
File Dimensione Formato  
ARTicle.pdf

non disponibili

Descrizione: Articolo
Tipologia: 1. Preprint / submitted version [pre- review]
Licenza: Non Pubblico - Accesso privato/ristretto
Dimensione 4.55 MB
Formato Adobe PDF
4.55 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
3D1DGEM.pdf

non disponibili

Descrizione: Articolo
Tipologia: 2a Post-print versione editoriale / Version of Record
Licenza: Non Pubblico - Accesso privato/ristretto
Dimensione 2.61 MB
Formato Adobe PDF
2.61 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
ARTicleBW.pdf

embargo fino al 09/12/2022

Descrizione: Articolo
Tipologia: 2. Post-print / Author's Accepted Manuscript
Licenza: PUBBLICO - Tutti i diritti riservati
Dimensione 4.56 MB
Formato Adobe PDF
4.56 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
Pubblicazioni consigliate

Caricamento pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11583/2943832