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 paper, an analysis of the superconvergence property of the semidiscrete discontinuous Galerkin (DG) method applied to one-dimensional time-dependent nonlinear scalar conservation laws is carried out. We prove that the error between the DG solution and a particular projection of the exact solution achieves $\left(k + \frac32\right)$-th order superconvergence when upwind fluxes are used. The results hold true for arbitrary nonuniform regular meshes and for piecewise polynomials of degree $k$ ($k \ge 1$), under the condition that $|f'(u)|$ possesses a uniform positive lower bound. Numerical experiments are provided to show that the superconvergence property actually holds true for nonlinear conservation laws with general flux functions, indicating that the restriction on $f(u)$ is artificial.

No abstract uploaded!

No abstract uploaded!

A spectral collocation scheme for the three-dimensional incompressible \(({\varvec{u}},p)\) formulation of the Navier鈥揝tokes equations, in domains \(\varOmega \) with a non-periodic boundary condition, is described. The key feature is the high order approximation, by means of a local Hermite interpolant, of a Neumann boundary condition for use in the numerical solution of the pressure Poisson system. The time updates of the velocity \({\varvec{u}}\) and pressure \(p\) are decoupled as a result of treating the pressure gradient in the momentum equation explicitly in time. The pressure update is computed from a pressure Poisson equation. Extension of the overall methodology to the Boussinesq system is also described. The uncoupling of the pressure and velocity time updates results in a highly efficient scheme that is simple to implement and well suited for simulating moderate to high Reynolds and Rayleigh number flows. Accuracy checks are presented, along with simulations of the lid-driven cavity flow and a differentially heated cavity flow, to demonstrate the scheme produces accurate three-dimensional results at a reasonable computational cost.