Source author record

Richard S. Falk

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

ResearcherUnclaimed source record

math.NA Numerical Analysis math.NT math.DG math.GT

Catalog footprint

What is connected

10works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Construction of polynomial preserving cochain extensions by blending

A classical technique to construct polynomial preserving extensions of scalar functions defined on the boundary of an $n$ simplex to the interior is to use so-called rational blending functions. The purpose of this paper is to generalize the construction by blending to the de Rham complex. More precisely, we define polynomial preserving extensions which map traces of $k$ forms defined on the boundary of the simplex to $k$ forms defined in the interior. Furthermore, the extensions are cochain maps, i.e., they commute with the exterior derivative.

preprint2022arXiv

The Bubble Transform and the de Rham Complex

The purpose of this paper is to discuss a generalization of the bubble transform to differential forms. The bubble transform was discussed in a previous paper by the authors for scalar valued functions, or zero-forms, and represents a new tool for the understanding of finite element spaces of arbitrary polynomial degree. The present paper contains a similar study for differential forms. From a simplicial mesh of the domain, we build a map which decomposes piecewise smooth $k$-forms into a sum of local bubbles supported on appropriate macroelements. The key properties of the decomposition are that it commutes with the exterior derivative and preserves the piecewise polynomial structure of the standard finite element spaces of $k$-forms. Furthermore, the transform is bounded in $L^2$ and also on the appropriate subspace consisting of $k$-forms with exterior derivatives in $L^2$.

preprint2021arXiv

Hidden Positivity and a New Approach to Numerical Computation of Hausdorff Dimension: Higher Order Methods

In [14], the authors developed a new approach to the computation of the Hausdorff dimension of the invariant set of an iterated function system or IFS. In this paper, we extend this approach to incorporate high order approximation methods. We again rely on the fact that we can associate to the IFS a parametrized family of positive, linear, Perron-Frobenius operators $L_s$, an idea known in varying degrees of generality for many years. Although $L_s$ is not compact in the setting we consider, it possesses a strictly positive $C^m$ eigenfunction $v_s$ with eigenvalue $R(L_s)$ for arbitrary $m$ and all other points $z$ in the spectrum of $L_s$ satisfy $|z| \le b$ for some constant $b < R(L_s)$. Under appropriate assumptions on the IFS, the Hausdorff dimension of the invariant set of the IFS is the value $s=s_*$ for which $R(L_s) =1$. This eigenvalue problem is then approximated by a collocation method at the extended Chebyshev points of each subinterval using continuous piecewise polynomials of arbitrary degree $r$. Using an extension of the Perron theory of positive matrices to matrices that map a cone $K$ to its interior and explicit a priori bounds on the derivatives of the strictly positive eigenfunction $v_s$, we give rigorous upper and lower bounds for the Hausdorff dimension $s_*$, and these bounds converge rapidly to $s_*$ as the mesh size decreases and/or the polynomial degree increases.

preprint2016arXiv

C^m Eigenfunctions of Perron-Frobenius Operators and a New Approach to Numerical Computation of Hausdorff Dimension

We develop a new approach to the computation of the Hausdorff dimension of the invariant set of an iterated function system or IFS. In the one dimensional case, our methods require only C^3 regularity of the maps in the IFS. The key idea, which has been known in varying degrees of generality for many years, is to associate to the IFS a parametrized family of positive, linear, Perron-Frobenius operators L_s. The operators L_s can typically be studied in many different Banach spaces. Here, unlike most of the literature, we study L_s in a Banach space of real-valued, C^k functions, k >= 2; and we note that L_s is not compact, but has a strictly positive eigenfunction v_s with positive eigenvalue lambda_s equal to the spectral radius of L_s. Under appropriate assumptions on the IFS, the Hausdorff dimension of the invariant set of the IFS is the value s=s_* for which lambda_s =1. This eigenvalue problem is then approximated by a collocation method using continuous piecewise linear functions (in one dimension) or bilinear functions (in two dimensions). Using the theory of positive linear operators and explicit a priori bounds on the derivatives of the strictly positive eigenfunction v_s, we give rigorous upper and lower bounds for the Hausdorff dimension s_*, and these bounds converge to s_* as the mesh size approaches zero.

preprint2013arXiv

The bubble transform: A new tool for analysis of finite element methods

