Numerical analysis of capillarity-driven thinning rheometry for polydisperse polymer solutions

기술 요약: 다분산 고분자 용액의 모세관 구동 가늘어짐 유변학에 대한 수치 해석

문제 정의
모세관 파단 신장 유변학(Capillary Breakup Extensional Rheometry, CaBER)은 고분자 용액의 신장 점도를 특성화하는 표준 기술이다. Oldroyd-B 모델로 설명되는 이상적인 단분산 시스템에서, 탄성 모세관(elastocapillary, EC) 영역에서의 필라멘트 가늘어짐 속도는 유체의 이완 시간($\tau$)에 반비례하며, 이를 통해 측정된 감쇄율($\dot{\epsilon}_{EC}$)로부터 특성 이완 시간 $\tau_{EC} = 2/(3\dot{\epsilon}_{EC})$를 정의할 수 있다. 그러나 실제 다분산 고분자 용액의 경우, $\tau_{EC}$의 해석은 모호하다. 기존 연구들은 $\tau_{EC}$와 전단 또는 선형 점탄성에서 측정된 이완 시간 사이의 불일치를 지적해 왔다. 특히, 최근의 이분산 폴리스티렌 혼합물(Calabrese et al., 2025)에 대한 실험은 고분자량(high-MW) 용액에 저분자량(low-MW) 종을 높은 농도로 첨가하더라도 측정된 $\tau_{EC}$에 미치는 영향이 미미함을 보여주었으며, 이는 가늘어짐 역학이 전체 앙상블 평균이 아닌 특정 분자량 분포의 하위 집합에 의해 지배됨을 시사한다. 분자량 분포, 흐름 이력, 그리고 그 결과인 $\tau_{EC}$를 연결하는 포괄적인 프레임워크가 결여되어 있었다.

방법론
저자들은 다분산성을 명시적으로 포함하여 원통형 고분자 필라멘트의 가늘어짐을 시뮬레이션하는 수치 모델을 개발하였다. 이 접근 방식은 다음을 결합한다:

1. 힘의 평형: 관성, 점성, 모세관 및 고분자 응력을 고려하여 필라멘트 반지름 $a(t)$에 대한 1차원 힘의 평형 방정식을 적용한다.
2. 구성 모델: 멀티 모드 FENE-PM (Finitely Extensible Nonlinear Elastic - Peterlin Mean) 모델을 사용한다. 고분자 용액은 서로 다른 "종"(서로 다른 분자량 분포를 나타내는 PS7 및 PS16 폴리스티렌 샘플)의 혼합물로 취급된다.
   • 선형 영역: 모델은 선형 점탄성을 기술하기 위해 사전 평균화된 유체역학적 상호작용을 포함하는 Zimm 모델을 활용하여 각 종에 대한 이완 시간 스펙트럼 $\tau_{is}$를 생성한다.
   • 비선형 영역: 모델은 FENE 스프링과 형태 의존적 항력(conformation-dependent drag) 항을 통해 유한 신장성을 통합한다. 이 항은 사슬이 코일에서 타원체로 늘어남에 따라 발생하는 유체역학적 마찰의 변화를 고려하여, 가장 긴 모드의 이완 시간($\hat{\tau}_{i1}$)을 수정한다.
3. 시뮬레이션: 필라멘트 반지름과 각 종/모드의 형태 텐서를 결정하는 상미분 방정식 시스템을 동시에 해결한다. 모델은 dioctyl phthalate (DOP) 내의 PS7 및 PS16 이분산 혼합물을 다룬 Calabrese et al. (2025)의 실험 데이터와 검증되었다.

주요 기여 및 결과

