  Local discontinuous Galerkin methods are developed for solving second order
and fourth order time-dependent partial differential equations defined on
static 2D manifolds. These schemes are second-order accurate with surfaces
triangulized by planar triangles and careful design of numerical fluxes. The
schemes are proven to be energy stable. Various numerical experiments are
provided to validate the new schemes.
