discretization.md

Theseus employs a nodal DGSEM formulation on tensor-product elements for spatial discretization [4, 5]. The solution strategy follows an entropy-conservative BR1 form [6], and uses explicit $s$-stage Runge-Kutta methods for time advancement.