Droplet Spreading on Nanowall-Textured Hydrophobic Surfaces: A Nonlinear Reduced-Order Model

기술 요약: 나노월(Nanowall) 질감을 가진 소수성 표면에서의 액적 퍼짐

문제 정의
비젖음성(non-wetting) 표면에 대한 액적 충돌은 병진 운동 에너지가 계면 에너지, 내부 운동 및 소산(dissipation)으로 전환되는 복잡한 과도 현상을 포함한다. 고전적인 스케일링 법칙(예: $\beta_{max} \sim We^{1/4}$)은 최대 변형 경향을 요약하지만, 표면 질감이 근본적인 힘과 에너지 균형을 어떻게 수정하는지는 유일하게 식별하지 못한다. 구체적으로, 표면 질감은 서로 경쟁하는 효과를 도입한다: 평형 젖음 기하학의 변화(Wenzel 대 Cassie–Baxter 상태), 접촉선을 고정시키는 에너지 장벽(pinning), 그리고 소산 및 회복 가능한 에너지를 변화시킬 수 있는 질감 내부로의 액체 침투 가능성 등이 그것이다.

본 연구에서 다루는 구체적인 설정은 Yilbas 등[11]의 실험 데이터를 기반으로 한, 기능화된 매끄러운 실리콘 표면과 나노월 질감 표면 위의 물 액적 충돌이다. 과제는 계산 비용이 많이 드는 완전 해상 다상 시뮬레이션(네 자릿수의 길이 스케일을 모두 분해해야 함)에 의존하지 않고, 관찰된 최대 퍼짐의 방향성과 조건 의존성을 재현할 수 있는 가장 단순한 방사형 동역학 시스템을 결정하는 것이다.

방법론
저자들은 응집된 방사형 대리 모델(lumped radial surrogate)로서 작동하는 비선형 1자유도(1-DOF) 차수 축소 모델을 공식화한다. 이 모델은 해상도가 높은 부피 유동(volume-of-fluid) 또는 위상장(phase-field) 시뮬레이션이 아니다.

1. 지배 동역학: 시스템은 참조 값 $q_0 = D_0/2$를 갖는 유효 반경 확장 $q(t)$를 추적한다. 운동 방정식은 다음과 같다:
   $$m\ddot{q} = -k_{cap}(q - q_0)\left[1 + \kappa_s \left(\frac{q - q_0}{q_0}\right)^2\right] - c_{visc}\dot{q} - F_{pin}\tanh\left(\frac{\dot{q}}{v_s}\right)$$

   • 질량 ($m$): 총 액적 질량과 동일하게 설정된다.
   • 모세관 강성 ($k_{cap}$): 관성-모세관 시간 척도에 의해 결정되며, 피팅(fitting)되지 않는다.
   • 비선형 강성: 3차 항($\kappa_s$)은 큰 변형 시 비현실적인 유연성을 방지하기 위해 진폭 의존적 강성 증가를 제공한다.
   • 점성 감쇠: 레이놀즈 수에 의존하는 계수로 모델링된다.
   • 운동 저항: 고속에서 $F_{pin}$의 크기에 도달하며 운동에 저항하는 정규화된 항이다.

2. 하이브리드 피닝 규칙: 정규화된 운동 항이 제로 속도에서 사라지므로(정적 피닝을 방지하기 위해), 명시적인 **이동-정지 이벤트(moving-to-stuck event)**가 부과된다. 궤적은 속도가 캡처 밴드($v_{stick}$) 미만으로 떨어지고 모세관 구동력이 정적 용량 임계값($\eta_s F_{pin}$)을 극복하기에 불충분할 때 종료된다.

3. 피닝 폐쇄(Pinning Closure): $F_{pin}$의 크기는 고정된 참조 조우 횟수(encounter count)로 스케일링된 기공 규모 힘 모델을 사용하여 추정된다. 특히, 수력 직경이 질감 간격($D_H = 2a$)과 같다는 가정하에, 무차원 피닝 크기는 간격 및 액적 크기 모두로부터 독립적이다.

