Solid adsorption: the missing mechanism for surfactant contact lines -- a phase-field approach

기술 요약: 고체 흡착: 계면활성제 접촉선의 누락된 메커니즘 – 위상장(Phase-Field) 접근법

문제 제기
계면활성제는 계면 장력을 변화시키고 젖음성을 조절함으로써 다상 시스템에서 결정적인 역할을 한다. 다양한 시스템에서 관찰되는 견고한 실험적 사실은 계면활성제가 계면 장력을 체계적으로 감소시켜, 친수성 및 소수성 표면 모두를 젖음성이 증가하는 방향으로 이동시킨다는 것이다. 가용성 계면활성제를 포함하는 이동 접촉선 문제에 대한 기존의 수치적 프레임워크들은 근본적인 한계를 공유하고 있는데, 바로 고체 벽에서의 계면활성제 흡착을 간과한다는 점이다. 유체-유체 계면 장력 감소만을 고려하는 기존 모델들은 계면활성제가 단순히 영의 방정식(Young's equation)을 재조정(rescale)할 뿐이라고 예측한다. 이는 실험 결과와 질적인 불일치를 초래한다. 즉, 기존 모델들은 친수성 기판은 더 친수성으로, 소수성 기판은 더 소수성으로 변한다고 예측한다. 또한, 이러한 모델들은 비대칭적인 고체-액체 및 고체-기체 계면 장력 변화를 필요로 하는 "오토포빙(autophobing, 자가 소수성 현상, 액적의 자발적 수축)"과 같은 현상을 설명하지 못한다. 저자들은 고체 표면 흡착의 부재가 이러한 불일치를 해결하기 위해 필요한 누락된 메커니즘임을 식별하였다.

방법론
저자들은 고체 표면 흡착과 계면 흡착을 모두 포함하여, 두 상 흐름 내 가용성 계면활성제에 대한 열역학적으로 일관된 위상장 모델을 개발하였다. 이 프레임워크는 제2법칙과 일치하는 변분 원리(variations principles)를 통해 도출되었다.

• 지배 방정식: 모델은 비압축성 나비에-스토크스(Navier-Stokes) 방정식과 위상장 변수($\phi$) 및 계면활성제 농도($\psi$)에 대한 칸-힐리야드(Cahn-Hilliard) 방정식을 결합한다. 벌크 자유 에너지는 혼합 에너지, 계면활성제 혼합 엔트로피, 계면 흡착 에너지, 그리고 계면활성제 용해도를 위한 벌크 페널티를 포함한다.
• 고체 흡착 메커니즘: 핵심적인 혁신은 벽면 자유 에너지를 정식화하는 것이다. 저자들은 국부적인 계면활성제 농도의 함수로서 고체-유체 표면 장력의 감소를 모델링하기 위해 랑뮈어-시슈코프스키(Langmuir-Szyszkowski) 상태 방정식을 도입하였다. 이를 통해 고체-액체 및 고체-기체 계면에서 비대칭적인 흡착 강도($\beta_{s1}$ 및 $\beta_{s2}$)를 설정할 수 있다.
• 경계 조건: 변분 유도는 벽면에서의 위상장 및 계면활성제 농도에 대한 운동학적 경계 조건을 도출한다. 결정적으로, 계면활성제 자유 에너지에 기울기 페널티 항($|\nabla \psi|^2$)을 포함함으로써 계면활성제 방정식을 공간에 대해 4차로 격상시켰다. 이러한 수학적 구조는 경계에서의 비자명한 계면활성제 농도의 법선 기울기를 규정하는 데 필수적이며, 이를 통해 고체 흡착 효과를 벌크 유체 분포로 전달한다.
• 수치 구현: 시스템은 2차 아담스-바슈포크(Adams-Bashforth) 스킴을 사용한 시차적 카테시안 격자(staggered Cartesian grid) 상에서 풀이되며, 나비에-스토크스 방정식을 위해 분할 단계법(fractional-step method)이 사용된다. 퇴화된 계면활성제 이동도(degenerate surfactant mobility)로 인해 발생하는 변수 계수 라플라시안을 처리하기 위해 방향 의존적인 전방-후방 분리(forward-backward splitting) 전략이 채택되었다.

주요 결과
모델은 깨끗한 액적 퍼짐에 대한 해석적 벤치마크를 통해 검증되었으며, 실험 데이터 및 기존 수치 연구와 비교되었다.

1. 젖음성 불일치의 해결: 고체 흡착을 포함했을 때, 모델은 모든 초기 접촉각에 대해 보편적인 친수성 증가(접촉각 감소)를 예측한다. 이는 계면활성제 첨가 시 친수성 및 소수성 표면 모두가 더 잘 젖게 된다는 실험적 관찰과 일치한다. 반면, 벽면 흡착이 없는 시뮬레이션은 기판의 기존 젖음 특성을 증폭시키는 오류 섞인 경향을 재현하였다.
2. 계면 곡률 및 마랑고니 응력: 고체 흡착을 포함하면 접촉선 근처에서 유체-유체 계면의 미세하지만 유의미한 굽힘(bending)이 나타난다. 이 곡률은 결합 메커니즘으로부터 발생한다: 고체 흡착은 벽 근처에서 계면활성제 농도 구배를 생성하여 벌크 흐름을 구동하는 마랑고니 응력을 발생시킨다. 비압축성 제약 조건은 이 흐름을 국부적인 압력장으로 변환하며, 이는 법선 응력 점프(normal stress jump)를 수정하여 계면이 곡률을 조정하도록 강제한다. 이 특징은 고체 흡착이 없는 모델에는 존재하지 않는다.
3. 오토포빙 예측: 흡착 비대칭성을 반전시킴으로써(고체-주변 환경 계면에 강한 흡착, 고체-액체 계면에 약한 흡착), 모델은 오토포빙을 성공적으로 예측한다. 이 영역에서 액적은 더 큰 평형 접촉각을 향해 자발적으로 수축한다. 시뮬레이션 결과, 접촉각이 약 $15^\circ$에서 $19^\circ$ 정도 증가하였으며, 이는 문헌에 보고된 실험 범위와 일치한다. 이 현상은 접촉선 앞의 "마른" 기판에 계면활성제가 우선적으로 축적되어, 고체-액체 계면 대비 고체-주변 환경 계면 장력을 급격히 감소시킴으로써 발생한다.
4. 스케일링 법칙: 저자들은 벽에서의 계면활성제 흡착층 두께에 대한 스케일링 법칙을 도출하였다. 이들은 두 가지 영역을 식별하였다: 하나는 층 두께가 계면활성제 확산 매개변수에 의존하는 영역(저기울기 영역)이고, 다른 하나는 층 두께가 벽에서의 계면활성제 강도에 의존하는 영역(고기울기 영역)이다.

의의 및 주장
본 논문은 고체 표면 흡착의 포함이 계면활성제가 포함된 접촉선 역학을 예측하는 데 필요한 "누락된 메커니즘"을 제공한다고 주장한다. 저자들은 이 메커니즘 없이는 수치 모델이 보편적인 젖음성 증가 경향을 질적으로 재현하거나 오토포빙과 같은 복잡한 현상을 설명할 수 없다고 단언한다.

본 연구는 다음을 확립한다:

• 고체 흡착은 단순히 부차적인 효과가 아니라, 접촉선에서의 힘의 근본적인 균형을 바꾸는 중심적인 동력이다.
• 제안된 모델의 열역학적 일관성, 특히 기울기 페널티를 포함한 자유 에너지 범함수로부터 경계 조건을 도출하는 것은 계면활성제 분포, 마랑고니 응력, 그리고 계면 곡률 사이의 결합을 포착하는 데 필수적이다.
• 본 모델은 벽면 흡착의 비대칭성만으로 향상된 퍼짐과 자발적 탈습(dewetting)을 모두 설명할 수 있는 통합된 프레임워크를 제공함으로써, 기존의 수치적 예측과 실험적 실재 사이의 간극을 성공적으로 메운다.

저자들은 본 연구가 기존 공식과 실험 사이의 질적인 불일치를 해결하였으며, 고체 흡착이 계면활성제가 포함된 젖음 현상을 정확하게 모델링하는 데 필요한 결정적인 요소임을 입증하였다고 결론짓는다.