Aqueous radiolysis

examples/aqueous_radiolysis.py is an example of a rather large system of reactions.

$ python analytic_diffusion.py --help

usage: aqueous_radiolysis.py [-h] [--t0 T0] [--tend TEND] [--x0 X0]
                             [--xend XEND] [-d DOSERATE] [-N N] [--nt NT]
                             [--nstencil NSTENCIL] [--logy] [--logt] [--logx]
                             [--name NAME] [--num-jacobian] [-s SAVEFIG] [-v]
                             [--plot] [--plot-jacobians]

    Integrates the reaction system defined by
    :download:`aqueous_radiolysis.json <examples/aqueous_radiolysis.json>`


optional arguments:
  -h, --help            show this help message and exit
  --t0 T0               1e-07
  --tend TEND           0.1
  --x0 X0               1e-09
  --xend XEND           0.1
  -d DOSERATE, --doserate DOSERATE
                        15
  -N N, --N N           1000
  --nt NT               512
  --nstencil NSTENCIL   3
  --logy                False
  --logt                False
  --logx                False
  --name NAME           u'aqueous_radiolysis'
  --num-jacobian        False
  -s SAVEFIG, --savefig SAVEFIG
                        u'None'
  -v, --verbose         False
  --plot                False
  --plot-jacobians      False

Here is an example generated by:

$ python analytic_diffusion.py --doserate 25 --plot     --savefig analytic_diffusion.png
../_images/aqueous_radiolysis.png
aqueous_radiolysis.integrate_rd(t0=1e-07, tend=0.1, x0=1e-09, xend=0.1, doserate=15, N=1000, nt=512, nstencil=3, logy=False, logt=False, logx=False, name=u'aqueous_radiolysis', num_jacobian=False, savefig=u'None', verbose=False, plot=False, plot_jacobians=False)[source]

Integrates the reaction system defined by aqueous_radiolysis.json