Gantumur Tsogtgerel

Douglas N. Arnold, Richard S. Falk, Johnny Guzmán et al.

In 1976, Dodziuk and Patodi employed Whitney forms to define a combinatorial codifferential operator on cochains, and they raised the question whether it is consistent in the sense that for a smooth enough differential form the combinatorial codifferential of the associated cochain converges to the exterior codifferential of the form as the triangulation is refined. In 1991, Smits proved this to be the case for the combinatorial codifferential applied to 1-forms in two dimensions under the additional assumption that the initial triangulation is refined in a completely regular fashion, by dividing each triangle into four similar triangles. In this paper we extend Smits's result to arbitrary dimensions, showing that the combinatorial codifferential on 1-forms is consistent if the triangulations are uniform or piecewise uniform in a certain precise sense. We also show that this restriction on the triangulations is needed, giving a counterexample in which a different regular refinement procedure, namely Whitney's standard subdivision, is used. Further, we show by numerical example that for 2-forms in three dimensions, the combinatorial codifferential is not consistent even for the most regular subdivision process.

Erick Schulz, Gantumur Tsogtgerel

Discrete exterior calculus (DEC) is a framework for constructing discrete versions of exterior differential calculus objects, and is widely used in computer graphics, computational topology, and discretizations of the Hodge-Laplace operator and other related partial differential equations. However, a rigorous convergence analysis of DEC has always been lacking; as far as we are aware, the only convergence proof of DEC so far appeared is for the scalar Poisson problem in two dimensions, and it is based on reinterpreting the discretization as a finite element method. Moreover, even in two dimensions, there have been some puzzling numerical experiments reported in the literature, apparently suggesting that there is convergence without consistency. In this paper, we develop a general independent framework for analyzing issues such as convergence of DEC without relying on theories of other discretization methods, and demonstrate its usefulness by establishing convergence results for DEC beyond the Poisson problem in two dimensions. Namely, we prove that DEC solutions to the scalar Poisson problem in arbitrary dimensions converge pointwise to the exact solution at least linearly with respect to the mesh size. We illustrate the findings by various numerical experiments, which show that the convergence is in fact of second order when the solution is sufficiently regular. The problems of explaining the second order convergence, and of proving convergence for general p-forms remain open.

Michael Holst, Gantumur Tsogtgerel, Yunrong Zhu

In this article we develop convergence theory for a general class of adaptive approximation algorithms for abstract nonlinear operator equations on Banach spaces, and use the theory to obtain convergence results for practical adaptive finite element methods (AFEM) applied to several classes of nonlinear elliptic equations. In the first part of the paper, we develop a weak-* convergence framework for nonlinear operators, whose Gateaux derivatives are locally Lipschitz and satisfy a local inf-sup condition. The framework can be viewed as extending the recent convergence results for linear problems of Morin, Siebert and Veeser to a general nonlinear setting. We formulate an abstract adaptive approximation algorithm for nonlinear operator equations in Banach spaces with local structure. The weak-* convergence framework is then applied to this class of abstract locally adaptive algorithms, giving a general convergence result. The convergence result is then applied to a standard AFEM algorithm in the case of several semilinear and quasi-linear scalar elliptic equations and elliptic systems, including: a semilinear problem with subcritical nonlinearity, the steady Navier-Stokes equations, and a quasilinear problem with nonlinear diffusion. This yields several new AFEM convergence results for these nonlinear problems. In the second part of the paper we develop a second abstract convergence framework based on strong contraction, extending the recent contraction results for linear problems of Cascon, Kreuzer, Nochetto, and Siebert and of Mekchay and Nochetto to abstract nonlinear problems. The contraction result is then applied to a standard AFEM algorithm for semilinear problems with sub- and super-critical nonlinearities and for the Hamiltonian constraint in general relativity.