Researcher profile

Andrea Bonito

Andrea Bonito contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2022arXiv

Approximation of the spectral fractional powers of the Laplace-Beltrami Operator

We consider numerical approximations of spectral fractional Laplace-Beltrami problems on closed surfaces. The proposed numerical algorithms rely on their Balakrishnan integral representation and consist of a sinc quadrature coupled with standard finite element methods for parametric surfaces. Possibly up to a log term, optimal rates of convergence are observed and derived analytically when the discrepancies between the exact solution and its numerical approximations are measured in $L^2$ and $H^1$. The performances of the algorithms are illustrated in different settings including the approximation of Gaussian fields on surfaces.

preprint2022arXiv

Error Estimates For A Linear Folding Model

An interior penalty discontinuous Galerkin method is devised to approximate minimizers of a linear folding model by discontinuous isoparametric finite element functions that account for an approximation of a folding arc. The numerical analysis of the discrete model includes an a priori error estimate in case of an accurate representation of the folding curve by the isoparametric mesh. Additional estimates show that geometric consistency errors may be controlled separately if the folding arc is approximated by piecewise polynomial curves. Various numerical experiments are carried out to validate the a priori error estimate for the folding model.

preprint2022arXiv

Modeling and simulation of thin sheet folding

The article addresses the mathematical modeling of the folding of a thin elastic sheet along a prescribed curved arc. A rigorous model reduction from a general hyperelastic material description is carried out under appropriate scaling conditions on the energy and the geometric properties of the folding arc in dependence on the small sheet thickness. The resulting two-dimensional model is a piecewise nonlinear Kirchhoff plate bending model with a continuity condition at the folding arc. A discontinuous Galerkin method and an iterative scheme are devised for the accurate numerical approximation of large deformations.

preprint2020arXiv

A divergence-conforming finite element method for the surface Stokes equation

The Stokes equation posed on surfaces is important in some physical models, but its numerical solution poses several challenges not encountered in the corresponding Euclidean setting. These include the fact that the velocity vector should be tangent to the given surface and the possible presence of degenerate modes (Killing fields) in the solution. We analyze a surface finite element method which provides solutions to these challenges. We consider an interior penalty method based on the well-known Brezzi-Douglas-Marini $H({\rm div})$-conforming finite element space. The resulting spaces are tangential to the surface, but require penalization of jumps across element interfaces in order to weakly maintain $H^1$ conformity of the velocity field. In addition our method exactly satisfies the incompressibility constraint in the surface Stokes problem. Secondly, we give a method which robustly filters Killing fields out of the solution. This problem is complicated by the fact that the dimension of the space of Killing fields may change with small perturbations of the surface. We first approximate the Killing fields via a Stokes eigenvalue problem and then give a method which is asymptotically guaranteed to correctly exclude them from the solution. The properties of our method are rigorously established via an error analysis and illustrated via numerical experiments.

preprint2020arXiv

Finite Element Approximation of a Strain-Limiting Elastic Model

We construct a finite element approximation of a strain-limiting elastic model on a bounded open domain in $\mathbb{R}^d$, $d \in \{2,3\}$. The sequence of finite element approximations is shown to exhibit strong convergence to the unique weak solution of the model. Assuming that the material parameters featuring in the model are Lipschitz-continuous, and assuming that the weak solution has additional regularity, the sequence of finite element approximations is shown to converge with a rate. An iterative algorithm is constructed for the solution of the system of nonlinear algebraic equations that arises from the finite element approximation. An appealing feature of the iterative algorithm is that it decouples the monotone and linear elastic parts of the nonlinearity in the model. In particular, our choice of piecewise constant approximation for the stress tensor (and continuous piecewise linear approximation for the displacement) allows us to compute the monotone part of the nonlinearity by solving an algebraic system with $d(d+1)/2$ unknowns independently on each element in the subdivision of the computational domain. The theoretical results are illustrated by numerical experiments.

preprint2020arXiv

Nonlinear Methods for Model Reduction

