Analysis of a Darcy-Cahn-Hilliard Diffuse Interface Model for the Hele-Shaw Flow and its Fully Discrete Finite Element Approximation

In this paper we present PDE and finite element analyses for a system of partial differential equations (PDEs) consisting of the Darcy equation and the Cahn-Hilliard equation, which arises as a diffuse interface model for the two phase Hele-Shaw flow. We propose a fully discrete implicit finite element method for approximating the PDE system, which consists of the implicit Euler method combined with a convex splitting energy strategy for the temporal discretization, the standard finite element discretization for the pressure and a split (or mixed) finite element discretization for the fourth order Cahn-Hilliard equation. It is shown that the proposed numerical method satisfies a mass conservation law in addition to a discrete energy law that mimics the basic energy law for the Darcy-Cahn-Hilliard phase field model and holds uniformly in the phase field parameter $\epsilon$. With help of the discrete energy law, we first prove that the fully discrete finite method is unconditionally energy stable and uniquely solvable at each time step. We then show that, using the compactness method, the finite element solution has an accumulation point that is a weak solution of the PDE system. As a result, the convergence result also provides a constructive proof of the existence of global-in-time weak solutions to the Darcy-Cahn-Hilliard phase field model in both two and three dimensions. Finally, we propose a nonlinear multigrid iterative algorithm to solve the finite element equations at each time step. Numerical experiments based on the overall solution method of combining the proposed finite element discretization and the nonlinear multigrid solver are presented to validate the theoretical results and to show the effectiveness of the proposed fully discrete finite element method for approximating the Darcy-Cahn-Hilliard phase field model.Comment: 30 pages, 4 tables, 2 figure

Dual virtual element method for discrete fractures networks

Discrete fracture networks is a key ingredient in the simulation of physical processes which involve fluid flow in the underground, when the surrounding rock matrix is considered impervious. In this paper we present two different models to compute the pressure field and Darcy velocity in the system. The first allows a normal flow out of a fracture at the intersections, while the second grants also a tangential flow along the intersections. For the numerical discretization, we use the mixed virtual finite element method as it is known to handle grid elements of, almost, any arbitrary shape. The flexibility of the discretization allows us to loosen the requirements on grid construction, and thus significantly simplify the flow discretization compared to traditional discrete fracture network models. A coarsening algorithm, from the algebraic multigrid literature, is also considered to further speed up the computation. The performance of the method is validated by numerical experiments

Virtual Element Methods for hyperbolic problems on polygonal meshes

In the present paper we develop the Virtual Element Method for hyperbolic problems on polygonal meshes, considering the linear wave equations as our model problem. After presenting the semi-discrete scheme, we derive the convergence estimates in H^1 semi-norm and L^2 norm. Moreover we develop a theoretical analysis on the stability for the fully discrete problem by comparing the Newmark method and the Bathe method. Finally we show the practical behaviour of the proposed method through a large array of numerical tests