Quantum lattice Boltzmann method via density-matrix encoding for fluid simulation with wall boundary conditions

기술 요약: 벽 경계 조건을 포함한 유체 시뮬레이션을 위한 밀도 행렬 인코딩 기반 양자 격자 볼츠만 방법(QLBM)

문제 정의
양자 격자 볼츠만 방법(QLBM)은 양자 병렬성과 지수적 메모리 효율성을 활용하여 효율적인 유체 역학 시뮬레이션을 제공하는 유망한 패러다임이지만, 그 실제 적용은 두 가지 주요 과제로 인해 제약을 받아왔다. 첫째는 비선형 충돌 항과 유니터리(unitary) 양자 연산 간의 부적합성이며, 둘째는 더 결정적으로, 현실적인 경계 조건을 부과하는 것의 어려움이다. Wang 등[37]에 의한 앙상블 기반 접근 방식과 같은 기존의 QLBM 프레임워크는 주로 주기적 경계 조건에 국한되어 있었다. 양자 컴퓨팅의 본질적인 유니터리 및 가역적 특성은 나비에-스토크스 방정식의 비선형적이고 비가역적인 성질, 그리고 임의의 기하학적 구조에 따른 고체 벽 경계(예: 노 슬립(no-slip) 조건) 및 입/출구 속도 프로파일의 특정 요구 사항과 충돌한다.

방법론
저자들은 세 가지 핵심적인 방법론적 발전을 통해 이러한 한계를 해결하는 일반화된 QLBM 프레임워크를 제안한다:

1. 밀도 행렬 인코딩(Density-Matrix Encoding): 순수 상태(pure-state) 인코딩을 넘어, 저자들은 격자 가스 앙상블의 확률 분포를 밀도 행렬($\varrho$)로 인코딩한다. 이 정식화는 크라우스 연산자(Krauss operators)를 통해 비유니터리 연산을 표현할 수 있게 한다. 완전 양의 보존적(Completely Positive Trace-Preserving, CP-TP) 맵 조건을 만족함으로써, 이 알고리즘은 클래식 LBM과 수학적으로 동일하면서도 스티른-체커 디라데이션 정리(Stinespring dilation theorem)를 통해 양자 하드웨어와 호환되는 확률적 과정(충돌 및 경계 반사 등)을 구현할 수 있다.
2. 벽 경계 조건을 위한 성분 교환(Component Exchange for Wall Boundaries): 스트리밍 연산자를 수정하지 않고 임의의 기하학적 구조에 대해 노 슬립 벽 조건을 강제하기 위해, 저자들은 스트리밍 단계 직전에 "성분 교환" 단계를 도입한다.
   • 격자점은 경계, 인접, 경계-내부, 유체-내부 범주로 분류된다.
   • 벽 근처의 인접 점들에 대해, 알고리즘은 입자 상태를 여러 브랜치(branch)로 분리한다. 유체를 향해 스트리밍하는 브랜치는 정상적으로 진행하는 반면, 벽을 향해 스트리밍하는 브랜치는 역방향 속도 구성을 가지고 동일한 격자점으로 다시 반사된다.
   • 이는 운송이 일어나기 전에 속도 성분을 교환하는 사전 계산된 치환(SWAP 연산)을 통해 수행되며, 자연스럽게 하프웨이 바운스 백(halfway bounce-back) 규칙을 만족시킨다.
3. 보조 레지스터를 통한 입/출구 조건(Inlet/Outlet Conditions via Ancilla Registers): 제어된 SWAP 연산을 사용하여 입/출구 속도 프로파일을 부과한다. 보조 레지스터(ancilla register)를 목표 속도 분포 상태로 준비한다. 제어된 SWAP은 대상 격자점의 속도 레지스터를 보조 레지스터와 교환하여, 국소 분포를 규정된 프로파일로 효과적으로 덮어쓴다. 그 후 보조 레지스터는 트레이스 아웃(trace out)된다.

전체 알고리즘은 (1) 성분 교환, (2) 스트리밍, (3) 충돌, (4) H-단계(속도 분포의 독립성을 복원하고 열역학적 평형을 유지하기 위함), (5) 입/출구 강제 실행의 순서로 반복된다.

