중성 부력 원통의 2차원 전단 흐름 내 움직임에 대한 수치 연구

초록

본 논문은 뉴턴 유체의 2차원 전단 흐름 속에서 중성 부력을 가진 원형·타원형 실린더의 운동을 직접 수치 시뮬레이션으로 조사한다. 시뮬레이션 결과는 기존 이론·실험·수치 자료와 비교해 검증했으며, 레이놀즈 수와 초기 위치에 따라 입자 중심이 두 개의 평형선(중심선·벽 근처) 중 하나로 이동한다는 것을 밝혔다.

상세 분석

이 연구는 유체-구조 상호작용을 정밀히 포착하기 위해 Immersed Boundary Method(IBM)를 기반으로 한 직접 수치 시뮬레이션(DNS)을 채택하였다. 유체는 비압축성 뉴턴 유체로 가정하고, 연속 방정식과 Navier‑Stokes 방정식을 비정규화하여 입자 직경·채널 높이·유체 점성 등으로 스케일링하였다. 입자는 중성 부력을 가지고 있어 중력·부력 효과가 상쇄되며, 따라서 관성·전단력만이 입자 운동을 지배한다.

시뮬레이션은 두 종류의 입자 형상(완전 원형, 장축 비율 2:1인 타원형)을 대상으로, 레이놀즈 수(Re) 0.1200, 채널-입자 비율(K) 0.50.8 범위에서 수행되었다. 초기 위치를 중심선(채널 중간)과 벽에서 일정 거리 떨어진 두 경우로 나누어, 입자 중심의 y‑좌표와 회전 각속도 ω(t)를 시간에 따라 추적하였다.

검증 단계에서는 Jeffery’s orbit(무한 희석, 무한 채널)과 Breugem et al. (2006)의 실험 데이터를 재현함으로써 수치 정확성을 확인하였다. 특히, 정규화된 각속도 Ω* = ω·a/U_max와 Re 사이에 Ω* ∝ Re^−0.5 형태의 전력 법칙이 관찰되어 기존 문헌과 일치하였다.

주요 물리적 메커니즘을 해석하면, 입자가 중심선에서 벗어나면 전단에 의해 발생하는 슬립 속도(v_slip = U_particle − U_fluid)가 비대칭적인 압력 분포를 유도한다. 이 압력 차는 Magnus force(회전으로 인한 양력)와 결합해 입자를 한쪽 벽으로 끌어당긴다. 동시에, 입자와 벽 사이의 유체가 압축되면서 발생하는 유동 저항(벽 반발력)이 입자를 다시 중심으로 복귀시키는 역할을 한다. Re가 낮을 때는 점성력이 지배적이어서 중심선이 전역적인 안정점이 된다. 그러나 Re가 일정 임계값(≈ 30–50) 이상으로 상승하면, 관성 효과가 강화되어 Magnus force와 슬립에 의한 비대칭 압력이 벽 반발력을 압도하고, 입자는 새로운 안정점(벽 근처)으로 이동한다.

형상 효과에서도 차이가 나타났다. 타원형 입자는 회전축이 길이 방향과 일치하면서 회전 각속도가 원형보다 크게 억제되며, 따라서 Magnus force가 약해진다. 결과적으로 동일 Re와 K 조건에서 타원형 입자는 원형보다 중심선에 더 오래 머무르며, 평형 위치 전이가 더 높은 Re에서 발생한다.

이러한 결과는 미세유체공학, 혈류 내 세포 이동, 그리고 입자 기반 혼합 공정 등에서 입자의 위치 제어와 회전 제어 전략을 설계하는 데 직접적인 시사점을 제공한다. 특히, 레이놀즈 수와 초기 위치를 조절함으로써 원하는 평형선(중심선 또는 벽 근처)으로 입자를 유도할 수 있다는 점은 채널 설계와 운영 조건 최적화에 활용 가능하다.