Finite element methods for surface PDEs

In this article we consider finite element methods for approximating the solution of partial differential equations on surfaces. We focus on surface finite elements on triangulated surfaces, implicit surface methods using level set descriptions of the surface, unfitted finite element methods and diffuse interface methods. In order to formulate the methods we present the necessary geometric analysis and, in the context of evolving surfaces, the necessary transport formulae. A wide variety of equations and applications are covered. Some ideas of the numerical analysis are presented along with illustrative numerical examples

Exhibiting cross-diffusion-induced patterns for reaction-diffusion systems on evolving domains and surfaces

The aim of this manuscript is to present for the first time the application of the finite element method for solving reaction-diffusion systems with cross-diffusion on continuously evolving domains and surfaces. Furthermore we present pattern formation generated by the reaction-diffusion systemwith cross-diffusion on evolving domains and surfaces. A two-component reaction-diffusion system with linear cross-diffusion in both u and v is presented. The finite element method is based on the approximation of the domain or surface by a triangulated domain or surface consisting of a union of triangles. For surfaces, the vertices of the triangulation lie on the continuous surface. A finite element space of functions is then defined by taking the continuous functions which are linear affine on each simplex of the triangulated domain or surface. To demonstrate the role of cross-diffusion to the theory of pattern formation, we compute patterns with model kinetic parameter values that belong only to the cross-diffusion parameter space; these do not belong to the standard parameter space for classical reaction-diffusion systems. Numerical results exhibited show the robustness, flexibility, versatility, and generality of our methodology; the methodology can deal with complicated evolution laws of the domain and surface, and these include uniform isotropic and anisotropic growth profiles as well as those profiles driven by chemical concentrations residing in the domain or on the surface

Finite element error analysis of wave equations with dynamic boundary conditions: L2L^2 estimates

$L^2$ norm error estimates of semi- and full discretisations, using bulk--surface finite elements and Runge--Kutta methods, of wave equations with dynamic boundary conditions are studied. The analysis resides on an abstract formulation and error estimates, via energy techniques, within this abstract setting. Four prototypical linear wave equations with dynamic boundary conditions are analysed which fit into the abstract framework. For problems with velocity terms, or with acoustic boundary conditions we prove surprising results: for such problems the spatial convergence order is shown to be less than two. These can also be observed in the presented numerical experiments

Error analysis for full discretizations of quasilinear parabolic problems on evolving surfaces

Convergence results are shown for full discretizations of quasilinear parabolic partial differential equations on evolving surfaces. As a semidiscretization in space the evolving surface finite element method is considered, using a regularity result of a generalized Ritz map, optimal order error estimates for the spatial discretization is shown. Combining this with the stability results for Runge--Kutta and BDF time integrators, we obtain convergence results for the fully discrete problems.Comment: -. arXiv admin note: text overlap with arXiv:1410.048

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

We present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C1-continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. The interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching

Quantization of systems with temporally varying discretization I: Evolving Hilbert spaces

A temporally varying discretization often features in discrete gravitational systems and appears in lattice field theory models subject to a coarse graining or refining dynamics. To better understand such discretization changing dynamics in the quantum theory, an according formalism for constrained variational discrete systems is constructed. While the present manuscript focuses on global evolution moves and, for simplicity, restricts to Euclidean configuration spaces, a companion article discusses local evolution moves. In order to link the covariant and canonical picture, the dynamics of the quantum states is generated by propagators which satisfy the canonical constraints and are constructed using the action and group averaging projectors. This projector formalism offers a systematic method for tracing and regularizing divergences in the resulting state sums. Non-trivial coarse graining evolution moves lead to non-unitary, and thus irreversible, projections of physical Hilbert spaces and Dirac observables such that these concepts become evolution move dependent on temporally varying discretizations. The formalism is illustrated in a toy model mimicking a `creation from nothing'. Subtleties arising when applying such a formalism to quantum gravity models are discussed.Comment: 45 pages, 1 appendix, 6 figures (additional explanations, now matches published version

Vistas in numerical relativity

Upcoming gravitational wave-experiments promise a window for discovering new physics in astronomy. Detection sensitivity of the broadband laser interferometric detectors LIGO/VIRGO may be enhanced by matched filtering with accurate wave-form templates. Where analytic methods break down, we have to resort to numerical relativity, often in Hamiltonian or various hyperbolic formulations. Well-posed numerical relativity requires consistency with the elliptic constraints of energy and momentum conservation. We explore this using a choice of gauge in the future and a dynamical gauge in the past. Applied to a polarized Gowdy wave, this enables solving {\em all} ten vacuum Einstein equations. Evolution of the Schwarzschild metric in 3+1 and, more generally, sufficient conditions for well-posed numerical relativity continue to be open challenges.Comment: invited talk, Asian Pacific CTP Winter School on black hole astrophysics, Pohang, Kore