[논문 리뷰] Quantifying the computability of the Lorenz system

Benjamin Kehlet, Anders Logg|arXiv (Cornell University)|2013. 06. 12.

Cellular Automata and Applications참고 문헌 18인용 수 7

한 줄 요약

이 논문은 반올림 오차가 시간이 지남에 따라 수치 정확도를 어떻게 제한하는지 분석하여 로렌츠 시스템의 계산 가능성을 정량화한다. 당접근 기반 후행 오차 분석을 사용하여, 계산 가능성은 기계 정밀도에 비례하여 선형적으로 증가하며, 해는 T ≈ 2.5 × n_mach까지 계산 가능하다는 모델을 유도한다. 여기서 n_mach는 유의 숫자의 개수이다. 예를 들어, 420자리 정밀도를 사용할 경우 T ≈ 1050까지 계산 가능하다.

ABSTRACT

It is well known that the computation of accurate trajectories of the Lorenz system is a difficult problem. Computed solutions are very sensitive to the discretization error determined by the time step size and polynomial order of the method, as well as round-off errors. In this work, we show how round-off errors limit the computability of the Lorenz system and quantify exactly the length of intervals over which solutions can be computed, expressed in terms of the floating point precision. Using adjoint-based \emph{a posteriori} error analysis techniques, we estimate the stability of computations with respect to initial data, discretization, and round-off errors, respectively. The analysis is verified by computing an accurate solution on the time interval $[0, 1000]$ using a very high order (order 200) finite element method and very high floating point precision 400 digits).

연구 동기 및 목표

• 주어진 기계 정밀도에서 로렌츠 시스템을 얼마나 오래 계산할 수 있는지에 대한 근본적인 질문을 해결하기 위해.
• 혼돈 시스템인 로렌츠 애트랙터와 같은 시스템의 수치적 해에서 반올림 오차가 차지하는 제한적 역할을 정량화하기 위해.
• 기존의 사전 오차 추정치가 과도하게 낙관적이며 진정한 계산 가능성의 한계를 포착하지 못한다는 것을 보여주기 위해.
• 초고정밀도(400자리 산술)를 사용한 고차수(q=200) 유한요소 방법을 통해 이론적 모델을 검증하기 위해.

제안 방법

• 초기 조건, 이산화, 반올림 오차에 대한 민감도를 추정하기 위해 당접근 기반 후행 오차 분석을 사용한다.
• 복합 오차 모델을 유도한다: E ≈ [0.5‖U(0)−u(0)‖ + 0.001Δt^{2q} + (0.002 + 0.0005q)Δt^{-1/2}ε_mach] × 10^{0.388T}로, 모든 주요 오차 원인을 포함한다.
• 다항식 차수 q=200인 고차수 연속 갈레르킨(cG(q)) 유한요소 방법을 사용하여 [0,1000]에서 기준 해를 계산한다.
• 반올림 오차를 최소화하고 이론적 오차 모델을 검증하기 위해 400자리 부동소수점 정밀도를 사용한다.
• 수치 실험에서 상수 C₂^{[q]}와 C₃^{[q]}를 경험적으로 적합하여 오차 모델을 校정한다.
• 이산화 오차와 반올림 오차 기여도를 균형 잡기 위해 최적의 시간 간격 Δt ≈ ε_mach^{1/(2q+0.5)}를 유도한다.

실험 결과

연구 질문

• RQ1 주어진 기계 정밀도에서 로렌츠 시스템을 얼마나 오래 정확하게 계산할 수 있는가?
• RQ2 반올림 오차는 이산화 오차에 비해 혼돈 시스템인 로렌츠 시스템의 계산 가능성을 어떻게 제한하는가?
• RQ3 후행 오차 분석은 기존의 사전 추정치보다 더 정확하고 현실적인 계산 가능성 추정치를 제공할 수 있는가?
• RQ4 기계 정밀도(유의 숫자의 개수)가 증가함에 따라 계산 가능성은 어떻게 변화하는가?

주요 결과

• 로렌츠 시스템의 계산 가능성은 기계 정밀도에 비례하여 선형적으로 증가한다: T ≈ 2.5 × n_mach, 여기서 n_mach는 유의 숫자의 개수이다.
• 16자리 정밀도(표준 더블 정밀도)를 사용할 경우 계산 가능성은 T ≈ 40으로 제한되며, 일반적으로 가정되는 T ≈ 1.1의 사전 추정치보다 훨씬 크다.
• 420자리 정밀도를 사용할 경우, 이 방법은 T = 1000까지 정확한 해를 계산할 수 있으며, 이는 이론적 예측인 T ≈ 1050을 확인한다.
• 오차 모델 E ≈ 0.002 × ε_mach × 10^{0.4T}는 오차 성장을 정확히 예측하며, ε_mach = 10^{-16}일 경우 T ≈ 50에서 E ≈ 0.001이 된다.
• 최적의 시간 간격은 Δt ≈ ε_mach^{1/(2q+0.5)}이며, q=100이고 ε_mach=10^{-420}일 경우 Δt ≈ 0.008로 실제 사용된 Δt=0.0037에 가깝다.
• T=1000일 때 안정성 인자 S_C(T) ≈ 2 × 10^{388}은 10^{0.388T}의 지수 오차 성장률을 확인하며, 오차 모델의 지수를 검증한다.