[논문 리뷰] Behavior of different numerical schemes for population genetic drift problems

Minxin Chen, Chun Liu외 3명|arXiv (Cornell University)|2014. 10. 21.

Arctic and Antarctic ice dynamics참고 문헌 14인용 수 3

한 줄 요약

이 논문은 열화된 대류지배 유전적 드리프트 방정식을 해결하기 위해 세 가지 유한체적 스킴—상류화 및 두 가지 중심 스킴—을 분석한다. 유한체적 스킴 중에서 유량의 구조를 유지하는 중심 스킴만이 인위적 수치 점탄성을 피하고, 총 확률과 기댓값을 정확히 보존하는 반면, 상류화 스킴과 다른 중심 스킴은 유사한 점탄성을 유도하여 잘못된 정상 상태를 초래한다.

ABSTRACT

In this paper, we focus on numerical methods for the genetic drift problems, which is governed by a degenerated convection-dominated parabolic equation. Due to the degeneration and convection, Dirac singularities will always be developed at boundary points as time evolves. In order to find a \emph{complete solution} which should keep the conservation of total probability and expectation, three different schemes based on finite volume methods are used to solve the equation numerically: one is a upwind scheme, the other two are different central schemes. We observed that all the methods are stable and can keep the total probability, but have totally different long-time behaviors concerning with the conservation of expectation. We prove that any extra infinitesimal diffusion leads to a same artificial steady state. So upwind scheme does not work due to its intrinsic numerical viscosity. We find one of the central schemes introduces a numerical viscosity term too, which is beyond the common understanding in the convection-diffusion community. Careful analysis is presented to prove that the other central scheme does work. Our study shows that the numerical methods should be carefully chosen and any method with intrinsic numerical viscosity must be avoided.

연구 동기 및 목표

• 열화된 대류지배 타원형 편미분방정식에 의해 지배되는 유전적 드리프트 시뮬레이션에서 총 확률과 기댓값을 정확히 보존하는 수치 스킴을 특정하는 것.
• 다양한 유한체적 스킴의 장기적 행동을 분석하며, 특히 수치 점탄성에 의해 유도되는 인위적 정상 상태에 초점하는 것.
• 모든 추가적인 무한소 확산이 동일한 잘못된 정상 상태 해를 유도함을 증명하여, 점탄성 없는 방법이 필수적임을 강조하는 것.
• 중심 스킴이 예상치 못한 수치 점탄성을 유도할 수 있음을 보여주어, 대류-확산 맥락에서 통상적인 이해에 도전하는 것.
• 특정 유량 기반 중심 스킴만이 디랙 특이성까지 포함한 완전한 해 구조를 유지함을 확립하는 것.

제안 방법

• 세 가지 유한체적 스킴을 사용: 하나의 상류화 스킴과 두 가지 중심 차분 방법으로, 모두 유전적 드리프트 방정식의 유량을 이산화하기 위해 설계됨.
• 유량이 분리된 대류 및 확산 항으로 나뉘는 경우(스킴 1 및 2) 또는 전체 유량으로 간주되는 경우(스킴 3)를 포함한 유한체적 설정을 적용.
• 모든 스킴에 대해 무조건적 안정성을 증명하기 위해 이산 최대원리와 에너지 추정을 적용.
• 이산 헬더 및 푸앵카레 유형 부등식을 사용하여 해의 감쇠에 대한 해석적 경계를 유도하여 장기적 행동을 연구.
• 모든 추가적인 무한소 확산이 동일한 인위적 정상 상태 해를 유도함을 증명하며, 스킴에 관계없이 동일한 결과를 얻음.
• 수치 점탄성 항의 구조를 분석하여, 스킴 2는 이차 수준의 수치 점탄성을 유도하는 반면, 스킴 1은 일차 수준의 점탄성을 유도함을 밝힘.

실험 결과

연구 질문

• RQ1 유전적 드리프트 방정식의 해에서 총 확률과 기댓값을 정확히 보존할 수 있는 수치 스킴는 무엇인가?
• RQ2 상류화 스킴과 일부 중심 스킴는 안정함에도 불구하고 왜 인위적 정상 상태를 생성하는가?
• RQ3 특히 중심 스킴에 의한 수치 점탄성은 열화된 대류지배 편미분방정식의 해에 장기적 행동에 어떤 영향을 미치는가?
• RQ4 중심 차분 스킴는 어떤 조건에서라도 수치 점탄성을 유도하지 않을 수 있는가?
• RQ5 디랙 특이성이 있는 문제에서 정확한 장기적 행동을 확보하기 위해 유량 이산화의 어떤 구조적 특성이 필수적인가?

주요 결과

• 모든 세 스킴이 무조건적 안정성을 보이며 총 확률을 보존하지만, 오직 유량 기반 중심 스킴(스킴 3)만이 시간에 따라 기댓값을 정확히 보존한다.
• 상류화 스킴은 일차 수준의 수치 점탄성을 유도하여 안정성에도 불구하고 잘못된 인위적 정상 상태를 초래한다.
• 중심 스킴 중 하나(스킴 2)는 일반적으로 중심 스킴가 점탄성이 없음을 믿는 관념과는 반대로 이차 수준의 수치 점탄성 항을 유도한다.
• 모든 추가적인 무한소 확산, 크기가 얼마든 간에 동일한 인위적 정상 상태 해를 유도하므로, 이러한 확산은 본질적으로 완전한 해와 호환되지 않음을 증명한다.
• 스킴 3는 중심 차분을 사용해 유량을 전체로 취급함으로써 유사한 점탄성을 피하고, 디랙 특이성이 포함된 해석적 해와 일치하는 정확한 장기적 행동을 도출한다.
• 본 연구는 표준 대류-확산 수치 방법이 이 문제에 대해 본질적으로 실패함을 드러내며, 이는 방법 선택에 신중함이 필요함을 시사한다.