A hybrid Volume of Fluid Phase-Field method for Direct Numerical Simulations of soluble surfactant-laden interfacial flows

본 논문은 용해성 계면활성제가 포함된 유동의 직접 수치 시뮬레이션을 위해 적응형 메쉬 세분화를 적용한 하이브리드 부피-유체 위상장 방법을 제시하며, 이는 벌크 및 계면 수송 간의 결합을 정확하게 포착하여 마랑고니 응력이 3 차원 기하학적 구조에서 기포 상승 역학을 어떻게 크게 변화시키는지 입증한다.

핵심

이 논문은 간단한 언어와 창의적인 비유를 사용하여 설명합니다.

큰 그림: "비눗물"을 시뮬레이션하는 새로운 방법

물이 담긴 유리잔 속으로 기포가 떠오르는 모습을 상상해 보세요. 물이 순수하다면 기포는 약간 흔들리지만 빠르게 곧바로 위로 솟아오릅니다. 하지만 물에 비누 (계면활성제) 를 넣으면 기포의 행동이 달라집니다. 속도가 느려지거나, 더 많이 흔들리거나, 심지어 경로를 바꾸기도 합니다.

이런 현상이 일어나는 이유는 비누 분자들이 기포 표면에 달라붙기를 좋아하기 때문입니다. 기포가 움직이면 이 분자들이 밀려다니면서 기포 막에 불균일한 "장력"을 만들어냅니다. 이 불균일한 장력은 기포를 다양한 방향으로 밀어붙이는 보이지 않는 손처럼 작용하여 기포의 움직임을 변화시킵니다.

문제점:
이를 컴퓨터로 시뮬레이션하는 것은 매우 어렵습니다. 마치 두 가지 상충되는 설정을 가진 카메라로 비눗방울을 촬영하려는 것과 같습니다.

1. 선명한 카메라: 압력과 모양을 계산하기 위해 기포의 가장자리를 칼날처럼 얇은 선으로 정확하게 보아야 합니다.
2. 흐릿한 카메라: 비누 분자가 그 가장자리를 따라 매끄럽게 퍼지는 모습을 보아야 비누의 움직임을 계산할 수 있습니다.

대부분의 컴퓨터 방법은 두 설정 중 하나를 선택하도록 강요합니다. 그 결과 시뮬레이션은 물리적으로 부정확해지거나, 복잡한 3 차원 모양을 실행하는 것이 계산상 불가능해집니다.

해결책:
이 논문의 저자들은 하이브리드 방법을 개발했습니다. 마치 두 가지 세계의 장점을 동시에 활용하는 "스플릿 스크린" 시뮬레이션과 같습니다.

• 선명한 가장자리 (부피-유체법): 기포의 가장자리를 선명하게 유지하고 액체의 양을 완벽하게 보존하는 방법을 사용합니다 (고화질 윤곽선과 같습니다).
• 매끄러운 비누 (위상장법): 비누 분자가 이동할 수 있는 매끄러운 고속도로 역할을 하는 두 번째 "흐릿한" 레이어를 사용합니다. 이를 통해 비누가 물과 기포 표면 사이에서 갇히거나 사라지지 않고 자연스럽게 이동할 수 있습니다.

작동 원리: "교통 통제관" 비유

이를 구현하기 위해 저자들은 비누 분자를 위한 디지털 교통 시스템을 만들었습니다.

1. 고속도로 (인터페이스): 기포 표면은 붐비는 고속도로입니다. 비누 분자는 자동차입니다.
2. 램프 (흡착/탈착):
   • 흡착: 물 (본체) 에서 온 비누 분자들이 고속도로 (기포 표면) 로 뛰어오르고자 합니다.
   • 탈착: 비누 분자들이 지쳐서 고속도로를 떠나 다시 물로 뛰어내립니다.
   • 새로운 방법은 매 순간 몇 대의 자동차가 오르거나 내리는지 정확히 계산하여, 자동차가 사라지거나 갑자기凭空 나타나는 일이 없도록 합니다.
3. 교통 체증 (마랑고니 응력): 기포의 한 곳에 비누 자동차들이 너무 많이 몰리면 그 부분이 "끈적해집니다" (높은 장력). 기포 막은 그 끈적한 부분에서 벗어나려고 당기는데, 이로 인해 기포가 느려지거나 흔들리는 힘이 발생합니다. 이 시뮬레이션은 이러한 줄다리기 현상을 완벽하게 포착합니다.

테스트 내용 ("운전 시험")