• 선택적 신장 메커니즘: 시뮬레이션은 EC 영역 동안 신장 속도 $\dot{\epsilon}$가 시간에 따라 감소함을 보여준다. 결과적으로, 유효 와이센베르크 수(Weissenberg number) $Wi = \dot{\epsilon}\tau > 0.5$ (코일-신장 전이 임계값)를 만족하는 고분자 사슬만이 계속 신장되어 탄성 응력에 유의미하게 기여한다. 이완 시간이 짧은(저분자량) 사슬들은 이 임계값 아래로 떨어져 이완되며, 더 이상 필라멘트를 지지하지 않는다.
• 고분자량 종의 지배력: 이완 시간은 분자량에 비례하므로($\tau \sim M^{1.5}$), 고분자량 종은 저분자량 종보다 $Wi > 0.5$를 더 오래 유지한다. 따라서 측정된$\tau_{EC}$는 고분자량 꼬리 부분(high-MW tail)에 크게 편향된 "응력 기여 하위 앙상블"의 사슬에 의해 결정된다. 이는 고분자량 용액에 저분자량 고분자를 추가하는 것이, 고분자량이 총 응력을 감당하기에 충분한 농도에 도달하면 $\tau_{EC}$를 변화시키지 못하는 이유를 설명한다.
• 형태 의존적 항력: 높은 농도에서의 실험 데이터를 일치시키기 위해서는 형태 의존적 항력의 포함이 필수적이었다. 이 항이 없으면 모델은 실험보다 너무 낮은 $\tau_{EC}$ 값을 예측하게 되며, 사슬이 늘어나고 상호작용함에 따라 관찰되는 이완 시간의 증가를 포착하지 못한다.
• 실험 매개변수 의존성: 본 연구는 $\tau_{EC}$가 유체의 고유한 물성이 아니라 실험 조건에 의존함을 확인한다. 구체적으로, 초기 필라멘트 직경($a_0$)과 사전 신장 이력(판 분리 속도)은 EC 영역이 확립되기 전 짧은 모드들이 이완될 수 있는 시간을 결정한다. 더 큰 초기 직경과 더 높은 사전 신장 속도는 더 긴 측정 $\tau_{EC}$ 값으로 이어진다.
• 평균의 재평가: 논문은 "신장성 가중 분자량"($M_L$)과 "응력 가중" 분자량을 대조한다. $M_L$이 응력에 기여할 수 있는 고분자량 사슬의 잠재력을 포착한다면, 시뮬레이션의 최소 신장 속도에서 도출된 응력 가중 평균은 EC 영역 동안 활성화된 사슬의 실제 하위 집합을 더 정확하게 반영한다.

의의 및 주장
본 논문은 다분산 유체의 CaBER 실험에서 추출된 특성 가늘어짐 시간 $\tau_{EC}$를 유체의 단일하고 고유한 이완 시간으로 해석해서는 안 된다고 주장한다. 대신, 이는 탄성 모세관 균형이 달성되는 순간 흐름에 의해 능동적으로 신장되는 특정 사슬 하위 집합의 이완 시간을 반영하는 실험 특이적 양이다.

저자들은 다분산 용액의 경우 $\tau_{EC}$가 다음의 함수라고 결론짓는다:

1. 분자량 분포 (특히 고분자량 꼬리 부분).
2. 총 고분자 농도.
3. 초기 필라멘트 직경 및 사전 신장 이력과 같은 실험 매개변수.

이러한 발견은 복잡한 유체에서 전단/선형 점탄성 이완 시간과 신장 가늘어짐 시간 사이의 불일치를 해결한다. 이는 고분자 필라멘트의 "방사성(spinnability)"과 안정성이 고분자량 분획에 의해 불균형하게 영향을 받는다는 중요한 통찰을 제공하며, 산업용 섬유 방사부터 생물학적 시스템의 액적 파쇄에 이르는 현상을 이해하는 데 핵심적인 역할을 한다. 본 연구는 단일 모드 Oldroyd-B 해석의 한계를 넘어, 다분산 유체의 CaBER 데이터를 분석하기 위한 계산 프레임워크를 제공한다.