Viscoelastic Droplet Impact on Surfaces with Sharp Wettability Contrast: Coupled Influence of Relaxation Time and Surface Tension

본 논문은 수치 시뮬레이션을 통해 점탄성 액적의 충돌 시 이완 시간과 표면 장력이 혼합 젖음성 표면을 가진 고체 표면에서의 확산, 수축, 반발 및 최종 형태에 미치는 복합적 영향을 규명하고, 이완 시간 증가가 탄성 에너지 저장을 통해 확산을 촉진하는 반면 표면 장력 증가는 확산을 억제하고 비대칭적인 형태 변화를 유도함을 보였습니다.

핵심

🍯 1. 연구의 주인공: "쫀득한 액적"과 "반쪽짜리 바닥"

• 액적 (방울): 일반적인 물방울은 그냥 퍼지거나 튕겨 나갑니다. 하지만 이 연구의 주인공은 고분자 (폴리머) 가 섞인 액체입니다.
  • 비유: 물방울이 아니라 치약이나 달걀 흰자처럼 늘어났다가 원래 모양으로 돌아오려는 성질 (탄성) 이 있는 액체라고 상상해 보세요.
• 표면: 바닥은 반쪽은 **물기를 아주 좋아하는 (친수성)**곳이고, 다른 반쪽은 **물을 아주 싫어하는 (소수성)**곳입니다.
  • 비유: 바닥의 왼쪽은 스펀지처럼 물을 빨아들이고, 오른쪽은 테플론 프라이팬처럼 물을 튕겨내는 상태입니다.

🎯 2. 연구의 핵심 질문

"이 쫀득한 액적이 스펀지 - 테플론 바닥에 떨어졌을 때, 액체의 **'쫀득함 (탄성)'**과 **물방울끼리 뭉치는 힘 (표면 장력)**이 어떻게 상호작용할까?"

🔍 3. 주요 발견 (세 가지 비유로 정리)

① "쫀득함"을 늘리면? (완화 시간의 영향)

액체가 얼마나 오래 '쫀득하게' 버틸 수 있는지를 완화 시간이라고 합니다. 이 시간을 길게 만들면 (액체가 더 끈적해질수록):

• 현상: 액체가 바닥에 닿았을 때 더 넓게 퍼집니다.
• 비유: 쫀득한 카라멜을 바닥에 떨어뜨리면, 일반 물방울보다 더 많이 퍼져나가면서 바닥을 감싸듯 넓게 깔립니다. 그 이유는 액체 내부에 탄성 에너지가 더 많이 저장되어, 퍼지는 힘을 더 오래 유지하기 때문입니다.
• 결과: 최대 퍼짐 크기가 약 13% 더 커지고, 액체 높이는 더 낮아집니다.

② "물기 뭉치는 힘"을 세게 하면? (표면 장력의 영향)

액체 분자들이 서로 더 강하게 붙어있도록 하는 힘 (표면 장력) 을 키우면:

• 현상: 액체가 바닥에 퍼지는 것을 막아줍니다.
• 비유: 방수 코팅을 한 것처럼 액체가 뭉쳐서 퍼지지 않으려 합니다. 퍼지는 범위는 약간 줄어들지만 (약 1%), 대신 높이가 더 두꺼워지고 튕겨 나가는 힘 (후퇴) 이 강해집니다.
• 결과: 넓게 퍼지는 대신, 더 둥글고 두툼하게 남습니다.

③ "반쪽짜리 바닥"의 마법: "주걱 모양"과 "신발 모양"

가장 재미있는 부분은 바닥이 반쪽씩 다를 때입니다.

• 현상: 액체가 **물기를 좋아하는 쪽 (스펀지)**으로 쏠리면서 모양이 기형적으로 변합니다.
• 비유:
  • 위에서 본 모습: 마치 **주걱 (Dustpan)**처럼 한쪽은 평평하고 다른 쪽은 올라간 모양이 됩니다.
  • 옆에서 본 모습: 마치 **신발 (Shoe)**처럼 앞코가 올라간 모양이 됩니다.
• 이유: 물을 좋아하는 쪽은 액체를 당겨서 넓게 퍼뜨리고, 물을 싫어하는 쪽은 액체를 밀어내서 튕겨냅니다. 이 서로 다른 힘이 액체를 한쪽으로 몰아넣어 독특한 모양을 만들어냅니다.

