Semilinear mixed problems on Hilbert complexes and their numerical approximation

Arnold, Falk, and Winther recently showed [Bull. Amer. Math. Soc. 47 (2010), 281-354] that linear, mixed variational problems, and their numerical approximation by mixed finite element methods, can be studied using the powerful, abstract language of Hilbert complexes. In another recent article [arXiv:1005.4455], we extended the Arnold-Falk-Winther framework by analyzing variational crimes (a la Strang) on Hilbert complexes. In particular, this gave a treatment of finite element exterior calculus on manifolds, generalizing techniques from surface finite element methods and recovering earlier a priori estimates for the Laplace-Beltrami operator on 2- and 3-surfaces, due to Dziuk [Lecture Notes in Math., vol. 1357 (1988), 142-155] and later Demlow [SIAM J. Numer. Anal., 47 (2009), 805-827], as special cases. In the present article, we extend the Hilbert complex framework in a second distinct direction: to the study of semilinear mixed problems. We do this, first, by introducing an operator-theoretic reformulation of the linear mixed problem, so that the semilinear problem can be expressed as an abstract Hammerstein equation. This allows us to obtain, for semilinear problems, a priori solution estimates and error estimates that reduce to the Arnold-Falk-Winther results in the linear case. We also consider the impact of variational crimes, extending the results of our previous article to these semilinear problems. As an immediate application, this new framework allows for mixed finite element methods to be applied to semilinear problems on surfaces.Comment: 22 pages; v2: major revision, particularly sharpening of error estimates in Section

Convergence and Optimality of Adaptive Mixed Methods on Surfaces

In a 1988 article, Dziuk introduced a nodal finite element method for the Laplace-Beltrami equation on 2-surfaces approximated by a piecewise-linear triangulation, initiating a line of research into surface finite element methods (SFEM). Demlow and Dziuk built on the original results, introducing an adaptive method for problems on 2-surfaces, and Demlow later extended the a priori theory to 3-surfaces and higher order elements. In a separate line of research, the Finite Element Exterior Calculus (FEEC) framework has been developed over the last decade by Arnold, Falk and Winther and others as a way to exploit the observation that mixed variational problems can be posed on a Hilbert complex, and Galerkin-type mixed methods can be obtained by solving finite dimensional subproblems. In 2011, Holst and Stern merged these two lines of research by developing a framework for variational crimes in abstract Hilbert complexes, allowing for application of the FEEC framework to problems that violate the subcomplex assumption of Arnold, Falk and Winther. When applied to Euclidean hypersurfaces, this new framework recovers the original a priori results and extends the theory to problems posed on surfaces of arbitrary dimensions. In yet another seemingly distinct line of research, Holst, Mihalik and Szypowski developed a convergence theory for a specific class of adaptive problems in the FEEC framework. Here, we bring these ideas together, showing convergence and optimality of an adaptive finite element method for the mixed formulation of the Hodge Laplacian on hypersurfaces.Comment: 22 pages, no figures. arXiv admin note: substantial text overlap with arXiv:1306.188

Adaptive Algorithms

Overwhelming empirical evidence in computational science and engineering proved that self-adaptive mesh-generation is a must-do in real-life problem computational partial differential equations. The mathematical understanding of corresponding algorithms concerns the overlap of two traditional mathematical disciplines, numerical analysis and approximation theory, with computational sciences. The half workshop was devoted to the mathematics of optimal convergence rates and instance optimality of the Dörfler marking or the maximum strategy in various versions of space discretisations and time-evolution problems with all kind of applications in the efficient numerical treatment of partial differential equations

3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.Comment: 39 pages, 14 figures. High resolution figures and supplemental movies available upon reques

Topological and Variational Methods for Differential Equations

These notes contain the extended abstracts of the talks presented at the workshop. The range of topics includes nonlinear Schrödinger equations, singularly perturbed equations, symmetry and nodal properties of solutions, long-time dynamics for parabolic equations, Morse theory

Constructing Reference Metrics on Multicube Representations of Arbitrary Manifolds

Reference metrics are used to define the differential structure on multicube representations of manifolds, i.e., they provide a simple and practical way to define what it means globally for tensor fields and their derivatives to be continuous. This paper introduces a general procedure for constructing reference metrics automatically on multicube representations of manifolds with arbitrary topologies. The method is tested here by constructing reference metrics for compact, orientable two-dimensional manifolds with genera between zero and five. These metrics are shown to satisfy the Gauss-Bonnet identity numerically to the level of truncation error (which converges toward zero as the numerical resolution is increased). These reference metrics can be made smoother and more uniform by evolving them with Ricci flow. This smoothing procedure is tested on the two-dimensional reference metrics constructed here. These smoothing evolutions (using volume-normalized Ricci flow with DeTurck gauge fixing) are all shown to produce reference metrics with constant scalar curvatures (at the level of numerical truncation error).Comment: 37 pages, 16 figures; additional introductory material added in version accepted for publicatio