[논문 리뷰] Finite Element Technique for Solving the Stream Function Form of a Linearized Navier-Stokes Equations Using Argyris Element

Faisal Fairag, Noah Al-Mulla|ArXiv.org|2004. 06. 03.

Advanced Numerical Methods in Computational Mathematics참고 문헌 12인용 수 3

한 줄 요약

이 논문은 선형화된 나비에-스토크스 방정식의 스트림 함수 형태를 해결하기 위해 아르기리스 원소를 사용하는 유한요소법을 제시하며, 계산 효율성을 향상시키기 위해 행렬의 폭을 최소화하기 위해 노드 재정렬을 중시한다. 이 방법은 최적의 수렴성과 안정성을 달성하며, 유닛 정사각형 영역에서 Re=1인 경우 제조된 정확한 해를 사용한 수치 결과는 정확한 해를 보여준다. 적응형 사다리꼴 적분 규칙을 사용한다.

ABSTRACT

The numerical implementation of finite element discretization method for the stream function formulation of a linearized Navier-Stokes equations is considered. Algorithm 1 is applied using Argyris element. Three global orderings of nodes are selected and registered in order to conclude the best banded structure of matrix and a fluid flow calculation is considered to test a problem which has a known solution. Visualization of global node orderings, matrix sparsity patterns and stream function contours are displayed showing the main features of the flow.

연구 동기 및 목표

• 선형화된 나비에-스토크스 방정식의 스트림 함수 수식을 위한 안정적이고 효율적인 유한요소 이산화를 개발한다.
• 다양한 전역 노드 순서가 유한요소 행렬의 폭과 흐린 구조에 미치는 영향을 평가한다.
• 비대칭이고 흐린 선형 시스템을 해결하기 위해 켈리지안 스테이블라이즈드(BiCGSTAB) 방법을 구현하고 테스트한다.
• 유닛 정사각형 영역에서 알려진 해석적 행동을 보이는 제조된 정확한 해를 사용하여 수치적 방법을 검증한다.
• 다양한 메쉬 정밀도와 적분 규칙에서 계산 비용, 수렴 속도 및 반복 솔버 성능을 분석한다.

제안 방법

• 스트림 함수 수식은 압력을 제거하고 비압축성을 자동으로 만족시켜, 네 번째 차수의 비하모닉형 편미분방정식으로 문제를 축소시킨다.
• 스옵레브 공간 $ H^2_0(\Omega) $ 를 사용하여 선형화된 문제의 약한 형태를 유도하며, 이중형 형태 $ a(\psi,\varphi) $ 와 $ b(\psi^*;\psi,\varphi) $ 는 라플라스 연산자와 대류 항을 나타낸다.
• 아르기리스 원소는 $ H^2_0(\Omega) $ 에서의 연속 유한요소로 사용되며, 네 번째 차수 문제에 필요한 $ C^1 $ 연속성을 보장한다.
• 유한요소 행렬의 폭을 최소화하고 강성 행렬의 조건수를 개선하기 위해 세 가지 전역 노드 재정렬 전략을 시험한다.
• BiCGSTAB 반복 솔버는 유한요소 수식에서 유도된 흐린 비대칭 선형 시스템에 적용되며, 잔차와 반복값 노름의 허용 오차를 통해 수렴을 모니터링한다.
• 요소 행렬을 평가하는 데 있어 정확도와 계산 비용을 평가하기 위해 네 점과 여섯 점의 가우스 적분 점을 사용한다.

실험 결과

연구 질문

• RQ1 노드 재정렬은 $ H^2 $-연속 이산화에서 유한요소 행렬의 폭과 흐린 패턴에 어떤 영향을 미치는가?
• RQ2 4점과 6점의 적분 점 선택이 아르기리스 유한요소 해의 정확도와 계산 비용에 어떤 영향을 미치는가?
• RQ3 메쉬 정밀도($ h = 1/3, 1/5, 1/9 $) 가 증가함에 따라 BiCGSTAB 반복 횟수와 CPU 시간은 어떻게 변화하는가?
• RQ4 제조된 정확한 해를 사용할 때 아르기리스 유한요소 방법은 스트림 함수 수식에서 최적의 수렴 속도를 달성할 수 있는가?
• RQ5 비선형 방정식의 초기 추정치 $ \psi^* $ 를 비하모닉 방정식을 풀어 얻었을 때 반복 솔버의 수렴에 어떤 영향을 미치는가?

주요 결과

• 노드 재정렬 1이 가장 좁은 행렬 폭과 가장 낮은 CPU 시간(28.45초)을 기록하여 다른 방법들보다 계산 효율성이 뛰어났다.
• $ h = 1/9 $ 인 경우, $ L^2 $-노름에서 상대 오차 $ 1.423 \times 10^{-4} $ 와 $ H^1 $-노름에서 $ 6.986 \times 10^{-4} $ 를 기록하였으며, BiCGSTAB 반복 횟수는 567.5회였다.
• 산술 연산 횟수는 메쉬 정밀도 증가에 따라 크게 증가하였으며, $ h=1/3 $ 에서 약 1750만 회에서 $ h=1/9 $ 에서 약 35억 회로 증가하였으며, 이는 $ H^2 $-연속 원소의 $ O(h^{-4}) $ 스케일링을 반영한다.
• 여섯 점 적분을 사용할 경우 네 점 적분에 비해 오차와 연산 횟수가 약간 증가했지만, 안정성과 수렴 행동이 향상되었다.
• 유선 시각화 결과 모든 메쉬 크기에서 물리적으로 일관된 유동 패턴을 확인하였으며, 이는 방법이 소용돌이 구조를 잘 포착할 수 있음을 검증한다.
• 모든 테스트된 메쉬 크기와 적분 규칙에서 방법은 강력한 수렴성을 보였으며, 잔차와 반복값 노름이 모두 $ 10^{-5} $ 이하로 감소하였다.