시간·공간 상관 난수 모델로 보는 입자와 입자쌍 확산

초록

본 논문은 동질·등방성 난류에서 입자와 입자쌍의 확산을 기술하기 위해, 시간적으로 상관된 난수 과정에 공간적 상관을 부여하는 Heisenberg‑유사 해밀토니안을 도입한 새로운 모델을 제시한다. 단일 입자에 대해 짧은 시간에서는 $t^{2}$의 탄성(ballistic) 거동을, 장시간에서는 정상 확산($t$)으로 전이함을 재현한다. 입자쌍의 평균 제곱 거리 역시 초기 $t^{2}$, 장기 $t$ 거동을 보이며, 특정 파라미터 구간에서 Richardson의 $t^{3}$ 법칙이 나타나는 중간 스케일을 확인한다. 또한 관성 입자의 질량 효과가 확산에 미치는 영향을 조사한다.

상세 분석

이 연구는 난류 입자 운반 현상을 확률적 모델링으로 접근하면서, 기존의 시간‑only 상관 난수 모델이 갖는 공간적 비연속성을 보완한다는 점에서 의미가 크다. 저자들은 유체 입자들의 위치와 속도를 각각 독립적인 Ornstein‑Uhlenbeck 과정으로 기술하고, 이들 사이에 Heisenberg‑형 해밀토니안 $H=-\sum_{i<j}J_{ij},\mathbf{u}{i}\cdot\mathbf{u}{j}$ 를 도입해 인접 입자 쌍 간의 에너지 교환을 강제한다. 여기서 $J_{ij}$는 거리 $r_{ij}$에 따라 감소하는 커플링 상수로, $J_{ij}\propto r_{ij}^{-p}$ 형태를 취한다. 이 구조는 물리적으로는 소용돌이 간 상호작용을 모사하며, 수학적으로는 공간 상관 함수를 직접 제어할 수 있게 한다.

시간 상관은 각 입자에 대한 OU 과정의 고유 시간 $\tau$ 로 정의되며, 이는 Kolmogorov 스케일에 맞추어 조정된다. 모델 파라미터 $p$, $\tau$, 그리고 난류 강도 $\epsilon$ 를 적절히 선택하면, 단일 입자의 평균 제곱 변위(MSD) $\langle\Delta x^{2}(t)\rangle$ 가 초기 $t^{2}$ (탄성) 구간을 지나, $t^{1}$ (정상 확산) 구간으로 부드럽게 전이함을 재현한다. 이는 실험 및 DNS 결과와 정량적으로 일치한다.

입자쌍에 대해서는 평균 제곱 거리 $\langle r^{2}(t)\rangle$ 를 추적한다. 짧은 시간에서는 두 입자가 동일한 흐름을 공유하므로 $r^{2}\sim t^{2}$ 를 보이며, 장기적으로는 독립적인 확산으로 돌아가 $r^{2}\sim t$ 로 수렴한다. 가장 흥미로운 것은 중간 스케일에서 $r^{2}\sim t^{3}$, 즉 Richardson 법칙이 나타나는 구간이 관측된 점이다. 이 현상은 $p$ 값이 2에 가까울 때, 즉 장거리 상관이 충분히 강할 때 두드러지며, 이는 전통적인 K41 이론의 스케일 간격을 초월하는 비선형 확산 메커니즘을 모델이 포착함을 의미한다.

관성 입자(밀도 차이와 크기 때문에 유체와 다른 응답을 보이는 입자)를 포함시키기 위해, 저자들은 Maxey‑Riley 방정식의 간소화 형태를 도입하고, 입자 질량 비 $\beta$ 와 응답 시간 $\tau_{p}$ 를 추가 파라미터로 설정한다. 결과적으로 관성 입자는 초기 탄성 구간이 짧아지고, 중간 스케일에서 $t^{3}$ 거동이 억제되며, 장기 확산 속도가 유체 입자보다 크게 변한다는 점을 확인한다. 이는 실제 환경(예: 대기 입자, 해양 플랑크톤)에서 관성 효과가 확산을 어떻게 변형시키는지를 정량적으로 설명한다.

수치 실험에서는 10⁴개의 입자를 3차원 주기적 경계조건 하에 시뮬레이션하고, 통계적 수렴을 위해 수천 개의 초기 조건을 평균한다. 모델의 민감도 분석을 통해 $p$, $\tau$, $\beta$ 의 변화가 MSD와 $r^{2}$ 스케일에 미치는 영향을 체계적으로 조사한다. 전반적으로 모델은 기존 DNS와 실험 데이터와 비교했을 때, 특히 중간 스케일에서 Richardson 법칙을 재현하는 능력에서 우수함을 보인다.

이 논문의 주요 공헌은 (1) 공간·시간 상관을 동시에 고려한 확률적 난류 모델을 제시, (2) 단일 입자와 입자쌍의 전이 스케일을 정량적으로 설명, (3) 관성 입자에 대한 확산 변형 메커니즘을 분석, (4) Richardson 스케일을 파라미터 공간에서 명시적으로 연결시킨 점이다. 다만, 해밀토니안의 커플링 함수 형태가 경험적으로 선택되었으며, 복잡한 비동질 난류(예: 전단 흐름)에는 직접 적용이 어려울 수 있다는 한계가 있다. 향후 연구에서는 커플링 함수를 물리적 측정값에 기반해 보정하고, 비등방성 및 경계 효과를 포함시키는 확장이 기대된다.