난류 모델의 변동성 탐구

초록

본 연구는 직접수치시뮬레이션(DNS)과 가우시안 필터링을 이용해 서브그리드 스케일 에너지 플럭스를 정량화하고, 필터링된 속도 구배 텐서의 불변량 Q·R에 조건화된 평균 및 확률분포를 분석한다. 클라크 모델은 평균 플럭스를 잘 재현하지만, 동일 (Q,R) 조건에서 실제 플럭스와 모델 예측 사이에 평균과 같은 규모의 강한 변동과 무거운 꼬리를 보인다. 스마그오르스키 모델 역시 비슷한 경향을 보이며, 이러한 결과는 확률적 서브그리드 모델링의 필요성을 강조한다.

상세 분석

본 논문은 고레인즈 수치실험을 기반으로, 가우시안 필터를 적용한 후 필터 스케일 ℓ 에서 정의되는 서브그리드 텐서 τₗᵢⱼ 와 에너지 플럭스 Πₗ = −τₗᵢⱼ ∂ⱼuᵢₗ 을 직접 계산한다. 필터링된 속도 구배 텐서 ∇u 의 두 불변량 Q = −½ Tr(∇u·∇u)와 R = −⅓ Tr(∇u·∇u·∇u) 를 이용해 (Q,R) 공간에 조건부 평균 ⟨Π|Q,R⟩ 를 구하고, 이를 전통적인 클라크 모델(τ≈ℓ²∂ₖuᵢ∂ₖuⱼ)과 스마그오르스키 모델(τ≈−ℓ²|S|S)과 비교한다. 결과는 다음과 같다. 첫째, DNS와 클라크 모델 모두 (Q,R) 평면에서 전방 에너지 전이(양의 Π)와 역전이(음의 Π)의 위치를 유사하게 재현한다. 그러나 정량적으로는 클라크 모델이 평균 플럭스를 전반적으로 과소평가한다. 둘째, 동일 (Q,R) 조건에서 실제 플럭스와 모델 예측 사이의 분포는 평균과 비슷한 규모의 표준편차를 보이며, 특히 양의 꼬리에서 비정규적, 뾰족한 형태의 ‘heavy‑tail’를 나타낸다. 이는 클라크 모델이 1차 ℓ² 항만을 포함하고, 고차 비국소 항(예: ℓ⁴ ∂⁴u 등)을 무시함에 따른 결과로 해석된다. 셋째, 스마그오르스키 모델은 구조적으로 Π≥0 를 보장하지만, 실제 DNS에서는 음의 플럭스가 존재함을 확인한다. 따라서 스마그오르스키 모델은 역전이 현상을 전혀 포착하지 못한다. 마지막으로, 조건부 PDF와 (Q,R) 공간에서의 공동 확률밀도함수(PDF)를 분석한 결과, 플럭스의 비대칭성(skewness)과 높은 첨도(kurtosis)가 두 모델 모두에서 충분히 재현되지 않으며, 이는 기존 결정론적 서브그리드 모델이 작은 스케일의 급격한 변동을 반영하지 못한다는 근거가 된다. 이러한 통계적 특성을 고려하면, 서브그리드 스케일의 랜덤성을 명시적으로 포함하는 확률적 모델링이 필요함을 강력히 시사한다.