정의
μ(I) 유변학은 조밀한 입자(그레인) 흐름을 기술하기 위해 고안된 경험적 유변학 모델이다. 이 모델에서는 전단 응력 σ와 정압 P의 비인 전단 마찰계수 μ를 관성수(I, inertial number)라는 무차원 변수의 함수로 정의한다.
$$ \mu(I) ;=; \frac{\sigma}{P};=;\mu_s ;+; \frac{\Delta\mu}{1 + I_0/I} $$
- μₛ : 정지 마찰계수(정체 상태에서의 μ)
- Δμ : 흐름이 시작될 때 μ가 증가하는 폭 (μ₂ – μₛ)
- I₀ : 전형적인 관성수 스케일
- I : 관성수 = $\frac{\dot\gamma d}{\sqrt{P/\rho_s}}$
- $\dot\gamma$ : 전단 변형률(전단 속도)
- $d$ : 입자 평균 직경
- $\rho_s$ : 입자 물질 밀도
이 식은 입자 간 마찰 및 충돌에 의해 발생하는 전단 저항을 압력에 대한 비율로 표현함으로써, 조밀한 비뉴턴 유동(예: 비산형 흐름, 토양, 모래)에서의 동적 거동을 간단히 설명한다.
역사·배경
- 초기 연구
- 1980‑1990년대, 입자 흐름을 연속체 모델로 해석하려는 시도가 늘어났지만, 전통적인 뉴턴 유체 모델은 입자 간 충돌과 마찰을 적절히 반영하지 못했다.
- 관성수 개념 도입
- 2004년 Jop, Forterre, and Pouliquen이 제안한 μ(I) 관계는 입자 흐름을 압력‑의존적인 전단 마찰계수로 기술하는 최초의 성공적인 경험법칙이었다.
- 후속 확장
- 2010년대 이후, 구체적인 물성(입자 형태, 입도 분포, 습도 등)과의 연관성을 포함한 μ(I) – ϕ(I)(밀도 ϕ) 모델, 그리고 비관성수 Bagnold number 등을 결합한 다중 변수 모델이 개발되었다.
핵심 식·요소
| 변수 | 의미 | 정의 |
|---|---|---|
| $\mu(I)$ | 전단 마찰계수 | $\mu = \sigma/P$ |
| $I$ | 관성수 | $I = \dot\gamma d / \sqrt{P/\rho_s}$ |
| $\mu_s$ | 정지 마찰계수 | 흐름이 정지할 때의 μ |
| $\mu_2$ | 고속 흐름 마찰계수 | $I \to \infty$에서의 μ |
| $\Delta\mu = \mu_2 - \mu_s$ | 마찰계수 변화량 | |
| $I_0$ | 전이 관성수 | μ가 $\mu_s$와 $\mu_2$ 사이에 전이하는 기준값 |
전단 응력은 압력과 μ(I)의 곱으로 주어지며, 전단 변형률은 관성수와 역으로 연관된다. 따라서 입자 흐름을 전산유체역학(CFD) 혹은 입자 기반 시뮬레이션(DEM)과 연계해 계산할 때, μ(I) 관계를 본질적인 constitutive law 로 사용한다.
실험적 검증
| 실험 방법 | 주요 결과 | 참고 논문 |
|---|---|---|
| 평면 전단 셀(Plane shear cell) | 전단 응력‑압력 비가 관성수에 따라 μ(I) 함수와 잘 맞음 | Jop et al., Nature 2006 |
| 경사면 흐름(Inclined plane flow) | 유동층 두께와 흐름 속도가 관성수에 의해 예측 가능 | Forterre & Pouliquen, J. Fluid Mech. 2008 |
| 원통형 콜럼버스(Annular shear) | 입도 분포와 입자 형태가 Δμ, I₀에 미치는 영향 보고 | GDR MiDi, Europhys. Lett. 2004 |
| 3D X‑ray CT + PIV | 내부 입자 속도장과 압력장이 동시에 측정되어 μ(I) 모델 검증 | Reddy et al., Phys. Rev. Lett. 2020 |
실험 결과는 입자 간 마찰계수, 탄성계수, 입자 형상 등에 따라 μₛ, Δμ, I₀ 값이 변한다는 점을 강조한다.
주요 적용 분야
- 산사태·토류·지반공학
- 사면 붕괴와 흐름을 예측하는 데 사용되며, 특히 비포화 상태의 모래 흐름 모델링에 적합.
- 산업 공정
- 제분, 제약, 식품공정 등 입자 혼합·이송 시스템에서 압력‑전단 관계를 설계에 활용.
- 천체 물리학
- 소행성 표면의 미세입자 흐름, 유성우 충돌 후 입자 재분포 등을 설명.
- 시뮬레이션
- CFD‑DEM 연계 모델에서 연속체 부분에 μ(I) 유변학을 적용해 계산 효율성을 높인다.
한계점·비판
| 한계 | 내용 |
|---|---|
| 희박·입자간 거리 큰 흐름 | 관성수 정의가 압력에 의존하므로, 압력이 거의 없는 경우(예: 자유 낙하)에는 적용 불가. |
| 습윤·점착성 입자 | 물 또는 점착제 존재 시, μ(I) 모델에 추가적인 캡시티브 힘이 필요. |
| 비균질 입도·비구형 입자 | Δμ, I₀ 파라미터가 입도 분포와 형상에 따라 크게 변해, 단일 파라미터로는 설명이 부족. |
| 전단 경계 조건 | 경계에서 발생하는 슬립·바운더리 레이어가 μ(I) 관계와 일치하지 않을 경우가 있음. |
| 동적 전이 영역 | 흐름이 정지→흐름, 혹은 흐름→정지 전이 시 히스테리시스가 나타나 μ(I) 식에 추가적인 히스테리시스 파라미터가 요구됨. |
이러한 한계는 μ(I)‑ϕ 모델, 비관성수(Bagnold number) 기반 확장, 혹은 시공간 가변 파라미터를 도입해 보완하려는 연구가 진행 중이다.
관련 연구·발전 동향
- μ(I)‑ϕ 모델 (Barker et al., Granular Matter 2019): 관성수와 동시에 입자 부피밀도 ϕ를 함수화하여 압축성 효과 포함.
- 다중 스케일 연계: DEM으로 미세입자 충돌을 직접 시뮬레이션하고, 얻은 μ(I) 파라미터를 CFD에 연동하는 Hybrid CFD‑DEM 접근법.
- 머신러닝 기반 파라미터 추정: 실험 데이터에서 μₛ, Δμ, I₀를 역학적으로 추정하기 위해 Gaussian Process Regression 등 활용.
- 습윤·점착 입자 모델: 물리적 점착력과 표면 장력을 고려한 μ(I, J)(J는 점착수) 형태 확장 연구가 활발.
참고문헌 (주요)
- Jop, P., Forterre, Y., & Pouliquen, O. (2006). A constitutive relation for dense granular flows. Nature, 441, 727‑730.
- Forterre, Y., & Pouliquen, O. (2008). Flows of dense granular media. Annual Review of Fluid Mechanics, 40, 1‑24.
- GDR MiDi. (2004). On dense granular flows. Eur. Phys. J. E, 14, 341‑364.
- Reddy, K. et al. (2020). X‑ray tomography of granular shear flows. Physical Review Letters, 124, 018001.
- Barker, L., Gray, J. M. N. T., & Hutter, K. (2019). A coupled μ(I)–ϕ rheology for granular materials. Granular Matter, 21, 70.
요약
μ(I) 유변학은 압력‑의존적인 전단 마찰계수를 관성수 I의 함수로 정의함으로써, 조밀한 입자 흐름을 간단하면서도 실험적으로 검증된 형태로 기술한다. 다양한 공학·과학 분야에 적용되고 있으나, 희박·습윤·비균질 입자 등 특수 상황에 대해서는 추가적인 모델링이 필요하다.