A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries

이 논문은 임의의 복잡한 3 차원 경계면에서 접촉선 역학과 히스테리시스를 정밀하게 모사하기 위해, 혼합 셀에서의 엄격한 질량 보존을 보장하고 작은 절단 셀로 인한 시간 간격 제한을 해결하며 3 차원 접촉각 조건을 정확히 부과하는 새로운 기하학적 보존적 VOF 방법을 제안합니다.

핵심

이 논문은 **"복잡한 모양의 물체 위를 흐르는 물방울을 컴퓨터로 얼마나 정교하게 시뮬레이션할 수 있을까?"**라는 질문에 대한 해답을 제시합니다.

연구진 (황충선, 한천양 등) 은 물방울이 구부러진 벽, 거친 표면, 혹은 구멍이 뚫린 판 위를 이동할 때 발생하는 **'접촉선 (물방울 가장자리)'**의 움직임을 매우 정확하게 계산하는 새로운 3 차원 컴퓨터 프로그램을 개발했습니다.

이 복잡한 과학적 내용을 일상적인 비유로 쉽게 설명해 드리겠습니다.

---

1. 핵심 문제: "물방울이 구석진 곳에 갇히면 어떻게 될까?"

컴퓨터 시뮬레이션은 보통 정사각형 격자 (블록) 를 쌓아 세상을 표현합니다.

• 평평한 바닥: 물방울이 평평한 바닥에 있으면 블록 모양과 잘 맞습니다.
• 복잡한 모양: 하지만 물방울이 구멍이 뚫린 판이나 구부러진 구 위를 움직일 때, 정사각형 블록과 물체의 경계가 부딪히면 문제가 생깁니다. 마치 레고 블록으로 둥근 공을 만들려고 할 때 생기는 틈처럼, 물이不该 들어갈 곳 (고체 부분) 으로 새어 들어갈 수도 있고, 반대로 물이 있어야 할 공간이 비어있을 수도 있습니다.

기존 프로그램들은 이 '틈새'를 처리할 때 물의 양이 조금씩 사라지거나 늘어나는 오류를 범했습니다. 마치 물을 옮길 때 컵 가장자리를 스쳐서 물이 조금씩 새는 것과 같습니다.

2. 이 연구의 해결책 1: "물방울을 옮기는 새로운 규칙" (보존적 이동)

이 연구진은 **"물이 새지 않도록 완벽하게 막는 새로운 이동 규칙"**을 만들었습니다.

• 비유: "물통을 옮길 때, 구멍이 난 바닥을 고려하다"
  기존 방법은 물이 있는 공간만 보고 물을 옮겼는데, 이 새로운 방법은 **"아, 이 공간에는 고체 (벽) 가 있어서 물이 들어갈 수 있는 공간이 줄어든구나"**라고 계산합니다.
  • 적응형 이동: 만약 물이 들어갈 공간이 매우 좁다면 (작은 조각), 물을 옮기는 속도를 조절하거나, 넘쳐난 물을 옆집에 재분배하는 방식을 썼습니다.
  • 효과: 이렇게 하면 물이 사라지거나 불어나는 일이 전혀 없습니다. 마치 물방울을 옮길 때 물이 한 방울도 새지 않도록 완벽하게 밀봉된 통을 사용하는 것과 같습니다. 또한, 작은 공간 때문에 시뮬레이션 속도가 느려지는 문제도 해결했습니다.

3. 이 연구의 해결책 2: "물방울의 각도를 정확히 맞추는 기술" (접촉각)

물방울이 벽에 닿을 때, 벽을 타고 올라가는 각도 (접촉각) 는 물방울의 모양을 결정하는 핵심입니다.