The usual approach to model reduction for parametric partial differential equations (PDEs) is to construct a linear space $V_n$ which approximates well the solution manifold $\mathcal{M}$ consisting of all solutions $u(y)$ with $y$ the vector of parameters. This linear reduced model $V_n$ is then used for various tasks such as building an online forward solver for the PDE or estimating parameters from data observations. It is well understood in other problems of numerical computation that nonlinear methods such as adaptive approximation, $n$-term approximation, and certain tree-based methods may provide improved numerical efficiency. For model reduction, a nonlinear method would replace the linear space $V_n$ by a nonlinear space $Σ_n$. This idea has already been suggested in recent papers on model reduction where the parameter domain is decomposed into a finite number of cells and a linear space of low dimension is assigned to each cell. Up to this point, little is known in terms of performance guarantees for such a nonlinear strategy. Moreover, most numerical experiments for nonlinear model reduction use a parameter dimension of only one or two. In this work, a step is made towards a more cohesive theory for nonlinear model reduction. Framing these methods in the general setting of library approximation allows us to give a first comparison of their performance with those of standard linear approximation for any general compact set. We then turn to the study these methods for solution manifolds of parametrized elliptic PDEs. We study a very specific example of library approximation where the parameter domain is split into a finite number $N$ of rectangular cells and where different reduced affine spaces of dimension $m$ are assigned to each cell. The performance of this nonlinear procedure is analyzed from the viewpoint of accuracy of approximation versus $m$ and $N$.

preprint2020arXiv

Numerical Simulations of Surface-Quasi Geostrophic Flows on Periodic Domains

We propose a novel algorithm for the approximation of surface-quasi geostrophic (SQG) flows modeled by a nonlinear partial differential equation coupling transport and fractional diffusion phenomena. The time discretization consists of an explicit strong-stability-preserving three-stage Runge-Kutta method while a flux-corrected-transport (FCT) method coupled with Dunford-Taylor representations of fractional operators is advocated for the space discretization. Standard continuous piecewise linear finite elements are employed and the algorithm does not have restrictions on the mesh structure nor on the computational domain. In the inviscid case, we show that the resulting scheme satisfies a discrete maximum principle property under a standard CFL condition and observe, in practice, its second-order accuracy in space. The algorithm successfully approximates several benchmarks with sharp transitions and fine structures typical of SQG flows. In addition, theoretical Kolmogorov energy decay rates are observed on a freely decaying atmospheric turbulence simulation.

preprint2020arXiv

Polynomial Approximation of Anisotropic Analytic Functions of Several Variables

Motivated by numerical methods for solving parametric partial differential equations, this paper studies the approximation of multivariate analytic functions by algebraic polynomials. We introduce various anisotropic model classes based on Taylor expansions, and study their approximation by finite dimensional polynomial spaces $\cal{P}_Λ$ described by lower sets $Λ$. Given a budget $n$ for the dimension of $\cal{P}_Λ$, we prove that certain lower sets $Λ_n$, with cardinality $n$, provide a certifiable approximation error that is in a certain sense optimal, and that these lower sets have a simple definition in terms of simplices. Our main goal is to obtain approximation results when the number of variables $d$ is large and even infinite, and so we concentrate almost exclusively on the case $d=\infty$. We also emphasize obtaining results which hold for the full range $n\ge 1$, rather than asymptotic results that only hold for $n$ sufficiently large. In applications, one typically wants $n$ small to comply with computational budgets.

preprint2020arXiv

Reduced Basis Approximations of the Solutions to Spectral Fractional Diffusion Problems

We consider the numerical approximation of the spectral fractional diffusion problem based on the so called Balakrishnan representation. The latter consists of an improper integral approximated via quadratures. At each quadrature point, a reaction-diffusion problem must be approximated and is the method bottle neck. In this work, we propose to reduce the computational cost using a reduced basis strategy allowing for a fast evaluation of the reaction-diffusion problems. The reduced basis does not depend on the fractional power $s$ for $0<s_{\min}\leq s \leq s_{\max}<1$. It is built offline once for all and used online irrespectively of the fractional power. We analyze the reduced basis strategy and show its exponential convergence. The analytical results are illustrated with insightful numerical experiments.

preprint2019arXiv

Electroconvection of Thin Liquid Crystals: Model Reduction and Numerical Simulations

We propose a finite element method for the numerical simulation of electroconvection of thin liquid crystals. The liquid is located in between two concentric circular electrodes which are either assumed to be of infinite height or slim. Each configuration results in a different nonlocal electro-magnetic model defined on a two dimensional bounded domain. The numerical method consists in approximating the surface charge density, the liquid velocity and pressure, and the electric potential in the two dimensional liquid region. Finite elements for the space discretization coupled with standard time stepping methods are put forward. Unlike for the infinite electrodes configuration, our numerical simulations indicate that slim electrodes are favorable for electroconvection to occur and are able to sustain the phenomena over long period of time. Furthermore, we provide a numerical study on the influence of the three main parameters of the system: the Rayleigh number, the Prandtl number and the electrodes aspect ratio.