Nonconforming Virtual Element Method for 2m2m-th Order Partial Differential Equations in Rn\mathbb R^n

A unified construction of the $H^m$-nonconforming virtual elements of any order $k$ is developed on any shape of polytope in $\mathbb R^n$ with constraints $m\leq n$ and $k\geq m$. As a vital tool in the construction, a generalized Green's identity for $H^m$ inner product is derived. The $H^m$-nonconforming virtual element methods are then used to approximate solutions of the $m$-harmonic equation. After establishing a bound on the jump related to the weak continuity, the optimal error estimate of the canonical interpolation, and the norm equivalence of the stabilization term, the optimal error estimates are derived for the $H^m$-nonconforming virtual element methods.Comment: 33page

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

Multigrid Methods for Hellan-Herrmann-Johnson Mixed Method of Kirchhoff Plate Bending Problems

A V-cycle multigrid method for the Hellan-Herrmann-Johnson (HHJ) discretization of the Kirchhoff plate bending problems is developed in this paper. It is shown that the contraction number of the V-cycle multigrid HHJ mixed method is bounded away from one uniformly with respect to the mesh size. The uniform convergence is achieved for the V-cycle multigrid method with only one smoothing step and without full elliptic regularity. The key is a stable decomposition of the kernel space which is derived from an exact sequence of the HHJ mixed method, and the strengthened Cauchy Schwarz inequality. Some numerical experiments are provided to confirm the proposed V-cycle multigrid method. The exact sequences of the HHJ mixed method and the corresponding commutative diagram is of some interest independent of the current context.Comment: 23 page

Nonconforming finite element Stokes complexes in three dimensions

Two nonconforming finite element Stokes complexes ended with the nonconforming $P_1$-$P_0$ element for the Stokes equation in three dimensions are constructed. And commutative diagrams are also shown by combining nonconforming finite element Stokes complexes and interpolation operators. The lower order $\boldsymbol H(\textrm{grad}~\textrm{curl})$-nonconforming finite element only has $14$ degrees of freedom, whose basis functions are explicitly given in terms of the barycentric coordinates. The $\boldsymbol H(\textrm{grad}~\textrm{curl})$-nonconforming elements are applied to solve the quad-curl problem, and optimal convergence is derived. By the nonconforming finite element Stokes complexes, the mixed finite element methods of the quad-curl problem is decoupled into two mixed methods of the Maxwell equation and the nonconforming $P_1$-$P_0$ element method for the Stokes equation, based on which a fast solver is developed.Comment: 20 page

OPR

The ability to reproduce a parallel execution is desirable for debugging and program reliability purposes. In debugging (13), the programmer needs to manually step back in time, while for resilience (6) this is automatically performed by the the application upon failure. To be useful, replay has to faithfully reproduce the original execution. For parallel programs the main challenge is inferring and maintaining the order of conflicting operations (data races). Deterministic record and replay (R&R) techniques have been developed for multithreaded shared memory programs (5), as well as distributed memory programs (14). Our main interest is techniques for large scale scientific (3; 4) programming models

Residual-Based A Posteriori Error Estimates for Symmetric Conforming Mixed Finite Elements for Linear Elasticity Problems

A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems of Dirichlet and mixed boundary conditions are proposed. Stability and efficiency of the estimators are proved. Finally, we provide numerical examples to verify the theoretical results