[논문 리뷰] Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximaTion (GIFT)

Elena Atroshchenko, Gang Xu외 2명|arXiv (Cornell University)|2017. 01. 01.

Advanced Numerical Analysis Techniques참고 문헌 12인용 수 1

한 줄 요약

이 논문은 등각 분석에서 기하 구조 파arameterization과 해 근사 사이의 분리를 가능하게 하는 기하 독립적 필드 근사(Geometry Independent Field approximation, GIFT)를 소개한다. NURBS를 기하에, PHT-splines를 필드에 사용하는 별도의 스퍼린 공간을 통해, CAD 기하를 변경하지 않고 국소 적응형 메esh 개선을 가능하게 하며, 비중첩 공간임에도 불구하고 최적 수렴 속도를 달성한다. 표준 패치 테스트의 실패에도 불구하고 이는 약한 일致성과 적절한 오차 지표 덕분에 일관성과 수렴성을 유지한다.

ABSTRACT

This paper presents an approach to generalize the concept of isogeometric analysis (IGA) by allowing different spaces for parameterization of the computational domain and for approximation of the solution field. The method inherits the main advantage of isogeometric analysis, i.e. preserves the original, exact CAD geometry (for example, given by NURBS), but allows pairing it with an approximation space which is more suitable/flexible for analysis, for example, T-splines, LR-splines, (truncated) hierarchical B-splines, and PHT-splines. This generalization offers the advantage of adaptive local refinement without the need to re-parameterize the domain, and therefore without weakening the link with the CAD model. We demonstrate the use of the method with different choices of the geometry and field splines, and show that, despite the failure of the standard patch test, the optimum convergence rate is achieved for non-nested spaces.

연구 동기 및 목표

• . 기존 등각 분석(IGA)의 강성 문제를 해결하고자 하며, 기하와 해 근사에 동일한 스퍼린 공간을 사용하는 전통적 접근의 제약을 극복하고자 한다.
• . NURBS 기반 IGA가 텐서 곱의 구조로 인해 국소 메쉬 개선이 어려운 점을 해결하고자 한다.
• . 정확한 CAD 기하를 유지하면서도 해 근사 공간을 국소적으로 풍부화시킬 수 있는 프레임워크를 개발하고자 한다.
• . 기하 매개변수화와 독립적인 해 근사 필드의 적응형 h-메쉬 개선이 가능하도록 하되, CAD 통합은 유지하고자 한다.
• . 해 근사 공간과 기하 공간이 비중첩일 경우에도 최적 수렴 속도를 달성할 수 있음을 검증하고자 한다.

제안 방법

• . GIFT는 기하 구조의 매개변수화(NURBS 또는 기타 기하 스퍼린)와 해 근사 공간(T-splines, LR-splines, PHT-splines, 또는 잘라낸 계층적 B-스플라인)을 분리한다.
• . 동일한 매개변수 영역에서 PDE를 풀기 위해 약한 일致성 있는 갈레르킨 형식을 사용하며, 해 변수의 근사 공간은 다르게 설정한다.
• . 국소 메쉬 개선은 해 근사 공간(예: PHT-splines)에만 수행되며, 기하 매개변수화를 변경하거나 재매개변수화가 필요하지 않다.
• . 각 매개변수 셀 단위로 PDE 잔여항을 이용해 잔여 오차 지표를 계산하고, 오차 임계값 기반의 마킹 전략을 통해 개선 대상 셀을 선정한다.
• . 개선은 PHT-spline 규칙에 따라 수행되며, 교차점과 T-정점에 새로운 기저 함수를 생성하고, 새로운 정점에서의 Bézier 계수를 0으로 재설정함으로써 기존 기저 함수를 갱신한다.
• . 새로운 해 근사 제어 변수는 이전 제어점에서 갱신하는 것이 아니라, 개선된 근사 공간에서 PDE를 재해결함으로써 확보된다.

실험 결과

연구 질문

• RQ1 . 등각 분석에서 기하와 필드 근사 간의 분리된 접근이 비중첩 근사 공간임에도 불구하고 최적 수렴 속도를 유지할 수 있는가?
• RQ2 . GIFT에서 표준 패치 테스트의 실패가 방법의 일관성과 수렴 성질을 손상시키는가?
• RQ3 . 원래 CAD 기하를 변경하거나 재매개변수화 없이도 해 근사 필드에 대해 국소 적응형 메쉬 개선을 수행할 수 있는가?
• RQ4 . 기하에 NURBS, 필드에 PHT-splines 등의 다른 스퍼린 유형을 사용할 경우 수치 해의 정확도와 효율성에 어떤 영향을 미치는가?
• RQ5 . 효과적이고 안정적인 국소 메쉬 개선을 보장하기 위해 어떤 오차 지표와 마킹 전략이 필요한가?

주요 결과

• . GIFT는 해 근사 공간과 기하 근사 공간이 비중첩일지라도 2D 및 3D 문제 모두에서 최적 수렴 속도를 달성한다.
• . 방법은 정확한 CAD 기하를 그대로 유지하여 등각 분석의 핵심 이점을 그대로 보존한다.
• . 표준 패치 테스트의 실패에도 불구하고, 약한 형식과 적절한 오차 지표 덕분에 방법은 일관성 있고 수렴성을 유지한다.
• . PHT-splines의 사용은 기하 매개변수화를 변경하지 않고도 해 근사 필드의 효과적인 국소 h-메쉬 개선을 가능하게 한다.
• . 잔여 오차 기반 지표 eK = ||f + Δuh||_L2(K) × hK 는 메쉬 개선 대상 셀을 마킹하는 데 신뢰할 수 있는 척도를 제공한다.
• . 하이브리드 마킹 전략(먼저 오차 기준 상위 ϵ% 셀을 마킹하고, 그 다음 평균값 기반 마킹을 적용)은 오차 분포가 매우 비균일한 경우의 안정성과 강건성을 향상시킨다.