On the deformation of a shear thinning viscoelastic drop in a steady electric field

이 논문은 오픈소스 솔버 Basilisk 를 이용한 수치 시뮬레이션을 통해 전도도 및 유전율 비율에 따라 정의된 다양한 영역에서 균일 전기장에 노출된 전단 박화 점탄성 (LPTT) 액적의 변형 및 파열 거동을 분석하고, 특히 탄성 증가가 변형을 억제하거나 비단조적 거동을 유발하는 등 뉴턴 유체와 구별되는 역학적 특성을 규명했습니다.

핵심

🧪 연구의 핵심: "끈적한 액방울 vs 전기장"

상상해 보세요. 전구 속의 전구처럼 두 개의 전극 사이에 끈적끈적한 액체 방울을 넣었습니다. 여기에 전기를 흘려보내면 어떻게 될까요?

1. 일반적인 물방울 (뉴턴 유체): 전기를 받으면 길쭉하게 늘어나다가, 전기 세기가 너무 강해지면 그냥 톡 하고 터져버립니다.
2. 이 연구의 주인공 (점탄성 액체): 이 액체는 치킨 스프링이나 쫄깃한 치약처럼 탄력이 있습니다. 당기면 늘어나지만, 다시 원래 모양으로 돌아오려는 성질도 있고, 너무 당기면 끊어지기도 합니다.

연구자들은 이 탄력 있는 액방울이 전기장 안에서 어떻게 변형되고, 언제 터지는지, 그리고 그 모양이 어떻게 바뀌는지 다양한 조건에서 실험했습니다.

---

🔍 주요 발견 3 가지 (비유로 설명)

연구는 전기의 세기와 액체의 성질에 따라 6 가지 다른 '무대' (영역) 로 나뉘어 관찰되었습니다. 그중 가장 흥미로운 세 가지 상황을 비유로 풀어보겠습니다.

1. "탄력 있는 스프링"의 저항 (PR+A 영역)

• 상황: 전기장이 약할 때는 방울이 길쭉한 타원형 (계란 모양) 으로 변합니다.
• 발견: 전기 세기를 조금 더 세게 하면, 보통 물방울은 바로 여러 조각으로 부서집니다. 하지만 탄력 있는 액체는 다릅니다.
  • 비유: 마치 튼튼한 고무줄을 당기는 것과 같습니다. 처음에는 잘 늘어나지만, 탄력 (Deborah 수, De) 이 강할수록 더 많이 당겨도 잘 끊어지지 않습니다.
  • 결과: 탄력이 강한 액체일수록 더 강한 전기를 받아야만 터집니다. 그리고 터질 때도 그냥 톡 터지는 게 아니라, 두 개의 방울이 붙어 있다가 나뉘거나, 세 개의 방울로 갈라지는 등 더 복잡한 모양을 보입니다. 탄력이 강할수록 터지는 시간이 늦어집니다.

2. "뾰족한 끝"이 생기는 마법 (PR+B 영역)

• 상황: 전기장의 방향과 액체의 성질에 따라 방울 끝이 뾰족하게 변하는 구간입니다.
• 발견: 여기서 재미있는 일이 일어납니다. 탄력 (De) 이 약할 때는 뾰족해지지 않다가, 탄력이 어느 정도 강해지면 갑자기 뾰족한 끝이 생깁니다. 그런데 탄력이 너무 강해지면 다시 뾰족해지지 않고 평평해지기도 합니다.
  • 비유: 악어 입을 상상해 보세요. 입을 벌리는 힘 (전기) 과 악어의 턱 근육 (탄력) 의 싸움입니다.
    • 근육이 약하면 입을 못 벌립니다.
    • 근육이 적당히 강하면 입을 크게 벌려 뾰족하게 만듭니다.
    • 하지만 근육이 너무 강하면 (너무 탄력이 있으면), 오히려 입을 다물고 싶어 하거나, 입 모양이 변하지 않고 유지됩니다.
  • 결과: 탄력이 너무 강해지면 오히려 변형이 줄어들기도 해서, 가장 잘 변형되는 탄력 정도가 따로 존재한다는 것을 발견했습니다.

