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

Reduced order modelling of nonaffine problems on parameterized NURBS multipatch geometries

基本信息

论文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)})$