High order ghost-FEM for incompressible Navier-Stokes equations on moving domains

1. 연구 배경 및 문제 정의 (Problem)

• 배경: 대기, 해양, 산업용 터빈 및 차량 등 다양한 환경에서 발생하는 비압축성 유체 흐름을 모델링하는 것은 중요합니다. 특히, 계면활성제 (surfactant) 가 포함된 2 상 유동 (예: 움직이는 포획물 주변의 계면활성제 흐름) 은 물리적 복잡성 (계면 역학, 표면 대류 - 확산, 표면 장력 결합 등) 으로 인해 수치 해석이 매우 어렵습니다.
• 핵심 문제:
  • 움직이는 도메인: 물체의 운동이 시간에 따라 변할 때, 기존의 적합 격자 (fitted mesh) 방법은 매 시간 단계마다 격자를 재생성 (remeshing) 해야 하므로 계산 비용이 매우 큽니다.
  • 기하학적 오차: 이동하는 경계를 단순한 다각형으로 근사할 때 발생하는 기하학적 오차로 인해 수치 해의 정확도가 떨어질 수 있습니다.
  • 고차 정확도: 비선형 대류 항을 처리하고 높은 시간 정확도를 유지하면서도 안정성을 확보하는 것은 여전히 도전 과제입니다.

2. 제안된 방법론 (Methodology)

이 논문은 이동 도메인에서 비압축성 나비에 - 스토크스 방정식을 풀기 위한 새로운 수치 기법을 제안하며, 다음과 같은 핵심 요소들을 통합합니다.

가. 공간 이산화: Ghost Finite Element Method (Ghost-FEM)

• Unfitted Mesh: 물체의 형상과 무관한 직교 격자 (Cartesian mesh) 를 사용합니다. 물체 내부와 경계를 잘라내는 (cut) 셀은 '활성 셀 (active cell)'로, 그 외는 '비활성 셀'로 처리합니다.
• Ghost Point 접근법: 경계 조건을 물리적 경계 ($\Gamma$) 가 아닌 근사된 경계 ($\Gamma_h$) 에 부과할 때 발생하는 문제를 해결하기 위해, '그리드 복원 (snapping-back-to-grid)' 알고리즘을 사용하여 내부 노드를 '유령 노드 (ghost nodes)'로 변환하고 해를 확장합니다. 이는 작은 절단 셀 (small cut cells) 로 인한 조건수 악화를 방지하고 구현을 단순화합니다.
• 유한 요소 공간: 속도에 대해 $Q_2$ (이차), 압력에 대해 $Q_1$ (선형) 을 사용하는 $Q_2/Q_1$ 통합 FEM 을 사용하여 inf-sup 조건을 만족시키고 3 차 정확도 (속도) 를 달성합니다.

나. 경계 조건 처리 및 기하학적 오차 보정

• Shifted Boundary Method (SBM): 경계 조건을 실제 경계가 아닌 근사 경계 ($\Gamma_h$) 에 부과할 때 발생하는 기하학적 오차를 보정하기 위해 테일러 급수 전개를 사용합니다. 이를 통해 고차 정확도를 유지하면서도 경계 조건을 약하게 (weakly) 부과할 수 있습니다.
• Nitsche's Method: 디리클레 (Dirichlet) 경계 조건을 약하게 부과하기 위해 Nitsche 기법을 사용하며, 안정화 파라미터 ($\lambda$) 는 일반 고유값 문제 (generalized eigenvalue problem) 를 풀어 최적화하여 수치적 안정성을 보장합니다.

다. 시간 이산화: IMEX (Implicit-Explicit) Scheme

• 비선형성 처리: 나비에 - 스토크스 방정식의 비선형 대류 항은 명시적 (Explicit) 으로, 점성 및 압력 항은 암시적 (Implicit) 으로 처리하는 IMEX 방식을 채택합니다.
• 고차 정확도: 2 단계 (2 차 정확도) 및 4 단계 (3 차 정확도) IMEX Butcher 테이블을 사용하여 시간 적분 정확도를 높였습니다.
• 단일화 (Monolithic) 접근: 속도와 압력을 분리하지 않고 하나의 행렬 시스템으로 연립하여 풀며, saddle-point 문제의 불안정성을 해결하기 위해 작은 안정화 항 ($\epsilon I$) 을 추가합니다.

라. 이동 도메인 처리: Aslam Extrapolation

