2025-11-10T03:15:01.076525

Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method

Adjerid, Lin, Meghaichi
The Frenet apparatus is a new framework for constructing high order geometry-conforming immersed finite element functions for interface problems. In this report, we present a procedure for constructing the local IFE bases in some detail as well as a new approach for constructing orthonormal bases using the singular value decomposition of the local generalized Vandermonde matrix. A sample implementation in MATLAB is provided to showcase the simplicity and extensionability of the framework.
academic

Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method

基本信息

  • 论文ID: 2510.12018
  • 标题: Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method
  • 作者: Slimane Adjerid (Virginia Tech), Tao Lin (Virginia Tech), Haroun Meghaichi (Ohio State University)
  • 分类: math.NA cs.NA
  • 发表时间: October 15, 2025
  • 论文链接: https://arxiv.org/abs/2510.12018

摘要

本文提出了基于Frenet装置的新框架，用于构造解决界面问题的高阶几何一致性浸入式有限元函数。论文详细介绍了构造局部IFE基函数的过程，以及使用局部广义Vandermonde矩阵奇异值分解构造正交归一化基的新方法。提供了MATLAB示例实现，展示了该框架的简洁性和可扩展性。

研究背景与动机

问题背景

本研究要解决椭圆界面问题： −∇⋅(β∇u)=f, on Ω−∪Ω+-\nabla \cdot (\beta\nabla u) = f, \text{ on } \Omega^- \cup \Omega^+ 其中计算域Ω⊂R2\Omega \subset \mathbb{R}^2被界面Γ\Gamma分为两个子域Ω−\Omega^-和Ω+\Omega^+，扩散系数β\beta是分段常数函数。

研究动机

  1. 现有IFE方法的局限性：传统浸入式有限元方法只能近似满足界面跳跃条件
  2. 高阶精度需求：需要构造任意阶多项式的IFE函数
  3. 几何一致性要求：需要精确满足界面跳跃条件的IFE函数
  4. 计算稳定性：需要条件数良好的基函数系统

重要性

界面问题广泛出现在多相流、材料科学、生物医学等领域，高精度数值求解方法对工程应用具有重要意义。

核心贡献

  1. 提出了基于Frenet装置的GC-IFE构造框架：利用微分几何中的Frenet标架简化界面跳跃条件
  2. 开发了两种基函数构造方法：
    • 基于特殊格式的初始构造方法
    • 基于扩展思想的广义构造方法
  3. 提出了两种正交化重构算法：
    • 基于质量矩阵SVD的方法（Approach 1）
    • 基于广义Vandermonde矩阵SVD的方法（Approach 2）
  4. 提供了完整的MATLAB实现：展示了方法的实用性和可扩展性

方法详解

Frenet变换

核心思想是利用Frenet变换将物理坐标(x,y)(x,y)转换为局部坐标(η,ξ)(\eta,\xi)： [ηξ]→PΓ[x(η,ξ)y(η,ξ)]=g(ξ)+ηn(ξ)\begin{bmatrix} \eta \\ \xi \end{bmatrix} \xrightarrow{P_\Gamma} \begin{bmatrix} x(\eta,\xi) \\ y(\eta,\xi) \end{bmatrix} = g(\xi) + \eta n(\xi)

其中g(ξ)g(\xi)是界面曲线的参数化，n(ξ)n(\xi)是法向量。

界面跳跃条件变换

在局部坐标系中，界面跳跃条件变为：

  1. ⟦u^⟧Γ^KF=0\llbracket \hat{u} \rrbracket_{\hat{\Gamma}_{K_F}} = 0
  2. ⟦β^u^η⟧Γ^KF=0\llbracket \hat{\beta}\hat{u}_\eta \rrbracket_{\hat{\Gamma}_{K_F}} = 0
  3. ⟦β^∂j∂ηjL(u^)⟧Γ^KF=0,j=0,1,…,m−2\llbracket \hat{\beta}\frac{\partial^j}{\partial\eta^j}\mathcal{L}(\hat{u}) \rrbracket_{\hat{\Gamma}_{K_F}} = 0, j = 0,1,\ldots,m-2

基函数构造方法

方法1：特殊格式构造

构造形如下式的基函数： ϕ^i,j(η,ξ)=1β^(η,ξ)qj(ηηh)pi(ξ−ξmidξh)\hat{\phi}_{i,j}(\eta,\xi) = \frac{1}{\hat{\beta}(\eta,\xi)}q_j\left(\frac{\eta}{\eta_h}\right)p_i\left(\frac{\xi-\xi_{mid}}{\xi_h}\right)

对于1≤j≤m,0≤i≤m1 \leq j \leq m, 0 \leq i \leq m的情况，这些函数自动满足跳跃条件。

方法2：广义扩展构造

