Search

We present a new mimetic finite difference method for diffusion problems that converges on grids with \textit{curved} (i.e., non-planar) faces. Crucially, it gives a symmetric discrete problem that uses only one discrete unknown per curved face. The principle at the core of our construction is to abandon the standard definition of local consistency of mimetic finite difference methods. Instead, we exploit the novel and global concept of $P_{0}$-consistency. Numerical examples confirm the consistency and the optimal convergence rate of the proposed mimetic method for cubic grids with randomly perturbed nodes as well as grids with curved boundaries., Comment: Accepted manuscript

We provide a novel framework to compute a discrete vector potential of a given discrete vector field on arbitrary polyhedral meshes. The framework exploits the concept of acyclic matching, a combinatorial tool at the core of discrete Morse theory. We introduce the new concept of complete acyclic matchings and we show that they give the same end result of Gaussian elimination. Basically, instead of doing costly row and column operations on a sparse matrix, we compute equivalent cheap combinatorial operations that preserve the underlying sparsity structure. Currently, the most efficient algorithms proposed in literature to find discrete vector potentials make use of tree-cotree techniques. We show that they compute a special type of complete acyclic matchings. Moreover, we show that the problem of computing them is equivalent to the problem of deciding whether a given mesh has a topological property called collapsibility. This fact gives a topological characterization of well-known termination problems of tree-cotree techniques. We propose a new recursive algorithm to compute discrete vector potentials. It works directly on basis elements of $1$- and $2$-chains by performing elementary Gaussian operations on them associated with acyclic matchings. However, the main novelty is that it can be applied recursively. Indeed, the recursion process allows us to sidetrack termination problems of the standard tree-cotree techniques. We tested the algorithm on pathological triangulations with known topological obstructions. In all tested problems we observe linear computational complexity as a function of mesh size. Moreover, the algorithm is purely graph-based so it is straightforward to implement and does not require specialized external procedures. We believe that our framework could offer new perspectives to sparse matrix computations.

We present an in-house modelling framework for Ferroelectric Tunnelling Junctions (FTJ), and an insightful study of the design of FTJs as synaptic devices. Results show that a moderately low-k tunnelling dielectric (e.g. SiO2) can increase the read current and the current dynamic range., Comment: 2021 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copy-righted component of this work in other works

The geometric reinterpretation of the Finite Element Method (FEM) shows that Raviart Thomas and Nedelec mass matrices map from degrees of freedoms (DoFs) attached to geometric elements of a tetrahedral grid to DoFs attached to the barycentric dual grid. The algebraic inverses of the mass matrices map DoFs attached to the barycentric dual grid back to DoFs attached to the corresponding primal tetrahedral grid, but they are of limited practical use since they are dense. In this paper we present a new geometric construction of sparse inverse mass matrices for arbitrary tetrahedral grids and possibly anisotropic materials, debunking the conventional wisdom that the barycentric dual grid prohibits a sparse representation for inverse mass matrices. In particular, we provide a unified framework for the construction of both edge and face mass matrices and their sparse inverses. Such a unifying principle relies on novel geometric reconstruction formulas, from which, according to a well established design strategy, local mass matrices are constructed as the sum of a consistent and a stabilization term. A major difference with the approaches proposed so far is that the consistent term is defined geometrically and explicitly, that is, without the necessity of computing the inverses of local matrices. This provides a sensible speedup and an easier implementation. We use these new sparse inverse mass matrices to discretize a three dimensional Poisson problem, providing the comparison between the results obtained by various formulations on a benchmark problem with analytical solution., Comment: 6 figures, submitted on 29 July 2020

An efficient numerical algorithm for the computation of linking number is presented. The algorithm keep tracks or rounding error so that it can ensure the correctness of the results., Comment: 19 pages 1 figure

This paper extends the $T$-$\Omega$ formulation for eddy currents based on higher order hierarchical basis functions so that it can deal with conductors of arbitrary topology. To this aim we supplement the classical hierarchical basis functions with non-local basis functions spanning the first de Rham cohomology group of the insulating region. Such non-local basis functions may be efficiently found in negligible time with the recently introduced D{\l}otko--Specogna (DS) algorithm.

We present an efficient algorithm for the construction of a basis of $H_2(\overline{\Omega},\partial\Omega;\mathbb Z)$ via the Poincar\'e--Lefschetz duality theorem. Denoting by $g$ the first Betti number of $\overline \Omega$ the idea is to find, first $g$ different $1$-boundaries of $\overline{\Omega}$ with supports contained in $\partial\Omega$ whose homology classes in $\mathbb R^3 \setminus \Omega$ form a basis of $H_1(\mathbb R^3 \setminus \Omega;\mathbb Z)$, and then to construct in $\overline{\Omega}$ a homological Seifert surface of each one of these $1$-boundaries. The Poincar\'e--Lefschetz duality theorem ensures that the relative homology classes of these homological Seifert surfaces in $\overline\Omega$ modulo $\partial\Omega$ form a basis of $H_2(\overline\Omega,\partial\Omega;\mathbb Z)$. We devise a simply procedure for the construction of the required set of $1$-boundaries of $\overline{\Omega}$ that, combined with a fast algorithm for the construction of homological Seifert surfaces, allows the efficient computation of a basis of $H_2(\overline{\Omega},\partial\Omega;\mathbb Z)$ via this very natural geometrical approach. Some numerical experiments show the efficiency of the method and its performance comparing with other algorithms., Comment: arXiv admin note: text overlap with arXiv:1409.5487

Numerical simulation of anechoic chambers is a hot topic since it can provide useful data about the performance of the EMC site. However, the mathematical nature of the problem, the physical dimensions of the simulated sites and the frequency ranges pose nontrivial challenges to the simulation. Computational requirements in particular will quickly become unmanageable if adequate techniques are not employed. In this work we describe a novel approach, based on equivalent elements, that enables the simulation of large chambers with modest computational resources. The method is then validated against real measurement results.

Let $\Omega$ be a bounded domain of $\mathbb{R}^3$ whose closure $\overline{\Omega}$ is polyhedral, and let $\mathcal{T}$ be a triangulation of $\overline{\Omega}$. Assuming that the boundary of $\Omega$ is sufficiently regular, we provide an explicit formula for the computation of homological Seifert surfaces of any $1$-boundary $\gamma$ of $\mathcal{T}$; namely, $2$-chains of $\mathcal{T}$ whose boundary is $\gamma$. It is based on the existence of special spanning trees of the complete dual graph of $\mathcal{T}$, and on the computation of certain linking numbers associated with those spanning trees. If the triangulation $\mathcal{T}$ is fine, the explicit formula is too expensive to be used directly. For this reason, making also use of a simple elimination procedure, we devise a fast algorithm for the computation of homological Seifert surfaces. Some numerical experiments illustrate the efficiency of this algorithm., Comment: 30 pages, 14 figures, 2 tables