C01 - Lithium-Cluster Dynamics

The Benchmark C01 Lithium Cluster Dynamics, the first benchmark in SNE, is taken from solid state physics. A balance model describes formation and decay of defect aggregates in alkali halides, produced by electron bombardment. The nonlinear ODE model consists of three states (defect aggregates) and is very stiff.

Tasks are simulation of the stiff system comparing ODE solvers, a parameter study for a critical decay parameter, and steady state calculation. Emphasis of educational aspects is simulation of stiff systems, and steady state analysis.