[논문 리뷰] Strong form mesh-free $hp$-adaptive solution of linear elasticity problem

Mitja Jančič, Gregor Kosec|arXiv (Cornell University)|2022. 10. 13.

Advanced Numerical Analysis Techniques인용 수 4

한 줄 요약

이 논문은 선형 탄성 문제를 해결하기 위한 메쉬 없는 hp-적응형 수치 방법을 제시한다. 이 방법은 반경 기저 함수를 이용한 유한차분(RBF-FD)과 혁신적인 암묵-수동(IMEX) 오차 지표를 결합한 것이다. 이 방법은 반복적으로 해를 구하고, 국소적 고차수 재평가를 통해 오차를 추정하며, 수정된 텍사스 3단계 전략을 사용해 정밀화 대상 노드를 표시하고, 노드를 재배치하며 근사 차수를 조정함으로써 2차원 및 3차원 접촉 문제에서 비적응형 방법보다 뛰어난 수렴성과 정확도를 달성한다.

ABSTRACT

We present an algorithm for $hp$-adaptive collocation-based mesh-free numerical analysis of partial differential equations. Our solution procedure follows a well-established iterative solve-estimate-mark-refine paradigm. The solve phase relies on the Radial Basis Function-generated Finite Differences (RBF-FD) using point clouds generated by advancing front node positioning algorithm that supports variable node density. In the estimate phase, we introduce an Implicit-Explicit (IMEX) error indicator, which assumes that the error relates to the difference between the implicitly obtained solution (from the solve phase) and a local explicit re-evaluation of the PDE at hand using a higher order approximation. Based on the IMEX error indicator, the modified Texas Three Step marking strategy is used to mark the computational nodes for $h$-, $p$- or $hp$-(de-)refinement. Finally, in the refine phase, nodes are repositioned and the order of the method is locally redefined using the variable order of the augmenting monomials according to the instructions from the mark phase. The performance of the introduced $hp$-adaptive method is first investigated on a two-dimensional Peak problem and further applied to two- and three-dimensional contact problems. We show that the proposed IMEX error indicator adequately captures the global behaviour of the error in all cases considered and that the proposed $hp$-adaptive solution procedure significantly outperforms the non-adaptive approach. The proposed $hp$-adaptive method stands for another important step towards a fully autonomous numerical method capable of solving complex problems in realistic geometries without the user intervention.

연구 동기 및 목표

• 실제 기하 형상에서 복잡한 선형 탄성 문제를 해결하기 위한 완전 자동화된 메쉬 없는 hp-적응형 수치 방법을 개발하는 것.
• 산산이 흩어진 노드 점군과 적응형 정밀화 전략을 활용하여 메쉬 생성 및 정밀화에 대한 사용자 간섭을 제거하는 것.
• h-정밀화(노드 재배치)와 p-정밀화(변동 근사 차수)를 결합하여 해의 정확도와 수렴 효율을 향상시키는 것.
• 피크 문제와 분석적 해가 없는 3차원 펄서 피로 접촉 문제를 포함한 기준 문제에 대해 제안된 방법을 검증하는 것.
• 고기울기 및 특이점 영역에서 전반적인 오차 행동을 효과적으로 포착할 수 있는 새로운 IMEX 기반 오차 지표를 도입하고 평가하는 것.

제안 방법

• 진행면 노드 배치 알고리즘을 사용해 생성된 점군 위에서 RBF-FD를 사용한 반복적 해결-추정-표시-정밀화 프로토콜을 적용한다.
• 해결 단계에서는 강형 형태의 편미분 연산자를 다항기저 함수와 가변 스텐실 크기를 사용한 RBF-FD를 통해 근사한다.
• 추정 단계에서는 암묵적으로 구해진 해와 국소적으로 명시적인 고차수 재평가를 통해 PDE를 재계산한 결과의 차이를 계산하는 IMEX 오차 지표를 도입한다.
• 표시 단계에서는 수정된 텍사스 3단계 전략을 사용해 IMEX 오차 지표에 기반해 h-, p-, 또는 hp-정밀화를 위한 노드를 분류한다.
• 정밀화 단계에서는 진행면 알고리즘을 사용해 노드를 재배치하고, 표시 결정에 따라 국소적으로 다항식을 추가함으로써 근사 차수를 조정한다.
• 이 방법은 메쉬 구조의 위상에 의존하지 않는 통합 프레임워크 내에서 h-정밀화(노드 추가/제거)와 p-정밀화(차수 증가)를 모두 지원한다.

실험 결과

연구 질문

• RQ1 IMEX 기반 오차 지표는 메쉬 없는 RBF-FD 해에서 선형 탄성 문제의 전반적 오차 행동을 효과적으로 포착할 수 있는가?
• RQ2 제안된 hp-적응형 프레임워크는 비적응형 RBF-FD 방법에 비해 수렴 속도와 정확도를 크게 향상시키는가?
• RQ3 특이점이나 국소적 응력 집중이 있는 문제, 예를 들어 접촉 문제에 대해 이 방법은 얼마나 잘 성능을 발휘하는가?
• RQ4 사용자 정의 메쉬 또는 수동 정밀화 없이도 이 방법은 복잡한 3차원 기하 형상에서 높은 정확도를 달성할 수 있는가?
• RQ5 변동 노드 밀도와 변동 근사 차수의 조합이 날카로운 기울기와 경계층을 효과적으로 해결하는 데 효과적인가?

주요 결과

• IMEX 오차 지표는 2차원 피크 문제와 3차원 펄서 피로 접촉 문제를 포함한 모든 시험 케이스에서 전반적 오차 행동을 성공적으로 포착하였다.
• 제안된 hp-적응형 방법은 비적응형 RBF-FD에 비해 유의미하게 빠른 수렴과 높은 정확도를 달성하였으며, 특히 고응력 기울기 및 특이점 영역에서 두드러진 성능을 보였다.
• 펄서 피로 접촉 문제에서 이 방법은 도메인 폭에 비해 200배 작아진 0.2067 mm 크기의 접촉 영역을 고해상도로 해결하여 소규모 특징을 모델링할 수 있는 능력을 입증하였다.
• 이론적 연구 결과와 일치하게, 초기 반복 단계에서 특이점 근처에서 근사 차수를 자동으로 증가시킨 후 국소적 필드 기술 향상을 위해 h-정밀화를 수행하였다.
• 변동 노드 밀도와 적응형 차수 선택의 활용은 계산 자원의 보다 효율적인 분포를 이끌어내어 총 계산 비용을 감소시켰다.
• 모서리 노드를 제거하여 수치적 안정성을 확보한 경우에도 이 방법은 안정성과 강건성을 유지하였으며, 실제 문제에 대한 실용적 타당성을 확인하였다.