In this paper, we propose a multi-scale discontinuous Galerkin (DG) method for second-order elliptic problems with curvilinear unidirectional rough coefficients by choosing a special non-polynomial approximation space. The key ingredient of the method lies in the incorporation of the local oscillatory features of the differential operators into the approximation space so as to capture the multi-scale solutions without having to resolve the finest scales. The unidirectional feature of the rough coefficients allows us to construct the basis functions of the DG non-polynomial approximation space explicitly, thereby greatly increasing the algorithm efficiency. Detailed error estimates for two-dimensional second-order DG methods are derived, and a general guidance on how to construct such non-polynomial basis is discussed. Numerical examples are also presented to validate and demonstrate the effectiveness of the algorithm.

In this work, we study the unsteady axisymmetric mixed convection boundary layer flow and heat transfer of non-Newtonian power law fluid over a cylinder. Different from most classical works, the temperature dependent variable fluid viscosity and thermal conductivity are taken into account in highly coupled velocity and temperature fields. The motion of the fluid can be modeled by a time-dependent nonlinear parabolic system in cylindrical coordinates, which is solved numerically by using Chebyshev spectral method along with the strong stability preserving (SSP) third order Runge-Kutta time discretization. We apply the numerical solver to problems with different power law indices, viscosity parameter, thermal conductivity parameter and Richardson numbers, and compute up to the steady state. The numerical solver is checked by testing the spectral convergence of the numerical approximation to a smooth exact solution of the PDEs with source terms. Moreover, the combined effects of pertinent physical parameters on the flow and heat transfer characteristics are analyzed in detail.

Solving wave propagation problems within heterogeneous media has been of great interest and has a wide range of applications in physics and engineering. The design of numerical methods for such general wave propagation problems is challenging because the energy conserving property has to be incorporated in the numerical algorithms in order to minimize the phase or shape errors after long time integration. In this paper, we focus on multi-dimensional wave problems and consider linear second-order wave equation in heterogeneous media. We develop and analyze an LDG method, in which numerical fluxes are carefully designed to maintain the energy preserving property and accuracy. Compatible high order energy conserving time integrators are also proposed. The optimal error estimates and the energy conserving property are proved for the semi-discrete methods. Our numerical experiments demonstrate optimal rates of convergence, and show that the errors of the numerical solutions do not grow significantly in time due to the energy conserving property.