cfdmfFTFoam: A front-tracking solver for multiphase flows on general unstructured grids in OpenFOAM

이 논문은 복잡한 알고리즘과 격자 제한으로 인해 부족했던 오픈소스 프론트-트래킹 솔버의 공백을 메우기 위해, 일반 비정렬 격자에서 작동하도록 OpenFOAM 에 통합된 새로운 솔버 'cfdmfFTFoam'을 개발하고 그 유효성을 검증한 내용을 다루고 있습니다.

핵심

🌟 핵심 비유: "유리잔 속의 비눗방울을 추적하는 마법사"

상상해 보세요. 유리잔 안에 물이 있고, 그 안에 비눗방울이 떠다니고 있습니다. 우리는 이 비눗방울이 어떻게 움직이고, 모양이 어떻게 변하는지 알고 싶습니다.

기존의 컴퓨터 프로그램들은 비눗방울을 **'투명한 페인트'**처럼 처리했습니다. (이것을 VOF 방법이라고 합니다.) 페인트가 번지듯 흐르기 때문에, 비눗방울이 아주 얇아지거나 찢어질 때 모양이 뭉개지거나 사라지는 오류가 자주 발생했습니다. 마치 물에 잉크를 떨어뜨려서 잉크의 경계를 정확히 따라가는 건 불가능한 것과 비슷합니다.

하지만 이 논문에서 소개된 cfdmfFTFoam은 완전히 다른 방식을 사용합니다.

1. "투명한 막" 대신 "실제 그물망"을 씌우다 (Front-Tracking Method)

이 프로그램은 비눗방울의 표면을 **실제 삼각형 조각들로 만든 유연한 그물망 (Front Mesh)**으로 감싸서 추적합니다.

• 비유: 비눗방울 위에 미세한 비즈로 만든 그물망을 씌운다고 상상해 보세요. 비눗방울이 움직일 때, 그 그물망도 함께 움직이며 비눗방울의 모양을 정확히 따라갑니다.
• 장점: 그물망이 비눗방울의 '진짜 경계'를 정확히 잡기 때문에, 비눗방울이 찢어지거나 합쳐질 때의 모양을 훨씬 더 선명하고 정확하게 보여줍니다.

2. "모든 모양의 유리잔"에 적용 가능 (General Unstructured Grids)

기존의 이 기술 (Front-Tracking) 은 컴퓨터가 계산하기 쉬운 **정사각형 격자 (Structured Grid)**에서만 작동했습니다. 마치 정사각형 타일만으로만 만든 방에서만 비눗방울을 추적할 수 있었던 거죠. 하지만 현실의 배관이나 기계는 구부러지고 복잡한 모양입니다.

• 혁신: 이 연구팀은 그물망 기술을 **어떤 모양의 유리잔 (불규칙한 격자)**에서도 작동하도록 개조했습니다. 마치 정사각형 타일뿐만 아니라, 삼각형, 육각형 등 모든 모양의 타일로 만든 복잡한 방에서도 비눗방울을 정확히 추적할 수 있게 된 것입니다.

3. "그물망"과 "물"의 대화 (Communication)

그물망 (비눗방울 표면) 은 물 (유체) 과 끊임없이 상호작용합니다.

• 상황: 비눗방울이 물속을 움직이면 물이 밀려나고, 물의 흐름이 비눗방울을 밀어냅니다.
• 해결책: 이 프로그램은 그물망의 위치 정보를 물의 흐름 계산에 전달하고, 물의 흐름 정보를 다시 그물망에 전달하는 **정교한 통역사 (RKPM 알고리즘)**를 개발했습니다. 두 세계가 오해 없이 소통하게 해주는 것이죠.

4. "그물망"이 찢어지지 않게 수리하기 (Remeshing & Volume Correction)

비눗방울이 심하게 늘어나거나 구부러지면 그물망이 찢어지거나 구멍이 날 수 있습니다.

• 수리공: 이 프로그램에는 자동 수리공이 붙어 있습니다. 그물망이 너무 늘어났으면 잘라내고 (Refining), 너무 뭉개졌으면 합쳐서 (Coarsening) 모양을 다듬습니다.
• 체중 조절: 비눗방울의 부피가 줄어들거나 늘어나면 안 되는데, 계산 오차 때문에 부피가 변할 수 있습니다. 이때 **체중 조절기 (Volume Correction)**가 작동하여 부피를 원래대로 되돌려줍니다.

---

