Guglielmo Scovazzi

Efthymios N. Karatzas, Giovanni Stabile, Leo Nouveau et al.

We propose a model order reduction technique integrating the Shifted Boundary Method (SBM) with a POD-Galerkin strategy. This approach allows to treat more complex parametrized domains in an efficient and straightforward way. The impact of the proposed approach is threefold. First, problems involving parametrizations of complex geometrical shapes and/or large domain deformations can be efficiently solved at full-order by means of the SBM, an unfitted boundary method that avoids remeshing and the tedious handling of cut cells by introducing an approximate surrogate boundary. Second, the computational effort is further reduced by the development of a reduced order model (ROM) technique based on a POD-Galerkin approach. Third, the SBM provides a smooth mapping from the true to the surrogate domain, and for this reason, the stability and performance of the reduced order basis are enhanced. This feature is the net result of the combination of the proposed ROM approach and the SBM. Similarly, the combination of the SBM with a projection-based ROM gives the great advantage of an easy and fast to implement algorithm considering geometrical parametrization with large deformations. The transformation of each geometry to a reference geometry (morphing) is in fact not required. These combined advantages will allow the solution of PDE problems more efficiently. We illustrate the performance of this approach on a number of two-dimensional Stokes flow problems.

David Michael Riley, Guglielmo Scovazzi, Ioannis Stefanou

Porous media containing cracks, fractures, or internal discontinuities arise throughout subsurface geomechanics, biomechanics, and materials science. Numerical simulation of the coupled hydromechanical response is inherently challenging because the pressure and displacement fields are tightly coupled through the Biot equations, requiring stable mixed formulations. These difficulties are compounded when cracks are present, because standard mesh-conforming approaches require costly, labor-intensive, body-fitted meshing, while unfitted methods often require cut-cell integration, enrichment functions, or additional stabilization. In this work, we use an alternative approach, we adapt the shifted interface method to coupled transient poroelasticity with embedded interfaces. The method replaces the true crack by a surrogate approximation where interface conditions are transferred through local expansions. A unified derivation yields shifted forms for both hydraulic transmission and mechanical traction coupling. Two enforcement strategies are extensively compared: a weak (integral) enforcement and a strong (pointwise) enforcement. Four test cases of increasing geometric complexity (offset mesh-aligned, boundary-intersecting angled, embedded angled, and multi-crack configurations) validate the formulation. Away from crack tips, interface residuals converge as O(h); near tips, localized post-processing artifacts degrade the global rate, but first-order convergence is recovered when a small tip region is excluded. A multi-crack demonstration with four simultaneously embedded cracks of distinct geometry and interface properties confirms the practical applicability of the framework. These results support the shifted interface method as a practical framework for poroelastic crack modeling on non-body-fitted meshes with geometrically complex embedded interfaces.