2025-11-18T02:58:13.538030

Reduced order modelling of nonaffine problems on parameterized NURBS multipatch geometries

Chasapi, Antolin, Buffa

This contribution explores the combined capabilities of reduced basis methods and IsoGeometric Analysis (IGA) in the context of parameterized partial differential equations. The introduction of IGA enables a unified simulation framework based on a single geometry representation for both design and analysis. The coupling of reduced basis methods with IGA has been motivated in particular by their combined capabilities for geometric design and solution of parameterized geometries. In most IGA applications, the geometry is modelled by multiple patches with different physical or geometrical parameters. In particular, we are interested in nonaffine problems characterized by a high-dimensional parameter space. We consider the Empirical Interpolation Method (EIM) to recover an affine parametric dependence and combine domain decomposition to reduce the dimensionality. We couple spline patches in a parameterized setting, where multiple evaluations are performed for a given set of geometrical parameters, and employ the Static Condensation Reduced Basis Element (SCRBE) method. At the common interface between adjacent patches a static condensation procedure is employed, whereas in the interior a reduced basis approximation enables an efficient offline/online decomposition. The full order model over which we setup the RB formulation is based on NURBS approximation, whereas the reduced basis construction relies on techniques such as the Greedy algorithm or proper orthogonal decomposition (POD). We demonstrate the developed procedure using an illustrative model problem on a three-dimensional geometry featuring a multi-dimensional geometrical parameterization.

academic

非アフィン問題の降次元化モデリング：パラメータ化NURBS多重パッチ幾何への応用

基本情報

論文ID: 2211.07348
タイトル: Reduced order modelling of nonaffine problems on parameterized NURBS multipatch geometries
著者: Margarita Chasapi, Pablo Antolin, Annalisa Buffa
分類: math.NA cs.NA
発表日: 2022年11月14日 (arXiv投稿)
論文リンク: https://arxiv.org/abs/2211.07348

要旨

本論文は、降次元基底法(Reduced Basis Methods)と等幾何解析(IsoGeometric Analysis, IGA)をパラメータ化偏微分方程式の求解に組み合わせた応用を探索している。等幾何解析の導入により、単一の幾何表現に基づく設計と分析の統一的なシミュレーション枠組みが実現される。特に複数のパッチでモデル化された複雑な幾何体に対して、本論文は高次元パラメータ空間の特性を持つ非アフィン問題に焦点を当てている。経験的補間法(EIM)を採用してアフィンパラメータ依存性を回復し、領域分解と組み合わせて次元削減を実現している。パラメータ化設定においてスプライン補丁を結合し、静的凝聚降次元基底要素(SCRBE)法を採用している。隣接する補丁の共通界面では静的凝聚手順を適用し、内部では降次元基底近似により効率的なオフライン/オンライン分解を実現している。

研究背景と動機

問題定義

中核問題: 複雑なパラメータ化幾何体上で偏微分方程式を求解する際、従来の全次元モデルは計算コストが過大であり、特に多重クエリと実時間応用シナリオにおいて問題となる。
技術的課題:
- 幾何パラメータ化に起因する非アフィンパラメータ依存性が降次元基底法の効率性を損なう
- 高次元パラメータ空間により標準的な降次元基底法の性能が低下する
- 多重パッチ幾何の複雑性がモデリング難度を増加させる
応用要件: 産業応用における複雑な幾何体は通常、多重パッチ表現を必要とし、多数の幾何パラメータを含むため、効率的なパラメータ化求解法が必要である。

研究動機

統一的枠組み: IGAを利用して幾何設計と数値解析の統一的表現を実現する
計算効率: 降次元基底法を通じてパラメータ化問題の高速求解を実現する
産業関連性: 産業的意義を持つ複雑な多重パッチ幾何問題に対応する

核心的貢献

方法の融合: IGA と降次元基底法を非アフィン多重パッチ幾何問題に対して初めて体系的に結合した
EIM-IGA統合: スプライン近似枠組みの下で経験的補間法を再定式化し、精密な幾何写像を利用してアフィン分解を回復した
SCRBE拡張: 静的凝聚降次元基底要素法をNURBS多重パッチ幾何に拡張した
領域分解戦略: パッチベースの領域分解法を提案し、高次元パラメータ空間に対応した
数値検証: 三次元多重パッチ幾何上で方法の有効性を検証した

方法の詳細

問題設定

パラメータ化楕円型偏微分方程式を求解する： $-\Delta u = f \text{ in } \Omega(\mu)$ $u = 0 \text{ on } \Gamma_D(\mu)$ $\frac{\partial u}{\partial n} = 0 \text{ on } \Gamma_N(\mu)$

ここで $\Omega(\mu) \subset \mathbb{R}^d$ はパラメータ化領域、 $\mu \in \mathcal{P} \subset \mathbb{R}^P$ は幾何パラメータである。

モデルアーキテクチャ

1. NURBS多重パッチ幾何表現

パラメータ領域: $\hat{\Omega} = [0,1]^{\hat{d}}$
幾何写像: $F(\xi;\mu) = \sum_i B_{i,p}(\xi)P_i(\mu)$
多重パッチ分解: $\Omega(\mu) = \bigcup_{k=1}^{n_p} \Omega_k(\mu^{(k)})$

2. 経験的補間法(EIM)

非アフィンパラメータ依存性を参照領域上の補間問題に変換する：

変換された弱形式: $\alpha(u_h, v_h;\mu) = \int_{\hat{\Omega}} (\nabla \hat{u}_h DF^{-1}) \cdot (\nabla \hat{v}_h DF^{-1})|\det DF| d\hat{\Omega}$

EIM近似: $g_M^\alpha(\xi;\mu) = I_{\xi}^M g^\alpha(\xi;\mu) = \sum_{m=1}^{M^\alpha} \theta_m^\alpha(\mu)\phi_m^\alpha(\xi)$

3. SCRBE法

静的凝聚分解:

バブル空間: $V_{h,b}^k = \{v_h \in V_h^k : v_h|_{\gamma_{k,j}} = 0\}$
ポート空間: $V_{h,\gamma}^{k,j} = \{v_h \in V_h^k : v_h|_{\gamma_{k,j}} \neq 0\}$

解の表現: $u_h(\mu)|_{\Omega_k} = b_{h,f}^k(\mu^{(k)}) + \sum_{j=1}^{n_\gamma^k} \sum_{r=1}^{N_{\gamma}^{k,j}} \hat{u}_{G_k(j),r}(\mu) \phi_{h}^{k,j,r}(\mu^{(k)})$