We construct a local Lax-Friedrichs type positivity-preserving flux for compressible Navier-Stokes equations, which can be easily extended to high dimensions for generic forms of equations of state, shear stress tensor and heat flux. With this positivity-preserving flux, any finite volume type schemes including discontinuous Galerkin (DG) schemes with strong stability preserving Runge-Kutta time discretizations satisfy a weak positivity property. With a simple and effi- cient positivity-preserving limiter, high order explicit Runge-Kutta DG schemes are rendered preserving the positivity of density and internal energy without los- ing local conservation or high order accuracy. Numerical tests suggest that the positivity-preserving flux and the positivity-preserving limiter do not induce excessive artificial viscosity, and the high order positivity-preserving DG schemes without other limiters can produce satisfying non-oscillatory solutions when the nonlinear diffusion in compressible Navier-Stokes equations is accurately resolved.

This is part II of our paper in which we propose and develop a dynamically bi-orthogonal method (DyBO) to study a class of time-dependent stochastic partial differential equations (SPDEs) whose solutions enjoy a low-dimensional structure. In part I of our paper [9], we derived the DyBO formulation and proposed numerical algorithms based on this formulation. Some important theoretical results regarding consistency and bi-orthogonality preservation were also established in the first part along with a range of numerical examples to illustrate the effectiveness of the DyBO method. In this paper, we focus on the computational complexity analysis and develop an effective adaptivity strategy to add or remove modes dynamically. Our complexity analysis shows that the ratio of computational complexities between the DyBO method and a generalized polynomial chaos method (gPC) is roughly of order O ((m/N p) 3) for a

Poisson equation on point cloud with Dirichlet boundary condition plays important role in many problems. In this paper, we use the volume constraint proposed by Du et.al to handle the Dirichlet boundary condition in the point integral method for Poisson equation on point cloud. We prove that the solution given by volume constraint converges to the true solution as the point cloud converges to the underlying smooth manifold.

We study semi-Lagrangian discontinuous Galerkin (SLDG) and Runge-Kutta discontinuous Galerkin (RKDG) schemes for some front propagation problems in the presence of an obstacle term, modeled by a nonlinear Hamilton-Jacobi equation of the form $\min(u_t + c u_x, u - g(x))=0$, in one space dimension. New convergence results and error bounds are obtained for Lipschitz regular data. These ``low regularity" assumptions are the natural ones for the solutions of the studied equations. Numerical tests are given to illustrate the behavior of our schemes.

For a class of convection-diffusion equations with variable diffusivity, we construct third order accurate discontinuous Galerkin (DG) schemes on both one and two dimensional rectangular meshes. The DG method with an explicit time stepping can well be applied to nonlinear convection–diffusion equations. It is shown that under suitable time step restrictions, the scaling limiter proposed in Liu and Yu (2014) [23]when coupled with the present DG schemes preserves the solution bounds indicated by the initial data, i.e., the maximum principle, while maintaining uniform third order accuracy. These schemes can be extended to rectangular meshes in three dimension. The crucial for all model scenarios is that an effective test set can be identified to verify the desired bounds of numerical solutions. This is achieved mainly by taking advantage of the flexible form of the diffusive flux and the adaptable decomposition of weighted cell averages. Numerical results are presented to validate the numerical methods and demonstrate their effectiveness.