Frenet Immersed Finite Element Spaces on Triangular Meshes

Tao Lin; Xu Zhang; Yuanhui Lin

arxiv: 2512.23238 · v2 · pith:JQ3CPHKYnew · submitted 2025-12-29 · 🧮 math.NA · cs.NA

Frenet Immersed Finite Element Spaces on Triangular Meshes

Yuanhui Lin , Xu Zhang , Tao Lin This is my paper

Pith reviewed 2026-05-25 07:11 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords immersed finite elementFrenet-Serret mappingtriangular mesheselliptic interface problemsdiscontinuous Galerkinjump conditionsgeometry-conforminghigh-order approximation

0 comments

The pith

Frenet-Serret mapping transforms curved interfaces into straight lines inside triangular elements, allowing exact jump conditions in high-order immersed finite element spaces.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs geometry-conforming immersed finite element spaces on triangular meshes for elliptic interface problems by applying a Frenet-Serret mapping to each cut element. This mapping straightens the smooth interface curve into a line segment, so that the jump conditions can be imposed exactly rather than approximately. Three variants are presented: monomial bases, orthogonal polynomial bases, and reconstructed bases that improve mass-matrix conditioning. These spaces are inserted into an interior-penalty discontinuous Galerkin scheme, where numerical tests recover the optimal convergence rates in both the H1 and L2 norms. A reader cares because the construction removes the need to align the mesh with the interface while preserving high-order accuracy on simple triangular grids.

Core claim

The central claim is that a Frenet-Serret mapping applied element-wise turns a smooth interface curve into a straight line segment inside each cut triangle; this geometry-conforming transformation permits three families of high-order IFE bases (monomial, orthogonal-polynomial, and reconstructed) whose approximation properties, when used inside an interior-penalty DG method, produce optimal H1 and L2 error rates for elliptic interface problems on unfitted triangular meshes.

What carries the argument

The Frenet-Serret mapping that straightens the interface curve into a line segment inside each cut element so jump conditions can be imposed exactly.

If this is right

The three IFE constructions deliver high-order approximation on unfitted triangular meshes for elliptic interface problems.
Optimal H1 and L2 convergence rates are obtained when the spaces are used inside interior-penalty discontinuous Galerkin schemes.
Reconstructed bases improve the conditioning of the local mass matrices relative to the monomial construction.
The triangular-mesh framework extends the earlier rectangular-mesh Frenet-IFE method while retaining the exact jump-condition property.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same mapping idea could be applied to other unfitted-mesh methods such as cut finite elements or extended finite elements on triangles.
Because the mapping is local to each cut element, the approach may combine naturally with adaptive mesh refinement near the interface.
Theoretical a-priori error bounds could be derived by transferring standard approximation theory across the mapping, once the mapping's smoothness properties are quantified.

Load-bearing premise

The interface curve must be smooth enough that the Frenet-Serret mapping is well-defined and turns the curve into a straight line segment inside every cut element.

What would settle it

Numerical tests on a sequence of successively refined triangular meshes that cut a smooth curved interface would falsify the claim if the observed H1 or L2 convergence rates fall below the predicted optimal orders.

read the original abstract

In this paper, we develop geometry-conforming immersed finite element (IFE) spaces on triangular meshes for elliptic interface problems. The construction is built on a Frenet-Serret mapping that transforms a smooth interface curve into a straight line, so that the interface jump conditions can be imposed exactly. Extending the framework of [9] from rectangular meshes to triangular meshes, we introduce three types of high-order Frenet-IFE constructions: an initial construction using monomial bases, a general construction using orthogonal polynomials, and reconstructed IFE bases designed to improve the conditioning of the mass matrix. The approximation properties of these new IFE spaces are investigated through extensive numerical experiments. We also incorporate the new IFE spaces into interior penalty discontinuous Galerkin methods for solving elliptic interface problems, and demonstrate optimal convergence rates in $H^1$- and $L^2$- norms.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Extends Frenet-IFE from rectangles to triangles with three basis constructions and reports optimal IPDG convergence rates, but the exact straightening claim for the mapping needs explicit verification.

read the letter

