Hybrid CutFEM-IFEM model for embedded and immersed interfaces with triply-cut cells for the 2D Poisson problem

기술 요약: 임베디드 및 침투 인터페이스를 위한 하이브리드 CutFEM-IFEM 모델

문제 정의
본 논문은 세 가지 서로 다른 상(phase)($\Omega_a, \Omega_b$ 두 개의 유체 상과 $\Omega_s$ 한 개의 고체 상)을 포함하는 2D Poisson 문제의 수치적 해법을 다룬다. 도메인은 두 가지 유형의 인터페이스를 포함한다:

1. 침투 인터페이스 ($\Gamma_f$): 두 유체 상을 분리하며, 재료 특성($\mu$), 해($u$), 그리고 플럭스($\mu \partial u / \partial n$)의 불연속성을 특징으로 한다.
2. 임베디드 인터페이스 ($\Gamma_g$): 고체 영역의 경계를 나타내며, 이 영역은 계산 도메인에서 제거되므로 경계 조건(Dirichlet 또는 Robin)을 부과해야 한다.

핵심 과제는 두 인터페이스가 동일한 요소(element)를 통과할 수 있는 단일 고정 배경 격자(background mesh) 위에서 이러한 상호작용을 시뮬레이션하는 것이다. 이 경우 "삼중 절단 셀(triply-cut cells)"이 발생한다. 전통적인 격자 적합(mesh-fitting) 방식은 인터페이스가 이동하거나 변형될 때 계산 비용이 많이 들고 취약하다. 따라서 언핏티드(unfitted) 방식이 선호되지만, 기존 기술들은 하나의 절단된 셀 내에서 강한 불연속성(IFEM을 통해 처리)과 임베디드 경계(CutFEM을 통해 처리)를 동시에 처리할 때, 특히 적분과 일관성 측면에서 어려움을 겪는 경우가 많다.

방법론
저자들은 단일 언핏티드 격자 위에서 두 방식의 강점을 활용하는 하이브리드 CutFEM–IFEM 정식화를 제안한다:

• 기저 함수 및 인리치먼트(Enrichment):

  • 유체-유체 인터페이스($\Gamma_f$)에 의해 절단된 셀의 경우, 이 방법은 IFEM 기저 함수를 사용한다. 이 함수들은 인터페이스를 가로질러 연속적이지만, 요소의 모서리(edge)를 가로질러서는 불연속적이며, 해와 플럭스의 점프 조건을 정확히 만족하도록 구성된다.
  • 비균질 점프(주어진 점프 $J_s$ 및 $J_f$)를 처리하기 위해, 시행 함수(trial function)는 특정 "리프트(lift)" 함수($j_h$는 해의 점프를 위해, $q_h$는 플럭스 점프를 위해)로 보강된다. 이러한 인리치먼트 함수는 정확한 점프를 재현하여, 오차가 점프의 크기에 의존하지 않도록 보장한다.
  • 임베디드 고체 경계($\Gamma_s$)의 경우, 자유도(degree of freedom)를 추가하지 않고 경계 조건(Dirichlet 또는 Robin)을 약하게 부과하기 위해 Nitsche 방식이 사용된다.

• 약형(Weak Formulation) 및 안정성:

  • 약형은 임베디드 경계를 위한 표준 Nitsche 항과 작은 절단 비율에 대한 시스템 안정화를 위한 고스트 페널티(ghost penalty) 항을 포함한다.
  • 삼중 절단 셀을 위해 중요한 수정 사항이 도입된다. 고체 경계가 두 유체 절단 요소 사이의 공유 모서리 일부를 제거하기 때문에, 표준 IFEM 모서리 적분은 일관성이 결여된다. 저자들은 남은 유체 부분의 플러스를 고려하는 특정 일관성 보정 항($c(u_h, v_h)$)을 유도하여 최적의 수렴성을 회복시킨다.

