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

Construcción de Funciones Base para el Método de Elementos Finitos Inmersos Conformes a la Geometría

Información Básica

  • ID del Artículo: 2510.12018
  • Título: Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method
  • Autores: Slimane Adjerid (Virginia Tech), Tao Lin (Virginia Tech), Haroun Meghaichi (Ohio State University)
  • Clasificación: math.NA cs.NA
  • Fecha de Publicación: 15 de octubre de 2025
  • Enlace del Artículo: https://arxiv.org/abs/2510.12018

Resumen

Este artículo propone un nuevo marco basado en el aparato de Frenet para construir funciones de elementos finitos inmersos (IFE) de alto orden conformes a la geometría que resuelven problemas de interfaz. El artículo detalla el proceso de construcción de funciones base IFE locales, así como un nuevo método para construir bases ortonormalizadas utilizando la descomposición en valores singulares (SVD) de matrices de Vandermonde generalizadas locales. Se proporcionan implementaciones de ejemplo en MATLAB que demuestran la simplicidad y escalabilidad del marco.

Antecedentes de Investigación y Motivación

Contexto del Problema

Esta investigación aborda problemas de interfaz elíptica: −∇⋅(β∇u)=f, en Ω−∪Ω+-\nabla \cdot (\beta\nabla u) = f, \text{ en } \Omega^- \cup \Omega^+ donde el dominio computacional Ω⊂R2\Omega \subset \mathbb{R}^2 se divide por la interfaz Γ\Gamma en dos subdominios Ω−\Omega^- y Ω+\Omega^+, siendo β\beta una función de coeficiente de difusión por partes constante.

Motivación de la Investigación

  1. Limitaciones de los métodos IFE existentes: Los métodos tradicionales de elementos finitos inmersos solo pueden satisfacer aproximadamente las condiciones de salto en la interfaz
  2. Requisitos de precisión de alto orden: Se necesita construir funciones IFE de polinomios de orden arbitrario
  3. Requisitos de conformidad geométrica: Se necesitan funciones IFE que satisfagan exactamente las condiciones de salto en la interfaz
  4. Estabilidad computacional: Se requiere un sistema de funciones base con números de condición favorables

Importancia

Los problemas de interfaz aparecen ampliamente en flujos multifásicos, ciencia de materiales, ingeniería biomédica y otros campos. Los métodos numéricos de solución de alta precisión tienen una importancia significativa para aplicaciones de ingeniería.

Contribuciones Principales

  1. Se propone un marco de construcción GC-IFE basado en el aparato de Frenet: Utiliza marcos de Frenet de geometría diferencial para simplificar las condiciones de salto en la interfaz
  2. Se desarrollan dos métodos de construcción de funciones base:
    • Método de construcción inicial basado en formato especial
    • Método de construcción generalizada basado en ideas de extensión
  3. Se proponen dos algoritmos de reconstrucción ortogonal:
    • Método basado en SVD de matriz de masa (Enfoque 1)
    • Método basado en SVD de matriz de Vandermonde generalizada (Enfoque 2)
  4. Se proporciona implementación completa en MATLAB: Demuestra la practicidad y escalabilidad del método

Explicación Detallada del Método

Transformación de Frenet

La idea central es utilizar la transformación de Frenet para convertir coordenadas físicas (x,y)(x,y) en coordenadas locales (η,ξ)(\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)

donde g(ξ)g(\xi) es la parametrización de la curva de interfaz y n(ξ)n(\xi) es el vector normal.

Transformación de Condiciones de Salto en la Interfaz

En el sistema de coordenadas local, las condiciones de salto en la interfaz se convierten en:

  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

Métodos de Construcción de Funciones Base

Método 1: Construcción de Formato Especial

Se construyen funciones base de la forma: ϕ^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)

Para 1≤j≤m,0≤i≤m1 \leq j \leq m, 0 \leq i \leq m, estas funciones satisfacen automáticamente las condiciones de salto.

Método 2: Construcción de Extensión Generalizada

Se adopta una forma más general:

undefined