Add DG discretization to Crystallization module

The crystallization PDE equations are currently discretized using higher order FV schemes, #112 and #241 .

The computational performance for simulating these models is very relevant since large systems and many terms need to be solved.

We should explore suitable higher/arbitrary order methods for the crystallization equations, such as DG.