[논문 리뷰] Improved approximation of eigenvalues in isogeometric methods for multi-patch geometries and Neumann boundaries

Thomas Horger, Alessandro Reali외 2명|arXiv (Cornell University)|2017. 01. 23.

Advanced Numerical Analysis Techniques인용 수 2

한 줄 요약

이 논문은 다패치 기하 구조에서 경계를 넘는 약한 연속성과 법선 도수의 페널티를 부과하는 고차수 페널티 기법을 제안한다. 모르타르 결합과 법선 도수 페널티를 조합함으로써, 이 방법은 제4차 문제와 고유값 문제에 대해 안정적이고 고차수의 이산화를 가능하게 하며, 제2차 고유값 문제에서의 스펙트럼 오염('이상치')을 효과적으로 제거하면서도 복잡한 기하 구조에서도 정확성을 유지한다.

ABSTRACT

We present a systematic study on higher-order penalty techniques for isogeometric mortar methods. In addition to the weak-continuity enforced by a mortar method, normal derivatives across the interface are penalized. The considered applications are fourth order problems as well as eigenvalue problems for second and fourth order equations. The hybrid coupling enables the discretization of fourth order problems in a multi-patch setting as well as a convenient implementation of natural boundary conditions. For second order eigenvalue problems, the pollution of the discrete spectrum - typically referred to as 'outliers' - can be avoided. Numerical results illustrate the good behaviour of the proposed method in simple systematic studies as well as more complex multi-patch mapped geometries for linear elasticity and Kirchhoff plates.

연구 동기 및 목표

• 이소지오메트릭 분석에서 제2차 문제의 이산 고유값 근사에서 발생하는 스펙트럼 오염('이상치')을 해결하기 위해.
• 복잡한 경계를 가진 다패치 기하 구조에서 제4차 문제에 대해 안정적이고 정확한 이산화를 가능하게 하기 위해.
• 연속성과 법선 도수 제약 조건을 약하게 부과하는 고차수 페널티 기법의 체계적인 프레임워크를 제공하기 위해.
• 이소지오메트릭 방법에서 자연 경계 조건(네이만 경계 조건)의 구현을 용이하게 하기 위해.

제안 방법

• 모르타르 방법을 통해 Nitsche의 방법과 페널티 항을 사용하여 패치 경계를 넘는 약한 연속성을 부과한다.
• 경계를 넘는 법선 도수를 고차수 페널티 항을 통해 페널티화하여 법선 허브의 연속성을 향상시킨다.
• 하이브리드 결합은 모르타르 연속성과 법선 도수 페널티를 조합하여 제4차 문제의 일致성 있는 이산화를 가능하게 한다.
• 이 방법은 제2차 및 제4차 고유값 문제에 적용되며, 특히 스펙트럼 오염을 방지하는 데 중점을 둔다.
• 이 방법은 다패치 기하 구조를 지원하며 선형 탄성 문제와 기르흐호프 플레이트 문제에 대해 검증된다.
• 형식은 경계 허브의 약한 부과를 통해 자연 경계 조건을 자연스럽게 통합한다.

실험 결과

연구 질문

• RQ1 고차수 페널티 항은 제2차 문제의 이소지오메트릭 고유값 근사에서 스펙트럼 오염을 효과적으로 제거할 수 있는가?
• RQ2 복잡한 경계를 가진 다패치 이소지오메트릭 설정에서 제4차 문제를 어떻게 정확하게 이산화할 수 있는가?
• RQ3 법선 도수를 페널티화하는 것은 이소지오메트릭 모르타르 방법의 정확성과 안정성에 어느 정도 기여하는가?
• RQ4 이 방법은 복잡한 매핑이 적용된 다패치 기하 구조에서도 고차수 수렴률을 유지할 수 있는가?
• RQ5 하이브리드 모르타르-페널티 접근법은 조건 수와 스펙트럼 성질 측면에서 기존의 이소지오메트릭 방법과 비교해 어떻게 다른가?

주요 결과

• 제안된 방법은 제2차 문제의 이산 고유값 근사에서 스펙트럼 오염('이상치')을 성공적으로 제거한다.
• 이 방법은 복잡한 매핑이 적용된 다패치 기하 구조에서 제4차 문제에 대해 안정적이고 정확한 해를 달성한다.
• 수치 실험을 통해 선형 탄성 문제와 기르흐호프 플레이트 문제에서 고차수 수렴률이 확인된다.
• 하이브리드 모르타르-페널티 형식은 네이만 경계 조건의 자연스러운 구현을 가능하게 한다.
• 체계적인 연구와 복잡한 기하 구조 테스트를 통해 이산 시스템의 강건한 성능과 양호한 조건 수를 입증한다.