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

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

基本信息

  • 论文ID: 2212.11934
  • 标题: A localized reduced basis approach for unfitted domain methods on parameterized geometries
  • 作者: Margarita Chasapi, Pablo Antolin, Annalisa Buffa (École Polytechnique Fédérale de Lausanne & CNR)
  • 分类: math.NA cs.NA
  • 发表期刊: Computer Methods in Applied Mechanics and Engineering (已接收)
  • 论文链接: https://arxiv.org/abs/2212.11934

摘要

本文提出了一个降阶建模(ROM)框架，用于求解定义在未拟合几何(unfitted geometries)上的参数化二阶线性椭圆偏微分方程。该方法旨在构建高效的基于投影的ROM，依赖于降基方法(reduced basis method)和离散经验插值等技术。由于几何参数在未拟合域离散化中的存在给标准ROM的应用带来了挑战，本文提出了基于(i)快照在背景网格上的扩展和(ii)局部化策略来减少降基函数数量的方法论。所得方法在计算上高效且准确，同时对底层离散化选择不可知。

研究背景与动机

核心问题

  1. 未拟合域方法的挑战: 在计算机辅助设计(CAD)和等几何分析(IGA)中，复杂几何通常通过修剪(trimming)操作处理，导致未拟合网格，其中几何与解的离散化解耦。
  2. 参数化几何的降阶建模困难: 当几何参数变化时，活跃基函数集合会发生变化，导致快照解向量可能具有不同长度，这阻碍了标准降基方法的直接应用。
  3. 非仿射参数依赖: 几何参数导致微分算子对参数的非仿射依赖，破坏了降基方法的基本假设。

研究重要性

  • 工业相关性: 修剪是CAD中表示复杂形状的主流技术，对实现统一的设计到分析框架至关重要
  • 多查询场景: 在设计优化、不确定性量化、控制等应用中需要多次求解参数化PDE
  • 计算效率需求: 需要高效的离线/在线程序来实现计算加速

现有方法局限性

  • 标准ROM需要参考域映射，高度依赖于具体问题
  • 现有未拟合域ROM工作未解决非仿射依赖问题，缺乏超降维技术
  • 全局降基对于扩展解的非线性流形近似效果差，需要大量基函数

核心贡献

  1. 提出了针对参数化未拟合几何的完整降维框架，结合超降维技术实现高效的离线/在线分解
  2. 开发了基于参数的局部化策略，通过k-means聚类构建多个小维度的局部子空间，显著减少降基函数数量
  3. 设计了与离散化方法无关的通用框架，适用于有限元和等几何方法
  4. 实现了基于径向基函数(RBF)插值的非侵入式在线系数计算，避免了对高保真度汇编例程的侵入性修改
  5. 在数值实验中证明了显著的计算加速：相比标准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}(‖μ-μⱼ‖₂)

3. 局部化策略

  • 参数空间聚类: 使用k-means算法将参数空间划分为N_c个子区域
  • 局部基构建: 为每个聚类分别构建DEIM近似和降基
  • 在线切换: 根据参数到聚类中心的距离选择相应的局部基

技术创新点

  1. 几何驱动的聚类策略: 基于边界距离的参数聚类，确保相似几何配置分组
  2. 多层次局部化: 同时对DEIM近似和降基进行局部化
  3. 非侵入式实现: 通过RBF插值避免对底层汇编例程的修改
  4. 通用框架设计: 与具体离散化方法解耦，支持多种未拟合域技术

实验设置

测试问题

  1. Poisson方程:
    • 1D参数化：带圆孔的正方形域，孔心沿对角线移动
    • 2D参数化：孔心位置和半径同时变化
  2. 线性弹性问题:
    • 多孔四分之一圆柱：多个圆孔半径参数化
    • 3D立方体球形包含物：球半径参数化

离散化设置

  • 基函数: 二次和三次C¹-C²连续B样条
  • 网格: 32×32 Cartesian网格(2D)，8×8×8网格(3D)
  • 自由度: 1225-3000个DOF

评价指标

  • 精度: L²和H¹范数相对误差
  • 效率: 在线CPU时间和加速比
  • 降维效果: 奇异值衰减和基函数数量

