High accuracy numerical methods for ordinary differential equations with discontinuous right-hand side

Stewart, David (1990). High accuracy numerical methods for ordinary differential equations with discontinuous right-hand side PhD Thesis, School of Physical Sciences. doi:10.14264/uql.2015.244

       
Attached Files (Some files may be inaccessible until you login with your UQ eSpace credentials)
Name Description MIMEType Size Downloads
THE6944.pdf Thesis full text application/pdf 3.61MB 1

Author Stewart, David
Thesis Title High accuracy numerical methods for ordinary differential equations with discontinuous right-hand side
School, Centre or Institute School of Physical Sciences
DOI 10.14264/uql.2015.244
Publication date 1990
Thesis type PhD Thesis
Supervisor Ludvik Bass
Total pages 186
Language eng
Subjects 010109 Ordinary Differential Equations, Difference Equations and Dynamical Systems
Formatted abstract
Ordinary Differential Equations (ODEs) with discontinuous right-hand sides arise in a number of applications in elementary mechanics in relation to dry friction problems; in control theory in relation to switching controls and also in solving optimal control problems; in problems involving elastic-plastic behaviour and other variational inequalities. The crucial difference between these problems and problems with discontinuous coefficients is that the discontinuity is in the "state" variables, rather than occuring in the "time" variable. To emphasise this distinction we consider autonomous problems 𝑥1 = 𝑓(𝑥) where 𝑓 is discontinuous.

The conventional theory regarding existence still works satisfactorily as long as the trajectory 𝑥(𝑡) is at a discontinuity for only finitely many 𝑡. However, there are many examples where the trajectory is forced to stay at the discontinuity. In such a case conventional theory and conventional numerical methods in general break down; in particular, the existence of local solutions cannot be guaranteed (e.g. 𝑥1 = — sgn(𝑥) + α if 0 < |α| < 1). The appropriate formulation of such a problem is to replace the ODE "𝑥1 = 𝑓(𝑥)" by a differential inclusion 𝑥1 ∈ 𝐹(𝑥) where 𝐹 is an appropriate set-valued function with nonempty compact convex values and is upper semi-continuous (in terms of the "⊆" ordering). This formulation is due to Filippov(1960). With this formulation solutions can readily be shown to exist, although they are not necessarily unique.

Numerical methods to compute solutions in the sense of Filippov have already been developed and convergence results proven by Taubert, Elliott and Niepage & Wendt. However, these methods are no better than 1st order accurate, and 1st order accuracy is established only with additional assumptions regarding the number of "switching points". Further, it is only guaranteed that a sub-sequence of the numerical solutions converge, unless solutions are unique. This limits the useful of these methods in solving problems where the desired solution is determined by other considerations, such as in two-point boundary value problems (TPBVPs).

The approach taken here is to assume that the discontinuous ODE has some particular structure, and to exploit that to obtain high order accuracy and other information. The structure assumed is that there is a collection of disjoint open sets with dense union in 𝐑𝑛, 𝑅1,...,𝑅𝑚 such that for 𝑥1 ∈ 𝑅𝑖 implies that 𝑓(𝑥) = 𝑓𝑖(𝑥) with 𝑓𝑖 smooth. Further it is assumed that the 𝑅𝑖 can be described by indicator functions ℎ𝑖: 𝑅𝑖 = {𝑥 | ℎ𝑖(𝑥) < ℎ𝑗(𝑥), 𝑗 ≠ 𝑖}. For each 𝑥 there is the active set 𝐼(𝑥) = {𝑖 | 𝑥 ∈ ̄𝑅𝑖}. The Filippov formulation then becomes

𝑥1 ∈ co{𝑓𝑖(𝑥)| 𝑖 ∈ 𝐼(𝑥)}.

As determining the 𝐼 from 𝑥 is not a well posed problem we suppose that 𝐼(𝑥(𝑡)) is a piecewise constant function. On the intervals on which 𝐼(𝑥(𝑡)) is constant and given, we can select elements of the above differential inclusion that keep the active set constant and moreover gives a smooth ODE to solve. At a switching point the set of possible new values for 𝐼(𝑥(𝑡)) are determined by means of a Linear Complementarity Problem. Given a new active set selected from the list of possibilities, we can continue by solving a new smooth ODE to maintain the new active set, until another switching point occurs. Switching points can be both determined by the ODE itself, or if they are not unique, by the user of the algorithm.

A complete algorithm for solving discontinuous ODEs by this approach has been developed and convergence results proven. These convergence results are of two types: Firstly, given only that the algorithm has produced a sequence of numerical solutions with step size approaching zero, then there is a converging subsequence and every converging subsequence converges to a solution. Secondly, if a particular solution is known which has only finitely many switching points and some non-degeneracy conditions are satisfied, then by suitable selections (in the algorithm) of new active sets a sequence of numerical solutions can be generated that converge to the given solution with the same order of accuracy as that of the smooth ODE solver used.

This algorithm has been implemented and numerical results obtained. These results confirm the expectation of high order accuracy of this method, and that it is competitive with previously published methods.

A number of extensions are proposed for this algorithm. The first is a "de-composition" type of approach that is particularly useful for friction problems with multiple friction surfaces. Convergence results are given for this extension. This extension has been implemented and numerical results obtained for a problem involving three friction surfaces.

Some further extensions are also proposed to avoid some of the non-degeneracy conditions; these are in effect "higher order" modifications to the main algorithm. Also, some difficulties in applying the main algorithm to Hamiltonian inclusions (which are closely related to optimal control theory) are noted.
Keyword Differential equations -- Numerical solutions
Discontinuous functions
Additional Notes Other Title: Numerical methods for ODEs with discontinuous RHS.

Document type: Thesis
Collection: UQ Theses (RHD) - UQ staff and students only
 
Citation counts: Google Scholar Search Google Scholar
Created: Tue, 27 Jan 2015, 16:35:54 EST by Mary-Anne Marrington on behalf of Scholarly Communication and Digitisation Service