가상 영역 기반 유체‑구조 상호작용의 안정성 및 조건수 분석

초록

본 논문은 가상 영역(Fictitious Domain) 접근법을 이용한 유체‑구조 상호작용(FSI) 문제에 대해, 분산 라그랑주 승수를 도입한 비맞춤형(Unfitted) 형식을 제안한다. 시간 이산화 후 얻어지는 정적 문제에 대해, 격자 정밀화와 인터페이스 위치에 관계없이 조건수가 일정하게 유지됨을 이론적으로 증명하고, 작은 절단 셀(small cut cells) 존재에도 불구하고 안정성이 보장됨을 수치 실험으로 확인한다.

상세 분석

이 연구는 기존의 피팅된 ALE 방식이 겪는 격자 변형 문제를 회피하고자, 유체와 구조를 각각 독립적인 격자로 전개한 뒤, 전체 영역을 가상 영역으로 확장한다. 핵심 아이디어는 구조 변형을 라그랑주 좌표계에서 기술하고, 그 변형 맵 X(s,t)를 통해 구조와 유체 사이의 운동 연속성을 분산 라그랑주 승수 λ∈Λ 로 강제한다는 점이다. 이때 λ와 변형 맵 사이의 쌍대 쌍을 정의하는 bilinear form c(·,·) 은 두 가지 선택(Λ₀와 Λ₁) 중 하나를 사용할 수 있으며, 두 경우 모두 c(μ,Z)=0 ⇒ Z=0 라는 강성 조건을 만족한다.

시간 이산화는 후진 오일러 스키마를 적용해 1차 정확도를 유지하면서, 비선형 대류항을 무시하고 선형화된 Stokes‑형 시스템을 얻는다. 이 시스템은 4개의 변수(u, p, X, λ) 로 구성된 saddle‑point 문제로 재구성되며, 연산상으로는 유체 영역 Ω와 구조 기준 영역 B 에 대한 각각의 변분식에 c(·,·) 로 연결된 coupling term 이 추가된다.

논문의 주요 수학적 기여는 다음과 같다. 첫째, 연속 및 이산 문제 모두에 대해 Lax‑Milgram‑type의 inf‑sup 조건을 만족함을 증명한다. 여기서 핵심은 c(·,·) 가 정의역 전체에 대해 강한 연속성을 가지고, λ와 X‑u(X) 사이의 제약을 정확히 반영한다는 점이다. 둘째, 조건수(Condition Number) 상한을 도출할 때, 격자 크기 h_Ω, h_B 와는 무관하게, 오직 물성 파라미터(밀도 차 δρ, 점성 ν, 강성 κ)와 시간 간격 Δt 에만 의존한다는 것을 보인다. 이는 작은 절단 셀(인터페이스가 격자를 거의 전부 가로지르는 경우) 이 존재하더라도 행렬 스펙트럼이 급격히 악화되지 않음을 의미한다.

수치 실험에서는 (i) 정밀한 교차 적분을 통한 정확한 coupling term 구현과 (ii) 단일 사다리꼴 적분을 이용한 근사 구현을 비교한다. 두 경우 모두 조건수가 동일하게 유지되며, 특히 인터페이스가 격자 셀을 거의 전부 차지하는 “sliver cell” 상황에서도 해의 수렴성 및 안정성이 손상되지 않는다. 또한, 시간 전진 테스트에서는 Δt 에 대한 제한 없이 무조건적인 안정성을 확인한다. 이는 기존 Nitsche‑type 방법에서 흔히 요구되는 ghost‑penalty 같은 인공적인 안정화 항이 필요 없음을 시사한다.

결론적으로, 이 논문은 가상 영역 기반 비맞춤형 FSI 모델이 작은 절단 셀에 민감하지 않으며, 조건수가 격자 정밀화와 인터페이스 위치에 독립적인 특성을 갖는다는 강력한 이론적·수치적 근거를 제공한다. 이는 복잡한 생체 조직 시뮬레이션이나 이동·변형이 큰 구조물에 대한 고효율 전산 모델링에 중요한 전진을 의미한다.