3. "납작하게 눌리는" 상황 (OB- 영역)

• 상황: 전기장이 방울을 납작하게 (원반 모양) 누르는 경우입니다.
• 발견: 탄력이 약하면 전기 세기가 조금만 강해져도 납작한 방울이 바로 쪼개져 버립니다. 하지만 탄력이 강해지면, 더 강한 전기를 받아도 납작한 모양을 유지하다가, 어느 순간에야 터집니다.
  • 비유: 젤리를 손으로 누르는 상황입니다.
    • 일반 젤리 (탄력 없음) 는 살짝만 누르면 찢어집니다.
    • 쫄깃한 젤리 (탄력 있음) 는 꾹꾹 눌러도 찢어지지 않고 납작하게 변형되다가, 한계점을 넘으면 톡 하고 터집니다.
  • 결과: 탄력이 강할수록 터지기까지 버티는 힘이 훨씬 강해집니다.

---

🆚 중요한 비교: "무한한 스프링" vs "한계가 있는 스프링"

이 논문은 기존에 많이 쓰이던 올드로드-B (Oldroyd-B) 모델과 새로 적용한 선형 판 - 테이 - 탠너 (LPTT) 모델을 비교했습니다.

• 올드로드-B 모델 (구식): 이 모델은 고분자 사슬이 무한히 늘어날 수 있다고 가정합니다. (현실적이지 않음)
  • 비유: 끝이 없는 고무줄처럼, 당기면 당길수록 계속 늘어나고 결국은 무한한 힘으로 터집니다.
• LPTT 모델 (현실적): 이 모델은 고분자 사슬에 **최대 길이 (한계)**가 있다고 봅니다. (현실의 플라스틱이나 고분자 용액과 비슷함)
  • 비유: 최대 길이까지만 늘어나는 스프링입니다. 일정 이상 당기면 더 이상 늘어나지 않고, 오히려 단단해지거나 (스트레인 하딩) 힘을 견디다가 끊어집니다.

결론적으로:
기존 모델 (올드로드-B) 은 탄력이 강해지면 액방울이 더 쉽게 변형되고 터진다고 예측했지만, **이 연구의 현실적인 모델 (LPTT)**은 탄력이 강해지면 오히려 더 잘 견디고, 터지는 시점이 늦어지며, 모양도 더 안정적임을 보여주었습니다.

---

💡 이 연구가 왜 중요할까요?

이 연구는 단순히 물리학 실험이 아닙니다. 우리 일상의 기술에 큰 영향을 줍니다.

1. 잉크젯 프린터: 잉크 방울이 종이에 떨어지기 전, 전기장으로 어떻게 모양을 잡는지 이해하면 더 선명한 인쇄가 가능합니다.
2. 마이크로 유체 칩: 작은 액적들을 전기로 조종할 때, 액체가 너무 빨리 터지지 않게 하거나 원하는 모양으로 변형시키는 데 이 지식이 필요합니다.
3. 기름 정제: 기름과 물이 섞인 액체를 전기로 분리할 때, 점성이 있는 액체들이 어떻게 행동하는지 알면 공정을 효율화할 수 있습니다.

한 줄 요약:

"전기장 안에서 탄력 있는 액체 방울은 우리가 생각했던 것보다 훨씬 튼튼하고 복잡한 모양을 만들며, 탄력이 강할수록 더 잘 버틴다는 것을 발견했습니다. 이는 더 정교한 미세 액적 제어 기술의 길을 열어줍니다."

기술 요약

논문 요약: 정상 전기장 하의 전단 박리 점탄성 액적 변형 연구

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