• 기존의 한계: 평평한 벽에서는 각도를 맞추기 쉬웠지만, 구부러진 벽이나 거친 표면에서는 컴퓨터가 "어느 각도로 올라가야 하지?"를 계산할 때 엉뚱한 값을 내뱉었습니다. 마치 구부러진 경사로에서 자동차가 미끄러지는 각도를 계산할 때, 평지 공식을 적용하는 오류와 비슷합니다.
• 새로운 기술 (포물선 피팅): 연구진은 물방울 가장자리의 모양을 단순한 직선이 아닌, 부드러운 포물선 (파라볼라) 곡선으로 추정하는 새로운 수학적 기법을 도입했습니다.
  • 비유: 그림을 그릴 때, 직선 자 대신 유연한 곡선 자를 사용하여 벽의 굴곡에 딱 맞게 물방울의 모양을 그리는 것입니다.
  • 히스테리시스 (기억력): 물방울이 미끄러질 때, 앞쪽은 넓어지고 뒤쪽은 좁아지는 현상 (접촉각 히스테리시스) 을 모방했습니다. 마치 끈적이는 테이프를 떼어낼 때, 붙어있는 힘과 떨어지는 힘이 다르듯이, 물방울이 벽을 미끄러질 때의 '붙어있음'과 '떨어짐'을 정교하게 계산합니다.

4. 검증된 결과: "어떤 모양에서도 완벽하게 작동"

이 새로운 프로그램은 다양한 난이도의 테스트를 통과했습니다.

1. 회전하는 물방울: 구멍이 뚫린 구형 껍질 안에서 물이 회전할 때, 물의 양이 절대 변하지 않았습니다.
2. 구멍이 뚫린 판에 떨어지는 물방울: 물방울이 구멍을 통과할 때, 구멍 가장자리에 물이 걸리거나 찢어지는 현상을 실제 실험 결과와 똑같이 재현했습니다.
3. 뾰족한 원뿔을 타고 오르는 물방울: 물방울이 원뿔 모양의 표면을 타고 스스로 올라가는 현상을 정확히 시뮬레이션했습니다.
4. 메쉬 (망) 위를 퍼지는 물방울: 복잡한 격자 모양의 표면을 따라 물방울이 퍼져나가는 모습을 정교하게 그렸습니다.

5. 결론: 왜 이것이 중요한가?

이 연구는 **"복잡한 3 차원 세상에서 물방울의 움직임을 가장 정확하게, 그리고 물이 새지 않도록 시뮬레이션하는 방법"**을 처음 제시했습니다.

• 실제 활용: 잉크젯 프린터의 잉크 방울, 미세 유체 칩, 나노 기술, 심지어 우주에서의 연료 관리 등 정교한 모양의 표면과 액체가 만나는 모든 분야에 이 기술이 적용될 수 있습니다.
• 핵심 메시지: 이제 컴퓨터는 더 이상 평평한 바닥만 다룰 수 있는 것이 아니라, 구불구불한 길, 구멍, 거친 표면 위에서도 물방울이 어떻게 움직일지 마치 실제 눈으로 보는 것처럼 정확하게 예측할 수 있게 되었습니다.

이 논문은 복잡한 수학 공식 뒤에 숨겨진 **"물방울의 정확한 춤"**을 컴퓨터가 완벽하게 따라 할 수 있게 만든 획기적인 기술입니다.

기술 요약

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

• 배경: 액체/기체/고체의 삼중 접선 (Triple junction) 이 움직이는 접촉선 (Contact line) 동역학은 미세유체, 잉크젯 프린팅, 박막 제어 등 다양한 공학 및 과학 분야에서 핵심적인 역할을 합니다.
• 문제점: 기존의 접촉선 모델링 연구는 주로 평탄한 정사각형 격자 (Regular boundary) 에 국한되어 있었습니다. 그러나 실제 응용에서는 곡면이나 임의의 복잡한 형상을 가진 고체 경계 (Complex geometric boundaries) 에서의 거동을 시뮬레이션해야 하는 요구가 증가하고 있습니다.
• 주요 난제:
  1. 혼합 셀 (Mixed Cells) 내의 질량 보존: 복잡한 고체 경계가 격자를 절단할 때 생성되는 '혼합 셀' (유체와 고체가 공존하는 셀) 에서 인터페이스를 이동시킬 때 국소 및 전역 질량 보존을 엄격하게 유지하는 것이 어렵습니다. 특히 작은 절단 셀 (Small cut cells) 로 인해 시간 간격 (Time-step) 제한 (CFL 조건) 이 매우 엄격해져 계산 효율이 떨어집니다.
  2. 복잡한 경계면에서의 접촉각 조건 부과: 평탄한 표면에서는 Height Function(HF) 방법을 통해 접촉각을 비교적 쉽게 부과할 수 있지만, 곡면이나 임의의 경계면에서는 접촉각 조건을 정확하게 적용하고 곡률을 정밀하게 추정하는 것이 어렵습니다. 기존 선형 보간법은 극단적인 접촉각이나 복잡한 기하학적 구조에서 큰 오차를 발생시킵니다.

