MPDFT 기반 고해상도 DNS에서 에너지 소산률 PDF 분석 및 부분 영역 분리

초록

4096³ DNS 데이터를 이용해 에너지 소산률 확률밀도함수(PDF)를 MPDFT 이론으로 정밀 분석하고, 전체 영역에서 추출한 w‑PDF를 기준으로 최대·최소 와전류(enstrophy) 영역(512³)에서 만든 max‑PDF와 min‑PDF를 각각 확대·축소하여 꼬리와 중심 부분의 일치를 확인하였다.

상세 분석

본 논문은 다중프랙탈 확률밀도함수 이론(MPDFT)을 이용해 완전 발달 난류의 에너지 소산률 PDF를 고해상도(4096³) 직접수치시뮬레이션(DNS) 데이터에서 추출하였다. MPDFT는 코히런트(비선형) 흐름과 비코히런트(점성 소산) 두 기여를 구분하고, PDF를 ‘꼬리 부분’과 ‘중심 부분’으로 나눈다. 꼬리 부분은 코히런트 구조가 지배한다고 가정하여, α‑지수의 레니 또는 HCT 최대엔트로피 분포 P⁽ⁿ⁾(α) 를 통해 직접 계산한다. 중심 부분은 코히런트와 비코히런트가 혼합되므로, 경험적 Tsallis‑형 함수 e^{‑g(ξ)} 를 도입해 연결점 ξ* (또는 ε* )에서 꼬리와 매끄럽게 이어지도록 한다.

핵심 파라미터는 인터미턴시 µ, 엔트로피 지수 q, 스케일링 파라미터 δ 이며, 저자들은 δ를 임의로 선택해도 물리량은 변하지 않도록 (1‑q)·ln δ 조합에만 의존하도록 설계하였다. DNS에서 얻은 전체 영역 PDF(w‑PDF)를 기준으로, 512³ 크기의 부분 영역을 최대 와전류(enstrophy)와 최소 와전류 영역으로 나누어 각각 max‑PDF와 min‑PDF를 만든다.

분석 결과, max‑PDF의 꼬리 부분은 w‑PDF에 적절한 확대 계수 χ (>1)를 곱하면 정확히 겹친다. 확대 계수와 연결점 ξ* 은 w‑PDF의 이론적 연결점과 일치하여, 부분 영역에서 코히런트 구조가 전체와 동일한 스케일링을 유지함을 보여준다. 반면 min‑PDF의 중심 부분은 w‑PDF를 동일한 축소 계수 σ (<1)로 스케일링하면 거의 완벽히 재현된다. 이는 최소 와전류 영역이 비코히런트(점성) 기여가 상대적으로 크게 작용함을 의미한다.

또한, 파라미터 µ=0.345, α₀=1.20, X=0.411 등이 δ=2^{1/4}, 2^{1/2}, 2 세 경우 모두 동일하게 추정되어, MPDFT가 스케일링 자유도에 강인함을 확인한다. 논문은 PDF의 꼬리와 중심을 각각 다른 물리적 메커니즘으로 해석함으로써, 난류의 다중프랙탈 구조와 점성 소산의 상호작용을 정량적으로 분리하는 새로운 방법론을 제시한다.