난류 채널 흐름에서 입자 혼합 시뮬레이션

초록

본 연구는 이산 입자와 난류 유체의 상호작용을 DEM과 확률적 난류 모델을 결합해 시뮬레이션하고, 두 종류 입자의 혼합 특성을 분석한다. 혼합 영역은 파이프 길이에 비례해 선형적으로 증가하며, 테일러 마이크로스케일 레이놀즈 수가 약 300 이하일 때는 레이놀즈 수에 민감하게 변하고, 이를 초과하면 혼합 속도가 거의 일정하게 유지된다.

상세 분석

이 논문은 전통적인 유한체적(DEM) 방식에 난류의 미세한 속도 변동을 재현하는 확률적 모델을 결합함으로써, 입자-유체 상호작용을 보다 현실적으로 구현한다. 입자는 구형이며, 입자 간 충돌은 Hertzian 접촉 모델과 마찰을 포함한 Cundall‑Strack 모델로 처리한다. 유체는 평균 흐름을 파이프 중심축에 평행한 포아송 흐름으로 가정하고, 난류의 고유 변동은 Kolmogorov‑Obukhov 스펙트럼을 기반으로 한 랜덤 가우시안 과정으로 재구성한다. 이렇게 얻어진 입자 가속도는 평균 유속과 난류 유동의 순간적 변동을 동시에 반영한다.

시뮬레이션은 두 종류의 입자를 동일한 질량·크기로 설정하고, 초기에는 완전히 분리된 상태에서 파이프 입구에 투입한다. 시간에 따라 입자들의 횡단면 분포를 측정해 혼합 폭을 정의하고, 파이프 길이에 따른 혼합 폭의 변화를 분석한다. 결과는 혼합 폭이 파이프 길이에 대해 거의 직선적인 증가를 보이며, 이는 입자들이 난류 소용돌이와 충돌을 통해 지속적으로 재배열된다는 물리적 메커니즘을 반영한다.

특히 레이놀즈 수(Re)와 테일러 마이크로스케일 레이놀즈 수(R_λ) 사이의 관계를 조사한 결과, R_λ≈300을 임계값으로 하여 두 구간으로 나뉜다. R_λ<300 구간에서는 난류 스케일이 입자 크기와 비슷하거나 더 작아, 입자들이 난류 구조에 강하게 종속되므로 혼합 속도가 Re에 민감하게 변한다. 반면 R_λ>300 구간에서는 난류 스케일이 충분히 커서 입자들이 다수의 독립적인 소용돌이에 노출되며, 이때 혼합 속도는 Re에 거의 의존하지 않는다. 이러한 현상은 Kolmogorov 스케일과 입자 직경 사이의 비율이 혼합 효율을 결정한다는 기존 이론과 일치한다.

또한, 입자 간 마찰계수와 탄성계수의 변화를 시험했을 때, 혼합 폭의 성장률에는 큰 영향을 미치지 않으며, 난류 유동 자체가 혼합을 주도한다는 결론을 도출한다. 이는 산업 현장에서 입자 물성보다는 흐름 조건을 최적화하는 것이 혼합 효율을 높이는 핵심임을 시사한다.