“…Besides, a key advantage of this scheme is the local solvability, that is, the auxiliary variables approximating the gradient of the solution can be locally eliminated. Over the past twenty years, there has been extensive study of the LDG methods for steady problems or in the semi-discrete framework, such as for elliptic problems [4], convection-diffusion problems [5], the Stokes system [7], the KdV type equations [18], Hamilton-Jacobi equations [17], time-dependent fourth order problems [10], and so on.…”