[논문 리뷰] Numerical method for Darcy flow derived using Discrete Exterior Calculus

Anil N. Hirani, K. B. Nakshatrala외 1명|arXiv (Cornell University)|2008. 10. 19.

Black Holes and Theoretical Physics참고 문헌 36인용 수 12

한 줄 요약

이 논문은 단순형 메esh에서 유량과 압력 변수를 통해 혼합형 포아송 방정식을 이산화하는 좌표 불변 수치 방법을 제시한다. 이 방법은 국소 질량 보존을 보장하며, 유량과 압력에 대해 최적 수렴률을 달성하고, 대각형이면서 양의 정부호인 허지 스타 이산화를 통해 공간적으로 변화하는 투수도를 처리한다. 2차원, 3차원 및 곡면 문제에서의 성능이 입증되었다.

ABSTRACT

We derive a numerical method for Darcy flow, hence also for Poisson's equation in mixed (first order) form, based on discrete exterior calculus (DEC). Exterior calculus is a generalization of vector calculus to smooth manifolds and DEC is one of its discretizations on simplicial complexes such as triangle and tetrahedral meshes. DEC is a coordinate invariant discretization, in that it does not depend on the embedding of the simplices or the whole mesh. We start by rewriting the governing equations of Darcy flow using the language of exterior calculus. This yields a formulation in terms of flux differential form and pressure. The numerical method is then derived by using the framework provided by DEC for discretizing differential forms and operators that act on forms. We also develop a discretization for spatially dependent Hodge star that varies with the permeability of the medium. This also allows us to address discontinuous permeability. The matrix representation for our discrete non-homogeneous Hodge star is diagonal, with positive diagonal entries. The resulting linear system of equations for flux and pressure are saddle type, with a diagonal matrix as the top left block. The performance of the proposed numerical method is illustrated on many standard test problems. These include patch tests in two and three dimensions, comparison with analytically known solution in two dimensions, layered medium with alternating permeability values, and a test with a change in permeability along the flow direction. We also show numerical evidence of convergence of the flux and the pressure. A convergence experiment is also included for Darcy flow on a surface. A short introduction to the relevant parts of smooth and discrete exterior calculus is included in this paper. We also include a discussion of the boundary condition in terms of exterior calculus.

연구 동기 및 목표

• 국소 질량 보존과 기하학적 불변성을 유지하는 Darcy 유량을 위한 수치 방법을 개발한다.
• 비균일하고 공간적으로 변화하는 투수도를 다룰 수 있도록 이산 외부 미분 계산(DEC)을 확장한다.
• 2차원, 3차원 및 곡면에서의 Darcy 유량을 위한 통합된 수식을 메쉬 기반 이산화로 제시한다.
• 표준 테스트 문제에서 유량과 압력 변수에 대해 최적 수렴률을 입증한다.
• 유한요소 외부 미분 계산과의 호환성을 확보하고, 향후 비등방성 투수도로의 확장을 위한 프레임워크를 수립한다.

제안 방법

• Darcy 유량의 근본 방정식을 외부 미분 기하학을 사용해 재구성하여, 유량을 미분형으로, 압력을 0-형으로 표현한다.
• 이산 외부 미분 계산(DEC)을 적용하여 미분형과 연산자를 이산화함으로써 스토크스 정리와 같은 위상적 및 기하학적 성질을 유지한다.
• 공간적으로 변화하는 투수도를 위한 허지 스타 연산자의 새로운 이산화를 개발하여, 양의 대각형 행렬을 얻는다.
• 유량은 변(또는 면)과 연결되고 압력은 정점(또는 세포)과 연결되는 스태그너 메쉬 쌍을 사용하여 국소 질량 보존을 보장한다.
• 유도된 선형 시스템은 대각형 상단 왼쪽 블록을 지닌 사다리꼴 형식이며, 효율적인 해법 전략을 가능하게 한다.
• 경계 조건은 외부 미분 기하학의 언어로 구현되어 기하학적 구조와의 일致를 유지한다.

실험 결과

연구 질문

• RQ1 DEC 기반 방법이 일반 단순형 메쉬에서 Darcy 유량의 유량과 압력에 대해 최적 수렴률을 달성할 수 있는가?
• RQ2 비균일하고 공간적으로 변화하는 투수도를 고려할 때, 허지 스타 연산자를 어떻게 이산화할 수 있으며, 양의 정부호성과 대각형 형식을 유지할 수 있는가?
• RQ3 DEC 프레임워크는 임베딩에 의존하지 않고도 곡면에서 Darcy 유량을 다룰 수 있도록 확장될 수 있는가?
• RQ4 불연속적이거나 비등방성 투수도가 존재하는 상황에서도 방법이 국소 질량 보존을 유지하는가?
• RQ5 표준 벤치마크 문제에서 DEC 기반 방법이 기존의 유한요소 또는 유한체적 방법과 비교해 어떻게 성능을 발휘하는가?

주요 결과

• 모든 표준 패치 테스트를 2차원 및 3차원에서 통과하여 일致성과 안정성을 확인하였다.
• 유량과 압력에 대한 수치적 수렴이 관찰되었으며, 표면에서의 수렴 실험을 통해 최적 수렴률이 확인되었다.
• 비균일한 투수도를 위한 이산 허지 스타는 대각형이면서 양의 성분을 지녀, 잘 조절된 선형 시스템을 보장한다.
• DEC에서 정확한 이산 스토크스 정리의 성질 덕분에 방법은 구성상 국소 질량 보존을 유지한다.
• 수식은 내재적이며 좌표 불변하므로, 임베딩에 의존하는 수정 없이도 반구와 같은 곡면에 직접 적용 가능하다.
• 층상 구조에서 번갈아가며 투수도가 변화하는 매체와 흐름 방향으로 투수도 기울기가 존재하는 유량 문제에서 뛰어난 성능을 보였다.