📊 이 프로그램이 왜 중요한가요? (결과 요약)

연구진은 이 프로그램이 얼마나 잘 작동하는지 여러 가지 실험으로 증명했습니다.

1. 비눗방울 변형 실험: 비눗방울을 강하게 늘려서 다시 원래 모양으로 돌아오게 했을 때, 기존 프로그램은 비눗방울이 찢어지거나 사라졌지만, 이 프로그램은 완벽하게 원래 모양을 유지했습니다.
2. 기포 상승 실험: 물속에서 기포가 올라가는 실험을 했을 때, 이 프로그램은 실험실에서 본 실제 기포의 모양과 속도와 거의 일치하는 결과를 보여주었습니다. 특히 물과 기름처럼 밀도가 아주 다른 액체가 섞일 때에도 잘 작동했습니다.

🚀 결론: 미래의 열쇠

이 cfdmfFTFoam 프로그램은 **오픈 소스 (누구나 무료로 사용할 수 있음)**로 공개되었습니다.

• 의미: 이제 전 세계의 연구자들이 이 도구를 이용해 약품 제조, 석유 추출, 엔진 연소, 심지어 인체 내 혈류 등 다양한 분야에서 액체와 기체의 복잡한 상호작용을 더 정밀하게 연구할 수 있게 되었습니다.

한 줄 요약:

"이 프로그램은 컴퓨터 안에서 비눗방울의 표면을 '그물망'으로 감싸서, 어떤 복잡한 모양의 용기 안에서도 비눗방울이 찢어지거나 사라지지 않고 정교하게 움직이는 모습을 완벽하게 재현하는 마법 같은 도구입니다."

기술 요약

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

• 배경: 다상 유동 (Multiphase flow) 의 정확한 예측을 위해 인터페이스를 명시적으로 추적하는 '전면 추적법 (Front-Tracking Method, FTM)'은 인터페이스의 파열/합체 물리, 정확한 법선 벡터 및 곡률 계산, 인터페이스 확산 오류 회피 등의 장점을 가집니다.
• 문제점:
  • 기존 오픈소스 FTM 솔버 (PARIS, FronTier++, TrioCFD, Basilisk 등) 는 대부분 **구조화된 직교 격자 (Structured Cartesian grids)**에만 국한되어 있습니다.
  • 복잡한 형상을 다루기 위해 필수적인 **일반 비정렬 격자 (General Unstructured grids)**를 지원하는 순수 FTM 솔버는 부재합니다.
  • OpenFOAM 기반의 기존 솔버 (lentFoam 등) 는 FTM 과 레벨셋 (Level-Set) 을 혼합한 하이브리드 방식이거나, 인터페이스를 순수한 라그랑지안 메쉬로 표현하지 않아 FTM 의 고유한 장점을 완전히 활용하지 못합니다.
• 목표: OpenFOAM 환경에서 일반 비정렬 격자를 지원하면서도 순수 FTM 알고리즘을 구현하여, 복잡한 다상 유동 문제를 해결할 수 있는 오픈소스 솔버를 개발하는 것입니다.

2. 방법론 (Methodology)

• 솔버 개발: 기존 FTM 코드인 Ftc3D(Gretar Tryggvason 등 개발) 를 OpenFOAM v9 기반의 cfdmfFTFoam 솔버로 통합했습니다.
• 핵심 알고리즘 및 구현:
  1. 격자 간 통신 (Front-to-Field Communication): 라그랑지안 전면 메쉬 (Front) 와 오일러 격자 (Eulerian grid) 간의 데이터 전달을 위해 **재현 커널 입자 방법 (Reproducing Kernel Particle Method, RKPM)**을 구현했습니다. 이를 통해 표면 장력을 격자 노드로 매핑하고, 밀도/점도 필드를 업데이트합니다.
  2. 부피 보정 (Volume Correction): FTM 의 고유한 단점인 부피 (질량) 비보존 문제를 해결하기 위해 VC2 알고리즘을 적용하여 전면 정점의 위치를 보정합니다.
  3. 재메싱 (Remeshing): 전면 메쉬의 품질을 유지하기 위해 세분화 (Refining), 축소 (Coarsening), 매끄러움 (Smoothing/Undulation Removal) 알고리즘 (TSUR3D, VCS 등) 을 구현했습니다.
  4. 표면 장력 계산: 요소 기반 직접 계산법 (Direct Element-Based Method) 을 사용하여 인터페이스에 작용하는 표면 장력을 계산합니다.
  5. 지시자 함수 (Indicator Function) 구성: Poisson 방정식 (PE) 해법과 최단 거리 변환 (CPT) 방법을 지원하며, 비정렬 격자에서는 PE 방법이 더 우수한 성능을 보입니다.
  6. 병렬 처리: OpenFOAM 의 도메인 분해 (Domain Decomposition) 아키텍처를 지원하여 병렬 연산이 가능합니다.

