Note on a weakly over-penalised symmetric interior penalty method on anisotropic meshes for the Poisson equation, Ver. 1

Hiroki Ishizaka

2024-03-17

Abstract:The purpose is to make an easy-to-understand note of "Special Topics in Finite Element Methods." There might be typos and mistakes. Therefore, I do not take any responsibility for unauthorised use.

Numerical Analysis

What problem does this paper attempt to address?

Hybrid weakly over-penalised symmetric interior penalty method on anisotropic meshes

Hiroki Ishizaka

2024-06-12

Abstract:In this study, we investigate a hybrid-type anisotropic weakly over-penalised symmetric interior penalty method for the Poisson equation on convex domains. Compared with the well-known hybrid discontinuous Galerkin methods, our approach is simple and easy to implement. Our primary contributions are the proposal of a new scheme and the demonstration of a proof for the consistency term, which allows us to estimate the anisotropic consistency error. The key idea of the proof is to apply the relation between the Raviart--Thomas finite element space and a discontinuous space. In numerical experiments, we compare the calculation results for standard and anisotropic mesh partitions.

Numerical Analysis
Anisotropic weakly over-penalised symmetric interior penalty method for the Stokes equation

Hiroki Ishizaka

2024-04-05

Abstract:In this study, we investigate an anisotropic weakly over-penalised symmetric interior penalty method for the Stokes equation {on convex domains}. Our approach is a simple discontinuous Galerkin method similar to the Crouzeix--Raviart finite element method. As our primary contribution, we show a new proof for the consistency term, which allows us to obtain an estimate of the anisotropic consistency error. The key idea of the proof is to apply the relation between the Raviart--Thomas finite element space and a discontinuous space. While inf-sup stable schemes of the discontinuous Galerkin method on shape-regular mesh partitions have been widely discussed, our results show that the Stokes element satisfies the inf-sup condition on anisotropic meshes. Furthermore, we provide an error estimate in an energy norm on anisotropic meshes. In numerical experiments, we compare calculation results for standard and anisotropic mesh partitions.

Numerical Analysis
An Approach to Numerically Solving the Poisson Equation

Zhichen Feng,Zheng-Mao Sheng

DOI: https://doi.org/10.1088/0031-8949/90/6/065603

2015-01-01

Physica Scripta

Abstract:We introduce an approach for numerically solving the Poisson equation by using a physical model, which is a way to solve a partial differential equation without the finite difference method. This method is especially useful for obtaining the solutions in very many free-charge neutral systems with open boundary conditions. It can be used for arbitrary geometry and mesh style and is more efficient comparing with the widely-used iterative algorithm with multigrid methods. It is especially suitable for parallel computing. This method can also be applied to numerically solving other partial differential equations whose Green functions exist in analytic expression.
Crouzeix–Raviart and Raviart–Thomas finite-element error analysis on anisotropic meshes violating the maximum-angle condition

Hiroki Ishizaka,Kenta Kobayashi,Takuya Tsuchiya

DOI: https://doi.org/10.1007/s13160-020-00455-7

2021-02-09

Japan Journal of Industrial and Applied Mathematics

Abstract:We investigate the piecewise linear nonconforming Crouzeix–Raviart and the lowest order Raviart–Thomas finite-element methods for the Poisson problem on three-dimensional anisotropic meshes. We first give error estimates of the Crouzeix–Raviart and the Raviart–Thomas finite-element approximate problems. We next present the equivalence between the Raviart–Thomas finite-element method and the enriched Crouzeix–Raviart finite-element method. We emphasize that we do not impose either shape-regular or maximum-angle condition during mesh partitioning. Numerical results confirm the results that we obtained.

mathematics, applied
Anisotropic Error Analysis of Weak Galerkin finite element method for Singularly Perturbed Biharmonic Problems

Aayushman Raina,Srinivasan Natesan,Şuayip Toprakseven

2024-09-11

Abstract:We consider the Weak Galerkin finite element approximation of the Singularly Perturbed Biharmonic elliptic problem on a unit square domain with clamped boundary conditions. Shishkin mesh is used for domain discretization as the solution exhibits boundary layers near the domain boundary. Error estimates in the equivalent $H^{2}-$ norm have been established and the uniform convergence of the proposed method has been proved. Numerical examples are presented corroborating our theoretical findings.

