수직 파이프에서 입자 흐름과 막힘 전이

초록

본 연구는 2차원 수직 파이프 내에서 입자들이 막히는 현상을 실험과 DEM 시뮬레이션으로 조사하고, 파이프‑입자 직경비 (D/d)와 마찰계수 (μ)의 조합이 막힘 확률 (J)에 미치는 영향을 규명한다. D/d ≈ 2.45에서 급격한 막힘‑흐름 전이가 나타나며, μ가 클수록 전이점이 낮아진다. 세 입자 아치 형성의 기하학·힘·토크 균형을 분석한 이론식이 실험·시뮬레이션 결과와 일치한다.

상세 분석

이 논문은 기존 3차원 원통형 파이프에서 보고된 비단조성(clogging probability)의 변화를 2차원 평면 파이프에서도 재현함으로써, 입자 배열의 차원적 제약이 막힘 현상에 미치는 영향을 명확히 한다. 핵심 변수는 파이프 직경 D와 입자 직경 d의 비율 (D/d)와 입자·벽 마찰계수 μ이며, 두 변수는 독립적으로가 아니라 상호보완적으로 작용한다. 실험에서는 스테인리스 스틸 원판(d = 10 mm, 12 mm)을 사용해 N = 4050개의 입자를 파이프 상단에 적재하고, 바인더 클립을 제거해 중력에 의해 흐르게 한다. 흐름이 완전히 종료되면 비클깅, 일부 입자가 남으면 클깅으로 판정한다. 각 조건에서 100300번의 반복 실험을 통해 J = N_c/N_t를 구한다.

DEM 시뮬레이션은 LMGC90 기반 접촉역학을 이용해 입자를 강체·비탄성으로 모델링하고, μ_p = μ_w = μ로 설정한다. 시뮬레이션 파라미터는 시간 단계 Δt = 5×10⁻⁴ s, 최대 2000 스텝(≈1 s)이며, 입자 속도가 10⁻⁶ m/s 이하로 감소하면 클깅으로 판단한다. 시뮬레이션 결과는 실험과 전반적으로 일치하지만, D/d > 2.73에서 약간의 차이가 나타난다. 이는 실제 입자 표면 거칠기와 시뮬레이션 마찰 모델 차이에서 기인한다.

J와 D/d의 관계는 비단조적이며, D/d ≈ 2와 D/d ≈ 2.45~2.8 구간에서 J가 최소값을 보인다. 특히 D/d ≈ 2.45에서 급격히 상승하는 전이점(D/d)_c는 μ에 따라 이동한다. μ가 증가할수록 (D/d)_c는 감소해, 마찰이 큰 경우 더 작은 파이프에서도 막힘이 발생한다는 점을 보여준다. 이는 마찰이 토크와 수평 힘을 보강해 아치가 더 쉽게 안정화되기 때문이다.

아치 형성 메커니즘을 해석하기 위해 저자들은 3입자 아치를 모델링한다. 아치의 기하학적 제약은 입자 중심 간 거리와 파이프 벽 사이의 접촉을 고려해 θ = arccos(0.5 D/d − 0.5) 로 표현된다. 실험·시뮬레이션에서 측정된 θ 분포는 D/d가 작을수록 좁고, D/d가 2.5 이상이면 폭넓어지며, D/d = 3에 가까워질수록 아치 형성 확률이 급격히 감소한다.

힘·토크 평형식은 다음과 같이 전개된다.

1. 수직 평형: F₂ sinθ + G = f₁ + f₂ cosθ
2. 수평 평형: F₂ cosθ + f₂ sinθ = F₁
3. 마찰 제한: f₁ ≤ μ F₁, f₂ ≤ μ F₂
4. 토크 평형: f₁ = f₂

이들 식을 결합하면 μ ≥ 1 − cosθ sinθ 이 도출되고, θ와 D/d의 관계를 대입하면 최종적인 전이 조건식
D/d ≥ max{ (3 − μ²)/(1 + μ²), 2 }
이 얻어진다. μ = 0.4일 때 D/d ≥ 2.45가 되며, 이는 실험·시뮬레이션에서 관찰된 (D/d)_c와 일치한다. 따라서 마찰이 충분히 크면 기하학적 제약보다 힘·토크 균형이 주된 막힘 메커니즘이 된다.

또한, 입자 간 접촉 수(N_f)와 벽 접촉 수를 D/d에 따라 조사한 결과, 2 < D/d < 2.73 구간에서는 벽에 접촉하는 입자가 한 개뿐이지만, D/d > 2.73에서는 두 개가 동시에 접촉해 토크 균형을 맞추기 위해 더 큰 전단력이 필요한다. 이 때문에 J가 급격히 감소한다.

결론적으로, 저자들은 마찰계수와 파이프‑입자 비율이 결합된 2차원 파이프에서의 막힘‑흐름 전이를 정량적으로 설명하는 이론적 프레임워크를 제시하고, 실험·시뮬레이션 결과와의 높은 일치를 통해 모델의 타당성을 검증하였다. 이 연구는 파이프 설계, 입자 운반 시스템, 핵연료 운송 등 다양한 산업 분야에서 막힘 위험을 예측하고 최소화하는 데 실질적인 가이드를 제공한다.