Source author record

Michael Neilan

Michael Neilan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis math.AP

Catalog footprint

What is connected

12works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

A Note on the Shape Regularity of Worsey-Farin Splits

We prove three-dimensional Worsey-Farin refinements inherit their parent triangulations' shape-regularity.

preprint2022arXiv

A CutFEM divergence-free discretization for the Stokes problem

We construct and analyze a CutFEM discretization for the Stokes problem based on the Scott-Vogelius pair. The discrete piecewise polynomial spaces are defined on macro-element triangulations which are not fitted to the smooth physical domain. Boundary conditions are imposed via penalization through the help of a Nitsche-type discretization, whereas stability with respect to small and anisotropic cuts of the bulk elements is ensured by adding local ghost penalty stabilization terms. We show stability of the scheme as well as a divergence--free property of the discrete velocity outside an $O(h)$ neighborhood of the boundary. To mitigate the error caused by the violation of the divergence-free condition, we introduce local grad-div stabilization. The error analysis shows that the grad-div parameter can scale like $O(h^{-1})$, allowing a rather heavy penalty for the violation of mass conservation, while still ensuring optimal order error estimates.

preprint2022arXiv

Convergence of Lagrange Finite Element Methods for Maxwell Eigenvalue Problem in 3D

We prove convergence of the Maxwell eigenvalue problem using quadratic or higher Lagrange finite elements on Worsey-Farin splits in three dimensions. To do this, we construct two Fortin-like operators to prove uniform convergence of the corresponding source problem. We present numerical experiments to illustrate the theoretical results.

preprint2021arXiv

Convergence of Lagrange finite elements for the Maxwell Eigenvalue Problem in 2D

We consider finite element approximations of the Maxwell eigenvalue problem in two dimensions. We prove, in certain settings, convergence of the discrete eigenvalues using Lagrange finite elements. In particular, we prove convergence in three scenarios: piecewise linear elements on Powell--Sabin triangulations, piecewise quadratic elements on Clough--Tocher triangulations, and piecewise quartics (and higher) elements on general shape-regular triangulations. We provide numerical experiments that support the theoretical results. The computations also show that, on general triangulations, the eigenvalue approximations are very sensitive to nearly singular vertices, i.e., vertices that fall on exactly two "almost" straight lines.

preprint2021arXiv

Low-order divergence-free approximations for the Stokes problem on Worsey-Farin and Powell-Sabin splits

We derive low-order, inf-sup stable and divergence-free finite element approximations for the Stokes problem using Worsey-Farin splits in three dimensions and Powell-Sabin splits in two dimensions. The velocity space simply consists of continuous, piecewise linear polynomials, where as the pressure space is a subspace of piecewise constants with weak continuity properties at singular edges (3D) and singular vertices (2D). We discuss implementation aspects that arise when coding the pressure space, and in particular, show that the pressure constraints can be enforced at an algebraic level.

preprint2020arXiv

Divergence--free Scott--Vogelius elements on curved domains

We construct and analyze an isoparametric finite element pair for the Stokes problem in two dimensions. The pair is defined by mapping the Scott-Vogelius finite element space via a Piola transform. The velocity space has the same degrees of freedom as the quadratic Lagrange finite element space, and therefore the proposed spaces reduce to the Scott-Vogelius pair in the interior of the domain. We prove that the resulting method converges with optimal order, is divergence--free, and is pressure robust. Numerical examples are provided which support the theoretical results.

preprint2016arXiv

Interior Penalty Discontinuous Galerkin Methods for Second Order Linear Non-Divergence Form Elliptic PDEs

This paper develops interior penalty discontinuous Galerkin (IP-DG) methods to approximate $W^{2,p}$ strong solutions of second order linear elliptic partial differential equations (PDEs) in non-divergence form with continuous coefficients. The proposed IP-DG methods are closely related to the IP-DG methods for advection-diffusion equations, and they are easy to implement on existing standard IP-DG software platforms. It is proved that the proposed IP-DG methods have unique solutions and converge with optimal rate to the $W^{2,p}$ strong solution in a discrete $W^{2,p}$-norm. The crux of the analysis is to establish a DG discrete counterpart of the Calderon-Zygmund estimate and to adapt a freezing coefficient technique used for the PDE analysis at the discrete level. As a byproduct of our analysis, we also establish broken $W^{1,p}$-norm error estimates for IP-DG approximations of constant coefficient elliptic PDEs. Numerical experiments are provided to gauge the performance of the proposed IP-DG methods and to validate the theoretical convergence results.

preprint2016arXiv

Numerical analysis of strongly nonlinear PDEs

