BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrea Bonito

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.APphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

20works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

20 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence and optimality of higher-order adaptive finite element methods for eigenvalue clusters

Proofs of convergence of adaptive finite element methods for the approximation of eigenvalues and eigenfunctions of linear elliptic problems have been given in a several recent papers. A key step in establishing such results for multiple and clustered eigenvalues was provided by Dai et. al. (2014), who proved convergence and optimality of AFEM for eigenvalues of multiplicity greater than one. There it was shown that a theoretical (non-computable) error estimator for which standard convergence proofs apply is equivalent to a standard computable estimator on sufficiently fine grids. Gallistl (2015) used a similar tool in order to prove that a standard adaptive FEM for controlling eigenvalue clusters for the Laplacian using continuous piecewise linear finite element spaces converges with optimal rate. When considering either higher-order finite element spaces or non-constant diffusion coefficients, however, the arguments of Dai et. al. and Gallistl do not yield equivalence of the practical and theoretical estimators for clustered eigenvalues. In this note we provide this missing key step, thus showing that standard adaptive FEM for clustered eigenvalues employing elements of arbitrary polynomial degree converge with optimal rate. We additionally establish that a key user-defined input parameter in the AFEM, the bulk marking parameter, may be chosen entirely independently of the properties of the target eigenvalue cluster. All of these results assume a fineness condition on the initial mesh in order to ensure that the nonlinearity is sufficiently resolved.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Diffusion Coefficients Estimation for Elliptic Partial Differential Equations

This paper considers the Dirichlet problem $$ -\mathrm{div}(a\nabla u_a)=f \quad \hbox{on}\,\,\ D, \qquad u_a=0\quad \hbox{on}\,\,\partial D, $$ for a Lipschitz domain $D\subset \mathbb R^d$, where $a$ is a scalar diffusion function. For a fixed $f$, we discuss under which conditions is $a$ uniquely determined and when can $a$ be stably recovered from the knowledge of $u_a$. A first result is that whenever $a\in H^1(D)$, with $0<λ\le a\le Λ$ on $D$, and $f\in L_\infty(D)$ is strictly positive, then $$ \|a-b\|_{L_2(D)}\le C\|u_a-u_b\|_{H_0^1(D)}^{1/6}. $$ More generally, it is shown that the assumption $a\in H^1(D)$ can be weakened to $a\in H^s(D)$, for certain $s<1$, at the expense of lowering the exponent $1/6$ to a value that depends on $s$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-Order AFEM for the Laplace-Beltrami Operator: Convergence Rates

We present a new AFEM for the Laplace-Beltrami operator with arbitrary polynomial degree on parametric surfaces, which are globally $W^1_\infty$ and piecewise in a suitable Besov class embedded in $C^{1,α}$ with $α\in (0,1]$. The idea is to have the surface sufficiently well resolved in $W^1_\infty$ relative to the current resolution of the PDE in $H^1$. This gives rise to a conditional contraction property of the PDE module. We present a suitable approximation class and discuss its relation to Besov regularity of the surface, solution, and forcing. We prove optimal convergence rates for AFEM which are dictated by the worst decay rate of the surface error in $W^1_\infty$ and PDE error in $H^1$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Approximation of Fractional Powers of Regularly Accretive Operators

We study the numerical approximation of fractional powers of accretive operators in this paper. Namely, if $A$ is the accretive operator associated with an accretive sesquilinear form $A(\cdot,\cdot)$ defined on a Hilbert space $\mathbb V$ contained in $L^2(Ω)$, we approximate $A^{-β}$ for $β\in (0,1)$. The fractional powers are defined in terms of the so-called Balakrishnan integral formula. Given a finite element approximation space $\mathbb V_h\subset \mathbb V$, $A^{-β}$ is approximated by $A_h^{-β}π_h$ where $A_h$ is the operator associated with the form $A(\cdot,\cdot)$ restricted to $\mathbb V_h$ and $π_h$ is the $L^2(Ω)$-projection onto $\mathbb V_h$. We first provide error estimates for $(A^β-A_h^βπ_h)f$ in Sobolev norms with index in [0,1] for appropriate $f$. These results depend on elliptic regularity properties of variational solutions involving the form $A(\cdot,\cdot)$ and are valid for the case of less than full elliptic regularity. We also construct and analyze an exponentially convergent sinc quadrature approximation to the Balakrishnan integral defining $A_h^βπ_h f$. Finally, the results of numerical computations illustrating the proposed method are given.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Approximation of Parabolic Equations Involving Fractional Powers of Elliptic Operators