• 삼중 절단 셀에서의 적분:

  • 본 논문은 여러 인터페이스에 의해 절단된 셀에서의 적분을 처리하기 위해 등가 다항식(equivalent polynomial) 기법을 확장한다.
  • 셀을 하위 삼각형으로 세분화하는 대신, 이 방법은 Heaviside 함수와 Dirac delta 함수를 등가 다항식으로 대체한다.
  • 삼중 절단 셀의 경우, 저자들은 절단에 의한 여덟 가지 가능한 위상적 구성(격리된 노드의 여부 및 삼중점 존재 여부에 기반)을 분류한다. 이들은 이러한 복잡한 기하학적 구조 내에서 각 상(고체, 유체 $a$, 유체 $b$)에 대한 모멘트 피팅(moment-fitting) 프로세스의 우변 벡터를 계산하기 위한 해석적 공식을 유도한다.

주요 기여

1. 하이브리드 정식화: 단일 격자 위에서 유체-유체 인터페이스를 위한 IFEM과 유체-고체 경계를 위한 Nitsche 기반 CutFEM을 결합한 통합 프레임워크를 제공하며, 선형($p=1$) 및 이차($p=2$) 요소 모두에 적용 가능하다.
2. 삼중 절단 적분: 두 개의 서로 다른 인터페이스에 의해 절단된 셀을 처리하기 위해 등가 다항식 적분법을 확장하였으며, 삼중점을 포함한 모든 위상적 구성을 다룬다.
3. 일관성 보정: 임베디드 경계에 의해 절단된 모서리를 고려하는 새로운 보정 항을 통해 삼중 절단 셀의 일관성을 확보하고 최적의 수렴성을 달성한다.
4. 정확한 점프 표현: 비균질 해 및 플럭스 점프를 처리하기 위해 인리치먼트 함수를 사용하여, 오차의 크기를 점프의 크기로부터 분리한다.

수치 결과
본 방법은 다양한 계수 대비(최대 $\mu_b/\mu_a = 1000$)와 점프 조건(균질 및 비균질)을 가진 제조된 해(manufactured solutions)를 통해 검증되었다.

• 수렴성: 이 방법은 선형 및 이차 요소 모두에 대해 $L^2$ 노름에서 $p+1$ 차수의 최적 수렴율을 달ates한다.
• 정확도: 단일 배경 격자를 사용함에도 불구하고, 두 인터페이스에 명시적으로 격자를 맞춘(fitted) 경우와 유사한 정확도를 보인다.
• 강건성: 높은 계수 대비에서도 안정적이고 정확한 결과를 유지한다.
• 한계점: 곡선형 임베디드 경계에 Robin 조건이 적용될 경우, 이차 수렴율($p=2$)이 2차로 저하된다. 저자들은 이를 레벨 셋 함수 내의 곡선 경계가 선형적으로 표현되었기 때문으로 보고하며, 최적의 수렴율을 회복하기 위해서는 더 높은 차수의 경계 표현이 필요함을 언급한다.

의의 및 주장
본 논문은 인터페이스가 이동하고 상호작용하는 복잡한 다상(multi-phase) 문제를 위한 견고하고 효율적인 수치 도구를 제공한다고 주장한다. IFEM과 CutFEM을 결합함으로써, 이 방법은 격자 재생성 없이도 높은 정확도를 유지하면서 격자 불일치 문제를 피할 수 있다. 저자들은 IFEM 처리를 통한 임머스드 인터페이스 처리와 Nitsche 처리를 통한 임베디드 경계 처리를 동일한 절단 셀 내에서 결합한 최초의 연구임을 강조하며, 특히 삼중 절단 셀의 적분 및 일관성 문제를 해결하였다. 결과는 이 하이브리드 접근 방식이 언핏티드 방식의 유연성을 갖추면서도 격자 적합 방식의 정확도를 제공함을 입증하며, 이는 움직이는 유체-유체 인터페이스 및 복잡한 고체 경계가 포함된 응용 분야에 적합함을 보여준다.