🚀 4. 왜 이 연구가 중요할까요?

이 연구는 단순히 물방울 실험이 아니라, 실제 산업에 큰 도움을 줍니다.

• 잉크젯 프린팅: 잉크가 종이나 플라스틱에 떨어졌을 때, 원하는 모양으로 정확히 퍼지게 하려면 액체의 '쫀득함'과 바닥의 '성질'을 어떻게 조절해야 하는지 알려줍니다.
• 스프레이 코팅: 페인트나 약품을 뿌릴 때, 특정 부분에만 골고루 묻히게 하거나, 원하지 않는 곳으로 퍼지지 않게 막을 수 있습니다.
• 마이크로 유체: 아주 작은 채널에서 액체를 움직일 때, 이 원리를 이용하면 액체를 원하는 방향으로만 흐르게 할 수 있습니다.

💡 요약

이 논문은 **"쫀득한 액체"**가 **"반쪽짜리 바닥"**에 떨어질 때, 액체의 탄성이 퍼짐을 도와주고, 표면 장력이 퍼짐을 막아주며, 바닥의 차이가 액체를 주걱 모양으로 변형시킨다는 것을 밝혀냈습니다.

이처럼 액체의 성질과 바닥의 성질을 잘 조절하면, 우리가 원하는 대로 액체를 조종할 수 있다는 귀중한 지식을 제공한 연구입니다.

기술 요약

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

• 배경: 잉크젯 프린팅, 스프레이 코팅, 미세유체공학 등 다양한 분야에서 액적의 확산, 수축, 반발 거동을 정밀하게 제어하는 것은 필수적입니다. 특히 고분자 용액, 현탁액, 생체 에어로졸과 같은 비뉴턴 유체 (점탄성 유체) 의 거동은 뉴턴 유체 모델로는 정확히 예측하기 어렵습니다.
• 문제점: 기존 연구들은 대부분 균일한 젖음성을 가정하거나, 중력과 모세관력의 상호작용을 체계적으로 고려하지 않았습니다. 또한, 점탄성 응력 완화, 접촉선 고정 (pinning), 그리고 급격한 젖음성 불연속성 (예: 초소수성과 초친수성의 경계) 간의 비선형 결합 효과에 대한 이해는 여전히 부족합니다.
• 연구 목표: 중력 - 모세관 평형 하에서 점탄성 유체의 완화 시간 (relaxation time) 과 표면 장력 (surface tension) 이 급격한 젖음성 대비를 가진 표면 (하이브리드 젖음성) 에 충돌하는 액적의 거동에 미치는 결합적 영향을 규명하는 것입니다.

2. 방법론 (Methodology)

• 수치 해석 도구: OpenFOAM 기반의 고충실도 3 차원 솔버 (RheoTool 패키지의 rheoInterFoam 사용) 를 개발하여 사용했습니다.
• 물리 모델:
  • 유체 모델: 점탄성 거동을 설명하기 위해 Oldroyd-B 구성 방정식을 적용했습니다.
  • 인터페이스 추적: VOF (Volume-of-Fluid) 방법을 사용하여 액 - 기 인터페이스를 추적했습니다.
  • 수치적 안정성: 높은 완화 시간에서의 수치적 불안정성을 해결하기 위해 로그 컨포메이션 (log-conformation) 형식을 도입하여 응력 텐서를 변환했습니다.
  • 접촉선 모델: 접촉선 속도에 의존하는 동적 접촉각 모델 (Kistler 의 Hoffman 함수 기반) 을 적용하여 젖음, 확산, 수축 거동을 정밀하게 모사했습니다.
• 시나리오 설정:
  • 직경 2 cm, 초기 낙하 속도 4 m/s 의 점탄성 액적이 중력 하에 표면에 충돌합니다.
  • 표면은 두 영역으로 나뉘어 급격한 젖음성 대비를 가집니다: Zone 1 (친수성, 정적 접촉각 0°) 과 Zone 2 (초소수성, 정적 접촉각 160°).
  • 변수: 완화 시간 ($\lambda$: 0.02 s ~ 0.12 s) 과 표면 장력 ($\sigma$: 0.05 N/m ~ 0.15 N/m) 을 변화시키며 Weber 수 (관성/모세관력) 와 Eötvös 수 (중력/모세관력) 의 변화를 체계적으로 분석했습니다.