• 배경: 전기장 하에서의 액적 변형 및 분열 역학은 미세유체공학, 잉크젯 프린팅, 전기수력학적 유체 조작 등 다양한 공학적 응용 분야에서 중요합니다. 기존 연구는 주로 뉴턴 유체나 Oldroyd-B 모델 (무한한 신장성 가정) 을 기반으로 수행되었습니다.
• 문제점: 많은 실제 유체 (고분자 용액, 생체 유체 등) 는 전단 박리 (shear thinning) 특성을 가지며, 고분자 사슬의 신장에는 물리적 한계 (finite extensibility) 가 존재합니다. Oldroyd-B 모델은 이러한 전단 박리 특성을 고려하지 못하며, 높은 변형률에서 비물리적인 무한한 응력 발산을 예측합니다.
• 연구 목표: 선형 Phan-Thien-Tanner (LPTT) 모델을 사용하여 전단 박리 특성과 유한한 고분자 신장성을 고려한 점탄성 액적의 변형 및 분열 거동을 정상 전기장 하에서 수치적으로 분석하고, 기존 Oldroyd-B 모델과의 차이점을 규명하는 것.

2. 방법론 (Methodology)

• 수치 해석 도구: 오픈소스 솔버인 Basilisk를 사용하여 축대칭 (axisymmetric) 2 차원 시뮬레이션을 수행했습니다.
• 물리 모델:
  • 유체 모델: 액적은 선형 Phan-Thien-Tanner (LPTT) 구성 방정식으로 모델링되었으며, 주변 유체는 뉴턴 유체로 가정했습니다.
  • 전기장: Maxwell 응력 텐서를 기반으로 한 전기체적력 (electric body force) 을 고려하며, 전하 대류 (charge convection) 효과를 포함했습니다.
  • 무차원 수: 전기 캐필러리 수 ($Ca_E$), 데보라 수 ($De$, 점탄성 효과의 척도), 전도도 비 ($\sigma_r$), 유전율 비 ($\epsilon_r$) 등을 주요 변수로 사용했습니다.
• 검증: Lid-driven cavity flow 문제를 통해 LPTT 유체의 수치 해석 정확도를 검증하고, 격자 무관성 (grid independence) 및 계산 영역 크기를 확인했습니다.
• 연구 범위: $(\sigma_r, \epsilon_r)$ 위상 공간의 6 가지 특성 영역 ($PR^+_A, PR^+_B, PR^-_A, PR^-_B, OB^+, OB^-$) 에서 대표 조합을 선정하여 다양한 $Ca_E$와 $De$ 조건에서 시뮬레이션을 수행했습니다.

3. 주요 결과 (Key Results)

가. 전단 박리 및 유한 신장성의 영향 (LPTT vs. Oldroyd-B)

• LPTT 모델의 특징: LPTT 모델은 높은 변형률에서 응력이 포화되는 현상 (strain hardening 후 응력 제한) 을 보여 물리적으로 더 현실적인 거동을 예측합니다. 반면 Oldroyd-B 모델은 무한한 신장성으로 인해 과도한 변형과 조기 분열을 예측하는 경향이 있습니다.

나. 위상 영역별 거동 분석

1. $PR^-_A, PR^-_B, OB^+$ 영역:

   • 1 차 및 2 차 변형 계수의 부호가 반대인 영역으로, 액적은 안정된 구형 (spheroidal) 또는 타원형 형태를 유지하며 분열되지 않습니다.
   • 이 영역에서는 점탄성 효과가 미미하여 뉴턴 유체와 거의 유사한 거동을 보입니다.

2. $PR^+_A$ 영역 (전도도 우세):

   • 변형: $Ca_E$가 임계값 이하일 때 안정된 장방형 (prolate) 형태를 유지하지만, 임계값을 초과하면 다엽 (multi-lobed) 형태를 형성하거나 분열합니다.
   • 점탄성 효과: $De$가 증가함에 따라 액적의 변형이 감소하고 임계 $Ca_E$가 증가하여 분열에 대한 저항이 커집니다.
   • 분열 양상: $De$가 높을수록 분열 시나리오가 변화하며 (예: 3 엽에서 2 엽으로 전환), 분열 시간이 지연되는 경향을 보입니다.

