Leaking Bucket Equation and Reversibility

How do reversible integrators work for nonreversible differential equations, such as the leaking bucket equation? This is investigated here for the asynchronous leapfrog integrator, a presumably novel variant of the well-known leapfrog method. For comparison, one can also choose the second-order nonreversible Runge–Kutta method, which is also known as Heun's method.
Both integrators empty the bucket only approximately, even in the long run. So the behavior of the reversed trajectories does not contradict the fact that any one-step method has a flat-line solution, given an autonomous rate function that is initially vanishing.



  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The leaking bucket equation for the water level is , where is the unit step function. The solution is easily seen to be and says that is zero for (because the bucket is empty then). In this Demonstration, we replace the exact solution by a numerical one, which is obtained by repeated application of a so-called integrator (actually the asynchronous leapfrog integrator) to the initial state of the system. We compute this solution well beyond the time that the bucket is empty. Because the integrator under consideration is reversible, we can take the final (empty) state and evolve it back in exact agreement with the forward evolution. Along this reversed trajectory, the bucket will start to get filled exactly at time . How can this be? It works because the state description according to the asynchronous leapfrog method contains, in addition to the water level, a second quantity. In most cases this quantity behaves like . However, when approaches the value 0 on a trajectory of the leaking bucket equation, the steepness of the right-hand side of this equation makes this second quantity oscillate. It is only the mean value over two consecutive steps that resembles . This oscillation together with a small linear trend of carries the information concerning the moment of emptying.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+