This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high

Local discontinuous Galerkin (LDG) methods are popular for convection-diffusion equations. In LDG methods, we introduce an auxiliary variable $p$ to represent the derivative of the primary variable $u$, and solve them on the same mesh. In this paper, we will introduce a new LDG method, and solve $u$ and $p$ on different meshes. The stability and error estimates will be investigated. The new algorithm is more flexible and flux-free for pure diffusion equations without introducing additional computational cost compared with the original LDG methods, since it is not necessary to solve each equation twice. Moreover, it is possible to construct third-order maximum-principle-preserving schemes based on the new algorithm. However, one cannot anticipate optimal accuracy in some special cases. In this paper, we will find out the reason for accuracy degeneration which further leads to several alternatives to obtain optimal convergence rates. Finally, several numerical experiments will be given to demonstrate the stability and optimal accuracy of the new algorithm.

A new residual-based stabilized finite element method is analyzed for solving the Stokes equations in terms of velocity and pressure, where the $H^{−1}$ norm is introduced in the measurement of the residuals to obtain a symmetric positive definite (SPD) method. The $H^{−1}$ norm is computable and can be always easily realized offline by the continuous linear finite element solution or the preconditioned counterpart of the Poisson Dirichlet problem. Although the $H^{−1}$ norm is computed in the linear element space, no matter what the finite element spaces for the velocity and the pressure are, optimal error bounds can be established when using continuous finite element pairs $R_l$−$R_m$ for velocity and pressure for any $l, m\ge 1$. Numerical experiments are performed to confirm the theoretical results obtained.

This paper deals with the numerical resolution of kinetic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. We focus on the kinetic model considered in \cite{DM, DLMP} where alignment is taken into account in addition of an attraction-repulsion interaction potential. We apply a discontinuous Galerkin method for the free transport and non-local drift velocity together with a spectral method for the velocity variable. Then, we analyse consistency and stability of the semi-discrete scheme. We propose several numerical experiments which provide a solid validation of the method and its underlying concepts.

We present two uniformly accurate numerical methods for discretizing the Zakharovsystem (ZS) with a dimensionless parameter 0< ε ≤ 1, which is inversely proportional to theacoustic speed. In the subsonic limit regime, i.e., 0< ε << 1, the solution of ZS propagates waves with O(ε)- andO(1)-wavelengths in time and space, respectively, and/or rapid outgoing initial layerswith speed O(1/ε) in space due to the singular perturbation of the wave operator in ZS and/or theincompatibility of the initial data. By adopting an asymptotic consistent formulation of ZS, wepresent a time-splitting exponential wave integrator (TS-EWI) method via applying a time-splittingtechnique and an exponential wave integrator for temporal derivatives in the nonlinear Schr ̈odingerequation and wave-type equation, respectively. By introducing a multiscale decomposition of ZS, wepropose a time-splitting multiscale time integrator (TS-MTI) method. Both methods are explicitand convergent exponentially in space for all kinds of initial data, which is uniformly for ε ∈ (0,1].The TS-EWI method is simpler to be implemented and it is only uniformly and optimally accuratein time for well-prepared initial data, while the TS-MTI method is uniformly and optimally accuratein time for any kind of initial data. Extensive numerical results are reported to show their efficiencyand accuracy, especially in the subsonic limit regime. Finally, the TS-MTI method is applied tostudy numerically convergence rates of ZS to its limiting models when ε→0+.