We review the construction and analysis of numerical methods for strongly nonlinear PDEs, with an emphasis on convex and nonconvex fully nonlinear equations and the convergence to viscosity solutions. We begin by describing a fundamental result in this area which states that stable, consistent, and monotone schemes converge as the discretization parameter tends to zero. We review methodologies to construct finite difference, finite element, and semi-Lagrangian schemes that satisfy these criteria, and, in addition, discuss some rather novel tools that have paved the way to derive rates of convergence within this framework.

preprint2015arXiv

$C^0$ discontinuous Galerkin finite element methods for second order linear elliptic partial differential equations in non-divergence form

This paper is concerned with finite element approximations of $W^{2,p}$ strong solutions of second-order linear elliptic partial differential equations (PDEs) in non-divergence form with continuous coefficients. A nonstandard (primal) finite element method, which uses finite-dimensional subspaces consisting globally continuous piecewise polynomial functions, is proposed and analyzed. The main novelty of the finite element method is to introduce an interior penalty term, which penalizes the jump of the flux across the interior element edges/faces, to augment a nonsymmetric piecewise defined and PDE-induced bilinear form. Existence, uniqueness and error estimate in a discrete $W^{2,p}$ energy norm are proved for the proposed finite element method. This is achieved by establishing a discrete Calderon-Zygmund-type estimate and mimicking strong solution PDE techniques at the discrete level. Numerical experiments are provided to test the performance of proposed finite element method and to validate the convergence theory.

preprint2013arXiv

Discontinuous Galerkin finite element differential calculus and applications to numerical solutions of linear and nonlinear partial differential equations

This paper develops a discontinuous Galerkin (DG) finite element differential calculus theory for approximating weak derivatives of Sobolev functions and piecewise Sobolev functions. By introducing numerical one-sided derivatives as building blocks, various first and second order numericaloperators such as the gradient, divergence, Hessian, and Laplacian operator are defined, and their corresponding calculus rules are established. Among the calculus rules are product and chain rules, integration by parts formulas and the divergence theorem. Approximation properties and the relationship between the proposed DG finite element numerical derivatives and some well-known finite difference numerical derivative formulas on Cartesian grids are also established. Efficient implementation of the DG finite element numerical differential operators is also proposed. Besides independent interest in numerical differentiation, the primary motivation and goal of developing the DG finite element differential calculus is to solve partial differential equations. It is shown that several existing finite element, finite difference and DG methods can be rewritten compactly using the proposed DG finite element differential calculus framework. Moreover, new DG methods for linear and nonlinear PDEs are also obtained from the framework.

preprint2012arXiv

Convergence of a fourth order singular perturbation of the $n$-dimensional radially symmetric Monge-Ampère equation

This paper concerns with the convergence analysis of a fourth order singular perturbation of the Dirichlet Monge-Ampère problem in the $n$-dimensional radial symmetric case. A detailed study of the fourth order problem is presented. In particular, various {\em a priori} estimates with explicit dependence on the perturbation parameter $\vepsi$ are derived, and a crucial convexity property is also proved for the solution of the fourth order problem. Using these estimates and the convexity property, we prove that the solution of the perturbed problem converges uniformly and compactly to the unique convex viscosity solution of the Dirichlet Monge-Ampère problem. Rates of convergence in the $H^k$-norm for $k=0,1,2$ are established, and illustrating numerical experiment results are also presented in the paper.

preprint2011arXiv

The Vanishing Moment Method for Fully Nonlinear Second Order Partial Differential Equations: Formulation, Theory, and Numerical Analysis

The vanishing moment method was introduced by the authors in [37] as a reliable methodology for computing viscosity solutions of fully nonlinear second order partial differential equations (PDEs), in particular, using Galerkin-type numerical methods such as finite element methods, spectral methods, and discontinuous Galerkin methods, a task which has not been practicable in the past. The crux of the vanishing moment method is the simple idea of approximating a fully nonlinear second order PDE by a family (parametrized by a small parameter $\vepsi$) of quasilinear higher order (in particular, fourth order) PDEs. The primary objectives of this book are to present a detailed convergent analysis for the method in the radial symmetric case and to carry out a comprehensive finite element numerical analysis for the vanishing moment equations (i.e., the regularized fourth order PDEs). Abstract methodological and convergence analysis frameworks of conforming finite element methods and mixed finite element methods are first developed for fully nonlinear second order PDEs in general settings. The abstract frameworks are then applied to three prototypical nonlinear equations, namely, the Monge-Ampère equation, the equation of prescribed Gauss curvature, and the infinity-Laplacian equation. Numerical experiments are also presented for each problem to validate the theoretical error estimate results and to gauge the efficiency of the proposed numerical methods and the vanishing moment methodology.