Metagraph-Based Domain-Decomposed Galerkin Reduced-Order Model

기술 요약: 메타그래프 기반 도메인 분할 갈레르킨 차수 축소 모델 (MBDD-G-ROM)

문제 정의
데이터 구동형 차수 축소 모델(ROM), 특히 적절 직교 분해(POD)에 기반한 갈레르킨 차수 축소 모델(G-ROM)은 매개변수 연구를 위한 상당한 계산 가속을 제공하며, 이는 저차원 근사 공간으로 지배 방정식을 투영함으로써 이루어진다. 그러나 기저가 되는 전계 모델(FOM)이 수백만 개의 자유도(예: 고레이놀즈 수 유동)로 확장됨에 따라, ROM의 구축 및 적용은 심각한 메모리 및 계산 병목 현상에 직면한다.

기존 방식들은 다음과 같은 구체적인 한계점을 가진다:

1. 전역 POD (Global POD): 전역 저차원 공간을 구축하는 것은 밀집된 축소 계수 행렬을 생성하여, 온라인 단계(조립 및 솔루션)의 분산 메모리 병렬화를 과도한 프로세스 간 통신으로 인해 비효율적으로 만든다.
2. 도메인 분할 G-ROM (DD-G-ROM): 이러한 방법들은 국부적 현상을 처리하기 위해 국부적 근사 공간을 구축하지만, 일반적으로 국부 기저를 구축하는 데 사용된 분할이 병렬 컴퓨팅에 사용되는 분할과 일치한다고 가정한다. 이러한 경직된 결합은 유연한 부하 분산(load balancing)을 방해하고, 가용 병렬 하드웨어와 독립적으로 로컬 서브도메인의 수를 조절하는 능력을 제한한다. 또한, 대부분의 기존 DD-G-ROM 프레임워크는 축소된 시스템의 선형 솔버를 포함하여 온라인 단계 전체를 병렬화하지 못한다.

방법론
저자들은 국부 근사 공간의 구축과 병렬 실행 전략을 분리하도록 설계된 **메타그래프 기반 도메인 분할 갈레르킨 차수 축소 모델(MBDD-G-ROM)**을 제안한다. 이 방법론은 계층적 그래프 표현에 의존한다:

1. 그래프 정의:

   • 계산점 그래프 (Calculation-Point Graph): FOM 자코비안 또는 강성 행렬의 희소 패턴에 대응하는 계산점(예: 유한 요소 노드) 간의 상호작용을 나타낸다.
   • 메타그래프 (Metagraph): "국부 근사 공간 서브도메인"(POD 계산 서브도메인) 간의 연결성을 나타낸다. 이 그래프에서 메타노드(metanode)는 서브도메인에 대응하며, 메타에지(metaedge)는 국부 POD 기저 함수의 지지 영역(support) 간 중첩에 의해 유도되는 블록 희소 패턴을 인코딩한다.

2. 2단계 분할 전략:

   • 단계 1 (POD 계산 서브도메인): 계산점 그래프를 분할하여 국부 POD 기저를 독립적으로 구축할 서브도메인을 정의한다(Local POD). 이를 통해 물리 법칙이나 기하학적 구조에 맞춘 임의의 도메인 분할이 가능하다.
   • 단계 2 (병렬 계산 서브도메인): 생성된 메타그래프를 다시 분할하여 병렬 프로세스(MPI 랭크)에 할당될 서브도메인을 정의한다. 이 두 번째 단계는 POD 서브도메인과 병렬 프로세스 사이의 결합을 해제하여, 병렬 프로세스 수와 상관없이 POD 서브도메인의 수와 기하학적 구조를 독립적으로 설정할 수 있게 한다.

3. 병렬 구현:

   • 본 방법은 국부 POD 서브도메인 간의 상호작용을 나타내는 밀집 블록 행렬을 처리하기 위해 중첩 도메인 분할 기법을 확장한다.
   • 벡터 합, 내적, 행렬-벡터 곱과 같은 병렬 연산은 메시지 패싱 인터페이스(MPI)를 사용하여 구현된다. 데이터 통신은 메타그래프 인접성에 의해 결정된, 축소된 시스템 내의 0이 아닌 블록에 해당하는 벡터 성분만을 교환함으로써 최적화된다.
   • 정적 부하 분산 (Static Load Balancing): 프레임워크는 추정된 계산 비용(예: 기저 벡터 수 또는 하이퍼-축소 요소 수)을 나타내는 메타노드 가중치를 포함한다. 가중치 그래프 분할 알고리즘(METIS)을 사용하여 이 가중치 메타노드를 병렬 프로세스에 분배함으로써 부하를 균형 있게 조절한다.