Numerical Analysis
Interior Penalty Mixed Finite Element Methods of Any Order in Any Dimension for Linear Elasticity with Strongly Symmetric Stress Tensor

Shuonan Wu,Shihua Gong,Jinchao Xu

DOI: https://doi.org/10.1142/s0218202517500567

2017-01-01

Mathematical Models and Methods in Applied Sciences

Abstract:We propose two classes of mixed finite elements for linear elasticity of any order, with interior penalty for nonconforming symmetric stress approximation. One key point of our method is to introduce some appropriate nonconforming face-bubble spaces based on the local decomposition of discrete symmetric tensors, with which the stability can be easily established. We prove the optimal L-2-error estimate for displacement and optimal Hh(div) error estimate for stress by adding an interior penalty term. The elements are easy to be implemented thanks to the explicit formulations of its basis functions. Moreover, the method can be applied to arbitrary simplicial grids for any spatial dimension in a unified fashion. Numerical tests for both 2D and 3D are provided to validate our theoretical results.
A study on the conditioning of finite element equations with arbitrary anisotropic meshes via a density function approach

Lennard Kamenski,Weizhang Huang

DOI: https://doi.org/10.4208/jms.v47n2.14.02

2013-09-18

Abstract:The linear finite element approximation of a general linear diffusion problem with arbitrary anisotropic meshes is considered. The conditioning of the resultant stiffness matrix and the Jacobi preconditioned stiffness matrix is investigated using a density function approach proposed by Fried in 1973. It is shown that the approach can be made mathematically rigorous for general domains and used to develop bounds on the smallest eigenvalue and the condition number that are sharper than existing estimates in one and two dimensions and comparable in three and higher dimensions. The new results reveal that the mesh concentration near the boundary has less influence on the condition number than the mesh concentration in the interior of the domain. This is especially true for the Jacobi preconditioned system where the former has little or almost no influence on the condition number. Numerical examples are presented.

Numerical Analysis
Weak discrete maximum principle of isoparametric finite element methods in curvilinear polyhedra

Buyang Li,Weifeng Qiu,Yupei Xie,Wenshan Yu

2023-03-16

Abstract:The weak maximum principle of the isoparametric finite element method is proved for the Poisson equation under the Dirichlet boundary condition in a (possibly concave) curvilinear polyhedral domain with edge openings smaller than $\pi$, which include smooth domains and smooth deformations of convex polyhedra. The proof relies on the analysis of a dual elliptic problem with a discontinuous coefficient matrix arising from the isoparametric finite elements. Therefore, the standard $H^2$ elliptic regularity which is required in the proof of the weak maximum principle in the literature does not hold for this dual problem. To overcome this difficulty, we have decomposed the solution into a smooth part and a nonsmooth part, and estimated the two parts by $H^2$ and $W^{1,p}$ estimates, respectively. As an application of the weak maximum principle, we have proved a maximum-norm best approximation property of the isoparametric finite element method for the Poisson equation in a curvilinear polyhedron. The proof contains non-trivial modifications of Schatz's argument due to the non-conformity of the iso-parametric finite elements, which requires us to construct a globally smooth flow map which maps the curvilinear polyhedron to a perturbed larger domain on which we can establish the $W^{1,\infty}$ regularity estimate of the Poisson equation uniformly with respect to the perturbation.

Numerical Analysis
Anisotropic mesh adaptation for high-order finite elements spaces with the log-simplex method. Application to discontinuous Galerkin methods

Olivier Coulaud,Adrien Loseille,Pierre Schrooyen

DOI: https://doi.org/10.1016/j.jcp.2024.112774

IF: 4.645

2024-01-24

Journal of Computational Physics

Abstract:In this article, a high-order solution-based mesh adaptation method is investigated. This later, which is called the log-simplex method, relies on the approximation of high-order differential form of the solution by a dedicated metric-field. The method is a natural extension of Hessian-based methods to the high-order case.

