Aarhus University Seal / Aarhus Universitets segl

Dieter Britz

Damping of Crank–Nicolson error oscillations

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review

  • Department of Chemistry
  • Department of Computer Science
The Crank–Nicolson (CN) simulation method has an oscillatory response to sharp initial transients. The technique is convenient but the oscillations make it less popular. Several ways of damping the oscillations in two types of electrochemical computations are investigated. For a simple one-dimensional system with an initial singularity, subdivision of the first time interval into a number of equal subintervals (the Pearson method) works rather well, and so does division with exponentially increasing subintervals, where however an optimum expansion parameter must be found. This method can be computationally more expensive with some systems. The simple device of starting with one backward implicit (BI, or Laasonen) step does damp the oscillations, but not always sufficiently. For electrochemical microdisk simulations which are two-dimensional in space and using CN, the use of a first BI step is much more effective and is recommended. Division into subintervals is also effective, and again, both the Pearson method and exponentially increasing subintervals methods are effective here. Exponentially increasing subintervals are often considerably more expensive computationally. Expanding intervals over the whole simulation period, although capable of satisfactory results, for most systems will require more cpu time compared with subdivision of the first interval only.
Original languageEnglish
JournalComputational Biology and Chemistry
Pages (from-to)253-263
Number of pages11
Publication statusPublished - 2003

See relations at Aarhus University Citationformats

ID: 36380592