BZPEERSource author record
Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.
Catalog footprint
Published work
In this work we present the mathematical foundation of an assembly code for finite element approximations of nonlocal models with compactly supported, weakly singular kernels. We demonstrate the code on a nonlocal diffusion model in various configurations and on a two-dimensional bond-based peridynamics model. The code nlfem is published under the MIT License and can be freely downloaded.
Shape optimization methods have been proven useful for identifying interfaces in models governed by partial differential equations. Here we consider a class of shape optimization problems constrained by nonlocal equations which involve interface-dependent kernels. We derive a novel shape derivative associated to the nonlocal system model and solve the problem by established numerical techniques.
The classical local Neumann problem is well studied and solutions of this problem lie, in general, in a Sobolev space. In this work, we focus on nonlocal Neumann problems with measurable, nonnegative kernels, whose solutions require less regularity assumptions. For kernels of this kind we formulate and study the weak formulation of the nonlocal Neumann problem and we investigate a nonlocal counterpart of the Sobolev space H1 as well as a resulting nonlocal trace space. We further establish, mainly for symmetric kernels, various existence results for the weak solution of the Neumann problem and we discuss related necessary conditions. Both, homogeneous and nonhomogeneous Neumann boundary conditions are considered. In addition to that, we present a new weak formulation of a Robin problem, where we reformulate the Robin problem into a regional problem.
The implementation of finite element methods (FEMs) for nonlocal models with a finite range of interaction poses challenges not faced in the partial differential equations (PDEs) setting. For example, one has to deal with weak forms involving double integrals which lead to discrete systems having higher assembly and solving costs due to possibly much lower sparsity compared to that of FEMs for PDEs. In addition, one may encounter non-smooth integrands. In many nonlocal models, nonlocal interactions are limited to bounded neighborhoods that are ubiquitously chosen to be Euclidean balls, resulting in the challenge of dealing with intersections of such balls with the finite elements. We focus on developing recipes for the efficient assembly of FEM stiffness matrices and on the choice of quadrature rules for the double integrals that contribute to the assembly efficiency and also posses sufficient accuracy. A major feature of our recipes is the use of approximate balls, e.g., several polygonal approximations of Euclidean balls, that, among other advantages, mitigate the challenge of dealing with ball-element intersections. We provide numerical illustrations of the relative accuracy and efficiency of the several approaches we develop.
A rigorous mathematical framework is provided for a substructuring-based domain-decomposition approach for nonlocal problems that feature interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains. Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal domain-decomposition methods.
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.