Colleague, the main new material is the move of the Frenet-Serret IFE construction to triangular meshes. They define three high-order versions: one built on monomial bases, one using orthogonal polynomials, and a reconstructed set meant to improve mass-matrix conditioning. These are then used inside an interior-penalty DG scheme for elliptic interface problems, with numerical tests showing optimal H1 and L2 rates under mesh refinement. That extension and the concrete basis options are the actual advance over the rectangular-mesh reference. The experiments are the part that lands: they check approximation properties across the three constructions and obtain the expected convergence, which is practical evidence for unstructured-mesh work. The soft spot is the mapping itself. The method stands on the claim that the local Frenet-Serret frame turns the curved interface into an exact straight segment inside every cut triangle, letting the jump conditions be imposed without approximation error. The abstract asserts this for smooth curves, but the stress-test concern is fair: without seeing the local derivation, it is not obvious that the mapping stays exact on general triangles rather than becoming approximate near vertices or under only C2 regularity. If the full paper supplies a clear geometric argument or explicit construction that confirms exactness, the rates can be attributed to the geometry-conforming property; otherwise they may partly reflect other factors. The smoothness assumption is standard and not a major issue by itself. This paper is aimed at people already working on immersed or DG methods for interface problems who need triangular-mesh options. A reader focused on high-order IFE constructions or conditioning fixes will get direct value from the three variants and the reported rates. It deserves peer review because the new constructions and the numerical evidence are specific enough to warrant referee scrutiny, even if the mapping details may need tightening in revision.

Referee Report

2 major / 0 minor

Summary. The paper develops geometry-conforming immersed finite element (IFE) spaces on triangular meshes for elliptic interface problems. The construction relies on a Frenet-Serret mapping that transforms a smooth interface curve into a straight line segment inside each cut element, enabling exact imposition of jump conditions. It extends prior rectangular-mesh work by introducing three high-order Frenet-IFE constructions (monomial bases, orthogonal polynomials, and reconstructed bases), investigates their approximation properties via numerical experiments, and incorporates the spaces into interior-penalty DG methods, reporting optimal H^1 and L^2 convergence rates.

Significance. If the Frenet-Serret mapping is shown to produce an exact linear image of the interface arc in arbitrary cut triangles, the work would provide a meaningful extension of IFE methods to unstructured triangular meshes while preserving the geometry-conforming property that supports exact jump enforcement. The numerical demonstration of optimal rates is a positive indicator of practical utility, though the absence of a supporting theoretical analysis limits the assessed impact.

major comments (2)

[Abstract] Abstract: the central claim that the Frenet-Serret mapping 'transforms a smooth interface curve into a straight line' so that jump conditions can be imposed exactly is load-bearing for attributing optimal convergence to the geometry-conforming construction. No derivation is supplied showing that the local Frenet frame on a general cut triangle maps the curved arc to a precise straight-line segment (as opposed to an approximation), which directly affects whether the observed rates can be credited to the claimed exactness rather than to the specific numerical tests.
[Abstract] The abstract states that 'numerical experiments demonstrate optimal convergence,' yet supplies no information on the mesh-refinement strategy, the precise interface regularity assumed, or how errors are computed in the IPDG experiments. This information is required to evaluate whether the reported H^1 and L^2 rates are robust and truly optimal.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed comments on our manuscript. We provide point-by-point responses below and indicate where revisions will be made to address the concerns.

read point-by-point responses

Referee: [Abstract] Abstract: the central claim that the Frenet-Serret mapping 'transforms a smooth interface curve into a straight line' so that jump conditions can be imposed exactly is load-bearing for attributing optimal convergence to the geometry-conforming construction. No derivation is supplied showing that the local Frenet frame on a general cut triangle maps the curved arc to a precise straight-line segment (as opposed to an approximation), which directly affects whether the observed rates can be credited to the claimed exactness rather than to the specific numerical tests.

Authors: The Frenet-Serret mapping is constructed locally in each cut element using the Frenet frame at a point on the interface arc, which by definition aligns the tangent vector and maps the curve to the local x-axis in the transformed coordinates. However, we agree that an explicit derivation confirming the exact straight-line image for arbitrary cut triangles would strengthen the presentation. In the revised manuscript, we will add a subsection or appendix providing this derivation based on the properties of the Frenet-Serret formulas. revision: yes
Referee: [Abstract] The abstract states that 'numerical experiments demonstrate optimal convergence,' yet supplies no information on the mesh-refinement strategy, the precise interface regularity assumed, or how errors are computed in the IPDG experiments. This information is required to evaluate whether the reported H^1 and L^2 rates are robust and truly optimal.

Authors: While the abstract is necessarily concise, the requested details on mesh refinement (uniform h-refinement), interface regularity (C^2 smooth curves), and error computation (standard L2 and H1 norms via quadrature) are fully specified in Sections 4 and 5 of the manuscript. To improve clarity, we will revise the abstract to briefly note the use of uniform triangular mesh refinement and the assumed interface smoothness. revision: partial

Circularity Check

