[논문 리뷰] Superconvergence of Galerkin variational integrators

Sina Ober‐Blöbaum, Mats Vermeeren|arXiv (Cornell University)|2021. 04. 25.

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

한 줄 요약

이 논문은 다항식의 차수 s 이하인 갈레르킨 변분 적분법과 충분히 정확한 구적법을 사용할 경우 초수렴(superconvergence)이 달성되며, 수렴 차수는 2s임을 증명한다. 증명은 변분법과 향상된 오차 추정을 활용하여 이전 결과를 일반 라그랑주 시스템으로 확장하며, 외력이 작용하는 시스템에 대한 일반화 가능성을 시사한다.

ABSTRACT

We study the order of convergence of Galerkin variational integrators for ordinary differential equations. Galerkin variational integrators approximate a variational (Lagrangian) problem by restricting the space of curves to the set of polynomials of degree at most $s$ and approximating the action integral using a quadrature rule. We show that, if the quadrature rule is sufficiently accurate, the order of the integrators thus obtained is $2s$.

연구 동기 및 목표

• 일반 미분방정식에 대한 갈레르킨 변분 적분법의 초수렴 성질을 확립하기.
• 충분히 정확한 구적법 조건 하에서 기존의 수렴 차수 min(s, u)를 2s로 확장하기.
• 변분법과 정밀한 오차 추정을 활용한 초수렴의 일반적 증명 제공하기.
• 확장된 라그랑주 형식을 통해 외력이 작용하는 시스템으로의 초수렴 확장 탐색하기.
• 기계 시스템에서 수정된 라그랑주 함수와 강건성 조건(coercivity conditions)에 대한 향후 연구의 기초 마련하기.

제안 방법

• 궤적을 차수 ≤ s인 다항식으로 제한하고 행동 적분을 구적법으로 근사하여 갈레르킨 변분 적분법을 구성한다.
• 행동 함수수의 오차를 추정하기 위해 변분법을 적용하며, Hall과 Leok(2015)의 이전 결과보다 개선된 오차 추정을 이룩한다.
• 라그랑주 함수의 비퇴도성(non-degeneracy) 조건을 활용하여 오일러-라그랑즈 방정식이 이阶 미분방정식임을 보장하고, 레지외르 변환의 가역성을 확보한다.
• 외력이 작용하는 시스템에 대해 상태 공간을 두 배로 늘리고 보조 변수를 도입하여 하밀턴 원리의 복원을 위한 수정된 라그랑주 함수 접근법을 사용한다.
• 조건 Q = q를 도입하여 원래의 외력 시스템의 역학을 복원하고, 외력이 작용하는 갈레르킨 적분법과 확장된 시스템의 갈레르킨 적분법 간의 동치성을 보여준다.
• 임계 곡선이 행동을 최소화한다고 가정하여 연속 해와 이산 최소화자 사이의 차이를 유계화한다.

실험 결과

연구 질문

• RQ1 구적법이 충분히 정확할 경우 갈레르킨 변분 적분법의 수렴 차수는 표준적인 min(s, u) 기준을 초월하는가?
• RQ2 변분 오차 분석과 향상된 오차 추정을 통해 초수렴 차수 2s에 대한 일반적 증명을 확립할 수 있는가?
• RQ3 외력 존재 조건 하에서도 초수렴이 유지되는가? 만약 그렇다면 어떤 조건에서 성립하는가?
• RQ4 위치에 따라 변화하는 질량 행렬을 갖는 기계 시스템의 라그랑주 함수에 대해, 임계 곡선을 최소화하기 위한 강건성 조건을 검증할 수 있는가?
• RQ5 확장된 라그랑주 형식을 통해 외력이 작용하는 시스템에 대한 초수렴을 엄밀하게 증명할 수 있는가?

주요 결과

• 구적법이 충분히 정확할 경우 갈레르킨 변분 적분법의 수렴 차수는 2s로 초수렴이 성립함을 입증함.
• 변분법 기법을 활용하여 행동 함수수 오차를 더 정밀하게 추정함으로써 이전 오차 추정의 한계를 초월함.
• 양의 정부호 질량 행렬을 갖는 기계 시스템의 라그랑주 함수에 대해서는, 임계 곡선을 최소화하기 위한 강건성 조건이 충족됨을 보여주며, 이는 주요 정리의 가정 조건을 만족함을 의미함.
• 수치적 증거와 저차수 예제(예: 중점 법칙)는 초수렴이 외력 시스템에도 성립할 것임을 지지하지만, 완전한 증명은 아직 열려 있음.
• 외력 시스템은 확장된 라그랑주 형식을 통해 재구성되어 하밀턴 원리를 복원할 수 있으며, 이로 인해 확장된 시스템에 대해 초수렴 결과를 적용할 수 있음.
• 현재 증명에서 외력 시스템에 대한 약간의 격차를 발견함: 식 (15)의 확장된 라그랑주 함수는 일반적으로 최소화자 가정을 만족하지 않음. 이는 향후 연구에서 이 가정 없이도 해의 편차 ˆq − ˜q가 작다는 것을 보장할 필요가 있음을 시사함.