2. 방법론 (Methodology)

이 논문은 Basilisk 오픈소스 솔버를 기반으로 하여, 기하학적 Volume-of-Fluid (VOF) 방법과 임베디드 경계 방법 (Embedded Boundary Method, EBM) 을 결합한 3 차원 수치 프레임워크를 제안합니다.

A. 혼합 셀을 위한 보존적 VOF 이송 (Conservative VOF Advection)

• 수정된 이송 및 재구성: 액체, 기체, 고체가 공존하는 혼합 셀 내에서 인터페이스를 재구성하기 위해 3 차원 Flood-fill 알고리즘을 개발했습니다. 이는 불규칙한 다면체 영역 내에서 액체 부피를 정확히 계산하고 PLIC(Piecewise Linear Interface Calculation) 방식을 적용합니다.
• 플럭스 보정 (Flux Correction): 기존 기하학적 이송 방식이 고체 경계로 인해 유체 부피가 차단되는 효과를 무시하여 발생하는 질량 오차를 보정하기 위해, 수정된 이송 폭 ($u_n^*$) 을 도입했습니다. 이는 고체로 막힌 부피를 가상으로 보정하여 국소 질량 보존을 엄격하게 달성합니다.
• 재분배 이송 전략 (Redistribution Advection Strategy): 작은 절단 셀로 인한 CFL 시간 간격 제한을 완전히 제거하기 위해 플럭스 재분배 (Flux Redistribution) 기법을 도입했습니다. 이 방법은 인터페이스 셀에서 발생하는 부피 과잉 또는 부족분을 하류 이웃 셀로 보존적으로 재분배하여, 시간 간격 제한 없이도 국소 및 전역 부피 보존을 유지합니다.

B. Height Function 기반의 3 차원 접촉각 부과 기법

• 포물면 피팅 (Paraboloid Fitting): 평탄한 표면에서의 선형 보간법을 넘어, 복잡한 경계면에서 접촉각 조건을 정확히 부과하기 위해 3 차원 포물면 (Paraboloid) 피팅 방식을 개발했습니다.
  • 예측 피팅 (Pre-fitting): 접촉선 셀 주변의 유효한 HF(Height Function) 데이터 포인트를 선별하기 위해 1 단계로 포물면을 미리 피팅하여 기하학적 일관성을 확보합니다.
  • 제약 피팅: 접촉각 조건 ($\theta$) 을 명시적으로 만족하도록 최종 포물면을 피팅하여 곡률 추정과 접촉각 부과를 동시에 수행합니다.
• 수직 Height Function (Vertical HF): 접촉각이 매우 작거나 ($<45^\circ$) 큰 ($>135^\circ$) 경우나 고체 표면이 수평면과 거의 평행한 경우, 기존 수평 HF 방법이 실패하는 문제를 해결하기 위해 수직 HF 방법을 도입하여 임의의 기하학적 구조에서도 접촉각을 부과할 수 있도록 했습니다.
• 이력 현상 (Hysteresis) 모델링: 접촉각 이력 현상 (접촉각이 전진각과 후퇴각 사이에서 변화하며 접촉선이 고정되는 현상) 을 모델링하기 위해, 부피 보존 원리에 기반한 이터레이션 (Bisection method) 을 통해 각 시간 단계에서 접촉각을 동적으로 결정하는 알고리즘을 구현했습니다.

3. 주요 기여 (Key Contributions)