0 steps flagged

No significant circularity; construction and convergence rest on independent geometric mapping and numerical tests

full rationale

The derivation introduces a Frenet-Serret mapping to straighten the interface inside cut triangles, defines three new IFE constructions (monomial, orthogonal, reconstructed), and validates approximation properties plus optimal IPDG convergence rates exclusively through numerical experiments. No equation or claim reduces by construction to a fitted parameter, self-definition, or unverified self-citation chain; the extension of [9] supplies only the rectangular-mesh starting point while the triangular-mesh mapping, basis choices, and observed rates are externally falsifiable via the reported experiments. The smoothness assumption on the interface is stated explicitly and does not smuggle in the target result.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the domain assumption that the interface is smooth and on standard Sobolev-space theory for elliptic problems; no free parameters or new entities are introduced in the abstract.

axioms (1)

domain assumption The interface curve is sufficiently smooth for the Frenet-Serret mapping to be well-defined and to map the curve to a straight line segment.
Invoked to justify exact imposition of jump conditions after the mapping.

pith-pipeline@v0.9.0 · 5674 in / 1255 out tokens · 54798 ms · 2026-05-25T07:11:58.814184+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

A Priori Error Analysis of a High-Order Selective Discontinuous Galerkin Method for Elliptic Interface Problems
math.NA 2026-05 unverdicted novelty 6.0

A high-order selective DG method with a new hybrid IFE space is introduced for elliptic interface problems on unfitted meshes, with proofs of optimal approximation, well-posedness, and a priori error estimates in ener...

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages · cited by 1 Pith paper

[1]

Marco Abate and Francesca Tovena,Curves and surfaces, Springer Science & Business Media, 2012

work page 2012
[2]

Abbena, S

E. Abbena, S. Salamon, and A. Gray,Modern differential geometry of curves and surfaces with mathematica, 3 ed., Chapman and Hall/CRC, 2006

work page 2006
[3]

Slimane Adjerid,A study of high-order immersed finite element spaces by pointwise interface conditions on curved interfaces, Comput. Math. Appl.128(2022), 331–353. MR 4512459

work page 2022
[4]

Methods Appl

Slimane Adjerid, Ivo Babuška, Ruchi Guo, and Tao Lin,An enriched immersed finite element method for inter- face problems with nonhomogeneous jump conditions, Comput. Methods Appl. Mech. Engrg.404(2023), Paper No. 115770, 37. MR 4518020

work page 2023
[5]

Slimane Adjerid, Mohamed Ben-Romdhane, and Tao Lin,Higher degree immersed finite element spaces con- structed according to the actual interface, Comput. Math. Appl.75(2018), no. 6, 1868–1881. MR 3775090

work page 2018
[6]

Slimane Adjerid, Ruchi Guo, and Tao Lin,High degree immersed finite element spaces by a least squares method, Int. J. Numer. Anal. Model.14(2017), no. 4-5, 604–626

work page 2017
[7]

Methods Appl

Slimane Adjerid, Tao Lin, and Haroun Meghaichi,A high order geometry conforming immersed finite ele- ment for elliptic interface problems, Comput. Methods Appl. Mech. Engrg.420(2024), Paper No. 116703, 21. MR 4682175

work page 2024
[8]

,Construction of basis functions for the geometry conforming immersed finite element method, arXiv preprint (2025), arXiv:2510.12018

work page arXiv 2025
[9]

Methods Appl

,The Frenet immersed finite element method for elliptic interface problems: an error analysis, Comput. Methods Appl. Mech. Engrg.438(2025), Paper No. 117829, 19. MR 4867707

work page 2025
[10]

Yuan Chen and Xu Zhang,A P 2-P1 partially penalized immersed finite element method for stokes interface prob- lems., Int. J. Numer. Anal. Model.18(2021), no. 1, 120–141

work page 2021
[11]

,An immersed c 0 interior penalty finite element method for biharmonic interface problems., submitted (2025), arXiv:2509.12555

work page arXiv 2025
[12]

Ruchi Guo and Tao Lin,A group of immersed finite element spaces for elliptic interface problems, IMA J. Numer. Anal.39(2019), no. 1, 482–511

work page 2019
[13]

Ruchi Guo and Tao Lin,A higher degree immersed finite element method based on a Cauchy extension for elliptic interface problems, SIAM J. Numer. Anal.57(2019), no. 4, 1545–1573. MR 3977109

work page 2019
[14]

Ruchi Guo and Xu Zhang,Solving three-dimensional interface problems with immersed finite elements: A-priori error analysis, J. Comput. Phys.441(2021), 110445. MR 4265655

