[논문 리뷰] $h-p$ spectral element methods for three dimensional elliptic problems on non-smooth domains using parallel computers

Akhlaq Husain|arXiv (Cornell University)|2011. 10. 11.

Advanced Numerical Methods in Computational Mathematics참고 문헌 37인용 수 9

한 줄 요약

이 논문은 비연속 영역에서 3차원 타원 경계값 문제를 지수 수렴성으로 해결하기 위한 비연속 h-p 스펙트럼 유한요소 방법을 제안한다. 특이점 근처에서 국소적인 구좌표계와 원통좌표계를 사용하고, 모서리와 가장자리 쪽으로 수축하는 기하학적 메쉬를 적용함으로써, 최소제곱 콜로케이션 방법을 통해 조건부 공액 경사법(Preconditioned Conjugate Gradient Method)으로 해를 구함으로써 명시적인 질량행렬과 강성행렬의 조립을 피함으로써 최적의 수렴 속도를 달성한다.

ABSTRACT

Elliptic partial differential equations arise in many fields of science and engineering such as steady state distribution of heat, fluid dynamics, structural/mechanical engineering, aerospace engineering and seismology etc. In three dimensions it is well known that the solutions of elliptic boundary value problems have singular behavior near the corners and edges of the domain. The singularities which arise are known as vertex, edge, and vertex-edge singularities. We propose a nonconforming h-p spectral element method to solve three dimensional elliptic boundary value problems on non-smooth domains to exponential accuracy. To overcome the singularities which arise in the neighbourhoods of the vertices, vertex-edges and edges we use local systems of coordinates. These local coordinates are modified versions of spherical and cylindrical coordinate systems in their respective neighbourhoods. Away from these neighbourhoods standard Cartesian coordinates are used. In each of these neighbourhoods we use a geometrical mesh which becomes finer near the corners and edges. We then derive differentiability estimates in these new set of variables and a stability estimate on which our method is based for a non-conforming h-p spectral element method. The Sobolev spaces in vertex-edge and edge neighbourhoods are anisotropic and become singular at the corners and edges. The method is essentially a least-squares collocation} method and a solution can be obtained using Preconditioned Conjugate Gradient Method (PCGM). To solve the minimization problem we need to solve the normal equations for the least-squares problem. The residuals in the normal equations can be obtained without computing and storing mass and stiffness matrices. Computational results for a number of model problems confirm the theoretical estimates obtained for the error and computational complexity.

연구 동기 및 목표

• 모서리, 가장자리, 정점-모서리 영역 근처에서 3차원 타원 PDE 해의 특이 행동을 다루기 위해.
• 해의 특이성에도 불구하고 지수 수렴성을 달성할 수 있는 고차수 수치 방법을 개발하기 위해.
• 대규모 문제를 효율적으로 병렬 계산할 수 있도록 비연속 스펙트럼 유한요소 설정을 활용하기 위해.
• 해 과정에서 질량행렬과 강성행렬의 조립 및 저장 비용을 피하기 위해.
• 모델 문제에 대한 수치 실험을 통해 이론적 오차 추정치와 계산 복잡도를 검증하기 위해.

제안 방법

• 특이 영역(정점, 가장자리, 정점-모서리) 근처에서 수정된 구좌표계와 원통좌표계를 사용하여 특이성을 더 정확히 해석한다.
• 이러한 국소 영역에 기하학적 메쉬를 적용하여 특이점 쪽으로 점점 더 미세해지도록 하여 해의 이방향성 행동을 포착한다.
• 특이점에서 벗어난 영역에서는 표준 카르테시안 좌표계를 사용하여 표준 스펙트럼 유한요소 절차와의 호환성을 확보한다.
• 수식은 최소제곱 콜로케이션 방법에 기반하며, 명시적인 행렬 조립 없이 약한 의미에서 잔차를 최소화한다.
• 특이성의 구조를 반영하는 이방향 소볼레프 공간에서 안정성과 미분 가능성 추정치를 도출한다.
• 최종 선형 시스템은 조건부 공액 경사법(Preconditioned Conjugate Gradient Method)을 사용하여 해석하며, 질량행렬이나 강성행렬을 저장하지 않고도 정규방정식에서 잔차를 직접 계산한다.

실험 결과

연구 질문

• RQ1 h-p 스펙트럼 유한요소 방법은 정점과 가장자리 특이성이 존재하는 비연속 영역에서 3차원 타원 문제에 대해 지수 수렴성을 달성할 수 있는가?
• RQ2 국소 좌표 변환은 3차원 영역에서 특이 해의 해상도를 어떻게 향상시킬 수 있는가?
• RQ3 특이점 근처의 이방향 메쉬 정밀화가 수렴 속도와 계산 비용에 어떤 영향을 미치는가?
• RQ4 최소제곱 콜로케이션 접근법은 큰 질량행렬과 강성행렬의 조립 및 저장 없이 효율적으로 구현될 수 있는가?
• RQ5 이론적 오차 추정치는 실질적인 대규모 계산에서 수치 결과와 어떻게 비교되는가?

주요 결과

• 이 방법은 정점, 가장자리, 정점-가장자리 특이성이 있는 문제에 대해 h-p 의미에서 지수 수렴 속도를 달성하여 이론적 예측을 확인한다.
• 국소적인 구좌표계와 원통좌표계의 사용은 기하학적 특이점 근처에서 특이 해의 해상도를 크게 향상시킨다.
• 계산 결과는 유도된 이론적 오차 추정치와 일치하는 최적의 수렴 행동을 보여준다.
• 이 방법은 질량행렬과 강성행렬의 저장 및 조립을 피함으로써 메모리 오버헤드를 줄이고 병렬 아키텍처에서의 확장성을 향상시킨다.
• 조건부 공액 경사법은 최소제곱 수식에서 유도된 대규모 선형 시스템을 효과적으로 해석한다.
• 모델 문제에 대한 수치 실험은 이론적 오차 한계와 계산 복잡도 추정치를 모두 검증하며, 다양한 특이 구성에 대해 안정적인 성능을 보여준다.