Off-site Enforcement of Natural Conditions on Smooth Boundaries for Finite Elements Upon Fitted Straight-edged Triangular Meshes

Jul 16, 2026

Off-site enforcement of natural boundary conditions addresses reduced accuracy when smooth domains are approximated by fitted straight-edged triangular meshes. Prior work suggested curved-edge elements to recover the intended theoretical approximation order. This study instead retains straight-edged elements and restores accuracy by adding terms to the bilinear form so boundary conditions on the approximating polytope better represent the true natural conditions. The method is applied to triangular Lagrange finite elements and supported by a rigorous reliability analysis for reaction–diffusion equations, with numerical experiments validating results.

LanguageEnglish

CategoryResearch & Report

FormatPDF

Pages25

Reading Time1 h 3 min

Abstract

Methodological background

Outline and section plan

arXiv :2607 .07963v1 [math .NA] 8 Jul 2026 Off-site enforcement of natural conditions on smooth boundaries for finite elements upon fitted straight-edged triangular meshes Vitoriano Ruas 1 * 1 Institut Jean Le Rond d’Alembert, CNRS UMR 7190, Sorbonne Universit, Paris, France. e-mail: [vitoriano.ruas@upmc.fr](vitoriano.ruas@upmc.fr) Abstract A few decades ago some possible remedies to an inaccurate enforcement of Neumann or Robin conditions prescribed on the boundary of a smooth domain, owing to the approximation of a curved domain by the union of straight-edged triangles or tetrahedra in a fitted mesh, were addressed in the literature. By that time authors such as Barrett & Elliott [5] advocated the use of elements with a single curved edge or face fitting the true boundary not only at two or three vertexes, but also at additional points on those curves or curved surfaces, so as to define a polynomial surface of a certain type compatible with the theoretical approximation order of the method in use. In this work we adopt a different approach, whose main feature is the use of a fitted mesh consisting of straight-edged elements only. The recovery of lost accuracy due to the domain’s approximation by a polytope is achieved by means of the addition of terms to the bilinear form, which account for natural boundary conditions of the same type to be prescribed on the approximating boundary, though much closer to the true ones. This technique is applied here to the case of triangular Lagrange finite elements, for which we give a rigorous reliability study in the solution of reaction-diffusion equations. Numerical experimentation is supplied in support of the theoretical results. Keywords: Boundary conditions; Curved domains; Finite elements; Lagrange; Modified bilinear form; Neumann; Robin; Second-order elliptic equations; Straight-edged triangles. AMS Subject Classification: 65N30, 74S05, 76M10, 78M10, 80M10 . 1 Methodological background In the framework of the non-affine simplicial finite-element solution of boundary value problems, the author and collaborators addressed a new approach to handle more accurately zero (resp. non zero) DOFs (degrees of freedom) prescribed on the boundary of a smooth N-dimensional domain for N = 2 , 3 ina series of publications finalized within the last ten years. In contrast to other widespread techniques such as the isoparametric method, in this approach both the shape-and test-function spaces consist only of polynomials and the computational domain is the polytope formed by the union of straight-edged N-simplexes of a fitting mesh. But its key point is the use of a trial space of polynomials (resp. manifold) different from the test space, in the sense that for the former prescribed DOFs on the boundary are enforced at their exact locations, while for the latter such values are enforced at shifted locationson the boundary of the approximating polytope instead. We refer to [20], [21][22], [24], [23], [6] and references therein for mote details on this methodology, as applied to different kinds of problems and formulations. It is noteworthy that similar principles have been exploited in parallel by other authors, mostly in connection with the finite volume method (see e.g. [11]) . More recently they were also applied in [26] to the DG (discontinuous Galerkin) method in the two-dimensional case. The reliability of the resulting method restricted to Dirichlet boundary conditions was formally established in [2], by adapting to the * Sorbonne Universit, Campus Pierre et Marie Curie, 4 place jussieu, Couloir 55-65, 4me tage, 75005 Paris, France. DG environment the mathematical analysis of this type of technique provided in [24] for the finite element method. We also observe that in these works such an approach is referred to as the ROD-method, where the acronym stands for reconstruction of off-site data. In contrast, in our own contributions such as [24], we call it a Petrov-Galerki

FAQs

What problem does the off-site enforcement approach aim to solve in finite element boundary treatments?

It targets the loss of accuracy when prescribing Neumann or Robin conditions on a smooth boundary that is approximated by unions of straight-edged triangles or tetrahedra.

How does the paper recover accuracy while using only straight-edged finite elements?

It adds corrective terms to the bilinear form so natural boundary conditions on the approximating polytope match the true ones more closely.

What is the main application and what kind of validation is provided?

The method is applied to triangular Lagrange finite elements for reaction–diffusion equations, with a rigorous reliability study and supporting numerical experiments.

Jul 16, 2026

Upload by Genevieve