• 도메인이 이동할 때, 이전 시간 단계의 해를 새로운 활성 영역으로 확장하기 위해 **Aslam의 2 차 외삽법 (quadratic extrapolation)**을 유한 요소 프레임워크에 적용했습니다. 이는 유체 영역 밖의 영역으로 해를 매끄럽게 확장하여 시간 적분 시 필요한 정보를 제공합니다.

3. 주요 기여 (Key Contributions)

1. 고차 Ghost-FEM 과 IMEX 의 결합: 이동 도메인에서 비압축성 유동 문제를 해결하기 위해 Ghost-FEM 공간 이산화와 고차 IMEX 시간 적분을 성공적으로 결합한 최초의 연구 중 하나입니다.
2. 안정성 보장: 안정화 파라미터를 일반 고유값 문제를 통해 동적으로 결정하여, 작은 절단 셀이 존재하는 경우에도 수치적 안정성을 확보했습니다.
3. 고차 정확도 달성: Shifted Boundary Method 와 Aslam 외삽법을 통해 기하학적 오차와 시간 적분 오차를 모두 고차 (공간 3 차, 시간 2~3 차) 로 제어할 수 있음을 입증했습니다.
4. 유연한 구현: 복잡한 형상과 이동 경계를 가진 도메인에 대해 재격자화 없이 효율적으로 시뮬레이션할 수 있는 프레임워크를 제시했습니다.

4. 수치 실험 및 결과 (Results)

논문은 다양한 벤치마크 문제를 통해 제안된 방법의 정확성과 효율성을 검증했습니다.

• 정확도 검증 (Stationary & Moving Domains):
  • 고정된 도메인 (원, 타원, 꽃 모양) 에서 Kim-Moin 해를 사용하여 오차를 분석한 결과, 2 단계 IMEX 는 2 차, 4 단계 IMEX 는 3 차 시간 정확도를 보였습니다.
  • 회전하는 도메인 (이동하는 타원 및 꽃) 에 대해서도 동일한 차수의 수렴성을 확인했습니다.
• 원판 주위 유동 (Flow around a disk, Re=100):
  • 카르만 와류 (Kármán vortex street) 형성 및 항력/양력 계수 ($C_D, C_L$) 를 계산하여 기존 연구 결과와 비교했습니다. 수치 결과가 참고 문헌의 범위와 잘 일치하며, 안정적인 진동 거동을 보였습니다.
• 구동 캐비티 (Driven Cavity):
  • 꽃 모양의 장애물이 있는 구동 캐비티 유동에서 레이놀즈 수 (1~1000) 에 따른 와류 구조 변화를 정확히 포착했습니다.
• 회전하는 꽃 모양 물체 및 진동하는 기포:
  • 회전하는 꽃 모양 물체 주위의 유동과 계면활성제 연구에 필요한 진동하는 기포 (타원형 변형 및 횡방향 진동) 주위의 유동을 시뮬레이션했습니다. 저 레이놀즈 수와 고 레이놀즈 수에서 각각 다른 와류 형성 패턴을 관찰했습니다.

5. 의의 및 결론 (Significance & Conclusion)

• 기술적 의의: 이 연구는 이동 경계 문제를 해결하는 데 있어 재격자화의 비효율성을 극복하고, 고차 정확도를 유지하면서도 구현이 간단한 강력한 수치 기법을 제시했습니다. 특히, Ghost-FEM 과 IMEX, Shifted Boundary Method 의 시너지는 복잡한 유체 - 구조 상호작용 (FSI) 문제 해결에 중요한 토대가 됩니다.
• 응용 가능성: 현재 연구는 물체의 운동이 유체에 미치는 영향 (일방향 결합) 만을 다루지만, 향후 **양방향 결합 (two-way coupling)**을 도입하여 진동하는 기포가 유체 흐름에 의해 다시 영향을 받는 상황을 시뮬레이션할 수 있게 됩니다. 이는 계면활성제 흡착 역학 연구 및 미세 유체 공학 분야에서 중요한 응용 가능성을 열어줍니다.
• 소프트웨어: 모든 실험은 Julia 기반의 오픈소스 라이브러리인 Gridap과 GridapEmbedded 플러그인을 사용하여 구현되었습니다.

요약하자면, 이 논문은 이동하는 경계를 가진 비압축성 유동 문제를 해결하기 위한 고차 정확도, 안정적이며 효율적인 수치 프레임워크를 제시하며, 복잡한 유체 역학 현상을 모델링하는 데 있어 중요한 진전을 이루었습니다.