采用更一般的形式： λ^j(η,ξ)={λ^j−(η,ξ)=∑i=1(m+1)2Ci,j−Ri(η,ξ),η<0λ^j+(η,ξ)=∑i=1(m+1)2Ci,j+Ri(η,ξ),η>0\hat{\lambda}_j(\eta,\xi) = \begin{cases} \hat{\lambda}_j^-(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^- R_i(\eta,\xi), & \eta < 0 \\ \hat{\lambda}_j^+(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^+ R_i(\eta,\xi), & \eta > 0 \end{cases}

通过求解线性方程组确定系数矩阵C±C^{\pm}。

正交化重构

Approach 1：质量矩阵SVD

设Mq=V1ΛV1TM_q = V_1\Lambda V_1^T，选择Q1=V1Λ−1/2Q_1 = V_1\Lambda^{-1/2}使得重构后的质量矩阵为单位矩阵。

Approach 2：Vandermonde矩阵SVD

对V~=WV\tilde{V} = \sqrt{W}V进行SVD分解：V~=U2ΣV2T\tilde{V} = U_2\Sigma V_2^T，选择Q2=V2Σ−1Q_2 = V_2\Sigma^{-1}。

理论分析表明Approach 2在数值实现中更稳定，因为避免了对条件数平方的矩阵进行SVD分解。

实验设置

测试问题

考虑域Ω=[−1,1]2\Omega = [-1,1]^2，界面为x2+y2=r02x^2 + y^2 = r_0^2（r0=1/3r_0 = 1/\sqrt{3}），测试函数： u(x,y)={1β+cos⁡(2πr2),r>r01β−cos⁡(2πr2)+cos⁡(2πr02)(1β+−1β−),r<r0u(x,y) = \begin{cases} \frac{1}{\beta^+}\cos(2\pi r^2), & r > r_0 \\ \frac{1}{\beta^-}\cos(2\pi r^2) + \cos(2\pi r_0^2)\left(\frac{1}{\beta^+} - \frac{1}{\beta^-}\right), & r < r_0 \end{cases}

评价指标

  • L2L^2投影误差的收敛阶
  • 质量矩阵条件数
  • 基函数构造的数值稳定性

实验结果

收敛性验证

对于不同多项式阶数m=1,2,3,4m = 1,2,3,4，L2L^2投影误差表现出最优收敛阶：

mmN=16N = 16N=32N = 32N=64N = 64N=128N = 128N=256N = 256收敛阶
18.14E-022.08E-025.23E-031.31E-033.28E-04≈2.0
29.29E-031.19E-031.50E-041.88E-052.36E-06≈3.0
38.80E-045.74E-053.62E-062.27E-071.42E-08≈4.0
47.55E-052.39E-067.49E-082.34E-097.33E-11≈5.0

条件数改善效果

重构前后质量矩阵条件数的比较显示：

  • 初始构造的条件数随多项式阶数指数增长
  • Approach 2能够稳定地将条件数降至1.0
  • Approach 1在高阶时出现数值不稳定

预处理效果

对于系数矩阵AA和A~\tilde{A}，Jacobi预处理器和行归一化预处理器都能显著改善条件数，特别是在高阶和小网格尺寸情况下。

相关工作

IFE方法发展

  • 传统IFE方法：只能近似满足界面条件
  • 扩展IFE方法：基于Cauchy扩展的高阶方法
  • 本文的GC-IFE：精确满足界面条件的几何一致性方法

正交化技术

  • 虚元方法中的正交多项式构造
  • DG方法中的正交基函数
  • 本文将正交化思想引入IFE框架

结论与讨论

主要结论

  1. 理论贡献：建立了基于Frenet装置的GC-IFE理论框架
  2. 算法创新：提出了两种构造方法和两种正交化策略
  3. 数值验证：证明了方法的最优收敛性和数值稳定性
  4. 实用价值：提供了完整的MATLAB实现

局限性

  1. 几何限制：要求界面曲线足够光滑（C3C^3）
  2. 计算复杂度：高阶情况下矩阵组装成本较高
  3. 扩展性：目前仅适用于二维椭圆界面问题

未来方向

  1. 扩展到三维问题和复杂几何
  2. 应用于时间相关的界面问题
  3. 扩展到弹性力学、Stokes系统等
  4. 开发更高效的实现算法

深度评价

优点

  1. 理论严谨：基于微分几何的Frenet装置提供了坚实的理论基础
  2. 方法创新：首次实现精确满足界面跳跃条件的高阶IFE方法
  3. 数值稳定：正交化重构显著改善了计算稳定性
  4. 实现完整：提供了详细的算法描述和MATLAB代码
  5. 验证充分：理论分析和数值实验相互印证

不足

  1. 适用范围有限：目前仅限于二维椭圆界面问题
  2. 计算成本：高阶情况下的计算复杂度分析不够深入
  3. 几何要求严格：对界面曲线的光滑性要求较高
  4. 实际应用案例：缺乏复杂工程问题的应用验证

影响力

  1. 学术价值：为IFE方法提供了新的理论框架和技术路径
  2. 实用价值：为高精度界面问题求解提供了有效工具
  3. 可复现性：完整的代码实现便于其他研究者使用和扩展
  4. 发展潜力：为后续研究奠定了坚实基础

适用场景

  1. 多相流动：不同流体界面的数值模拟
  2. 材料科学：复合材料界面问题的高精度计算
  3. 生物医学：组织界面的传热传质问题
  4. 地球物理：地层界面的波传播问题

参考文献

论文引用了22篇相关文献，涵盖了微分几何、数值分析、有限元方法等多个领域的重要工作，为本研究提供了坚实的理论基础。