HOME    »    SCIENTIFIC RESOURCES    »    Volumes
Abstracts and Talk Materials
Finite Element Circus featuring a Scientific Celebration of Falk, Pasciak, and Wahlbin
November 5-6, 2010

Douglas N. Arnold

Plenary talk: Canonical families of finite elements
November 6, 2010

The most familiar family of finite elements is the Lagrange family, which provide the canonical finite element approximation of H1 on simplicial meshes in any dimension. In this talk we discuss families of simplicial and cubical finite elements—some previously known and some new—which are natural extensions of the Lagrange family in various ways. Even for some of the long known elements, a modern viewpoint based on the finite element exterior calculus provides new properties and insights.

James H. Bramble

Plenary talk: Analysis of a Cartesian PML approximation to an acoustic scattering problem
November 6, 2010

We consider a Cartesian PML approximation to solutions of acoustic scattering problems on an unbounded domain in ℝ2 and ℝ3. The perfectly matched layer (PML) technique modifies the equations outside of a bounded domain containing the region of interest. This is done in such a way that the new problem (still on an unbounded domain) has a solution which agrees with the solution of the original problem. The new problem has a solution which decays much faster, thus suggesting replacing it by a problem on a bounded domain. The perfectly matched layer (PML) technique, in a curvilinear coordinate system and in Cartesian coordinates, has been studied for acoustic scattering applications both in theory and computation. Using a different approach we extend the results of Kim and Pasciak concerning the PML technique in Cartesian coordinates. The exponential convergence of approximate solutions as a function of domain size and/or the PML "strength" parameter, σ0 is also shown. We note that once the stability and convergence of the (continuous) truncated problem has been established, the analysis of the resulting finite element approximations is then classical. Finally, the results of numerical computations illustrating the theory, in terms of efficiency and parameter dependence of the Cartesian PML approach will be given.

Vidar Thomée

Plenary talk: On the lumped mass finite element method for parabolic problems
November 5, 2010

We study the lumped mass method for the model homogeneous heat equation with homogeneous Dirichlet boundary conditions. We first recall that the maximum principle for the heat equation does not carry over to the the spatially semidiscrete standard Galerkin finite element method, using continuous, piecewise linear approximating functions. However, for the lumped mass variant the situation is more advantageous. We present necessary and sufficient conditions on the triangulation, expressed in terms of properties of the stiffness matrix, for the semidiscrete lumped mass solution operator to be a positive operator or a contraction in the maximum-norm.

We then turn to error estimates in the L2-norm. Improving earlier results we show that known optimal order smooth initial data error estimates for the standard Galerkin method carry over to the lumped mass method, whereas nonsmooth initial data estimates require special assumptions on the triangulations.

We also discuss the application to time discretization by the backward Euler and Crank-Nicolson methods.