Source author record

Gianmarco Manzini

Gianmarco Manzini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton's principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee $C^0$-continuity of the solution at inter-element boundaries, and enforces the required $C^1$ regularity with the aid of additional variational and interior penalty terms in the weak form. The phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

preprint2022arXiv

Mesh Quality Agglomeration algorithm for the Virtual Element Method applied to Discrete Fracture Networks

We propose a quality-based optimization strategy to reduce the total number of degrees of freedom associated to a discrete problem defined over a polygonal tessellation with the Virtual Element Method. The presented Quality Agglomeration algorithm relies only on the geometrical properties of the problem polygonal mesh, agglomerating groups of neighboring elements. We test this approach in the context of fractured porous media, in which the generation of a global conforming mesh on a Discrete Fracture Network leads to a considerable number of unknowns, due to the presence of highly complex geometries and the significant size of the computational domains. We show the efficiency and the robustness of our approach, applied independently on each fracture for different network configurations, exploiting the flexibility of the Virtual Element Method in handling general polygonal elements.

preprint2022arXiv

Nonnegative Canonical Tensor Decomposition with Linear Constraints: nnCANDELINC

There is an emerging interest in tensor factorization applications in big-data analytics and machine learning. To speed up the factorization of extra-large datasets, organized in multidimensional arrays (aka tensors), easy to compute compression-based tensor representations, such as Tucker and Tensor Train formats, are used to approximate the initial large-tensor. Further, tensor factorization is used to extract latent features that can facilitate discoveries of new mechanisms and signatures hidden in the data, where the explainability of the latent features is of principal importance. Nonnegative tensor factorization extracts latent features that are naturally sparse and parts of the data, which makes them easily interpretable. However, to take into account available domain knowledge and subject matter expertise, additional constraints often need to be imposed, which lead us to Canonical decomposition with linear constraints (CANDELINC), a Canonical Polyadic Decomposition with rank deficient factors. In CANDELINC, Tucker compression is used as a pre-processing step, which leads to a larger residual error but to more explainable latent features. Here, we propose a nonnegative CANDELINC (nnCANDELINC) accomplished via a specific nonnegative Tucker decomposition; we refer to as minimal or canonical nonnegative Tucker. We derive several results required to understand the specificity of nnCANDELINC, focusing on the difficulties of preserving the nonnegative rank to its Tucker core and comparing the real-valued to the nonnegative case. Finally, we demonstrate nnCANDELINC performance on synthetic and real-world examples.

preprint2021arXiv

Extended virtual element method for two-dimensional linear elastic fracture

In this paper, we propose an eXtended Virtual Element Method (X-VEM) for two-dimensional linear elastic fracture. This approach, which is an extension of the standard Virtual Element Method (VEM), facilitates mesh-independent modeling of crack discontinuities and elastic crack-tip singularities on general polygonal meshes. For elastic fracture in the X-VEM, the standard virtual element space is augmented by additional basis functions that are constructed by multiplying standard virtual basis functions by suitable enrichment fields, such as asymptotic mixed-mode crack-tip solutions. The design of the X-VEM requires an extended projector that maps functions lying in the extended virtual element space onto a set spanned by linear polynomials and the enrichment fields. An efficient scheme to compute the mixed-mode stress intensity factors using the domain form of the interaction integral is described. The formulation permits integration of weakly singular functions to be performed over the boundary edges of the element. Numerical experiments are conducted on benchmark mixed-mode linear elastic fracture problems that demonstrate the sound accuracy and optimal convergence in energy of the proposed formulation.

preprint2021arXiv

Stability and conservation properties of Hermite-based approximations of the Vlasov-Poisson system

Spectral approximation based on Hermite-Fourier expansion of the Vlasov-Poisson model for a collisionless plasma in the electro-static limit is provided, by including high-order artificial collision operators of Lenard-Bernstein type. These differential operators are suitably designed in order to preserve the physically-meaningful invariants (number of particles, momentum, energy). In view of time-discretization, stability results in appropriate norms are presented. In this study, necessary conditions link the magnitude of the artificial collision term, the number of spectral modes of the discretization, as well as the time-step. The analysis, carried out in full for the Hermite discretization of a simple linear problem in one-dimension, is then partly extended to cover the complete nonlinear Vlasov-Poisson model.

preprint2021arXiv

Stabilization of the nonconforming virtual element method

We address the issue of designing robust stabilization terms for the nonconforming virtual element method. To this end, we transfer the problem of defining the stabilizing bilinear form from the elemental nonconforming virtual element space, whose functions are not known in closed form, to the dual space spanned by the known functionals providing the degrees of freedom. By this approach, we manage to construct different bilinear forms yielding optimal or quasi-optimal stability bounds and error estimates, under weaker assumptions on the tessellation than the ones usually considered in this framework. In particular, we prove optimality under geometrical assumptions allowing a mesh to have a very large number of arbitrarily small edges per element. Finally, we numerically assess the performance of the VEM for several different stabilizations fitting with our new framework on a set of representative test cases.