physics, mathematical,computer science, interdisciplinary applications
A study on anisotropic mesh adaptation for finite element approximation of eigenvalue problems with anisotropic diffusion operators

Jingyue Wang,Weizhang Huang

DOI: https://doi.org/10.1137/140958554

2014-02-25

Abstract:Anisotropic mesh adaptation is studied for the linear finite element solution of eigenvalue problems with anisotropic diffusion operators. The M-uniform mesh approach is employed with which any nonuniform mesh is characterized mathematically as a uniform one in the metric specified by a metric tensor. Bounds on the error in the computed eigenvalues are established for quasi-M-uniform meshes. Numerical examples arising from the Laplace-Beltrami operator on parameterized surfaces, nonlinear diffusion, thermal diffusion in a magnetic field in plasma physics, and the Laplacian operator on an L-shape domain are presented. Numerical results show that anisotropic adaptive meshes can lead to more accurate computed eigenvalues than uniform or isotropic adaptive meshes. They also confirm the second order convergence of the error that is predicted by the theoretical analysis. The effects of approximation of curved boundaries on the computation of eigenvalue problems is also studied in two dimensions. It is shown that the initial mesh used to define the geometry of the physical domain should contain at least \sqrt{N} boundary points to keep the effects of boundary approximation at the level of the error of the finite element approximation, where N is the number of the elements in the final adaptive mesh. Only about N^{1/3} boundary points in the initial mesh are needed for boundary value problems. This implies that the computation of eigenvalue problems is more sensitive to the boundary approximation than that of boundary value problems.

Numerical Analysis
Anisotropic modified Crouzeix-Raviart finite element method for the stationary Navier-Stokes equation

Hiroki Ishizaka

2024-09-14

Abstract:We studied an anisotropic modified Crouzeix--Raviart finite element method for the rotational form of a stationary incompressible Navier--Stokes equation with large irrotational body forces. We present an anisotropic $H^1$ error estimate for the velocity of the modified Crouzeix--Raviart finite element method for the Navier--Stokes equation. The modified Crouzeix--Raviart finite element scheme was obtained using a lifting operator that mapped the velocity test functions to $H(÷;\Omega)$-conforming finite element spaces. Because no shape-regularity mesh conditions are imposed, anisotropic meshes can be used for the analysis. The core idea of the proof involves using the relation between the Raviart--Thomas and Crouzeix--Raviart finite element spaces. Furthermore, we present a discrete Sobolev inequality under semi-regular mesh conditions to estimate the stability of the proposed method, and confirm the results obtained through numerical experiments.

Numerical Analysis
Sharp Penalty Term and Time Step Bounds for the Interior Penalty Discontinuous Galerkin Method for Linear Hyperbolic Problems

S. Geevers,J.J.W. van der Vegt

DOI: https://doi.org/10.1137/16M1091290

2018-07-06

Abstract:We present sharp and sufficient bounds for the interior penalty term and time step size to ensure stability of the Symmetric Interior Penalty Discontinuous Galerkin (SIPDG) method combined with an explicit time-stepping scheme. These conditions hold for generic meshes, including unstructured non-conforming heterogeneous meshes of mixed element types, and apply to a large class of linear hyperbolic problems, including the acoustic wave equation, the (an)isotropic elastic wave equations and Maxwell's equations. The penalty term bounds are computed element-wise, while bounds for the time step size are computed at weighted submeshes requiring only a small number of elements and faces. Numerical results illustrate the sharpness of these bounds.

Numerical Analysis
Anisotropic Adaptive Finite Elements for a p -Laplacian Problem

Paride Passelli,Marco Picasso

DOI: https://doi.org/10.1515/cmam-2022-0205

2024-06-26

Computational Methods in Applied Mathematics

Abstract:The p -Laplacian problem - ∇ ⋅ ( ( μ + | ∇ ⁡ u | p - 2 ) ⁢ ∇ ⁡ u ) = f is considered, where μ is a given positive number. An anisotropic a posteriori residual-based error estimator is presented. The error estimator is shown to be equivalent, up to higher order terms, to the error in a quasi-norm. The involved constants being independent of μ, the solution, the mesh size and aspect ratio. An adaptive algorithm is proposed and numerical results are presented when p = 3 . From this model problem, we propose a simplified error estimator and use it in the framework of an industrial application, namely a nonlinear Navier–Stokes problem arising from aluminium electrolysis.