새로운 차를 고속도로에 내보내기 전에, 저자들은 세 가지 특정 시험으로 운전 학교를 통과시켰습니다.

1. 신장 시험 (팽창하는 구): 비누가 덮인 기포를 불어 올렸습니다. 기포가 커짐에 따라 비누가 고르게 퍼지는지 확인했습니다. 시뮬레이션은 수학적 계산과 완벽하게 일치했습니다.
2. 회전 시험 (회전하는 기포): 비누가 붙은 기포를 회전시켰습니다. 비누가 원 주변을 올바르게 이동하며 새지 않는지 확인했습니다. 역시 시뮬레이션은 정확했습니다.
3. 교환 시험 (평평한 벽): 비누가 물에서 평평한 벽으로 이동했다가 다시 돌아오는 모습을 관찰했습니다. 세 가지 시나리오를 테스트했습니다.
   • 오직 붙기만 함: 비누가 달라붙는가? 네.
   • 오직 떨어지기만 함: 비누가 떠나는가? 네.
   • 둘 다: 균형이 잡히는가? 네.

메인 이벤트: 떠오르는 기포

마지막으로, 저자들은 새로운 방법으로 3 차원 물탱크 속을 떠오르는 기포를 시뮬레이션했습니다.

• 깨끗한 기포: 상대적으로 빠르고 곧게 떠올랐습니다.
• "불용성" 비누 기포: 비누가 표면에 고정되어 떠날 수 없었습니다. 이로 인해 기포 뒤쪽에 강력한 "교통 체증"이 생겨 기포가 크게 느려졌습니다.
• "가용성" 비누 기포 (실제 상황): 여기서 새로운 방법이 빛을 발합니다. 비누가 기포에 오르내릴 수 있었습니다.
  • 비누가 쉽게 떨어졌다면 (높은 "탈착"), 기포는 깨끗한 기포와 거의 비슷하게 행동했습니다.
  • 비누가 쉽게 붙었다면 (높은 "흡착"), 기포는 비누가 고정된 버전처럼 행동했습니다.
  • 중간 상태에서는 기포가 복잡한 춤을 추었습니다: 속도가 느려지고, 경로를 바꾸며, 떠오르는 동안 물속에 비누의 "흔적"을 남겼습니다.

왜 중요한가 (논문에 따르면)

저자들은 이 방법이 견고하고, 확장 가능하며, 정확하다고 주장합니다.

• 견고함: 기포가 이상한 모양이 되거나 부서져도 시스템이 충돌하지 않습니다.
• 확장성: 슈퍼컴퓨터에서 실행하여 거대하고 복잡한 3 차원 시뮬레이션을 효율적으로 처리할 수 있습니다.
• 정확성: 기포가 얼마나 빠르게 떠오르고 어떻게 흔들리는지 정확하게 예측하여 실제 물리 현상과 일치합니다.

간단히 말해: 그들은 마침내 3 차원 공간에서 비눗방울이 어떻게 행동하는지 시뮬레이션할 수 있는 새로운 디지털 엔진을 구축했습니다. 기포의 모양과 기포 막 위를 오르내리는 비누 분자 사이의 까다로운 춤을 처리하면서도 정확도를 잃거나 컴퓨터를 충돌시키지 않습니다.

기술 요약

기술 요약: 용해성 계면활성제가 포함된 계면 유동을 위한 하이브리드 부피-유체 및 위상장 방법

문제 제기
용해성 계면활성제가 포함된 계면 유동의 정확한 수치 시뮬레이션은 문제의 다중 규모 및 다중 물리학적 특성으로 인해 formidable한 도전 과제를 제시합니다. 불용성 계면활성과 달리, 용해성 종은 벌크 대류 - 확산, 계면 수송, 그리고 동적 교환 (흡착/탈착) 과정을 포함하는 결합된 시스템에 의해 지배됩니다. 이러한 상호작용은 공간적 및 시간적으로 변화하는 표면 장력장을 생성하여, 마랑고니 응력을 유발하고 유동 위상, 안정성 및 운동을 크게 변화시킵니다. 기존 수치 방법들은 다음과 같은 명확한 한계에 직면해 있습니다:

