Related papers: A primal finite element scheme of the Hodge Laplac…
A nonlinear Helmholtz equation (NLH) with high wave number and Sommerfeld radiation condition is approximated by the perfectly matched layer (PML) technique and then discretized by the linear finite element method (FEM).…
Due to the nonlocal feature of fractional differential operators, the numerical solution to fractional partial differential equations usually requires expensive memory and computation costs. This paper develops a fast scheme for fractional…
In this work, a method for solving the constraints of general relativity is presented, where first all geometrical objects are written in terms of a set of orthonormal triads and a flat Weitzenbock connection, which depends on the triads…
A high-order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the $Q_k$ finite element with mass lumping on rectangular grids, the second-order…
In this paper, we study the mixed element schemes of the Reissner-Mindlin plate model and the Kirchhoff plate model in multiply-connected domains. By a regular decomposition of $H_0({\rm rot},\Omega)$ and a Helmholtz decomposition of its…
By time discretization of a second-order primal-dual dynamical system with damping $\alpha/t$ where an inertial construction in the sense of Nesterov is needed only for the primal variable, we propose a fast primal-dual algorithm for a…
In this paper, we develop two fully nonconforming (both H(grad curl)-nonconforming and H(curl)-nonconforming) finite elements on cubical meshes which can fit into the Stokes complex. The newly proposed elements have 24 and 36 degrees of…
In this paper we consider an initial-boundary value problem with a Caputo time derivative of order $\alpha\in(0,1)$. The solution typically exhibits a weak singularity near the initial time and this causes a reduction in the orders of…
In this paper, we propose an inexact proximal Newton-type method for nonconvex composite problems. We establish the global convergence rate of the order $\mathcal{O}(k^{-1/2})$ in terms of the minimal norm of the KKT residual mapping and…
We analyze a divergence based first order system least squares method applied to a second order elliptic model problem with homogeneous boundary conditions. We prove optimal convergence in the $L^2(\Omega)$ norm for the scalar variable.…
We propose and analyze a general framework for space-time finite element methods that is based on least-squares finite element methods for solving a first-order reformulation of the thick parabolic obstacle problem. Discretizations based on…
We study the most general class of linear boundary-value problems for systems of $r$-th order ordinary differential equations whose solutions range over the complex H\"older space $C^{n+r,\alpha}$, with $0\leq n\in\mathbb{Z}$ and…
We study the variational problem $$\inf \{\lambda_k(\Omega): \Omega\ \textup{open in}\ \R^m,\ |\Omega| < \infty, \ \h(\partial \Omega) \le 1 \},$$ where $\lambda_k(\Omega)$ is the $k$'th eigenvalue of the Dirichlet Laplacian acting in…
The paper considers the initial-boundary value problem for equation $D^\rho_t u(x,t)+ (-\Delta)^\sigma u(x,t)=0$, $\rho\in (0,1)$, $\sigma>0$, in an N-dimensional domain $\Omega$ with a homogeneous Dirichlet condition. The fractional…
We develop a novel unified randomized block-coordinate primal-dual algorithm to solve a class of nonsmooth constrained convex optimization problems, which covers different existing variants and model settings from the literature. We prove…
We hybridize the methods of finite element exterior calculus for the Hodge-Laplace problem on differential $k$-forms in $\mathbb{R}^n$. In the cases $k = 0$ and $k = n$, we recover well-known primal and mixed hybrid methods for the scalar…
We consider time-harmonic linear elasticity equations in domains containing two-dimensional semi-infinite strips. Since for such problems there exist modes with different signs of group and phase velocity, standard perfectly matched layer…
It is shown that the h-adaptive mixed finite element method for the discretization of eigenvalue clusters of the Laplace operator produces optimal convergence rates in terms of nonlinear approximation classes. The results are valid for the…
Due to its highly oscillating solution, the Helmholtz equation is numerically challenging to solve. To obtain a reasonable solution, a mesh size that is much smaller than the reciprocal of the wavenumber is typically required (known as the…
We consider fourth order singularly perturbed eigenvalue problems in one-dimension and the approximation of their solution by the $h$ version of the Finite Element Method (FEM). In particular, we use piecewise Hermite polynomials of degree…