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.

Using the convex integration technique for the three-dimensional Navier-Stokes equations introduced by T. Buckmaster and V. Vicol, it is shown the existence of non-unique weak solutions for the 3D Navier-Stokes equations with fractional hyperviscosity $(-\Delta)^{\theta}$, whenever the exponent $\theta$ is less than J.-L. Lions' exponent $5/4$, i.e., when $\theta < 5/4$.

Adapting the convex integration technique introduced in \cite{dLSz5} and subsequently developed in \cite{IsettVicol,IsettOh16,Buckmaster2013transporting}, we construct H\"{o}lder continuous weak solutions to the three dimensional Prandtl system and some other models with vertical viscosity.

We consider admissible weak solutions to the compressible Euler system with source terms, which include rotating shallow water system and the Euler system with damping as special examples. In the case of anti-symmetric sources such as rotations, for general piecewise Lipschitz initial densities and some suitably constructed initial momentum, we obtain infinitely many global admissible weak solutions. Furthermore, we construct a class of finite-states admissible weak solutions to the Euler system with anti-symmetric sources. Under the additional smallness assumption on the initial densities, we also obtain multiple global-in-time admissible weak solutions for more general sources including damping. The basic framework are based on the convex integration method developed by De~Lellis and Sz\'{e}kelyhidi \cite{dLSz1,dLSz2} for the Euler system. One of the main ingredients of this paper is the construction of specified localized plane wave perturbations which are compatible with a given source term.

No abstract uploaded!