• 경계 적분/요소 방법: Stokes 유동 영역으로 제한되며, 위상 변화 (분열/병합) 에 어려움을 겪습니다.
• 프론트 추적 (Front-Tracking): 높은 정확도를 제공하지만 3 차원에서의 구현이 복잡하며, 명시적인 위상 처리가 필요하고 병렬 컴퓨팅에서 부하 균형 문제를 겪습니다.
• 레벨-셋 방법: 위상 변화를 자연스럽게 처리하지만, 재초기화 및 사후 보정이 필요한 질량 보존 문제를 종종 겪습니다.
• 위상장 방법: 위상 변화를 자연스럽게 처리하지만, 종종 강한 대류에 의해 위반될 수 있는 열역학적 평형 가정에 의존하여 질량 비보존을 초래합니다.
• 부피-유체 (VoF) 방법: 질량 보존과 날카로운 계면에 탁월하지만, 전통적으로 0 두께 계면에서의 계면활성제 수송에 필요한 표면 미분 연산자를 정의하는 데 어려움을 겪습니다.

적응형 메쉬 정제 (AMR) 를 지원하고, 동적 벌크 - 계면 교환을 갖는 용해성 계면활성제를 처리하며, 확장 가능한 오픈소스 프레임워크 내에서 작동하는 기하학적 VoF 기반 방법에 대한 특정 필요성이 존재합니다.

방법론
저자들은 오픈소스 코드 Basilisk에 구현된 하이브리드 수치 프레임워크를 제시하며, 이는 기하학적 부피 - 유체 (VoF) 방법과 확산 위상장 접근법을 결합합니다.

1. 계면 표현:

   • VoF ($c$): 한 유체의 부피 분율을 나타내는 스칼라 장으로, 보존적 계면 추적 및 운동량 해를 위해 사용됩니다. 이는 정확한 질량 보존을 보장하고 기하학적 재구성을 지원합니다.
   • 위상장 ($\phi$): Allen-Cahn 기반 정밀 보존 확산 계면 (ACDI) 형식에 따라 진화하는 보조 변수입니다. 이는 계면활성제 수송을 위한 매끄러운 운반체 역할을 하며, 날카로운 계면에서의 명시적 표면 미분 계산 없이 정규화된 법선 벡터, 곡률 및 표면 기울기를 제공합니다.
   • 결합: 위상장은 일관성을 유지하기 위해 부호화된 거리 함수를 사용하여 VoF 장에서 주기적으로 재초기화됩니다.

2. 지배 방정식:

   • 유체 운동: 공간적으로 변화하는 밀도와 점도를 갖는 비압축성 Navier-Stokes 방정식입니다. 표면 장력 힘은 체적 체적력으로 통합됩니다: 등방성 모세관 힘 ($\sigma \kappa \mathbf{n} \delta_s$) 및 접선 방향 마랑고니 힘 ($\nabla_s \sigma \delta_m$).
   • 계면활성제 수송: 두 개의 스칼라 장이 해결됩니다: 벌크 농도 ($F$) 및 계면 농도 ($f$, $\Gamma$의 체적 표현).
     • 계면 수송 방정식은 표면 미분을 피하기 위해 정규화되어, 표면 농도 $\Gamma$를 확산 계면 내에 국소화된 체적 장 $f = \Gamma \delta_\phi$로 변환합니다.
     • 동역학: 흡착 및 탈착은 벌크 및 계면 농도를 연결하는 계면에서 국소화된 정규화된 소스/싱크 항 ($J$) 을 통해 모델링됩니다. 이 모델은 표면 장력을 농도와 연관시키기 위해 일반적인 상태 방정식 (예: Langmuir 등온선) 을 지원합니다.

3. 수치 이산화:

   • 공간: 동적으로 정제된 4 분나무/8 분나무 격자 (AMR) 에서의 유한 체적 이산화입니다. 농도 및 유동의 급격한 기울기를 해결하기 위해 해상도가 계면 근처에 집중됩니다.
   • 시간: 연산자 분할이 사용됩니다. 대류는 명시적으로 처리됩니다 (기하학적 언스플릿 Bell–Collela–Glaz 방식), 반면 확산, 반확산 및 소스 항은 수정된 다중 그리드 Poisson 솔버를 사용하여 암시적으로 처리됩니다.
   • 안정성: 강한 계면 기울기 하에서의 안정성을 보장하기 위해 대류 및 반확산 기여를 모두 포함하는 유효 속도에 기반한 Courant–Friedrichs–Lewy (CFL) 조건이 강제됩니다.

주요 기여

