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 de Fonctions de Base pour la Méthode d'Éléments Finis Immergés Conformes à la Géométrie

Informations Fondamentales

  • ID de l'article: 2510.12018
  • Titre: Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method
  • Auteurs: Slimane Adjeril (Virginia Tech), Tao Lin (Virginia Tech), Haroun Meghaichi (Ohio State University)
  • Classification: math.NA cs.NA
  • Date de publication: 15 octobre 2025
  • Lien de l'article: https://arxiv.org/abs/2510.12018

Résumé

Cet article propose un nouveau cadre basé sur l'appareil de Frenet pour la construction de fonctions d'éléments finis immergés (IFE) géométriquement conformes d'ordre élevé résolvant les problèmes d'interface. L'article détaille le processus de construction des fonctions de base IFE locales, ainsi qu'une nouvelle méthode de construction de bases orthonormalisées utilisant la décomposition en valeurs singulières (SVD) de matrices de Vandermonde généralisées locales. Des implémentations d'exemples MATLAB sont fournies, démontrant la concision et l'extensibilité du cadre.

Contexte et Motivation de la Recherche

Contexte du Problème

Cette recherche aborde les problèmes d'interface elliptique: −∇⋅(β∇u)=f, sur Ω−∪Ω+-\nabla \cdot (\beta\nabla u) = f, \text{ sur } \Omega^- \cup \Omega^+ où le domaine de calcul Ω⊂R2\Omega \subset \mathbb{R}^2 est divisé par l'interface Γ\Gamma en deux sous-domaines Ω−\Omega^- et Ω+\Omega^+, et le coefficient de diffusion β\beta est une fonction constante par morceaux.

Motivation de la Recherche

  1. Limitations des méthodes IFE existantes: Les méthodes traditionnelles d'éléments finis immergés ne peuvent que satisfaire approximativement les conditions de saut à l'interface
  2. Exigences de précision d'ordre élevé: Nécessité de construire des fonctions IFE pour des polynômes d'ordre arbitraire
  3. Exigences de conformité géométrique: Nécessité de fonctions IFE satisfaisant exactement les conditions de saut à l'interface
  4. Stabilité numérique: Nécessité de systèmes de fonctions de base avec des nombres de condition favorables

Importance

Les problèmes d'interface apparaissent largement dans les écoulements multiphasiques, la science des matériaux, la biomédecine et autres domaines. Les méthodes numériques de résolution haute précision sont d'une importance capitale pour les applications d'ingénierie.

Contributions Principales

  1. Proposition d'un cadre de construction GC-IFE basé sur l'appareil de Frenet: Utilisation du repère de Frenet en géométrie différentielle pour simplifier les conditions de saut à l'interface
  2. Développement de deux méthodes de construction de fonctions de base:
    • Méthode de construction initiale basée sur des formats spéciaux
    • Méthode de construction généralisée basée sur l'idée d'extension
  3. Proposition de deux algorithmes de réorthogonalisation:
    • Approche 1: Méthode basée sur la SVD de la matrice de masse
    • Approche 2: Méthode basée sur la SVD de la matrice de Vandermonde généralisée
  4. Fourniture d'une implémentation MATLAB complète: Démonstration de la praticité et de l'extensibilité de la méthode

Détails de la Méthode

Transformation de Frenet

L'idée centrale consiste à utiliser la transformation de Frenet pour convertir les coordonnées physiques (x,y)(x,y) en coordonnées 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)

où g(ξ)g(\xi) est la paramétrisation de la courbe d'interface et n(ξ)n(\xi) est le vecteur normal.

Transformation des Conditions de Saut à l'Interface

Dans le système de coordonnées locales, les conditions de saut à l'interface deviennent:

  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éthodes de Construction de Fonctions de Base

Méthode 1: Construction par Format Spécial

Construction de fonctions de base de la forme: ϕ^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)

Pour 1≤j≤m,0≤i≤m1 \leq j \leq m, 0 \leq i \leq m, ces fonctions satisfont automatiquement les conditions de saut.

Méthode 2: Construction Généralisée par Extension

Adoption d'une forme plus générale:

\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}$$ Les coefficients de la matrice $C^{\pm}$ sont déterminés par la résolution d'un système d'équations linéaires. ### Réorthogonalisation #### Approche 1: SVD de la Matrice de Masse Soit $M_q = V_1\Lambda V_1^T$, on choisit $Q_1 = V_1\Lambda^{-1/2}$ de sorte que la matrice de masse reconstruite soit la matrice identité. #### Approche 2: SVD de la Matrice de Vandermonde Décomposition en valeurs singulières de $\tilde{V} = \sqrt{W}V$: $\tilde{V} = U_2\Sigma V_2^T$, on choisit $Q_2 = V_2\Sigma^{-1}$. L'analyse théorique montre que l'Approche 2 est plus stable numériquement, car elle évite la décomposition en valeurs singulières d'une matrice dont le nombre de condition est au carré. ## Configuration Expérimentale ### Problèmes de Test Considération du domaine $\Omega = [-1,1]^2$, interface $x^2 + y^2 = r_0^2$ ($r_0 = 1/\sqrt{3}$), fonction de test: $$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}$$ ### Indicateurs d'Évaluation - Ordre de convergence de l'erreur de projection $L^2$ - Nombre de condition de la matrice de masse - Stabilité numérique de la construction des fonctions de base ## Résultats Expérimentaux ### Vérification de la Convergence Pour différents ordres polynomiaux $m = 1,2,3,4$, l'erreur de projection $L^2$ présente un ordre de convergence optimal: | $m$ | $N = 16$ | $N = 32$ | $N = 64$ | $N = 128$ | $N = 256$ | Ordre de convergence | |-----|----------|----------|----------|-----------|-----------|----------------------| | 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 | ### Effet d'Amélioration du Nombre de Condition La comparaison des nombres de condition de la matrice de masse avant et après reconstruction montre: - Le nombre de condition de la construction initiale augmente exponentiellement avec l'ordre polynomial - L'Approche 2 peut réduire le nombre de condition de manière stable à 1.0 - L'Approche 1 présente une instabilité numérique aux ordres élevés ### Effet du Préconditionnement Pour les matrices de coefficients $A$ et $\tilde{A}$, les préconditionneurs de Jacobi et de normalisation par ligne améliorent tous deux significativement le nombre de condition, particulièrement dans les cas d'ordre élevé et de petite taille de maille. ## Travaux Connexes ### Développement des Méthodes IFE - Méthodes IFE traditionnelles: Satisfaction approximative des conditions d'interface - Méthodes IFE étendues: Méthodes d'ordre élevé basées sur l'extension de Cauchy - GC-IFE de cet article: Méthode de conformité géométrique satisfaisant exactement les conditions d'interface ### Techniques d'Orthogonalisation - Construction de polynômes orthogonaux dans les méthodes d'éléments virtuels - Fonctions de base orthogonales dans les méthodes DG - Introduction de l'idée d'orthogonalisation dans le cadre IFE dans cet article ## Conclusions et Discussion ### Conclusions Principales 1. **Contribution théorique**: Établissement d'un cadre théorique GC-IFE basé sur l'appareil de Frenet 2. **Innovation algorithmique**: Proposition de deux méthodes de construction et deux stratégies d'orthogonalisation 3. **Vérification numérique**: Démonstration de la convergence optimale et de la stabilité numérique de la méthode 4. **Valeur pratique**: Fourniture d'une implémentation MATLAB complète ### Limitations 1. **Restrictions géométriques**: Exigence que la courbe d'interface soit suffisamment lisse ($C^3$) 2. **Complexité de calcul**: Coût d'assemblage de matrice élevé dans les cas d'ordre élevé 3. **Extensibilité**: Actuellement applicable uniquement aux problèmes d'interface elliptique bidimensionnels ### Directions Futures 1. Extension aux problèmes tridimensionnels et géométries complexes 2. Application aux problèmes d'interface dépendant du temps 3. Extension aux systèmes de mécanique élastique, Stokes, etc. 4. Développement d'algorithmes d'implémentation plus efficaces ## Évaluation Approfondie ### Avantages 1. **Rigueur théorique**: L'appareil de Frenet basé sur la géométrie différentielle fournit une base théorique solide 2. **Innovation méthodologique**: Première réalisation d'une méthode IFE d'ordre élevé satisfaisant exactement les conditions de saut à l'interface 3. **Stabilité numérique**: La réorthogonalisation améliore significativement la stabilité numérique 4. **Implémentation complète**: Fourniture de descriptions d'algorithmes détaillées et de code MATLAB 5. **Vérification suffisante**: L'analyse théorique et les expériences numériques se corroborent mutuellement ### Insuffisances 1. **Portée d'application limitée**: Actuellement limitée aux problèmes d'interface elliptique bidimensionnels 2. **Coût de calcul**: L'analyse de la complexité de calcul dans les cas d'ordre élevé est insuffisante 3. **Exigences géométriques strictes**: Exigences élevées concernant la régularité de la courbe d'interface 4. **Cas d'application pratique**: Manque de vérification sur des problèmes d'ingénierie complexes ### Impact 1. **Valeur académique**: Fourniture d'un nouveau cadre théorique et d'une voie technique pour les méthodes IFE 2. **Valeur pratique**: Fourniture d'outils efficaces pour la résolution haute précision de problèmes d'interface 3. **Reproductibilité**: L'implémentation complète du code facilite l'utilisation et l'extension par d'autres chercheurs 4. **Potentiel de développement**: Établissement d'une base solide pour les recherches ultérieures ### Scénarios d'Application 1. **Écoulements multiphasiques**: Simulation numérique d'interfaces entre différents fluides 2. **Science des matériaux**: Calcul haute précision de problèmes d'interface de matériaux composites 3. **Biomédecine**: Problèmes de transfert thermique et de masse aux interfaces tissulaires 4. **Géophysique**: Problèmes de propagation d'ondes aux interfaces géologiques ## Références L'article cite 22 références pertinentes, couvrant plusieurs domaines importants tels que la géométrie différentielle, l'analyse numérique et les méthodes d'éléments finis, fournissant une base théorique solide pour cette recherche.