p‑워셔슈타인 거리와 네트워크, 3차원에서 1차원으로의 수렴

초록

본 논문은 가스 네트워크를 모델링하는 메트릭 그래프 위에서 동적·정적 p‑워셔슈타인 거리를 정의하고, 정점에 질량 저장 여부에 따른 두 가지 연속성 방정식(키르히호프 조건과 저장 모델)을 제시한다. 또한, 파이프 직경을 0으로 보내 3차원 도메인에서 메트릭 그래프로 수렴할 때 정적 워셔슈타인 거리의 수렴을 c‑사이클리컬 단조성 이론을 이용해 증명한다. 마지막으로 수치 실험을 통해 이론적 결과를 검증한다.

상세 분석

논문은 먼저 최적 수송 이론의 동적 형식인 Benamou‑Brenier 접근을 메트릭 그래프에 확장한다. 그래프는 각 간선에 길이 ℓₑ를 부여한 1차원 구간들의 합집합으로 정의되며, 정점은 파이프의 연결점으로 취급한다. 이때 거리 d(x,y)는 그래프 상의 최단 경로 길이로 정의되고, 이를 기반으로 정적 p‑워셔슈타인 거리 Wₚ(μ,ν) 를 (1)식과 (2)식 형태로 기술한다. 동적 형식에서는 연속 방정식 ∂ₜρₜ+∇·(ρₜvₜ)=0 를 각 간선에 적용하고, 정점에서의 질량 흐름을 두 가지 방식으로 모델링한다. 첫 번째는 전통적인 키르히호프 조건으로, 모든 입·출 흐름의 합이 0이 되도록 하여 정점에 질량이 저장되지 않음을 가정한다. 두 번째는 정점에 별도의 질량 변수 γᵥ(t)를 도입하고, ∂ₜγᵥ(t)=∑ₑ∈E(v) ĵₑᵥ(t) 로 표현되는 저장 모델이다. 두 경우 모두 연속 방정식은 약한 형태로 정의되며, 라그랑주 승수를 이용한 이중화 과정을 통해 동적 최적화 문제를 제시한다. 특히, 저장 모델에서는 흐름과 질량 변수를 결합한 복합 변수 j=(ĵ, j) 를 도입해 전체 연속 방정식 ∂ₜμₜ+div jₜ=0 를 만족하도록 만든다.

정적 거리의 수렴 결과는 3차원 파이프 네트워크 Ω_ε (ε는 파이프 반경) 가 ε→0 일 때 메트릭 그래프 G 로 위상·측도적으로 수렴한다는 가정 하에 전개된다. 저자는 최적 수송 계획 π_ε 가 c‑사이클리컬 단조성(c‑cyclical monotonicity) 성질을 만족한다는 점을 이용해, π_ε 가 G 위의 최적 계획 π* 로 수렴함을 보인다. 핵심은 비용 함수 c_ε(x,y)=|x−y|ᵖ가 ε에 따라 그래프 거리 d_G(x,y)ᵖ 로 점근적으로 변환된다는 사실이며, 이를 통해 Wₚ(μ_ε,ν_ε)→Wₚ(μ,ν) 를 엄밀히 증명한다. 증명 과정에서는 Γ‑수렴 기법과 함께 측도 이론적 도구(밀도 함수의 약한 수렴, 연속성 등)를 활용한다.

마지막으로, 저자는 수치 실험으로 (i) 키르히호프 조건을 적용한 그래프 상의 흐름, (ii) 정점 저장 모델을 적용한 흐름, (iii) 3차원 파이프 모델을 미세 메쉬로 해석한 결과를 비교한다. 실험은 유한요소법(FEM) 기반의 연속 방정식 해석과 Sinkhorn 알고리즘을 이용한 정적 워셔슈타인 거리 계산을 포함한다. 결과는 ε가 충분히 작을 때 3차원 모델과 1차원 그래프 모델이 거의 동일한 거리와 흐름 패턴을 보이며, 정점 저장 모델이 키르히호프 모델에 비해 질량 보존을 더 유연하게 처리함을 확인한다. 전체적으로 논문은 메트릭 그래프를 통한 네트워크 최적 수송 이론을 체계화하고, 실제 물리적 파이프 네트워크와의 연계성을 수학적으로 뒷받침한다.