[논문 리뷰] Lie group computation of finite difference schemes

Emma Hoarau, Claire David|ArXiv.org|2006. 11. 29.

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

한 줄 요약

이 논문은 유한 차분 스킴의 미분 근사에 기반하여 Mathematica 기반의 기호 계산 패키지를 제시한다. 이는 유한 차분 스킴의 리 군 대칭을 결정하는 데 사용된다. 스킴의 미분 표현에 리 군 분석을 적용함으로써, 유지되거나 상실되는 대칭을 식별한다. 이 방법은 FTCS 스킴이 시간 및 공간 이동 대칭과 수정된 확대군을 유지하지만, 갈릴레이 및 사영 불변성을 상실함을 보여준다.

ABSTRACT

A Mathematica based program has been elaborated in order to determine the symmetry group of a finite difference equation, by means of its differential representation. The package provides functions which enable us to solve the determining equations of the related Lie group

연구 동기 및 목표

• 유한 차분 스킴의 리 군 분석을 자동화하기 위한 기호 계산 도구를 개발하는 것.
• 원래 PDE의 연속 대칭 중에서 수치 이산화 과정에서 유지되거나 깨지는 대칭을 규명하는 것.
• 군 이론적 방법을 사용하여 대칭을 유지하는 유한 차분 스킴을 체계적으로 구성할 수 있도록 하는 것.
• 미분 근사 차수에 따라 대칭 결정의 계산 비용이 어떻게 변하는지 평가하는 것.

제안 방법

• 유한 차분 스킴의 미분 근사를 기반으로 리 군 이론을 적용하여, 이를 PDE 시스템으로 간주한다.
• 벡터 장의 임계 확장에 기반한 결정 방정식을 유도하고, 이를 통해 무한소 생성자를 계산한다.
• 독립 변수, 종속 변수 및 스텝 크기 변수(x, t, u, h, τ, ν)의 다항식 전개를 사용하여 무한소를 매개변수화한다.
• 패키지는 Mathematica의 기호 연산 능력을 활용하여 대수적 결정 방정식의 해를 자동으로 구한다.
• 결정 시스템의 해로부터 대칭 군을 재구성하며, 결과는 표준 스킴에서 검증된다.
• 실행 시간 측정을 통해 미분 근사 차수에 따라 변하는 계산 성능을 평가한다.

실험 결과

연구 질문

• RQ1 주어진 유한 차분 스킴이 원래 PDE의 어떤 연속 대칭을 유지하는가?
• RQ2 미분 근사의 차수는 대칭 결정의 계산 복잡도에 어떤 영향을 미치는가?
• RQ3 체계적인 기호 패키지가 유한 차분 스킴의 리 대칭 계산을 자동화할 수 있는가?
• RQ4 특히 FTCS 및 라크스-벤드로프와 같은 표준 스킴에서 이산화 과정에서 상실되는 대칭은 무엇인가?
• RQ5 스텝 크기 변수(h, τ, ν)는 대칭 군에 의해 어떻게 변환되며, 이는 스킴의 불변성에 어떤 함의를 갖는가?

주요 결과

• FTCS 스킴은 시간 및 공간 이동 대칭과 x, t, u, h, τ를 포함한 두 개의 확대 대칭으로 생성되는 4차원 리 군을 유지한다.
• FTCS 스킴의 대칭 대수는 A_d ⊕ A_{h,τ}와 동형이며, 여기서 A_d는 버거스 방정식의 확대 대수이고, A_{h,τ}는 스텝 크기에 작용하는 대칭이다.
• FTCS 스킴은 이산화 과정으로 인해 갈릴레이 대칭과 사영 불변성을 상실하여 대칭 깨짐이 발생한다.
• 라크스-벤드로프 및 크랭크-니콜슨 스킴도 동일한 4매개변수 대칭 군을 가지며, 유사한 대칭 구조를 가짐을 시사한다.
• 미분 근사 차수에 따라 실행 시간이 급격히 증가한다: 차수 4일 때 FTCS는 395.2초, 차수 6일 때 라크스-벤드로프는 972.2초, 크랭크-니콜슨은 906.6초이다.
• 기호 패키지는 표준 PC에서 Mathematica 5.2를 사용하여 테스트된 모든 스킴의 대칭을 성공적으로 계산하였으며, 이는 실현 가능성과 확장성의 가능성을 보여준다.