4. 하이퍼-축소 (Hyper-Reduction): 효율성을 더욱 높이기 위해, 본 방법은 일부 요소에 대한 적분을 근사하여 축소 시스템의 조립 비용을 줄이는 경험적 큐버처 방법(ECM)을 통합한다.

주요 기여

• 새로운 그래프 표현: 메타그래프 개념의 도입은 임의의 도메인 분할을 수용할 수 있는 DD-G-ROM을 위한 일반적인 그래프 이론적 정식화를 제공한다.
• 분리된 분할 (Decoupled Partitioning): 2단계 분할 전략을 통해 국부 근사 공간의 서브도메인을 병렬 계산 서브도메인과 독립적으로 지정할 수 있다. 이는 병렬 프로세스와 1대 1 대응 관계를 요구하지 않으면서도 많은 수의 국부 서브도메인을 사용할 수 있게 한다.
• 전체 온라인 단계 병렬화: 본 프레임워크는 축소 시스템의 조립과 반복 솔버를 통한 축소 선형 시스템의 솔루션을 포함한 전체 온라인 단계의 분산 메모리 병렬화를 달성한다.
• 부하 분산 메커니즘: 메타노드에 계산 비용 추정치(예: 기저 벡터 수 또는 하이퍼-축소 요소 수)를 기준으로 가중치를 할당함으로써 자연스러운 정적 부하 분산 메커니즘을 제공한다.

수치 결과
본 방법은 비정상 확산 방정식과 3D 실린더 주위의 비압축성 나비에-스토크스 방정식이라는 두 가지 테스트 케이스를 통해 검증되었다.

• 정확도: 확산 방정식 케이스에서 POD 계산 서브도메인(및 그에 따른 총 기저 벡터 수)의 수를 늘릴수록 근사 정확도가 향상되었으며, 이는 국부 POD 구현의 타당성을 확인시켜 준다.
• 병렬 효율성: 강력 스케일링(Strong-scaling) 테스트 결과, 본 방법은 온라인 단계에서 높은 병렬 효율성을 유지함을 보여주었다. ROM 효율성(FEM 대비 가속도)은 특히 FOM 자유도가 큰 문제에서 병렬 프로세스 수가 증가함에 따라 안정적이거나 증가하는 경향을 보였다.
• 확장성: 축소된 선형 솔버는 양호한 강력 스케일링 동작을 보였으나, 프로세스 수가 매우 많아져 프로세스당 축소 차원이 통신 오버헤드를 상쇄할 만큼 작아질 경우 확장성이 저하되었다.
• 부하 분산: 정적 부하 분산 테스트에서 하이퍼-축소 요소 수($n^{(i)}_{HROM}$)를 기준으로 메타노드 가중치를 할당했을 때, 가중치 없는 분할에 비해 계산 효율이 최대 약 15% 개선되었으며, 이는 가중치 그래프 접근 방식의 유용성을 입증한다.

의의 및 주장
본 논문은 MBDD-G-ROM이 대규모 DD-G-ROM 분석의 메모리 제약을 완화하기 위한 일반적인 프레임워크를 제공한다고 주장한다. POD 구축 서브도메인과 병렬 실행 서브도메인을 분리함으로써, 이 방법은 기존 방식보다 고성능 컴퓨팅 자원을 더 유연하게 사용할 수 있게 한다. 저자들은 본 프레임워크가 전역 POD의 밀집된 축소 행렬 문제와 기존 DD-G-ROM의 경직된 결합 문제를 성공적으로 해결하여, 오프라인(기저 구축) 및 온라인(시뮬레이션) 단계 모두에서 높은 병렬 효율성을 달성했다고 단언한다. 연구의 결론은 본 방법이 솔루션의 정확도를 유지하면서도 복잡하고 대규모인 문제에 대해 확장 가능한 병렬 분석을 가능하게 한다는 것이다.