参数设置

  • POD容差: ε_ = 10⁻⁵, ε_^d = 10⁻⁷
  • 训练样本: N_s = 100-1000, N_s^d = 250-1000
  • 聚类数: N_c = 4-16, N_c^d = 8-16
  • RBF类型: 三次径向基函数

实验结果

主要结果

1. DEIM近似效果

  • 全局方法: 需要349个矩阵项和124个向量项(1D情况)
  • 局部方法: 最大31-38个矩阵项和11-15个向量项(16个聚类)
  • 精度: 达到10⁻⁵级别误差

2. 降基维数对比

问题全局ROM局部ROM聚类数
Poisson 1DN=182N=354
Poisson 2DN=201N=1716
弹性多孔N=198N=1610

3. 计算加速比

问题全局ROM加速局部ROM加速
Poisson 2D9.1×17.6×
弹性多孔16.5×130×

消融实验

  • 聚类数影响: 增加聚类数显著减少基函数数量，提高精度
  • 训练样本大小: 需要足够丰富的训练集(N_s^d ≥ 100)保证DEIM精度
  • RBF插值精度: 三次RBF提供良好的光滑性而无虚假振荡

3D验证

在3D立方体球形包含物问题中：

  • DEIM基函数从180个减少到35个
  • 降基维数从35个减少到14个
  • 验证了方法在高维几何中的有效性

相关工作

等几何降阶建模

  • 已有工作主要集中在保形几何的IGA-ROM结合
  • 在流体力学、抛物型问题、心脏电生理等领域的成功应用
  • 缺乏对未拟合域的系统研究

未拟合域降阶建模

  • 早期工作: 虚拟域方法+PGD、低秩近似用于时变界面
  • 近期进展: CutFEM-ROM、Shifted Boundary Method-ROM
  • 局限性: 缺乏超降维技术，未充分解决非仿射依赖问题

局部降基方法

  • 聚类方法: 基于解快照聚类的局部ROM
  • hp-降基: 参数空间自适应划分
  • 本文贡献: 首次针对未拟合几何的参数驱动局部化

结论与讨论

主要结论

  1. 有效性验证: 所提方法在保持精度的同时显著降低了在线计算成本
  2. 通用性: 框架与底层离散化方法无关，适用于多种未拟合域技术
  3. 扩展性: 成功扩展到3D复杂几何，展现了良好的可扩展性

局限性

  1. 聚类数选择: 需要先验选择聚类数，虽然提供了方差准则但仍需经验判断
  2. 参数维数: 高维参数空间可能导致RBF插值条件数问题
  3. 几何复杂性: 对于极其复杂的几何变化，可能需要更多聚类

未来方向

  1. Greedy算法: 应用Greedy算法构建局部降基和误差认证
  2. 四阶PDE: 扩展到更复杂的PDE问题
  3. 优化应用: 在设计、形状和拓扑优化中的应用

深度评价

优点

  1. 方法创新性强: 首次系统地将局部化策略应用于未拟合域ROM，解决了重要的技术难题
  2. 理论框架完整: 从快照扩展到局部化策略，形成了完整的方法论体系
  3. 实验验证充分: 涵盖2D/3D、不同PDE类型、多种参数化场景的全面验证
  4. 实用价值高: 显著的计算加速和与离散化方法无关的特性使其具有广泛应用前景

不足

  1. 理论分析有限: 缺乏收敛性和误差估计的严格理论分析
  2. 参数依赖性: 方法效果高度依赖于聚类策略和参数选择
  3. 复杂性权衡: 虽然在线效率高，但离线阶段的计算和存储成本显著增加

影响力

  1. 学术贡献: 为未拟合域ROM开辟了新的研究方向，预期引发后续研究
  2. 工程价值: 对CAD-分析一体化和工业几何优化具有重要意义
  3. 可复现性: 提供了详细的算法描述和参数设置，便于复现和应用

适用场景

  • 涉及复杂几何变化的参数化设计优化
  • 修剪几何的多查询分析
  • 实时仿真和控制应用
  • 不确定性量化中的几何参数研究

参考文献

论文引用了丰富的相关文献，包括：

  • 未拟合域方法的经典工作(Peskin, Hughes等)
  • 等几何分析的基础文献
  • 降基方法和经验插值的理论基础
  • 局部ROM的前沿研究

该论文为未拟合域上的参数化PDE降阶建模提供了重要的方法论贡献，在理论创新和实际应用之间取得了良好的平衡。