The why and how of finite elements

The development of the finite element method is traced to show how ideas from structural and fluid mechanics, the calculus of variations, functional analysis and the calculus of finite differences have been forged to provide a tool which minimizes the mismatch between the behaviour of a continuous system and that of a discrete model of the system assembled from finite elements. Geometrical flexibility of the model is achieved by the use of polygonal and curved elements. The behaviour of any point of an element is described in terms of its behaviour at discrete points or nodes of the element. In treating neutron transport the finite element method can be applied to phase-space, or the spatial dependence can be treated by the use of finite elements in conjuction with expansions in orthogonal functions for the directional dependence. The maximum principle for the second-order even-parity Boltzmann equations is used to demonstrate the precision and flexibility of the finite element method by solving the problems of a doglegged duct in a shield and a cylindrical fuel element in a square lattice cell. The geometrical interpretation of the boundary-free maximum principle with the aid of a suitable Hilbert space then leads to completely boundary-free weighted residual or Galerkin schemes for both the first- and second-order forms of the Boltzmann equation. Imposing essential boundary conditions leads to classical schemes, a sketch of finite element treatments of the multigroup Boltzmann equation is given. (author)