[논문 리뷰] Multigrid Methods for Constrained Minimization Problems and Application to Saddle Point Problems

Long Chen|arXiv (Cornell University)|2016. 01. 15.

Advanced Numerical Methods in Computational Mathematics참고 문헌 45인용 수 7

한 줄 요약

이 논문은 편미분방정식(Poisson, Darcy, Stokes 방정식 포함)의 혼합 유한요소 해석에서 유래하는 제약 최소화 문제와 안장점 시스템을 위한 다중격자 방법을 제시한다. 제약 공간의 다중레벨 분해를 기반으로 하는 순차 하위공간 최적화로 간주하고, 스토케르-샤르스형 다중격자 스무스러이터(다중 겹침 슈워츠 스무스러이터)를 사용함으로써, 안정적 분해와 강화된 코시-슈와르츠 조건을 가정할 경우, 전체 정규성 조건이 없더라도 균일 수렴성을 증명한다. 또한, 정확한 해가 아닌 하나의 경사계산 단계만을 사용하는 경우에도 비일치 유한요소 방법으로의 프레임워크 확장을 성공적으로 수행한다.

ABSTRACT

The first order condition of the constrained minimization problem leads to a saddle point problem. A multigrid method using a multiplicative Schwarz smoother for saddle point problems can thus be interpreted as a successive subspace optimization method based on a multilevel decomposition of the constraint space. Convergence theory is developed for successive subspace optimization methods based on two assumptions on the space decomposition: stable decomposition and strengthened Cauchy-Schwarz inequality, and successfully applied to the saddle point systems arising from mixed finite element methods for Poisson and Stokes equations. Uniform convergence is obtained without the full regularity assumption of the underlying partial differential equations. As a byproduct, a V-cycle multigrid method for non-conforming finite elements is developed and proved to be uniform convergent with even one smoothing step.

연구 동기 및 목표

• 혼합 유한요소 방법에서 유래하는 안장점 시스템에서 발생하는 제약 최소화 문제를 위한 다중격자 방법을 개발한다.
• 안정적 분해(SD)와 강화된 코시-슈와르츠(SCS)라는 두 가지 핵심 가정을 기반으로 하는 순차 하위공간 최적화의 수렴 이론을 수립한다.
• Poisson, Darcy, Stokes 방정식의 혼합 유한요소에 대해 전체 정규성 조건 없이도 V-사이클 다중격자 방법의 균일 수렴성을 증명한다.
• 비일치 유한요소 방법으로의 프레임워크 확장을 수행하고, 오직 한 번의 스무스링 단계만으로도 균일 수렴성을 증명한다.
• 다중격자 방법에서 곱셈형 슈워츠 스무스러이터의 이론적 근거를 제공하여, 이전 분석이 덧셈형에 국한되어 있던 격차를 메운다.

제안 방법

• 방법은 제약 공간 $\mathcal{K}$의 다중레벨 분해를 사용하며, $\mathcal{K} = \sum_{i=1}^{N} \mathcal{K}_i$로 표현된다. 각 하위공간 $\mathcal{K}_i$에서 국소적 제약 최소화 문제를 풀게 된다.
• 각 국소 문제는 작은 안장점 시스템을 푸는 것과 동치이며, 다중 겹침 슈워츠 스무스러이터(반카 유형)로 해석된다.
• 에너지 안정적 분해(SD)와 강화된 코시-슈와르츠(SCS) 조건을 가정할 경우, 수렴성이 증명되며, 수축 인자는 $1 - \frac{1}{1 + C_A C_S}$로 나타난다.
• 국소 문제의 정확한 해가 아닌, 근사해를 허용하는 프레임워크를 제공하며, 정확한 해 대신 한 번의 경사계산 단계를 사용한다. 또한, 최적의 완화 파rameter는 최소화 관점에서 유도된다.
• 비일치 방법의 경우, 겹침 슈워츠 스무스러이터와 네스팅된 플럭스 공간을 통한 격자 간 전달 연산자를 사용하며, 혼합 방법과의 등가성을 활용한다.
• 분석은 Poisson 및 Darcy 방정식에 대한 H(div) 유한요소와 Stokes 방정식에 대한 발산이 없는 내재된 유한요소에 적용된다.

실험 결과

연구 질문

• RQ1 순차 하위공간 최적화 기반의 다중격자 방법이 혼합 유한요소 방법에서 유래하는 안장점 시스템에 대해 균일 수렴성을 달성할 수 있는가?
• RQ2 기존의 다중격자 분석에서 일반적으로 요구하는 전체 정규성 조건 없이도 수렴 이론이 유지되는가?
• RQ3 비일치 유한요소 방법으로의 프레임워크 확장이 오직 한 번의 스무스링 단계만으로도 가능한가?
• RQ4 이전 분석이 덧셈형에 국한되어 있었기 때문에, 다중격자 방법에서 곱셈형 슈워츠 스무스러이터의 이론적 근거는 어떻게 마련할 수 있는가?
• RQ5 크루아제-라비에르 비일치 방법과 혼합 방법 간의 등가성을 활용하여, 포isson 방정식에 대해 균일 수렴성을 보장하는 V-사이클 다중격자 방법을 설계할 수 있는가?

주요 결과

• 안정적 분해(SD)와 강화된 코시-슈와르츠(SCS) 조건을 가정할 경우, Poisson 및 Darcy 방정식의 혼합 유한요소 해석에 대해 V-사이클 다중격자 방법이 전체 정규성 조건 없이도 균일 수렴성을 확보한다.
• 비일치 유한요소 방법의 경우, 오직 한 번의 스무스링 단계만으로도 균일 수렴성을 달성하며, 이는 이전의 다중격자 수렴 이론이 다루지 못한 결과이다.
• 3차원에서 라비아르트-라비에르 요소에 대해 새로운 다중레벨 분해를 제안하여, 이중성 논증이 필요 없어지며, 따라서 전체 정규성 조건을 완화한다.
• 내재된 발산이 없는 유한요소를 사용한 스토크스 방정식에 대해서도 수렴 이론을 확장하여, 전체 정규성 조건 없이도 균일 수렴성을 증명하였으며, 한 번의 스무스링 단계만으로도 V-사이클에 대해 성립한다.
• Vanka 스무스러이터의 에너지 최소화 해석이 명확하게 드러나며, 최소화 원리에 기반한 최적의 완화 파rameter 선택이 정당화된다.
• 다중레벨 H(div) 공간에서 강화된 코시-슈와르츠 부등식과 안정적 분해 조건이 성립함을 확인하여, 혼합 방법의 수렴 증명이 가능해졌다.