[논문 리뷰] Longitudinal and transverse velocity scaling exponents from merging of the Vortex filament and Multifractal models

K. P. Zybin, V. A. Sirota|arXiv (Cornell University)|2012. 04. 06.

Fluid Dynamics and Turbulent Flows참고 문헌 1인용 수 3

한 줄 요약

이 논문은 비틀림 필라멘트 모델과 다중분포(MF) 가설을 융합하여 난류에서 종방향 및 횡방향 속도 스케일링 지수 간의 관측된 차이를 통합적으로 설명한다. 물리적 추론과 특이성 스펙트럼 D(h)의 포물선 형태를 사용하여 피팅 매개변수 없이 ζₙ∥ 및 ζₙ⊥의 해석적 표현을 유도하며, 고차수 모멘트 및 종방향 대 횡방향 구조함수 간의 상이한 스케일링 행동에 대해 DNS 및 실험 데이터와 뛰어난 일치를 이룬다.

ABSTRACT

We suggest a simple explanation of the difference between transverse and longitudinal scaling exponents observed in experiments and simulations. Based on the Vortex filament model and Multifractal conjecture, we calculate both scaling exponents for any n without any fitting parameters and ESS anzatz. The results are in very good agreement with the data of simulations.

연구 동기 및 목표

• 모델링 및 실험에서 관측된 종방향 및 횡방향 속도 구조함수 스케일링 지수 간 오랫동안 지속된 모순을 해결하기 위해.
• 비틀림 필라멘트의 역학 및 기하학적 성질에 뿌리를 둔 물리적 설명을 제공하기 위해.
• 비틀림 필라멘트 모델과 다중분포 가설을 융합하여 스케일링 지수를 일관되게 유도하기 위해.
• ESS 가설이나 피팅 매개변수에 의존하지 않고 임의의 n에 대해 ζₙ∥ 및 ζₙ⊥을 예측하기 위해.
• 종방향과 횡방향에서 ζₙ∥ 및 ζₙ⊥가 다른 이유, 특히 고차수에서의 차이를 축축한 필라멘트 형태(축대칭 대 강한 굽은 형태)에 기반해 설명하기 위해.

제안 방법

• 비틀림 필라멘트 모델을 사용하여 난류 속도 구조함수를 모델링하며, 이 때 회전율은 얇은 필라멘트에 집중된다.
• 지역적 스케일링 행동을 기술하기 위해 특이성 스펙트럼 D(h)를 사용하는 다중분포 가설을 적용하며, P ∝ l^{3−D(h)}이다.
• 큰 변동성 정리(Large Fluctuations Theorem)를 적용하여 구조함수를 h에 대한 적분으로 표현한다: ⟨Δvⁿ⟩ ∝ ∫ l^{nh} l^{3−D(h)} dμ(h).
• D(h)를 오목한 포물선 D(h) = 3 − b(h − h₀)²로 가정하여 ζₙ = min_h (nh + 3 − D(h))를 통해 스케일링 지수를 유도한다.
• 필라멘트 기하학에 기반해 횡방향과 종방향 기여를 구분하며, 각각 D⊥(h)와 D∥(h) 함수를 할당한다.
• 경계 조건 D⊥(0) = 3, D⊥(h₀) = 3, D∥(0) = 0, D∥(h₀) = 3을 사용하여 D(h)의 매개변수를 구하고, 각 경우에 대해 별개의 h₀ 및 b를 도출한다.

실험 결과

연구 질문

• RQ1 이sovotropic 및 균일한 난류에서 종방향 및 횡방향 속도 구조함수 스케일링 지수가 다른 이유는 무엇인가요?
• RQ2 비틀림 필라멘트 내에서 종방향 대 횡방향 속도 증분의 상이한 스케일링 행동을 유도하는 물리적 메커니즘은 무엇인가요?
• RQ3 ESS 가설이나 피팅 매개변수 없이 ζₙ∥ 및 ζₙ⊥ 간의 차이를 해석적으로 도출할 수 있을까요?
• RQ4 비틀림 필라멘트의 기하학적 성질—특히 축대칭성 대 강한 굽음—은 스케일링 지수에 어떻게 영향을 미치나요?
• RQ5 종방향 및 횡방향 구조함수의 관측된 스케일링 지수를 재현할 수 있는 특이성 스펙트럼 D(h)의 함수 형태는 무엇인가요?

주요 결과

• 모델은 n ≤ 16.3일 때 ζₙ∥ = 0.367n − 1.12×10⁻²n², n > 16.3일 때 ζₙ∥ = 3으로 예측하며, 실험 및 DNS 데이터와 일치한다.
• 모델은 n ≤ 10.2일 때 ζₙ⊥ = 0.391n − 1.91×10⁻²n², n > 10.2일 때 ζₙ⊥ = 2로 예측하며, Benzi 및 Gotoh의 DNS 결과와 뛰어난 일치를 보인다.
• 스케일링 지수의 차이는 축대칭 필라멘트와 강하게 굽은 필라멘트의 별개의 필라멘트 형태에 기인한다: 축대칭 필라멘트는 횡방향 구조함수에 기여하고, 강하게 굽은 필라멘트는 종방향 구조함수를 지배한다.
• n=2일 때 정확한 조건 ζ₂ = ζ₂⊥를 만족하며, 작은 편차(1.6×10⁻²)가 존재하지만, 이는 고차수 D(h) 항을 추가함으로써 수정 가능하며, 이는 시각적 피팅에 영향을 주지 않는다.
• 횡방향 및 종방향의 경우에 대해 도출된 D(h) 함수는 물리적으로 구별된다: D⊥(0) = 3 및 D∥(0) = 0으로, h=0에서의 다른 스케일링 행동을 반영한다.
• 모델은 비틀림 필라멘트 역학 및 다중분포 스케일링의 물리 원리에 기반하여 피팅 매개변수 없이 모든 n에 대해 ζₙ∥ 및 ζₙ⊥의 완전한 해석적 프레임워크를 제공한다.