3. 주요 기여 (Key Contributions)

• 첫 번째 통합 프레임워크: Oldroyd-B 점탄성, 동적 접촉각 이력 현상, 그리고 급격한 젖음성 불연속성을 중력 하에서 동시에 해결하는 완전 결합 3 차원 수치 프레임워크를 최초로 제시했습니다.
• 비선형 결합 효과 규명: 내부 점탄성 응력과 극단적인 이종 젖음성 패턴 간의 비선형 결합 메커니즘을 규명하여, 기존 뉴턴 유체 기반의 예측을 넘어서는 통찰을 제공했습니다.
• 새로운 평형 형태 발견: 급격한 젖음성 대비와 점탄성 효과의 상호작용으로 인해 발생하는 독특한 '호일 (dustpan) 모양' (상부 뷰) 과 '신발 (shoe) 모양' (측면 뷰) 의 평형 형태를 발견하고 그 물리적 기작을 설명했습니다.

4. 주요 결과 (Results)

A. 완화 시간 ($\lambda$) 의 영향

• 확산 증가: 완화 시간을 0.02 s 에서 0.12 s 로 증가시키면 탄성 에너지 저장량이 증가하여 최대 확산 직경이 약 12.9% 증가했습니다 (약 24.97 mm → 28.09~28.17 mm).
• 수축 지연: 탄성 기억 효과로 인해 점성 소산이 지연되어 확산 단계가 길어지고, 반발 (recoil) 이 지연되었습니다.
• 비대칭성 증폭: 하이브리드 표면에서 점탄성 유체는 친수성 영역으로의 유체 이동을 더욱 촉진하여 비대칭 확산을 강화했습니다.
• 최종 형태: 완화 시간이 길어질수록 액적은 더 두껍고 비대칭적인 형태를 유지하며, 친수성 쪽으로 더 많이 퍼지는 경향을 보였습니다.

B. 표면 장력 ($\sigma$) 의 영향

• 확산 억제: 표면 장력을 0.05 N/m 에서 0.15 N/m 로 증가시키면 모세관 복원력이 강해져 최대 확산 직경이 약 1.1% 감소했습니다 (27.21 mm → 26.90 mm).
• 수직 회복 촉진: 높은 표면 장력은 액적의 최소 높이를 약 3.3% 증가시켰으며 (2.12 mm → 2.20 mm), 수직 방향의 빠른 회복을 유도했습니다.
• 응력 재분배: 높은 표면 장력은 액적의 확산을 제한하고 곡률 구동 수축을 촉진하여, 점탄성 응력을 더 좁은 영역에 집중시키는 효과를 낳았습니다.

C. 하이브리드 표면에서의 거동

• 방향성 이동: 급격한 젖음성 대비 (0° vs 160°) 는 친수성 영역으로의 방향성 유체 이동을 유발했습니다.
• 고유한 형태: 친수성 영역은 액체를 끌어당겨 확산을 촉진하고, 소수성 영역은 반발을 유도하여 액적이 호일 모양 (top view) 과 신발 모양 (side view) 의 평형 상태에 도달하게 되었습니다. 이는 친수성 쪽의 접촉선이 고정되고 소수성 쪽에서 수축이 일어나면서 발생하는 비대칭적인 모세관 압력 구배 때문입니다.

5. 의의 및 결론 (Significance)

• 물리적 통찰: 탄성 응력, 모세관 작용, 그리고 표면 젖음성 구배 간의 복잡한 상호작용을 정량적으로 규명했습니다. 특히 완화 시간이 길수록 내부 유동 진동이 지속되고 비대칭 응력 분포가 유지됨을 확인했습니다.
• 공학적 적용: 잉크젯 프린팅, 스프레이 코팅, 미세유체 장치 등에서 점탄성 잉크나 유체의 확산, 반발, 최종 형태를 정밀하게 제어하기 위한 설계 지침을 제공합니다.
• 표면 공학: 비균일 젖음성 표면을 설계할 때 점탄성 유체의 특성을 고려해야만 원하는 패턴 형성 및 유체 제어가 가능함을 시사합니다.

이 연구는 점탄성 액적의 충돌 역학을 이해하는 데 있어 중력, 모세관력, 탄성력의 결합적 역할을 체계적으로 분석한 선구적인 작업으로 평가됩니다.