To easily generalize the maximum-principle-satisfying schemes for scalar conservation laws to convection diffusion equations, we propose a non-conventional high order finite volume weighted essentially non-oscillatory (WENO) scheme which can be proven maximum-principle-satisfying. Two-dimensional extensions are straightforward. We also show that the same idea can be used to construct high order schemes preserving the maximum principle for two-dimensional incompressible Navier-Stokes equations in the vorticity stream-function formulation. Numerical tests for the fifth order WENO schemes are reported.

This is the first of a series of papers on the subject of projection methods for viscous incompressible flow calculations. The purpose of these papers is to provide a thorough understanding of the numerical phenomena involved in the projection methods, particularly when boundaries are present, and point to ways of designing more efficient, robust, and accurate numerical methods based on the primitive variable formulation. This paper contains the following topics: 1. convergence and optimal error estimates for both velocity and pressure up to the boundary; 2. explicit characterization of the numerical boundary layers in the pressure approximations and the intermediate velocity fields; 3. the effect of choosing different numerical boundary conditions at the projection step. We will show that a different choice of boundary conditions gives rise to different boundary layer structures. In particular, the straightforward Dirichlet boundary condition for the pressure leads to O(1) numerical boundary layers in the pressure and deteriorates the accuracy in the interior; and 4. postprocessing the numerical solutions to get more accurate approximations for the pressure.

In this paper, we extend the high order finite-difference method with subcell resolution (SR) for two-species stiff one-reaction models to multispecies and multireaction chemical reactive flows, which are significantly more difficult because of the multiple scales generated by different reactions. For reaction problems, when the reaction time scale is very small, the reaction zone scale is also small and the governing equations become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present SR method for reactive Euler system is a fractional step method. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with certain computed flow variables in the shock region modified by the Harten subcell resolution idea. Several numerical examples of multispecies and multireaction reactive flows are performed in both one and two dimensions. Studies demonstrate that the SR method can capture the correct propagation speed of discontinuities in very coarse meshes.

The computation of compressible flows becomes more challenging when the Mach number has different orders of magnitude. When the Mach number is of order one, modern shock capturing methods are able to capture shocks and other complex structures with high numerical resolutions. However, if the Mach number is small, the acoustic waves lead to stiffness in time and exces-sively large numerical viscosity, thus demanding much smaller time step and mesh size than normally needed for incompressible flow simulation. In this paper, we develop an all-speed asymptotic preserving (AP) numerical scheme for the compressible isentropic Euler and Navier-Stokes equations that is uni-formly stable and accurate for all Mach numbers. Our idea is to split the system into two parts: one involves a slow, nonlinear and conservative hyper-bolic system adequate for the use of modern shock capturing methods, and the other a linear hyperbolic system which contains the stiff acoustic dynam-ics, to be solved implicitly. This implicit part is reformulated into a standard pressure Poisson projection system, and thus possesses sufficient structure for efficient fast Fourier transform solution techniques. In the zero Mach number limit, the scheme automatically becomes a projection method-like incompress-ible solver. We present numerical results in one and two dimensions in both compressible and incompressible regimes.

We investigate the angular distribution of Ly$\alpha$ photons transferring in or emergent from an optically thick medium. Since the evolutions of specific intensity $I$ in the frequency space and the angular space are coupled with each other, we first develop the WENO numerical solver in order to find the time-dependent solutions of the integro-differential equation of $I$ in frequency and angular space simultaneously. We first show that the solutions with the Eddington approximation, which assume $I$ to be linearly dependent on the angular variable $\mu$, yield similar frequency profiles of the photon flux as that without the Eddington approximation. However, the solutions of the $\mu$ distribution evolution are significantly different from that given by Eddington approximation. First, the angular distribution of $I$ are found to be substantially dependent on the frequency of photons. For photons with the resonant frequency $\nu_0$, $I$ contains only a linear term of $\mu$. For photons with frequency at the double peaks of the flux, the $\mu$-distribution is highly anisotropic, in which most photons are in the direction of radial forward. Moreover, either at $\nu_0$ or at the double peaks, the $\mu$ distributions actually are independent of the initial $\mu$ distribution of photons of the source. This is because the photons with frequency either of $\nu_0$ or of the double peaks undergo the process of forgetting their initial conditions due to the resonant scattering. We also show that the optically thick medium is a collimator of photons at the double peaks. Photons from the double peaks form a forward beam with very small spread angle.