Local Discontinuous Galerkin Methods for One-Dimensional Second Order Fully Nonlinear Elliptic and Parabolic Equations

• Authors: Xiaobing Feng;Thomas Lewis
• Venue: Journal of Scientific Computing
• Year: 2014

This paper is concerned with developing accurate and efficient nonstandard discontinuous Galerkin methods for fully nonlinear second order elliptic and parabolic partial differential equations (PDEs) in the case of one spatial dimension. The primary goal of the paper to develop a general framework for constructing high order local discontinuous Galerkin (LDG) methods for approximating viscosity solutions of these fully nonlinear PDEs which are merely continuous functions by definition. In order to capture discontinuities of the first order derivative $$u_x$$ux of the solution $$u$$u, two independent functions $$q^-$$q- and $$q^+$$q+ are introduced to approximate one-sided derivatives of $$u$$u. Similarly, to capture the discontinuities of the second order derivative $$u_{xx}$$uxx, four independent functions $$p^{- -}, p^{- +}, p^{+ -}$$p--,p-+,p+-, and $$p^{+ +}$$p++ are used to approximate one-sided derivatives of $$q^-$$q- and $$q^+$$q+. The proposed LDG framework, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a given fully nonlinear problem into a mostly linear system of equations where the given nonlinear differential operator must be replaced by a numerical operator which allows multiple value inputs of the first and second order derivatives $$u_x$$ux and $$u_{xx}$$uxx. An easy to verify set of criteria for constructing "good" numerical operators is also proposed. It consists of consistency and generalized monotonicity. To ensure such a generalized monotonicity property, the crux of the construction is to introduce the numerical moment in the numerical operator, which plays a critical role in the proposed LDG framework. The generalized monotonicity gives the LDG methods the ability to select the viscosity solution among all possible solutions. The proposed framework extends a companion finite difference framework developed by Feng and Lewis (J Comp Appl Math 254:81---98, 2013) and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes. Numerical experiments are also presented to demonstrate the accuracy, efficiency and utility of the proposed LDG methods.