[논문 리뷰] Preconditioning a coupled model for reactive transport in porous media

Laila Amir, Michel Kern|arXiv (Cornell University)|2017. 10. 04.

Numerical methods for differential equations참고 문헌 66인용 수 4

한 줄 요약

이 논문은 다공성 매질 내 결합된 반응성 이동 문제를 해결하기 위해 이동 가능한 종의 농도를 제거하고 총 고정 농도만을 풀이 대상으로 삼는 전처리가 적용된 뉴턴–크릴로프 방법을 제시한다. 이 방법은 블록 가우스–자이델 또는 블록 재귀 전처리를 사용하여 비선형 및 선형 해법에서 메쉬에 독립적인 수렴성을 달성하며, 기존 방법에 비해 효율성이 크게 향상된다.

ABSTRACT

We study numerical methods for solving reactive transport problems in porous media that allow a separation of transport and chemistry at the software level, while keeping a tight numerical coupling between both subsystems. After recalling how each individual subsystem can be solved, we give a formulation that eliminates the local chemical concentrations, and keeps the total concentrations as unknowns. The block structure of the model is exploited both at the nonlinear level, by eliminating some unknowns, and at the linear level by using block Gauss-Seidel or block Jacobi preconditioning. The methods are applied to the easy 1D case of the MoMaS benchmark.

연구 동기 및 목표

• 전체 음성 방법의 계산 비용을 줄이고 수렴 속도를 높임으로써 반응성 이동의 효율성을 향상시키기.
• 이동 가능한 종의 농도를 제거하고 총 고정 농도만을 풀이 대상으로 삼는 수식을 개발하여 시스템을 단순화하기.
• 메쉬에 독립적인 수렴 성질을 유지하는 블록 전처리 기법을 분석하고 적용하기.
• 기본 1차원 반응성 이동 문제에 대해 제안된 방법의 강건성과 확장성을 입증하기.
• 소프트웨어 결합 최소화와 높은 정확도를 확보하면서도 복잡한 반응성 이동을 효율적으로 시뮬레이션할 수 있도록 하기.

제안 방법

• 이 방법은 운반 및 화학 반응 방정식의 완전히 결합된 비선형 시스템을 해결하기 위해 뉴턴–크릴로프 솔버를 사용한다.
• 지역 화학 농도를 제거하고 총 농도만을 미지수로 유지함으로써 문제를 재구성한다.
• 뉴턴 방법에서 유도된 선형 시스템에 대해 블록 가우스–자이델 및 블록 재귀 전처리를 적용한다.
• 전처리 전략이 비선형 수준에서 변수 제거와 동일한 효과를 가지며, 조건 수를 향상시킴을 보여준다.
• 벡터 곱을 통해 야코비안 행렬을 정확히 계산하며, 이는 운반 및 화학 솔버의 야코비안에 접근이 가능해야 한다.
• mass action 법칙을 직접 대입하지 않아 모듈성을 유지하면서도 정확도를 향상시킨다.

실험 결과

연구 질문

• RQ1 블록 전처리 전략이 반응성 이동 시뮬레이션에서 메쉬에 독립적인 수렴을 달성할 수 있는가?
• RQ2 이동 가능한 종의 농도를 제거하면 더 나은 조건을 갖춘 선형 시스템이 되어 수렴 속도가 빨라지는가?
• RQ3 전처리가 적용된 뉴턴–크릴로프 방법의 반복 횟수와 CPU 시간은 기존 방법과 비교해 어떻게 되는가?
• RQ4 실제 반응성 이동 기준 문제에 적용했을 때 이 방법이 강건성과 효율성을 유지할 수 있는가?
• RQ5 제안된 수식은 정확도를 유지하면서 계산 비용을 얼마나 줄일 수 있는가?

주요 결과

• 다양한 메쉬 크기에서 비선형 및 선형 반복 횟수가 거의 일정하게 유지되어 메쉬에 독립적인 수렴이 이루어짐을 나타낸다.
• 블록 가우스–자이델 및 블록 재귀 전처리기는 전처리 없을 경우에 비해 크릴로프 반복 횟수를 크게 감소시킨다.
• 메쉬 정밀도를 높여도 야코비안 행렬의 값 도메인은 0에서 멀리 떨어져 있어 메쉬에 독립적인 수렴을 뒷받침한다.
• 이동 가능한 종의 농도를 제거하는 대안적 수식은 다른 전처리 전략에 비해 CPU 시간을 줄인다.
• 이 방법은 큰 시간 간격으로도 수렴을 달성하여 연산자 분할 방법에서 흔히 발생하는 질량 오류를 피할 수 있다.
• MoMaS 기준 문제에 대한 수치 결과는 기존 문헌 데이터와 강한 일치를 보이며, 이 방법의 정확도를 검증한다.