무작위 등방성 텐서의 새로운 확률 항등식

초록

본 논문은 3차원 무작위 2차 텐서장에 대해, 회전 대칭성으로부터 유도되는 다섯 개의 비자명한 확률 항등식을 제시한다. LDU 분해와 주요 소행렬식(minor)을 이용해 식을 전개하고, 축대칭 경우에도 적용 가능한 변형 형태를 제시한다. 직접 수치 시뮬레이션(DNS)으로 확인한 결과, 등방성 난류와 채널 흐름에서 항등식이 1에 수렴함을 보이며, 벽 근처에서는 축대칭 항등식이 더 잘 유지됨을 확인한다.

상세 분석

논문은 먼저 무작위 텐서 A의 양의 정부호 이차형식 Γ=AᵀA를 LDU 분해(Γ=ZᵀD²Z) 형태로 표현한다. 여기서 D=diag(D₁,…,D_d)는 상삼각 행렬 Z와 연관된 대각 성분이며, D_i²는 Γ의 k차 주요 소행렬식 s_k와 s_{k‑1}의 비(s_k/s_{k‑1})로 정의된다(식 3). 중요한 점은 D_i가 고유값이 아니라 선택된 좌표계에 의존한다는 점이다.

다음으로, 텐서 A가 회전군 O(d) 아래에서 확률분포 ρ(A)가 불변일 때, 임의의 정수 집합 {m_i}에 대해 ⟨∏i D_i^{m_i+1}⟩가 인접한 두 지수 m_p, m{p+1} 사이의 교환에 대해 대칭임을 보인다(식 5). 이는 SO(2) 부분군에 대한 적분이 갖는 비자명한 대칭성에서 비롯된다. 전부 대칭군 O(d)까지 확대하면, 모든 순열에 대해 대칭성이 유지된다.

이 대칭성을 이용해 m_i=−i를 대입하면, ⟨∏i D_i^{−i}⟩=1이라는 기본 항등식(식 7)을 얻는다. D_i를 s_k로 다시 쓰면, 주요 소행렬식들의 곱에 대한 기대값이 1이 되는 식(8)이 도출된다. d=3인 경우, 3!−1=5개의 서로 다른 순열에 대응하는 다섯 개의 비자명한 항등식이 표 1에 정리된다. 각 항등식은 ⟨s₁ s{2}^{‑1} s₃⟩, ⟨s_{2}^{‑1} s₁ s₂⟩ 등으로 나타나며, 모두 1에 수렴한다.

축대칭(axial symmetry) 경우에는 회전군이 한 축을 중심으로 제한되므로, 표 1의 일부 항등식만 유지된다. 좌표축을 회전시켜 새로운 좌표계에서 동일한 항등식을 재작성하면, 주요 소행렬식이 아닌 일반적인 2‑rank 소행렬식(m_{ij})이 등장한다. 이는 표 2에 정리된 바와 같이, 축대칭을 갖는 텐서에 대해 두 개의 연속적인 1‑파라미터 항등식 군이 존재함을 의미한다.

실제 검증을 위해 저자들은 Johns Hopkins Turbulence Database의 DNS 데이터를 사용했다. 등방성 난류와 고레시스수 채널 흐름에서 속도 구배 텐서 A_{ij}=∂u_i/∂x_j를 계산하고, Γ=AᵀA를 구성한 뒤 표 1·2의 좌변을 공간 평균(또는 시간 평균)으로 추정하였다. 결과는 그림 1에 나타나듯이, 등방성 영역에서는 모든 평균값이 1에 매우 가깝게 수렴하고, 벽에 가까워질수록 일부 항등식은 크게 벗어나지만 축대칭 항등식은 상대적으로 작은 편차를 보였다. 또한 인위적으로 상수 전단 A_{12}=a를 추가한 실험에서도, 항등식들의 편차가 전단 강도에 민감하게 반응함을 확인하였다. 이는 제시된 항등식이 통계적 등방성·축대칭을 정량적으로 검증하는 강력한 지표임을 시사한다.

이러한 결과는 기존의 단순한 2‑점 차분 통계(예: ⟨δu_i δu_i⟩/3)보다 훨씬 풍부한 구조 정보를 제공한다. 특히, 소행렬식 기반의 비선형 조합은 고차 모멘트와 연관된 비정상성(intermittency) 현상을 포착할 수 있으며, 실험·시뮬레이션 데이터에서의 작은 비대칭성도 민감하게 드러낸다. 따라서 제안된 확률 항등식은 난류 이론에서 등방성 가정의 검증, 모델링의 제약조건 설정, 그리고 복합 물리 현상(예: 마그네토수치류, 스칼라 운반)에서의 텐서 구조 분석에 광범위하게 활용될 수 있다.