Internally discontinuous finite elements for moving interface problems

All physical process change with time and although physicists and engineers often freeze them to simplify their model studies, this is not always possible. Add to this the existence of several subregions within the modelling domain with different physical properties, the size and number of which change with time, and one has produced what might appear to be an intractable problem. Previous studies have shown that adaptive grids are an inconvenient way of modelling such moving interface problems and suggested that the use of a fixed grid with special discontinuous elements affords an efficient solution technique. This work presents and develops the characteristic matrices for two such elements, namely a quadratic triangle with an internal interface modelled by two straight lines and a quadratic isoparametric element. Both these elements are tested against an analytical solution to a simple Poisson equation. Such tests reveal the performance of the proposed elements to be satisfactory except for the definition of gradients near the interface.