preprint2021arXiv

The conforming virtual element method for polyharmonic and elastodynamics problems: a review

In this paper, we review recent results on the conforming virtual element approximation of polyharmonic and elastodynamics problems. The structure and the content of this review is motivated by three paradigmatic examples of applications: classical and anisotropic Cahn-Hilliard equation and phase field models for brittle fracture, that are briefly discussed in the first part of the paper. We present and discuss the mathematical details of the conforming virtual element approximation of linear polyharmonic problems, the classical Cahn-Hilliard equation and linear elastodynamics problems.

preprint2021arXiv

The role of mesh quality and mesh quality indicators in the Virtual Element Method

Since its introduction, the Virtual Element Method (VEM) was shown to be able to deal with a large variety of polygons, while achieving good convergence rates. The regularity assumptions proposed in the VEM literature to guarantee the convergence on a theoretical basis are therefore quite general. They have been deduced in analogy to the similar conditions developed in the Finite Element Methods (FEMs) analysis. In this work, we experimentally show that the VEM still converges with almost optimal rates and low errors in the L2 and H1 norms even if we significantly break the regularity assumptions that are used in the literature. These results suggest that the regularity assumptions proposed so far might be overestimated. We also exhibit examples on which the VEM sub-optimally converges or diverges. Finally, we introduce a mesh quality indicator that experimentally correlates the entity of the violation of the regularity assumptions and the performance of the VEM solution, thus predicting if a dataset is potentially critical for VEM.

preprint2021arXiv

VEM and the Mesh

In this work we report some results, obtained within the framework of the ERC Project CHANGE, on the impact on the performance of the virtual element method of the shape of the polygonal elements of the underlying mesh. More in detail, after reviewing the state of the art, we present a) an experimental analysis of the convergence of the VEM under condition violating the standard shape regularity assumptions, b) an analysis of the correlation between some mesh quality metrics and a set of different performance indexes, and c) a suitably designed mesh quality indicator, aimed at predicting the quality of the performance of the VEM on a given mesh.

preprint2016arXiv

Annotations on the virtual element method for second-order elliptic problems

This document contains working annotations on the Virtual Element Method (VEM) for the approximate solution of diffusion problems with variable coefficients. To read this document you are assumed to have familiarity with concepts from the numerical discretization of Partial Differential Equations (PDEs) and, in particular, the Finite Element Method (FEM). This document is not an introduction to the FEM, for which many textbooks are available. Eventually, this document will evolve into a tutorial introduction to the VEM (but this is really a long-term goal).

preprint2016arXiv

Discretization of mixed formulations of elliptic problems on polyhedral meshes

We review basic design principles underpinning the construction of mimetic finite difference and a few finite volume and finite element schemes for mixed formulations of elliptic problems. For a class of low-order mixed-hybrid schemes, we show connections between these principles and prove that the consistency and stability conditions must lead to a member of the mimetic family of schemes regardless of the selected discretization framework. Finally, we give two examples of using flexibility of the mimetic framework: derivation of higher-order schemes and convergent schemes for nonlinear problems with small diffusion coefficients.

preprint2016arXiv

The non-conforming virtual element method for the Stokes equations

We present the non-conforming Virtual Element Method (VEM) for the numerical approximation of velocity and pressure in the steady Stokes problem. The pressure is approximated using discontinuous piecewise polynomials, while each component of the velocity is approximated using the nonconforming virtual element space. On each mesh element the local virtual space contains the space of polynomials of up to a given degree, plus suitable non-polynomial functions. The virtual element functions are implicitly defined as the solution of local Poisson problems with polynomial Neumann boundary conditions. As typical in VEM approaches, the explicit evaluation of the non-polynomial functions is not required. This approach makes it possible to construct nonconforming (virtual) spaces for any polynomial degree regardless of the parity, for two-and three-dimensional problems, and for meshes with very general polygonal and polyhedral elements. We show that the non-conforming VEM is inf-sup stable and establish optimal a priori error estimates for the velocity and pressure approximations. Numerical examples confirm the convergence analysis and the effectiveness of the method in providing high-order accurate approximations.

preprint2015arXiv

Conforming and nonconforming virtual element methods for elliptic problems

We present in a unified framework new conforming and nonconforming Virtual Element Methods (VEM) for general second order elliptic problems in two and three dimensions. The differential operator is split into its symmetric and non-symmetric parts and conditions for stability and accuracy on their discrete counterparts are established. These conditions are shown to lead to optimal $H^1$- and $L^2$-error estimates, confirmed by numerical experiments on a set of polygonal meshes. The accuracy of the numerical approximation provided by the two methods is shown to be comparable.