주파수 영역 변수분리 기반 영역분할을 이용한 파라메트릭 동적 시스템 차원 축소 모델링

초록

본 논문은 시간‑Fourier 변환을 통해 파라메트릭 동적 PDE를 복소수값 타원식으로 변환하고, 복소수 전용 변수분리(VS)와 비중첩 영역분할(DD‑VS) 기법을 결합한 차원 축소 방법을 제안한다. 오프라인 단계에서 인터페이스와 서브문제에 대한 별도 표현을 사전 구축하고, 온라인 단계에서는 파라미터만 입력해 즉시 해를 복원한다. 공간 격자와 무관한 온라인 연산으로 높은 효율성을 확보한다.

상세 분석

본 연구는 파라메트릭 동적 시스템을 시간 영역에서 직접 해석하는 대신, 전체 시간 구간에 대해 Fourier 변환을 적용함으로써 복소수값 주파수‑파라메터 의존 타원식(2.2)으로 전환한다. 이 과정에서 원래의 비선형 연산자는 선형 복소수 연산자로 변환되며, 주파수 ω는 1차원 파라메터로 취급된다. 핵심은 복소수값 문제에 대해 기존 실수‑값 VS 방법을 그대로 적용할 수 없다는 점이다. 저자들은 복소수 해를 실수와 허수 부분으로 분리하고, 각각에 대해 동일한 VS 절차를 수행함으로써 ‘복소수‑VS’ 모델을 구축한다. 여기서 affine 분해(3.2)가 필수적인데, 파라메터 μ=(ω,ξ) 에 대한 bilinear 형 a와 linear 형 b를 μ‑의존 함수와 μ‑독립 행렬·벡터의 선형 결합 형태로 표현한다. 이러한 affine 구조는 오프라인‑온라인 분리를 가능하게 하며, 온라인 단계에서는 μ‑의 함수값만 평가하면 된다.

다음으로 도입된 DD‑VS는 비중첩 영역분할 프레임워크에 VS를 결합한다. 전체 도메인을 서브도메인으로 분할하고, 각 서브도메인의 내부 자유도와 인터페이스 자유도를 구분한다. 블록 가우시안 소거를 통해 전역 시스템을 인터페이스 Schur 보완식으로 축소하고, 이 인터페이스 문제에 대해 VS‑surrogate 모델을 구축한다. 서브도메인 내부 문제 역시 VS를 적용해 별도 surrogate를 만든다. 복소수 문제에 대해 실·허수 쌍을 동시에 다루는 방식은 기존 실수‑전용 DD‑VS와 구조적으로 동일하지만, 복소수 행렬을 실수 블록 형태로 재구성함으로써 기존 알고리즘을 그대로 활용한다.

오프라인 단계에서는 (i) 파라메터 샘플링을 통한 μ_k 선택, (ii) 각 μ_k 에 대해 residual r_k(μ) 를 계산하고, (iii) VS‑greedy 절차로 새로운 분리 항 ζ_j(μ)·c_j 를 추가한다. 이 과정은 인터페이스와 서브문제 각각에 독립적으로 수행되며, 최종적으로 N개의 분리 항을 갖는 저차원 모델이 완성된다. 온라인 단계에서는 주어진 새로운 파라메터 ξ에 대해 μ=(ω,ξ) 의 함수값을 평가하고, 사전 구축된 ζ_j(μ)·c_j 를 선형 결합해 복소수 해 ˆu(·,ω;ξ)를 즉시 얻는다. 이후 Gaussian quadrature 기반 LGL 포인트와 가중치를 이용해 역 Fourier 변환(2.4)을 수행하면 시간‑도메인 해 u(x,t;ξ)를 복원한다. 중요한 점은 온라인 연산이 전혀 FEM 격자 크기에 의존하지 않으며, 복소수 연산만으로 충분하다는 것이다.

이 방법은 기존 Reduced Basis(POD) 방식과 달리 전역 행렬을 직접 압축하지 않고, 인터페이스와 서브문제 각각을 독립적으로 압축한다는 장점이 있다. 또한, 복소수‑VS와 DD‑VS를 결합함으로써 고주파수 다중 파라메터 문제에서도 병렬 처리와 메모리 절감 효과를 기대할 수 있다. 다만, affine 분해가 성립하지 않는 비선형 파라메터 의존성이나, 강한 비정형 경계조건에 대해서는 추가적인 비선형 사전 처리나 DEIM(Discrete Empirical Interpolation Method)과 같은 기법이 필요할 수 있다. 또한, 복소수 부분을 실수‑허수 쌍으로 변환하는 과정에서 행렬 차원이 두 배가 되므로, 매우 큰 규모의 3D 문제에서는 메모리 관리가 여전히 도전 과제로 남는다.