Search

The bubble transform is a procedure to decompose differential forms, which are piecewise smooth with respect to a given triangulation of the domain, into a sum of local bubbles. In this paper, an improved version of a construction in the setting of the de Rham complex previously proposed by the authors is presented. The major improvement in the decomposition is that unlike the previous results, in which the individual bubbles were rational functions with the property that groups of local bubbles summed up to preserve piecewise smoothness, the new decomposition is strictly space-preserving in the sense that each local bubble preserves piecewise smoothness. An important property of the transform is that the construction only depends on the given triangulation of the domain and is independent of any finite element space. On the other hand, all the standard piecewise polynomial spaces are invariant under the transform. Other key properties of the transform are that it commutes with the exterior derivative, is bounded in L^2, and satisfies the stable decomposition property.

Solution methods for the nonlinear partial differential equation of the Rudin-Osher-Fatemi (ROF) and minimum-surface models are fundamental for many modern applications. Many efficient algorithms have been proposed. First order methods are common. They are popular due to their simplicity and easy implementation. Some second order Newton-type iterative methods have been proposed like Chan-Golub-Mulet method. In this paper, we propose a new Newton-Krylov solver for primal-dual finite element discretization of the ROF model. The method is so simple that we just need to use some diagonal preconditioners during the iterations. Theoretically, the proposed preconditioners are further proved to be robust and optimal with respect to the mesh size, the penalization parameter, the regularization parameter, and the iterative step, essentially it is a parameter independent preconditioner. We first discretize the primal-dual system by using mixed finite element methods, and then linearize the discrete system by Newton\textquoteright s method. Exploiting the well-posedness of the linearized problem on appropriate Sobolev spaces equipped with proper norms, we propose block diagonal preconditioners for the corresponding system solved with the minimum residual method. Numerical results are presented to support the theoretical results., Comment: 22 pages, 24 figures

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.

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$., Comment: Some typos and other minor errors corrected and additional explanations included