Convergence of a CutFEM for fluid--structure interaction with a deforming interface

기술 요약: 변형되는 인터페이스를 가진 유체-구조 상호작용을 위한 CutFEM의 수렴성

문제 설정
본 연구는 비압축성 점성 유체와 선형 탄성 고체를 포함하는 유체-구조 상호작용(FSI) 문제의 수치적 근사를 다룬다. 이 시스템은 고체의 변형에 의해 유체 영역이 정의되는 움직이는 인터페이스를 특징으로 한다. 지배 방정식은 유체 영역에서의 비압축성 나비에-스토크스 방정식(오일러 기술)과 고체 영역에서의 선형 탄성 방정식(라그랑주 기술)으로 구성된다. 결합은 운동학적 연속성(속도 일치)과 진화하는 인터페이스 상의 응력 동적 평형을 통해 강제된다.

이러한 문제를, 특히 큰 변형이 발생하는 경우 수치적으로 해결하는 데 있어 주요한 과제는 계산 격량에 대한 인터페이스의 상대적 움직임이다. 전통적인 적합 격자 방법(예: ALE)은 큰 회전이나 변형 하에서 격자 저하 문제를 겪는다. 미적합 방법(Immersed Boundary 또는 Fictitious Domain 방법 등)은 격자 재생성을 피할 수 있으나, 불연속성이나 경계 조건을 해상하지 못해 인터페이스 근처에서의 정확도가 떨어지는 경우가 많다.

방법론
저자들은 이 완전 결합된 FSI 시스템을 위해 반이산(semi-discrete) 임머스드 CutFEM(Cut Finite Element Method)을 제안하고 분석한다. 이 방법은 다음과 같은 핵심 요소들을 사용한다:

1. 미적합 격자(Unfitted Mesh): 유체는 고체 경계와 일치하지 않는 고정된 배경 격자 상에서 이산화된다. 고체는 그 참조 도메인 상에서 등매개변수(isoparametric) 적합 격자를 사용하여 이산화된다.
2. 라그랑주 승수 정식화: 유체-고체 결합 조건은 라그랑주 승수를 사용하여 약하게 강제된다. 이 승수는 참조 고체 인터페이스로 역추적(pull back)된 유체 트랙션(traction)을 나타낸다. 운동학적 조건($u = \partial_t d \circ X^{-1}$)은 현재의 인터페이스 상에서 강제되며, 트랙션 균형은 참조 형상을 포함하는 적분 관계로 정식화된다.
3. 고스트 페널티 안정화(Ghost Penalty Stabilization): 임의로 작은 컷 요소들에 대해 안정성과 강건성을 보장하기 위해, 속도 및 압력 정식화에 전역 고스트 페널티 항이 추가된다. 이 항들은 배경 격자의 내부 면을 가로지르는 법선 미분값의 점프를 페널티로 부과하여, 확장된 유체 영역까지 유한 요소 공간의 안정성을 효과적으로 확장한다.
4. 반이산 스킴(Semi-Discrete Scheme): 본 방법은 시간 연속적인 공간 반이산 방식이다. 유체 속도, 압력, 고체 변위, 그리고 라그랑주 승수는 각각 차수가 $k$(속도 및 변위)와 $k-1$(압력)인 유한 요소 공간을 사용하여 근사된다.

주요 기여 및 이론적 프레임워크
본 논문의 핵심 기여는 진정으로 변형되는 인터페이스를 가진 완전 결합된 FSI 시스템에 적용된 미적합 유한 요소 방법의 첫 번째 엄밀한 사전(a priori) 오차 분석이다. 이전의 엄밀한 분석들은 주로 선형화된 모델, 작은 변위 근사, 또는 적합 격자 접근법에 국한되어 있었다.

수렴성 증명은 몇 가지 새로운 분석적 요소들에 의존한다:

• 보간된 형상(Interpolated Configuration): 저자들은 정확한 인터페이스 $\Gamma[X]$, 보간된 인터페이스 $\Gamma[X^*_H]$, 그리고 수치적 인터페이스 $\Gamma[X_H]$를 비교하기 위해 정확한 해의 보간된 형상 $X^*_H$를 도입한다. 이를 통해 인터페이스 변위에 의해 발생하는 오차를 체계적으로 정량화할 수 있다.
• 균등 추정(Uniform Estimates): 분석은 움직이는 컷 도메인 상에서 균등 노름 등가성, 트레이스 부등식, 역 추정 및 inf-sup 조건을 확립한다. 결정적으로, 이러한 추정치들은 인터페이스가 배경 격자 내에서 위치하는 방식에 관계하여, 그리고 수치적 변형이 정확한 변형과 충분히 가깝다는 조건 하에 시간과 관계없이 균등하다.
• 오차의 결합(Coupling of Errors): 해결해야 할 중요한 분석적 어려움은 유체 오차와 기하학적 오차(고체 변위) 사이의 양방향 의존성이다. 고체 변위의 오차는 수치적 유체 영역과 인터페이스를 변화시키며, 이는 다시 유체 해에 영향을 미친다. 저자들은 인터페이스 상의 투영된 운동학적 불일치에 대한 안정성 추정을 도출하고 연속성-시간(continuation-in-time) 논증을 활용함으로써 이 루프를 닫는다.
• 연속성 논증(Continuation Argument): 증명은 최대 시간 $t^*$까지의 오차에 대한 사전 경계(특히 $W^{1,\infty}$ 및 $L^\infty$ 경계)를 가정한다. $t^*$와 무관한 상수를 갖는 에너지 추정을 증명하고, $k \ge 3$일 때 도출된 수렴 속도가 가정된 경계를 함의함을 보여줌으로써, 전체 시간 구간 $[0, T]$에 대한 추정치의 유효성을 확장한다.

주요 결과
정확한 해의 충분한 정칙성(regularity)과 격자 크기가 $c_1 H \le h \le c_2 H$ (여기서 $h$와 $H$는 각각 유체와 고체의 격자 크기)를 만족한다는 가정 하에, 저자들은 차수가 $k \ge 3$인 유한 요소에 대해 다음의 수렴 속도를 증명한다:

• 유체 속도 ($u_h$) 및 고체 변위 ($d_H$): 자연 에너지 노름에서 $k$ 차수의 수렴.
• 유체 압력 ($p_h$): ($L^2$ 노름과 관련된) 에너지 노름에서 $k-1$ 차수의 수렴.
• 라그랑주 승수 (인터페이스 트랙션): $H^{-1/2}$ 노름에서 $k-1$ 차수의 수렴.

구체적인 오차 경계는 다음과 같다:
$$\max_{t \in [0,T]} \left( \|u_h - u\|_{L^2}^2 + \|d_H - d \circ \Phi_H\|_{H^1}^2 + \|\partial_t d_H - \partial_t d \circ \Phi_H\|_{L^2}^2 \right) + \int_0^T \|u_h - u\|_{H^1}^2 dt + H^2 \int_0^T (\|p_h - p\|_{L^2}^2 + \|\sigma_H - \sigma \circ \Phi_H\|_{H^{-1/2}}^2) dt \le C(h^{2k} + H^{2k})$$

의의
본 논문은 변형되는 인터페이스를 가진 완전 결합된 FSI 시스템에 대해 미적합 유한 요소 방법의 첫 번째 엄밀한 수렴 분석을 수행했다고 주장한다. 이전 연구들은 주어진 움직이는 도메인이나 단순화된 결합 문제에 대해서는 분석하였으나, 완전 결합된 FSI 시스템에서의 비선형성과 움직이는 기하 구조와 유체 해 사이의 피드백 루프는 미적합 맥락에서 엄밀하게 다뤄진 바 없다. 이 결과는 복잡한 FSI 응용 분야에서 CutFEM을 사용할 수 있는 이론적 토대를 제공하며, 작은 변위나 고정 도메인 근사에 의존하지 않고도 이 방법의 정확성과 안정성을 검증한다. $k \ge 3$이라는 제한은 증명에 사용된 특정 역 추정(inverse estimates)과 관련된 기술적인 사항이며, 저자들은 향-개선을 통해 이 결과가 $k=2$까지 확장될 수 있음을 시사한다.