No abstract

In this work, we investigate the numerical approximation of the compressible Navier-Stokes equations under the framework of discontinuous Galerkin methods. For discretization of the viscous and heat fluxes, we extend and apply the symmetric direct discontinuous Galerkin (SDDG) method which is originally introduced for scalar diffusion problems. The original compressible Navier-Stokes equations are rewritten into an equivalent form via homogeneity tensors. Then, the numerical diffusive fluxes are constructed from the weak formulation of primal equations directly without converting the second-order equations to a first-order system. Additional numerical flux functions involving the jump of second order derivative of test functions are added to the original direct discontinuous Galerkin (DDG) discretization. A number of numerical tests are carried out to assess the practical performance of the SDDG method for the two dimensional compressible Navier-Stokes equations. These numerical results obtained demonstrate that the SDDG method can achieve the optimal order of accuracy. Especially, compared with the well-established symmetric interior penalty (SIP) method [18], the SDDG method can maintain the expected optimal order of convergence with a smaller penalty coefficient.

An adaptive mesh refinement algorithm based on a continuous adjoint approach is developed. Both the primal equation and the adjoint equation are approximated with the discontinuous Galerkin (DG) method. The proposed adaptive algorithm is used in compressible Euler equations. Numerical tests are made to show the superiority of the proposed adaptive algorithm.

No abstract

The aim of this work is to develop a hybridizable discontinuous Galerkin method for elliptic problems. In the proposed method, the numerical flux functions are constructed from the weak formulation of primal equation directly without converting the second-order equation to a first-order system. In order to guarantee the stability and convergence of the method, we derive a computable lower bound for the constant in numerical flux functions. We also establish a prior error estimation and give some theoretical analysis for the proposed method. Finally, a numerical experiment is presented to verify the theoretical results.