The purpose of this paper is to discuss the construction of a linear operator, referred to as the bubble transform, which maps scalar functions defined on a bounded domain $Ω$ in $\mathbb{R}^n$ into a collection of functions with local support. In fact, for a given simplicial triangulation of $Ω$, the associated bubble transform produces a decomposition of functions on $Ω$ into a sum of functions with support on the corresponding macroelements. The transform is bounded in both $L^2$ and the Sobolev space $H^1$, it is local, and it preserves the corresponding continuous piecewise polynomial spaces. As a consequence, this transform is a useful tool for constructing local projection operators into finite element spaces such that the appropriate operator norms are bounded independently of polynomial degree. The transform is basically constructed by two families of operators, local averaging operators and rational trace preserving cut--off operators.

preprint2012arXiv

Local bounded cochain projection

We construct projections from the space of differential k-forms which belong to L2 and whose exterior derivative also belongs to L2, to finite dimensional subspaces of piecewise polynomial differential forms defined on a simplicial mesh. These projections have the properties that they commute with the exterior derivative and are bounded independent of the mesh size. Unlike some other recent work in this direction, the projections are also locally defined in the sense that they are defined by local operators on overlapping macroelements, in the spirit of the Clement interpolant.

preprint2011arXiv

Mixed finite element approximation of the vector Laplacian with Dirichlet boundary conditions

We consider the finite element solution of the vector Laplace equation on a domain in two dimensions. For various choices of boundary conditions, it is known that a mixed finite element method, in which the rotation of the solution is introduced as a second unknown, is advantageous, and appropriate choices of mixed finite element spaces lead to a stable, optimally convergent discretization. However, the theory that leads to these conclusions does not apply to the case of Dirichlet boundary conditions, in which both components of the solution vanish on the boundary. We show, by computational example, that indeed such mixed finite elements do not perform optimally in this case, and we analyze the suboptimal convergence that does occur. As we indicate, these results have implications for the solution of the biharmonic equation and of the Stokes equations using a mixed formulation involving the vorticity.

preprint2009arXiv

Finite element exterior calculus: from Hodge theory to numerical stability

This article reports on the confluence of two streams of research, one emanating from the fields of numerical analysis and scientific computation, the other from topology and geometry. In it we consider the numerical discretization of partial differential equations that are related to differential complexes so that de Rham cohomology and Hodge theory are key tools for the continuous problem. After a brief introduction to finite element methods, the discretization methods we consider, we develop an abstract Hilbert space framework for analyzing stability and convergence. In this framework, the differential complex is represented by a complex of Hilbert spaces and stability is obtained by transferring Hodge theoretic structures from the continuous level to the discrete. We show stable discretization is achieved if the finite element spaces satisfy two hypotheses: they form a subcomplex and there exists a bounded cochain projection from the full complex to the subcomplex. Next, we consider the most canonical example of the abstract theory, in which the Hilbert complex is the de Rham complex of a domain in Euclidean space. We use the Koszul complex to construct two families of finite element differential forms, show that these can be arranged in subcomplexes of the de Rham complex in numerous ways, and for each construct a bounded cochain projection. The abstract theory therefore applies to give the stability and convergence of finite element approximations of the Hodge Laplacian. Other applications are considered as well, especially to the equations of elasticity. Background material is included to make the presentation self-contained for a variety of readers.

preprint2008arXiv

Geometric decompositions and local bases for spaces of finite element differential forms

We study the two primary families of spaces of finite element differential forms with respect to a simplicial mesh in any number of space dimensions. These spaces are generalizations of the classical finite element spaces for vector fields, frequently referred to as Raviart-Thomas, Brezzi-Douglas-Marini, and Nedelec spaces. In the present paper, we derive geometric decompositions of these spaces which lead directly to explicit local bases for them, generalizing the Bernstein basis for ordinary Lagrange finite elements. The approach applies to both families of finite element spaces, for arbitrary polynomial degree, arbitrary order of the differential forms, and an arbitrary simplicial triangulation in any number of space dimensions. A prominent role in the construction is played by the notion of a consistent family of extension operators, which expresses in an abstract framework a sufficient condition for deriving a geometric decomposition of a finite element space leading to a local basis.

preprint2007arXiv

Mixed finite element methods for linear elasticity with weakly imposed symmetry

In this paper, we construct new finite element methods for the approximation of the equations of linear elasticity in three space dimensions that produce direct approximations to both stresses and displacements. The methods are based on a modified form of the Hellinger--Reissner variational principle that only weakly imposes the symmetry condition on the stresses. Although this approach has been previously used by a number of authors, a key new ingredient here is a constructive derivation of the elasticity complex starting from the de Rham complex. By mimicking this construction in the discrete case, we derive new mixed finite elements for elasticity in a systematic manner from known discretizations of the de Rham complex. These elements appear to be simpler than the ones previously derived. For example, we construct stable discretizations which use only piecewise linear elements to approximate the stress field and piecewise constant functions to approximate the displacement field.