mathematics, applied
Anisotropic finite elements with high aspect ratio for an Asymptotic Preserving method for highly anisotropic elliptic equation

Jacek Narski

DOI: https://doi.org/10.48550/arXiv.1302.4269

2013-02-18

Abstract:The concern of this work is the generalization of an Asymptotic Preserving method for the highly anisotropic elliptic equations presented in [P. Degond, A. Lozinski, J. Narski, and C. Negulescu. An asymptotic-preserving method for highly anisotropic elliptic equations based on a micro-macro decomposition. J. Comput. Phys., 231(7):2724{2740, 2012]. The limitations of the method introduced there in are omitted by the introduction of a stabilization term in the Asymptotic Reformulation. Furthermore, anisotropic error indicators and mesh adaptation algorithms are proposed and tested allowing to reduce considerably the number of mesh points required to achieve prescribed precision. Reported meshes have maximum aspect ratio greater than 500.

Numerical Analysis
A Posteriori Error Analysis for Finite Element Solution of One-Dimensional Elliptic Differential Equations Using Equidistributing Meshes

Yinnian He,Weizhang Huang

DOI: https://doi.org/10.1016/j.cam.2015.10.033

IF: 2.872

2015-01-01

Journal of Computational and Applied Mathematics

Abstract:The paper is concerned with the adaptive linear finite element solution of linear one-dimensional elliptic differential equations using equidistributing meshes. A strategy is developed for defining such meshes based on a residual-based a posteriori error estimate. The mesh and the finite element solution are determined by the coupled system of the finite element equation and the equidistribution relation (i.e., the mesh equation). An iterative algorithm is proposed for solving this system for the mesh and the finite element solution. The existence of an equidistributing mesh is proven for a given sufficiently large number of the points with help of a result on the continuous dependence of the finite element solution on the mesh, which is also established in the current work. Error bounds for the finite element solution are obtained for the equidistributing and quasi-equidistributing meshes. They show that adaptive meshes can lead to more accurate solutions than a uniform mesh and it is unnecessary to compute the equidistribution relation accurately for the equidistributing meshes. The departure from the equidistributing meshes has only a mild effect on the finite element error. Numerical examples are given to illustrate the convergence of the iterative algorithm and the theoretical findings.
Very weak finite element methods: discretisation and applications

Douglas Ramalho Queiroz Pacheco

DOI: https://doi.org/10.1108/ec-10-2023-0699

2024-03-21

Engineering Computations

Abstract:Purpose This study aims to propose and numerically assess different ways of discretising a very weak formulation of the Poisson problem. Design/methodology/approach We use integration by parts twice to shift smoothness requirements to the test functions, thereby allowing low-regularity data and solutions. Findings Various conforming discretisations are presented and tested, with numerical results indicating good accuracy and stability in different types of problems. Originality/value This is one of the first articles to propose and test concrete discretisations for very weak variational formulations in primal form. The numerical results, which include a problem based on real MRI data, indicate the potential of very weak finite element methods for tackling problems with low regularity.

mechanics,computer science, interdisciplinary applications,engineering, multidisciplinary,mathematics
An Interface Penalty Finite Element Method for Elliptic Interface Problems on Piecewise Meshes

Xiaoxiao He,Weibing Deng,Haijun Wu

DOI: https://doi.org/10.1016/j.cam.2019.112473

IF: 2.872

2019-01-01

Journal of Computational and Applied Mathematics

Abstract:An interface penalty finite element method (IPFEM) is proposed for elliptic interface problems, which allows to use different meshes in different subdomains separated by the interface. The transmission conditions across the interface are treated by the Nitsche’s method (or penalty technique) with some harmonic-weighted averages. Both symmetric IPFEM and non-symmetric IPFEM are analyzed. Optimal order error estimates in energy norm and L2-norm and the flux error estimate in L2-norm are proved for the symmetric IPFEM. In particular, the relative error estimate in energy norm and the flux error estimate in L2-norm are independent of the ratio of mesh sizes and the contrast of the discontinuous coefficients across the interface. Error estimates for the non-symmetric IPFEM are also obtained. Furthermore, optimal mesh sizes with respect to the error estimate in energy norm and the flux error estimate in L2-norm are discussed, respectively. Numerical examples are also provided to confirm the theoretical results and show that the IPFEM on a mesh with optimal mesh sizes may give much better approximations than the one on the quasi-uniform mesh with the same number of nodal points.
Nonconforming Dirichlet boundary conditions in implicit material point method by means of penalty augmentation

Bodhinanda Chandra,Veronika Singer,Tobias Teschemacher,Roland Wüchner,Antonia Larese

DOI: https://doi.org/10.1007/s11440-020-01123-3

2021-01-24

Acta Geotechnica

Abstract:Abstract In many geomechanics applications, material boundaries are subjected to large displacements and deformation. Under these circumstances, the application of boundary conditions using particle methods, such as the material point method (MPM), becomes a challenging task since material boundaries do not coincide with the background mesh. This paper presents a formulation of penalty augmentation to impose nonhomogeneous, nonconforming Dirichlet boundary conditions in implicit MPM. The penalty augmentation is implemented utilizing boundary particles, which can move either according to or independently from the material deformation. Furthermore, releasing contact boundary condition, as well as the capability to accommodate slip boundaries, is introduced in the current work. The accuracy of the proposed method is assessed in both 2D and 3D cases, by convergence analysis reaching the analytical solution and by comparing the results of nonconforming and classical grid-conforming simulations.

engineering, geological
Anisotropic mesh adaptation for finite element solution of Anisotropic Porous Medium Equation

Xianping Li

DOI: https://doi.org/10.1016/j.camwa.2017.08.005

2018-03-01

Abstract:Anisotropic Porous Medium Equation (APME) is developed as an extension of the Porous Medium Equation (PME) for anisotropic porous media. A special analytical solution is derived for APME for time-independent diffusion. Anisotropic mesh adaptation for linear finite element solution of APME is discussed and numerical results for two dimensional examples are presented. The solution errors using anisotropic adaptive meshes show second order convergence.

mathematics, applied

Note on a weakly over-penalised symmetric interior penalty method on anisotropic meshes for the Poisson equation, Ver. 1

Hybrid weakly over-penalised symmetric interior penalty method on anisotropic meshes

Anisotropic weakly over-penalised symmetric interior penalty method for the Stokes equation

An Approach to Numerically Solving the Poisson Equation

Crouzeix–Raviart and Raviart–Thomas finite-element error analysis on anisotropic meshes violating the maximum-angle condition

Anisotropic Error Analysis of Weak Galerkin finite element method for Singularly Perturbed Biharmonic Problems

Interior Penalty Mixed Finite Element Methods of Any Order in Any Dimension for Linear Elasticity with Strongly Symmetric Stress Tensor

A study on the conditioning of finite element equations with arbitrary anisotropic meshes via a density function approach

Weak discrete maximum principle of isoparametric finite element methods in curvilinear polyhedra

Anisotropic mesh adaptation for high-order finite elements spaces with the log-simplex method. Application to discontinuous Galerkin methods

A study on anisotropic mesh adaptation for finite element approximation of eigenvalue problems with anisotropic diffusion operators

Anisotropic modified Crouzeix-Raviart finite element method for the stationary Navier-Stokes equation

Sharp Penalty Term and Time Step Bounds for the Interior Penalty Discontinuous Galerkin Method for Linear Hyperbolic Problems

Anisotropic Adaptive Finite Elements for a p -Laplacian Problem

Anisotropic finite elements with high aspect ratio for an Asymptotic Preserving method for highly anisotropic elliptic equation

A Posteriori Error Analysis for Finite Element Solution of One-Dimensional Elliptic Differential Equations Using Equidistributing Meshes

Very weak finite element methods: discretisation and applications

An Interface Penalty Finite Element Method for Elliptic Interface Problems on Piecewise Meshes

Nonconforming Dirichlet boundary conditions in implicit material point method by means of penalty augmentation

Anisotropic mesh adaptation for finite element solution of Anisotropic Porous Medium Equation