두 상 흐름을 위한 파라메트릭 유한요소와 XFEM 기반 비뉴턴·탄성 시뮬레이션

초록

본 논문은 Oldroyd‑B 모델을 이용한 두 상 Navier‑Stokes 흐름에 대해 파라메트릭 유한요소와 XFEM을 결합한 비맞춤(unfitted) 방법을 제안한다. 인터페이스는 저차원 파라메트릭 메쉬로 직접 추적하고, 압력은 XFEM 보강으로 불연속을 정확히 포착한다. 시간·공간 이산화는 무조건적인 안정성과 해 존재성을 보장하며, 체적 보존과 에너지 소산 법칙을 만족한다. 수치 실험을 통해 2차원·3차원 사례에서 방법의 효율성과 정확성을 확인한다.

상세 분석

이 연구는 두 상 유체 사이의 자유계면을 명시적으로 추적하는 파라메트릭 유한요소 기법을 기반으로, Oldroyd‑B 모델에 의해 기술되는 점탄성 효과를 동시에 고려한다. 기존의 레벨셋이나 위상장(field) 방식과 달리, 인터페이스는 저차원( d‑1 차원) 메쉬를 직접 이동시켜 구현되며, 이는 자유계면의 기하학적 정확성을 크게 향상시킨다. 그러나 파라메트릭 메쉬는 변형이 심해질 경우 격자 퇴화가 발생할 위험이 있다. 이를 해결하기 위해 저자들은 Dziuk‑style의 정상속도 기반 이동과 접선 자유도를 이용한 재배열 없는 메쉬 품질 유지 전략을 채택하였다.

본 논문의 핵심 수치적 공헌은 두 가지이다. 첫째, 압력과 속도 사이의 불연속을 정확히 포착하기 위해 XFEM(확장 유한요소법)을 적용하였다. 압력 공간에 인터페이스를 가로지르는 특수한 풍선형 함수 하나를 추가함으로써, 반대상(phase) 간의 표면장력에 의한 압력 점프를 정확히 재현하고, 반면에 체적 보존을 반정밀 수준에서 완전하게 유지한다. 둘째, Oldroyd‑B 방정식에 포함된 비선형 대류항과 응력 확산항을 다루기 위해 반명시적(sem-implicit) 시간 스키마와 Schur 보완을 이용한 고정점 반복을 설계하였다. 이 과정에서 연속적인 에너지 소산 법칙을 이산화 수준에서도 유지하도록 설계했으며, 이를 통해 “무조건적인(stable for any time step) 안정성”과 “해 존재성”을 엄격히 증명하였다.

수학적 분석에서는 먼저 연속 모델에 대한 에너지 식을 도출하고, 이를 기반으로 이산 에너지 불등식(energy inequality)을 구축한다. 특히, Oldroyd‑B의 탄성 에너지 G/2·Tr(B−ln B−I) 를 사용함으로써 B 텐서가 양정(positive‑definite)임을 보장하고, α·ΔB 확산항을 통해 정규화된 해의 정규성을 확보한다. 인터페이스 조건은 속도 연속성, 응력 점프, 그리고 표면장력에 의한 평균곡률 항을 포함한다. 이러한 조건들을 변분 형태로 변환한 뒤, 파라메트릭 메쉬와 XFEM 보강을 적용한 불맞춤(fitted‑unfitted) 공간에 투영한다.

알고리즘 흐름은 다음과 같다. (1) 현재 시각의 인터페이스와 물리량(u, p, B)을 이용해 선형화된 Navier‑Stokes‑Oldroyd‑B 시스템을 구성한다. (2) Schur 보완을 통해 압력·속도와 탄성 텐서 B를 분리하고, 고정점 반복으로 비선형성을 처리한다. (3) 인터페이스 메쉬는 정상속도 V·ν에 의해 이동하고, 재파라메트라이징 없이 접선 자유도를 이용해 격자 품질을 유지한다. (4) 압력 공간에 XFEM 풍선 함수를 삽입해 체적 보존을 강화한다.

수치 실험에서는 2D “스위스 치즈”와 3D “스피어‑드롭” 문제를 다루며, 높은 Weissenberg 수(>100)에서도 수렴성과 안정성을 확인한다. 특히, 압력 점프와 인터페이스 위치가 정확히 재현되며, 체적 오차는 10⁻⁶ 수준으로 억제된다. 또한, α>0 인 경우 응력 확산이 수치적 진동을 감소시켜 고Weissenberg 수치에서도 안정적인 시뮬레이션이 가능함을 보여준다.

전반적으로, 이 논문은 파라메트릭 인터페이스 추적, XFEM 보강, 그리고 Oldroyd‑B 점탄성 모델을 통합한 최초의 무조건 안정적인 두 상 흐름 시뮬레이션 프레임워크를 제공한다는 점에서 학술적·실용적 의의가 크다.