미세역학과 거시유변학의 상호작용 베시클 현탁액 점도 비단조성

초록

본 연구는 격자볼츠만-침투경계법을 이용해 유동 내 베시클의 미세운동과 매크로 점도 사이의 연관성을 조사한다. 점도 대비(내부/외부 점도비) Λ가 변할 때, 베시클은 탱크‑트레딩(TT)에서 텀블링(TB)으로 전이하며, 이 전이점에서 전체 현탁액의 내재 점도 η가 최소값을 보이는 비단조적 거동을 확인한다. 농도, 캡릴러리 수(막힘 강도), 팽창도 Δ 등을 변화시켜도 이 현상이 유지되지만, 고농도·고팽창 경우에는 효과가 약해진다.

상세 분석

이 논문은 2차원 격자볼츠만(LBM)과 침투경계법(IBM)을 결합한 수치 모델을 구축하여, 점도 대비 Λ가 1‒20 범위에서 베시클 현탁액의 점성 거동을 정량적으로 분석한다. 먼저, 단일 베시클(극히 희석된 경우)에서 TT와 TB 상태를 구분하기 위해 기울기 각 θ를 관찰한다. TT 상태에서는 θ가 일정하게 유지되며, 이에 따라 유동에 가해지는 전단 응력이 일정해 내재 점도 η(t)가 시간에 따라 거의 변하지 않는다. 반면 TB 상태에서는 θ가 주기적으로 변동하고, 이와 동조하여 벽면 전단 응력 σ_xy도 주기적으로 변하므로 η(t)도 진동한다. 이러한 미시적 동역학 차이가 매크로 점도 특성에 직접적인 영향을 미친다.

Λ가 증가하면 TT 영역에서는 베시클이 흐름 방향에 더 정렬되어 유동 저항이 감소하고, η는 감소한다. Λ가 임계값 Λ_c≈7.8(본 연구의 제한조건 χ=0.2) 에 도달하면 TT‑TB 전이가 일어나며, η는 최소값을 기록한다. 전이 이후 TB 영역에서는 베시클이 회전하면서 유동에 더 큰 교란을 일으키고, 내부·외부 점도 차이가 커질수록 전단 손실이 증가해 η가 다시 상승한다. 이 비단조적 η(Λ) 곡선은 Danker‑Misbah 이론과 Ghigliotti 등(2D 경계적분법) 결과와 정성적으로 일치한다.

다음으로 농도 φ를 7.5%, 15.1%, 22.6% 로 변화시킨 다베시클 시스템을 조사한다. φ가 증가하면 베시클 간 충돌 빈도가 높아지고, 충돌 시 각 베시클의 θ가 순간적으로 최대값을 취한다. 이는 전단 저항을 일시적으로 증가시켜 η(t)의 변동폭을 확대한다. 평균 η는 φ가 증가함에 따라 전반적으로 상승하지만, Λ에 대한 비단조성은 여전히 유지된다. 다만 고농도에서는 TT‑TB 전이점 주변에서 η의 최소값이 완만해져, 실험적으로 관측하기 어려운 수준으로 완화된다.

캡릴러리 수 Ca(막탄성)와 팽창도 Δ 역시 η(Λ) 곡선에 영향을 미친다. Ca가 작을수록 베시클이 쉽게 변형되어 TT‑TB 전이가 낮은 Λ에서 발생하고, 이에 따라 η 최소점이 왼쪽으로 이동한다. 반대로 Ca가 크면 베시클이 강직해 전이점이 높은 Λ로 이동한다. Δ가 1에 가까운 거의 구형 베시클(2D에서는 원형)에서는 형태 변화가 제한되어 TT‑TB 전이가 거의 사라지고, η는 Λ에 대해 거의 일정한 값을 보인다. 이는 기존 연구에서 보고된 “점도 대비에 대한 단조적 증가” 현상과 일맥상통한다.

또한, 본 연구는 Reynolds 수 Re=0.5(비정상적인 저속 흐름)와 제한된 채널(χ=0.2)에서도 비단조적 η(Λ) 현상이 재현됨을 보여, 이전 연구(Lamura‑Gompper)의 부정적 결과가 수치적 구현 차이(예: 충돌 모델, 경계 조건)에서 비롯될 가능성을 제시한다. 전체적으로, 미세동역학(베시클의 회전·변형)과 매크로 점도 사이의 직접적인 연결 고리를 명확히 규명했으며, 농도·탄성·팽창도 등 실험적 파라미터가 이 연결을 어떻게 조절하는지 체계적으로 제시한다.