Investigation and analysis of the numerical approach to solve the multi-term time-fractional advection-diffusion model

In this paper, a methodical approach is presented to approximate the multi-term fractional advection-diffusion model (MT-FAD). The Lagrange squared interpolation is used to discretize temporal fractional derivatives, and Legendre polynomials are shifted as an operator to discretize the spatial fractional derivatives. The advantage of these numerical techniques lies in the orthogonality of Legendre polynomials and its matrix operations. A quadratic implicit design as well as its stability and convergence analysis are evaluated. It should be noted that the theoretical proof obtained from this study represents the first results for these numerical schemes. Finally, we provide three numerical examples to verify the validity of the proposed methods and demonstrate their accuracy and effectiveness in comparison with previous studies shown in [W. P. Bu, X. T. Liu, Y. F. Tang, J. Y. Yang, Finite element multigrid method for multi-term time fractional advection diffusion equations, Int. J. Model. Simul. Sci. Comput., 6 (2015), 1540001]

A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure

In this paper we formulate and test numerically a fully-coupled discontinuous Galerkin (DG) method for incompressible two-phase flow with discontinuous capillary pressure. The spatial discretization uses the symmetric interior penalty DG formulation with weighted averages and is based on a wetting-phase potential / capillary potential formulation of the two-phase flow system. After discretizing in time with diagonally implicit Runge-Kutta schemes the resulting systems of nonlinear algebraic equations are solved with Newton's method and the arising systems of linear equations are solved efficiently and in parallel with an algebraic multigrid method. The new scheme is investigated for various test problems from the literature and is also compared to a cell-centered finite volume scheme in terms of accuracy and time to solution. We find that the method is accurate, robust and efficient. In particular no post-processing of the DG velocity field is necessary in contrast to results reported by several authors for decoupled schemes. Moreover, the solver scales well in parallel and three-dimensional problems with up to nearly 100 million degrees of freedom per time step have been computed on 1000 processors