1. 3 차원 보존적 인터페이스 이송 체계: 임베디드 경계가 있는 3 차원 공간에서 혼합 셀 내의 부피 이송을 정확히 계산하고, 작은 절단 셀로 인한 시간 간격 제한을 제거하면서도 국소 및 전역 질량 보존을 엄격하게 유지하는 완전 기하학적 (Fully geometric) 알고리즘을 최초로 제안했습니다.
2. 3 차원 Height Function 기반 접촉각 부과: 임의의 복잡한 3 차원 고체 표면에서 접촉각 조건을 정확하게 부과하고, 접촉각 이력 현상을 모델링할 수 있는 새로운 알고리즘을 개발했습니다. 이는 기존 선형 보간법의 오차를 해결하고 곡률 추정의 안정성을 크게 향상시켰습니다.

4. 검증 및 결과 (Results)

제안된 방법의 정확성, 견고성, 보존성을 검증하기 위해 다양한 벤치마크 테스트를 수행했습니다.

• 부피 보존 및 시간 간격 제한:
  • 회전하는 액체 껍질 (Concentric sphere) 시뮬레이션에서 기존 방식은 국소 질량 오차가 발생했으나, 제안된 보정 방식은 기계 정밀도 (Machine precision) 수준의 오차 ($L_\infty \approx 10^{-15}$) 를 달성했습니다.
  • 원뿔형 컷아웃이 있는 회전 액체 껍질 테스트에서 재분배 전략은 작은 절단 셀로 인한 시간 간격 제한을 제거하면서도 기존 방식 대비 형상 오차를 10 배 이상 감소시켰습니다.
• 접촉각 정확도:
  • 평평한 판과 구형 표면에서의 액적 확산 테스트에서, 제안된 포물면 피팅 방식은 접촉각 ($\theta$) 과 고체 표면의 위치/방향 ($d$) 에 무관하게 이론값과 매우 잘 일치하는 평형 반경을 보여주었습니다. 반면, 기존 선형 피팅 방식은 표면 깊이 ($d$) 에 따라 오차가 크게 증가하고 접촉선 대칭성이 깨지는 현상을 보였습니다.
• 이력 현상 및 복잡한 기하학:
  • 전단 유동 하의 액적 운동 시뮬레이션에서 접촉각 이력 현상을 정확히 재현하여 정적 및 동적 핀닝 (Pinning) 및 디핀닝 (Depinning) 거동을 잘 묘사했습니다.
  • 사인파형 기저면, 구멍이 있는 판, 원뿔형 표면, 메쉬 구조 등 매우 복잡한 3 차원 기하학에서 액적의 충돌, 확산, 자발적 이동 (Spontaneous motion) 을 성공적으로 시뮬레이션하여 방법론의 범용성과 안정성을 입증했습니다.

5. 의의 및 결론 (Significance)

이 연구는 임의의 복잡한 3 차원 표면에서 이동 접촉선 동역학을 정확하게 해석할 수 있는 완전히 기하학적이고 보존적인 VOF 기반 프레임워크를 최초로 확립했다는 점에서 의의가 큽니다.

• 기술적 진보: 기존 방법들이 직면했던 '작은 절단 셀로 인한 계산 비효율성'과 '복잡한 경계면에서의 접촉각/곡률 추정 오차'라는 두 가지 근본적인 난제를 동시에 해결했습니다.
• 실용성: 이 방법은 미세유체 장치, 코팅 공정, 3D 프린팅 등 복잡한 형상의 고체 표면과 유체 상호작용이 중요한 다양한 공학 및 과학 분야에서 정밀한 시뮬레이션을 가능하게 하는 강력한 도구로 활용될 수 있습니다.
• 확장성: 오픈소스 솔버 (Basilisk) 에 통합되어 향후 다양한 3 차원 다상 유동 문제에 적용 및 확장될 수 있는 기반을 마련했습니다.

결론적으로, 본 논문은 3 차원 다상 유동 시뮬레이션 분야에서 접촉선 모델링의 정확도와 효율성을 획기적으로 향상시킨 중요한 연구 성과입니다.