3. 주요 기여 (Key Contributions)

• 최초의 오픈소스 순수 FTM 솔버: 일반 비정렬 격자 (Hexahedral, Tetrahedral, Polyhedral 등) 를 지원하는 최초의 오픈소스 FTM 솔버 (cfdmfFTFoam) 를 공개했습니다.
• 알고리즘 통합 및 개선: Ftc3D 의 핵심 로직을 OpenFOAM 프레임워크에 완벽하게 통합하고, RKPM 기반의 통신, 부피 보정, 다양한 재메싱 기법 등을 추가하여 솔버의 안정성과 정확도를 높였습니다.
• 검증 데이터 제공: 다양한 벤치마크 테스트를 통해 솔버의 정확성을 입증하고, 사용자를 위한 튜토리얼 케이스 및 소스 코드 (GitHub, CPC Library) 를 공개했습니다.

4. 검증 결과 (Results)

논문에서는 5 가지 주요 벤치마크를 통해 솔버를 검증했습니다:

1. 3 차원 변형 벤치마크 (3D Deformation): Leveque 의 유속장에서 드롭렛의 변형을 시뮬레이션. FTM 은 VOF(VOF-MULES, isoAdvector) 에 비해 인터페이스의 불필요한 파열 없이 큰 변형을 정확하게 추적하며, 부피 보존 오차가 1% 미만으로 매우 낮음을 확인했습니다.
2. 정지 드롭렛 벤치마크 (Stagnant Droplet): 정적 상태의 드롭렛에서 표면 장력과 압력 차이를 검증. cfdmfFTFoam 은 VOF 솔버에 비해 압력 오차 (0.046% vs 12.6%) 와 기생 전류 (Parasitic Current) 가 현저히 낮음을 보였습니다.
3. Poiseuille 유동 내 드롭렛 (Stokes Regime): 다양한 격자 형태 (육면체, 프리즘, 다면체, 사면체) 에서의 드롭렛 변형을 시뮬레이션. 해석해와 잘 일치하며, 4 개 프로세서에서 최적의 병렬 효율을 보였습니다.
4. 테일러 드롭렛 상승 (Taylor Droplet Rising): 실험 데이터 (Hayashi et al.) 와 비교. 초기 가속도 단계에서 발생하는 복잡한 유동 (후면의 공동 형성) 을 재현하기 위해 언듈레이션 제거 (Undulation removal) 빈도 조정이 필요함을 보였으며, 최종 형태가 실험과 잘 일치했습니다.
5. 자유 상승 기포 (Free Rising Bubble): 높은 밀도/점도 비율 ($10^3 \sim 10^4$) 과 큰 관성 효과를 가진 복잡한 기포 상승 시뮬레이션. 상승 속도와 기포 형태를 정확히 예측했으나, 매우 큰 물성비에서는 운동량 이산화의 불일치로 인한 오차가 발생하여 향후 개선이 필요함을 지적했습니다.

5. 의의 및 결론 (Significance)

• 연구 도구 제공: 복잡한 다상 유동 현상 (합체, 파열, 접촉선 동역학 등) 을 연구하는 연구자들에게 일반 비정렬 격자를 지원하는 강력한 오픈소스 도구를 제공합니다.
• 알고리즘 발전의 기반: cfdmfFTFoam 은 향후 FTM 알고리즘의 개선 (접촉선 동역학, 위상 변화, 열/물질 전달 등) 및 하이브리드 전략 개발을 위한 기반 솔버로 활용될 수 있습니다.
• OpenFOAM 생태계 확장: OpenFOAM 이 기존에 주로 VOF 나 Level-Set 에 의존하던 다상 유동 해법에서, 정확도가 높은 순수 FTM 기반 해법을 지원할 수 있게 함으로써 CFD 소프트웨어의 범위를 확장했습니다.

이 솔버는 GitHub와 CPC Library를 통해 공개되어 있으며, 다상 유동 연구의 정확도와 효율성을 높이는 데 기여할 것으로 기대됩니다.