4. 교정 및 검증:

   • $\kappa_s$는 단일 그래픽 데이터 포인트(매끄러운 표면, $V=20\,\mu$L, $h=20\,$mm)에 맞춰 피팅된다.
   • 질감 표면 응답은 추가적인 피팅 없이 고정된 파라미터를 사용하여 예측된다.
   • 민감도 분석: 폐쇄 파라미터($\kappa_s, F_{pin}, \zeta_0, v_s$)와 수치적 설정에 대해 국부 결정론적 민감도 분석을 수행하여 물리적 제어 요인과 수치적 아티팩트를 구분한다.

주요 기여

• 모델 계층 구조: 본 논문은 비선형 강성과 피닝 저항의 구체적인 역할을 격리하기 위해 모델 계층(선형 비피닝, 비선선형 비피닝, 비선형 피닝)을 구축한다.
• 명시적 하이브리드 동역학: 정적 피닝을 처리하기 위해 수학적으로 명시된 이동-정지 규칙을 도입하며, 단순한 정규화만으로는 피닝된 상태를 유지할 수 없음을 인정한다.
• 식별 가능성 및 투명성: 본 연구는 파라미터를 문헌 값, 재구성된 입력값, 피팅된 계수, 또는 가정된 폐쇄 값으로 명확히 분류한다. 또한, 한 조건에서의 일치가 모든 폐쇄를 검증할 수는 없음을 보여준다.
• 범위 정의: 모델은 최대 반경 확장 계수($\beta_{max}$)를 예측하는 데 엄격히 제한된다. 이는 수직 반등, 접촉 시간, 젖음 상태 전이 또는 물리적 이탈(lift-off) 예측을 명시적으로 제외한다.

결과

• 기본 케이스 ($V=20\,\mu$L, $h=20\,$mm):
  • 선형 비피닝 모델은 확장을 과대 예측한다 ($\beta_{max} = 3.955$).
  • 비선형 비피닝 모델(교정됨)은 매끄러운 표면 목표치와 일치한다 ($\beta_{max} = 2.253$).
  • 비선형 질감 모델은 감소된 확장을 예측한다 ($\beta_{max} = 1.810$), 이는 Yilbas 등이 관찰한 나노월 표면에서의 감소된 확산의 질적 경향과 일치한다.
• 조건 전이: 파라미터가 고정된 상태에서, 모델은 재구성된 웨버 수($We$)에 따라 두 표면 모두에서 최대 확장이 단조 증가함을 예측한다. 질감 표면은 조사된 부피 및 방출 높이 범위 전반에서 매끄러운 표면보다 일관되게 낮은 $\beta_{max}$를 나타낸다.
• 민감도:
  • **비선형 강성($\kappa_s$)**과 **피닝 크기($F_{pin}$)**가 최대 확장의 주요 제어 요인으로 식별되었다.
  • 피닝 크기는 질감 표면에 특화된 지배적 제어 요인이다.
  • 수치적 파라미터(솔버 허용 오차, 이벤트 국소화, 운동 평활화 속도)는 보고된 정밀도 내에서 예측된 최대 확장에 미치는 영향이 미미하다.
  • 종료 시간(Endpoint Time): 궤적이 종료되는 시간은 스위칭 임계값 근처의 정적 용량 계수($\eta_s$)와 피닝 크기에 매우 민감하며, 이는 종료 시간이 물리적 예측이라기보다는 상태 분류 진단임을 나타낸다.

의의 및 주장
본 논문은 가용한 최대 확산 정보로부터 무엇을 추론할 수 있고 무엇을 할 수 없는지를 정의하는 재현 가능하고 계산 비용이 적은 방사형 대리 모델을 제공한다는 점을 겸허히 주장한다.

• 나노월 질감이 매끄러운 표면에 비해 최대 확산을 감소시킨다는 질적 경향을 성공적으로 재현한다.
• 확산 감소가 외부 피닝 작업과 하이브리드 정적 용량 규칙의 결합된 효과에 의해 구동됨을 명확히 한다.
• 저자들은 모델이 수직 반등, 법선 복원 계수 또는 물리적 접촉 시간을 예측하지 않음을 강조한다. 질감 케이스의 "제로" 반경 회복 프록시는 반경 방향 비회복에 대한 분류 관례이지, 실험적 반등 높이가 0임을 예측하는 것이 아니다.
• 본 연구는 복잡한 젖음 상태 전이 및 에너지 회복 메커니즘을 포착하는 데 있어 단일 좌표 접근법의 한계를 강조하며, 향후 2좌표 또는 해상도가 높은 다상 모델을 위한 명시적인 베이스라인 역할을 한다.