METaL Seminars

Gabriele Loli

Efficient solvers for Isogeometric Analysis

We discuss efficient linear solvers for PDEs discretized through Isogeometric Analysis (IGA). In IGA, spline and Non-Uniform Rational B-Spline (NURBS) bases are used both to describe the geometry and to build the approximation spaces, allowing one to work with high-order smooth functions while keeping an exact representation of the domain. The methods presented here rely on the tensor-product structure of multivariate B-splines. In particular, we exploit the Kronecker-product structure of the matrices arising from Galerkin discretizations to construct preconditioners for several classes of PDEs.

For elliptic and parabolic problems, we use Fast Diagonalization (together with some extensions) to invert preconditioners expressed as sums of Kronecker products of univariate matrices. For hyperbolic problems, in particular the wave equation, we instead turn to Schur decompositions, which allow us to design solvers that avoid assembling very large space–time matrices.

A set of numerical experiments will be shown to illustrate the performance of the proposed strategies.

Antonio Maria D'Altri

Computational Analysis of Masonry and Heritage Structures:

recent advances across aifferent scales

This seminar provides an overview of computational modeling techniques applied to masonry and heritage structures, emphasizing recent advances across different scales, i.e., from a full-scale monument (e.g., medieval castle) down to the micrometer scale of a single pore in porous building material. The presentation begins by addressing the structural scale, focusing on the numerical modeling and structural analysis of full-scale heritage structures. A key challenge discussed is finite element (FE) mesh generation from point clouds, featuring procedures such as Cloud2FEM. This is followed by an overview of analysis approaches for conducting structural assessments on large-scale historical structures. Next, the focus narrows to a high-fidelity block-based mechanical model for masonry, where damaging blocks interact through cohesive-frictional contact constraints. Applications, including those in the dynamic regime under earthquake-like actions, will be shown. A recent advancement is the "continualization" of this block-based model through deriving a mechanism-based strength domain and implementing it within plasticity and damage computational frameworks. Finally, the discussion shifts to the micrometer scale with the simulation of salt crystallization-induced damage in porous materials. This involves multiphase modeling of water and salt transport within porous building materials, coupled with micromechanics-based damage estimations. The seminar will conclude with a brief mention of advances in the Virtual Element Method (VEM) and an overview of the Marie Skłodowska-Curie Postdoctoral Fellowships (MSCA PF) call, highlighting this opportunity for future research funding for young investigators.

Matteo Negri

From Energy Functionals to Cracks: 

bridging analysis and computation in phase-field fracture

We present a general approach for quasi-static and dynamic fracture propagation based on energy functionals. Quasi-static (steady-state) evolutions are governed by energy identity and equilibrium, with respect to displacement and crack variations. We present these basic ingredients and the resulting system of evolution equations for a toy sharp-crack model and then for a phase-field model. In the latter case the evolution is computed by means of constrained incremental problems, which we solve using a staggered algorithm. Consistency between sharp and phase-field approaches is supported by Γ-convergence and numerical simulations, however the overall theoretical and computational picture is far from being complete and raises open problems on convergence of energies, algorithms and evolutions. For dynamic problems, evolutions are governed by energy identity and equilibrium (as in quasi-static) together with stationary action. Again, we compute phase-field evolutions by a staggered scheme where the wave equation is solved by a Newmark method, in a way that numerical dissipation is negligible; as a result simulations feature energy conservation and a precise representation of P- and S-waves. In this setting, high order phase-field energies turn out to be crucial to reduce the computational cost without compromising accuracy.