[논문 리뷰] Numerical computations of split Bregman method for fourth order total variation flow

Yoshikazu Giga, Yuki Ueda|arXiv (Cornell University)|2019. 06. 11.

Fluid Dynamics and Turbulent Flows참고 문헌 51인용 수 10

한 줄 요약

이 논문은 후진 오일러 시간 이산화와 조각별 상수 공간 이산화를 사용하여 제4차 총변동 흐름과 스폰의 제4차 모델을 해결하기 위한 새로운 분할 브레그만 프레임워크를 제안한다. 스폰의 모델을 위해 새로운 셸링 연산자를 도입하고, 주기적 경계 조건 하에서 효과적인 수치적 해를 도출하여 2D 시뮬레이션에서 정확한 평평한 면과 표면 완화를 달성하였으며, 매개변수 스케일링은 O(h⁻⁴) 및 O(h⁻²)를 사용하였다.

ABSTRACT

The split Bregman framework for Osher-Sol\'e-Vese (OSV) model and fourth order total variation flow are studied. We discretize the problem by piecewise constant function and compute $ abla(-\Delta_{\mathrm{av}})^{-1}$ approximately and exactly. Furthermore, we provide a new shrinkage operator for Spohn's fourth order model. Numerical experiments are demonstrated for fourth order problems under periodic boundary condition.

연구 동기 및 목표

• 제4차 총변동 흐름과 스폰의 제4차 모델을 효율적으로 계산할 수 있는 수치적 방법을 개발한다. 이는 매우 특이하고 계산하기 어려운 문제이다.
• 두 번째 차수 문제에 잘 정립된 분할 브레그만 프레임워크를 후진 오일러 시간 스텝을 사용하여 제4차 영역으로 확장한다.
• 결정성 있는 결정 표면 완화의 안정적이고 정확한 수치적 해를 가능하게 하는 스폰의 모델에 특화된 새로운 셸링 연산자를 도입한다.
• 평탄한 면과 불연속성과 같은 핵심 물리적 특성을 유지하는 강력한 수치적 방법을 제공한다. 이는 이론적 기대와 일치한다.

제안 방법

• 1D 또는 2D 토러스에서 조각별 상수 유한요소를 사용하여 제4차 PDE의 특이성을 다루기 위해 문제를 이산화한다.
• 하위미분 공식화에 후진 오일러 방법을 적용하여 PDE를 총변동과 H⁻¹ 적합성 항을 포함한 최소화 문제로 변환한다.
• 분할 브레그만 프레임워크를 사용하여 L1 및 L2 부분을 분리하여 효율적인 교대 최소화를 가능하게 한다.
• 두 가지 방법을 구현한다: 하나는 ∇(−Δₐᵥ)⁻¹의 근사치를 사용하고, 다른 하나는 조각별 상수 함수에 대해 두 차수 B-스플라인을 이용한 정확한 역라플라시안을 사용한다.
• 분할 브레그만 프레임워크 하에서 오일러-라그랑주 방정식을 풀어 스폰의 모델을 위한 새로운 셸링 연산자를 유도하며, 이는 ℓ¹ 및 ℓᵖ 노름을 모두 포함한다.
• 일致한 매개변수 스케일링을 적용한다: 2D에서 λ = O(h⁻⁴), µ = O(h⁻²)이며, 이는 수치 실험에서 안정성과 수렴성을 보장한다.

실험 결과

연구 질문

• RQ1 특이성이 두 번째 차수 대비 더 높은 제4차 총변동 흐름과 스폰의 모델에 대해 분할 브레그만 프레임워크를 성공적으로 확장할 수 있는가?
• RQ2 조각별 상수 함수의 맥락에서 제4차 PDE의 맥락에서 역라플라시안 연산자를 효율적이고 정확하게 계산할 수 있는가?
• RQ3 스폰의 제4차 모델의 조합된 ℓ¹ 및 ℓᵖ 구조를 다루기 위해 필요한 새로운 셸링 연산자는 무엇이며, 이는 수치적 안정성을 어떻게 향상시키는가?
• RQ4 수치적 해가 이론적으로 예측한 바와 같이 평탄한 면과 불연속성과 같은 예상되는 물리적 특성을 재현하는가?

주요 결과

• 제안된 분할 브레그만 방법은 주기적 경계 조건 하에서 제4차 총변동 흐름을 성공적으로 계산하여 이론적 기대와 일치하는 평탄한 면과 불연속성을 가진 해를 도출한다.
• 등방성 및 비등방성 제4차 총변동 흐름에 대한 수치 결과는 날카운 전이를 보이며 안정적인 진화를 보이며, 비등방성 경우는 조각별 상수 구조를 유지하여 이론적 추측을 지지한다.
• 스폰의 모델을 위한 새로운 셸링 연산자는 결정 표면 완화의 효과적인 시뮬레이션을 가능하게 하여 면 형성과 부드러운 진화를 만들어내며 이는 이전 연구와 일치한다.
• 2D에서 λ = 1.25h⁻⁴ 및 µ = 5h⁻²를 사용할 경우 뛰어난 성능을 보이며, 복잡한 제4차 PDE에 대해 강력성과 수렴성을 입증한다.
• 정확한 역라플라시안 계산을 위해 두 차수 B-스플라인을 사용하는 것은 근사 방법에 비해 정확도를 향상시키며, 특히 비국소 효과를 잘 포착하는 데 유리하다.