work page 2021
[15]

Methods Partial Differential Equations24(2008), no

Xiaoming He, Tao Lin, and Yanping Lin,Approximation capability of a bilinear immersed finite element space, Numer. Methods Partial Differential Equations24(2008), no. 5, 1265–1300. MR 2427191 (2009c:65307)

work page 2008
[16]

Hesthaven and Tim Warburton,Nodal discontinuous Galerkin methods, Texts in Applied Mathematics, vol

Jan S. Hesthaven and Tim Warburton,Nodal discontinuous Galerkin methods, Texts in Applied Mathematics, vol. 54, Springer, New York, 2008, Algorithms, analysis, and applications. MR 2372235

work page 2008
[17]

Kafafy, T

R. Kafafy, T . Lin, Y. Lin, and J. Wang,Three-dimensional immersed finite element methods for electric field sim- ulation in composite materials, Internat. J. Numer. Methods Engrg.64(2005), no. 7, 940–972. MR 2172214 (2006e:78008)

work page 2005
[18]

Z. Li, T . Lin, Y. Lin, and R. C. Rogers,An immersed finite element space and its approximation capability, Numer. Methods Partial Differential Equations20(2004), no. 3, 338–367. MR 2046521 (2005f:65153)

work page 2004
[19]

1, 61–98

Zhilin Li, Tao Lin, and Xiaohui Wu,New Cartesian grid methods for interface problems using finite element formulation, Numerische Mathematik96(2003), no. 1, 61–98

work page 2003
[20]

Tao Lin, Yanping Lin, and Xu Zhang,Partially penalized immersed finite element methods for elliptic interface problems, SIAM J. Numer. Anal.53(2015), no. 2, 1121–1144. MR 3338673

work page 2015
[21]

Tao Lin, Dongwoo Sheen, and Xu Zhang,A nonconforming immersed finite element method for elliptic interface problems, J. Sci. Comput.79(2019), no. 1, 442–463

work page 2019
[22]

Barret O’Neill,Elementary differential geometry, second ed., Academic Press, 2006

work page 2006
[23]

FR35, SIAM, Philadelphia, 2008

Beatrice Rivière,Discontinuous galerkin methods for solving elliptic and parabolic equations, Frontiers in Ap- plied Mathematics, vol. FR35, SIAM, Philadelphia, 2008. 20 TAO LIN, YUANHUI LIN, AND XU ZHANG

work page 2008
[24]

Qiao Zhuang, Zhongqiang Zhang, Marcus Sarkis, and Tao Lin,Higher-degree rectangular immersed finite ele- ments discontinuous Galerkin methods for elliptic interface problems, J. Comput. Appl. Math.476(2026), Paper No. 117126, 26. MR 4968689 DEPARTMENT OFMATHEMATICS, VIRGINIATECH, BLACKSBURG, VA 24060, USA. Email address:tlin@vt.edu SCHOOL OFMATHEMATICS, S...

work page 2026

[1] [1]

Marco Abate and Francesca Tovena,Curves and surfaces, Springer Science & Business Media, 2012

work page 2012

[2] [2]

Abbena, S

E. Abbena, S. Salamon, and A. Gray,Modern differential geometry of curves and surfaces with mathematica, 3 ed., Chapman and Hall/CRC, 2006

work page 2006

[3] [3]

Slimane Adjerid,A study of high-order immersed finite element spaces by pointwise interface conditions on curved interfaces, Comput. Math. Appl.128(2022), 331–353. MR 4512459

work page 2022

[4] [4]

Methods Appl

Slimane Adjerid, Ivo Babuška, Ruchi Guo, and Tao Lin,An enriched immersed finite element method for inter- face problems with nonhomogeneous jump conditions, Comput. Methods Appl. Mech. Engrg.404(2023), Paper No. 115770, 37. MR 4518020

work page 2023

[5] [5]

Slimane Adjerid, Mohamed Ben-Romdhane, and Tao Lin,Higher degree immersed finite element spaces con- structed according to the actual interface, Comput. Math. Appl.75(2018), no. 6, 1868–1881. MR 3775090

work page 2018

[6] [6]

Slimane Adjerid, Ruchi Guo, and Tao Lin,High degree immersed finite element spaces by a least squares method, Int. J. Numer. Anal. Model.14(2017), no. 4-5, 604–626

work page 2017

[7] [7]

Methods Appl

Slimane Adjerid, Tao Lin, and Haroun Meghaichi,A high order geometry conforming immersed finite ele- ment for elliptic interface problems, Comput. Methods Appl. Mech. Engrg.420(2024), Paper No. 116703, 21. MR 4682175

