[논문 리뷰] Numerical solution of the incompressible Navier-Stokes equations for chemical mixers via quantum-inspired Tensor Train Finite Element Method

Egor Kornev, Sergey Dolgov외 4명|arXiv (Cornell University)|2023. 05. 18.

Computational Physics and Python Applications인용 수 4

한 줄 요약

이 논문은 복잡한 T-믹서 기하구조에서 비압축성 나비에-스토크스 방정식을 해결하기 위한 양자 기반 텐서 트레이스 유한요소법(TetraFEM)을 제안한다. 유한요소 행렬을 양자화된 텐서 트레이스(QTT) 형식으로 표현함으로써, 기존 FEM 대비 지수적 메모리 압축과 속도 향상을 달성하며, TT-반올림 허용오차를 조절하여 정확도를 제어할 수 있다.

ABSTRACT

The solution of computational fluid dynamics problems is one of the most computationally hard tasks, especially in the case of complex geometries and turbulent flow regimes. We propose to use Tensor Train (TT) methods, which possess logarithmic complexity in problem size and have great similarities with quantum algorithms in the structure of data representation. We develop the Tensor train Finite Element Method -- TetraFEM -- and the explicit numerical scheme for the solution of the incompressible Navier-Stokes equation via Tensor Trains. We test this approach on the simulation of liquids mixing in a T-shape mixer, which, to our knowledge, was done for the first time using tensor methods in such non-trivial geometries. As expected, we achieve exponential compression in memory of all FEM matrices and demonstrate an exponential speed-up compared to the conventional FEM implementation on dense meshes. In addition, we discuss the possibility of extending this method to a quantum computer to solve more complex problems. This paper is based on work we conducted for Evonik Industries AG.

연구 동기 및 목표

• 복잡한 기하구조(예: T-믹서)에서 기존 FEM을 사용한 유체 흐름 시뮬레이션의 높은 계산 비용을 해결한다.
• 순시 유동 시뮬레이션에서 조밀한 메esh 이산화로 인한 메모리 및 실행 시간 병목 현상을 극복한다.
• 유한요소 행렬의 지수적 압축과 더 빠른 해법 시간을 가능하게 하는 텐서 기반 수치적 방법을 개발한다.
• 고계수 해를 처리하기 위해 이론적으로 양자 컴퓨팅으로의 확장을 탐색한다.
• 실험적 검증을 통해 비틀린 T-믹서 기준 테스트 케이스에서 방법의 정확성과 효율성을 입증한다.

제안 방법

• T-믹서 도메인을 곡선형 사각형 요소로 이산화하고, 국소 행렬 조립을 위해 기준 요소로 매핑한다.
• 유한요소 방법을 사용하여 국소 강성, 질량, 일阶 도함수 행렬을 구성하고, 이를 전역 행렬로 조립한다.
• 모든 FEM 행렬과 해 벡터를 양자화된 텐서 트레이스(QTT) 형식으로 표현하여 낮은 질서 구조를 활용한다.
• TT-형식 연산을 사용하여 비압축성 나비에-스토크스 방정식의 명시적 시간 스텝 스킴을 구현한다.
• TT-반올림 및 TT-크로스 근사법을 적용하여 행렬을 압축하고 해의 정확도를 제어한다.
• TT 형식의 행렬-벡터 곱(MatVec) 및 선형 시스템 해법기를 사용하여 시간 진전 해를 효율적으로 계산한다.

실험 결과

연구 질문

• RQ1 텐서 트레이스 방법은 T-믹서와 같은 복잡하고 비직사각형 유동 도메인에 효과적으로 적용될 수 있는가?
• RQ2 나비에-스토크스 시뮬레이션에서 FEM 행렬의 TT-형식 표현이 얼마나 메모리 사용을 줄이고 해법 시간을 가속화하는가?
• RQ3 TT-반올림 허용오차는 정확도와 계산 속도 사이의 상호 보완적 관계에 어떻게 영향을 미치는가?
• RQ4 TT 기반 수치적 방법은 양자 컴퓨터에서 효율적으로 실행될 수 있는가?
• RQ5 기존 FEM 및 실험 데이터와 비교했을 때, 이 방법은 물리적 정확성을 유지하는가?

주요 결과

• 제안된 TetraFEM 방법은 모든 FEM 행렬에 대해 지수적 압축을 달성하여 기존 FEM 대비 메모리 사용을 크게 줄였다.
• 특히 고해상도 시뮬레이션에서 조밀한 메쉬에서 기존 FEM 대비 계산 시간에 지수적 속도 향상이 관찰되었다.
• TetraFEM과 기존 FEM 해 사이의 상대 오차는 제어 가능하며, 더 엄격한 TT-반올림 허용오차로 갈수록 감소한다. 32×32 메쉬에서 ε = 1e-4일 경우 오차는 1% 미만이었다.
• 유동 복잡성이 증가함에 따라 압력 벡터의 TT-랭크가 시간이 지남에 따라 증가하여 해의 상관관계가 증가함을 나타낸다.
• 이 방법은 T-믹서에서 물리적 흐름 거동을 성공적으로 재현하였으며, 실험 데이터와 일치하는 해를 도출하였다.
• 이 프레임워크는 양자 컴퓨팅으로 확장 가능하며, 양자 병렬성과 효율적인 상태 준비 덕분에 고계수 TT 해에 대한 연산이 유리하게 스케일링된다.