[논문 리뷰] Flux Recovery and Superconvergence of Quadratic Immersed Interface Finite Elements

So‐Hsiang Chou, Champike Attanayake|arXiv (Cornell University)|2015. 12. 14.

Lattice Boltzmann Simulation Studies참고 문헌 7인용 수 4

한 줄 요약

이 논문은 연속하지 않은 확산 계수를 가진 1차원 인터페이스 문제에 적용된 2차 임mersed 유한요소 방법을 위한 플럭스 복원 기법을 제시한다. 이 방법은 흡수 항이 존재하지 않을 경우 플럭스와 압력에 대해 끝 노드와 인터페이스 점에서 초수렴을 달성하며, 이 경우 정확한 복원(오차 0)을 보이며, 흡수 항이 존재할 경우 플럭스 오차의 정밀도가 3차가 된다. 이는 수치 실험을 통해 확인되었다.

ABSTRACT

We introduce a flux recovery scheme for the computed solution of a quadratic immersed finite element method. The recovery is done at nodes and interface point first and by interpolation at the remaining points. We show that the end nodes are superconvergence points for both the primary variable $p$ and its flux $u$. Furthermore, in the case of piecewise constant diffusion coefficient without the absorption term the errors at end nodes and interface point in the approximation of $u$ and $p$ are zero. In the general case, flux error at end nodes and interface point is third order. Numerical results are provided to confirm the theory.

연구 동기 및 목표

• 추가 시스템을 풀지 않고도 정확한 플럭스 근사가 가능한 2차 임머스드 유한요소 해를 위한 플럭스 복원 기법을 개발하기 위해.
• 끝 노드와 인터페이스 점에서 복원된 플럭스와 압력의 초수렴 성질을 조사하기 위해.
• 복원된 플럭스에 대해 연속적인 첫 번째 기본정리와 유사한 이산 보존 법칙을 수립하기 위해.
• 다양한 메쉬 크기와 인터페이스 조건 하에서 이론적 결과를 수치 실험으로 확인하기 위해.

제안 방법

• 요소 간격에서 유한요소 해를 가중치 적분을 통해 플럭스 복원 기법을 구성하였으며, 국소 지지 집합을 가지며 특정 노드 값이 주어진 시험 함수를 활용하였다.
• 이 기법은 2차 임머스드 유한요소 기저 함수를 사용하여 인터페이스를 넘는 연속성과 점프 조건을 강제함으로써 약한 형태와의 일致성을 확보하였다.
• 부분적분과 점프 조건을 활용한 오차 분석을 통해 주요 점에서 압력과 플럭스에 대한 점별 오차 추정치를 도출하였다.
• 흡수 항 q(x) = 0일 경우, 기초 갈레르킨 사영의 정확성 덕분에 끝 노드에서 플럭스와 압력의 정확한 복원이 보장된다.
• 이산 보존 법칙을 유도하였으며, 이는 인접한 노드에서 복원된 플럭스의 차이가 요소 내에서 소스 항의 적분과 같다는 것을 보여준다.
• 다양한 메쉬 크기와 다항식 차수를 가진 두 가지 벤치마크 문제를 사용하여 수치적 검증을 수행하였으며, 노드와 인터페이스 점에서의 오차를 비교하였다.

실험 결과

연구 질문

• RQ1 2차 임머스드 유한요소에 대해 끝 노드와 인터페이스 점에서 초수렴을 달성할 수 있는 플럭스 복원 기법을 구성할 수 있는가?
• RQ2 흡수 항 q(x)가 0이 아닐 경우, 끝 노드와 인터페이스 점에서 플럭스 오차의 수렴 차수는 얼마인가?
• RQ3 복원된 플럭스는 연속적인 첫 번째 기본정리와 유사한 이산 보존 법칙을 만족하는가?
• RQ4 q(x) = 0일 경우 복원된 플럭스가 끝 노드와 인터페이스 점에서 정확한가, 그리고 만약 그렇다면 어떤 조건에서 성립하는가?

주요 결과

• q(x) = 0일 경우, 끝 노드와 인터페이스 점에서 플럭스 오차는 정확히 0이 되며, 이는 플럭스의 정확한 복원을 의미한다.
• 일반적인 경우 q(x) ≠ 0일 경우, 끝 노드와 인터페이스 점에서 플럭스 오차는 3차 정밀도를 가지며, 즉 O(h³)이다.
• q(x) = 0일 경우, 끝 노드에서 압력은 초수렴하며, 이 점에서의 오차는 정확히 0이다.
• 수치 결과는 q(x) ≠ 0일 경우 끝 노드와 인터페이스 점에서 플럭스 오차의 수렴률이 약 O(h³)임을 확인한다.
• 중점 노드와 인터페이스 점에서의 압력에 대해서는 q(x) = 0일 경우 약 O(h⁴)의 수렴률을 보이며, q(x) ≠ 0일 경우 O(h²)의 수렴률을 보이며, 이는 이론적 예측과 일치한다.
• 이산 보존 법칙이 성립한다: uₕ(tᵢ) - uₕ(tᵢ₋₁) = ∫_{Iᵢ} f(x) dx. 이는 연속적인 보존 성질을 그대로 반영한다.