No abstract uploaded!

In this paper, we develop high-order bound-preserving (BP) discontinuous Galerkin (DG) methods for the coupled system of compressible miscible displacements on triangular meshes. We consider the problem with multi-component fluid mixture and the (volumetric) concentration of the $j$th component, $c_j$, should be between 0 and 1. There are three main difficulties. Firstly, $c_j$ does not satisfy a maximum-principle. Therefore, the numerical techniques introduced in (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 3091-3120) cannot be applied directly. The main idea is to apply the positivity-preserving techniques to all $c_j's$ and enforce $\sum_jc_j=1$ simultaneously to obtain physically relevant approximations. By doing so, we have to treat the time derivative of the pressure $dp/dt$ as a source in the concentration equation and choose suitable fluxes in the pressure and concentration equations. Secondly, it is not easy to construct first-order numerical fluxes for interior penalty DG methods on triangular meshes. One of the key points in the high-order BP technique applied in this paper is the combination of high-order and lower-order numerical fluxes. We will construct second-order BP schemes and use the second-order numerical fluxes as the lower-order one. Finally, the classical slope limiter cannot be applied to $c_j$. To construct the BP technique, we will not approximate $c_j$ directly. Therefore, a new limiter will be introduced. Numerical experiments will be given to demonstrate the high-order accuracy and good performance of the numerical technique.

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.

In this paper, a new type of finite difference Hermite weighted essentially non-oscillatory (HWENO) schemes are constructed for solving Hamilton-Jacobi (HJ) equations. Point values of both the solution and its first derivatives are used in the HWENO reconstruction and evolved via time advancing. While the evolution of the solution is still through the classical numerical fluxes to ensure convergence to weak solutions, the evolution of the first derivatives of the solution is through a simple dimension-by-dimension non-conservative procedure to gain efficiency. The main advantages of this new scheme include its compactness in the spatial field and its simplicity in the reconstructions. Extensive numerical experiments in one and two dimensional cases are performed to verify the accuracy, high resolution and efficiency of this new scheme.

We show that the gradient mapping of the squared norm of FischerBurmeister function is globally Lipschitz continuous and semismooth, which provides a theoretical basis for solving nonlinear second-order cone complementarity problems via the conjugate gradient method and the semismooth Newtons method.