무거운 콜모고로프 입자가 관성 난류 연쇄를 억제한다

초록

**
본 연구는 직경이 콜모고로프 길이와 동일한 구형 입자를 고밀도로(밀도비 100–1500, 질량비 0.1–0.6) 분산시켰을 때, 레이놀즈 수 150의 동질·등방성 난류에 미치는 영향을 입자‑해결 직접수치시뮬레이션으로 조사한다. 입자 관성이 커질수록 에너지 스펙트럼이 전통적인 κ⁻⁵ᐟ³ 에서 κ⁻¹ 으로 변형되고, 비선형 에너지 전달이 크게 억제된다. 두 번째 구조함수는 입자 직경보다 큰 거리에서 로그 스케일을 보이며, 입자 주변에서는 축방향 변형과 와류 압축이 강화된다. 입자 군집은 밀도비와 스토크스 수가 증가함에 따라 약화되지만, 남은 군집은 여전히 고변형·저와류 영역을 선호한다.

**

상세 분석

**
본 논문은 콜모고로프 길이와 동일한 직경 D ≈ η 을 갖는 구형 입자를 고밀도(부피분율 Φₚ = 10⁻³)로 분산시킨 경우, 입자 관성(밀도비 Ψₚ = 100–1500, 질량비 Mₚ = 0.1–0.6)이 난류의 에너지 전달 메커니즘에 미치는 영향을 정량적으로 규명한다. 우선, 입자‑해결 DNS(IBM 기반)와 2048³ 격자, Re_λ ≈ 150이라는 충분히 넓은 스케일 분리를 확보한 설정을 사용함으로써, 입자와 가장 작은 난류 구조 사이의 상호작용을 정확히 포착한다.

에너지 스펙트럼 분석에서, Ψₚ가 증가함에 따라 전통적인 Kolmogorov κ⁻⁵ᐟ³ 스케일이 급격히 사라지고, κ⁻¹ 스펙트럼이 지배적인 새로운 정규 상태가 나타난다. κ⁻¹ 스펙트럼은 비선형 전이항 T(k) 가 거의 0에 가까워짐을 의미하며, 이는 입자‑유체 상호작용 Πₚ 와 점성 소산 ε 이 에너지 균형을 주도한다는 것을 시사한다. 구조함수 S₂(r) 에 대해서는 r > D 영역에서 S₂(r) ∝ ln r 형태의 로그 스케일이 관측되며, 이는 입자에 의해 흐름이 급격히 탈동조화(decorrelation)된 결과로 해석된다.

또한, 입자 주변의 속도 구배 텐서 분석을 통해, 입자 표면 근처에서 축방향 변형 S₁₁ 이 크게 증가하고, 회전 성분 Ω 이 억제되는 ‘축압축‑와류 압축’ 현상이 확인된다. 이는 무거운 입자가 주변 유체를 끌어당겨 고변형 영역을 형성하고, 동시에 와류를 억제함으로써 전체 난류 강도를 감소시키는 메커니즘과 일치한다.

입자 군집 현상에 대해서는, 질량비가 커질수록 클러스터의 부피와 입자 수가 감소하고, 군집 간 평균 거리와 클러스터의 차원수가 감소한다. 그럼에도 불구하고, 남은 군집은 높은 변형·낮은 와류 영역을 선호한다는 Maxey‑type 선호 샘플링이 유지된다. 이는 입자 관성이 클수록 원심력에 의해 고변형 영역에 머무르는 경향이 강화된 결과로 볼 수 있다.

마지막으로, 스케일별 에너지 예산을 정량화한 결과, 입자‑유체 상호작용 항 Πₚ(k) 가 전체 스펙트럼에 걸쳐 일정하게 양(energy extraction)이며, 비선형 전이 T(k) 는 거의 0에 수렴한다. 이는 ‘에너지 차단(energy bottleneck)’ 현상이 입자에 의해 직접 대체된다는 새로운 물리적 통찰을 제공한다. 전체적으로, 본 연구는 콜모고로프‑크기 입자가 난류 연쇄를 근본적으로 재구성하고, 고관성 입자 집합이 에너지 흐름을 억제하는 메커니즘을 최초로 수치적으로 입증하였다.

**