전기장 하 구형 점성탄성 드롭의 대변형 및 파괴 메커니즘

초록

본 연구는 Oldroyd‑B 모델을 이용해 균일 직류 전기장에 노출된 점성탄성 액적의 변형과 파괴를 수치적으로 조사한다. 전도도 비(σ_r)와 유전율 비(ε_r)의 여섯 가지 대표 조합을 선택해 전기 모자볼 캡릴러리 수(Ca_E)와 Deborah 수(De)의 영향을 체계적으로 분석하였다. 결과는 전도도·유전율 비에 따라 변형 형태가 크게 달라지며, 탄성은 주로 변형을 억제하고 임계 전기 캡릴러리 수를 상승시키는 역할을 함을 보여준다.

상세 분석

본 논문은 전기장 하에서의 점성탄성 액적 변형을 Oldroyd‑B 모델로 기술하고, 이를 Basilisk 오픈소스 솔버로 3차원 전산유체역학 시뮬레이션을 수행함으로써 기존의 레키‑다이얼렉트릭(LDM) 이론이 다루지 못한 비선형·대변형 현상을 정량적으로 파악한다. 연구는 (σ_r, ε_r) 파라미터 공간을 여섯 개 영역(PR_A^+, PR_B^+, PR_A^−, PR_B^−, OB^+, OB^−)으로 구분하고, 각 영역에서 대표적인 비율쌍을 선택하였다.

1. 구역별 변형 특성

   • PR_A^−, PR_B^−, OB^+: 1차·2차 변형 계수가 동일 부호를 가져 뉴턴 유체와 거의 구분되지 않는다. 즉, 탄성 효과가 미미하고 전기 응력에 의한 변형이 전형적인 구형→타원형 전이를 보인다.
   • PR_A^+: 전기 캡릴러리 수 Ca_E가 임계값을 초과하면 다엽(lobed) 형태가 발생한다. De가 증가할수록 탄성 저항이 강화돼 변형이 억제되고, 임계 Ca_E가 상승한다. 이는 탄성 응력이 전기 응력에 비해 상대적으로 크게 작용함을 의미한다.
   • PR_B^+: 임계 Ca_E 이상에서 원뿔형(콘) 끝을 가진 형태가 나타난다. De가 낮을 때는 변형이 감소하지만, 일정 범위에서 비단조적(non‑monotonic) 변화를 보이며, 임계 Ca_E 자체도 De에 따라 상승·하강을 반복한다. 고 Ca_E·고 De 영역에서는 변형이 일시적으로 최대·최소값을 갖는 과도 현상이 나타나며, 구형과 뾰족한 형태 사이를 반복적으로 진동한다.
   • OB^−: 전기장에 의해 평평한(oblates) 형태로 압축되며, 일정 Ca_E를 초과하면 파열(breakup)된다. De가 증가하면 변형 크기가 커지고 임계 Ca_E는 감소한다. 저 Ca_E·고 De 조건에서는 드롭 표면에 움푹 들어간(dimpling) 현상이 발생하고, 중심 위치가 주기적으로 진동한다.

2. 탄성·전기 상호작용 메커니즘

   • Oldroyd‑B 모델의 고정 점성(η_s)과 탄성(η_p) 구성은 전기 응력에 대한 응답 시간을 Deborah 수(De=λ/τ_c)로 통합한다. De가 클수록 응력 이완이 느려져 전기장에 의해 가해지는 변형을 즉시 따라가지 못하고, 내부 탄성 응력이 축적된다. 이 축적된 탄성 응력은 전기 응력과 반대 방향으로 작용해 변형을 억제하고, 특히 PR_A^+, PR_B^+ 구역에서 임계 Ca_E를 상승시킨다.
   • 전도도 비와 유전율 비가 변함에 따라 전기장 내부와 외부에서 발생하는 정상·전단 전기 응력이 달라진다. 예를 들어 PR_B^+에서는 높은 전도도 비가 전기 전류를 강화해 전기 전단 응력이 크게 작용, 이에 따라 원뿔형 끝이 형성된다. 반면 OB^−에서는 높은 유전율 비가 전기 압축 응력을 강화해 평평한 형태를 유도한다.

3. 수치적 검증 및 그리드 독립성

   • 저자들은 Basilisk의 적응형 격자(AMR) 기능을 활용해 인터페이스 근처를 고해상도로 재구성하고, 격자 수렴 테스트를 수행해 결과의 수치적 신뢰성을 확보하였다. 특히 고 De·고 Ca_E 상황에서 발생하는 뾰족한 끝(conical tip)이나 다엽 구조는 매우 미세한 인터페이스 곡률을 필요로 하므로, 2^−10 수준의 최소 셀 크기를 사용해 정확도를 검증하였다.

4. 실용적 함의

   • 마이크로플루이딕 디스펜서, 잉크젯 프린팅, 전기적 스프레이 등에서 액적의 형태 제어는 품질과 효율에 직접적인 영향을 미친다. 본 연구는 De를 조절(예: 폴리머 용액의 점탄성 조절)함으로써 원하는 변형 모드(예: 다엽형 vs 원뿔형)와 파괴 임계값을 설계할 수 있음을 제시한다. 특히 PR_A^+ 구역에서는 De를 높여 다엽형 파열을 억제하고, PR_B^+ 구역에서는 De를 최적화해 원뿔형 끝을 이용한 정밀 분사(nozzle‑free) 기술에 활용할 수 있다.