[논문 리뷰] $h$ and $hp$-adaptive Interpolation by Transformed Snapshots for Parametric and Stochastic Hyperbolic PDEs

Gerrit Welper|arXiv (Cornell University)|2017. 10. 31.

Probabilistic and Robust Engineering Design참고 문헌 35인용 수 15

한 줄 요약

이 논문은 충격 구조의 위상 변화가 발생하는 파rametric 및 확률적 초음속 PDE를 해결하기 위해 변환된 스냅샷을 이용한 $h$ 및 $hp$-적응형 보간법을 제안한다. 국소적인 매개변수 공간의 메esh 조정과 좌표축에 따라 정렬된 변환을 활용함으로써, 비연속성에도 불구하고 지수 수렴 속도를 달성하며, 충격 구조가 변화하는 문제에서 기존의 축소 기저 및 POD 방법보다 뛰어난 성능을 보인다.

ABSTRACT

The numerical approximation of solutions of parametric or stochastic hyperbolic PDEs is still a serious challenge. Because of shock singularities, most methods from the elliptic and parabolic regime, such as reduced basis methods, POD or polynomial chaos expansions, show a poor performance. Recently, Welper [Interpolation of functions with parameter dependent jumps by transformed snapshots. SIAM Journal on Scientific Computing, 39(4):A1225-A1250, 2017] introduced a new approximation method, based on the alignment of the jump sets of the snapshots. If the structure of the jump sets changes with parameter, this assumption is too restrictive. However, these changes are typically local in parameter space, so that in this paper, we explore $h$ and $hp$-adaptive methods to resolve them. Since local refinements do not scale to high dimensions, we introduce an alternative "tensorized" adaption method.

연구 동기 및 목표

• 충격 비연속성으로 인해 전통적인 축소 기저, POD 및 다항식 혼합 방법의 성능이 열 劣하는 문제를 해결하기 위해.
• 변환된 스냅샷 보간법(TSI)의 한계를 극복하기 위해, 매개변수 간에 절단 집합의 위상이 고정되어 있음을 가정하는 것을 넘어서기 위해.
• 충격 구조의 국소적 위상 변화를 해소하기 위해 매개변수 공간에서 $h$ 및 $hp$-적응형 메쉬 조정 전략을 개발하기 위해.
• 고차원 매개변수 공간에서의 확장성 문제를 완화하기 위해 텐서형 적응 전략을 도입하기 위해.
• 변형된 충격 위상 변화를 보이는 모델 문제에 대해 지수 수렴 속도를 달성하고, 벤치마크 PDE에 대한 수치 실험을 통해 검증하기 위해.

제안 방법

• 매개변수에 따라 변하는 영역 변환 $\phi(\mu,\eta)(x)$를 사용하여, 다양한 매개변수 값 간에 절단 비연속성을 정렬하는 변환된 스냅샷 보간법(TSI)을 적용한다.
• 충격 충돌이나 충격에서 희석파로의 전이와 같은 충격 위상 변화가 발생하는 국소 영역을 해소하기 위해 매개변수 공간에서 $h$- 및 $hp$-적응형 메쉬 조정을 수행한다.
• 국소적 메쉬 조정을 좌표축에 따라 정렬함으로써 텐서형 적응 전략을 구성하여, 고차원 매개변수 공간에서도 확장성을 확보한다.
• 변환된 스냅샷을 사용하여 매개변수 공간에서 라그랑주 보간을 수행한다: $ u(x,\mu) \approx \sum_{\eta \in \mathcal{P}_n} \ell_\eta(\mu) \, u(\phi(\mu,\eta)(x), \eta) $.
• 차수 2의 유한요소 기반 변환 $\phi(\mu,\eta)(x) = x + \varphi(\mu,\eta)(x)$를 사용하며, 경계 조건 $\varphi \cdot \nu = 0$를 적용하여 미분동형성과 영역 보존을 보장한다.
• 공간 해상도 향상을 위해 두 번째 차수의 정적 적분 및 다항식 근사법을 변환 및 해 평가에 활용한다.

실험 결과

연구 질문

• RQ1 $h$ 및 $hp$-적응형 메쉬 조정 전략이 파rametric 초음속 PDE의 충격 구조에서 국소적 위상 변화를 효과적으로 해소할 수 있는가?
• RQ2 좌표축에 따라 정렬된 변환을 기반으로 한 텐서형 적응 전략이 표준 국소적 메쉬 조정 대비 고차원 매개변수 공간에서의 확장성에 얼마나 기여하는가?
• RQ3 변형된 충격 위상 변화를 보이는 문제에 대해 제안된 적응형 TSI 방법이 달성할 수 있는 수렴 속도는 무엇인가?
• RQ4 비연속 해를 갖는 문제에서 기존의 축소 기저 및 POD 방법 대비 오차 및 강인성 측면에서 이 방법은 어떻게 성능을 발휘하는가?
• RQ5 이 방법은 충격 충돌 및 충격에서 희석파로의 전이와 같은 중요한 전이 현상을 최소한의 메쉬 조정으로도 잘 포착할 수 있는가?

주요 결과

• 모델 문제에 대해 충격 위상 변화가 발생하는 경우 지수 수렴 속도를 달성하며, 기존의 축소 기저 및 POD 방법보다 뛰어난 성능을 보였다.
• 초기 조건에 두 개의 절단이 포함된 비점성 버거스 방정식에 대한 수치 실험 결과, 국소적 메쉬 조정이 충격 충돌 및 희석파 전이 현상을 효과적으로 포착하였다.
• $\mu = -0.5$인 경우, 충격 충돌 시점에 메쉬 조정이 집중되었고, $\mu = 0.5$인 경우는 $t=0$에서 초기 절단이 희석파로 전이되는 시점에 메쉬 조정이 이루어졌다.
• 시험 케이스에서 전체 오차는 이미 공간 해상도 $h = 0.01$을 갖는 스냅샷의 오차 수준에 가까워졌으며, 이는 최소한의 메쉬 조정으로도 높은 정확도를 달성했음을 시사한다.
• 앞면 단계 흐름 예시에서는 고마하 수에서 새로운 충격이 형성되는 것을 포착하였고, 메쉬 조정은 전이점 근처에 국소화되었다.
• 이 방법은 해상에서 계단 효과를 피하고 오차를 공간 메쉬 크기의 약 두 배 내외로 유지함으로써 강인성과 정확성을 입증하였다 ($h = 0.025$).