In this article we survey r-adaptive (or moving grid) methods for solving time-dependent partial differential equations (PDEs). Although these methods have received much less attention than their h- and p-adaptive counterparts, particularly within the finite element community, we review the substantial progress that has been made in developing more robust and reliable algorithms and in understanding the basic principles behind these methods, and we give some numerical examples illustrative of the wide classes of problems for which these methods are suitable alternatives to the traditional ones.More specifically, we first examine the basic geometric properties of moving meshes in both one and higher spatial dimensions, and discuss the discretization process for PDEs on such moving meshes (both structured and unstructured). In particular, we consider the issues of mesh regularity, equidistribution, alignment, and associated variational methods. An overview is given of the general interpolation error analysis for a function or a truncation error on such an adaptive mesh. Guided by these principles, we show how to design effective moving mesh strategies. We then examine in more detail how these strategies can be implemented in practice. The first class of methods which we consider are based upon controlling mesh density and hence are called position-based methods. These make use of a so-called moving mesh PDE (MMPDE) approach and variational methods, as well as optimal transport methods. This is followed by an analysis of methods which have a more Lagrange-like interpretation, and due to this focus are called velocity-based methods. These include the moving finite element method (MFE), the geometric conservation law (GCL) methods, and the deformation map method. Finally, we present a number of specific types of examples for which the use of a moving mesh method is particularly effective in applications. These include scale-invariant problems, blow-up problems, problems with moving fronts and problems in meteorology. We conclude that, whilst r-adaptive methods are still in their relatively early stages of development, with many outstanding questions remaining, they have enormous potential and indeed can produce an optimal form of adaptivity for many problems.

Abstract.A review is presented of the one-parameter, nonsmooth bifurcations that occur in a variety of continuous-time piecewise-smooth dynamical systems. Motivated by applications, a pragmatic approach is taken to defining a discontinuity-induced bifurcation (DIB) as a nontrivial interaction of a limit set with respect to a codimension-one discontinuity boundary in phase space. Only DIBs that are local are considered, that is, bifurcations involving equilibria or a single point of boundary interaction along a limit cycle for flows. Three classes of systems are considered, involving either state jumps, jumps in the vector field, or jumps in some derivative of the vector field. A rich array of dynamics are revealed, involving the sudden creation or disappearance of attractors, jumps to chaos, bifurcation diagrams with sharp corners, and cascades of period adding. For each kind of bifurcation identified, where possible, a kind of "normal form" or discontinuity mapping (DM) is given, together with a canonical example and an application. The goal is always to explain dynamics that may be observed in simulations of systems which include friction oscillators, impact oscillators, DC-DC converters, and problems in control theory.

In this paper we consider the numerical solution of PDEs with blow-up for which scaling invariance plays a natural role in describing the underlying solution structures. It is a challenging numerical problem to capture the qualitative behaviour in the blow-up region, and the use of nonuniform meshes is essential. We consider moving mesh methods for which the mesh is determined using so-called moving mesh partial differential equations (MMPDEs). Specifically, the underlying PDE and the MMPDE are solved for the blow-up solution and the computational mesh simultaneously. Motivated by the desire for the MMPDE to preserve the scaling invariance of the underlying problem, we study the effect of different choices of MMPDEs and monitor functions. It is shown that for suitable ones the MMPDE solution evolves towards a (moving) mesh which close to the blow-up point automatically places the mesh points in such a manner that the ignition kernel, which is well known to be a natural coordinate in describing the behaviour of blow-up, approaches a constant as-T (the blow-up time). Several numerical examples are given to verify the theory for these MMPDE methods and to illustrate their efficacy.

No abstract

No abstract

No abstract

No abstract

This paper constructs and analyses an adaptive moving mesh scheme for the numerical simulation of singular PDEs in one or more spatial dimensions. The scheme is based on computing a Legendre transformation from a regular to a spatially non-uniform mesh via the solution of a relaxed form of the Monge-Ampere equation. The method is shown to preserve the inherent scaling properties of the PDE and to identify natural computational coordinates. Numerical examples are presented in one and two dimensions which demonstrate the effectiveness of this approach.