Marie E. Rognes

Martin S. Alnaes, Anders Logg, Kristian B. Oelgaard et al.

We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies for multi-field problems, general differential operators and flexible tensor algebra. With these features, UFL has been used to effortlessly express finite element methods for complex systems of partial differential equations in near-mathematical notation, resulting in compact, intuitive and readable programs. We present in this work the language and its construction. An implementation of UFL is freely available as an open-source software library. The library generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations. Some application examples are presented and libraries that support UFL are highlighted.

Andre Massing, Mats G. Larson, Anders Logg et al.

We develop a Nitsche fictitious domain method for the Stokes problem starting from a stabilized Galerkin finite element method with low order elements for both the velocity and the pressure. By introducing additional penalty terms for the jumps in the normal velocity and pressure gradients in the vicinity of the boundary, we show that the method is inf-sup stable. As a consequence, optimal order a priori error estimates are established. Moreover, the condition number of the resulting stiffness matrix is shown to be bounded independently of the location of the boundary. We discuss a general, flexible and freely available implementation of the method in three spatial dimensions and present numerical examples supporting the theoretical results.

Jeonghun J. Lee, Eleonora Piersanti, Kent-Andre Mardal et al.

In this paper, we present and analyze a new mixed finite element formulation of a general family of quasi-static multiple-network poroelasticity (MPET) equations. The MPET equations describe flow and deformation in an elastic porous medium that is permeated by multiple fluid networks of differing characteristics. As such, the MPET equations represent a generalization of Biot's equations, and numerical discretizations of the MPET equations face similar challenges. Here, we focus on the nearly incompressible case for which standard mixed finite element discretizations of the MPET equations perform poorly. Instead, we propose a new mixed finite element formulation based on introducing an additional total pressure variable. By presenting energy estimates for the continuous solutions and a priori error estimates for a family of compatible semi-discretizations, we show that this formulation is robust in the limits of incompressibility, vanishing storage coefficients, and vanishing transfer between networks. These theoretical results are corroborated by numerical experiments. Our primary interest in the MPET equations stems from the use of these equations in modelling interactions between biological fluids and tissues in physiological settings. So, we additionally present physiologically realistic numerical results for blood and tissue fluid flow interactions in the human brain.

André Massing, Mats G. Larson, Anders Logg et al.

We develop a Nitsche-based formulation for a general class of stabilized finite element methods for the Stokes problem posed on a pair of overlapping, non-matching meshes. By ex- tending the least-squares stabilization to the overlap region, we prove that the method is stable, consistent, and optimally convergent. To avoid an ill-conditioned linear algebra system, the scheme is augmented by a least-squares term measuring the discontinuity of the solution in the overlap region of the two meshes. As a consequence, we may prove an estimate for the condition number of the resulting stiffness matrix that is independent of the location of the interface. Finally, we present numerical examples in three spatial dimensions illustrating and confirming the theoretical results.

Marie E. Rognes, Anders Logg

This article presents a general and novel approach to the automation of goal-oriented error control in the solution of nonlinear stationary finite element variational problems. The approach is based on automated linearization to obtain the linearized dual problem, automated derivation and evaluation of a posteriori error estimates, and automated adaptive mesh refinement to control the error in a given goal functional to within a given tolerance. Numerical examples representing a variety of different discretizations of linear and nonlinear partial differential equations are presented, including Poisson's equation, a mixed formulation of linear elasticity, and the incompressible Navier-Stokes equations.

Matteo Croci, Michael B. Giles, Marie E. Rognes et al.

When solving stochastic partial differential equations (SPDEs) driven by additive spatial white noise, the efficient sampling of white noise realizations can be challenging. Here, we present a new sampling technique that can be used to efficiently compute white noise samples in a finite element method and multilevel Monte Carlo (MLMC) setting. The key idea is to exploit the finite element matrix assembly procedure and factorize each local mass matrix independently, hence avoiding the factorization of a large matrix. Moreover, in a MLMC framework, the white noise samples must be coupled between subsequent levels. We show how our technique can be used to enforce this coupling even in the case of non-nested mesh hierarchies. We demonstrate the efficacy of our method with numerical experiments. We observe optimal convergence rates for the finite element solution of the elliptic SPDEs of interest in 2D and 3D and we show convergence of the sampled field covariances. In a MLMC setting, a good coupling is enforced and the telescoping sum is respected.

Patrick E. Farrell, Johan E. Hake, Simon W. Funke et al.

Mathematical models that couple partial differential equations (PDEs) and spatially distributed ordinary differential equations (ODEs) arise in biology, medicine, chemistry and many other fields. In this paper we discuss an extension to the FEniCS finite element software for expressing and efficiently solving such coupled systems. Given an ODE described using an augmentation of the Unified Form Language (UFL) and a discretisation described by an arbitrary Butcher tableau, efficient code is automatically generated for the parallel solution of the ODE. The high-level description of the solution algorithm also facilitates the automatic derivation of the adjoint and tangent linearization of coupled PDE-ODE solvers. We demonstrate the capabilities of the approach on examples from cardiac electrophysiology and mitochondrial swelling.

Andre Massing, Mats G. Larson, Anders Logg et al.

We present a new composite mesh finite element method for fluid--structure interaction problems. The method is based on surrounding the structure by a boundary-fitted fluid mesh which is embedded into a fixed background fluid mesh. The embedding allows for an arbitrary overlap of the fluid meshes. The coupling between the embedded and background fluid meshes is enforced using a stabilized Nitsche formulation which allows us to establish stability and optimal order \emph{a priori} error estimates, see~\cite{MassingLarsonLoggEtAl2013}. We consider here a steady state fluid--structure interaction problem where a hyperelastic structure interacts with a viscous fluid modeled by the Stokes equations. We evaluate an iterative solution procedure based on splitting and present three-dimensional numerical examples.

Douglas N. Arnold, Marie E. Rognes

The stability properties of simple element choices for the mixed formulation of the Laplacian are investigated numerically. The element choices studied use vector Lagrange elements, i.e., the space of continuous piecewise polynomial vector fields of degree at most r, for the vector variable, and the divergence of this space, which consists of discontinuous piecewise polynomials of one degree lower, for the scalar variable. For polynomial degrees r equal 2 or 3, this pair of spaces was found to be stable for all mesh families tested. In particular, it is stable on diagonal mesh families, in contrast to its behaviour for the Stokes equations. For degree r equal 1, stability holds for some meshes, but not for others. Additionally, convergence was observed precisely for the methods that were observed to be stable. However, it seems that optimal order L2 estimates for the vector variable, known to hold for r>3, do not hold for lower degrees.