Stabilized nonconforming finite element methods for data assimilation in incompressible flows

We consider a stabilized nonconforming finite element method for data assimilation in incompressible flow subject to the Stokes' equations. The method uses a primal dual structure that allows for the inclusion of nonstandard data. Error estimates are obtained that are optimal compared to the conditional stability of the ill-posed data assimilation problem

Stabilized mixed finite element methods for linear elasticity on simplicial grids in Rn\mathbb{R}^{n}

In this paper, we design two classes of stabilized mixed finite element methods for linear elasticity on simplicial grids. In the first class of elements, we use $\boldsymbol{H}(\mathbf{div}, \Omega; \mathbb{S})$-$P_k$ and $\boldsymbol{L}^2(\Omega; \mathbb{R}^n)$-$P_{k-1}$ to approximate the stress and displacement spaces, respectively, for $1\leq k\leq n$, and employ a stabilization technique in terms of the jump of the discrete displacement over the faces of the triangulation under consideration; in the second class of elements, we use $\boldsymbol{H}_0^1(\Omega; \mathbb{R}^n)$-$P_{k}$ to approximate the displacement space for $1\leq k\leq n$, and adopt the stabilization technique suggested by Brezzi, Fortin, and Marini. We establish the discrete inf-sup conditions, and consequently present the a priori error analysis for them. The main ingredient for the analysis is two special interpolation operators, which can be constructed using a crucial $\boldsymbol{H}(\mathbf{div})$ bubble function space of polynomials on each element. The feature of these methods is the low number of global degrees of freedom in the lowest order case. We present some numerical results to demonstrate the theoretical estimates.Comment: 16 pages, 1 figur

Stable finite element pair for Stokes problem and discrete Stokes complex on quadrilateral grids

In this paper, we first construct a nonconforming finite element pair for the incompressible Stokes problem on quadrilateral grids, and then construct a discrete Stokes complex associated with that finite element pair. The finite element spaces involved consist of piecewise polynomials only, and the divergence-free condition is imposed in a primal formulation. Combined with some existing results, these constructions can be generated onto grids that consist of both triangular and quadrilateral cells

Simulation of rock salt dissolution and its impact on land subsidence

Extensive land subsidence can occur due to subsurface dissolution of evaporites such as halite and gypsum. This paper explores techniques to simulate the salt dissolution forming an intrastratal karst, which is embedded in a sequence of carbonates, marls, anhydrite and gypsum. A numerical model is developed to simulate laminar flow in a subhorizontal void, which corresponds to an opening intrastratal karst. The numerical model is based on the laminar steady-state Stokes flow equation, and the advection dispersion transport equation coupled with the dissolution equation. The flow equation is solved using the nonconforming Crouzeix-Raviart (CR) finite element approximation for the Stokes equation. For the transport equation, a combination between discontinuous Galerkin method and multipoint flux approximation method is proposed. The numerical effect of the dissolution is considered by using a dynamic mesh variation that increases the size of the mesh based on the amount of dissolved salt. The numerical method is applied to a 2D geological cross section representing a Horst and Graben structure in the Tabular Jura of northwestern Switzerland. The model simulates salt dissolution within the geological section and predicts the amount of vertical dissolution as an indicator of potential subsidence that could occur. Simulation results showed that the highest dissolution amount is observed near the normal fault zones, and, therefore, the highest subsidence rates are expected above normal fault zones

The nonconforming virtual element method for eigenvalue problems

We analyse the nonconforming Virtual Element Method (VEM) for the approximation of elliptic eigenvalue problems. The nonconforming VEM allow to treat in the same formulation the two- and three-dimensional case.We present two possible formulations of the discrete problem, derived respectively by the nonstabilized and stabilized approximation of the L^2-inner product, and we study the convergence properties of the corresponding discrete eigenvalue problem. The proposed schemes provide a correct approximation of the spectrum, in particular we prove optimal-order error estimates for the eigenfunctions and the usual double order of convergence of the eigenvalues. Finally we show a large set of numerical tests supporting the theoretical results, including a comparison with the conforming Virtual Element choice