[논문 리뷰] Reminiscences about numerical schemes

С. К. Годунов|ArXiv.org|2008. 10. 03.

Computational Fluid Dynamics and Aerodynamics참고 문헌 11인용 수 5

한 줄 요약

이 논문은 1953년에서 1969년 사이에 초순응력 시스템의 보존법칙을 위한 고두노프 스킴의 발전에 대해 S. K. 고두노프가 기록한 개인적 회고이다. 인공 점성에 기반한 1차 스킴에서 출발하여 고차수 방법, 격자 적응, 충격 불연속성의 모델링에서 수치적 점성이 차지하는 핵심적 역할까지의 발전 과정을 기술하며, 최종적으로 엔트로피 조건을 만족하는 일반화된 해의 형식화와 스킴에 따라 달라지는 수치적 소산 효과의 인식으로 이어진다.

ABSTRACT

This preprint appeared firstly in Russian in 1997. Some truncated versions of this preprint were published in English and French, here a fully translated version is presented. The translation in English was done by O. V. Feodoritova and V. Deledicque to whom I express my gratitude. In the present paper I will describe how the first variant of the Godunov's scheme has been elaborated in 1953-1954 and tell about all modifications realized by myself (until 1969) and the group of scientists from the Institute of Applied Mathematics in Moscow (which has become the M.V.Keldysh Institute of Applied Mathematics).

연구 동기 및 목표

• 모스크바 응용수학연구소에서의 초기 년도 동안 경험한 바를 바탕으로 고두노프의 초순응력 보존법칙 수치 스킴의 역사적 발전을 재구성하는 것.
• 충격파 근처에서 2차 스킴의 안정화에 직면한 과제와 그로 인해 1차 정확도를 가진 보존적 방법으로의 전환을 설명하는 것.
• 수치적 점성이 물리적 소산을 어떻게 근사하는지와 이와의 연관성, 이로 인해 같은 문제에 대해 서로 다른 스킴이 다른 해를 도출할 수 있음을 밝혀내는 것.
• 수치 스킴, 일반화된 해, 열역학적 엔트로피 조건 사이의 연결 고리를 설정하는 것. 특히 이동파 해와 소점성 극한을 통해 설명하는 것.
• 물리학자들과 수학자들 간의 다학제적 논의가 계산 유체역학의 이론적 및 계산적 기초를 어떻게 다듬는지 강조하는 것.

제안 방법

• A. I. 주코프가 사용한 기존 알고리즘을 수정하여, 잘라내기 오차에 기인한 수치적 점성으로 충격 프로파일을 안정화시키는 1차 유한차분을 사용하였다.
• 고두노프는 오일러 기준의 보존적 유한부피 형식을 도입하여, 셀 경계에서의 유량을 정확한 리만 해법을 사용해 계산하였다. 이후 이는 고두노프 스킴로 체계화되었다.
• 수치적 점성을 분석하기 위해 첫 번째 미분근사 개념을 활용하였으며, 잘라내기 오차와 차분 스킴에서의 인위적 소산 간의 연관성을 설정하였다.
• 작은 점성 하에서의 해의 거동을 분석하기 위해 형태가 $ q_i = q_i((x - \rho t)/\varepsilon) $ 인 이동파 해를 사용하였으며, 이로 인해 소산 항의 형태에 따라 해가 의존함을 밝혀냈다.
• A. V. 자브로딘과 G. P. 프로코포프가 개발한 이동 격자 및 도메인 분할 기법을 통해 복잡한 다스케일 문제와 이동 불연속성을 가진 문제의 해를 구할 수 있었다.
• 격자 생성에서 발생하는 타원형 시스템의 이론적 분석은 수치 선형대수학의 발전과 함께, 행렬 펜슬에 대한 스펙트럼 이분 기준의 제안으로 이어졌다.

실험 결과

연구 질문

• RQ1 왜 2차 스킴은 충격파 근처에서 안정된 해를 도출하지 못했으며, 이러한 불안정성의 배경이 되는 물리적 또는 수학적 원리는 무엇이었는가?
• RQ2 유한차분 스킴에서 유래하는 수치적 점성은 물리적 소산과 어떻게 관련되어 있으며, 이를 초순응력 보존법칙의 일반화된(엔트로피 조건을 만족하는) 해를 근사하는 데 사용할 수 있는가?
• RQ3 다른 수치적 스킴이 동일한 초순응력 시스템에 대해 서로 다른 해를 도출할 수 있는 정도는 어느 정도이며, 점성이 0으로 수렴할 때의 극한 해의 유일성은 무엇에 의해 결정되는가?
• RQ4 불연속 흐름과 엔트로피 조건의 맥락에서 일반화된 해의 구조를 어떻게 엄밀하게 정의하고 정당화할 수 있는가?
• RQ5 실험 물리학자들로부터의 물리적 직관과 비판적 피드백은 수치 알고리즘과 그 이론적 기초를 어떻게 다듬는 데 기여하는가?

주요 결과

• 정확한 리만 해법과 보존적 유량을 기반으로 한 1차 고두노프 스킴은 충격 포착에 있어 강력하고 안정적이었으며, 물리적으로 비합리적인 진동을 유발하는 2차 스킴보다 뛰어난 성능을 보였다.
• 유한차분 스킴에서 유래하는 수치적 점성은 작은 인위적 소산을 도입하는 것과 동일시되었지만, 그로 인한 해는 점성이 0으로 수렴하는 것 외에도 소산 항의 형태에 따라 달라졌다.
• B. F. 디아첸코가 제시한 예시는 동일한 초순응력 시스템에 대해 서로 다른 수치적 스킴이 다른 해를 도출할 수 있음을 보여주었으며, 이는 해의 스킴 불변성 가정을 무너뜨렸다.
• 이동파 해의 사용은 오직 특정 형태의 점성만이 엔트로피 조건을 만족하는 불연속 해를 도출할 수 있음을 드러내었으며, 이로 인해 수치 스킴과 열역학 원리 사이의 연결 고리를 확립하였다.
• 이동 격자 및 도메인 분할의 개발은 이동하는 파면과 다중 불연속성을 가진 복잡한 흐름의 시뮬레이션을 가능하게 하였으며, 이로 인해 방법의 적용 범위가 확장되었다.
• 격자 생성에 관한 이론적 연구는 스펙트럼 이분 기준의 제안과 특이 벡터 계산을 위한 고정밀 알고리즘 개발으로 이어졌으며, 이는 수치 선형대수학 전반에 걸쳐 광범위한 영향을 미쳤다.