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

Ragnar Winther

Ragnar Winther 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.APmath.DGmath.GT

Catalog footprint

What is connected

14works
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

14 published item(s)

preprint2022arXiv

A uniform preconditioner for a Newton algorithm for total-variation minimization and minimum-surface problems

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.

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Well-conditioned frames for high order finite element methods

The purpose of this paper is to discuss representations of high order $C^0$ finite element spaces on simplicial meshes in any dimension. When computing with high order piecewise polynomials the conditioning of the basis is likely to be important. The main result of this paper is a construction of representations by frames such that the associated $L^2$ condition number is bounded independently of the polynomial degree. To our knowledge, such a representation has not been presented earlier. The main tools we will use for the construction is the bubble transform, introduced previously in [Falk and Winther, Found Comput Math (2016) 16: 297], and properties of Jacobi polynomials on simplexes in higher dimensions. We also include a brief discussion of preconditioned iterative methods for the finite element systems in the setting of representations by frames.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Parameter-robust discretization and preconditioning of Biot's consolidation model

Biot's consolidation model in poroelasticity has a number of applications in science, medicine, and engineering. The model depends on various parameters, and in practical applications these parameters ranges over several orders of magnitude. A current challenge is to design discretization techniques and solution algorithms that are well behaved with respect to these variations. The purpose of this paper is to study finite element discretizations of this model and construct block diagonal preconditioners for the discrete Biot systems. The approach taken here is to consider the stability of the problem in non-standard or weighted Hilbert spaces and employ the operator preconditioning approach. We derive preconditioners that are robust with respect to both the variations of the parameters and the mesh refinement. The parameters of interest are small time-step sizes, large bulk and shear moduli, and small hydraulic conductivity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonconforming tetrahedral mixed finite elements for elasticity

This paper presents a nonconforming finite element approximation of the space of symmetric tensors with square integrable divergence, on tetrahedral meshes. Used for stress approximation together with the full space of piecewise linear vector fields for displacement, this gives a stable mixed finite element method which is shown to be linearly convergent for both the stress and displacement, and which is significantly simpler than any stable conforming mixed finite element method. The method may be viewed as the three-dimensional analogue of a previously developed element in two dimensions. As in that case, a variant of the method is proposed as well, in which the displacement approximation is reduced to piecewise rigid motions and the stress space is reduced accordingly, but the linear convergence is retained.

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A uniform inf--sup condition with applications to preconditioning

A uniform inf-sup condition related to a parameter dependent Stokes problem is established. Such conditions are intimately connected to the construction of uniform preconditioners for the problem, i.e., preconditioners which behave uniformly well with respect to variations in the model parameter as well as the discretization parameter. For the present model, similar results have been derived before, but only by utilizing extra regularity ensured by convexity of the domain. The purpose of this paper is to remove this artificial assumption. As a byproduct of our analysis, in the two dimensional case we also construct a new projection operator for the Taylor-Hood element which is uniformly bounded in $L^2$ and commutes with the divergence operator. This construction is based on a tight connection between a subspace of the Taylor-Hood velocity space and the lowest order Nedelec edge element.

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On variational eigenvalue approximation of semidefinite operators

Eigenvalue problems for semidefinite operators with infinite dimensional kernels appear for instance in electromagnetics. Variational discretizations with edge elements have long been analyzed in terms of a discrete compactness property. As an alternative, we show here how the abstract theory can be developed in terms of a geometric property called the vanishing gap condition. This condition is shown to be equivalent to eigenvalue convergence and intermediate between two different discrete variants of Friedrichs estimates. Next we turn to a more practical means of checking these properties. We introduce a notion of compatible operator and show how the previous conditions are equivalent to the existence of such operators with various convergence properties. In particular the vanishing gap condition is shown to be equivalent to the existence of compatible operators satisfying an Aubin-Nitsche estimate. Finally we give examples demonstrating that the implications not shown to be equivalences, indeed are not.

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Finite elements for symmetric tensors in three dimension

We construct finite element subspaces of the space of symmetric tensors with square-integrable divergence on a three-dimensional domain. These spaces can be used to approximate the stress field in the classical Hellinger--Reissner mixed formulation of the elasticty equations, when standard discontinous finite element spaces are used to approximate the displacement field. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and there is one for each positive value of the polynomial degree used for the displacements. For each degree, these provide a stable finite element discretization. The construction of the spaces is closely tied to discretizations of the elasticity complex, and can be viewed as the three-dimensional analogue of the triangular element family for plane elasticity previously proposed by Arnold and Winther.

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

0 citations0 reviewsNo code linkedOpen access