• 하이브리드 프레임워크: VoF 의 질량 보존 및 날카로운 계면 장점을 유지하면서, 계면활성제 수송을 위한 매끄러운 운반체로만 위상장 변수를 활용하고 조각별 선형 계면에서의 기하학적 복잡성을 피하는 방법 개발.
• 용해성 계면활성제 능력: 불용성 계면활성제로 제한되었던 이전 VoF-위상장 작업을 완전히 용해성 사례로 확장하여, 결합된 벌크 대류 - 확산 및 동적 교환 항을 통합.
• 확장성 및 적응성: 이 프레임워크는 완전히 적응적 (AMR) 이며 2 차원, 축대칭 및 3 차원 구성을 지원하고 높은 병렬 확장성을 가지며, 벌크 및 계면 장 모두에서 다중 규모 기울기를 해결할 수 있습니다.
• 오픈소스 구현: 이 방법은 Basilisk 에 구현되어 커뮤니티에 재현 가능하고 접근 가능한 도구를 제공합니다.

결과 및 검증
이 논문은 일련의 벤치마크 테스트 및 적용 사례를 통해 방법을 검증합니다:

1. 유동 없는 검증:

   • 확장하는 구 (3D): 등방성 팽창 중 불용성 계면활성제 희석에 대한 해석적 해의 회복을 검증했습니다.
   • 회전하는 계면 (2D): 이동하는 계면에서의 불용성 계면활성제 대류 - 확산에 대한 2 차 수렴을 확인했습니다.
   • 평평한 계면에서의 흡착/탈착: 일정 플럭스 흡착, 순수 탈착 및 완전 Langmuir 동역학에 대한 해석적 해에 대해 결합된 벌크 - 계면 솔버를 검증했습니다. 이 방법은 평형으로의 전환과 벌크 내 고갈/농축 층의 형성을 정확하게 포착합니다.

2. 축대칭 상승 기포:

   • 불용성 한계: 마랑고니 응력으로 인한 종단 상승 속도의 감소를 시연했으며, 계면활성제 농도가 감소함에 따라 결과가 깨끗한 계면 한계로 수렴하는 것을 보였습니다.
   • 용해성 동역학: Biot 수 (탈착) 및 Damköhler 수 (흡착) 를 체계적으로 변화시켰습니다.
     • 높은 Biot (빠른 탈착): 계면활성제가 계면에서 빠르게 제거됨에 따라 시스템은 깨끗한 계면 한계에 접근합니다.
     • 높은 Damköhler (빠른 흡착): 계면이 빠르게 포화됨에 따라 시스템은 표면 비압축성 (불용성) 한계에 접근합니다.
   • 유동 구조: 시각화는 기포 후미에 계면활성제의 "정체 캡" 형성을 보여주었으며, 이는 내부 순환을 억제하고 와류 구조를 수정합니다.

3. 3 차원 상승 기포:

   • 용해성 계면활성제가 있는 3 차원 상승 기포의 시뮬레이션은 모양 변형, 측면 궤적 편차 (나선 운동) 및 탈착으로 인한 후미 플룸 형성 등 복잡한 현상을 포착했습니다.
   • 이 방법은 3 차원에서 깨끗한 및 불용성 점근 한계 모두를 성공적으로 회복하여, 비정상 유동 패턴 및 변형 기하학을 처리하는 데 있어 견고함을 입증했습니다.

의의 및 주장
이 논문은 고전적 불용성 계면활성제 모델과 현실적인 용해성 공식 사이의 간극을 연결하는 견고하고 다재다능한 수치 프레임워크를 제공한다고 주장합니다. 그 의의는 다음과 같습니다:

• 통합 처리: 단일 오픈소스 플랫폼 내에서 비선형 결합, 진화하는 계면 및 완전한 3 차원 구성을 처리합니다.
• 물리적 충실도: 난류 폐쇄 모델이나 경험적 근사에 의존하지 않고 유체역학, 벌크/계면 수송 및 마랑고니 응력 간의 복잡한 상호작용을 포착합니다.
• 실용적 적용성: 이 방법은 산업적 및 생물학적으로 관련된 용해성 계면활성제가 존재하는 계면 재연결 (예: 제트 분열, 파도 깨짐) 이 포함된 자유 표면 유동에 적합합니다.

저자들은 현재 프레임워크가 이산화의 고유한 수치 확산으로 인해 중간 Péclet 수로 제한된다고 겸손하게 지적합니다. 그들은 높은 Péclet (대류 지배) 영역에서 날카로운 계면 기울기를 정확하게 해결하는 것을 향후 과제로 식별하며, 수치 확산을 최소화하는 형식을 개발하는 것을 목표로 합니다.