[논문 리뷰] A new symmetric linearly implicit exponential integrator preserving polynomial invariants or Lyapunov functions for conservative or dissipative systems

Lu Li|arXiv (Cornell University)|2021. 04. 25.

Numerical methods for differential equations참고 문헌 30인용 수 9

한 줄 요약

이 논문은 보존성 및 소산성 강성 시스템에 대해 다항형 불변량 또는 리아푸노프 함수를 유지하는 대칭적이고 선형적으로 암시적인 지수 적분기를 제안한다. 극성화된 이산 기울기와 지수 적분기를 조합함으로써, 전적으로 암시적인 방법에 비해 상당히 감소된 계산 비용으로 에너지 보존과 뛰어난 장기 안정성을 달성한다.

ABSTRACT

We present a new linearly implicit exponential integrator that preserves the polynomial first integrals or Lyapunov functions for the conservative and dissipative stiff equations, respectively. The method is tested by both oscillated ordinary differential equations and partial differential equations, e.g., an averaged system in wind-induced oscillation, the Fermi-Pasta-Ulam systems, and the polynomial pendulum oscillators. The numerical simulations confirm the conservative properties of the proposed method and demonstrate its good behavior in superior running speed when compared with fully implicit schemes for long-time simulations.

연구 동기 및 목표

• 강성 있는 보존성 및 소산성 시스템에 대해 다항형 에너지 함수를 갖는 구조 보존 적분기를 개발한다.
• 전적으로 암시적인 방법의 한계를 해결하기 위해 계산 비용을 감소시키기 위해 선형적으로 암시적인 방법을 도입한다.
• 기하학적 적분 기법을 통해 대칭성과 장기 안정성을 확보한다.
• 극성화된 이산 기울기와 지수 시간 적분을 사용하여 에너지 또는 리아푸노프 함수를 정확히 유지한다.
• 전적으로 암시적 및 명시적 방법과 비교해 장기 시뮬레이션에서 뛰어난 성능을 보여준다.

제안 방법

• 변수 상수 공식과 지수 시간 스텝을 사용하여 대칭적 선형적으로 암시적인 지수 적분기를 구성한다.
• 비선형성을 다중 시간 단계에 걸쳐 분할하기 위해 극성화된 이산 기울기를 적용하여 선형적으로 암시적인 방법을 가능하게 한다.
• 다항형 위치 에너지 함수의 제곱 극성화(예: 세차, 사차)를 사용하여 에너지 또는 리아푸노프 함수를 유지하는 이산 기울기를 정의한다.
• 식 yn+2 = exp(2hJM)yn + 2hφ(2hJM)J∇¯U(yn, yn+1, yn+2)로 방법을 구성하며, 여기서 ∇¯U는 극성화된 이산 기울기이다.
• 다중점 극성화와 이산 기울기 조건을 사용하여 다항형 잠재 에너지 함수의 고차항으로 일반화한다.
• 일致한 평균화와 기울기 재구성 기반으로 2h 시간 단계에 걸쳐 방법을 구성함으로써 대칭성을 확보한다.

실험 결과

연구 질문

• RQ1 선형적으로 암시적인 지수 적분기는 강성 있는 보존성 및 소산성 시스템에서 다항형 불변량 또는 리아푸노프 함수를 유지할 수 있는가?
• RQ2 장기 시뮬레이션에서 전적으로 암시적 및 명시적 방법과 비교해 제안된 방법은 정확도와 계산 비용 측면에서 어떻게 비교되는가?
• RQ3 적분기의 대칭성은 선형 오차 증가 및 불변량의 근사 보존과 같은 향상된 장기 행동을 이끌어내는가?
• RQ4 제한된 자유 에너지나 추가 스칼라 변수 없이도 고차 다항형 잠재 에너지 함수로 일반화할 수 있는가?
• RQ5 다항형 에너지 함수를 갖는 시스템에 대해 이 방법의 정밀도 순서와 수렴 행동은 무엇인가?

주요 결과

• 제안된 LIEEP 방법은 보존성 시스템, 예를 들어 α-Fermi–Pasta–Ulam 시스템과 다항형 진자 진동자에서 극성화된 에너지를 정확히 유지한다.
• 소산성 시스템에서는 리아푸노프 함수의 단조 감소를 유지하여 물리적 소산 행동와 일치한다.
• 수렴도를 통해 시간 순서가 2임을 확인하였다. h = 1/2^i (i = 1에서 5까지)를 사용한 수렴 플롯을 통해 검증하였다.
• 수치 결과는 LIEEP 방법이 전적으로 암시적인 EAVF 방법보다 상당히 낮은 계산 비용을 기록함을 보여주며, 특히 장기 시뮬레이션에서 뚜렷한 이점이 있다.
• 적절한 극성화된 에너지 공식을 사용할 경우 일부 시스템에서 초수렴 행동를 보인다.
• h = 0.3일 때 다항형 진자 진동자의 수치 해는 [0, 1000] 동안 진짜 운동을 잘 근사하며, 원통형 위상공간의 구조를 유지한다.