work page 2024

[8] [8]

,Construction of basis functions for the geometry conforming immersed finite element method, arXiv preprint (2025), arXiv:2510.12018

work page arXiv 2025

[9] [9]

Methods Appl

,The Frenet immersed finite element method for elliptic interface problems: an error analysis, Comput. Methods Appl. Mech. Engrg.438(2025), Paper No. 117829, 19. MR 4867707

work page 2025

[10] [10]

Yuan Chen and Xu Zhang,A P 2-P1 partially penalized immersed finite element method for stokes interface prob- lems., Int. J. Numer. Anal. Model.18(2021), no. 1, 120–141

work page 2021

[11] [11]

,An immersed c 0 interior penalty finite element method for biharmonic interface problems., submitted (2025), arXiv:2509.12555

work page arXiv 2025

[12] [12]

Ruchi Guo and Tao Lin,A group of immersed finite element spaces for elliptic interface problems, IMA J. Numer. Anal.39(2019), no. 1, 482–511

work page 2019

[13] [13]

Ruchi Guo and Tao Lin,A higher degree immersed finite element method based on a Cauchy extension for elliptic interface problems, SIAM J. Numer. Anal.57(2019), no. 4, 1545–1573. MR 3977109

work page 2019

[14] [14]

Ruchi Guo and Xu Zhang,Solving three-dimensional interface problems with immersed finite elements: A-priori error analysis, J. Comput. Phys.441(2021), 110445. MR 4265655

work page 2021

[15] [15]

Methods Partial Differential Equations24(2008), no

Xiaoming He, Tao Lin, and Yanping Lin,Approximation capability of a bilinear immersed finite element space, Numer. Methods Partial Differential Equations24(2008), no. 5, 1265–1300. MR 2427191 (2009c:65307)

work page 2008

[16] [16]

Hesthaven and Tim Warburton,Nodal discontinuous Galerkin methods, Texts in Applied Mathematics, vol

Jan S. Hesthaven and Tim Warburton,Nodal discontinuous Galerkin methods, Texts in Applied Mathematics, vol. 54, Springer, New York, 2008, Algorithms, analysis, and applications. MR 2372235

work page 2008

[17] [17]

Kafafy, T

R. Kafafy, T . Lin, Y. Lin, and J. Wang,Three-dimensional immersed finite element methods for electric field sim- ulation in composite materials, Internat. J. Numer. Methods Engrg.64(2005), no. 7, 940–972. MR 2172214 (2006e:78008)

work page 2005

[18] [18]

Z. Li, T . Lin, Y. Lin, and R. C. Rogers,An immersed finite element space and its approximation capability, Numer. Methods Partial Differential Equations20(2004), no. 3, 338–367. MR 2046521 (2005f:65153)

work page 2004

[19] [19]

1, 61–98

Zhilin Li, Tao Lin, and Xiaohui Wu,New Cartesian grid methods for interface problems using finite element formulation, Numerische Mathematik96(2003), no. 1, 61–98

work page 2003

[20] [20]

Tao Lin, Yanping Lin, and Xu Zhang,Partially penalized immersed finite element methods for elliptic interface problems, SIAM J. Numer. Anal.53(2015), no. 2, 1121–1144. MR 3338673

work page 2015

[21] [21]

Tao Lin, Dongwoo Sheen, and Xu Zhang,A nonconforming immersed finite element method for elliptic interface problems, J. Sci. Comput.79(2019), no. 1, 442–463

work page 2019

[22] [22]

Barret O’Neill,Elementary differential geometry, second ed., Academic Press, 2006

work page 2006

[23] [23]

FR35, SIAM, Philadelphia, 2008

Beatrice Rivière,Discontinuous galerkin methods for solving elliptic and parabolic equations, Frontiers in Ap- plied Mathematics, vol. FR35, SIAM, Philadelphia, 2008. 20 TAO LIN, YUANHUI LIN, AND XU ZHANG

work page 2008

[24] [24]

Qiao Zhuang, Zhongqiang Zhang, Marcus Sarkis, and Tao Lin,Higher-degree rectangular immersed finite ele- ments discontinuous Galerkin methods for elliptic interface problems, J. Comput. Appl. Math.476(2026), Paper No. 117126, 26. MR 4968689 DEPARTMENT OFMATHEMATICS, VIRGINIATECH, BLACKSBURG, VA 24060, USA. Email address:tlin@vt.edu SCHOOL OFMATHEMATICS, S...

work page 2026