A linear, fully decoupled, and unconditionally energy-stable SAV-FEM for the Cahn--Hilliard--Navier--Stokes model

문제 정의
본 논문은 두 상의 비압축성 유동을 기술하는 결합된 비선형 편미분 방정식 모델인 Cahn–Hilliard–Navier–Stokes (CHNS) 시스템의 수치 시뮬레이션을 다룬다. 이 시스템은 다음과 같은 수치 이산화 측면에서의 상당한 과제를 포함한다:

1. 비선형 결합: 속도장(velocity field)과 상 변수(phase variable)는 상 방정식의 대류 항과 운동량 방정식의 응력 항을 통해 서로 결합되어 있다.
2. 비선형성 및 강성(Stiffness): 시스템은 비선형 에너지 포텐셜(Ginzburg–Landau 자유 에너지)을 특징으로 하며, 작은 계면 폭 파라미터로 인해 발생하는 강성을 가진다.
3. 에너지 안정성: 선형적이고 완전히 분리된 정식화(formulation)를 추구할 때, 이산 수준에서 물리적인 에너지 소산 법칙을 보존하는 스킴을 설계하는 것은 어렵다.

기존 방법들은 이러한 비선형성을 처리하기 위해 다중 스칼라 보조 변수(MSAV) 또는 제로 에너지 기여(ZEC) 기법이 결합된 보조 변수 활용 방식에 의존하는 경우가 많다. 이러한 접근 방식들은 효과적이기는 하지만, 비선형 시스템을 여러 번 분해하거나 여러 개의 보조 상 미분 방정식을 풀어야 하므로 계산 복잡도를 증가시킨다.

방법론
저자들은 스칼라 보조 변수(SAV) 접근법과 Implicit–Explicit (IMEX) Euler 시간 이산화를 결합한, 선형적이고, 완전히 분리되며, 무조건적으로 에너지 안정적인 완전 이산 유한 요소 스킴을 제안한다.

• 재정식화(Reformulation): $r(t) = \sqrt{E_1[\phi(t)]}$인 단일 스칼라 보조 변수를 도입하여, Cahn–Hilliard와 Navier–Stokes 방정식 모두에서 발생하는 모든 비선형 항을 동시에 재정식화한다. 표준 SAV 접근 방식과 달리, 저자들은 특정 보정 항을 포함하는 이 단일 변수에 대한 새로운 업데이트 방정식을 도입한다. 이를 통해 여러 개의 보조 변수로 분리하지 않고도 다양한 비선형성을 처리할 수 있다.
• 이산화:
  • 시간(Temporal): IMEX Euler 스킴은 선형 항은 암시적(implicit)으로, 비선형 항은 명시적(explicit)으로 처리한다.
  • 공간(Spatial): 표준 유한 요소법(FEM)이 사용된다. 속도와 압력 공간은 inf-sup 조건을 만족하며(예: MINI 요소 또는 Taylor-Hood 요소), 상 필드는 연속적인 조각별 선형 요소를 사용한다.
• 알고리즘 구조: 결과적인 완전 이산 스킴은 각 시간 단계마다 세 개의 독립적인 부문제(subproblem)로 분해된다:
  1. 선형 Cahn–Hilliard 유형의 시스템.
  2. 선형 Stokes 유형의 시스템.
  3. 보조 변수를 위한 스칼라 이차 대수 방정식.
     결정적으로, 속도와 압력은 각 단계에서 결합된 Stokes 시스템으로 직접 풀리며(압력 보정 투영법을 사용하는 대신), 이는 오차 분석을 단순화한다.

주요 기여 및 결과

1. 단일 보조 변수의 효율성: 본 연구의 주요 기여는 특정 업데이트 전략과 결합된 단일 스칼라 보조 변수가 CHNS 시스템의 모든 비선형성을 처리하기에 충분하다는 것을 입증한 것이다. 이는 여러 개의 보조 변수를 요구하는 MSAV 또는 SAV-ZEC 방법들에 비해 계산 오버헤드를 줄여준다.
2. 무조건적 에너지 안정성: 저자들은 제안된 스킴이 무조건적 이산 에너지 소산 법칙을 만족함을 증축한다. 수정된 이산 에너지는 시간 단계 크기에 관계없이 각 시간 단계에서 단조 감소한다.
3. 최적 오차 추정(Optimal Error Estimates): 표준 정칙성 가정($u, \phi \in H^2$) 하에서, 본 논문은 완전 이산 근사에 대한 최적 차수의 $L^2$ 오차 추정치를 확립한다.
   • 상 필드, 화학 포텐셜, 그리고 압력은 최적의 수렴율을 달enc한다.
   • 속도는 최적 차수의 $L^2$ 수렴을 달성한다.
   • 중요한 기술적 성과는 비전역적 리프시츠(non-globally Lipschitz) 특성을 가진 비선형 포텐셜 $g(\phi) = F'(\phi)$를 다룬 것이다. 저자들은 수치 해의 **균등 $L^\infty$ 유계성(uniform $L^\infty$ boundedness)**을 증명하였으며, 이를 통해 자유 에너지에 대한 제한적인 가정 없이 오차 분석을 진행할 수 있었다.
4. 수치 검증: 수치 실험을 통해 제안된 방법의 이론적 수렴율(공간 이산화에 대해 $L^2$에서 $O(h^2)$, $H^1$에서 $O(h)$ 관찰)을 확인하였으며, 상 분리 및 조립(coarsening) 역학을 포함한 다양한 테스트 케이스에서 수정된 이산 에너지가 단조 감소함을 입증하였다.

의의
본 논문은 제안된 접근 방식이 CHNS 시스템에 대한 기존 SAV 기반 방법들보다 더 효율적이고 단순한 대안을 제공한다고 주장한다. 단일 보조 변수와 통합된 업데이트 전략을 활용함으로써, 이 방법은 시스템을 여러 번 분해해야 하는 복잡성을 피한다. 또한, 투영 단계 없이 Stokes 부체계를 직접 해결함으로써 이론적 오차 분석을 단순화하며, 높은 정칙성을 요구할 수 있는 일부 투영 기반 방법들과 비교하여 더 약한 정칙성 가정($\phi \in H^2$) 하에서도 최적의 오차 추정을 가능하게 한다. 본 연구는 계산적으로 효율적(선형 및 분리형)이면서도 이론적으로 엄밀한(무조건적 안정성 및 최적의 오차 경계) 두 상 유동 시뮬레이션을 위한 견고한 프레임워크를 제공한다.