주요 기여

• 정적 벽을 위한 일반 경계 스킴: 본 논문은 이전의 앙상블 기반 QLBM이 가졌던 주기적 경계 제한을 극복하고, 임의의 기하학적 구조를 가진 정적 고체 벽 경계와 입/출구 조건을 처리할 수 있는 QLBM 알고리즘을 제시한다. 저자들은 현재 프레임워크에서 움직이는 경계 문제는 완전히 다루어지지 않았으며 아직 구현되지 않은 특별한 처리가 필요함을 명시적으로 언급한다.
• 비유니터리 구현: 밀도 행렬 인코딩과 크라우스 연산자가 양자 장치에서 필요한 비유니터리 유체 역학 연산(충돌 및 경계 반사)을 어떻게 구현할 수 있는지 보여준다.
• 알고리즘 설계: "성분 교환" 단계의 구체적인 설계는 표준 스트리밍 회로를 수정하지 않고도 유지할 수 있게 하여, 복잡한 기하학적 구조에 대한 회로 구조를 단순화한다.
• 복잡도 분석: 저자들은 경계 처리를 포함하더라도 스트리밍 단계가 지배적인 항($O(2^q n^2)$)이므로 원래 QLBM의 점근적 게이트 복잡도가 변하지 않는다는 자원 분석을 제공한다.

시뮬레이션 결과
알고리즘은 클래식 컴퓨터 상의 양자 에뮬레이터를 사용하여 몇 가지 2D 벤치마크 흐름에 대해 검증되었다:

• 감쇠 푸아죄유 흐름(Decaying Poiseuille Flow): 시뮬레이션은 시간 경과에 따른 삼각 함수 형태의 속도 프로파일을 정확하게 복구하였고 벽에서의 노 슬립 조건을 강제하였으며, 격자 세분화에 따른 수렴성을 보여주었다.
• 후방 단계 흐름(Backward-Facing Step Flow): 이 방법은 흐름 박리와 재순환 영역을 성공적으로 포착하였다. 무차원 재순환 길이($L_r/h$)는 OpenFOAM 참조 솔루션(1.375 대 1.1563)과 밀접한 일치를 보였으며, 레이놀즈 수가 증가함에 따라 상대 오차가 감소하였다.
• 실린더를 지나는 흐름(Flow Past a Cylinder): 알고리즘은 곡선형 경계를 효과적으로 처리하였으며, 실린더 뒤쪽의 저속 영역을 참조 데이터 대비 12.53%에서 17.34%의 상대 오차로 포착하였다.
• 복잡한 기하학적 구조 ("PKU" 장애물): 방법론은 "PKU"라는 글자 모양의 복잡한 형상을 지나는 흐름에 대해 테스트되었으며, 틈새와 후류에서의 유선 및 와류 형성을 성공적으로 포착하여 비정형 기하학적 구조로의 확장성을 입증하였다.

의의 및 주장
본 논문은 이 연구가 실제 전산유체역학(CFD)을 위한 실질적인 QLBM 배포를 향한 중요한 단계라고 주장한다. 정적 벽에 대한 경계 조건 병목 현상을 해결함으로써, 제안된 방법은 양자 유체 시뮬레이션의 적용 범위를 이상적인 주기적 영역에서 고체 벽 및 복잡한 기하학적 구조를 포함하는 문제로 확장한다. 저자들은 제안된 방식이 물리적 경계의 비유니터리 요구 사항을 처리할 수 있는 유연한 회로 설계를 도입하면서도, 양자 컴퓨팅의 이론적 이점(병렬성 및 메모리 효율성)을 유지한다고 강조한다.

그러나 저자들은 즉각적인 실용적 배포에 대해서는 신중한 태도를 보인다. 그들은 H-단계가 상당한 복잡성과 회로 깊이를 초래하며, 이것이 현재로서는 중규모 문제에 대한 양자 가속도를 저해할 수 있음을 인정한다. 또한, 현재 프레임워크는 정적 경계만을 다루며, 움직이거나 변형되는 경계에 대한 연구는 향후 과제로 남겨두고 있다. 이 작업은 완전히 실현된 하드웨어 데모라기보다는 기초적인 알고리즘 확장으로서 제시되었다.