An adaptive moving grid method for one-dimensional systems of partial differential equations

Abstract

We describe a fully adaptive, moving grid method for solving initial-boundary value problems for systems of one-space dimensional partial differential equations whose solutions exhibit rapid variations in space and time. The method, based on finite-differences, is of the Lagrangian type and has been derived through a co-ordinate transformation which leads to equidistribution in space of the second derivative. Our technique is "intermediate" between static regridding methods, where nodes remain fixed for intervals of time, and continuously moving grid methods, where the node movement and the PDE integration are fully coupled. In our approach, the computation of the moving grids and the solution on these grids are carried out separately, while the nodes are moved at each time-step. Two error monitors have been implemented, one to govern the time-step selection and the other to eventually adapt the number of moving nodes. The method allows the use of different moving grids for different components in the PDE system. Numerical experiments are presented for a set of five sample problems from the literature, including two problems from combu~tion.,,.,

A numerical study of three moving-grid methods for one-dimensional partial differential equations which are based on the method of lines

“…Observe that the numerical front is ahead of the true one. Concerning the quality of the reference solution we note that in [27] it is claimed that the reference solution is "exact up to plotting accuracy," except perhaps in the neighbourhood of x = 0 at the first output time t = 0.26. All experiments with the present flame problem, including those with Methods II and III, show a deviation here.…”

“…This obviously depends on the nature of the solution sought, which can be nicely illustrated by examining Problem I of Section 4 (cf. [27,Section 5.3] ). Let us consider its solution near the left boundary, while the steep front is forming (the ignition phase).…”

“…Our third example problem is a two-component, semi-linear hyperbolic system, the solution of which is constituted by two waves travelling in opposite directions (copied from Madsen [15], see also [27]). The system is given by Note that these are functions with a mere C 1 continuity, which represent wave pulses located at x = -0.2 and x = 0.2, respectively.…”

An evaluation of the gradient-weighted moving-finite-element method in one space dimension

“…Our third example problem is a two-component, semilinear hyperbolic system, the solution of which is given by two pulses traveling in opposite directions (copied from [10], see also [8,20,21 ]). The system is given by…”

An adaptive moving mesh method for thin film flow equations with surface tension

“…[24] use general purpose publicly available solvers for parabolic PDEs which have built-in adaptive mesh capabilities (PDECOL [28] based on collocation methods; TOMS731 [29][30][31][32] based on a monitor function to adapt mesh and a Lagrangian method for moving the mesh points). Sun et al [33] use a h-adaptive finite element mesh refinement method based on an optimal interpolation error estimate for a two dimensional thin film equation of gravity driven flow down an inclined plane.…”