We study the numerical approximation of a time dependent equation involving fractional powers of an elliptic operator $L$ defined to be the unbounded operator associated with a Hermitian, coercive and bounded sesquilinear form on $H^1_0(Ω)$. The time dependent solution $u(x,t)$ is represented as a Dunford Taylor integral along a contour in the complex plane. The contour integrals are approximated using sinc quadratures. In the case of homogeneous right-hand-sides and initial value $v$, the approximation results in a linear combination of functions $(z_qI-L)^{-1}v\in H^1_0(Ω)$ for a finite number of quadrature points $z_q$ lying along the contour. In turn, these quantities are approximated using complex valued continuous piecewise linear finite elements. Our main result provides $L^2(Ω)$ error estimates between the solution $u(\cdot,t)$ and its final approximation. Numerical results illustrating the behavior of the algorithms are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Interior Penalty Method with $C^0$ Finite Elements for the Approximation of the Maxwell Equations in Heterogeneous Media: Convergence Analysis with Minimal Regularity

The present paper proposes and analyzes an interior penalty technique using $C^0$-finite elements to solve the Maxwell equations in domains with heterogeneous properties. The convergence analysis for the boundary value problem and the eigenvalue problem is done assuming only minimal regularity in Lipschitz domains. The method is shown to converge for any polynomial degrees and to be spectrally correct.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bilayer Plates: Model Reduction, $Γ$-Convergent Finite Element Approximation and Discrete Gradient Flow

The bending of bilayer plates is a mechanism which allows for large deformations via small externally induced lattice mismatches of the underlying materials. Its mathematical modeling, discussed herein, consists of a nonlinear fourth order problem with a pointwise isometry constraint. A discretization based on Kirchhoff quadrilaterals is devised and its $Γ$-convergence is proved. An iterative method that decreases the energy is proposed and its convergence to stationary configurations is investigated. Its performance, as well as reduced model capabilities, are explored via several insightful numerical experiments involving large (geometrically nonlinear) deformations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical Simulations of Bouncing Jets

Bouncing jets are fascinating phenomenons occurring under certain conditions when a jet impinges on a free surface. This effect is observed when the fluid is Newtonian and the jet falls in a bath undergoing a solid motion. It occurs also for non-Newtonian fluids when the jets falls in a vessel at rest containing the same fluid. We investigate numerically the impact of the experimental setting and the rheological properties of the fluid on the onset of the bouncing phenomenon. Our investigations show that the occurrence of a thin lubricating layer of air separating the jet and the rest of the liquid is a key factor for the bouncing of the jet to happen. The numerical technique that is used consists of a projection method for the Navier-Stokes system coupled with a level set formulation for the representation of the interface. The space approximation is done with adaptive finite elements. Adaptive refinement is shown to be very important to capture the thin layer of air that is responsible for the bouncing.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adaptive Finite Element Methods for Elliptic Problems with Discontinuous Coefficients

Elliptic partial differential equations (PDEs) with discontinuous diffusion coefficients occur in application domains such as diffusions through porous media, electro-magnetic field propagation on heterogeneous media, and diffusion processes on rough surfaces. The standard approach to numerically treating such problems using finite element methods is to assume that the discontinuities lie on the boundaries of the cells in the initial triangulation. However, this does not match applications where discontinuities occur on curves, surfaces, or manifolds, and could even be unknown beforehand. One of the obstacles to treating such discontinuity problems is that the usual perturbation theory for elliptic PDEs assumes bounds for the distortion of the coefficients in the $L_\infty$ norm and this in turn requires that the discontinuities are matched exactly when the coefficients are approximated. We present a new approach based on distortion of the coefficients in an $L_q$ norm with $q<\infty$ which therefore does not require the exact matching of the discontinuities. We then use this new distortion theory to formulate new adaptive finite element methods (AFEMs) for such discontinuity problems. We show that such AFEMs are optimal in the sense of distortion versus number of computations, and report insightful numerical results supporting our analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical Approximation of Fractional Powers of Elliptic Operators

We present and study a novel numerical algorithm to approximate the action of $T^β:=L^{-β}$ where $L$ is a symmetric and positive definite unbounded operator on a Hilbert space $H_0$. The numerical method is based on a representation formula for $T^{-β}$ in terms of Bochner integrals involving $(I+t^2L)^{-1}$ for $t\in(0,\infty)$. To develop an approximation to $T^β$, we introduce a finite element approximation $L_h$ to $L$ and base our approximation to $T^β$ on $T_h^β:= L_h^{-β}$. The direct evaluation of $T_h^β$ is extremely expensive as it involves expansion in the basis of eigenfunctions for $L_h$. The above mentioned representation formula holds for $T_h^{-β}$ and we propose three quadrature approximations denoted generically by $Q_h^β$. The two results of this paper bound the errors in the $H_0$ inner product of $T^β-T_h^βπ_h$ and $T_h^β-Q_h^β$ where $π_h$ is the $H_0$ orthogonal projection into the finite element space. We note that the evaluation of $Q_h^β$ involves application of $(I+(t_i)^2L_h)^{-1}$ with $t_i$ being either a quadrature point or its inverse. Efficient solution algorithms for these problems are available and the problems at different quadrature points can be straightforwardly solved in parallel. Numerical experiments illustrating the theoretical estimates are provided for both the quadrature error $T_h^β-Q_h^β$ and the finite element error $T^β-T_h^βπ_h$.

0 citations0 reviewsNo code linkedOpen access