[논문 리뷰] On the stability of robust dynamical low-rank approximations for hyperbolic problems

Jonas Kusch, Lukas Einkemmer외 1명|arXiv (Cornell University)|2021. 01. 01.

Model Reduction and Neural Networks참고 문헌 35인용 수 2

한 줄 요약

이 논문은 초순수성과 시간 이산화 이전에 동적 저질서 근사(DLRA)를 적용하는 새로운 프로젝터 분할 통합기법을 제안하며, 고전적 CFL 조건을 복원하고 쌍곡 문제에 대한 안정성을 향상시킨다. 일阶 모멘트에 대해 뛰어난 안정성을 보이며, 운동량 수송 방정식에서 산산각 항에 대한 효율적이고 안정적인 갱신을 제시한다. 수치 실험을 통해 검증되었다.

ABSTRACT

The dynamical low-rank approximation (DLRA) is used to treat high-dimensional problems that arise in such diverse fields as kinetic transport and uncertainty quantification. Even though it is well known that certain spatial and temporal discretizations when combined with the DLRA approach can result in numerical instability, this phenomenon is poorly understood. In this paper we perform a L2 stability analysis for the corresponding nonlinear equations of motion. This reveals the source of the instability for the projector splitting integrator when first discretizing the equations and then applying the DLRA. Based on this we propose a projector splitting integrator, based on applying DLRA to the continuous system before performing the discretization, that recovers the classic CFL condition. We also show that the unconventional integrator has more favorable stability properties and explain why the projector splitting integrator performs better when approximating higher moments, while the unconventional integrator is generally superior for first order moments. Furthermore, an efficient and stable dynamical low-rank update for the scattering term in kinetic transport is proposed. Numerical experiments for kinetic transport and uncertainty quantification, which confirm the results of the stability analysis, are presented.

연구 동기 및 목표

• 쌍곡 문제의 공간 및 시간 이산화 이후에 DLRA를 적용할 때 발생하는 수치적 불안정성의 원인을 이해하기 위해.
• 공간 및 시간 이산화 이전에 DLRA를 적용함으로써 고전적 CFL 조건을 유지하는 안정화된 프로젝터 분할 통합기법을 개발하기 위해.
• 새로운 통합기법과 기존 프로젝터 분할 방법 간의 안정성과 정확도를 비교하여, 일阶 및 고계 모멘트에 대해 분석하기 위해.
• 운동량 수송 방정식의 산산각 항에 대해 효율적이고 안정적인 동적 저질서 갱신을 설계하기 위해.
• 이론적 안정성 분석 결과를 운동량 수송 및 불확실성 정량화 문제에서의 수치 실험을 통해 검증하기 위해.

제안 방법

• DLRA를 이산화 이후 적용할 때 발생하는 비선형 운동 방정식의 L2 안정성 분석을 수행하여 불안정성의 원인을 규명한다.
• 공간 및 시간 이산화 이전에 연속 시스템에 대해 DLRA를 적용하는 새로운 프로젝터 분할 통합기법을 제안한다.
• 조기 저질서 투영을 통해 안정성을 확보함으로써 고전적 CFL 조건을 복원하는 수정된 통합기법을 유도한다.
• 일阶 모멘트와 고계 모멘트에 대해 기존 통합기법과 비표준 통합기법의 안정성 특성을 분석하고 비교한다.
• 운동량 수송 방정식의 산산각 항에 대해 안정성과 효율성을 유지하는 특수화된 동적 저질서 갱신 방법을 설계한다.
• 운동량 수송 및 불확실성 정량화 문제에 대해 제안된 통합기법을 구현하고 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 DLRA를 공간 및 시간 이산화 이후에 적용할 때 프로젝터 분할 통합기법에서 발생하는 수치적 불안정성의 원인은 무엇인가?
• RQ2 공간 및 시간 이산화 이전에 저질서 근사를 적용함으로써 고전적 CFL 조건을 복원할 수 있는 DLRA 기반 통합기법을 설계할 수 있는가?
• RQ3 새로운 통합기법과 기존 프로젝터 분할 방법 간의 안정성과 정확도는 일阶 모멘트와 고계 모멘트에 대해 어떻게 비교되는가?
• RQ4 운동량 수송 방정식의 산산각 항에 대해 효율적이고 안정적인 동적 저질서 갱신을 구성할 수 있는가?
• RQ5 수치 실험은 제안된 통합기법의 이론적 안정성 예측을 확인하는가?

주요 결과

• 기존 프로젝터 분할 통합기법의 불안정성은 공간 및 시간 이산화 이후에 DLRA를 적용함으로써 안정성에 필요한 기초 구조를 파괴하기 때문이다.
• 이산화 이전에 DLRA를 적용하는 제안된 통합기법은 고전적 CFL 조건을 성공적으로 복원하여 안정성을 보장한다.
• 비표준 통합기법은 전반적으로 더 유리한 안정성 특성을 보이며, 특히 일阶 모멘트에 대해 뛰어나다.
• 기존 프로젝터 분할 통합기법은 비선형 운동 방정식의 특정 구조 덕분에 고계 모멘트 근사에 있어 더 우수한 성능을 보인다.
• 운동량 수송 방정식의 산산각 항에 대해 효율적이고 안정적인 동적 저질서 갱신 방법을 성공적으로 유도하고 수치적으로 검증하였다.
• 운동량 수송 및 불확실성 정량화 문제에서의 수치 실험은 이론적 안정성 분석을 확인하였으며, 제안된 통합기법의 우수성을 입증하였다.