[논문 리뷰] A higher-order finite-volume discretization method for Poisson's equation in cut cell geometries

Dharshi Devendran, Daniel Graves외 1명|arXiv (Cornell University)|2014. 11. 16.

Advanced Numerical Methods in Computational Mathematics참고 문헌 15인용 수 4

한 줄 요약

이 논문은 가중치가 부여된 최소 제곱 재구성 기반의 면 유량 재구성 기법을 사용하여 카르테시안 컷 셀 격자에서 푸아송 방정식에 대해 고차수 유한체적 방법을 제안한다. 이 방법은 제2차 및 제4차 정확도를 달성하며, 신중하게 선택된 가중치를 통해 안정성을 확보한다. 이는 다양한 2차원 및 3차원 기하 구조에서 이산 라플라스 연산자가 엄격히 음수인 고유값을 가지며 최적의 수렴성을 보여준다.

ABSTRACT

We present a method for generating higher-order finite volume discretizations for Poisson's equation on Cartesian cut cell grids in two and three dimensions. The discretization is in flux-divergence form, and stencils for the flux are computed by solving small weighted least-squares linear systems. Weights are the key in generating a stable discretization. We apply the method to solve Poisson's equation on a variety of geometries, and we demonstrate that the method can achieve second and fourth order accuracy in both truncation and solution error for these examples. We also show that the Laplacian operator has only stable eigenvalues for each of these examples.

연구 동기 및 목표

• 복잡한 기하 구조에서 컷 셀을 사용한 푸아송 방정식에 대해 고차수 정확도를 달성하는 데 도전하는 것.
• 기존의 컷 셀 유한체적 방법에서 임베디드 경계 근처에서 제1차 절단 오차를 초래하는 한계를 극복하는 것.
• 작은 컷 셀 부피와 복잡한 인터페이스 조건에서도 고차수 정확도를 유지하는 안정적이고 보존적인 이산화 방법을 개발하는 것.
• 다양한 기하 구조에서 이산 라플라스 연산자가 엄격히 음수인 고유값만을 가지도록 보장하여 안정성을 확보하는 것.
• 메시팅 복잡도를 최소화하면서도 보존이 필요한 문제, 예를 들어 열 및 질량 이동 문제에 대한 강력한 적용을 가능하게 하는 것.

제안 방법

• 이산 보존성을 체계적으로 확보하기 위해 푸아송 방정식을 유량 발산 형태로 재구성한다.
• 각 컷 셀에서 해의 다항식 재구성을 통해 가중치가 부여된 최소 제곱 시스템을 사용해 면 유량을 계산한다.
• 해의 안정성과 조건수 제어를 위해 가중치를 사용하여 과소결정 선형 시스템을 해결한다.
• 지역 기하 구조와 부피 분율에 맞게 스텐실을 조정하면서 2차원 및 3차원 카르테시안 컷 셀 격자에 이 방법을 적용한다.
• 보존성과 고차수 정확도를 동시에 확보하기 위해 일관되고 보존적인 유량 재구성 기법을 사용한다.
• 가상 셀과 일관된 재구성 기법을 사용해 임베디드 경계에서 디리클레 및 뉴먼 경계 조건을 적용한다.

실험 결과

연구 질문

• RQ1 커트 셀 격자에서 푸아송 방정식에 대해 제2차 및 제4차 정확도를 달성하는 고차수 유한체적 방법을 구성할 수 있는가? 이는 절단 오차와 해 오차 양쪽에서 수렴 정도가 제2차 및 제4차 수준에 도달하는가?
• RQ2 적절한 가중치를 사용한 과소결정 최소 제곱 시스템을 통해 작은 컷 셀에 대해 안정적인 유량 재구성을 어떻게 생성할 수 있는가?
• RQ3 결과적으로 도출된 이산 라플라스 연산자가 다양한 복잡한 기하 구조에서 엄격히 음수인 고유값만을 가지며 안정성을 유지하는가?
• RQ4 소규모 기하학적 변형, 예를 들어 임베디드 경계 위치의 약간의 이동에 대해 이 방법은 얼마나 강인한가?
• RQ5 이 방법은 최소한의 수정으로 변수 계수 푸아송 방정식 및 기타 발산-유량 형태 문제로 확장 가능한가?

주요 결과

• 2차원 정사각형 기하 구조에 내부로 네 개의 원이 있는 경우, 해 오차에서 제4차 수렴을 달성한다. 그리드 해상도가 $h=1/32$ 에서 $h=1/64$로 증가함에 따라 $L_\text{infty}$ 오차는 $1.04\times10^{-6}$ 에서 $6.30\times10^{-8}$ 로 감소한다.
• 3차원 구 기하 구조의 경우, 해 오차는 약 제4차로 수렴하며, $h$ 가 $1/64$ 에서 $1/128$ 로 반으로 줄어들면서 $L_\infty$ 오차는 $9.16\times10^{-5}$ 에서 $6.27\times10^{-6}$ 로 감소한다.
• 모든 테스트된 기하 구조에서 이산 라플라스 연산자가 엄격히 음수인 고유값만을 가지며, 디리클레 및 뉴먼 경계 조건 모두에서 수치적 안정성이 확인된다.
• 소규모 기하학적 변형에 대해서도 방법은 강인하다: 2차원에서 원의 중심을 $0.001$ 만큼 이동하거나 반지름을 $0.005$ 만큼 변경해도 수렴 속도가 떨어지지 않는다.
• 라플라스 연산자의 스펙트럼에는 비영인 허수부를 가진 작은 고유값 군이 존재하지만, 주요 고유값은 실수이면서 음수로 유지되어 안정성을 유지한다.
• 2차원 정사각형 기하 구조에 내부로 네 개의 원이 있는 경우, 절단 오차에서 제2차 수렴을 달성한다. $h$ 가 반으로 줄어들면서 $L_1$ 절단 오차는 $1.12\times10^{-3}$ 에서 $4.62\times10^{-5}$ 로 감소한다.