2025-11-12T12:07:10.593714

A localized reduced basis approach for unfitted domain methods on parameterized geometries

Chasapi, Antolin, Buffa

This work introduces a reduced order modeling (ROM) framework for the solution of parameterized second-order linear elliptic partial differential equations formulated on unfitted geometries. The goal is to construct efficient projection-based ROMs, which rely on techniques such as the reduced basis method and discrete empirical interpolation. The presence of geometrical parameters in unfitted domain discretizations entails challenges for the application of standard ROMs. Therefore, in this work we propose a methodology based on i) extension of snapshots on the background mesh and ii) localization strategies to decrease the number of reduced basis functions. The method we obtain is computationally efficient and accurate, while it is agnostic with respect to the underlying discretization choice. We test the applicability of the proposed framework with numerical experiments on two model problems, namely the Poisson and linear elasticity problems. In particular, we study several benchmarks formulated on two-dimensional, trimmed domains discretized with splines and we observe a significant reduction of the online computational cost compared to standard ROMs for the same level of accuracy. Moreover, we show the applicability of our methodology to a three-dimensional geometry of a linear elastic problem.

academic

非適合領域法における局所化縮約基底法：パラメータ化幾何への応用

基本情報

論文ID: 2212.11934
タイトル: A localized reduced basis approach for unfitted domain methods on parameterized geometries
著者: Margarita Chasapi, Pablo Antolin, Annalisa Buffa (スイス連邦工科大学ローザンヌ校 & CNR)
分類: math.NA cs.NA
掲載誌: Computer Methods in Applied Mechanics and Engineering (採択済み)
論文リンク: https://arxiv.org/abs/2212.11934

要旨

本論文は、非適合幾何上で定義されたパラメータ化二階線形楕円型偏微分方程式を求解するための縮約秩モデル化(ROM)フレームワークを提案している。本手法は、縮約基底法および離散経験的補間などの技術に依存する、効率的な投影ベースのROMの構築を目指している。幾何パラメータが非適合領域離散化に存在することが標準ROMの適用に課題をもたらすため、本論文は(i)背景メッシュ上のスナップショット拡張と(ii)縮約基底関数の数を削減するための局所化戦略に基づく方法論を提案する。得られた手法は計算効率が高く正確であり、同時に基礎となる離散化選択に対して不可知的である。

研究背景と動機

核心的課題

非適合領域法の課題: コンピュータ支援設計(CAD)および等幾何解析(IGA)において、複雑な幾何は通常トリミング操作によって処理され、非適合メッシュをもたらす。ここで幾何は解の離散化から切り離される。
パラメータ化幾何における縮約秩モデル化の困難: 幾何パラメータが変化する場合、活性基底関数の集合が変化し、スナップショット解ベクトルが異なる長さを持つ可能性があり、これが標準縮約基底法の直接的な適用を阻害する。
非アフィン型パラメータ依存性: 幾何パラメータは微分作用素のパラメータに対する非アフィン型依存性をもたらし、縮約基底法の基本的な仮定を破壊する。

研究の重要性

産業的関連性: トリミングはCADにおける複雑な形状表現の主流技術であり、統一的な設計から解析へのフレームワークの実現に不可欠である
多クエリシナリオ: 設計最適化、不確実性定量化、制御などの応用において、パラメータ化PDEの多重求解が必要である
計算効率の必要性: 計算加速を実現するための効率的なオフライン/オンライン手順が必要である

既存手法の限界

標準ROMは参照領域マッピングを必要とし、問題に高度に依存している
既存の非適合領域ROM研究は非アフィン型依存性を解決しておらず、超縮約技術が不足している
グローバル縮約基底は拡張解の非線形多様体近似に対して効果が低く、多くの基底関数を必要とする

核心的貢献

パラメータ化非適合幾何に対する完全な縮約化フレームワークを提案し、超縮約技術を組み込んで効率的なオフライン/オンライン分解を実現
パラメータベースの局所化戦略を開発し、k-means クラスタリングを通じて複数の小次元局所部分空間を構築し、縮約基底関数の数を大幅に削減
離散化手法に依存しない汎用フレームワークを設計し、有限要素法および等幾何法に適用可能
径向基底関数(RBF)補間に基づく非侵襲的なオンライン係数計算を実装し、高忠実度アセンブリルーチンへの侵襲的な修正を回避
数値実験において顕著な計算加速を実証：標準ROMと比較して、同じ精度レベルで17.6倍から130倍のオンライン計算加速を実現

方法論の詳細

問題設定

パラメータ化非適合領域Ω(μ)上で定義された線形楕円型PDEを求解する：

Lu = f in Ω(μ)

ここでμ ∈ P ⊂ R^M は幾何パラメータベクトルであり、領域はカット操作により定義される：

Ω(μ) = Ω̂₀(μ) \ ⋃ᵢ₌₁ᴷ Ω̂ᵢ(μ)

核心的方法論アーキテクチャ

1. 非適合領域離散化

背景領域設定: パラメータ化領域Ω(μ)をパラメータに依存しない背景領域Ω₀に埋め込む
離散空間構築: 背景領域上で離散関数空間V_{h,0}を定義し、活性空間V_h(μ)はΩ(μ)と交差するサポートを持つ基底関数により張られる
拡張戦略: 異なるパラメータ下の解ベクトルを背景領域に拡張し、非活性領域でゼロに設定

2. 離散経験的補間法(DEIM)

アフィン分解: 剛性行列と右辺ベクトルをアフィン形式に近似する：
```
Â(μ) ≈ ∑_{q=1}^{Q_a} θ_q^a(μ)Âq
f̂(μ) ≈ ∑_{q=1}^{Q_f} θ_q^f(μ)f̂q
```
RBF補間: 径向基底関数補間を使用してオンラインで高速にパラメータ依存係数を計算する：
```
θ_q^a(μ) ≈ ∑_{j=1}^{N_s^d} ω_{q,j}^a φ_{q,j}(‖μ-μⱼ‖₂)
```