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：广义扩展构造

采用更一般的形式：

\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^{\pm}$。 ### 正交化重构 #### Approach 1：质量矩阵SVD 设$M_q = V_1\Lambda V_1^T$，选择$Q_1 = V_1\Lambda^{-1/2}$使得重构后的质量矩阵为单位矩阵。 #### Approach 2：Vandermonde矩阵SVD 对$\tilde{V} = \sqrt{W}V$进行SVD分解：$\tilde{V} = U_2\Sigma V_2^T$，选择$Q_2 = V_2\Sigma^{-1}$。 理论分析表明Approach 2在数值实现中更稳定，因为避免了对条件数平方的矩阵进行SVD分解。 ## 实验设置 ### 测试问题 考虑域$\Omega = [-1,1]^2$，界面为$x^2 + y^2 = r_0^2$（$r_0 = 1/\sqrt{3}$），测试函数： $$u(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}$$ ### 评价指标 - $L^2$投影误差的收敛阶 - 质量矩阵条件数 - 基函数构造的数值稳定性 ## 实验结果 ### 收敛性验证 对于不同多项式阶数$m = 1,2,3,4$，$L^2$投影误差表现出最优收敛阶： | $m$ | $N = 16$ | $N = 32$ | $N = 64$ | $N = 128$ | $N = 256$ | 收敛阶 | |-----|----------|----------|----------|-----------|-----------|--------| | 1 | 8.14E-02 | 2.08E-02 | 5.23E-03 | 1.31E-03 | 3.28E-04 | ≈2.0 | | 2 | 9.29E-03 | 1.19E-03 | 1.50E-04 | 1.88E-05 | 2.36E-06 | ≈3.0 | | 3 | 8.80E-04 | 5.74E-05 | 3.62E-06 | 2.27E-07 | 1.42E-08 | ≈4.0 | | 4 | 7.55E-05 | 2.39E-06 | 7.49E-08 | 2.34E-09 | 7.33E-11 | ≈5.0 | ### 条件数改善效果 重构前后质量矩阵条件数的比较显示： - 初始构造的条件数随多项式阶数指数增长 - Approach 2能够稳定地将条件数降至1.0 - Approach 1在高阶时出现数值不稳定 ### 预处理效果 对于系数矩阵$A$和$\tilde{A}$，Jacobi预处理器和行归一化预处理器都能显著改善条件数，特别是在高阶和小网格尺寸情况下。 ## 相关工作 ### IFE方法发展 - 传统IFE方法：只能近似满足界面条件 - 扩展IFE方法：基于Cauchy扩展的高阶方法 - 本文的GC-IFE：精确满足界面条件的几何一致性方法 ### 正交化技术 - 虚元方法中的正交多项式构造 - DG方法中的正交基函数 - 本文将正交化思想引入IFE框架 ## 结论与讨论 ### 主要结论 1. **理论贡献**：建立了基于Frenet装置的GC-IFE理论框架 2. **算法创新**：提出了两种构造方法和两种正交化策略 3. **数值验证**：证明了方法的最优收敛性和数值稳定性 4. **实用价值**：提供了完整的MATLAB实现 ### 局限性 1. **几何限制**：要求界面曲线足够光滑（$C^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篇相关文献，涵盖了微分几何、数值分析、有限元方法等多个领域的重要工作，为本研究提供了坚实的理论基础。