No abstract

An efficient nonlinear multigrid method for a mixed finite element method of the Darcy-Forchheimer model is constructed in this paper. A Peaceman-Rachford type iteration is used as a smoother to decouple the nonlinearity from the divergence constraint. The nonlinear equation can be solved element-wise with a closed formulae. The linear saddle point system for the constraint is reduced into a symmetric positive definite system of Poisson type. Furthermore an empirical choice of the parameter used in the splitting is proposed and the resulting multigrid method is robust to the so-called Forchheimer number which controls the strength of the nonlinearity. By comparing the number of iterations and CPU time of different solvers in several numerical experiments, our multigrid method is shown to convergent with a rate independent of the mesh size and the Forchheimer number and with a nearly linear computational cost.

A finite difference scheme for the one-dimensional space fractional diffusion equation is presented and analysed. The scheme is constructed by modifying the shifted Grünwald approximation to the spatial fractional derivative and using an asymmetric discretisation technique. By calculating the unknowns in differential nodal point sequences at the odd and even time levels, the discrete solution of the scheme can be obtained explicitly. We prove that the scheme is uniformly stable. The error between the discrete solution and the analytical solution in the discrete l2 norm is optimal in some cases. Numerical results for several examples are consistent with the theoretical analysis.

In this paper, the solution of the Darcy-Forchheimer model in high contrast heterogeneous media is studied. This problem is solved by a mixed finite element method (MFEM) on a fine grid (the reference solution), where the pressure is approximated by piecewise constant elements; meanwhile, the velocity is discretized by the lowest order Raviart-Thomas elements. The solution on a coarse grid is performed by using the mixed generalized multiscale finite element method (mixed GMsFEM). The nonlinear equation can be solved by the well known Picard iteration. Several numerical experiments are presented in a two-dimensional heterogeneous domain to show the good applicability of the proposed multiscale method.

No abstract