3. $PR^+_B$ 영역 (유전율 우세):

   • 형태 전이: 임계 $Ca_E$ 이하에서는 장방형 형태를 유지하지만, 초과 시 뾰족한 끝 (pointed ends) 을 가진 형태로 변형됩니다.
   • 비단조적 거동 (Non-monotonic behavior): 변형 정도와 임계 $Ca_E$가 $De$에 대해 비단조적으로 변화합니다. 즉, 낮은 $De$에서는 변형이 증가하다가 높은 $De$에서는 감소하는 경향을 보입니다. 이는 중간 $De$ 영역에서 점탄성이 변형을 억제하지만, 매우 높은 $De$에서는 그 저항이 감소하여 뉴턴 유체보다 더 큰 변형을 허용함을 의미합니다.
   • Oldroyd-B 와의 차이: LPTT 액적은 Oldroyd-B 액적보다 더 큰 변형을 보이며, 뾰족한 끝 형태가 더 쉽게 형성됩니다.

4. $OB^-$ 영역 (구형/평평한 형태):

   • 변형: 액적은 평평한 형태 (oblate) 로 변형하다가 임계 $Ca_E$를 초과하면 분열합니다.
   • 비단조적 거동: 변형 크기는 $De$가 증가함에 따라 처음에는 증가하다가 (strain hardening 효과) 이후 감소합니다 (stress saturation 효과).
   • 안정성: 높은 $De$에서 LPTT 액적은 Oldroyd-B 액적에 비해 분열에 대한 저항이 훨씬 크며, 더 안정적입니다. Oldroyd-B 는 높은 $De$에서 오히려 변형이 커지고 분열이 빨라지는 반면, LPTT 는 변형이 억제됩니다.

4. 주요 기여 및 결론 (Contributions & Significance)

• 모델의 현실성 제고: 전단 박리 특성과 유한한 고분자 신장성을 고려한 LPTT 모델을 전기수력학적 (EHD) 액적 변형 연구에 적용함으로써, 고분자 용액 등 실제 유체의 거동을 더 정확하게 예측할 수 있음을 입증했습니다.
• 점탄성 효과의 복잡성 규명: 점탄성 ($De$) 이 액적의 변형과 분열에 미치는 영향이 단순하지 않으며, 전기장 조건 ($\sigma_r, \epsilon_r$) 과 전기 캐필러리 수 ($Ca_E$) 에 따라 변형 저항이 증가하거나 감소하는 비단조적 거동을 발견했습니다.
• Oldroyd-B 모델의 한계 지적: 기존에 널리 사용되던 Oldroyd-B 모델이 높은 전기장 하에서 액적의 안정성을 과소평가하거나 (분열을 너무 일찍 예측), 변형 거동을 정확히 재현하지 못함을 보였습니다. 특히 $OB^-$ 영역에서 LPTT 모델은 높은 $De$에서 액적 안정성을 유지하는 반면, Oldroyd-B 는 불안정해지는 등 근본적인 차이를 보였습니다.
• 공학적 시사점: 전기장을 이용한 미세 유체 조작, 잉크젯 프린팅, 유체 분무 등 다양한 응용 분야에서 점탄성 유체의 거동을 제어하기 위해 전단 박리 특성과 유한 신장성을 고려한 설계가 필수적임을 시사합니다.

결론적으로, 본 연구는 전기장 하에서 전단 박리 점탄성 액적의 거동이 Oldroyd-B 모델로 예측되는 바와 질적으로 다를 수 있음을 보여주었으며, 특히 높은 데보라 수 ($De$) 영역에서 LPTT 모델이 예측하는 유한 신장성 효과가 액적의 안정성과 분열 역학에 결정적인 영향을 미친다는 점을 강조했습니다.