시간 스텝에 강인한 PolyDG 의사응력 스토크스 해법 검증

초록

본 논문은 다형격자(PolyDG) 기반 의사응력 형태의 비정상 스토크스 방정식을 시간-공간 전역적으로 이산화하고, 시간 스텝이 작아질 때 발생하는 행렬 조건 악화를 완화하기 위해 결함제거 Conjugate Gradient와 집합형 Block‑Jacobi 전처리를 제안한다. 수치 실험을 통해 두 방법 모두 Δt 에 독립적인 반복 횟수를 보이며, 향후 공간 격자 크기 h 에 대한 다중격자 결합 전략을 제시한다.

상세 분석

이 연구는 비정상 스토크스 방정식의 의사응력(pseudo‑stress) 변수를 도입함으로써 비뉴턴 흐름 및 계면 결합 문제에서 응력 텐서의 정확한 재현이 가능하도록 설계하였다. 공간 이산화는 폴리토피(다각형) 격자 위에서의 불연속 Galerkin 방법(PolyDG)을 사용하고, 시간 적분은 완전 암시적 Euler 스킴을 적용해 전형적인 대칭 양정(positive‑definite) 시스템을 만든다. 핵심 난제는 시간 스텝 Δt 가 감소하면 질량 행렬 M 의 영특이성 때문에 시스템 행렬 A* = M + Δt·A 의 조건수가 1/Δt 비례로 급격히 악화된다는 점이다. 기존 CG 혹은 표준 Block‑Jacobi 전처리는 이 스케일링을 억제하지 못해 반복 횟수가 급증한다.

첫 번째 해결책인 Deflated CG는 M 의 영공간을 정확히 포착하는 기저 V 를 구성하고, A* 에 대한 V‑정규 직교 투영 π_A*(S) 를 이용해 문제를 두 부분으로 분해한다. 이렇게 하면 가장 작은 고유값(또는 영특이값)들이 CG 과정에서 “숨겨져” 효과적으로 조건수가 감소한다. 구현상 V^T A* V 의 역을 직접 해석적으로 구하거나, Δt^2 (B1+B3) 형태의 라플라시안 연산자로 근사해 다중격자 전처리와 결합할 여지를 남긴다.

두 번째 전략인 집합형(Block‑Jacobi) 전처리는 각 요소 κ 에 대해 네 개의 텐서 성분(σ_11, σ_12, σ_21, σ_22)을 하나의 블록으로 묶어 블록 대각선 행렬을 구성한다. 이 블록은 요소 내부에서 강한 결합을 유지하면서도 블록 대각선 우세성을 확보해, 전통적인 스칼라형 Block‑Jacobi 보다 Δt 에 대한 민감도가 크게 감소한다.

수치 실험에서는 다항 차수 p=3, 2차 단위 정사각형 영역을 50800개의 다형격자로 분할해 h 를 점진적으로 감소시켰다. 조건수 측정 결과, Δt 가 10⁻⁴ 이하로 작아질 때 A* 의 조건수가 ≈1/Δt 로 증가했으며, 집합형 전처리 적용 시 조건수가 약 23배 정도 낮아졌다. 반복 횟수 표에서는 표준 CG가 Δt=10⁻⁶, h≈0.09 일 때 3000회 이상 필요했지만, Deflated CG는 60100회, 집합형 Block‑Jacobi 전처리와 CG 결합은 5565회 수준으로 Δt 와 h 에 거의 무관한 수렴 특성을 보였다.

결론적으로, 의사응력 PolyDG 이산화에서 발생하는 시간 스텝 의존성 문제는 두 제안된 선형 솔버가 효과적으로 해결한다. 향후 연구는 영공간 전처리와 다중격자 방법을 결합해 공간 격자 크기 h 에 대해서도 조건수와 반복 횟수를 독립적으로 유지하는 완전한 시간‑공간 강인 솔버를 구축하는 방향으로 진행될 예정이다.